We propose the meshfree-based physics-informed neural networks for solving the unsteady Oseen equations.Firstly,based on the ideas of meshfree and small sample learning,we only randomly select a small number of spatio...We propose the meshfree-based physics-informed neural networks for solving the unsteady Oseen equations.Firstly,based on the ideas of meshfree and small sample learning,we only randomly select a small number of spatiotemporal points to train the neural network instead of forming a mesh.Specifically,we optimize the neural network by minimizing the loss function to satisfy the differential operators,initial condition and boundary condition.Then,we prove the convergence of the loss function and the convergence of the neural network.In addition,the feasibility and effectiveness of the method are verified by the results of numerical experiments,and the theoretical derivation is verified by the relative error between the neural network solution and the analytical solution.展开更多
Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be ...Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be satisfied exactly without increasing the computational time. The novel scheme is useful for computing transonic, strong disturbance, unsteady flows with high reduced frequencies. The scheme makes use of curvefitted orthogonal meshes and the lattice control technique to obtain the optimal grid distribution. The numerical results are satisfactory.展开更多
This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse th...This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse the behaviour of axial velocity and shear stress for unsteady flow of nun-Newtonian visco-elastie fluids in tube, and to provide a theoretical base for the projection of pipe-line engineering.展开更多
The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for...The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.展开更多
In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split i...In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.展开更多
The forward flight of a model butterfly was stud- ied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational f...The forward flight of a model butterfly was stud- ied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational forces, where the aerodynamic forces were generated by flapping wings which moved with the body, allowing the body os- cillations of the model butterfly to be simulated. The main results are as follows: (1) The aerodynamic force produced by the wings is approximately perpendicular to the long-axis of body and is much larger in the downstroke than in the up- stroke. In the downstroke the body pitch angle is small and the large aerodynamic force points up and slightly backward, giving the weight-supporting vertical force and a small neg- ative horizontal force, whilst in the upstroke, the body an- gle is large and the relatively small aerodynamic force points forward and slightly downward, giving a positive horizon- tal force which overcomes the body drag and the negative horizontal force generated in the downstroke. (2) Pitching oscillation of the butterfly body plays an equivalent role of the wing-rotation of many other insects. (3) The body-mass- specific power of the model butterfly is 33.3 W/kg, not very different from that of many other insects, e.g., fruitflies and dragonflies.展开更多
The aim of this paper is to study the unsteady korteweg-de vries equation that plays an important role in describing the shallow water.Two analytical techniques namely the Sardar-subequation method and the energy bala...The aim of this paper is to study the unsteady korteweg-de vries equation that plays an important role in describing the shallow water.Two analytical techniques namely the Sardar-subequation method and the energy balance method are employed to seek the abundant traveling wave solutions for the first time.By these two methods,plenty of traveling wave solutions such as the bright solitary wave solutions,dark solitary wave solutions,singular periodic wave solutions and perfect periodic wave solution that expressed in terms of the generalized hyperbolic functions,generalized trigonometric functions and the cosine function are obtained.Finally,the dynamic behaviors of the solutions are described through the 3D plot and 2D curve.The results in this paper demonstrate that the proposed methods are powerful and effective to construct the traveling wave solutions of the nonlinear evolution equations in ocean engineering and science.展开更多
Stability related to theoretical model for catastrophic weather prediction, which includes non-hydrostatic perfect elastic model and anelastic model, is discussed and analyzed in detail. It is proved that non-hydrosta...Stability related to theoretical model for catastrophic weather prediction, which includes non-hydrostatic perfect elastic model and anelastic model, is discussed and analyzed in detail. It is proved that non-hydrostatic perfect elastic equations set is stable in the class of infinitely differentiable function. However, for the anelastic equations set, its continuity equation is changed in form because of the particular hypothesis for fluid, so "the matching consisting of both viscosity coefficient and incompressible assumption" appears, thereby the most important equations set of this class in practical prediction shows the same instability in topological property as Navier-Stokes equation, which should be avoided first in practical numerical prediction. In light of this, the referenced suggestions to amend the applied model are finally presented.展开更多
In this paper, the Crank-Nicholson + component-consistent pressure correction method for the numerical solution of the unsteady incompressible Navier-Stokes equation of [1] on the rectangular half-Staggered mesh has b...In this paper, the Crank-Nicholson + component-consistent pressure correction method for the numerical solution of the unsteady incompressible Navier-Stokes equation of [1] on the rectangular half-Staggered mesh has been extended to the curvilinear half-Staggered mesh. The discrete projection, both for the projection step in the solution procedure and for the related differential-algebraic equations, has been carefully studied and verified. It is proved that the proposed method is also unconditionally (in t) nonlinearly stable on the curvilinear mesh, provided the mesh is not too skewed. It is seen that for problems with an outflow boundary, the half-Staggered mesh is especially advantageous. Results of preliminary numerical experiments support these claims.展开更多
Fourier-Legendre spectral approximation for the unsteady Navier-Stokes equations is analyzed. The generalized stability and convergence are proved respectively.
The numerical methods for computing the stability derivatives of the aircraft by solving unsteady sensitivity equations which was proposed in our previous papers was extended to solve three-dimensional problems in thi...The numerical methods for computing the stability derivatives of the aircraft by solving unsteady sensitivity equations which was proposed in our previous papers was extended to solve three-dimensional problems in this paper.Both the static and dynamic derivatives of the hypersonic blunt cone undergoing pitching oscillation around a fixed point were computed using the new methods.The predicted static derivative and dynamic derivative were found to be in reasonable agreement with the experimental data.For the present method,it is possible to distinguish the components of dynamic derivatives caused by different state parameters.It is found that C_(m_α) and C_(mq) are usually of opposite signs and tend to eliminate each other,which makes C_(m_α)+C_(mq) being much smaller than its individual components.Another feature of this method is that the moment of pressure derivatives proposed in the present paper can be used to predict the contribution of each part of the blunt cone to the overall stability quantitatively.It is found that the head region is crucial for the static stability and the body region contributes most to the dynamic stability.展开更多
To find intrinsically different symmetry reductions and inequivalent group invariant solutions of the 2D unsteady incompressible boundary-layer equations, a two-dimensional optimal system is constructed which attribut...To find intrinsically different symmetry reductions and inequivalent group invariant solutions of the 2D unsteady incompressible boundary-layer equations, a two-dimensional optimal system is constructed which attributed to the classification of the corresponding Lie subalgebras. The comprehensiveness and inequivalence of the optimal system are shown clearly under different values of invariants. Then by virtue of the optimal system obtained, the boundary-layer equations are directly reduced to a system of ordinary differential equations(ODEs) by only one step. It has been shown that not only do we recover many of the known results but also find some new reductions and explicit solutions, which may be previously unknown.展开更多
In this paper, a full discrete local projection stabilized (LPS) method is proposed to solve the optimal control problems of the unsteady Navier-Stokes equations with equal order elements. Convective effects and pre...In this paper, a full discrete local projection stabilized (LPS) method is proposed to solve the optimal control problems of the unsteady Navier-Stokes equations with equal order elements. Convective effects and pressure are both stabilized by using the LPS method. A priori error estimates uniformly with respect to the Reynolds number are obtained, providing the true solutions are sufficient smooth. Numerical experiments are implemented to illustrate and confirm our theoretical analysis.展开更多
In this paper, numerical prediction of ship motion responses in long-crest irregular waves by the URANS-VOF method is presented. A white noise spectrum is applied to generate the incoming waves to evaluate the motion ...In this paper, numerical prediction of ship motion responses in long-crest irregular waves by the URANS-VOF method is presented. A white noise spectrum is applied to generate the incoming waves to evaluate the motion responses. The procedure can replace a decade of simulations in regular wave with one single run to obtain a complete curve of linear motion response, considerably reducing computation time. A correction procedure is employed to adjust the wave generation signal based on the wave spectrum and achieves fairly better results in the wave tank. Three ship models with five wave conditions are introduced to validate the method. The computations in this paper are completed by using the solver naoe-FOAM-SJTU, a solver developed for ship and ocean engineering based on the open source code OpenFOAM. The computational motion responses by the irregular wave procedure are compared with the results by regular wave, experiments and strip theory. Transfer functions by irregular wave closely agree with the data obtained in the regular waves, showing negligible difference. The comparison between computational results and experiments also show good agreements. The results better predicted by CFD method than strip theories indicate that this method can compensate for the inaccuracy of the strip theories. The results confirm that the irregular wave procedure is a promising method for the accurate prediction of motion responses with less accuracy loss and higher efficiency compared with the regular wave procedure.展开更多
Based on the field velocity method,a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper.Responses to pitching and plunging motions of the airfoil are simulated...Based on the field velocity method,a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper.Responses to pitching and plunging motions of the airfoil are simulated under different conditions.The obtained results are compared with those of moving grid method and good agreement is achieved.In the conventional field velocity method,the Euler solver is usually used to simulate the movement of the airfoil.However,when viscous effect is considered,unsteady Navier-Stokes equations have to be solved and the viscous flux correction must be taken into account.In this work,the viscous flux correction is introduced into the conventional field velocity method when non-uniform grid speed distribution is occurred.Numerical experiments for the flow around NACA0012 airfoil showed that the proposed approach can well simulate viscous moving boundary flow problems.展开更多
基金Project supported in part by the National Natural Science Foundation of China(Grant No.11771259)Shaanxi Provincial Joint Laboratory of Artificial Intelligence(GrantNo.2022JCSYS05)+1 种基金Innovative Team Project of Shaanxi Provincial Department of Education(Grant No.21JP013)Shaanxi Provincial Social Science Fund Annual Project(Grant No.2022D332)。
文摘We propose the meshfree-based physics-informed neural networks for solving the unsteady Oseen equations.Firstly,based on the ideas of meshfree and small sample learning,we only randomly select a small number of spatiotemporal points to train the neural network instead of forming a mesh.Specifically,we optimize the neural network by minimizing the loss function to satisfy the differential operators,initial condition and boundary condition.Then,we prove the convergence of the loss function and the convergence of the neural network.In addition,the feasibility and effectiveness of the method are verified by the results of numerical experiments,and the theoretical derivation is verified by the relative error between the neural network solution and the analytical solution.
文摘Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be satisfied exactly without increasing the computational time. The novel scheme is useful for computing transonic, strong disturbance, unsteady flows with high reduced frequencies. The scheme makes use of curvefitted orthogonal meshes and the lattice control technique to obtain the optimal grid distribution. The numerical results are satisfactory.
文摘This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse the behaviour of axial velocity and shear stress for unsteady flow of nun-Newtonian visco-elastie fluids in tube, and to provide a theoretical base for the projection of pipe-line engineering.
基金The project is supported by the National Nature Science Foundation of China(NNSFC)
文摘The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.
基金financially supported by the National Natural Science Foundation of China(Grant No.51349011)the Foundation of Si’chuan Educational Committee(Grant No.17ZB0452)+1 种基金the Innovation Team Project of Si’chuan Educational Committee(Grant No.18TD0019)the Longshan Academic Talent Research Support Program of the Southwest of Science and Technology(Grant Nos.18LZX715 and 18LZX410)
文摘In order to solve unsteady incompressible Navier–Stokes(N–S) equations, a new stabilized finite element method,called the viscous-splitting least square FEM, is proposed. In the model, the N–S equations are split into diffusive and convective parts in each time step. The diffusive part is discretized by the backward difference method in time and discretized by the standard Galerkin method in space. The convective part is a first-order nonlinear equation.After the linearization of the nonlinear part by Newton’s method, the convective part is also discretized by the backward difference method in time and discretized by least square scheme in space. C0-type element can be used for interpolation of the velocity and pressure in the present model. Driven cavity flow and flow past a circular cylinder are conducted to validate the present model. Numerical results agree with previous numerical results, and the model has high accuracy and can be used to simulate problems with complex geometry.
基金supported by the National Natural Science Foundation of China(11232002)the Ph.D.Student Foundation of Chinese Ministry of Education(30400002011105001)
文摘The forward flight of a model butterfly was stud- ied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational forces, where the aerodynamic forces were generated by flapping wings which moved with the body, allowing the body os- cillations of the model butterfly to be simulated. The main results are as follows: (1) The aerodynamic force produced by the wings is approximately perpendicular to the long-axis of body and is much larger in the downstroke than in the up- stroke. In the downstroke the body pitch angle is small and the large aerodynamic force points up and slightly backward, giving the weight-supporting vertical force and a small neg- ative horizontal force, whilst in the upstroke, the body an- gle is large and the relatively small aerodynamic force points forward and slightly downward, giving a positive horizon- tal force which overcomes the body drag and the negative horizontal force generated in the downstroke. (2) Pitching oscillation of the butterfly body plays an equivalent role of the wing-rotation of many other insects. (3) The body-mass- specific power of the model butterfly is 33.3 W/kg, not very different from that of many other insects, e.g., fruitflies and dragonflies.
基金the Key Programs of Universities in Henan Province of China(22A140006)the Fundamental Research Funds for the Universities of Henan Province(NSFRF210324)+1 种基金Pro-gram of Henan Polytechnic University(B2018-40)Innovative Sci-entists and Technicians Team of Henan Provincial High Education(21IRTSTHN016).
文摘The aim of this paper is to study the unsteady korteweg-de vries equation that plays an important role in describing the shallow water.Two analytical techniques namely the Sardar-subequation method and the energy balance method are employed to seek the abundant traveling wave solutions for the first time.By these two methods,plenty of traveling wave solutions such as the bright solitary wave solutions,dark solitary wave solutions,singular periodic wave solutions and perfect periodic wave solution that expressed in terms of the generalized hyperbolic functions,generalized trigonometric functions and the cosine function are obtained.Finally,the dynamic behaviors of the solutions are described through the 3D plot and 2D curve.The results in this paper demonstrate that the proposed methods are powerful and effective to construct the traveling wave solutions of the nonlinear evolution equations in ocean engineering and science.
基金Project supported by the National Natural Science Foundation of China (Major Program of the Tenth Five-Year Plan) (No.90411006)the Post-Doctoral Science Foundation of Jiangsu Province of China(No.0602024C)
文摘Stability related to theoretical model for catastrophic weather prediction, which includes non-hydrostatic perfect elastic model and anelastic model, is discussed and analyzed in detail. It is proved that non-hydrostatic perfect elastic equations set is stable in the class of infinitely differentiable function. However, for the anelastic equations set, its continuity equation is changed in form because of the particular hypothesis for fluid, so "the matching consisting of both viscosity coefficient and incompressible assumption" appears, thereby the most important equations set of this class in practical prediction shows the same instability in topological property as Navier-Stokes equation, which should be avoided first in practical numerical prediction. In light of this, the referenced suggestions to amend the applied model are finally presented.
基金Supported by Projects 19472068 and 19772056 of the National Natural Science Foundation ofChina and the Laboratory of Scientifi
文摘In this paper, the Crank-Nicholson + component-consistent pressure correction method for the numerical solution of the unsteady incompressible Navier-Stokes equation of [1] on the rectangular half-Staggered mesh has been extended to the curvilinear half-Staggered mesh. The discrete projection, both for the projection step in the solution procedure and for the related differential-algebraic equations, has been carefully studied and verified. It is proved that the proposed method is also unconditionally (in t) nonlinearly stable on the curvilinear mesh, provided the mesh is not too skewed. It is seen that for problems with an outflow boundary, the half-Staggered mesh is especially advantageous. Results of preliminary numerical experiments support these claims.
文摘Fourier-Legendre spectral approximation for the unsteady Navier-Stokes equations is analyzed. The generalized stability and convergence are proved respectively.
基金This work is supported by national numerical wind tunnel project under contract number 2018-ZT4A072016YFA0401200 of national key research and development program of China.
文摘The numerical methods for computing the stability derivatives of the aircraft by solving unsteady sensitivity equations which was proposed in our previous papers was extended to solve three-dimensional problems in this paper.Both the static and dynamic derivatives of the hypersonic blunt cone undergoing pitching oscillation around a fixed point were computed using the new methods.The predicted static derivative and dynamic derivative were found to be in reasonable agreement with the experimental data.For the present method,it is possible to distinguish the components of dynamic derivatives caused by different state parameters.It is found that C_(m_α) and C_(mq) are usually of opposite signs and tend to eliminate each other,which makes C_(m_α)+C_(mq) being much smaller than its individual components.Another feature of this method is that the moment of pressure derivatives proposed in the present paper can be used to predict the contribution of each part of the blunt cone to the overall stability quantitatively.It is found that the head region is crucial for the static stability and the body region contributes most to the dynamic stability.
基金Supported by the Global Change Research Program of China under Grant No.2015CB953904National Natural Science Foundation of China under Grant Nos.11275072,11435005,11675054Shanghai Collaborative Innovation Center of Trustworthy Software for Internet of Things under Grant No.ZF1213
文摘To find intrinsically different symmetry reductions and inequivalent group invariant solutions of the 2D unsteady incompressible boundary-layer equations, a two-dimensional optimal system is constructed which attributed to the classification of the corresponding Lie subalgebras. The comprehensiveness and inequivalence of the optimal system are shown clearly under different values of invariants. Then by virtue of the optimal system obtained, the boundary-layer equations are directly reduced to a system of ordinary differential equations(ODEs) by only one step. It has been shown that not only do we recover many of the known results but also find some new reductions and explicit solutions, which may be previously unknown.
基金This work is supported by the Natural Science Foundation of China (No. 11271273) and the Scientific Research Foundation of the Education Department of Sichuan Province of China (No.16ZB0300). The authors would like to thank the associate editor and anonymous referees comments to improve the quality of the manuscript.
文摘In this paper, a full discrete local projection stabilized (LPS) method is proposed to solve the optimal control problems of the unsteady Navier-Stokes equations with equal order elements. Convective effects and pressure are both stabilized by using the LPS method. A priori error estimates uniformly with respect to the Reynolds number are obtained, providing the true solutions are sufficient smooth. Numerical experiments are implemented to illustrate and confirm our theoretical analysis.
基金Y.Guo's research was supported in part by NSF grant 1810868supported in part by NSF DMS-1501031,DMS-1900923+1 种基金the University of WisconsinMadison Graduate School with funding from the Wisconsin Alumni Research Foundationpartially supported by MIUR-Prin。
文摘The goal of this paper is to study the important diffusive expansion via an alternative mathematical approach other than that in [21].
基金supported by the National Natural Science Foundation of China(Grant Nos.51379125,11272120)the National Key Basic Research Development Program of China(973Program,Grant No.2013CB036103)the High Technology of Marine Research Project of the Ministry of Industry and Information Technology of China
文摘In this paper, numerical prediction of ship motion responses in long-crest irregular waves by the URANS-VOF method is presented. A white noise spectrum is applied to generate the incoming waves to evaluate the motion responses. The procedure can replace a decade of simulations in regular wave with one single run to obtain a complete curve of linear motion response, considerably reducing computation time. A correction procedure is employed to adjust the wave generation signal based on the wave spectrum and achieves fairly better results in the wave tank. Three ship models with five wave conditions are introduced to validate the method. The computations in this paper are completed by using the solver naoe-FOAM-SJTU, a solver developed for ship and ocean engineering based on the open source code OpenFOAM. The computational motion responses by the irregular wave procedure are compared with the results by regular wave, experiments and strip theory. Transfer functions by irregular wave closely agree with the data obtained in the regular waves, showing negligible difference. The comparison between computational results and experiments also show good agreements. The results better predicted by CFD method than strip theories indicate that this method can compensate for the inaccuracy of the strip theories. The results confirm that the irregular wave procedure is a promising method for the accurate prediction of motion responses with less accuracy loss and higher efficiency compared with the regular wave procedure.
基金This work was supported by The National Basic Research Program of China(Grant No.2007CB714600)Funding of Jiangsu Innovation Program for Graduate Education(Grant No.CXLX110170).
文摘Based on the field velocity method,a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper.Responses to pitching and plunging motions of the airfoil are simulated under different conditions.The obtained results are compared with those of moving grid method and good agreement is achieved.In the conventional field velocity method,the Euler solver is usually used to simulate the movement of the airfoil.However,when viscous effect is considered,unsteady Navier-Stokes equations have to be solved and the viscous flux correction must be taken into account.In this work,the viscous flux correction is introduced into the conventional field velocity method when non-uniform grid speed distribution is occurred.Numerical experiments for the flow around NACA0012 airfoil showed that the proposed approach can well simulate viscous moving boundary flow problems.