In this paper,we describe how to construct a finite-difference shockcapturing method for the numerical solution of the Euler equation of gas dynamics on arbitrary two-dimensional domainΩ,possibly with moving boundary...In this paper,we describe how to construct a finite-difference shockcapturing method for the numerical solution of the Euler equation of gas dynamics on arbitrary two-dimensional domainΩ,possibly with moving boundary.The boundaries of the domain are assumed to be changing due to the movement of solid objects/obstacles/walls.Although the motion of the boundary could be coupled with the fluid,all of the numerical tests are performed assuming that such a motion is prescribed and independent of the fluid flow.The method is based on discretizing the equation on a regular Cartesian grid in a rectangular domainΩ_(R)⊃Ω.Ωe identify inner and ghost points.The inner points are the grid points located insideΩ,while the ghost points are the grid points that are outsideΩbut have at least one neighbor insideΩ.The evolution equations for inner points data are obtained from the discretization of the governing equation,while the data at the ghost points are obtained by a suitable extrapolation of the primitive variables(density,velocities and pressure).Particular care is devoted to a proper description of the boundary conditions for both fixed and time dependent domains.Several numerical experiments are conducted to illustrate the validity of themethod.Ωe demonstrate that the second order of accuracy is numerically assessed on genuinely two-dimensional problems.展开更多
We propose a finite-difference ghost-point approach for the numerical solution of Cauchy-Navier equations in linear elasticity problems on arbitrary unbounded domains.The technique is based on a smooth coordinate tran...We propose a finite-difference ghost-point approach for the numerical solution of Cauchy-Navier equations in linear elasticity problems on arbitrary unbounded domains.The technique is based on a smooth coordinate transformation,which maps an unbounded domain into a unit square.Arbitrary geometries are defined by suitable level-set functions.The equations are discretized by classical nine-point stencil on interior points,while boundary conditions and high order reconstructions are used to define the field variables at ghost-points,which are grid nodes external to the domain with a neighbor inside the domain.The linear system arising from such discretization is solved by a multigrid strategy.The approach is then applied to solve elasticity problems in volcanology for computing the displacement caused by pressure sources.The method is suitable to treat problems in which the geometry of the source often changes(explore the effects of different scenarios,or solve inverse problems in which the geometry itself is part of the unknown),since it does not require complex re-meshing when the geometry is modified.Several numerical tests are successfully performed,which asses the effectiveness of the present approach.展开更多
基金The work of A.Chertock was supported in part by the NSF Grants DMS-1216974 and DMS-1521051The work of A.Kurganov was supported in part by the NSF Grants DMS-1216957 and DMS-1521009The work of G.Russo was supported partially by the University of Catania,Project F.I.R.Charge Transport in Graphene and Low Dimensional Systems,and partially by ITN-ETN Horizon 2020 Project Mod Comp Shock,Modeling and Computation on Shocks and Interfaces,Project Reference 642768.
文摘In this paper,we describe how to construct a finite-difference shockcapturing method for the numerical solution of the Euler equation of gas dynamics on arbitrary two-dimensional domainΩ,possibly with moving boundary.The boundaries of the domain are assumed to be changing due to the movement of solid objects/obstacles/walls.Although the motion of the boundary could be coupled with the fluid,all of the numerical tests are performed assuming that such a motion is prescribed and independent of the fluid flow.The method is based on discretizing the equation on a regular Cartesian grid in a rectangular domainΩ_(R)⊃Ω.Ωe identify inner and ghost points.The inner points are the grid points located insideΩ,while the ghost points are the grid points that are outsideΩbut have at least one neighbor insideΩ.The evolution equations for inner points data are obtained from the discretization of the governing equation,while the data at the ghost points are obtained by a suitable extrapolation of the primitive variables(density,velocities and pressure).Particular care is devoted to a proper description of the boundary conditions for both fixed and time dependent domains.Several numerical experiments are conducted to illustrate the validity of themethod.Ωe demonstrate that the second order of accuracy is numerically assessed on genuinely two-dimensional problems.
基金the OTRIONS project under the European Territorial Cooperation Programme Greece-Italy 2007-2013,and by PRIN 2009“Innovative numerical methods for hyperbolic problems with applications to fluid dynamics,kinetic theory and computational biology”.
文摘We propose a finite-difference ghost-point approach for the numerical solution of Cauchy-Navier equations in linear elasticity problems on arbitrary unbounded domains.The technique is based on a smooth coordinate transformation,which maps an unbounded domain into a unit square.Arbitrary geometries are defined by suitable level-set functions.The equations are discretized by classical nine-point stencil on interior points,while boundary conditions and high order reconstructions are used to define the field variables at ghost-points,which are grid nodes external to the domain with a neighbor inside the domain.The linear system arising from such discretization is solved by a multigrid strategy.The approach is then applied to solve elasticity problems in volcanology for computing the displacement caused by pressure sources.The method is suitable to treat problems in which the geometry of the source often changes(explore the effects of different scenarios,or solve inverse problems in which the geometry itself is part of the unknown),since it does not require complex re-meshing when the geometry is modified.Several numerical tests are successfully performed,which asses the effectiveness of the present approach.
