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.展开更多
A novel two-dimensional （2D） simulation method of positive corona current pulses is proposed. A control-volume- based finite element method （CV-FEM） is used to solve continuity equations, and the Galerkin finite e...A novel two-dimensional （2D） simulation method of positive corona current pulses is proposed. A control-volume- based finite element method （CV-FEM） is used to solve continuity equations, and the Galerkin finite element method （FEM） is used to solve Poisson＇s equation. In the proposed method, photoionization is considered by adopting an exact Helmholtz photoionization model. Furthermore, fully implicit discretization and variable time step are used to ensure the time-efficiency of the present method. Finally, the method is applied to a positive rod-plane corona problem. The numerical results are in agreement with the experimental results, and the validity of the proposed method is verified.展开更多
In the research, changes of apple chemistry, and molecule, under stresses, are n terms of morphology, physiology, bio- illustrated and research and identifica- tion methods of apple resistance are explored involving ...In the research, changes of apple chemistry, and molecule, under stresses, are n terms of morphology, physiology, bio- illustrated and research and identifica- tion methods of apple resistance are explored involving drought-resistance, flood-re- sistance, salt-stress resistance, cold-hardiness and heat-resistance. In addition prospects of apple resistance research are proposed, as well.展开更多
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.展开更多
In this paper, we establish travelling wave solutions for some nonlinear evolution equations. The first integral method is used to construct the travelling wave solutions of the modified Benjamin-Bona-Mahony and the c...In this paper, we establish travelling wave solutions for some nonlinear evolution equations. The first integral method is used to construct the travelling wave solutions of the modified Benjamin-Bona-Mahony and the coupled Klein-Gordon equations. The obtained results include periodic and solitary wave solutions. The first integral method presents a wider applicability to handling nonlinear wave equations.展开更多
In this paper, we use the fractional complex transform and the （G＇/G）-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is prop...In this paper, we use the fractional complex transform and the （G＇/G）-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is proposed to convert a partial fractional differential equation with Jumarie＇s modified Riemann-Liouville derivative into its ordinary differential equation. It is shown that the considered transform and method are very efficient and powerful in solving wide classes of nonlinear fractional order equations.展开更多
The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF me...The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.展开更多
A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variable...A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.展开更多
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.展开更多
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.展开更多
A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some resul...A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some results on the related Jacobi interpolation are established. A pseudospectral scheme is proposed for the Kuramoto-Sivashisky equation. A skew symmetric decomposition is used for dealing with the nonlinear convection term. The stability and convergence of the proposed scheme are proved. The error estimates are obtained. Numerical results show the efficiency of this approach.展开更多
This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification o...This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.展开更多
In this paper, a non-Newtonian third-grade blood in coronary and femoral arteries is simulated analytically and numerically. The blood is considered as the third- grade non-Newtonian fluid under the periodic body acce...In this paper, a non-Newtonian third-grade blood in coronary and femoral arteries is simulated analytically and numerically. The blood is considered as the third- grade non-Newtonian fluid under the periodic body acceleration motion and the pulsatile pressure gradient. The hybrid multi-step differential transformation method （Hybrid-Ms- DTM） and the Crank-Nicholson method （CNM） are used to solve the partial differential equation （PDE）, and a good agreement between them is observed in the results. The effects of the some physical parameters such as the amplitude, the lead angle, and the body acceleration frequency on the velocity and shear stress profiles are considered. The results show that increasing the amplitude, Ag, and reducing the lead angle of body acceleration, 9, make higher velocity profiles on the center line of both arteries. Also, the maximum wall shear stress increases when Ag increases.展开更多
This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singula...This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.展开更多
Efficient optimization strategy of multibody systems is developed in this paper. Aug- mented Lagrange method is used to transform constrained optimal problem into unconstrained form firstly. Then methods based on seco...Efficient optimization strategy of multibody systems is developed in this paper. Aug- mented Lagrange method is used to transform constrained optimal problem into unconstrained form firstly. Then methods based on second order sensitivity are used to solve the unconstrained problem, where the sensitivity is solved by hybrid method. Generalized-α method and generalized-α projection method for the differential-algebraic equation, which shows more efficient properties with the lager time step, are presented to get state variables and adjoint variables during the optimization procedure. Numerical results validate the accuracy and efficiency of the methods is presented.展开更多
基金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.
基金Project supported by the National Basic Research Program of China(Grant No.2011CB209402)the National Natural Science Foundation of China(Grant No.51177041)the Fundamental Research Funds for the Central Universities,China(Grant No.12QX01)
文摘A novel two-dimensional （2D） simulation method of positive corona current pulses is proposed. A control-volume- based finite element method （CV-FEM） is used to solve continuity equations, and the Galerkin finite element method （FEM） is used to solve Poisson＇s equation. In the proposed method, photoionization is considered by adopting an exact Helmholtz photoionization model. Furthermore, fully implicit discretization and variable time step are used to ensure the time-efficiency of the present method. Finally, the method is applied to a positive rod-plane corona problem. The numerical results are in agreement with the experimental results, and the validity of the proposed method is verified.
基金Supported by Shandong Provincial Natural Science Foundation in China(ZR2011CM034)~~
文摘In the research, changes of apple chemistry, and molecule, under stresses, are n terms of morphology, physiology, bio- illustrated and research and identifica- tion methods of apple resistance are explored involving drought-resistance, flood-re- sistance, salt-stress resistance, cold-hardiness and heat-resistance. In addition prospects of apple resistance research are proposed, as well.
基金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.
文摘In this paper, we establish travelling wave solutions for some nonlinear evolution equations. The first integral method is used to construct the travelling wave solutions of the modified Benjamin-Bona-Mahony and the coupled Klein-Gordon equations. The obtained results include periodic and solitary wave solutions. The first integral method presents a wider applicability to handling nonlinear wave equations.
文摘In this paper, we use the fractional complex transform and the （G＇/G）-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is proposed to convert a partial fractional differential equation with Jumarie＇s modified Riemann-Liouville derivative into its ordinary differential equation. It is shown that the considered transform and method are very efficient and powerful in solving wide classes of nonlinear fractional order equations.
文摘The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.
基金supported by the National Special Fund for Major Research Instrument Development(2011YQ140145)111 Project (B07009)+1 种基金the National Natural Science Foundation of China(11002013)Defense Industrial Technology Development Program(A2120110001 and B2120110011)
文摘A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.
文摘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 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.
文摘A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some results on the related Jacobi interpolation are established. A pseudospectral scheme is proposed for the Kuramoto-Sivashisky equation. A skew symmetric decomposition is used for dealing with the nonlinear convection term. The stability and convergence of the proposed scheme are proved. The error estimates are obtained. Numerical results show the efficiency of this approach.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10771019 and 10826107)
文摘This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics. This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional. Using this method, a rapid convergent sequence is produced which converges to the exact solutions of the problem. Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.
文摘In this paper, a non-Newtonian third-grade blood in coronary and femoral arteries is simulated analytically and numerically. The blood is considered as the third- grade non-Newtonian fluid under the periodic body acceleration motion and the pulsatile pressure gradient. The hybrid multi-step differential transformation method （Hybrid-Ms- DTM） and the Crank-Nicholson method （CNM） are used to solve the partial differential equation （PDE）, and a good agreement between them is observed in the results. The effects of the some physical parameters such as the amplitude, the lead angle, and the body acceleration frequency on the velocity and shear stress profiles are considered. The results show that increasing the amplitude, Ag, and reducing the lead angle of body acceleration, 9, make higher velocity profiles on the center line of both arteries. Also, the maximum wall shear stress increases when Ag increases.
基金supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.
基金supported by the National Natural Science Foundation of China (11002075 and 10972110)
文摘Efficient optimization strategy of multibody systems is developed in this paper. Aug- mented Lagrange method is used to transform constrained optimal problem into unconstrained form firstly. Then methods based on second order sensitivity are used to solve the unconstrained problem, where the sensitivity is solved by hybrid method. Generalized-α method and generalized-α projection method for the differential-algebraic equation, which shows more efficient properties with the lager time step, are presented to get state variables and adjoint variables during the optimization procedure. Numerical results validate the accuracy and efficiency of the methods is presented.
