Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems,whose basic concept is to embed physical laws to constrain/inform neural networks,with the need of l...Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems,whose basic concept is to embed physical laws to constrain/inform neural networks,with the need of less data for training a reliable model.This can be achieved by incorporating the residual of physics equations into the loss function.Through minimizing the loss function,the network could approximate the solution.In this paper,we propose a mixed-variable scheme of physics-informed neural network(PINN)for fluid dynamics and apply it to simulate steady and transient laminar flows at low Reynolds numbers.A parametric study indicates that the mixed-variable scheme can improve the PINN trainability and the solution accuracy.The predicted velocity and pressure fields by the proposed PINN approach are also compared with the reference numerical solutions.Simulation results demonstrate great potential of the proposed PINN for fluid flow simulation with a high accuracy.展开更多
For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process...For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element(NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.展开更多
A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses parti...A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses particles to replace the mesh in traditional methods, the governing equations are discretized by virtue of the relationship of particles, and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method （ICCG）, the free surface is tracked by the change of numerical density. A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF, and an additional modification step was added to make the simulation more stable. The results show that the MPS method is suitable for the simulation of viscous liquid sloshing, with the advantage in arranging the particles easily, especially on some complex curved surface.展开更多
Smoothed particle hydrodynamics （SPH） is a Lagrangian, meshfree particle method and has been widely applied to diffe- rent areas in engineering and science. Since its original extension to modeling free surface flow...Smoothed particle hydrodynamics （SPH） is a Lagrangian, meshfree particle method and has been widely applied to diffe- rent areas in engineering and science. Since its original extension to modeling free surface flows by Monaghan in 1994, SPH has been gradually developed into an attractive approach for modeling viscous incompressible fluid flows. This paper presents an overview on the recent progresses of SPH in modeling viscous incompressible flows in four major aspects which are closely related to the computational accuracy of SPH simulations. The advantages and disadvantages of different SPH particle approximation sche- mes, pressure field solution approaches, solid boundary treatment algorithms and particle adapting algorithms are described and analyzed. Some new perspectives and fuRtre trends in SPH modeling of viscous incompressible flows are discussed.展开更多
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.展开更多
Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is c...Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model（MPS or HMPS） in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model（Enhanced MPS）, thus, the developed coupled solvers for fluid structure interaction（FSI） are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics（SPH）-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH（ISPH） fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.展开更多
In this article a finite volume method is proposed to solve viscous incompressible Navier-Stokes equations in two-dimensional regions with corners and curved boundaries. A hybrid collocated-grid variable arrangement i...In this article a finite volume method is proposed to solve viscous incompressible Navier-Stokes equations in two-dimensional regions with corners and curved boundaries. A hybrid collocated-grid variable arrangement is adopted, in which the velocity and pressure are stored at the centroid and the circumcenters of the triangular control cell, respectively. The cell flux is defined at the mid-point of the cell face. Second-order implicit time integration schemes are used for convection and diffusion terms. The second-order upwind scheme is used for convection fluxes. The present method is validated by results of several viscous flows.展开更多
Dynamic modeling for incompressible hyperelastic materials with large deformation is an important issue in biomimetic applications. The previously proposed lower-order fully parameterized absolute nodal coordinate for...Dynamic modeling for incompressible hyperelastic materials with large deformation is an important issue in biomimetic applications. The previously proposed lower-order fully parameterized absolute nodal coordinate formulation（ANCF） beam element employs cubic interpolation in the longitudinal direction and linear interpolation in the transverse direction, whereas it cannot accurately describe the large bending deformation. On this account, a novel modeling method for studying the dynamic behavior of nonlinear materials is proposed in this paper. In this formulation, a higher-order beam element characterized by quadratic interpolation in the transverse directions is used in this investigation. Based on the Yeoh model and volumetric energy penalty function, the nonlinear elastic force matrices are derived within the ANCF framework. The feasibility and availability of the Yeoh model are verified through static experiment of nonlinear incompressible materials. Furthermore,dynamic simulation of a silicone cantilever beam under the gravity force is implemented to validate the superiority of the higher-order beam element. The simulation results obtained based on the Yeoh model by employing three different ANCF beam elements are compared with the result achieved from a commercial finite element package as the reference result. It is found that the results acquired utilizing a higher-order beam element are in good agreement with the reference results,while the results obtained using a lower-order beam element are different from the reference results. In addition, the stiffening problem caused by volumetric locking can be resolved effectively by applying a higher-order beam element. It is concluded that the proposed higher-order beam element formulation has satisfying accuracy in simulating dynamic motion process of the silicone beam.展开更多
The method of semi-false transient was used to numerically compute the incompressible steady flow into the porous medium in this paper. The fundamental equations were established and numerically solved in the united f...The method of semi-false transient was used to numerically compute the incompressible steady flow into the porous medium in this paper. The fundamental equations were established and numerically solved in the united flow field, which included the space region and the porous region. The non-equidistant non-orthogonal semi-staggered mesh system was used in the method of semi-false transient. The computational results of two problems concerning the flow into porous medium from space region, in which there was the backward flow besides main flow, were obtained and discussed. It is seen from the results that the backward flow is generally not present in the porous medium as the osmotic resistance is very large.展开更多
Three-dimensional numerical simulation of a uniform incompressible viscous flow around a stationary circular cylinder was conducted. The CFX-4 software was used to calculate the hydrodynamic characteristics of the flo...Three-dimensional numerical simulation of a uniform incompressible viscous flow around a stationary circular cylinder was conducted. The CFX-4 software was used to calculate the hydrodynamic characteristics of the flow and the finite volume method for incompressible Navier-Stokes equations was employed in the program. The simulation of the flow was performed for Re=1013 and RE=104 respectively within the sub-critical region. In order to overcome numerical instability for the high Reynolds number flows, a quadratic upwind scheme was incorporated for the Navier-stokes equations. The periodicity boundary condition was used at the ends of the cylinder. It was found that the evolution of the lift and drag coefficients in each plane along the cylinder span is different. Comparison between the predicted results based on the three-dimensional and the two-dimensional analysis was also given. It is concluded that at high Reynolds number the effect of three-dimensionality of the flow around the circular cylinder is remarkable, and in addition hydrodynamic coefficients with of 3-D simulation are less than those given by 2-D simulation.展开更多
A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the gov...A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.展开更多
This paper is joint with [27]. The authors prove in this article the existence and reveal its structure of uniform attractor for a two-dimensional nonautonomous incompressible non-Newtonian fluid with a new class of e...This paper is joint with [27]. The authors prove in this article the existence and reveal its structure of uniform attractor for a two-dimensional nonautonomous incompressible non-Newtonian fluid with a new class of external forces.展开更多
Menasco showed that a closed incompressible surface in the complement of a non-split prime alternating link in S^(3) contains a circle isotopic in the link complement to a meridian of the links.Based on this result,he...Menasco showed that a closed incompressible surface in the complement of a non-split prime alternating link in S^(3) contains a circle isotopic in the link complement to a meridian of the links.Based on this result,he was able to argue the hyperbolicity of non-split prime alternating links in S3.Adams et al.showed that if F■S×I\L is an essential torus,then F contains a circle which is isotopic in S×I\L to a meridian of L.The author generalizes his result as follows:Let S be a closed orientable surface,L be a fully alternating link in S×I\If F ■ S×I\L is a closed essential surface,then F contains a circle which is isotopic in S×I\L to a meridian of L.展开更多
The radial symmetric motion problem was examined for a spherical shell composed of a class of imperfect incompressible hyper-elastic materials, in which the materials may be viewed as the homogeneous incompressible is...The radial symmetric motion problem was examined for a spherical shell composed of a class of imperfect incompressible hyper-elastic materials, in which the materials may be viewed as the homogeneous incompressible isotropic neo-Hookean material with radial perturbations. A second-order nonlinear ordinary differential equation that describes the radial motion of the inner surface of the shell was obtained. And the first integral of the equation was then carded out. Via analyzing the dynamical properties of the solution of the differential equation, the effects of the prescribed imperfection parameter of the material and the ratio of the inner and the outer radii of the underformed shell on the motion of the inner surface of the shell were discussed, and the corresponding numerical examples were carded out simultaneously. In particular, for some given parameters, it was proved that, there exists a positive critical value, and the motion of the inner surface with respect to time will present a nonlinear periodic oscillation as the difference between the inner and the outer presses does not exceed the critical value. However, as the difference exceeds the critical value, the motion of the inner surface with respect to time will increase infinitely. That is to say, the shell will be destroyed ultimately.展开更多
In this paper,we prove that there exists a unique local solution for the Cauchy problem of a system of the incompressible Navier-Stokes-Landau-Lifshitz equations with the Dzyaloshinskii-Moriya interaction and V-flow t...In this paper,we prove that there exists a unique local solution for the Cauchy problem of a system of the incompressible Navier-Stokes-Landau-Lifshitz equations with the Dzyaloshinskii-Moriya interaction and V-flow term inR^(2) and R^(3).Our methods rely upon approximating the system with a perturbed parabolic system and parallel transport.展开更多
Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of ...Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.展开更多
文摘Physics-informed deep learning has drawn tremendous interest in recent years to solve computational physics problems,whose basic concept is to embed physical laws to constrain/inform neural networks,with the need of less data for training a reliable model.This can be achieved by incorporating the residual of physics equations into the loss function.Through minimizing the loss function,the network could approximate the solution.In this paper,we propose a mixed-variable scheme of physics-informed neural network(PINN)for fluid dynamics and apply it to simulate steady and transient laminar flows at low Reynolds numbers.A parametric study indicates that the mixed-variable scheme can improve the PINN trainability and the solution accuracy.The predicted velocity and pressure fields by the proposed PINN approach are also compared with the reference numerical solutions.Simulation results demonstrate great potential of the proposed PINN for fluid flow simulation with a high accuracy.
基金supported by the National Natural Science Foundations of China (Grant 11502286)
文摘For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element(NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.
基金the National Natural Science Foundation under Grant No.50579035
文摘A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses particles to replace the mesh in traditional methods, the governing equations are discretized by virtue of the relationship of particles, and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method （ICCG）, the free surface is tracked by the change of numerical density. A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF, and an additional modification step was added to make the simulation more stable. The results show that the MPS method is suitable for the simulation of viscous liquid sloshing, with the advantage in arranging the particles easily, especially on some complex curved surface.
基金Project supported by the National Natural Science Foun-dation of China(Grant Nos.11172306,U1530110)the Institu-te of Systems Engineering,China Academy of Engineering Physics(Grant No.2013KJZ01)
文摘Smoothed particle hydrodynamics （SPH） is a Lagrangian, meshfree particle method and has been widely applied to diffe- rent areas in engineering and science. Since its original extension to modeling free surface flows by Monaghan in 1994, SPH has been gradually developed into an attractive approach for modeling viscous incompressible fluid flows. This paper presents an overview on the recent progresses of SPH in modeling viscous incompressible flows in four major aspects which are closely related to the computational accuracy of SPH simulations. The advantages and disadvantages of different SPH particle approximation sche- mes, pressure field solution approaches, solid boundary treatment algorithms and particle adapting algorithms are described and analyzed. Some new perspectives and fuRtre trends in SPH modeling of viscous incompressible flows are discussed.
基金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.
文摘Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model（MPS or HMPS） in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model（Enhanced MPS）, thus, the developed coupled solvers for fluid structure interaction（FSI） are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics（SPH）-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH（ISPH） fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.
基金Project supported by the National Natural Science Foundation of China(Grant No.10771134).
文摘In this article a finite volume method is proposed to solve viscous incompressible Navier-Stokes equations in two-dimensional regions with corners and curved boundaries. A hybrid collocated-grid variable arrangement is adopted, in which the velocity and pressure are stored at the centroid and the circumcenters of the triangular control cell, respectively. The cell flux is defined at the mid-point of the cell face. Second-order implicit time integration schemes are used for convection and diffusion terms. The second-order upwind scheme is used for convection fluxes. The present method is validated by results of several viscous flows.
基金supported by the National Natural Science Foundation of China (11772186 and 11272203)
文摘Dynamic modeling for incompressible hyperelastic materials with large deformation is an important issue in biomimetic applications. The previously proposed lower-order fully parameterized absolute nodal coordinate formulation（ANCF） beam element employs cubic interpolation in the longitudinal direction and linear interpolation in the transverse direction, whereas it cannot accurately describe the large bending deformation. On this account, a novel modeling method for studying the dynamic behavior of nonlinear materials is proposed in this paper. In this formulation, a higher-order beam element characterized by quadratic interpolation in the transverse directions is used in this investigation. Based on the Yeoh model and volumetric energy penalty function, the nonlinear elastic force matrices are derived within the ANCF framework. The feasibility and availability of the Yeoh model are verified through static experiment of nonlinear incompressible materials. Furthermore,dynamic simulation of a silicone cantilever beam under the gravity force is implemented to validate the superiority of the higher-order beam element. The simulation results obtained based on the Yeoh model by employing three different ANCF beam elements are compared with the result achieved from a commercial finite element package as the reference result. It is found that the results acquired utilizing a higher-order beam element are in good agreement with the reference results,while the results obtained using a lower-order beam element are different from the reference results. In addition, the stiffening problem caused by volumetric locking can be resolved effectively by applying a higher-order beam element. It is concluded that the proposed higher-order beam element formulation has satisfying accuracy in simulating dynamic motion process of the silicone beam.
文摘The method of semi-false transient was used to numerically compute the incompressible steady flow into the porous medium in this paper. The fundamental equations were established and numerically solved in the united flow field, which included the space region and the porous region. The non-equidistant non-orthogonal semi-staggered mesh system was used in the method of semi-false transient. The computational results of two problems concerning the flow into porous medium from space region, in which there was the backward flow besides main flow, were obtained and discussed. It is seen from the results that the backward flow is generally not present in the porous medium as the osmotic resistance is very large.
文摘Three-dimensional numerical simulation of a uniform incompressible viscous flow around a stationary circular cylinder was conducted. The CFX-4 software was used to calculate the hydrodynamic characteristics of the flow and the finite volume method for incompressible Navier-Stokes equations was employed in the program. The simulation of the flow was performed for Re=1013 and RE=104 respectively within the sub-critical region. In order to overcome numerical instability for the high Reynolds number flows, a quadratic upwind scheme was incorporated for the Navier-stokes equations. The periodicity boundary condition was used at the ends of the cylinder. It was found that the evolution of the lift and drag coefficients in each plane along the cylinder span is different. Comparison between the predicted results based on the three-dimensional and the two-dimensional analysis was also given. It is concluded that at high Reynolds number the effect of three-dimensionality of the flow around the circular cylinder is remarkable, and in addition hydrodynamic coefficients with of 3-D simulation are less than those given by 2-D simulation.
基金Supported by the National Natural Science Foundation of China(11272153)
文摘A lattice Boltzmann flux solver(LBFS)is presented for simulation of fluid flows.Like the conventional computational fluid dynamics(CFD)solvers,the new solver also applies the finite volume method to discretize the governing differential equations,but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers.Instead,it is evaluated from local solution of lattice Boltzmann equation(LBE)at cell interface.Two versions of LBFS are presented in this paper.One is to locally apply one-dimensional compressible lattice Boltzmann(LB)model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves.The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows.The present solver removes the drawbacks of conventional lattice Boltzmann method(LBM)such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.
基金Sponsored by the NSFC (10901121,10826091 and 10771139)NSF for Postdoctors of China (20090460952)+2 种基金NSF of Zhejiang Province (Y6080077)NSF of Wenzhou University (2008YYLQ01)by the Zhejiang Youth Teacher Training Project and Wenzhou 551 Project
文摘This paper is joint with [27]. The authors prove in this article the existence and reveal its structure of uniform attractor for a two-dimensional nonautonomous incompressible non-Newtonian fluid with a new class of external forces.
文摘Menasco showed that a closed incompressible surface in the complement of a non-split prime alternating link in S^(3) contains a circle isotopic in the link complement to a meridian of the links.Based on this result,he was able to argue the hyperbolicity of non-split prime alternating links in S3.Adams et al.showed that if F■S×I\L is an essential torus,then F contains a circle which is isotopic in S×I\L to a meridian of L.The author generalizes his result as follows:Let S be a closed orientable surface,L be a fully alternating link in S×I\If F ■ S×I\L is a closed essential surface,then F contains a circle which is isotopic in S×I\L to a meridian of L.
基金国家自然科学基金，Municipal Key Subject Program of Shanghai
文摘The radial symmetric motion problem was examined for a spherical shell composed of a class of imperfect incompressible hyper-elastic materials, in which the materials may be viewed as the homogeneous incompressible isotropic neo-Hookean material with radial perturbations. A second-order nonlinear ordinary differential equation that describes the radial motion of the inner surface of the shell was obtained. And the first integral of the equation was then carded out. Via analyzing the dynamical properties of the solution of the differential equation, the effects of the prescribed imperfection parameter of the material and the ratio of the inner and the outer radii of the underformed shell on the motion of the inner surface of the shell were discussed, and the corresponding numerical examples were carded out simultaneously. In particular, for some given parameters, it was proved that, there exists a positive critical value, and the motion of the inner surface with respect to time will present a nonlinear periodic oscillation as the difference between the inner and the outer presses does not exceed the critical value. However, as the difference exceeds the critical value, the motion of the inner surface with respect to time will increase infinitely. That is to say, the shell will be destroyed ultimately.
文摘In this paper,we prove that there exists a unique local solution for the Cauchy problem of a system of the incompressible Navier-Stokes-Landau-Lifshitz equations with the Dzyaloshinskii-Moriya interaction and V-flow term inR^(2) and R^(3).Our methods rely upon approximating the system with a perturbed parabolic system and parallel transport.
文摘Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.
