Abstract: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.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a Volume Penalization (VP) method for handling inhomogeneous Neumann boundary conditions in complex geometries. Specifically, the author aims to extend the existing Volume Penalization method for homogeneous Neumann boundary conditions so as to be able to simulate scalar fluxes through walls. This method allows the use of simple computational domains (such as Cartesian domains) to solve the governing equations, thereby simplifying numerical simulations in complex geometries. ### Main Problems and Solutions 1. **Limitations of Existing Methods**: - Numerical modeling of multi - physical problems in complex geometries remains a challenge in computational fluid dynamics. - In particular, for problems such as heat and mass exchange, traditional numerical methods are difficult to handle accurate simulations under complex boundary conditions. 2. **Advantages of the Volume Penalty Method**: - The volume penalty method treats the solid as a porous medium and introduces a penalty parameter \(\eta\) to handle boundary conditions, making it possible to impose complex boundary conditions in a simple computational domain. - This method can effectively simulate scalar fluxes through walls while maintaining flux continuity. 3. **Specific Implementation**: - The paper first verifies the effectiveness of this method on the one - dimensional Poisson equation and analyzes the penalty error and discretization error. - Then, the method is extended to the two - dimensional Poisson equation and applied to more complex geometries. - Finally, the method is applied to the two - dimensional steady - state free convection problem coupling the Navier - Stokes equations and the convection - diffusion equations. ### Mathematical Formula Representation - **One - Dimensional Poisson Equation**: \[ -\frac{d^2 w(x)}{dx^2}=f(x) \quad \text{for } x \in \Omega_f \] where \(\Omega_f = \{x|0 < x < \pi\}\), and inhomogeneous Neumann boundary conditions are imposed: \[ \left.\frac{dw}{dx}\right|_{x = 0}=\left.\frac{dw}{dx}\right|_{x=\pi}=\alpha \] - **Volume Penalty Representation**: \[ -\frac{d}{dx}\left(\theta\frac{dv}{dx}+\chi\alpha\right)=(1 - \chi)f \] where, \[ \theta(x)=1-\chi(x)+\eta\chi(x) \] \(\chi(x)\) is a mask function, defined as: \[ \chi(x)= \begin{cases} 0 & \text{for } 0 < x < \pi\\ \frac{1}{2} & \text{at } x = 0,\pi\\ 1 & \text{for } \pi < x < 2\pi \end{cases} \] ### Summary This paper successfully extends the volume penalty method so that it can handle inhomogeneous Neumann boundary conditions, thus providing an effective method for scalar flux simulation in complex geometries. Through detailed mathematical derivations and numerical experiments, the feasibility and accuracy of this method are proved.

Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry

Pressure-tight and non-stiff volume penalization for compressible flows

An effective preconditioning strategy for volume penalized incompressible/low Mach multiphase flow solvers

Penalty method for the Navier-Stokes-Fourier system with Dirichlet boundary conditions: convergence and error estimates

Immersed boundary method for high-order flux reconstruction based on volume penalization

A combined volume penalization / selective frequency damping approach for immersed boundary methods: application to moving geometries

Penalization of stationary Navier-Stokes equations and applications in topology optimization

Generalization of the Penalized Wall Function Method for Modeling of Turbulent Flows with Adverse Pressure Gradient

Convergence and error estimates of a penalization finite volume method for the compressible Navier-Stokes system

Compressible flow simulation with moving geometries using the Brinkman penalization in high-order Discontinuous Galerkin

A combined volume penalization / selective frequency damping approach for immersed boundary methods applied to high-order schemes

Penalization of stationary Navier–Stokes equations and applications in topology optimization

A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

On the Calculation of the Brinkman Penalization Term in Density-Based Topology Optimization of Fluid-Dependent Problems

An Interior Penalty coupling strategy for Isogeometric non-conformal Kirchhoff-Love shell patches

$d$-dimensional extension of a penalization method for Neumann or Robin boundary conditions: a boundary layer approach and numerical experiments

Scalar conservation laws with stochastic discontinuous flux function

Maxwell boundary condition for discrete velocity methods: macroscopic physical constraints and Lagrange multiplier-based implementation

Numerical Study of Lattice Boltzmann Methods for a Convection–diffusion Equation Coupled with Navier–Stokes Equations

A simple diffuse interface approach for compressible flows around moving solids of arbitrary shape based on a reduced Baer-Nunziato model

Finite element discretization of the steady, generalized Navier-Stokes equations with inhomogeneous Dirichlet boundary conditions