5 research outputs found
Analysis and discretization of the volume penalized Laplace operator with Neumann boundary conditions
We study the properties of an approximation of the Laplace operator with
Neumann boundary conditions using volume penalization. For the one-dimensional
Poisson equation we compute explicitly the exact solution of the penalized
equation and quantify the penalization error. Numerical simulations using
finite differences allow then to assess the discretisation and penalization
errors. The eigenvalue problem of the penalized Laplace operator with Neumann
boundary conditions is also studied. As examples in two space dimensions, we
consider a Poisson equation with Neumann boundary conditions in rectangular and
circular domains
Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry
We develop a volume penalization method for inhomogeneous Neumann boundary
conditions, generalizing the flux-based volume penalization method for
homogeneous Neumann boundary condition proposed by Kadoch et al. [J. Comput.
Phys. 231 (2012) 4365]. The generalized method allows us to model scalar flux
through walls in geometries of complex shape using simple, e.g. Cartesian,
domains for solving the governing equations. We examine the properties of the
method, by considering a one-dimensional Poisson equation with different
Neumann boundary conditions. The penalized Laplace operator is discretized by
second order central finite-differences and interpolation. The discretization
and penalization errors are thus assessed for several test problems.
Convergence properties of the discretized operator and the solution of the
penalized equation are analyzed. The generalized method is then applied to an
advection-diffusion equation coupled with the Navier-Stokes equations in an
annular domain which is immersed in a square domain. The application is
verified by numerical simulation of steady free convection in a concentric
annulus heated through the inner cylinder surface using an extended square
domain.Comment: 32 pages, 19 figure
An Efficient Numerical Approach for Solving Nonlinear Coupled Hyperbolic Partial Differential Equations with Nonlocal Conditions
One of the most important advantages of collocation method is the possibility of dealing with nonlinear partial differential equations (PDEs) as well as PDEs with variable coefficients. A numerical solution based on a Jacobi collocation method is extended to solve nonlinear coupled hyperbolic PDEs with variable coefficients subject to initial-boundary nonlocal conservation conditions. This approach, based on Jacobi polynomials and Gauss-Lobatto quadrature integration, reduces solving the nonlinear coupled hyperbolic PDEs with variable coefficients to a system of nonlinear ordinary differential equation which is far easier to solve. In fact, we deal with initial-boundary coupled hyperbolic PDEs with variable coefficients as well as initial-nonlocal conditions. Using triangular, soliton, and exponential-triangular solutions as exact solutions, the obtained results show that the proposed numerical algorithm is efficient and very accurate
A pseudo-spectral method with volume penalisation for magnetohydrodynamic turbulence in confined domains
International audienceWe present a Fourier pseudo-spectral method for solving the resistive magnetohydrodynamic equations of incompressible flow in confined domains. A volume penalisation method allows to take into account boundary conditions and the geometry of the domain. A code validation is presented for the z-pinch test case. Numerical simulations of decaying MHD turbulence in two space dimensions show spontaneous spin-up of the flow in non-axisymmetric geometries, which is reflected by the generation of angular momentum. First results of decaying MHD turbulence in a cylinder illustrate the potential of the new method for three-dimensional simulations