The Anjialing No. 1 Coal Mine in Shanxi Province, China, contains a complicated old goaf and an unknown water distribution that hold high potential for serious water hazards. Due to poor detection resolution, previous...The Anjialing No. 1 Coal Mine in Shanxi Province, China, contains a complicated old goaf and an unknown water distribution that hold high potential for serious water hazards. Due to poor detection resolution, previous attempts have failed to determine the scope of the old goal and the water distribution in the mine by separate use of various exploration methods such as seismic method, direct current resistivity, audio magnetotellurics, controlled-source audio-frequency magnetotellurics, and transient electromag-netics. To solve this difficult problem, a combination of the wide-field electromagnetic method and the flow field fitting method with three-dimensional resistivity data inversion was applied to determine the precise scope of the goal and the locations where water is present, and to identify the hydraulic con- nection between the water layers so as to provide reliable technical support for safe coal production. Reasonable results were achieved, with all these goals being met. As a result, a mining area of nearly 4 km^2 has been released for operation.展开更多
Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with...Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.展开更多
Fluid flow can significantly change the evolution of microstructural morphology. However, relatively little is known how the fluid flow, concentration and microstructure affect each other quantitatively, which is esse...Fluid flow can significantly change the evolution of microstructural morphology. However, relatively little is known how the fluid flow, concentration and microstructure affect each other quantitatively, which is essential to optimize processing parameters. A quantitative simulation study of Al-Cu solidified equiaxed dendrite evolution under forced flow based on phase field-lattice Boltzmann method(PF-LBM) is carried out. Results obtained are validated by Gibbs-Thomson relation at the dendrite tip. Compared with the equiaxed dendrite evolution without flow, the upstream dendrite arm is enhanced while the downstream arm is inhibited. Besides, as the inlet flow rate increases, the secondary arms attached onto the upstream primary arm and the upstream side of the primary arm normal to the inflow velocity has been well developed. Results show that sidewise instabilities of the primary dendrite arm and onset of secondary arm is caused by the local concentration perturbation and will be enhanced or inhabited by the flow. It is believed that the coupled PF-LBM method is able to handle dendrite evolution under forced flow quantitatively, which helps in investigating the solidified dendrite morphology evolution.展开更多
In this paper, on the basis of experimental data of two kinds of chemical explosions, the piston-pushing model of spherical blast-waves and the second-order Godunov-type scheme of finite difference methods with high i...In this paper, on the basis of experimental data of two kinds of chemical explosions, the piston-pushing model of spherical blast-waves and the second-order Godunov-type scheme of finite difference methods with high identification to discontinuity are used to the numerical reconstruction of part of an actual hemispherical blast-wave flow field by properly adjusting the moving bounary conditions of a piston. This method is simple and reliable. It is suitable to the evaluation of effects of the blast-wave flow field away from the explosion center.展开更多
The in-cylinder flow field of the internal combustion engine is an important factor affecting the quality and combustion quality of the fuel mixture in the cylinder. In order to calculate the high-precision flow field...The in-cylinder flow field of the internal combustion engine is an important factor affecting the quality and combustion quality of the fuel mixture in the cylinder. In order to calculate the high-precision flow field, the paper presents a flow field calculation method based on the optical flow algorithm. The motion of the point was calculated using the change in pixel intensity within two temporally adjacent frame images. The results show the high accuracy and resolution of the flow field at small displacement conditions.展开更多
In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Rey...In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k - ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propagates over breakwater. When wave crest propagates over breakwater, the anticlockwise vortex may generate. On the contrary, when wave hollow propagates over breakwater, the clockwise vortex may generate. Meanwhile, the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy, turbulent dissipation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.展开更多
The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at th...The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.展开更多
An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient....An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient. The cross-correlation algorithm with sub-pixel accuracy is introduced and used to calculate the background-element displacement of a disturbed image and a fourth-order difference scheme is also developed to solve the Poisson equation. An experiment for a disturbed flow field caused by a burning candle is performed to validate the built BOS system and the results indicate that density and temperature distribution of the disturbed flow field can be reconstructed accurately. A notable conclusion is that in order to make the reconstructed results have a satisfactory accuracy, the inquiry step length should be less than the size of the interrogation window.展开更多
Submerged horizontal plate can be considered as a new concept breakwater. In order to reveal the wave elimination mechanism of this type breakwater, boundary element method is utilized to investigate the velocity fiel...Submerged horizontal plate can be considered as a new concept breakwater. In order to reveal the wave elimination mechanism of this type breakwater, boundary element method is utilized to investigate the velocity field around plate carefully. The flow field analysis shows that the interaction between incident wave and reverse flow caused by submerged plate will lead to the formation of wave elimination area around both sides of the plate. The velocity magnitude of flow field has been reduced and this is the main reason of wave elimination.展开更多
Numerical simulation based on phase field method is performed to describe solidification process of pure material in a free or forced flow. The evolution of the interface is showed, and the effects of mesh grid and fl...Numerical simulation based on phase field method is performed to describe solidification process of pure material in a free or forced flow. The evolution of the interface is showed, and the effects of mesh grid and flow velocity on succinonitrite shape are studied. These results indicate that crystal grows into an equiaxial dendrite in a free flow and into an asymmetrical dendritic in a forced flow. With increasing flow velocity, the upstream dendritic arm tip grows faster and the downstream arm grows slower. However, the evolution of the perpendicular tip has no significant change. In addition, mesh grid has no influence on dendritic growth shape when mesh grid is above 300×300.展开更多
The phase-field model coupled with a flow field was used to simulate the solidification of pure materials by the finite difference method. The effects of initial crystal radius, the space step and the interface thickn...The phase-field model coupled with a flow field was used to simulate the solidification of pure materials by the finite difference method. The effects of initial crystal radius, the space step and the interface thickness on the dendrite growth were studied. Results indicate that the grain grows into an equiaxial dendrite during free flow and into a typical branched structure under forced flow. The radius of an initial crystal can affect the growth of side-branches but not the stability of the dendrite＇s tip when an appropriate value is assigned to it. With an increase in space steps, side-branches appear at the upstream of the longitudinal principal branch and they grow rapidly. With an increase in the interface thickness, the trunk of the longitudinal upstream and lateral principal branches grow longer and become more slender while the number of secondary branches increases.展开更多
In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum alo...In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.展开更多
Based on an analogy among the governing equations among the electrostatic field,static magnetic field and fluid flow,all these equations can be expressed by a general differential equation,and the control volume metho...Based on an analogy among the governing equations among the electrostatic field,static magnetic field and fluid flow,all these equations can be expressed by a general differential equation,and the control volume method,which has been widely used in computational fluid dynamics,was developed to solve all these equations.Numerical results showed that control volume method is a successful numerical in the computational fluid dynamics and the electromagnetics,and the upwind(or downwind)scheme have the clear mathematical meaning and can be applied to calculate the charges' movement in the electrostatic field.展开更多
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.展开更多
文摘The Anjialing No. 1 Coal Mine in Shanxi Province, China, contains a complicated old goaf and an unknown water distribution that hold high potential for serious water hazards. Due to poor detection resolution, previous attempts have failed to determine the scope of the old goal and the water distribution in the mine by separate use of various exploration methods such as seismic method, direct current resistivity, audio magnetotellurics, controlled-source audio-frequency magnetotellurics, and transient electromag-netics. To solve this difficult problem, a combination of the wide-field electromagnetic method and the flow field fitting method with three-dimensional resistivity data inversion was applied to determine the precise scope of the goal and the locations where water is present, and to identify the hydraulic con- nection between the water layers so as to provide reliable technical support for safe coal production. Reasonable results were achieved, with all these goals being met. As a result, a mining area of nearly 4 km^2 has been released for operation.
文摘Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.
基金financially supported by the National Key Research and Development Program of China(Grant No.2017YFB1103700)the National Science Foundation of China(Grant Nos.51320105003 and 51674153)
文摘Fluid flow can significantly change the evolution of microstructural morphology. However, relatively little is known how the fluid flow, concentration and microstructure affect each other quantitatively, which is essential to optimize processing parameters. A quantitative simulation study of Al-Cu solidified equiaxed dendrite evolution under forced flow based on phase field-lattice Boltzmann method(PF-LBM) is carried out. Results obtained are validated by Gibbs-Thomson relation at the dendrite tip. Compared with the equiaxed dendrite evolution without flow, the upstream dendrite arm is enhanced while the downstream arm is inhibited. Besides, as the inlet flow rate increases, the secondary arms attached onto the upstream primary arm and the upstream side of the primary arm normal to the inflow velocity has been well developed. Results show that sidewise instabilities of the primary dendrite arm and onset of secondary arm is caused by the local concentration perturbation and will be enhanced or inhabited by the flow. It is believed that the coupled PF-LBM method is able to handle dendrite evolution under forced flow quantitatively, which helps in investigating the solidified dendrite morphology evolution.
文摘In this paper, on the basis of experimental data of two kinds of chemical explosions, the piston-pushing model of spherical blast-waves and the second-order Godunov-type scheme of finite difference methods with high identification to discontinuity are used to the numerical reconstruction of part of an actual hemispherical blast-wave flow field by properly adjusting the moving bounary conditions of a piston. This method is simple and reliable. It is suitable to the evaluation of effects of the blast-wave flow field away from the explosion center.
文摘The in-cylinder flow field of the internal combustion engine is an important factor affecting the quality and combustion quality of the fuel mixture in the cylinder. In order to calculate the high-precision flow field, the paper presents a flow field calculation method based on the optical flow algorithm. The motion of the point was calculated using the change in pixel intensity within two temporally adjacent frame images. The results show the high accuracy and resolution of the flow field at small displacement conditions.
基金The National Natural Science Foundation of China under contract Nos 50979008 and 50909009Program for Hunan Province Key Laboratory of WaterSediment Sciences & Flood Hazard Prevention and Open Research Fund Program of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University No.2008490911
文摘In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k - ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propagates over breakwater. When wave crest propagates over breakwater, the anticlockwise vortex may generate. On the contrary, when wave hollow propagates over breakwater, the clockwise vortex may generate. Meanwhile, the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy, turbulent dissipation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.
文摘The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.
基金supported by the Key Program of the National Natural Science Foundation of China(Grant No.NSFC 91441205)
文摘An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient. The cross-correlation algorithm with sub-pixel accuracy is introduced and used to calculate the background-element displacement of a disturbed image and a fourth-order difference scheme is also developed to solve the Poisson equation. An experiment for a disturbed flow field caused by a burning candle is performed to validate the built BOS system and the results indicate that density and temperature distribution of the disturbed flow field can be reconstructed accurately. A notable conclusion is that in order to make the reconstructed results have a satisfactory accuracy, the inquiry step length should be less than the size of the interrogation window.
基金supported by the Fundamental Research Funds for the Dalian University of Technology(Grant No.DUT10LK43)the National Key Basic Research Program of China(Grant No.2013CB036101)
文摘Submerged horizontal plate can be considered as a new concept breakwater. In order to reveal the wave elimination mechanism of this type breakwater, boundary element method is utilized to investigate the velocity field around plate carefully. The flow field analysis shows that the interaction between incident wave and reverse flow caused by submerged plate will lead to the formation of wave elimination area around both sides of the plate. The velocity magnitude of flow field has been reduced and this is the main reason of wave elimination.
基金supported by the Major Program of the National Natural Science Foundation of China (Nos.50331040 and 60171034).
文摘Numerical simulation based on phase field method is performed to describe solidification process of pure material in a free or forced flow. The evolution of the interface is showed, and the effects of mesh grid and flow velocity on succinonitrite shape are studied. These results indicate that crystal grows into an equiaxial dendrite in a free flow and into an asymmetrical dendritic in a forced flow. With increasing flow velocity, the upstream dendritic arm tip grows faster and the downstream arm grows slower. However, the evolution of the perpendicular tip has no significant change. In addition, mesh grid has no influence on dendritic growth shape when mesh grid is above 300×300.
基金supported by a Program of the Lanzhou Science and Technology Bureau (No.2009-1-9)the Doctoral Fund of Lanzhou University of Technology(No.SB01200606)
文摘The phase-field model coupled with a flow field was used to simulate the solidification of pure materials by the finite difference method. The effects of initial crystal radius, the space step and the interface thickness on the dendrite growth were studied. Results indicate that the grain grows into an equiaxial dendrite during free flow and into a typical branched structure under forced flow. The radius of an initial crystal can affect the growth of side-branches but not the stability of the dendrite＇s tip when an appropriate value is assigned to it. With an increase in space steps, side-branches appear at the upstream of the longitudinal principal branch and they grow rapidly. With an increase in the interface thickness, the trunk of the longitudinal upstream and lateral principal branches grow longer and become more slender while the number of secondary branches increases.
文摘In this paper, unsteady free convection heat transfer flow over a vertical plate in the presence of a magnetic field is discussed in detail. The dimensionless partial differential equations of continuity, momentum along energy are analyzed with suitable transformations. For numerical calculation, an implicit finite difference method is applied to solve a set of nonlinear dimensionless partial differential equations. Dimensionless velocity and temperature profile are also investigated due to the effects of assumed parameters in the concerned problem. An explicit finite difference technique is used to compute velocity and temperature profiles. The stability conditions are also examined.
基金Item Sponsored by 111 project (B07015) the National Natural Science Foundation of China and Shanghai Baosteel (No.50834010) +1 种基金the Fundamental Research Funds for the Central Universities (N100409007) the Doctor Startup Foundation of Liaoning Province (20111009)
文摘Based on an analogy among the governing equations among the electrostatic field,static magnetic field and fluid flow,all these equations can be expressed by a general differential equation,and the control volume method,which has been widely used in computational fluid dynamics,was developed to solve all these equations.Numerical results showed that control volume method is a successful numerical in the computational fluid dynamics and the electromagnetics,and the upwind(or downwind)scheme have the clear mathematical meaning and can be applied to calculate the charges' movement in the electrostatic field.
文摘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.
