Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assu...Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assumptions and approximations. In this kind of groundwater systems, studying contamination spreading and analyzing risk are challenging tasks. The main difficulty stems from determining both the travel times and the maximum distances covered by pollutants. In this context, the risk of contamination in the deep carbonate aquifer of the Augusta coastal area is presented. We used a geostatistical approach and numerical codes (MODFLOW-2000, MT3DMS) to reconstruct the complex geological framework of the study area where several contamination scenarios of hypothetical point source in a risk assessment framework were simulated. Results on the contaminant spreading are discussed and the effect of the geological structures characterizing the zone under study, namely horst and graben, are described.展开更多
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.展开更多
In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) nu...In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) numerical simulation,the mean and fluctuating forces,Strouhal number(St) and vortex shedding pattern in the wake for each cylinder were analyzed with the spacing ratio(L /D) ranging from 1.5 to 6.0.It was found that,four different vortex modes(viz.,flip-flopping,shielding anti-phase-synchronized,in-phasesynchronized and anti-phase-synchronized) gradually appear with the increase of the L/D ratio.The average drag coefficient of the upstream cylinders is larger than that of the downstream cylinders,while the downstream cylinders usually undergo serious fluctuating forces.When the L/D ratio ranges from 3.0 to 4.0,the dominant frequency of the drag coefficient is equal to the value of St of upstream cylinders.This indicates that a simultaneous resonance in the in-flow and cross-flow directions may occur for some single structures of a multi-body oscillating system.For the 3-D numerical simulation,the L/D and aspect ratios are kept constant as 5.0 and 10,respectively.It was found that some vortices are formed in the wake of the upstream cylinders.Besides,with the same spacing ratio,the calculated drag coefficient and lift coefficient fluctuation are slightly larger than the 2-D results,but with a phase difference.展开更多
Numerical simulation of three-dimensional flow field and film cooling effectiveness in film-cooled turbine rotor and stationary turbine cascade were carried out by using the k- ε turbulence model, and the predictions...Numerical simulation of three-dimensional flow field and film cooling effectiveness in film-cooled turbine rotor and stationary turbine cascade were carried out by using the k- ε turbulence model, and the predictions of the three-dimensional velocities were compared with the measured results by Laser-Doppler Velocimetry （LDV）. Results reveal the secondary flow near the blade surface in the wake region behind the jet hole. Compared with the stationary cascade, there are the centrifugal force and Coriolis force existing in the flow field of the turbine rotor, and these forces make the three-dimensional flow field change in the turbine rotor, especially for the radial velocity. The effect of rotation on the flow field and the film cooling effectiveness on the pressure side is more apparent than that on the suction side as is shown in the computational and measured results, and the low film cooling effectiveness appears on the pressure surface of the turbine rotor blade compared with that of the stationary cascade.展开更多
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.展开更多
文摘Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assumptions and approximations. In this kind of groundwater systems, studying contamination spreading and analyzing risk are challenging tasks. The main difficulty stems from determining both the travel times and the maximum distances covered by pollutants. In this context, the risk of contamination in the deep carbonate aquifer of the Augusta coastal area is presented. We used a geostatistical approach and numerical codes (MODFLOW-2000, MT3DMS) to reconstruct the complex geological framework of the study area where several contamination scenarios of hypothetical point source in a risk assessment framework were simulated. Results on the contaminant spreading are discussed and the effect of the geological structures characterizing the zone under study, namely horst and graben, are described.
基金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.
基金Sponsored by the National Natural Science Foundation of China(Grant No.50538020)
文摘In the present paper,two-and three-dimensional numerical simulations of the flow interference between four cylinders in an in-line square arrangement at Re = 200 are performed.Assisted with the two-dimensional(2-D) numerical simulation,the mean and fluctuating forces,Strouhal number(St) and vortex shedding pattern in the wake for each cylinder were analyzed with the spacing ratio(L /D) ranging from 1.5 to 6.0.It was found that,four different vortex modes(viz.,flip-flopping,shielding anti-phase-synchronized,in-phasesynchronized and anti-phase-synchronized) gradually appear with the increase of the L/D ratio.The average drag coefficient of the upstream cylinders is larger than that of the downstream cylinders,while the downstream cylinders usually undergo serious fluctuating forces.When the L/D ratio ranges from 3.0 to 4.0,the dominant frequency of the drag coefficient is equal to the value of St of upstream cylinders.This indicates that a simultaneous resonance in the in-flow and cross-flow directions may occur for some single structures of a multi-body oscillating system.For the 3-D numerical simulation,the L/D and aspect ratios are kept constant as 5.0 and 10,respectively.It was found that some vortices are formed in the wake of the upstream cylinders.Besides,with the same spacing ratio,the calculated drag coefficient and lift coefficient fluctuation are slightly larger than the 2-D results,but with a phase difference.
基金the National Natural Science Foundation of China (Grant No. 50406017).
文摘Numerical simulation of three-dimensional flow field and film cooling effectiveness in film-cooled turbine rotor and stationary turbine cascade were carried out by using the k- ε turbulence model, and the predictions of the three-dimensional velocities were compared with the measured results by Laser-Doppler Velocimetry （LDV）. Results reveal the secondary flow near the blade surface in the wake region behind the jet hole. Compared with the stationary cascade, there are the centrifugal force and Coriolis force existing in the flow field of the turbine rotor, and these forces make the three-dimensional flow field change in the turbine rotor, especially for the radial velocity. The effect of rotation on the flow field and the film cooling effectiveness on the pressure side is more apparent than that on the suction side as is shown in the computational and measured results, and the low film cooling effectiveness appears on the pressure surface of the turbine rotor blade compared with that of the stationary cascade.
文摘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.
