A simplified simulation method based on the FDTD technique that can handle active devices is proposed. This method well suits the electrical crosstalk analysis of multi-channel integrated, opto-electronic mixed module...A simplified simulation method based on the FDTD technique that can handle active devices is proposed. This method well suits the electrical crosstalk analysis of multi-channel integrated, opto-electronic mixed modules. We apply this method to an 8-channel integrated super-compact high-sensitivity optical module. The results show good agreement between simulations and measurements.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through poro...In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.展开更多
The aim of this work is to build a 3D numerical model to study the characteristics of passive drag on competitive swimmers taking into account the impact of the free surface. This model solves the 3D incompressible Na...The aim of this work is to build a 3D numerical model to study the characteristics of passive drag on competitive swimmers taking into account the impact of the free surface. This model solves the 3D incompressible Navier-Stokes equations using RNG k-ε turbulence closure. The volume of fluid（VOF） method is used to locate the free surface. The 3D virtual model is created by Computer Aided Industrial Design（CAID） software, Rhinoceros. Firstly, a specific posture of swimming is studied. The simulation results are in good agreement with the data from mannequin towing experiments. The effects of a swimmer＇s arms and legs positions on swimming performance are then studied. Finally, it is demonstrated that the present method is capable of simulating gliding near the free surface.展开更多
Many studies revealed that the Earth medium＇s lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method ar...Many studies revealed that the Earth medium＇s lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method are adopted to quantify the effects of lateral heterogeneity caused by material parameters and fault dip angle on the co- and postseismic deformation in the near- and far-field. Our results show that： 1） the medium＇s lateral heterogeneity does affect the co-seismic deformation, with the effects increasing with the medium＇s lateral heterogeneity caused by material parameters; 2） the Lame parameters play a more dominant role than density in the effects caused by lateral heterogeneity; 3） when a fault＇s dip angle is smaller than 90, the effects of the medium＇s lateral heterogeneity on the hanging wall are greater than on the footwall; 4） the impact of lateral heterogeneity caused by the viscosity coefficient on the post-seismic deformation can affect a large area, including the near- and far-field.展开更多
The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in...The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in a mountain railway,we established an aerodynamic model involving a train exiting the tunnel,and verified it in the Fluent environment.Overset mesh technology was adopted to characterize the train’s movement.The flow field involving the train,tunnel,and crosswinds was simulated using the Reynolds-averaged turbulence model.Then,we built a comprehensive train-track coupled dynamic model considering the influences of ADLs,to investigate the vehicles’dynamic responses.The aerodynamics and dynamic behaviors of the train when affected by crosswinds with different velocities and directions are analyzed and discussed.The results show that the near-wall side crosswind leads to sharper variations in ADLs than the far-wall side crosswind.The leading vehicle suffers from more severe ADLs than other vehicles,which worsens the wheel-rail interaction and causes low-frequency vibration of the car body.When the crosswind velocity exceeds 20 m/s,significant wheel-rail impacts occur,and the running safety of the train worsens rapidly.展开更多
Grouting is a widely applied technique for reinforcing fractured zones in deep soft rock tunnels.By infiltrating rock fissures,slurry materials enhance structural integrity and improve the overall stability of the sur...Grouting is a widely applied technique for reinforcing fractured zones in deep soft rock tunnels.By infiltrating rock fissures,slurry materials enhance structural integrity and improve the overall stability of the surrounding rock.The performance of grouting is primarily governed by the flow behavior and diffusion extent of the slurry.This review considers recent advances in the theory and methodology of slurry flow and diffusion in fractured rock.It examines commonly used grout materials,including cement-based,chemical,and composite formulations,each offering distinct advantages for specific geological conditions.Themechanisms of reinforcement vary significantly across materials,requiring tailored application strategies.The rheological properties of grouting slurries,particularly cement-based types,have been widely modeled using classical constitutive approaches.However,the influence of time-and space-dependent viscosity evolution on slurry behavior remains underexplored.Experimental studies have provided valuable insights into slurry diffusion,yet further research is needed to capture real-time behavior under multi-scale and multi-physics coupling conditions.Similarly,current numerical simulations are largely limited to twoand three-dimensional models of single-fracture flow.These models often neglect the complexity of fracture networks and geological heterogeneity,highlighting a need for more realistic and integrated simulation frameworks.Future research should focus on:(1)fine-scale modeling of slurry hydration and mechanical reinforcement processes;(2)cross-scale analysis of slurry flow under coupled thermal,hydraulic,andmechanical fields;and(3)development of realtime,three-dimensional dynamic simulation tools to capture the full grouting process.These efforts will strengthen the theoretical foundation and practical effectiveness of grouting in complex underground environments.展开更多
A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck ...A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.展开更多
For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is develo...For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.展开更多
A computational modeling for the sheet cavitating flows is presented. The cavitation model is implemented in a viscous Navier-Stokes solver. The cavity interface and shape are determined using an iterative procedure m...A computational modeling for the sheet cavitating flows is presented. The cavitation model is implemented in a viscous Navier-Stokes solver. The cavity interface and shape are determined using an iterative procedure matching the cavity surface to a constant pressure boundary. The pressure distribution, as well as its gradient on the wall, is taken into account in updating the cavity shape iteratively. Numerical computations are performed for the sheet cavitating flows at a range of cavitation numbers across the hemispheric headform/cylinder body with different grid numbers. The influence of the relaxation factor in the cavity shape updating scheme for the algorithm accuracy and reliability is conducted through comparison with other two cavity shape updating numerical schemes. The results obtained are reasonable and the iterative procedure of cavity shape updating is quite stable, which demonstrate the superiority of the proposed cavitation model and algorithms.展开更多
For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or f...For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.展开更多
In this paper,we present an overview of numerical simulation methods for the flow around typical underwater vehicles at high Reynolds numbers,which highlights the dominant flow structures in different regions of inter...In this paper,we present an overview of numerical simulation methods for the flow around typical underwater vehicles at high Reynolds numbers,which highlights the dominant flow structures in different regions of interest.This overview covers the forebody,midbody,stern,wake region,and appendages and summarizes flow phenomena,including laminar-to-turbulent transition,turbulent boundary layers,flow under the influence of curvatures,wake interactions,and all associated complex vortex structures.Furthermore,the current issues and challenges of capturing these flow structures are addressed.This overview provides a deep insight into the use of numerical simulation methods,including the Reynolds-averaged Navier–Stokes(RANS)method,large eddy simulation(LES)method,and the hybrid RANS/LES method,and evaluates their applicability in capturing detailed flow features.展开更多
A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by t...A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method （LBM）, Particle-fluid coupling is realized by an immersed boundary method （IBM）, which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.展开更多
The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and...The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and most importantly, no paradoxical parameters appeared in the model. Moreover, it can simulate different microstructure patterns owing to the use of a new different free energy function for the simulation of morphologies of polymer. The new free energy function considers both the cases of T〈Tm and T≥Tm, which is more reasonable than that in published literatures that all ignored the T≥Tm case. In order to show the validity of the modified model, the finite difference method is used to solve the model and different crystallization morphologies during the solidification process of isotactic polystyrene are obtained under different conditions. Numerical results show that the growth rate of the initial secondary arms is obviously increased as the anisotropy strength increases. But the anisotropy strength seems to have no apparent effect on the global growth rate. The whole growth process of the dendrite depends mainly upon the latent heat and the latent heat has a direct effect on the tip radius and tip velocity of side branches.展开更多
Three kinds of models based on the same flexible pipe with 8 layers have been separately created to investigate the effects of different modeling approaches on numerical simulation results of finite element(FE)models ...Three kinds of models based on the same flexible pipe with 8 layers have been separately created to investigate the effects of different modeling approaches on numerical simulation results of finite element(FE)models for unbonded flexible pipes.Then the mechanical property of the unbonded flexible pipe under tension,torsion and bending load has been analyzed and compared via ABAQUS software on the basis of three created models.The research shows that different modeling methods of flexible pipes make a great difference in the results.Especially,modeling simplifications of the carcass and pressure armor have a great impact on the accuracy of the results.Model 3,in which the carcass is simulated by spiral isot ropic shell and other layers are Simula ted by solid element,possesses good adaptability,which has been proved by comparing the experiment data and other models.This paper can offer a reference for the FE modeling methods,selection and mechanical property analysis of unbonded flexible pipe.展开更多
The thin-walled tube flexure（TWTF） hinges have important potential application value in the deployment mechanisms of satellite and solar array, but the optimal design of the TWTF hinges haven＇t been completely solv...The thin-walled tube flexure（TWTF） hinges have important potential application value in the deployment mechanisms of satellite and solar array, but the optimal design of the TWTF hinges haven＇t been completely solved, which restricts their applications. An optimal design method for the qusai-static folding and deploying of TWTF hinges with double slots is presented based on the response surface theory. Firstly, the full factorial method is employed to design of the experiments. Then, the finite element models of the TWTF hinges with double slots are constructed to simulate the qusai-static folding and deploying non-linear analysis. What＇s more, the mathematical model of the TWTF flexure hinge quasi-static folding and deploying properties are derived by the response surface method. Considering of small mass and high stability, the peak moment of quasi-static folding and deploying as well as the lightless are set as the objectives to get the optimal performances. The relative errors of the objectives between the optimal design results and the FE analysis results are less than 7%, which demonstrates the precision of the surrogate models. Lastly, the parameter study shows that both the slots length and the slots width both have significant effects to the peak moment of quasi-static folding and deploying of TWTF hinges with double slots. However, the maximum Mises stress of quasi-static folding is more sensitive to the slots length than the slots width. The proposed research can be applied to optimize other thin-walled flexure hinges under quasi-static folding and deploying, which is of great importance to design of flexure hinges with high stability and low stress.展开更多
This paper focuses on the analytical derivation and the numerical simulation analyses to predict the interaction influences between a landslide and a new tunnel in mountain areas. Based on the slip-line theory, the di...This paper focuses on the analytical derivation and the numerical simulation analyses to predict the interaction influences between a landslide and a new tunnel in mountain areas. Based on the slip-line theory, the disturbance range induced by tunneling and the minimum safe distance between the tunnel vault and the sliding belt are obtained in consideration of the mechanical analyses of relaxed rocks over the tunnel opening. The influence factors for the minimum safe crossing distance are conducted,including the tunnel radius, the friction angle of surrounding rocks, the inclination angle of sliding belt,and the friction coefficient of surrounding rocks. Secondly, taking account of the compressive zone and relaxed rocks caused by tunneling, the Sarma method is employed to calculate the safety factor of landslide. Finally, the analytical solutions for interaction between the tunnel and the landslide are compared with a series of numerical simulations, considering the cases for different perpendicular distances between the tunnel vault and the sliding belt. Resultsshow that the distance between the tunnel vault and the slip zone has significant influence on the rock stress and strain. For the case of the minimum crossing distance, a plastic zone in the landslide traversed by tunneling would be formed with rather large range, which seriously threatens the stability of landslide. This work demonstrates that the minimum safe crossing distance obtained from numerical simulation is in a good agreement with that calculated by the proposed analytical solutions.展开更多
The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the...The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.展开更多
As one of the main reasons causing leakage heat load in a refrigerator,mass and heat transfer through refrigerator door seal is of great importance to be studied.In this paper,a model is presented for numerical simula...As one of the main reasons causing leakage heat load in a refrigerator,mass and heat transfer through refrigerator door seal is of great importance to be studied.In this paper,a model is presented for numerical simulation of mass and heat transfer process through refrigerator door seal,and an experiment apparatus is designed and set up as well for comparison.A two-dimensional model and tracer gas method are used in simulation and experiment,respectively.It can be found that the relative deviations of air infiltration rate between the simulated results and experimental results were less than 1%,and the temperature difference errors at two special points of the door seal were less than 2.03℃.In conclusion,the simulated results are in good agreement with the experimental results.This paper initially sets up a model that can accurately simulate the heat and mass transfer through the refrigerator door seal,and the model can be used in refrigerator door seal optimization research in the follow-up study.展开更多
文摘A simplified simulation method based on the FDTD technique that can handle active devices is proposed. This method well suits the electrical crosstalk analysis of multi-channel integrated, opto-electronic mixed modules. We apply this method to an 8-channel integrated super-compact high-sensitivity optical module. The results show good agreement between simulations and measurements.
基金supported by the Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金Projects 50534090 and 50674090 supported by the National Natural Science Foundation of China2005CB221503 by the National Key Basic ResearchDevelopment Program (973 Program)
文摘In order to simulate field distribution rules,mathematical models for 3-D air flows and gas transfer in the goaf of the coal mining are established,based on theories of permeability and dynamic dispersion through porous media. A gas dispersion equation in a 3-D field is calculated by use of numerical method on a weighted upstream multi-element balance. Based on data of an example with a U type ventilation mode,surface charts of air pressure distribution and gas concentration are drawn by Graphtool software. Finally,a comparison between actually measured results in the model test and the numerical simulation results is made to proves the numerical implementation feasible.
基金financially supported by the Fundamental Research Funds for the Central Universities of China
文摘The aim of this work is to build a 3D numerical model to study the characteristics of passive drag on competitive swimmers taking into account the impact of the free surface. This model solves the 3D incompressible Navier-Stokes equations using RNG k-ε turbulence closure. The volume of fluid（VOF） method is used to locate the free surface. The 3D virtual model is created by Computer Aided Industrial Design（CAID） software, Rhinoceros. Firstly, a specific posture of swimming is studied. The simulation results are in good agreement with the data from mannequin towing experiments. The effects of a swimmer＇s arms and legs positions on swimming performance are then studied. Finally, it is demonstrated that the present method is capable of simulating gliding near the free surface.
基金co-supported by the National Natural Science Foundation of China (41431069)the State Key Development Program for Basic Research of China (2013CB733304, 2013CB733303)+1 种基金the Doctoral Fund of Ministry of Education of China (20110141130010)China Postdoctoral Science Foundation funded project (2013M542062)
文摘Many studies revealed that the Earth medium＇s lateral heterogeneity can cause considerable effects on the co- and post-seismic deformation field. In this study, the threedimensional finite element numerical method are adopted to quantify the effects of lateral heterogeneity caused by material parameters and fault dip angle on the co- and postseismic deformation in the near- and far-field. Our results show that： 1） the medium＇s lateral heterogeneity does affect the co-seismic deformation, with the effects increasing with the medium＇s lateral heterogeneity caused by material parameters; 2） the Lame parameters play a more dominant role than density in the effects caused by lateral heterogeneity; 3） when a fault＇s dip angle is smaller than 90, the effects of the medium＇s lateral heterogeneity on the hanging wall are greater than on the footwall; 4） the impact of lateral heterogeneity caused by the viscosity coefficient on the post-seismic deformation can affect a large area, including the near- and far-field.
基金National Natural Science Foundation of China(No.52388102)New Cornerstone Science Foundation through the Xplorer Prize.
文摘The dynamic performance of high-speed trains is significantly influenced by sudden changes in aerodynamic loads(ADLs)when exiting a tunnel in a windy environment.Focusing on a double-track tunnel under construction in a mountain railway,we established an aerodynamic model involving a train exiting the tunnel,and verified it in the Fluent environment.Overset mesh technology was adopted to characterize the train’s movement.The flow field involving the train,tunnel,and crosswinds was simulated using the Reynolds-averaged turbulence model.Then,we built a comprehensive train-track coupled dynamic model considering the influences of ADLs,to investigate the vehicles’dynamic responses.The aerodynamics and dynamic behaviors of the train when affected by crosswinds with different velocities and directions are analyzed and discussed.The results show that the near-wall side crosswind leads to sharper variations in ADLs than the far-wall side crosswind.The leading vehicle suffers from more severe ADLs than other vehicles,which worsens the wheel-rail interaction and causes low-frequency vibration of the car body.When the crosswind velocity exceeds 20 m/s,significant wheel-rail impacts occur,and the running safety of the train worsens rapidly.
基金funded by the National Natural Science Foundation of China[U22A20234]Hubei Province key research and development project[2023BCB121]Wuhan innovation supporting projects[2023020201010079].
文摘Grouting is a widely applied technique for reinforcing fractured zones in deep soft rock tunnels.By infiltrating rock fissures,slurry materials enhance structural integrity and improve the overall stability of the surrounding rock.The performance of grouting is primarily governed by the flow behavior and diffusion extent of the slurry.This review considers recent advances in the theory and methodology of slurry flow and diffusion in fractured rock.It examines commonly used grout materials,including cement-based,chemical,and composite formulations,each offering distinct advantages for specific geological conditions.Themechanisms of reinforcement vary significantly across materials,requiring tailored application strategies.The rheological properties of grouting slurries,particularly cement-based types,have been widely modeled using classical constitutive approaches.However,the influence of time-and space-dependent viscosity evolution on slurry behavior remains underexplored.Experimental studies have provided valuable insights into slurry diffusion,yet further research is needed to capture real-time behavior under multi-scale and multi-physics coupling conditions.Similarly,current numerical simulations are largely limited to twoand three-dimensional models of single-fracture flow.These models often neglect the complexity of fracture networks and geological heterogeneity,highlighting a need for more realistic and integrated simulation frameworks.Future research should focus on:(1)fine-scale modeling of slurry hydration and mechanical reinforcement processes;(2)cross-scale analysis of slurry flow under coupled thermal,hydraulic,andmechanical fields;and(3)development of realtime,three-dimensional dynamic simulation tools to capture the full grouting process.These efforts will strengthen the theoretical foundation and practical effectiveness of grouting in complex underground environments.
基金financial support from the National Natural Science Foundation of China (No. 51204181)the Research Fund for the Doctoral Program of Higher Education of China (No. 20110095120004)+2 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe Fundamental Research Funds for the Central Universities (Nos. 2011QNA10 and 2010QNB17)the China Postdoctoral Science Foundation (No. 20110491485)
文摘A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.
基金financially supported by the Program for New Century Excellent Talents in University(No.NCET-13-0229,NCET-09-0396)the National Science & Technology Key Projects of Numerical Control(No.2012ZX04010-031,2012ZX0412-011)the National High Technology Research and Development Program("863"Program)of China(No.2013031003)
文摘For the sake of a more accurate shell boundary and calculation of radiation heat transfer in the Directional Solidification(DS) process, a radiation heat transfer model based on the Finite Element Method(FEM)is developed in this study. Key technologies, such as distinguishing boundaries automatically, local matrix and lumped heat capacity matrix, are also stated. In order to analyze the effect of withdrawing rate on DS process,the solidification processes of a complex superalloy turbine blade in the High Rate Solidification(HRS) process with different withdrawing rates are simulated; and by comparing the simulation results, it is found that the most suitable withdrawing rate is determined to be 5.0 mm·min^(-1). Finally, the accuracy and reliability of the radiation heat transfer model are verified, because of the accordance of simulation results with practical process.
基金This project is supported by the Grants from Specialized Research Fund for the Doctoral Program of Higher Education(SRFDP), China(No.20040698049)Natural Science Foundation of Xi'an Jiaotong University, China(No.2004).
文摘A computational modeling for the sheet cavitating flows is presented. The cavitation model is implemented in a viscous Navier-Stokes solver. The cavity interface and shape are determined using an iterative procedure matching the cavity surface to a constant pressure boundary. The pressure distribution, as well as its gradient on the wall, is taken into account in updating the cavity shape iteratively. Numerical computations are performed for the sheet cavitating flows at a range of cavitation numbers across the hemispheric headform/cylinder body with different grid numbers. The influence of the relaxation factor in the cavity shape updating scheme for the algorithm accuracy and reliability is conducted through comparison with other two cavity shape updating numerical schemes. The results obtained are reasonable and the iterative procedure of cavity shape updating is quite stable, which demonstrate the superiority of the proposed cavitation model and algorithms.
基金supported in part by the Higher Education Commission of Pakistan under PPCR programsupported by the National Magnetic Confinement Fusion Program under Grant No.2013GB104004Fundamental Research Fund for Chinese Central Universities
文摘For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.
基金Supported by the National Natural Science Foundation of China under Grant No.52131102.
文摘In this paper,we present an overview of numerical simulation methods for the flow around typical underwater vehicles at high Reynolds numbers,which highlights the dominant flow structures in different regions of interest.This overview covers the forebody,midbody,stern,wake region,and appendages and summarizes flow phenomena,including laminar-to-turbulent transition,turbulent boundary layers,flow under the influence of curvatures,wake interactions,and all associated complex vortex structures.Furthermore,the current issues and challenges of capturing these flow structures are addressed.This overview provides a deep insight into the use of numerical simulation methods,including the Reynolds-averaged Navier–Stokes(RANS)method,large eddy simulation(LES)method,and the hybrid RANS/LES method,and evaluates their applicability in capturing detailed flow features.
基金sponsored by Ministry of Finance under the grant ZDYZ2008-2National Key Science and Technology Project under the grant 2008ZX05014-003-006HZthe Chinese Academy of Sciences under the grant KGCX2-YW-124
文摘A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method （LBM）, Particle-fluid coupling is realized by an immersed boundary method （IBM）, which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.
基金This work is supported by the National Natural Science Foundation of China （No.11402210）, the Natural Science Foundation of Shanxi Province （No.2012011019-2）, and the Doctoral Fund of Taiyuan University of Science and Technology （No.20152024）.
文摘The irrationality of existing phase field model is analyzed and a modified phase-field model is proposed for polymer crystal growth, in which the parameters are obtained from real materials and very simple to use, and most importantly, no paradoxical parameters appeared in the model. Moreover, it can simulate different microstructure patterns owing to the use of a new different free energy function for the simulation of morphologies of polymer. The new free energy function considers both the cases of T〈Tm and T≥Tm, which is more reasonable than that in published literatures that all ignored the T≥Tm case. In order to show the validity of the modified model, the finite difference method is used to solve the model and different crystallization morphologies during the solidification process of isotactic polystyrene are obtained under different conditions. Numerical results show that the growth rate of the initial secondary arms is obviously increased as the anisotropy strength increases. But the anisotropy strength seems to have no apparent effect on the global growth rate. The whole growth process of the dendrite depends mainly upon the latent heat and the latent heat has a direct effect on the tip radius and tip velocity of side branches.
基金the National Natural Science Foundation of China(Nos.11627802 and 51709118)the China Postdoctoral Science Foundation(No.2017M612669)+2 种基金the Science and Technology Program of Guangzhou(No.201804010482)the Fundamental Research Funds for the Central Universities of China(No.2017BQ089)the Fund of State Key Laboratory of Ocean Engineering of Shanghai Jiao Tong University(No.1708)
文摘Three kinds of models based on the same flexible pipe with 8 layers have been separately created to investigate the effects of different modeling approaches on numerical simulation results of finite element(FE)models for unbonded flexible pipes.Then the mechanical property of the unbonded flexible pipe under tension,torsion and bending load has been analyzed and compared via ABAQUS software on the basis of three created models.The research shows that different modeling methods of flexible pipes make a great difference in the results.Especially,modeling simplifications of the carcass and pressure armor have a great impact on the accuracy of the results.Model 3,in which the carcass is simulated by spiral isot ropic shell and other layers are Simula ted by solid element,possesses good adaptability,which has been proved by comparing the experiment data and other models.This paper can offer a reference for the FE modeling methods,selection and mechanical property analysis of unbonded flexible pipe.
基金supported by National Natural Science Foundation ofChina(Grant No.50935002)
文摘The thin-walled tube flexure（TWTF） hinges have important potential application value in the deployment mechanisms of satellite and solar array, but the optimal design of the TWTF hinges haven＇t been completely solved, which restricts their applications. An optimal design method for the qusai-static folding and deploying of TWTF hinges with double slots is presented based on the response surface theory. Firstly, the full factorial method is employed to design of the experiments. Then, the finite element models of the TWTF hinges with double slots are constructed to simulate the qusai-static folding and deploying non-linear analysis. What＇s more, the mathematical model of the TWTF flexure hinge quasi-static folding and deploying properties are derived by the response surface method. Considering of small mass and high stability, the peak moment of quasi-static folding and deploying as well as the lightless are set as the objectives to get the optimal performances. The relative errors of the objectives between the optimal design results and the FE analysis results are less than 7%, which demonstrates the precision of the surrogate models. Lastly, the parameter study shows that both the slots length and the slots width both have significant effects to the peak moment of quasi-static folding and deploying of TWTF hinges with double slots. However, the maximum Mises stress of quasi-static folding is more sensitive to the slots length than the slots width. The proposed research can be applied to optimize other thin-walled flexure hinges under quasi-static folding and deploying, which is of great importance to design of flexure hinges with high stability and low stress.
基金financial support provided by Natural Science Foundation of China (Grant No. 51008188)by Shanghai Natural Science Foundation (Grant No. 15ZR1429400)+2 种基金by Open Project Program of State Key Laboratory Breeding Base of Mountain Bridge and Tunnel Engineering (Grant No. CQSLBF-Y15-1)by Open Project Program of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection (Grant No. SKLGP2015K015)by the Open Project Program of Key Laboratory of Geohazard Prevention of Hilly Mountains, Ministry of Land and Resources (Grant No. 2015k005)
文摘This paper focuses on the analytical derivation and the numerical simulation analyses to predict the interaction influences between a landslide and a new tunnel in mountain areas. Based on the slip-line theory, the disturbance range induced by tunneling and the minimum safe distance between the tunnel vault and the sliding belt are obtained in consideration of the mechanical analyses of relaxed rocks over the tunnel opening. The influence factors for the minimum safe crossing distance are conducted,including the tunnel radius, the friction angle of surrounding rocks, the inclination angle of sliding belt,and the friction coefficient of surrounding rocks. Secondly, taking account of the compressive zone and relaxed rocks caused by tunneling, the Sarma method is employed to calculate the safety factor of landslide. Finally, the analytical solutions for interaction between the tunnel and the landslide are compared with a series of numerical simulations, considering the cases for different perpendicular distances between the tunnel vault and the sliding belt. Resultsshow that the distance between the tunnel vault and the slip zone has significant influence on the rock stress and strain. For the case of the minimum crossing distance, a plastic zone in the landslide traversed by tunneling would be formed with rather large range, which seriously threatens the stability of landslide. This work demonstrates that the minimum safe crossing distance obtained from numerical simulation is in a good agreement with that calculated by the proposed analytical solutions.
基金supported by Scientific Research Foundation for Returned Scholars,Ministry of Education of China
文摘The dry-gas seal has been widely used in different industries. With increased spin speed of the rotator shaft, turbulence occurs in the gas film between the stator and rotor seal faces. For the micro-scale flow in the gas film and grooves, turbulence can change the pressure distribution of the gas film. Hence, the seal performance is influenced. However, turbulence effects and methods for their evaluation are not considered in the existing industrial designs of dry-gas seal. The present paper numerically obtains the turbulent flow fields of a spiral-groove dry-gas seal to analyze turbulence effects on seal performance. The direct numerical simulation （DNS） and Reynolds-averaged Navier-Stokes （RANS） methods are utilized to predict the velocity field properties in the grooves and gas film. The key performance parameter, open force, is obtained by integrating the pressure distribution, and the obtained result is in good agreement with the experimental data of other researchers. Very large velocity gradients are found in the sealing gas film because of the geometrical effects of the grooves. Considering turbulence effects, the calculation results show that both the gas film pressure and open force decrease. The RANS method underestimates the performance, compared with the DNS. The solution of the conventional Reynolds lubrication equation without turbulence effects suffers from significant calculation errors and a small application scope. The present study helps elucidate the physical mechanism of the hydrodynamic effects of grooves for improving and optimizing the industrial design or seal face pattern of a dry-gas seal.
基金Supported by the National Science Fund for Distinguished Young Scholar(51525604)111 project B16038
文摘As one of the main reasons causing leakage heat load in a refrigerator,mass and heat transfer through refrigerator door seal is of great importance to be studied.In this paper,a model is presented for numerical simulation of mass and heat transfer process through refrigerator door seal,and an experiment apparatus is designed and set up as well for comparison.A two-dimensional model and tracer gas method are used in simulation and experiment,respectively.It can be found that the relative deviations of air infiltration rate between the simulated results and experimental results were less than 1%,and the temperature difference errors at two special points of the door seal were less than 2.03℃.In conclusion,the simulated results are in good agreement with the experimental results.This paper initially sets up a model that can accurately simulate the heat and mass transfer through the refrigerator door seal,and the model can be used in refrigerator door seal optimization research in the follow-up study.
