Discrete element method can effectively simulate the discontinuity,inhomogeneity and large deformation and failure of rock and soil.Based on the innovative matrix computing of the discrete element method,the highperfo...Discrete element method can effectively simulate the discontinuity,inhomogeneity and large deformation and failure of rock and soil.Based on the innovative matrix computing of the discrete element method,the highperformance discrete element software MatDEM may handle millions of elements in one computer,and enables the discrete element simulation at the engineering scale.It supports heat calculation,multi-field and fluidsolid coupling numerical simulations.Furthermore,the software integrates pre-processing,solver,postprocessing,and powerful secondary development,allowing recompiling new discrete element software.The basic principles of the DEM,the implement and development of the MatDEM software,and its applications are introduced in this paper.The software and sample source code are available online(http://matdem.com).展开更多
Numerical simulation is dominant in solving partial differential equations(PDEs),but balancing fine-grained grids with low computational costs is challenging.Recently,solving PDEs with neural networks(NNs)has gained i...Numerical simulation is dominant in solving partial differential equations(PDEs),but balancing fine-grained grids with low computational costs is challenging.Recently,solving PDEs with neural networks(NNs)has gained interest,yet cost-effectiveness and high accuracy remain a challenge.This work introduces a novel paradigm for solving PDEs,called multi-scale neural computing(MSNC),considering spectral bias of NNs and local approximation properties in the finite difference method(FDM).The MSNC decomposes the solution with a NN for efficient capture of global scale and the FDM for detailed description of local scale,aiming to balance costs and accuracy.Demonstrated advantages include higher accuracy(10 times for 1D PDEs,20 times for 2D PDEs)and lower costs(4 times for 1D PDEs,16 times for 2D PDEs)than the standard FDM.The MSNC also exhibits stable convergence and rigorous boundary condition satisfaction,showcasing the potential for hybrid of NN and numerical method.展开更多
Three-dimensional(3 D)reconstruction of icosahedral viruses has played a crucial role in the development of cryoelectron microscopy single-particle reconstruction,with many cryo-electron microscopy techniques first es...Three-dimensional(3 D)reconstruction of icosahedral viruses has played a crucial role in the development of cryoelectron microscopy single-particle reconstruction,with many cryo-electron microscopy techniques first established for structural studies of icosahedral viruses,owing to their high symmetry and large mass.This review summarizes the computational methods for icosahedral and symmetry-mismatch reconstruction of viruses,as well as the likely challenges and bottlenecks in virus reconstruction,such as symmetry mismatch reconstruction,contrast transformation function(CTF)correction,and particle distortion.展开更多
To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Un...To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks....Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.展开更多
An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multi...An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.展开更多
This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed...This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.展开更多
The main aim in this research is comparison the parameters of some storm events in the watershed using two loss models in Unit hydrograph method by HEC-HMS. SCS Curve Number and Green-Ampt methods by developing loss m...The main aim in this research is comparison the parameters of some storm events in the watershed using two loss models in Unit hydrograph method by HEC-HMS. SCS Curve Number and Green-Ampt methods by developing loss model as a major component in runoff and flood modeling. The study is conducted in the Kuala Lumpurwatershed with674 km2 area located in Klang basin inMalaysia. The catchment delineation is generated for the Klang watershed to get sub-watershed parameters by using HEC-GeoHMS extension in ARCGIS. Then all the necessary parameters are assigned to the models applied in this study to run the runoff and flood model. The results showed that there was no significant difference between the SCS-CN and Green-Ampt loss method applied in the Klang watershed. Estimated direct runoff and Peak discharge (r = 0.98) indicates a statistically positive correlations between the results of the study. And also it has been attempted to use objective functions in HEC-HMS (percent error peaks and percent error volume) to classify the methods. The selection of best method is on the base of considering least difference between the results of simulation to observed events in hydrographs so that it can address which model is suit for runoff-flood simulation in Klang watershed. Results showed that SCS CN and Green-Ampt methods, in three events by fitting with percent error in peak and percent error in volume had no significant difference.展开更多
Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the productio...Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.展开更多
Among the most basic challenges of hydrology are the prediction and quantification of catchment surface runoff. The runoff curve number （CN） is a key factor in determining runoff in the SCS （Soil Conservation Serv...Among the most basic challenges of hydrology are the prediction and quantification of catchment surface runoff. The runoff curve number （CN） is a key factor in determining runoff in the SCS （Soil Conservation Service） based hydrologic modeling method. The traditional SCS-CN method for calculating the composite curve number is very tedious and consumes a major portion of the hydrologic modeling time. Therefore, geographic information systems （G/S） are now being used in combination with the SCS-CN method. This paper assesses the modeling of flow in West Bank catchments using the GIS-based SCS-CN method. The West Bank, Palestine, is characterized as an arid to semi-arid region with annual rainfall depths ranging between 100 mm in the vicinity of the Jordan River to 700 mm in the mountains extending across the central parts of the region. The estimated composite curve number for the entire West Bank is about 50 assuming dry conditions. This paper clearly demonstrates that the integration of GIS with the SCS-CN method provides a powerful tool for estimating runoff volumes in West Bank catchments, representing arid to semi-arid catchments of Palestine.展开更多
A scheme for general purposed FDTD visual scientific computing software is introduced in this paper using object-oriented design (OOD) method. By abstracting the parameters of FDTD grids to an individual class and sep...A scheme for general purposed FDTD visual scientific computing software is introduced in this paper using object-oriented design (OOD) method. By abstracting the parameters of FDTD grids to an individual class and separating from the iteration procedure, the visual software can be adapted to more comprehensive computing problems. Real-time gray degree graphic and wave curve of the results can be achieved using DirectX technique. The special difference equation and data structure in dispersive medium are considered, and the peculiarity of parameters in perfectly matched layer are also discussed.展开更多
In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and di...In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.展开更多
The computational methods of a typical dynamic mathematical model that can describe the differential element and the inertial element for the system simulation are researched. The stability of numerical solutions of t...The computational methods of a typical dynamic mathematical model that can describe the differential element and the inertial element for the system simulation are researched. The stability of numerical solutions of the dynamic mathematical model is researched. By means of theoretical analysis, the error formulas, the error sign criteria and the error relationship criterion of the implicit Euler method and the trapezoidal method are given, the dynamic factor affecting the computational accuracy has been found, the formula and the methods of computing the dynamic factor are given. The computational accuracy of the dynamic mathematical model like this can be improved by use of the dynamic factor.展开更多
A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-unifor...A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.展开更多
Attempting to find a fast computing method to DHT （distinguished hyperbolic trajectory）, this study first proves that the errors of the stable DHT can be ignored in normal direction when they are computed as the tra...Attempting to find a fast computing method to DHT （distinguished hyperbolic trajectory）, this study first proves that the errors of the stable DHT can be ignored in normal direction when they are computed as the trajectories extend. This conclusion means that the stable flow with perturbation will approach to the real trajectory as it extends over time. Based on this theory and combined with the improved DHT computing method, this paper reports a new fast computing method to DHT, which magnifies the DHT computing speed without decreasing its accuracy.展开更多
Based on the efficient hybrid methods for solving initial value problems of stiff ODEs, this paper derives a parallel scheme that can be used to solve the problems on parallel computers with N processors, and discusse...Based on the efficient hybrid methods for solving initial value problems of stiff ODEs, this paper derives a parallel scheme that can be used to solve the problems on parallel computers with N processors, and discusses the iteratively B-convergence of the Newton iterative process, finally, the paper provides some numberical results which show that the parallel scheme is highly efficient as N is not too large.展开更多
In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hover...In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.展开更多
基金Financial supports from the Natural Science Foundation of China(41761134089,41977218)Six Talent Peaks Project of Jiangsu Province(RJFW-003)the Fundamental Research Funds for the Central Universities(14380103)are gratefully acknowledged.
文摘Discrete element method can effectively simulate the discontinuity,inhomogeneity and large deformation and failure of rock and soil.Based on the innovative matrix computing of the discrete element method,the highperformance discrete element software MatDEM may handle millions of elements in one computer,and enables the discrete element simulation at the engineering scale.It supports heat calculation,multi-field and fluidsolid coupling numerical simulations.Furthermore,the software integrates pre-processing,solver,postprocessing,and powerful secondary development,allowing recompiling new discrete element software.The basic principles of the DEM,the implement and development of the MatDEM software,and its applications are introduced in this paper.The software and sample source code are available online(http://matdem.com).
基金supported by the National Natural Science Foundation of China(Grant No.92152301).
文摘Numerical simulation is dominant in solving partial differential equations(PDEs),but balancing fine-grained grids with low computational costs is challenging.Recently,solving PDEs with neural networks(NNs)has gained interest,yet cost-effectiveness and high accuracy remain a challenge.This work introduces a novel paradigm for solving PDEs,called multi-scale neural computing(MSNC),considering spectral bias of NNs and local approximation properties in the finite difference method(FDM).The MSNC decomposes the solution with a NN for efficient capture of global scale and the FDM for detailed description of local scale,aiming to balance costs and accuracy.Demonstrated advantages include higher accuracy(10 times for 1D PDEs,20 times for 2D PDEs)and lower costs(4 times for 1D PDEs,16 times for 2D PDEs)than the standard FDM.The MSNC also exhibits stable convergence and rigorous boundary condition satisfaction,showcasing the potential for hybrid of NN and numerical method.
基金Project supported by the National Key R&D Program of China(Grant No.2016YFA0501100)the National Natural Science Foundation of China(Grant Nos.91530321,31570742,and 31570727)Science and Technology Planning Project of Hunan Province,China(Grant No.2017RS3033)
文摘Three-dimensional(3 D)reconstruction of icosahedral viruses has played a crucial role in the development of cryoelectron microscopy single-particle reconstruction,with many cryo-electron microscopy techniques first established for structural studies of icosahedral viruses,owing to their high symmetry and large mass.This review summarizes the computational methods for icosahedral and symmetry-mismatch reconstruction of viruses,as well as the likely challenges and bottlenecks in virus reconstruction,such as symmetry mismatch reconstruction,contrast transformation function(CTF)correction,and particle distortion.
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.2023YJS053)the National Natural Science Foundation of China(Grant No.52278386).
文摘To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
文摘Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.
基金supported by the National Natural Science Foundation(10721062,11072051,90715037,10728205,91015003, 51021140004)the Program of Introducing Talents of Discipline to Universities(B08014)the National Key Basic Research Special Foundation of China(2010CB832704).
文摘An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.
基金National Natural Science Foundation of China(51076006)
文摘This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.
文摘The main aim in this research is comparison the parameters of some storm events in the watershed using two loss models in Unit hydrograph method by HEC-HMS. SCS Curve Number and Green-Ampt methods by developing loss model as a major component in runoff and flood modeling. The study is conducted in the Kuala Lumpurwatershed with674 km2 area located in Klang basin inMalaysia. The catchment delineation is generated for the Klang watershed to get sub-watershed parameters by using HEC-GeoHMS extension in ARCGIS. Then all the necessary parameters are assigned to the models applied in this study to run the runoff and flood model. The results showed that there was no significant difference between the SCS-CN and Green-Ampt loss method applied in the Klang watershed. Estimated direct runoff and Peak discharge (r = 0.98) indicates a statistically positive correlations between the results of the study. And also it has been attempted to use objective functions in HEC-HMS (percent error peaks and percent error volume) to classify the methods. The selection of best method is on the base of considering least difference between the results of simulation to observed events in hydrographs so that it can address which model is suit for runoff-flood simulation in Klang watershed. Results showed that SCS CN and Green-Ampt methods, in three events by fitting with percent error in peak and percent error in volume had no significant difference.
基金National Natural Science Foundation of China(U2004176,22008055)Technology Research Project of Henan Province(232102240034)are gratefully acknowledged.
文摘Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.
基金supported by the GLOWA-JR Project of the German Federal Ministry of Education and Research (BMBF)
文摘Among the most basic challenges of hydrology are the prediction and quantification of catchment surface runoff. The runoff curve number （CN） is a key factor in determining runoff in the SCS （Soil Conservation Service） based hydrologic modeling method. The traditional SCS-CN method for calculating the composite curve number is very tedious and consumes a major portion of the hydrologic modeling time. Therefore, geographic information systems （G/S） are now being used in combination with the SCS-CN method. This paper assesses the modeling of flow in West Bank catchments using the GIS-based SCS-CN method. The West Bank, Palestine, is characterized as an arid to semi-arid region with annual rainfall depths ranging between 100 mm in the vicinity of the Jordan River to 700 mm in the mountains extending across the central parts of the region. The estimated composite curve number for the entire West Bank is about 50 assuming dry conditions. This paper clearly demonstrates that the integration of GIS with the SCS-CN method provides a powerful tool for estimating runoff volumes in West Bank catchments, representing arid to semi-arid catchments of Palestine.
基金This project was supported by the National Natural Science Foundation (No. 69831020).
文摘A scheme for general purposed FDTD visual scientific computing software is introduced in this paper using object-oriented design (OOD) method. By abstracting the parameters of FDTD grids to an individual class and separating from the iteration procedure, the visual software can be adapted to more comprehensive computing problems. Real-time gray degree graphic and wave curve of the results can be achieved using DirectX technique. The special difference equation and data structure in dispersive medium are considered, and the peculiarity of parameters in perfectly matched layer are also discussed.
基金national natural science foundation natural science foundation of Gansu province.
文摘In this paper, based on the implicit Runge-Kutta(IRK) methods, we derive a class of parallel scheme that can be implemented on the parallel computers with Ns(N is a positive even number) processors efficiently, and discuss the iteratively B-convergence of the Newton iterative process for solving the algebraic equations of the scheme, secondly we present a strategy providing initial values parallelly for the iterative process. Finally, some numerical results show that our parallel scheme is higher efficient as N is not so large.
文摘The computational methods of a typical dynamic mathematical model that can describe the differential element and the inertial element for the system simulation are researched. The stability of numerical solutions of the dynamic mathematical model is researched. By means of theoretical analysis, the error formulas, the error sign criteria and the error relationship criterion of the implicit Euler method and the trapezoidal method are given, the dynamic factor affecting the computational accuracy has been found, the formula and the methods of computing the dynamic factor are given. The computational accuracy of the dynamic mathematical model like this can be improved by use of the dynamic factor.
基金Project supported by the National Natural Science Foundation of China (Grant No. 51075311)
文摘A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.
基金supported by the National Natural Science Foundation of China (Grant No. 60872159)
文摘Attempting to find a fast computing method to DHT （distinguished hyperbolic trajectory）, this study first proves that the errors of the stable DHT can be ignored in normal direction when they are computed as the trajectories extend. This conclusion means that the stable flow with perturbation will approach to the real trajectory as it extends over time. Based on this theory and combined with the improved DHT computing method, this paper reports a new fast computing method to DHT, which magnifies the DHT computing speed without decreasing its accuracy.
文摘Based on the efficient hybrid methods for solving initial value problems of stiff ODEs, this paper derives a parallel scheme that can be used to solve the problems on parallel computers with N processors, and discusses the iteratively B-convergence of the Newton iterative process, finally, the paper provides some numberical results which show that the parallel scheme is highly efficient as N is not too large.
文摘In order to overcome the efficiency problem of the conventional gradient-based optimal design method,a highly-efficient viscous adjoint-based RANS equations method is applied to the aerodynamic optimal design of hovering rotor airfoil.The C-shaped body-fitted mesh is firstly automatically generated around the airfoil by solving the Poisson equations,and the Navier-Stokes(N-S)equations combined with Spalart-Allmaras(S-A)one-equation turbulence model are used as the governing equations to acquire the reliable flowfield variables.Then,according to multi-constrained characteristics of the optimization of high lift/drag ratio for hovering rotor airfoil,its corresponding adjoint equations,boundary conditions and gradient expressions are newly derived.On these bases,two representative rotor airfoils,NACA0012 airfoil and SC1095 airfoil,are selected as numerical examples to optimize their synthesized aerodynamic characteristics about lift/drag ratio in hover,and better aerodynamic performance of optimal airfoils are obtained compared with the baseline.Furthermore,the new designed rotor with the optimized rotor airfoil has better hover aerodynamic characteristics compared with the baseline rotor.In contrast to the baseline airfoils optimized by the finite difference method,it is demonstrated that the adjoint optimal algorithm itself is practical and highly-efficient for the aerodynamic optimization of hover rotor airfoil.
