Related papers: Casimir preserving numerical method for global mul…
A proof of optimal-order error estimates is given for the full discretization of the bulk--surface Cahn--Hilliard system with dynamic boundary conditions in a smooth domain. The numerical method combines a linear bulk--surface finite…
Wave propagation problems for heterogeneous media are known to have many applications in physics and engineering. Recently, there has been an increasing interest in stochastic effects due to the uncertainty, which may arise from impurities…
This thesis deals with the investigation of a H(div)-conforming hybrid discontinuous Galerkin discretization for incompressible turbulent flows. The discretization method provides many physical and solving-oriented properties, which may be…
In this work, we develop a novel numerical scheme to solve the classical Keller--Segel (KS) model which simultaneously preserves its intrinsic mathematical structure and achieves optimal accuracy. The model is reformulated into a gradient…
This paper focuses on the numerical approximation of the linearized shallow water equations using hybridizable discontinuous Galerkin (HDG) methods, leveraging the Hamiltonian structure of the evolution system. First, we propose an…
This paper concerns preservation of velocity and pressure equilibria in smooth, compressible, multicomponent flows in the inviscid limit. First, we derive the velocity-equilibrium and pressure-equilibrium conditions of a standard…
A simple Hamiltonian modeling framework for general models in nonlinear optics is given. This framework is specialized to describe the Hamiltonian structure of electromagnetic phenomena in cubicly nonlinear optical media. The model has a…
In this paper, we present a new computational framework to approximate a Cahn-Hilliard-Navier-Stokes model with variable density and degenerate mobility that preserves the mass of the mixture, the pointwise bounds of the density and the…
We present a novel method for guaranteeing linear momentum in learned physics simulations. Unlike existing methods, we enforce conservation of momentum with a hard constraint, which we realize via antisymmetrical continuous convolutional…
We propose a physics-constrained machine learning method-based on reservoir computing- to time-accurately predict extreme events and long-term velocity statistics in a model of turbulent shear flow. The method leverages the strengths of two…
We develop several efficient numerical schemes which preserve exactly the global constraints for constrained gradient flows. Our schemes are based on the SAV approach combined with the Lagrangian multiplier approach. They are as efficient…
Simulating infiltration in porous media using Richards' equation remains computationally challenging due to its parabolic structure and nonlinear coefficients. While a wide range of numerical methods for differential equations have been…
This paper focuses on the numerical simulation of geothermal systems in complex geological settings. The physical model is based on two-phase Darcy flows coupling the mass conservation of the water component with the energy conservation and…
This work proposes a hyper-reduction method for nonlinear parametric dynamical systems characterized by gradient fields such as Hamiltonian systems and gradient flows. The gradient structure is associated with conservation of invariants or…
A Lagrangian numerical scheme for solving nonlinear degenerate Fokker-Planck equations in space dimensions $d\ge2$ is presented. It applies to a large class of nonlinear diffusion equations, whose dynamics are driven by internal energies…
The dynamics of an ideal fluid or plasma is constrained by topological invariants such as the circulation of (canonical) momentum or, equivalently, the flux of the vorticity or magnetic fields. In the Hamiltonian formalism, topological…
We formulate a novel numerical method suitable for the solution of topology optimization problems in solid mechanics. The most salient feature of the new approach is that the space and time discrete equations of the numerical method can be…
The Navier-Stokes-Korteweg (NSK) equations are a classical diffuse-interface model for compressible two-phase flow. As direct numerical simulations based on the NSK system are quite expensive and in some cases even impossible, we consider a…
The combined effectiveness of model reduction and the quasilinear approximation for the reproduction of the low-order statistics of oceanic surface boundary-layer turbulence is investigated. Idealized horizontally homogeneous problems of…
A nonhydrostatic dynamical core has been developed by using the multi-moment finite volume method that ensures the rigorous numerical conservation. To represent the spherical geometry free of polar problems, the cubed-sphere grid is…