A Moving Mesh Method for Porous Medium Equation by the Onsager Variational Principle

Read original: arXiv:2403.20030 - Published 4/1/2024 by Si Xiao, Xianmin Xu

🔮

Introduction

The porous medium equation (PME) is a prominent mathematical model used to describe various physical and biological phenomena, including gas flow, nonlinear heat transport, and groundwater movement. The PME is a nonlinear partial differential equation with a finite speed of propagation, in contrast to the linear heat equation where propagation speed is infinite.

The PME exhibits several interesting properties, including the formation of a free boundary and a waiting time phenomenon where the free boundary remains stationary until a critical time threshold is reached. Solving the PME numerically presents several challenges, such as capturing the free boundary, dealing with singularities, and computing the waiting time accurately.

Early numerical methods often transformed the PME into an alternative pressure equation, which exhibited better regularity. More recent approaches have utilized high-order local discontinuous Galerkin methods, positivity-preserving limiters, and moving mesh methods to address the complexities of the PME. Lagrangian-type methods, formulated in a reference domain, have also been developed, leveraging the energy variational method and Wasserstein gradient flow theory.

The paper aims to devise a numerical method that combines the strengths of both conventional moving mesh methods and Lagrangian-type methods by utilizing the Onsager variational principle as an approximation tool. The authors derive the PME from the Onsager principle and then formulate semi-discrete numerical schemes in mixed form, demonstrating their energy dissipation properties and exploring the impact of mass conservation.

Derivation of the PME by the Onsager principle

This section derives a model for gas flow in a homogeneous porous medium using the Onsager variational principle. The free energy of the gas is given by an integral over the domain, with a free energy density function f(ρ) that depends on the gas mass density ρ. Two approaches are presented to derive the Porous Medium Equation (PME) using the Onsager principle.

In the first approach, the Rayleighian functional is minimized with respect to the mass flux J, leading directly to the PME. In the second approach, the Rayleighian is minimized with respect to the time derivative of density ∂tρ and the gas velocity v, subject to the mass conservation equation. This results in a mixed form of the PME, where the Lagrange multiplier λ represents the gas pressure.

The second approach also shows that the system exhibits an energy decay property, where the time derivative of the free energy is less than or equal to 0. Finally, the PME is rewritten in a closed form that includes boundary conditions, and the well-posedness of the equation is discussed.

A moving mesh finite element method

This paper presents a moving mesh finite element method for the porous media equation (PME) derived using the Onsager variational principle. The method is first presented in one dimension and then generalized to higher dimensions.

The key aspects of the method are:

Semi-discretization: The domain is partitioned into N+1 time-dependent knots, and the finite element space Vh is used to approximate the density function ρ. This leads to a system of differential-algebraic equations for the nodal values of ρ and the mesh positions x.
Discrete energy and dissipation: The discrete energy functional Eh and dissipation functional Φh are defined. It is shown that the time derivative of Eh is non-positive, indicating the energy stability of the method.
Full discretization: Explicit and implicit time discretization schemes are derived. The implicit scheme is proven to be energy stable.
Implementation: The fully discrete schemes are formulated as decoupled linear systems that can be efficiently solved. The implicit scheme requires a fixed-point iteration.

The derivation and results can be generalized to higher dimensional cases, as demonstrated in the appendix and numerical examples in 2D.

Modified numerical schemes

The provided text discusses an alternative semi-discrete scheme to preserve mass conservation in the numerical simulation of a free boundary problem. The key points are:

The previous semi-discrete scheme (30) did not guarantee mass conservation due to the homogeneous boundary condition on the Lagrange multiplier.
To address this, the authors propose a new semi-discrete scheme that uses a piecewise linear approximation for the density and allows the grid nodes to move in a Lagrangian manner. This avoids imposing a boundary condition on the Lagrange multiplier.
The new variational problem and corresponding Euler-Lagrange equations are derived for this semi-discrete scheme. An algebraic formulation is also provided.
It is shown that the total mass is conserved in the semi-discrete scheme, in contrast to the previous one.
Explicit and implicit time discretization schemes are also developed, but the author notes that these fully discrete schemes do not satisfy the mass conservation property. Instead, the simpler schemes from the previous section are recommended for the numerical examples.

