LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation ha...LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.展开更多
The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/...The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/LS-DYNA 3D Nonlinear Dynamic Finite Element Software. In this package ANSYS is the preprocessor and LS-DYNA is the postprocessor. Numerical models in the paper to actual were l：10 and the element mesh was dissected in scanning mode utilizing the symmetry characteristics of the numerical model. Five different advance rates were studied. Parameters, such as the time required to maximum stress, the action time of the available stress, the maximum velocity of the nodes, the stress penetration time, the magnitude of the stress peak and the time duration for high stress were numerically simulated. The 2.2 m advance appeared optimum from an analysis of the simulation results. The results from numerical simulation have been validated by tests with physical models.展开更多
The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock...The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.展开更多
A numerical model based on a Continuum-based Distinct Element Method(CDEM) was used to carry out a dynamic simulation of the interaction between shield and rock strata movement in longwall mining. In Northern China, t...A numerical model based on a Continuum-based Distinct Element Method(CDEM) was used to carry out a dynamic simulation of the interaction between shield and rock strata movement in longwall mining. In Northern China, the Ordos coal field geological conditions and operational characteristics were used as a case example. The CDEM was constructed on Ordos coal field shield's operation characteristics and geological conditions. Numerical modelling was carried out to investigate the effects of different mining heights on the caving process, movement characteristics, equilibrium and stability conditions of overburden as the interaction between shield and surrounding rocks. With the numerical model, the internal factors for changes in shield resistance under different mining heights was found. The quantitative relationship between mining heights and shield resistance was also obtained by the numerical simulation.展开更多
It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numer...It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numerical simulation study method,boundary condition and governing equations are presented.It is easy to study complicated problems from simple conditions,so the program concerning boundary layer condition is compiled to solve the problem.Here the spectral method is introduced,and the results are tested by Dorod’s results.展开更多
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.展开更多
This paper presents the simulation results of the wind environment around a single high-rise building and that around two tall buildings in tandem arrangement by using the lattice Boltzmann method with an aim to under...This paper presents the simulation results of the wind environment around a single high-rise building and that around two tall buildings in tandem arrangement by using the lattice Boltzmann method with an aim to understand the ventilation issues around high-rise buildings in an urban environment.We analyzed the velocity distribution around the buildings and performed numericl simulations to reveal the formation and evolution law of the complex vortex system around the high-rise buildings.Numerical simulation results manifest a periodicity phenamenon in the process of the vortex evolution.For the case of two high-rise buildings,wind velocity in the space between the two buildings is very small,which is nearly a silent regime.Wind velocity above the front building is relatively larger and the maximum wind velocity is approximately 2.5 times the incoming wind velocity.The numerical results can be used in layout planning of high-rise residential buildings to create better environment for ventilation purpose in an urban area.展开更多
In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the sof...In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.展开更多
Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-materia...Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.展开更多
The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal min...The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal mine safety. Traditional deformation monitoring methods are mostly based on single parameter, in this paper, multiple approaches are integrated: firstly, both electric and elastic models are established,from which electric field distribution and seismic wave recording are calculated and finally, the resistivity profiles and source position information are determined using inversion methods, from which then the deformation and failure of mine floor are evaluated. According to the inversion results of both electric and seismic field signals, multiple-parameter dynamic monitoring of surrounding rock deformation in deep mine can be performed. The methodology is validated using numerical simulation results which shows that the multi-parameter dynamic monitoring methods have better results for surrounding rock deformation in deep mine monitoring than single parameter methods.展开更多
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.展开更多
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.展开更多
In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numer...In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.展开更多
Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eul...Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.展开更多
The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to ...The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to simulate the transient internal flow field and the external performance of the pump during starting and stopping periods.The terms of accelerations due to variable speeds in the flow governing equations were analyzed in a multiple reference of frame(MRF).A transient CFD simulation was performed for a typical centrifugal pump by using ANSYS-CFX with the standard k-εturbulence model.The entire simulation process was composed of four stages:start-up,normal run,shutdown and post-shutdown.The function of rotating speed with regard to time was set by CEL language directly into the impeller domain in the pre-processor of the software to conduct variable speed simulation.The variations of the flow field in the centrifugal pump were obtained from the transient simulation.The changing laws of flow rate,head and other performance parameters over time were also analyzed and summarized.展开更多
Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,compara...Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.展开更多
The neutron shielding component of ITER （International Thermonuclear Experimental Reactor） vacuum vessel is a kind of structure resembling a wall in appearance. A FE （finite element） model is set up by using ANSYS...The neutron shielding component of ITER （International Thermonuclear Experimental Reactor） vacuum vessel is a kind of structure resembling a wall in appearance. A FE （finite element） model is set up by using ANSYS code in terms of its structural features. Static analysis, thermal expansion analysis and dynamic analysis are performed. The static results show that the stress and displacement distribution are allowable, but the high stress appears in the junction between the upper and lower parts. The modal analysis indicates that the biggest deformation exists in the port area. Through modal superposition, the single-point response has been found with the lower rank frequency of the acceleration seismic response spectrum. But the deformation and the stress values are within the permissible limit. The analysis results would benefit the work in the next step and provide some reference for the implementation of the engineering plan in the future.展开更多
This article presents a modeling and simulation method for transient thermal analyses of integrated circuits(ICs)using the original and voltage-in-current(VinC)latency insertion method(LIM).LIM-based algorithms are a ...This article presents a modeling and simulation method for transient thermal analyses of integrated circuits(ICs)using the original and voltage-in-current(VinC)latency insertion method(LIM).LIM-based algorithms are a set of fast transient simulation methods that solve electrical circuits in a leapfrog updating manner without relying on large matrix operations used in conventional Simulation Program with Integrated Circuit Emphasis(SPICE)-based methods which can significantly slow down the solution process.The conversion from the thermal to electrical model is performed first by using the analogy between heat and electrical conduction.Since electrical inductance has no thermal equivalence,a modified VinC LIM formulation is presented which removes the requirement of the insertion of fictitious inductors.Numerical examples are presented,which show that the modified VinC LIM formulation outperforms the basic LIM formulation,in terms of both stability and accuracy in the transient thermal simulation of ICs.展开更多
A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed b...A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed by replacing the internal gas with an equivalent solid in the modeling,which can make it convenient to simulate the thermal-mechanical coupling effects in the solid research objects with gases in them.The applied thermal expansion coefficient,Young's modulus and Poisson's ratio of the equivalent solid material are derived.A series of tests have been conducted;and the proposed equivalent solid method to simulate the gas effects is validated.展开更多
文摘LS-DYNA program and the principle of ALE method were introduced, and the target features of the reinforced concrete penetration were analyzed by using the D material model and the ALE method. A numerical simulation has been done to show the penetration visually and veritably. The simulation results are analyzed carefully and explicitly prove their significance to the research of reinforced concrete penetration.
文摘The mechanism of stress generation and propagation by detonation loading in five separate independent advance of ore breaking patterns is discussed in the paper. An elastic numerical model was developed using AN- SYS/LS-DYNA 3D Nonlinear Dynamic Finite Element Software. In this package ANSYS is the preprocessor and LS-DYNA is the postprocessor. Numerical models in the paper to actual were l：10 and the element mesh was dissected in scanning mode utilizing the symmetry characteristics of the numerical model. Five different advance rates were studied. Parameters, such as the time required to maximum stress, the action time of the available stress, the maximum velocity of the nodes, the stress penetration time, the magnitude of the stress peak and the time duration for high stress were numerically simulated. The 2.2 m advance appeared optimum from an analysis of the simulation results. The results from numerical simulation have been validated by tests with physical models.
基金sponsored by the project No.50404001 from the National Natural Science Foundation of Chinathe National Key Fundamental Research & Development Project(Grant No.2007CB209601)+1 种基金the China National PetroleumCorporation Fundamental Research Program (Grant No.06A30102)the China Postdoctoral Science Foundation(Project No.2004035350)
文摘The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method （LBM） was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation （I-Sw）. The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.
基金funded by the National Natural Science Foundation of China (Nos. U1261207 and 51274086)
文摘A numerical model based on a Continuum-based Distinct Element Method(CDEM) was used to carry out a dynamic simulation of the interaction between shield and rock strata movement in longwall mining. In Northern China, the Ordos coal field geological conditions and operational characteristics were used as a case example. The CDEM was constructed on Ordos coal field shield's operation characteristics and geological conditions. Numerical modelling was carried out to investigate the effects of different mining heights on the caving process, movement characteristics, equilibrium and stability conditions of overburden as the interaction between shield and surrounding rocks. With the numerical model, the internal factors for changes in shield resistance under different mining heights was found. The quantitative relationship between mining heights and shield resistance was also obtained by the numerical simulation.
文摘It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numerical simulation study method,boundary condition and governing equations are presented.It is easy to study complicated problems from simple conditions,so the program concerning boundary layer condition is compiled to solve the problem.Here the spectral method is introduced,and the results are tested by Dorod’s results.
基金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.
基金Funded by the Natural Science Foundation of Tianjin Municipality (No. 06YFJMJC05300)the Science and Technology Development Foundation for Universities of Tianjin Municipality (No. 20060823)
文摘This paper presents the simulation results of the wind environment around a single high-rise building and that around two tall buildings in tandem arrangement by using the lattice Boltzmann method with an aim to understand the ventilation issues around high-rise buildings in an urban environment.We analyzed the velocity distribution around the buildings and performed numericl simulations to reveal the formation and evolution law of the complex vortex system around the high-rise buildings.Numerical simulation results manifest a periodicity phenamenon in the process of the vortex evolution.For the case of two high-rise buildings,wind velocity in the space between the two buildings is very small,which is nearly a silent regime.Wind velocity above the front building is relatively larger and the maximum wind velocity is approximately 2.5 times the incoming wind velocity.The numerical results can be used in layout planning of high-rise residential buildings to create better environment for ventilation purpose in an urban area.
文摘In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.
基金Supported by the National Natural Science Foundation of China(10476011)~~
文摘Numerical simulations are performed on the interface with large deformation induced by the interaction between a moving shock and two consecutive bubbles. The high performance of the level set method for multi-material interfaces is demonstrated. Discontinuous Galerkin finite element method is used to solve Euleri- an equations. And the fifth-order weighted essentially non-oscillatory （WENO） scheme is used to solve the level set equation for capturing multi-material interfaces. The ghost fluid method is used to deal with the interfacial boundary condition. Results are obtained for two bubble interacting with a moving shock. The contours of the constant density and the pressure at different time are given. In the computational domain, three different cases are considered, i.e. two helium bubbles, a helium bubble followed by an R22 bubble in the direction of the moving shock, and an R22 bubble followed by a helium bubble. Computational results indicate that multi-mate- rial interfaces can be properly captured by the level set method. Therefore, for problems involving the flow of three different materials with two different interfaces, each interface separating two different materials can be similarly handled.
基金financial support from the Fundamental Research Funds for the Central Universities of China (No. 2015QNB19)the financial support from the Open Fund of Key Laboratory of Safety and High-efficiency Coal Mining, Ministry of Education of China (No. JYBSYS2015107)+2 种基金the National Natural Science Foundation of China (Nos. 51404254, 41430317 and U1261202)the China Postdoctoral Science Foundation of China (No. 2014M560465)the Jiangsu Planned Projects for Postdoctoral Research Funds of China (No. 1302050B)
文摘The level of deformation development of surrounding rocks is a vital predictor to evaluate impending coal mine disasters and it is important to establish accurate measurements of the deformed status to ensure coal mine safety. Traditional deformation monitoring methods are mostly based on single parameter, in this paper, multiple approaches are integrated: firstly, both electric and elastic models are established,from which electric field distribution and seismic wave recording are calculated and finally, the resistivity profiles and source position information are determined using inversion methods, from which then the deformation and failure of mine floor are evaluated. According to the inversion results of both electric and seismic field signals, multiple-parameter dynamic monitoring of surrounding rock deformation in deep mine can be performed. The methodology is validated using numerical simulation results which shows that the multi-parameter dynamic monitoring methods have better results for surrounding rock deformation in deep mine monitoring than single parameter methods.
基金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.
基金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.
文摘In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.
文摘Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.
文摘The instantaneous variations of the hydraulic characteristics take place in centrifugal pumps during their start-up,shutdown and other variable speed operations.In this paper,the variable speed method was proposed to simulate the transient internal flow field and the external performance of the pump during starting and stopping periods.The terms of accelerations due to variable speeds in the flow governing equations were analyzed in a multiple reference of frame(MRF).A transient CFD simulation was performed for a typical centrifugal pump by using ANSYS-CFX with the standard k-εturbulence model.The entire simulation process was composed of four stages:start-up,normal run,shutdown and post-shutdown.The function of rotating speed with regard to time was set by CEL language directly into the impeller domain in the pre-processor of the software to conduct variable speed simulation.The variations of the flow field in the centrifugal pump were obtained from the transient simulation.The changing laws of flow rate,head and other performance parameters over time were also analyzed and summarized.
基金National Natural Science Foundation of China(No.51674280,51774308,51704033,51722406,51950410591)Shandong Provincial Natural Science Foundation(ZR2019JQ21,JQ201808)+3 种基金the Fundamental Research Funds for the Central Universities(No.20CX02113A)National Science and Technology Major Project(2016ZX05014-000407)Program for Changjiang Scholars and Innovative Research Team in University(IRT_16R69)PetroChina Innovation Foundation(No.2018D-5007-0210)。
文摘Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.
基金the National 973 program of China(No.2004CB720704)
文摘The neutron shielding component of ITER （International Thermonuclear Experimental Reactor） vacuum vessel is a kind of structure resembling a wall in appearance. A FE （finite element） model is set up by using ANSYS code in terms of its structural features. Static analysis, thermal expansion analysis and dynamic analysis are performed. The static results show that the stress and displacement distribution are allowable, but the high stress appears in the junction between the upper and lower parts. The modal analysis indicates that the biggest deformation exists in the port area. Through modal superposition, the single-point response has been found with the lower rank frequency of the acceleration seismic response spectrum. But the deformation and the stress values are within the permissible limit. The analysis results would benefit the work in the next step and provide some reference for the implementation of the engineering plan in the future.
基金This work was supported by the Fundamental Research Grant Scheme(FRGS)sponsored by the Ministry of Higher Education,Malaysia under Grant No.FRGS/1/2020/TK0/USM/02/7.
文摘This article presents a modeling and simulation method for transient thermal analyses of integrated circuits(ICs)using the original and voltage-in-current(VinC)latency insertion method(LIM).LIM-based algorithms are a set of fast transient simulation methods that solve electrical circuits in a leapfrog updating manner without relying on large matrix operations used in conventional Simulation Program with Integrated Circuit Emphasis(SPICE)-based methods which can significantly slow down the solution process.The conversion from the thermal to electrical model is performed first by using the analogy between heat and electrical conduction.Since electrical inductance has no thermal equivalence,a modified VinC LIM formulation is presented which removes the requirement of the insertion of fictitious inductors.Numerical examples are presented,which show that the modified VinC LIM formulation outperforms the basic LIM formulation,in terms of both stability and accuracy in the transient thermal simulation of ICs.
基金Supported by Natural Science Foundation of China(Nos.10772049 and 11072062)Research and Development Program of China(863 Program, No.2009AA04Z408)+1 种基金Natural Science Foundation of Shanghai(No.06ZR14009)Pujiang Scholar Program and the Wangdao Scholar Program(No.08076) of Fudan University
文摘A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed by replacing the internal gas with an equivalent solid in the modeling,which can make it convenient to simulate the thermal-mechanical coupling effects in the solid research objects with gases in them.The applied thermal expansion coefficient,Young's modulus and Poisson's ratio of the equivalent solid material are derived.A series of tests have been conducted;and the proposed equivalent solid method to simulate the gas effects is validated.
