In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solut...In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.展开更多
Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mec...Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems.展开更多
This study presents the first step of a research project that aims at using a three-dimensional （3D） hybridfinite-discrete element method （FDEM） to investigate the development of an excavation damaged zone（EDZ） ...This study presents the first step of a research project that aims at using a three-dimensional （3D） hybridfinite-discrete element method （FDEM） to investigate the development of an excavation damaged zone（EDZ） around tunnels in a clay shale formation known as Opalinus Clay. The 3D FDEM was first calibratedagainst standard laboratory experiments, including Brazilian disc test and uniaxial compression test. Theeffect of increasing confining pressure on the mechanical response and fracture propagation of the rockwas quantified under triaxial compression tests. Polyaxial （or true triaxial） simulations highlighted theeffect of the intermediate principal stress （s2） on fracture directions in the model： as the intermediateprincipal stress increased, fractures tended to align in the direction parallel to the plane defined by themajor and intermediate principal stresses. The peak strength was also shown to vary with changing s2. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the ...A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.展开更多
A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fractu...A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.展开更多
Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering appl...Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering applications. A three-dimensional anisotropic heterogeneous steady random seepage model was developed on the basis of the finite element method. A statistical analysis of the distribution characteristics of soil parameters sampled from the main embankment of the Yangtze River in the Southern Jingzhou zone of China was conducted. The Kolomogorov-Smirnov test verified the statistical hypothesis that the permeability coefficient tensor has a Gaussian distribution. With the help of numerical analysis of the stochastic seepage field using the developed model, various statistical and random characteristics of the stochastic seepage field of the main embankment of the Yangtze River in the Southern Jingzhou zone of China were investigated. The model was also examined with statistical testing. Through the introduction of random variation of the upstream and downstream water levels into the model, the effects of the boundary randomness due to variation of the downstream and upstream water levels on the variation of simulated results presented with a vector series of the random seepage field were analyzed. Furthermore, the combined influence of the variation of the soil permeability coefficient and such seepage resistance measures as the cut-off wall and relief ditch on the hydraulic head distribution was analyzed and compared with the results obtained by determinate analysis. Meanwhile, sensitivities of the hydraulic gradient and downstream exit height to the variation of boundary water level were studied. The validity of the simulated results was verified by stochastic testing and measured data. The developed model provides more detail and a full stochastic algorithm to characterize and analyze three-dimensional stochastic seepage field problems.展开更多
A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevo...A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.展开更多
The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-die...The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.展开更多
To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models...To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.展开更多
By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A specia...By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A special numerical code is introduced, in whicha lattice model and statistical approaches are used to simulate the initial heterogeneity ofmaterial properties. The size of displacement-load step is adap-tively determined so that only fewelements would fail in each load step. When the tensile principal strain in an element exceeds theultimate strain of this element, the element is considered broken and its Young's modulus is set tobe very low. Some important behaviors of heterogeneous brittle materials are indicated using thiscode. Load-displacement curves and figures of three-dimensional fracture patterns are alsonumerically obtained, which are similar to those observed in laboratory tests.展开更多
This present paper has proved the theorem of the Point Optimal Variable Successive Over Relaxation (OVSOR) method of the three-dimensional unsteady flow in the reservoir, and has put forward a formu- la for calculatin...This present paper has proved the theorem of the Point Optimal Variable Successive Over Relaxation (OVSOR) method of the three-dimensional unsteady flow in the reservoir, and has put forward a formu- la for calculating optimal parameters for OVSOR which vary with space points and time points. Using this method, internal memory of computer is the smallest, calculating work is the smallest, and calculating funds are the smallest. It is very easy to operate on microcomputers for three-dimensional res- ervoir simulation. The method is stable and convergent even if the time steps are taken to be large (for example, one year). The same applies for space steps. It is applicable both for homogeneous, isotropic porous mediums and for heterogeneous, anisotropic porous medium. On IBM microcomputers with internal memory of 512 thousand bytes, 8000 grid points may be cal- culated for three-dimensional simulation. It takes only two minutes to get convergence for one time step. It may be extended to three-dimensional heat conduction equation and three-dimensional simulation of the ground water flow. It looks much more advantageous for two-dimensional simulation.展开更多
We simulated three-dimensional heat transfer inside a horizontal rotating drum using the discrete element method and a thermal conduction model.The aim was to determine the effect of end-wall heating on thermal behavi...We simulated three-dimensional heat transfer inside a horizontal rotating drum using the discrete element method and a thermal conduction model.The aim was to determine the effect of end-wall heating on thermal behavior of a granular bed.The simulation showed that the end-wall heating significantly affects the axial temperature profile of the bed,particularly when the length-to-diameter ratio is low.Particles near the wall heated faster and became more thermally uniform than those in the center of the drum.The region affected by the end heating gradually increased over time.Increasing the rotation speed enhanced the heat conduction rate,and increasing the fill level reduced the mean temperature and thermal uniformity of the granular bed.Heat transfer was also simulated for drums with different length-to-diameter ratios.展开更多
Direct numerical simulation of coherent structures in the three-dimensional transitional jet with a moderate Reynolds number of 5000 was conducted. The finite volume method was used to discretize the governing equatio...Direct numerical simulation of coherent structures in the three-dimensional transitional jet with a moderate Reynolds number of 5000 was conducted. The finite volume method was used to discretize the governing equations in space; the low-storage, three-order Runge-Kutta scheme was used for time integration. The comparisons between the statistical results of the flow field; the related experimental data were performed to validate the reliability of the present numerical schemes. The emphasis was placed on the study of the spatial evolution of the three-dimensional coherent vortex structures as well as their interactions. It is found that the evolution of the spanwise vortex structures in three-dimensional space is similar to that in two-dimensional jet. The spanwise vortex structures are subject to three-dimensional instability; induce the formation of the streamwise; lateral vortex structures. Going with the breakup; mixing of the spanwise vortex structures, the streamwise; transverse vortex tubes also fall to pieces; the mixing arranged small-scale structures are formed in the flow field. Finally, the arrangement relationship among the spanwise, the streamwise; the lateral vortex structures was analyzed; their interactions were also discussed.展开更多
To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters for nasal structure, thirty three-dimensional, anatomically accurate representations of adult nasal cavity m...To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters for nasal structure, thirty three-dimensional, anatomically accurate representations of adult nasal cavity models were reconstructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated using the fluid dynamics with the finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to varying airflow distribution in the nasal cavities and the main airflow passed through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than a half of overall resistance. The characteristic model of nasal cavity was extracted based on the characteristic points and dimensions deducted from the original models. It showed that either the geometric structure or the air-flow field of the two kinds of model was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that properly represented the original model in research for nasal cavity.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10732010,10972010,and 11332002)
文摘In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.
基金Supported by the Research Starting Funds for Imported Talents of Ningxia University under Grant No BQD2012011
文摘Three-dimensional simulations of ferroelectric hysteresis and butterfly loops are carried out based on solving the time dependent Ginzburg-Landau equations using a finite volume method. The influence of externally mechanical loadings with a tensile strain and a compressive strain on the hysteresis and butterfly loops is studied numerically. Different from the traditional finite element and finite difference methods, the finite volume method is applicable to simulate the ferroelectric phase transitions and properties of ferroelectric materials even for more realistic and physical problems.
文摘This study presents the first step of a research project that aims at using a three-dimensional （3D） hybridfinite-discrete element method （FDEM） to investigate the development of an excavation damaged zone（EDZ） around tunnels in a clay shale formation known as Opalinus Clay. The 3D FDEM was first calibratedagainst standard laboratory experiments, including Brazilian disc test and uniaxial compression test. Theeffect of increasing confining pressure on the mechanical response and fracture propagation of the rockwas quantified under triaxial compression tests. Polyaxial （or true triaxial） simulations highlighted theeffect of the intermediate principal stress （s2） on fracture directions in the model： as the intermediateprincipal stress increased, fractures tended to align in the direction parallel to the plane defined by themajor and intermediate principal stresses. The peak strength was also shown to vary with changing s2. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.
基金Supported by the National Natural Science Foundation of China(51827804)CNPC Strategic Cooperation Science and Technology Major Project(ZLZX2020-01-05)Open Fund of State Key Laboratory of Rock Mechanics and Engineering(SKLGME021024).
文摘A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.
基金supported by the National Natural Science Foundation of China (Grant No. 50379046)the Doctoral Fund of the Ministry of Education of China (Grant No. A50221)
文摘Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering applications. A three-dimensional anisotropic heterogeneous steady random seepage model was developed on the basis of the finite element method. A statistical analysis of the distribution characteristics of soil parameters sampled from the main embankment of the Yangtze River in the Southern Jingzhou zone of China was conducted. The Kolomogorov-Smirnov test verified the statistical hypothesis that the permeability coefficient tensor has a Gaussian distribution. With the help of numerical analysis of the stochastic seepage field using the developed model, various statistical and random characteristics of the stochastic seepage field of the main embankment of the Yangtze River in the Southern Jingzhou zone of China were investigated. The model was also examined with statistical testing. Through the introduction of random variation of the upstream and downstream water levels into the model, the effects of the boundary randomness due to variation of the downstream and upstream water levels on the variation of simulated results presented with a vector series of the random seepage field were analyzed. Furthermore, the combined influence of the variation of the soil permeability coefficient and such seepage resistance measures as the cut-off wall and relief ditch on the hydraulic head distribution was analyzed and compared with the results obtained by determinate analysis. Meanwhile, sensitivities of the hydraulic gradient and downstream exit height to the variation of boundary water level were studied. The validity of the simulated results was verified by stochastic testing and measured data. The developed model provides more detail and a full stochastic algorithm to characterize and analyze three-dimensional stochastic seepage field problems.
基金Supported by Marie Curie International Incoming Fellowship (No. PIIF-GA-2009-253453)
文摘A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.
基金supported by the Foundation for Returned Scholars,the Ministry of Education of China
文摘The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.
基金the National Natural Science Foundation of China (1047202510672036)the Natural Science Foundation of Liaoning Province,China (20032109)
文摘To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters of nasal structure, three-dimensional, anatomically accurate representations of 30 adult nasal cavity models were recons- tructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated by fluid dynamics with finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to different airflow distribution in the nasal cavities and variation of the main airstream passing through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than half of the overall resistance. The characteristic model of nasal cavity was extracted on the basis of characteristic points and dimensions deduced from the original models. It showed that either the geometric structure or the airflow field of the two kinds of models was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that could properly represent the original model in model studies on nasal cavity.
文摘By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A special numerical code is introduced, in whicha lattice model and statistical approaches are used to simulate the initial heterogeneity ofmaterial properties. The size of displacement-load step is adap-tively determined so that only fewelements would fail in each load step. When the tensile principal strain in an element exceeds theultimate strain of this element, the element is considered broken and its Young's modulus is set tobe very low. Some important behaviors of heterogeneous brittle materials are indicated using thiscode. Load-displacement curves and figures of three-dimensional fracture patterns are alsonumerically obtained, which are similar to those observed in laboratory tests.
文摘This present paper has proved the theorem of the Point Optimal Variable Successive Over Relaxation (OVSOR) method of the three-dimensional unsteady flow in the reservoir, and has put forward a formu- la for calculating optimal parameters for OVSOR which vary with space points and time points. Using this method, internal memory of computer is the smallest, calculating work is the smallest, and calculating funds are the smallest. It is very easy to operate on microcomputers for three-dimensional res- ervoir simulation. The method is stable and convergent even if the time steps are taken to be large (for example, one year). The same applies for space steps. It is applicable both for homogeneous, isotropic porous mediums and for heterogeneous, anisotropic porous medium. On IBM microcomputers with internal memory of 512 thousand bytes, 8000 grid points may be cal- culated for three-dimensional simulation. It takes only two minutes to get convergence for one time step. It may be extended to three-dimensional heat conduction equation and three-dimensional simulation of the ground water flow. It looks much more advantageous for two-dimensional simulation.
文摘We simulated three-dimensional heat transfer inside a horizontal rotating drum using the discrete element method and a thermal conduction model.The aim was to determine the effect of end-wall heating on thermal behavior of a granular bed.The simulation showed that the end-wall heating significantly affects the axial temperature profile of the bed,particularly when the length-to-diameter ratio is low.Particles near the wall heated faster and became more thermally uniform than those in the center of the drum.The region affected by the end heating gradually increased over time.Increasing the rotation speed enhanced the heat conduction rate,and increasing the fill level reduced the mean temperature and thermal uniformity of the granular bed.Heat transfer was also simulated for drums with different length-to-diameter ratios.
基金Supported by the National Natural Science Foundation of China (Grant No. 50506027)
文摘Direct numerical simulation of coherent structures in the three-dimensional transitional jet with a moderate Reynolds number of 5000 was conducted. The finite volume method was used to discretize the governing equations in space; the low-storage, three-order Runge-Kutta scheme was used for time integration. The comparisons between the statistical results of the flow field; the related experimental data were performed to validate the reliability of the present numerical schemes. The emphasis was placed on the study of the spatial evolution of the three-dimensional coherent vortex structures as well as their interactions. It is found that the evolution of the spanwise vortex structures in three-dimensional space is similar to that in two-dimensional jet. The spanwise vortex structures are subject to three-dimensional instability; induce the formation of the streamwise; lateral vortex structures. Going with the breakup; mixing of the spanwise vortex structures, the streamwise; transverse vortex tubes also fall to pieces; the mixing arranged small-scale structures are formed in the flow field. Finally, the arrangement relationship among the spanwise, the streamwise; the lateral vortex structures was analyzed; their interactions were also discussed.
文摘To study the airflow distribution in human nasal cavity during respiration and the characteristic parameters for nasal structure, thirty three-dimensional, anatomically accurate representations of adult nasal cavity models were reconstructed based on processed tomography images collected from normal people. The airflow fields in nasal cavities were simulated using the fluid dynamics with the finite element software ANSYS. The results showed that the difference of human nasal cavity structure led to varying airflow distribution in the nasal cavities and the main airflow passed through the common nasal meatus. The nasal resistance in the regions of nasal valve and nasal vestibule accounted for more than a half of overall resistance. The characteristic model of nasal cavity was extracted based on the characteristic points and dimensions deducted from the original models. It showed that either the geometric structure or the air-flow field of the two kinds of model was similar. The characteristic dimensions were the characteristic parameters of nasal cavity that properly represented the original model in research for nasal cavity.
