Second Order Multigrid Methods for Elliptic Problems with Discontinuous Coefficients on an Arbitrary Interface, I: One Dimensional Problems

In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface.Second order accuracy for the first derivative is obtained as well.The method is based on the Ghost Fluid Method,making use of ghost points on which the value is defined by suitable interface conditions.The multi-domain formulation is adopted,where the problem is split in two sub-problems and interface conditions will be enforced to close the problem.Interface conditions are relaxed together with the internal equations(following the approach proposed in[10]in the case of smooth coefficients),leading to an iterative method on all the set of grid values(inside points and ghost points).A multigrid approach with a suitable definition of the restriction operator is provided.The restriction of the defect is performed separately for both sub-problems,providing a convergence factor close to the one measured in the case of smooth coefficient and independent on the magnitude of the jump in the coefficient.Numerical tests will confirm the second order accuracy.Although the method is proposed in one dimension,the extension in higher dimension is currently underway[12]and it will be carried out by combining the discretization of[10]with the multigrid approach of[11]for Elliptic problems with non-eliminated boundary conditions in arbitrary domain.

A Second Order Finite-Difference Ghost-Point Method for Elasticity Problems on Unbounded Domains with Applications to Volcanology

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.

ASecond-Order Finite-Difference Method for Compressible Fluids in Domains with Moving Boundaries

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.