The effects of interface shape on stress wave distribution and attenuation were investiga- ted using finite element method （ FEM ）. The simulation results indicate that when the stress wave propagates from SiC ceram...The effects of interface shape on stress wave distribution and attenuation were investiga- ted using finite element method （ FEM ）. The simulation results indicate that when the stress wave propagates from SiC ceramic to A1 alloy, the tensile stress decreases and the attenuation coefficient of the stress wave increases with increasing central angle of the concave interface between SiC and A1. But for the convex interface, the tensile stress increases and attenuation coefficient decreases with increasing central angle. As the stress wave propagates from A1 alloy to SiC ceramic, the atten- uation coefficient of stress wave decreases with increasing the central angle of the concave interface. For the convex interface, the attenuation coefficient increases with increasing central angle.展开更多
An Arbitrary Lagrangian-Eulerian(ALE)approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation.A two-way interaction between a nonlinear finite element method and a fin...An Arbitrary Lagrangian-Eulerian(ALE)approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation.A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished.In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact,a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces.In addition,the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method.Based on this method,fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained.Numerical results such as maximum Root Mean Square(RMS)drag,general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results.This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.展开更多
Aiming at the high angle of attack pull-up and multi-channel roll pull-up coupling problems of high maneuvering aircraft, this paper establishes the flight attitude control rate by means of unsteady flow numerical sol...Aiming at the high angle of attack pull-up and multi-channel roll pull-up coupling problems of high maneuvering aircraft, this paper establishes the flight attitude control rate by means of unsteady flow numerical solution, dynamic unstructured nested mesh assembly method and numerical solution method of flight mechanics equation. On this basis, a virtual flight simulation platform integrating pneumatics, motion and control is established. Based on this virtual flight simulation platform, F-16 aircraft is simulated by high angle of attack pull-up flight mode and multi-channel roll pull-up coupling flight mode. Finally, the influence of rudder on the yaw control channel is investigated. The results show that the numerical virtual flight simulation platform established in this paper has the ability to simulate maneuvering flight of aircraft.展开更多
Temporal mode direct numerical simulation (DNS) has been done for a supersonic turbulent boundary layer on a flat plate with Mach number 4.5. It was found that the mean flow profile, the normal-wise distribution of tu...Temporal mode direct numerical simulation (DNS) has been done for a supersonic turbulent boundary layer on a flat plate with Mach number 4.5. It was found that the mean flow profile, the normal-wise distribution of turbulent Mach number and the root mean square (RMS) of the fluctuations of various variables, as well as the Reynolds stresses, bore similarity in nature, when the turbulence reached a fully developed state. But the compressibility effect was strong and must be considered. The strong Reynolds analogy (SRA) and the Morkovin hypothesis were no longer valid. From the end of transition to the fully developed state of turbulence, it was in the transient period, for which the similarity did not hold.展开更多
The near field dynamics of buoyant reactive jets with adjacent sidewalls is investigated by time-dependent three-dimensional direct simulations. The physical problem is a fuel jet issuing vertically into an oxidant am...The near field dynamics of buoyant reactive jets with adjacent sidewalls is investigated by time-dependent three-dimensional direct simulations. The physical problem is a fuel jet issuing vertically into an oxidant ambient envi- ronment in a corner configuration with sidewall boundaries. Simulation results are presented for two cases with different jet nozzle geometries: a corner-round reactive jet and a cor- ner-square reactive jet with the same cross-sectional area on the nozzle plane. Buoyancy-induced large vortical structures evolve spatially in the flow field and transition to turbulence occurs downstream. Calculation of the mean flow properties shows that entrainment of the corner-round jet is stronger than that of the corner-square jet due to the stronger vortex deformation in the corner-round case.展开更多
We analyze the electromagnetic interaction between local surface plasmon polaritons （SPPs） and an atmospheric surface wave plasma jet （ASWPJ） in combination with our designed discharge device. Before discharge, th...We analyze the electromagnetic interaction between local surface plasmon polaritons （SPPs） and an atmospheric surface wave plasma jet （ASWPJ） in combination with our designed discharge device. Before discharge, the excitation of the SPPs and the spatial distribution of the enhanced electric field are analyzed. During discharge, the critical breakdown electric field of the gases at atmospheric gas pressure and the surface wave of the SPPs converted into electron plasma waves at resonant points are studied. After discharge, the ionization development process of the ASWPJ is simulated using a two- dimensional fluid model. Our results suggest that the local enhanced electric field of SPPs is merely the precondition of gas breakdown, and the key mechanism in maintaining the discharge development of a low-power ASWPJ is the wave-mode conversion of the local enhanced electric field at the resonant point.展开更多
Discrete element method （DEM） is used to study the factors affecting agglomeration in three-dimensional copper particle systems during solid-state sintering. A new parameter is proposed to characterize agglomeration...Discrete element method （DEM） is used to study the factors affecting agglomeration in three-dimensional copper particle systems during solid-state sintering. A new parameter is proposed to characterize agglomeration. The effects of a series of factors are studied, including particle size, size distribution, inter-particle tangential viscosity, tem- perature, initial density and initial distribution of particles on agglomeration. We find that the systems with smaller particles, broader particle size distribution, smaller viscos- ity, higher sintering temperature and smaller initial density have stronger particle agglomeration and different distribu- tions of particles induce different agglomerations. This study should be very useful for understanding the phenomenon of agglomeration and the micro-structural evolution during sin- tering and guiding sintering routes to avoid detrimental ag- glomeration.展开更多
The formation mechanism of an EFP(explosively formed projectile)using a double curvature liner under the overpressure effect generated by a regular oblique reflection was investigated in this paper.Based on the detona...The formation mechanism of an EFP(explosively formed projectile)using a double curvature liner under the overpressure effect generated by a regular oblique reflection was investigated in this paper.Based on the detonation wave propagation theory,the change of the incident angle of the detonation wave collision at different positions and the distribution area of the overpressure on the surface of the liner were calculated.Three dimensional numerical simulations of the formation process of the EFP with tail.as well as the ability to penetrate 45#steel were performed using LS-DYNA software,and the EFP ve locity,the penetration ability,and the forming were assessed via experiments and x_ray photographs.The experimental results coincides with those of the simulations.Results indicate that the collision of the detonation wave was controlled to be a regular oblique reflection acting on the liner by setting the di-mensions of the unit charge and maintai ning the pressure at the collision point region at more than 2.4 times the CJ detonation when the incident angle approached the cnitical angle.The distance from the liner midline to the boundary of the area within which the pressure ratio of the regular oblique reflection pressure to the qJ detonation pressure was greater than 2.5,2,and 15was approximately 0.66 mm,132 mm,and 3.3 mm,respectively.Itis noted that pressure gradient caused the liner to turn inside out in the middle to form the head of the EFP and close the two tails of the EFP at approximately 120μs.The penetration depth of the EFP into a 45#steel target exceeded 30 mm,and there was radial expansion between the head and tail of the EFP,increasing the penetration resistance of the EFP.Therefore,the structural size of the unit charge and the liner can be further optimized to reduce resist ance to increase the penetration ability of the EFP.展开更多
This paper is concerned with establishing a reduced-order extrapolating fi- nite volume element （FVE） format based on proper orthogonal decomposition （POD） for two-dimensional （2D） hyperbolic equations. For this...This paper is concerned with establishing a reduced-order extrapolating fi- nite volume element （FVE） format based on proper orthogonal decomposition （POD） for two-dimensional （2D） hyperbolic equations. For this purpose, a semi discrete variational format relative time and a fully discrete FVE format for the 2D hyperbolic equations are built, and a set of snapshots from the very few FVE solutions are extracted on the first very short time interval. Then, the POD basis from the snapshots is formulated, and the reduced-order POD extrapolating FVE format containing very few degrees of freedom but holding sufficiently high accuracy is built. Next, the error estimates of the reduced-order solutions and the algorithm procedure for solving the reduced-order for- mat are furnished. Finally, a numerical example is shown to confirm the correctness of theoretical conclusions. This means that the format is efficient and feasible to solve the 2D hyperbolic equations.展开更多
Single-pulse chaos are studied for a functionally graded materials rectangular plate. By means of the global perturbation method, explicit conditions for the existence of a SiZnikov-type homoclinic orbit are obtained ...Single-pulse chaos are studied for a functionally graded materials rectangular plate. By means of the global perturbation method, explicit conditions for the existence of a SiZnikov-type homoclinic orbit are obtained for this sys- tem, which suggests that chaos are likely to take place. Then, numerical simulations are given to test the analytical predic- tions. And from our analysis, when the chaotic motion oc- curs, there are a quasi-period motion in a two-dimensional subspace and chaos in another two-dimensional supplemen- tary subspace.展开更多
The purpose of this overview is to provide a concise introduction to the methodology and current advances in molecular dynamics(MD) simulations. MD simulations emerged as a powerful and popular tool to study dynamic b...The purpose of this overview is to provide a concise introduction to the methodology and current advances in molecular dynamics(MD) simulations. MD simulations emerged as a powerful and popular tool to study dynamic behavior of proteins and macromolecule complexes at the atomic resolution. This approach can extend static structural data, such as X-ray crystallography, into dynamic domains with realistic timescales(up to millisecond) and high precision, therefore becoming a veritable computational microscope. This perspective covers current advances and methodology in the simulation of protein folding and drug design as illustrated by several important published examples Overall, recent progress in the simulation field points to the direction that MD will have significant impact on molecular biology and pharmaceutical science.展开更多
In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservation...In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservations principle at the groove/land boundary in each method. The block-weight approach was implemented to deal with the non-coincidence of mesh and radial groove pattern in nu- merical method. It was observed that the numerical solutions had higher precision as mesh number exceed 70 x 70, and the relaxation iteration of differential scheme presented the fastest convergence speed when relaxation factor was close to 1.94.展开更多
基金Supported by the National Basic Research Program of China("973" Program)(613135)
文摘The effects of interface shape on stress wave distribution and attenuation were investiga- ted using finite element method （ FEM ）. The simulation results indicate that when the stress wave propagates from SiC ceramic to A1 alloy, the tensile stress decreases and the attenuation coefficient of the stress wave increases with increasing central angle of the concave interface between SiC and A1. But for the convex interface, the tensile stress increases and attenuation coefficient decreases with increasing central angle. As the stress wave propagates from A1 alloy to SiC ceramic, the atten- uation coefficient of stress wave decreases with increasing the central angle of the concave interface. For the convex interface, the attenuation coefficient increases with increasing central angle.
基金co-supported by National Nature Sciences Foundation of China(Nos 11972192,11172137)A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)。
文摘An Arbitrary Lagrangian-Eulerian(ALE)approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation.A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished.In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact,a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces.In addition,the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method.Based on this method,fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained.Numerical results such as maximum Root Mean Square(RMS)drag,general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results.This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.
文摘Aiming at the high angle of attack pull-up and multi-channel roll pull-up coupling problems of high maneuvering aircraft, this paper establishes the flight attitude control rate by means of unsteady flow numerical solution, dynamic unstructured nested mesh assembly method and numerical solution method of flight mechanics equation. On this basis, a virtual flight simulation platform integrating pneumatics, motion and control is established. Based on this virtual flight simulation platform, F-16 aircraft is simulated by high angle of attack pull-up flight mode and multi-channel roll pull-up coupling flight mode. Finally, the influence of rudder on the yaw control channel is investigated. The results show that the numerical virtual flight simulation platform established in this paper has the ability to simulate maneuvering flight of aircraft.
基金supported by the National Natural Science Foundation of China(Grant No.90205021)the Foundation for the Author of National Excellent Doctoral Dissertation of China(Grant No.200328)partially supported by Liu-Hui Center of Applied Mathematics,Nankai University and Tianjin University.
文摘Temporal mode direct numerical simulation (DNS) has been done for a supersonic turbulent boundary layer on a flat plate with Mach number 4.5. It was found that the mean flow profile, the normal-wise distribution of turbulent Mach number and the root mean square (RMS) of the fluctuations of various variables, as well as the Reynolds stresses, bore similarity in nature, when the turbulence reached a fully developed state. But the compressibility effect was strong and must be considered. The strong Reynolds analogy (SRA) and the Morkovin hypothesis were no longer valid. From the end of transition to the fully developed state of turbulence, it was in the transient period, for which the similarity did not hold.
文摘The near field dynamics of buoyant reactive jets with adjacent sidewalls is investigated by time-dependent three-dimensional direct simulations. The physical problem is a fuel jet issuing vertically into an oxidant ambient envi- ronment in a corner configuration with sidewall boundaries. Simulation results are presented for two cases with different jet nozzle geometries: a corner-round reactive jet and a cor- ner-square reactive jet with the same cross-sectional area on the nozzle plane. Buoyancy-induced large vortical structures evolve spatially in the flow field and transition to turbulence occurs downstream. Calculation of the mean flow properties shows that entrainment of the corner-round jet is stronger than that of the corner-square jet due to the stronger vortex deformation in the corner-round case.
基金Project supported by the National Natural Science Foundation of China(Grant No.11105002)the Open-end Fund of State Key Laboratory of Structural Analysis for Industrial Equipment,China(Grant No.GZ1215)+1 种基金the Natural Science Foundation for University in Anhui Province of China(Grant No.KJ2013A106)the Doctoral Scientific Research Funds of Anhui University of Science and Technology,China
文摘We analyze the electromagnetic interaction between local surface plasmon polaritons （SPPs） and an atmospheric surface wave plasma jet （ASWPJ） in combination with our designed discharge device. Before discharge, the excitation of the SPPs and the spatial distribution of the enhanced electric field are analyzed. During discharge, the critical breakdown electric field of the gases at atmospheric gas pressure and the surface wave of the SPPs converted into electron plasma waves at resonant points are studied. After discharge, the ionization development process of the ASWPJ is simulated using a two- dimensional fluid model. Our results suggest that the local enhanced electric field of SPPs is merely the precondition of gas breakdown, and the key mechanism in maintaining the discharge development of a low-power ASWPJ is the wave-mode conversion of the local enhanced electric field at the resonant point.
基金supported by the National Natural Science Foundation of China (10972220,11125211 and 11021262)973 Project(2012CB937500)
文摘Discrete element method （DEM） is used to study the factors affecting agglomeration in three-dimensional copper particle systems during solid-state sintering. A new parameter is proposed to characterize agglomeration. The effects of a series of factors are studied, including particle size, size distribution, inter-particle tangential viscosity, tem- perature, initial density and initial distribution of particles on agglomeration. We find that the systems with smaller particles, broader particle size distribution, smaller viscos- ity, higher sintering temperature and smaller initial density have stronger particle agglomeration and different distribu- tions of particles induce different agglomerations. This study should be very useful for understanding the phenomenon of agglomeration and the micro-structural evolution during sin- tering and guiding sintering routes to avoid detrimental ag- glomeration.
基金The work presented in this paper has been supported by the science foundation(YT20-01-02)of Nanjing Vocational University of Industry Technology and the National Science Foundation of China under NO.11802141.
文摘The formation mechanism of an EFP(explosively formed projectile)using a double curvature liner under the overpressure effect generated by a regular oblique reflection was investigated in this paper.Based on the detonation wave propagation theory,the change of the incident angle of the detonation wave collision at different positions and the distribution area of the overpressure on the surface of the liner were calculated.Three dimensional numerical simulations of the formation process of the EFP with tail.as well as the ability to penetrate 45#steel were performed using LS-DYNA software,and the EFP ve locity,the penetration ability,and the forming were assessed via experiments and x_ray photographs.The experimental results coincides with those of the simulations.Results indicate that the collision of the detonation wave was controlled to be a regular oblique reflection acting on the liner by setting the di-mensions of the unit charge and maintai ning the pressure at the collision point region at more than 2.4 times the CJ detonation when the incident angle approached the cnitical angle.The distance from the liner midline to the boundary of the area within which the pressure ratio of the regular oblique reflection pressure to the qJ detonation pressure was greater than 2.5,2,and 15was approximately 0.66 mm,132 mm,and 3.3 mm,respectively.Itis noted that pressure gradient caused the liner to turn inside out in the middle to form the head of the EFP and close the two tails of the EFP at approximately 120μs.The penetration depth of the EFP into a 45#steel target exceeded 30 mm,and there was radial expansion between the head and tail of the EFP,increasing the penetration resistance of the EFP.Therefore,the structural size of the unit charge and the liner can be further optimized to reduce resist ance to increase the penetration ability of the EFP.
基金Project supported by the National Natural Science Foundation of China(Nos.11271127 and11671106)
文摘This paper is concerned with establishing a reduced-order extrapolating fi- nite volume element （FVE） format based on proper orthogonal decomposition （POD） for two-dimensional （2D） hyperbolic equations. For this purpose, a semi discrete variational format relative time and a fully discrete FVE format for the 2D hyperbolic equations are built, and a set of snapshots from the very few FVE solutions are extracted on the first very short time interval. Then, the POD basis from the snapshots is formulated, and the reduced-order POD extrapolating FVE format containing very few degrees of freedom but holding sufficiently high accuracy is built. Next, the error estimates of the reduced-order solutions and the algorithm procedure for solving the reduced-order for- mat are furnished. Finally, a numerical example is shown to confirm the correctness of theoretical conclusions. This means that the format is efficient and feasible to solve the 2D hyperbolic equations.
基金supported by the National Natural Science Foundation of China(11172125,11202095 and 11201226)Natural Science Foundation of Henan,China(2009B110009,B2008-56 and 649106)
文摘Single-pulse chaos are studied for a functionally graded materials rectangular plate. By means of the global perturbation method, explicit conditions for the existence of a SiZnikov-type homoclinic orbit are obtained for this sys- tem, which suggests that chaos are likely to take place. Then, numerical simulations are given to test the analytical predic- tions. And from our analysis, when the chaotic motion oc- curs, there are a quasi-period motion in a two-dimensional subspace and chaos in another two-dimensional supplemen- tary subspace.
文摘The purpose of this overview is to provide a concise introduction to the methodology and current advances in molecular dynamics(MD) simulations. MD simulations emerged as a powerful and popular tool to study dynamic behavior of proteins and macromolecule complexes at the atomic resolution. This approach can extend static structural data, such as X-ray crystallography, into dynamic domains with realistic timescales(up to millisecond) and high precision, therefore becoming a veritable computational microscope. This perspective covers current advances and methodology in the simulation of protein folding and drug design as illustrated by several important published examples Overall, recent progress in the simulation field points to the direction that MD will have significant impact on molecular biology and pharmaceutical science.
基金Supported by the Ministerial Level Foundation(2220060029)
文摘In order to calculate the pressure distribution of radial grooved thrust bearing, analytical and numerical methods were applied respectively. Grooved region and land region were linked by u- sing the mass conservations principle at the groove/land boundary in each method. The block-weight approach was implemented to deal with the non-coincidence of mesh and radial groove pattern in nu- merical method. It was observed that the numerical solutions had higher precision as mesh number exceed 70 x 70, and the relaxation iteration of differential scheme presented the fastest convergence speed when relaxation factor was close to 1.94.
