In this paper, we develop an accelerated conservative sharp-interface method for multiphase flows simulations. Traditional multiphase simulation methods use the minimum time step of all fluids obtained according to the CFL conditions to evolve the fluid states, which limits the computational efficiency, as the sound speed of one fluid may be much larger than of others. To address this issue, based

In this paper we introduce constitutive artificial neural networks (CANNs), a novel machine learning architecture for data-driven modeling of the mechanical constitutive behavior of materials. CANNs are able to incorporate by their very design information from three different sources, namely stress-strain data, theoretical knowledge from materials theory, and diverse additional information (e.g., about

When parameter estimation is solved in a high-dimensional space, the dimensionality reduction strategy becomes the primary consideration for alleviating the tremendous computational cost. In the present study, the discrete empirical interpolation method (DEIM) is explored to retrieve the initial condition (IC) by combining the polynomial chaos (PC) based ensemble Kalman filter (i.e. PC-EnKF), where

Flux-corrected transport (FCT) is one of the flux limiter methods. Unlike the total variation diminishing methods, obtaining the known FCT formulas for computing flux limiters is not quite transparent, and their transformation is not obvious when the original differential operator changes. We propose a novel formal mathematical approach to design flux correction for weighted hybrid difference schemes

The Discontinuous Petrov-Galerkin (DPG) method and the exponential integrators are two well stablished numerical methods for solving Partial Differential Equations (PDEs) and stiff systems of Ordinary Differential Equations (ODEs), respectively. In this work, we apply the DPG method in the time variable for linear parabolic problems and we calculate the optimal test functions analytically. We show

We present a new intersection-distribution-based remapping method between arbitrary polygonal meshes for indirect staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics. All cell-centered material quantities are conservatively remapped using intersections between the Lagrangian (old, source) mesh and the rezoned (new, target) mesh. The new nodal masses are obtained by conservative distribution

In large-eddy simulations (LES) a computational-domain translation velocity can be used to improve performance by allowing longer time-step intervals. The continuous equations are Galilean invariant, however, standard finite-difference-based discretizations are not discretely invariant with the error being proportional to the product of the local translation velocity and the truncation error. Even

The influence of different boundary conditions on the spectral properties of the incompressible Navier-Stokes equations is investigated. By using the Fourier-Laplace transform technique, we determine the spectra, extract the decay rate in time, and investigate the dispersion relation. In contrast to an infinite domain, where only diffusion affects the convergence, we show that also the propagation

This review summarizes the rigorous mathematical theory behind the lattice Boltzmann equation (LBE). Relevant properties of the Boltzmann equation and a derivation of the LBE from the Boltzmann equation are presented. A summary of some important LBE models is provided. Focus is given to results from the numerical analysis of the LBE as a solver for the nearly incompressible Navier-Stokes equations

Dilatancy plays a key role in mixtures of grains and fluid but is poorly investigated in dry granular flows. These flows may however dilate by more than 10% in granular column collapses. We investigate here dilatancy effects in dry flows with a shallow depth-averaged model designed to be further applied to simulate natural landslides. We use a compressible μ(I), ϕeq(I) rheology with a dilatancy law

Current multi-component, multiphase pseudo-potential lattice Boltzmann models have thermodynamic inconsistencies that prevent them to correctly predict the thermodynamic phase behavior of partially miscible multi-component mixtures, such as hydrocarbon mixtures. This paper identifies these inconsistencies and attempts to design a thermodynamically consistent multi-component, multiphase pseudo-potential

Developing accurate, robust and efficient implicit time-marching method is crucial for the rapid convergence of high Reynolds number wall-bounded flows. In the presence that a massive structured grid is partitioned into multiple multiblocks for parallelization, a real-time implicit boundary condition treatment taking the multiblock communications into account affects the convergence acceleration dramatically

A numerical approach for the modeling and simulation of fluid-structure interaction (FSI) in multi-material and multi-phase systems with potential phase-changes dynamics is presented. The boundary conditions at the interface between the fluids and structures are enforced using an immersed boundary technique to couple the Eulerian multi-material solver to the Lagrangian structural solver and maintain

The numerical solutions to nonlinear hyperbolic balance laws at (or near) steady state may develop spurious oscillations due to the imbalance between flux and source terms. In the present article, we study a high order well-balanced discontinuous Galerkin (DG) scheme for balance law with subsonic flow, which preserves equilibrium solutions of the flow exactly, and also provides non-oscillatory solutions

Partial differential equations are frequently solved using a global basis, such as the Fourier series, due to excellent convergence. However, convergence becomes impaired when discontinuities are present due to the Gibbs phenomenon, negatively impacting simulation speed and possibly generating spurious solutions. We resolve this by supplementing the smooth global basis with an inherently discontinuous

In this paper, a novel semi-implicit lubrication dynamics method that can efficiently simulate dense non-colloidal suspensions is proposed. To reduce the computational cost in the presented methodology, inter-particle lubrication-based forces and torques alone are considered together with a short-range repulsion to enforce finite inter-particle separation due to surface roughness, Brownian forces or

A novel multiple-relaxation-time (MRT) lattice Boltzmann model is proposed for the radiative transfer equation (RTE). In this paper, the discussion and implementation are restricted to the grey (frequency-independent) radiative transfer equation. We establish this model by regarding the RTE as a particular convection-diffusion equation without the diffusion term. The equilibrium distribution function

Uncertainty is embedded in groundwater flow modeling because of the heterogeneity of hydraulic conductivity and the scarcity of measurements. To quantify the uncertainty of the modeled hydraulic head, this study proposes an improved version of the meshless generalized finite difference method (GFDM) for solving the statistical moment equation (ME). The proposed GFDM adopts a new support sub-domain

Domain Decomposition Methods (DDM) have been widely used in the Computational Electromagnetics (CEM) community in the last years to tackle large-scale problems with Finite Element Methods (FEM). Non-conformal DDM is more flexible (e.g., independently created meshes for different parts of the problem under analysis are supported) but may introduce an approximation error. In this communication, a thorough

In this paper, a low dissipation fifth-order finite difference nested multi-resolution weighted essentially non-oscillatory (WENO) scheme is presented for solving Euler/Navier-Stokes equations in multi-dimensions on structured meshes. In the reconstruction procedures of this finite difference WENO, a series of nested unequal-sized central stencils are used to design spatial reconstruction polynomials

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical

We present an immersed boundary projection method for flow-structure interaction problems involving rigid bodies with complex geometries. Dynamics of a rigid body interacting with fluid flow and kinematics of other rigid bodies undergoing prescribed motions are coupled implicitly with the incompressible vorticity equations. In particular, the method is formulated in a frame of reference fixed on the

Data-driven surrogate models are widely used for applications such as design optimization and uncertainty quantification, where repeated evaluations of an expensive simulator are required. For most partial differential equation (PDE) simulations, the outputs of interest are often spatial or spatial-temporal fields, leading to very high-dimensional outputs. Despite the success of existing data-driven

This paper presents mathematical formulations and methods to predict the effect of forced-flight variance reduction on Monte Carlo tally variance and calculation time. This includes deducing biasing operators that are then used to construct a history-score probability density function (HSPDF), which represents all possible Monte Carlo random walks and gives the probability of a Monte Carlo history

In this paper we study finite element method for three-dimensional incompressible resistive magnetohydrodynamic equations, in which the velocity, the current density, and the magnetic induction are divergence-free. It is desirable that the discrete solutions should also satisfy divergence-free conditions exactly especially for the momentum equations. Inspired by constrained transport method, we devise

We introduce DeepMoD, a Deep learning based Model Discovery algorithm. DeepMoD discovers the partial differential equation underlying a spatio-temporal data set using sparse regression on a library of possible functions and their derivatives. A neural network is used as function approximator and its output is used to constructs the function library, allowing to perform the sparse regression within

Exponential and IMEX integrators are considered for wave turbulence problems, which are stiff systems that often require long-time simulations. Applying a multi-scale asymptotic analysis to the difference equations for a simple wave turbulence problem reveals that Exponential Runge-Kutta methods incorrectly model interactions between resonant waves unless they are specifically designed otherwise, and

We consider the problem of electrowetting on dielectric. The system involves the dynamics of a conducting droplet, which is immersed in another dielectric fluid, on a dielectric substrate under an applied voltage. The fluid dynamics is modeled by the two-phase incompressible Navier-Stokes equations with the standard interface conditions, the Navier slip condition on the substrate and a contact angle

Drops on a free-flow/porous-medium-flow interface have a strong influence on the exchange of mass, momentum and energy between the two macroscopic flow regimes. Modeling droplet-related pore-scale processes in a macro-scale context is challenging due to the scale gap, but might be rewarding due to relatively low computational costs. We develop a three-domain approach to model drop formation, growth

We develop a nested hybridizable discontinuous Galerkin (HDG) method to numerically solve the Maxwell's equations coupled with a hydrodynamic model for the conduction-band electrons in metals. The HDG method leverages static condensation to eliminate the degrees of freedom of the approximate solution defined in the elements, yielding a linear system in terms of the degrees of freedom of the approximate

An improved immersed boundary method coupled to turbulence wall models on Cartesian grids is proposed, for producing smooth wall surface pressure and skin friction at high Reynolds numbers. Spurious oscillations are frequently observed on these quantities with most immersed boundary wall modeling methods, especially for the skin friction which is found to be very sensitive to the solid surface's position

In this work we develop a discretisation method for the mixed formulation of the magnetostatic problem supporting arbitrary orders and polyhedral meshes. The method is based on a global discrete de Rham (DDR) sequence, obtained by patching the local spaces constructed in [22] by enforcing the single-valuedness of the components attached to the boundary of each element. The first main contribution of

Knowing how the solution to time-harmonic wave scattering problems depends on medium properties and boundary conditions is pivotal in wave-based inverse problems, e.g. for imaging. This paper is devoted to the exposition of a computationally efficient method, called the adjoint state method, that allows to quantify the influence of media properties, directly and through boundary conditions, in the

An energetically balanced, implicit integrator for non-hydrostatic vertical atmospheric dynamics on the sphere is presented. The integrator allows for the exact balance of energy exchanges in space and time for vertical atmospheric motions by preserving the skew-symmetry of the non-canonical Hamiltonian formulation of the compressible Euler equations. The performance of the integrator is accelerated

A system of mass balance equations is set up for plasma-chemical simulations, we derive the diffusion velocities from Stefan-Maxwell theory. An algorithm to linearly transform the system of mass balance equations in the context of non-LTE plasmas is described. This transformation is derived from the reaction set, and eliminates part of the chemical source term, enforces invariants such as conservation

We present a method for the data-driven learning of physical phenomena whose evolution in time depends on history terms. It is well known that a Mori-Zwanzig-type projection produces a description of the physical phenomena that depends on history, and also incorporates noise. If the data stream is sampled from the projected Mori-Zwanzig manifold, the description of the phenomenon will always depend

The time fractional Klein–Kramers model (TFKKM) is obtained by incorporating the subdiffusive mechanisms into the Klein–Kramers formalism. The TFKKM can efficiently express subdiffusion while an external force field is present in the phase space. The model describes the escape of a particle over a barrier and has a significant role in examining a variety of systems including slow (subdiffusion) dynamics

This work introduces a numerical fixed point method to approximate the solutions of a Vlasov-Poisson-Boltzmann boundary value problem which arises when modeling a bi-species collisional sheath. Our method relies on the exact integration of the transport equations by means of the characteristic curves. A special care is given about the choice of a suitable phase space discretization together with the

In this paper, new Locally Linearized Runge-Kutta formulas of Dormand and Prince are introduced with the purpose of integrating large systems of initial value problems. The new formulas are obtained by replacing the Padé approximation for matrix exponential by a Krylov-Padé approximation in the embedded formulas proposed in a previous work. Unlike other high order exponential integrators, these new

This paper presents a machine-learning approach for determining the optimal anisotropy in a computational mesh, in the context of an output-based adaptive solution procedure. Artificial neural networks are used to predict the desired element aspect ratio from readily accessible features of the primal and adjoint solutions. Whereas the sizing of the element is still based on an adjoint-weighted residual

We consider the initial-boundary value Euler equations with the aim to derive boundary conditions that yield an entropy bound for the physical (Navier-Stokes) entropy. We begin by reviewing the entropy bound obtained for standard no-penetration wall boundary conditions and propose a numerical implementation. The main results is the derivation of full-state boundary conditions (far-field, inlet, outlet)

The recently developed quadrature by expansion (QBX) technique [24] accurately evaluates the layer potentials with singular, weakly or nearly singular, or even hyper singular kernels in the integral equation reformulations of partial differential equations. The idea is to form a local complex polynomial or partial wave expansion centered at a point away from the boundary to avoid the singularity in

A conservative sequential fully implicit method is derived for compositional reservoir simulation. Multi-phase flow in porous media comprises coupled complex processes: i.e. elliptic flow equation, hyperbolic transport equation and highly nonlinear phase equilibrium equation. These processes contain very different mathematical characteristics that cannot be efficiently solved by one numerical method

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover

Classical methods of calibration usually imply the minimisation of an objective function with respect to some control parameters. This function measures the error between some observations and the results obtained by a numerical model. In the presence of uncontrollable additional parameters that we model as random inputs, the objective function becomes a random variable, and notions of robustness have

A decoupled lattice Boltzmann method (LBM) is developed from the conventional multiple relaxation time (MRT) model by eliminating the second-order mutual effects of viscosity modes. The dependency between relaxation times and viscosity is decoupled to remove the O(u3) errors in the original model, and thus a stabilization strategy is proposed by adjusting relaxation times without changing the local

We present a multiphase level set method with local volume conservation for capillary-controlled displacement in porous structures. Standard numerical formulations of the level set method for capillary-controlled (or, curvature-driven) motions assume phase pressures and interface properties are spatially uniform and disregard the fact that separate phase ganglia typically have distinct pressures. This

We present a super-time-stepping scheme for numerically solving parabolic partial differential equations with Dirichlet boundary conditions (BC). Using the general Forward Euler scheme, one can show that by taking varying step sizes there is the potential of propagating the solution forward in time by a greater amount than with uniform step sizes, while maintaining the same order of accuracy. As shown

Over the past two decades, there has been much development in discontinuous Galerkin methods for incompressible flows and for compressible flows with a positive Mach number, but almost no attention has been paid to variable-density flows at low speeds. This paper presents a pressure-based discontinuous Galerkin method for flow in the low-Mach number limit. We use a variable-density pressure correction

In this paper, high order well-balanced finite difference weighted essentially non-oscillatory methods to solve general systems of balance laws are presented. Two different families are introduced: while the methods in the first one preserve every stationary solution, those in the second family only preserve a given set of stationary solutions that depend on some parameters. The accuracy, well-balancedness

This paper discusses energy-conserving time-discretizations for finite element particle-in-cell discretizations of the Vlasov–Maxwell system. A geometric spatially discrete system can be obtained using a standard particle-in-cell discretization of the particle distribution and compatible finite element spaces for the fields to discretize the Poisson bracket of the Vlasov–Maxwell model (see Kraus et

We aim at analyzing a novel structure preserving difference scheme for the high dimensional nonlinear space-fractional Schrödinger equation. The temporal direction is discretized by the modified Crank-Nicolson scheme, while the spacial variable is approximated by formulas from the shifted convolution quadrature. The energy and mass preserving properties are proved rigorously for the scheme based on

Unresolved gradients produce numerical oscillations and inaccurate results. The most straightforward solution to such a problem is to increase the resolution of the computational grid. However, this is often prohibitively expensive and may lead to ecessive execution times. By training a neural network to predict the shape of the solution, we show that it is possible to reduce numerical oscillations

Simulation of multiphase flows, which are ubiquitous in nature and engineering applications, require coupled capturing or tracking of the interfaces in conjunction with the solution of fluid motion often occurring at multiple scales. In this contribution, we will present unified cascaded LB methods based on central moments for the solution of the incompressible two-phase flows at high density ratios

Partitioned solutions to fluid-structure interaction problems often employ a Dirichlet-Neumann decomposition, where the fluid equations are solved subject to Dirichlet boundary conditions on velocity from the structure, and the structure equations are solved subject to forces from the fluid. In some scenarios, such as an elastic balloon filling with air, an incompressible fluid domain may have pure

This article presents a high order conservative flux optimization (CFO) finite element method that preserve the important mass conservation property locally for the convection-diffusion equations. The numerical scheme is an extension of the first order CFO scheme proposed in [25], it is based on the classical Galerkin finite element method enhanced by a flux approximation on the boundary of a prescribed

We propose an Eulerian approach to compute the expected finite time Lyapunov exponent (FTLE) of uncertain flow fields. The definition extends the usual FTLE for deterministic dynamical systems. Instead of performing Monte Carlo simulations as in typical Lagrangian computations, our approach associates each initial flow particle with a probability density function (PDF) which satisfies an advection-diffusion

We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with periodic

For numerical schemes to the incompressible Navier-Stokes equations with variable density, it is a critical property to preserve the bounds of density. A bound-preserving high order accurate scheme can be constructed by using high order discontinuous Galerkin (DG) methods or finite volume methods with a bound-preserving limiter for the density evolution equation, with any popular numerical method for