Numerical examples

This section presents numerical results to test the accuracy and effectiveness of the proposed numerical methods for solving the porous medium equation (PME) in one and two dimensions.

For the one-dimensional case, the authors use the Barenblatt-Pattle solution as the test problem. They find that the implicit numerical scheme achieves second-order convergence on a uniform mesh when the exponent m=2, but shows slower convergence when m=5 due to the more singular solution. Using a non-uniform initial mesh generated by a least squares method improves the convergence rate for both cases.

The authors also demonstrate that the numerical method can accurately capture the waiting time phenomenon, where the support of the solution does not change until a critical time.

For the two-dimensional case, similar convergence results are observed - optimal rates are obtained on non-uniform initial meshes, especially for larger m. The authors show examples of the two-dimensional Barenblatt-Pattle solution as well as more general initial data, including a case where two separate solution regions merge together over time.

Overall, the numerical results indicate the proposed methods can effectively solve the PME in one and two dimensions, accurately capturing key features of the solutions such as free boundaries and topology changes.

Conclusion

This paper presents a new moving mesh method for solving the porous medium equation, which is derived using the Onsager principle. The method can be applied to both the continuous and semi-discrete formulations of the porous medium equation, ensuring that the scheme maintains the same energy dissipation structure as the continuous problem.

The paper introduces a fully discrete explicit decoupled scheme and an implicit scheme. Numerical examples demonstrate the effectiveness of both schemes, showing that optimal convergence rates for the L2 error can be achieved when the initial meshes are appropriately selected.

The method is able to naturally capture the waiting time phenomena and can be extended to higher-dimensional problems and higher-order approximations. While the derivation of the method is intuitive, the error estimate for the moving mesh method remains an open question. Optimal convergence estimates will require suitable assumptions regarding the initial meshes and the regularity properties of the system.

Appendix A Numerical scheme in two dimensions

This section presents the discrete numerical scheme for the two-dimensional case of the Porous Medium Equation (PME). The key steps are:

Define a polygonal domain Ω_h and a finite element space V_h^t on a mesh T_h. The basis functions in V_h^t are piecewise linear.
Express the approximate solution ρ_h as a linear combination of the basis functions. This allows computing the time and spatial derivatives of ρ_h.
Derive the discrete Rayleighian functional R̃_h that combines the discrete energy functional E_h, the dissipation Φ_h, and the constraints from the continuum equation.
Compute the Euler-Lagrange equations by taking derivatives of R̃_h with respect to the unknowns (ρ, v_x, v_y, λ). This results in a system of algebraic equations.
The discrete system can be further discretized in time using an explicit or implicit scheme, though the details are omitted for simplicity.

The overall approach is to derive a finite element discretization of the PME, expressed in terms of the solution variables and their derivatives, and then obtain the corresponding Euler-Lagrange equations that define the discrete scheme.

This summary was produced with help from an AI and may contain inaccuracies - check out the links to read the original source documents!

Follow @aimodelsfyi on 𝕏 →

Related Papers

🔮

A Moving Mesh Method for Porous Medium Equation by the Onsager Variational Principle

Si Xiao, Xianmin Xu

In this paper, we introduce a new approach to solving the porous medium equation using a moving mesh finite element method that leverages the Onsager variational principle as an approximation tool. Both the continuous and discrete problems are formulated based on the Onsager principle. The energy dissipation structure is maintained in the semi-discrete and fully implicit discrete schemes. We also develop a fully decoupled explicit scheme by which only a few linear equations are solved sequentially in each time step. The numerical schemes exhibit an optimal convergence rate when the initial mesh is appropriately selected to ensure accurate approximation of the initial data. Furthermore, the method naturally captures the waiting time phenomena without requiring any manual intervention.

4/1/2024

Portable, Massively Parallel Implementation of a Material Point Method for Compressible Flows

Paolo Joseph Baioni, Tommaso Benacchio, Luigi Capone, Carlo de Falco

The recent evolution of software and hardware technologies is leading to a renewed computational interest in Particle-In-Cell (PIC) methods such as the Material Point Method (MPM). Indeed, provided some critical aspects are properly handled, PIC methods can be cast in formulations suitable to the requirements of data locality and fine-grained parallelism of modern hardware accelerators as Graphics Processing Units (GPUs). Such a rapid and continuous technological development increases also the importance of generic and portable implementations. While continuum mechanics simulations have already shown the capabilities of MPM on a wide range of phenomena, the use of the method in compressible fluid dynamics is less frequent, especially in the supersonic regime. In this paper we present a portable, highly parallel, GPU based MPM solver for compressible gas dynamics. The implementation aims to reach a good compromise between portability and efficiency and to give a first assessment of the potential of this approach in reproducing high speed gas flows, also taking into account solid obstacles. The proposed model constitutes a new step towards the realization of a monolithic MPM solver for Fluid-Structure Interaction (FSI) problems at all Mach numbers up to the supersonic regime.

4/29/2024

📊

A finite element-based physics-informed operator learning framework for spatiotemporal partial differential equations on arbitrary domains

Yusuke Yamazaki, Ali Harandi, Mayu Muramatsu, Alexandre Viardin, Markus Apel, Tim Brepols, Stefanie Reese, Shahed Rezaei

We propose a novel finite element-based physics-informed operator learning framework that allows for predicting spatiotemporal dynamics governed by partial differential equations (PDEs). The proposed framework employs a loss function inspired by the finite element method (FEM) with the implicit Euler time integration scheme. A transient thermal conduction problem is considered to benchmark the performance. The proposed operator learning framework takes a temperature field at the current time step as input and predicts a temperature field at the next time step. The Galerkin discretized weak formulation of the heat equation is employed to incorporate physics into the loss function, which is coined finite operator learning (FOL). Upon training, the networks successfully predict the temperature evolution over time for any initial temperature field at high accuracy compared to the FEM solution. The framework is also confirmed to be applicable to a heterogeneous thermal conductivity and arbitrary geometry. The advantages of FOL can be summarized as follows: First, the training is performed in an unsupervised manner, avoiding the need for a large data set prepared from costly simulations or experiments. Instead, random temperature patterns generated by the Gaussian random process and the Fourier series, combined with constant temperature fields, are used as training data to cover possible temperature cases. Second, shape functions and backward difference approximation are exploited for the domain discretization, resulting in a purely algebraic equation. This enhances training efficiency, as one avoids time-consuming automatic differentiation when optimizing weights and biases while accepting possible discretization errors. Finally, thanks to the interpolation power of FEM, any arbitrary geometry can be handled with FOL, which is crucial to addressing various engineering application scenarios.

8/7/2024

🧠

Finite Operator Learning: Bridging Neural Operators and Numerical Methods for Efficient Parametric Solution and Optimization of PDEs

Shahed Rezaei, Reza Najian Asl, Kianoosh Taghikhani, Ahmad Moeineddin, Michael Kaliske, Markus Apel

We introduce a method that combines neural operators, physics-informed machine learning, and standard numerical methods for solving PDEs. The proposed approach extends each of the aforementioned methods and unifies them within a single framework. We can parametrically solve partial differential equations in a data-free manner and provide accurate sensitivities, meaning the derivatives of the solution space with respect to the design space. These capabilities enable gradient-based optimization without the typical sensitivity analysis costs, unlike adjoint methods that scale directly with the number of response functions. Our Finite Operator Learning (FOL) approach uses an uncomplicated feed-forward neural network model to directly map the discrete design space (i.e. parametric input space) to the discrete solution space (i.e. finite number of sensor points in the arbitrary shape domain) ensuring compliance with physical laws by designing them into loss functions. The discretized governing equations, as well as the design and solution spaces, can be derived from any well-established numerical techniques. In this work, we employ the Finite Element Method (FEM) to approximate fields and their spatial derivatives. Subsequently, we conduct Sobolev training to minimize a multi-objective loss function, which includes the discretized weak form of the energy functional, boundary conditions violations, and the stationarity of the residuals with respect to the design variables. Our study focuses on the steady-state heat equation within heterogeneous materials that exhibits significant phase contrast and possibly temperature-dependent conductivity. The network's tangent matrix is directly used for gradient-based optimization to improve the microstructure's heat transfer characteristics. ...

7/8/2024