-
A consistent and conservative model and its scheme for N-phase-M-component incompressible flows J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-25 Ziyang Huang; Guang Lin; Arezoo M. Ardekani
In the present work, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including
-
Horizontal approximation schemes for the staggered reduced latitude-longitude grid J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-25 Gordey S. Goyman; Vladimir V. Shashkin
In this work, we consider discretization of the linearized shallow water equations on a reduced latitude-longitude grid with an analogue of Arakawa C-type variables staggering. The resulting schemes are based on the use of longitudinal interpolation procedures and can be of arbitrary order of accuracy. We also present the analysis of conservation and wave propagation properties of these schemes for
-
A Method for Representing Periodic Functions and Enforcing Exactly Periodic Boundary Conditions with Deep Neural Networks J. Comput. Phys. (IF 2.985) Pub Date : 2021-03-01 Suchuan Dong; Naxian Ni
We present a simple and effective method for representing periodic functions and enforcing exactly the periodic boundary conditions for solving differential equations with deep neural networks (DNN). The method stems from some simple properties about function compositions involving periodic functions. It essentially composes a DNN-represented arbitrary function with a set of independent periodic functions
-
Fixed Inducing Points Online Bayesian Calibration for Computer Models with an Application to a Scale-Resolving CFD Simulation J. Comput. Phys. (IF 2.985) Pub Date : 2021-03-01 Yu Duan; Matthew D. Eaton; Michael J. Bluck
This paper proposes a novel fixed inducing points online Bayesian calibration (FIPO-BC) algorithm to efficiently learn the model parameters using a benchmark database. The standard Bayesian calibration (STD-BC) algorithm provides a statistical method to calibrate the parameters of computationally expensive models. However, the STD-BC algorithm does not scale well with regard to the number of data points
-
A locally implicit time-reversible sonic point processing algorithm for one-dimensional shallow-water equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-23 Nikita Afanasiev; Vasily Goloviznin
In this paper, we introduce a new locally implicit sonic point processing algorithm for one-dimensional shallow-water equations. The algorithm is based on transferring the shallow-water invariants along the characteristics around the sonic point, and it can be coupled with conservative-characteristic methods, which have problems modelling the transonic flows. For each sonic point, a system of two non-linear
-
A probabilistic generative model for semi-supervised training of coarse-grained surrogates and enforcing physical constraints through virtual observables J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-22 Maximilian Rixner; Phaedon-Stelios Koutsourelakis
The data-centric construction of inexpensive surrogates for fine-grained, physical models has been at the forefront of computational physics due to its significant utility in many-query tasks such as uncertainty quantification. Recent efforts have taken advantage of the enabling technologies from the field of machine learning (e.g., deep neural networks) in combination with simulation data. While such
-
Very high-order Cartesian-grid finite difference method on arbitrary geometries J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-23 S. Clain; D. Lopes; R.M.S. Pereira
An arbitrary order finite difference method for curved boundary domains with Cartesian grid is proposed. The technique handles in a universal manner Dirichlet, Neumann or Robin conditions. We introduce the Reconstruction Off-site Data (ROD) method, that transfers in polynomial functions the information located on the physical boundary to the computational domain. Three major advantages are: (1) a simple
-
Parallel load-balancing for combustion with spray for large-scale simulation J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 A. Thari; N.C.W. Treleaven; M. Staufer; G.J. Page
An asynchronous task-based Eulerian-Lagrangian approach for efficient parallel multi-physics simulations that can scale for arbitrary large number of particles and non-uniformly distributed particles is presented. The parallel methodology is based on a task-based partitioning of the multi-physics problem, where each single-physics problem is considered as a task and carried out using its own set of
-
A parallel-in-time two-sided preconditioning for all-at-once system from a non-local evolutionary equation with weakly singular kernel J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-23 Xue-lei Lin; Michael K. Ng; Yajing Zhi
In this paper, we study a parallel-in-time (PinT) algorithm for all-at-once system from a non-local evolutionary equation with weakly singular kernel where the temporal term involves a non-local convolution with a weakly singular kernel and the spatial term is the usual Laplacian operator with variable coefficients. Such a problem has been intensively studied in recent years thanks to the numerous
-
High order pressure-based semi-implicit IMEX schemes for the 3D Navier-Stokes equations at all Mach numbers J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-26 Walter Boscheri; Lorenzo Pareschi
This article aims at developing a high order pressure-based solver for the solution of the 3D compressible Navier-Stokes system at all Mach numbers. We propose a cell-centered discretization of the governing equations that splits the fluxes into a fast and a slow scale part, that are treated implicitly and explicitly, respectively. A novel semi-implicit discretization is proposed for the kinetic energy
-
An improved Immersed Boundary Method for turbulent flow simulations on Cartesian grids J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-26 Benjamin Constant; Stéphanie Péron; Héloïse Beaugendre; Christophe Benoit
In this paper, we present recent improvements of an Immersed Boundary Method (IBM) for the simulation of turbulent compressible flows on Cartesian grids. The proposed approach enables to remove spurious oscillations at the wall on skin pressure and friction coefficients. Results are compared to a body-fitted approach using the same wall function, showing that the stair-step immersed boundary provides
-
Particle integrator for particle-in-cell simulations of ultra-high intensity laser-plasma interactions J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-26 Kavin Tangtartharakul; Guangye Chen; Alexey Arefiev
Particle-in-cell codes are the most widely used simulation tools for kinetic studies of ultra-intense laser-plasma interactions. Using the motion of a single electron in a plane electromagnetic wave as a benchmark problem, we show surprising deterioration of the numerical accuracy of the PIC algorithm with increasing normalized wave amplitude for typical time-step and grid sizes. Two significant sources
-
Structure Probing Neural Network Deflation J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-26 Yiqi Gu; Chunmei Wang; Haizhao Yang
Deep learning is a powerful tool for solving nonlinear differential equations, but usually, only the solution corresponding to the flattest local minimizer can be found due to the implicit regularization of stochastic gradient descent. This paper proposes a network-based structure probing deflation method to make deep learning capable of identifying multiple solutions that are ubiquitous and important
-
Hermite spectral method for Fokker-Planck-Landau equation modeling collisional plasma J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-26 Ruo Li; Yinuo Ren; Yanli Wang
We propose an Hermite spectral method for the Fokker-Planck-Landau (FPL) equation. Both the distribution functions and the collision terms are approximated by series expansions of the Hermite functions. To handle the complexity of the quadratic FPL collision operator, a reduced collision model is built by adopting the quadratic collision operator for the lower-order terms and the diffusive Fokker-Planck
-
Quantum Fokker-Planck Modeling of Degenerate Electrons J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-25 Hai P. Le
An implicit and conservative numerical scheme is proposed for the isotropic quantum Fokker-Planck equation describing the evolution of degenerate electrons subject to elastic collisions with other electrons and ions. The electron-ion and electron-electron collision operators are discretized using a discontinuous Galerkin method, and the electron energy distribution is updated by an implicit time integration
-
Energy-consistent finite difference schemes for compressible hydrodynamics and magnetohydrodynamics using nonlinear filtering J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-25 Haruhisa Iijima
In this paper, an energy-consistent finite difference scheme for the compressible hydrodynamic and magnetohydrodynamic (MHD) equations is introduced. For the compressible magnetohydrodynamics, an energy-consistent finite difference formulation is derived using the product rule for the spatial difference. The conservation properties of the internal, kinetic, and magnetic energy equations can be satisfied
-
A comparison of the finite difference and multiresolution method for the elliptic equations with Dirichlet boundary conditions on irregular domains J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Ping Yin; Jacques Liandrat; Wanqiang Shen
We make a comparison of the finite difference and multiresolution method for solving the elliptic equations on irregular domains. The Dirichlet boundary condition is treated by the ghost fluid method (GFM) for the finite difference method and the Lagrange multiplier for the multiresolution method. Numerical results illustrate the improved convergence rate of errors and their gradients with the multiresolution
-
A mathematical model for thermal single-phase flow and reactive transport in fractured porous media J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-18 Alessio Fumagalli; Anna Scotti
In this paper we present a mathematical model and a numerical workflow for the simulation of a thermal single-phase flow with reactive transport in porous media, in the presence of fractures. The latter are thin regions which might behave as high or low permeability channels depending on their physical parameters, and are thus of paramount importance in underground flow problems. Chemical reactions
-
An implicit HDG method for linear convection-diffusion with dual time stepping J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Ruben Sevilla
This work presents, for the first time, a dual time stepping (DTS) approach to solve the global system of equations that appears in the hybridisable discontinuous Galerkin (HDG) formulation of convection-diffusion problems. A proof of the existence and uniqueness of the steady state solution of the HDG global problem with DTS is presented. The stability limit of the DTS approach is derived using a
-
Peridynamics enabled learning partial differential equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Ali C. Bekar; Erdogan Madenci
This study presents an approach to discover the significant terms in partial differential equations (PDEs) that describe particular phenomena based on the measured data. The relationship between the known field data and its continuous representation of PDEs is achieved through a linear regression model. It specifically employs the peridynamic differential operator (PDDO) and sparse linear regression
-
Data-driven discovery of coarse-grained equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Joseph Bakarji; Daniel M. Tartakovsky
Statistical (machine learning) tools for equation discovery require large amounts of data that are typically computer generated rather than experimentally observed. Multiscale modeling and stochastic simulations are two areas where learning on simulated data can lead to such discovery. In both, the data are generated with a reliable but impractical model, e.g., molecular dynamics simulations, while
-
Deep learning of the spanwise-averaged Navier–Stokes equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Bernat Font; Gabriel D. Weymouth; Vinh-Tan Nguyen; Owen R. Tutty
Simulations of turbulent fluid flow around long cylindrical structures are computationally expensive because of the vast range of length scales, requiring simplifications such as dimensional reduction. Current dimensionality reduction techniques such as strip-theory and depth-averaged methods do not take into account the natural flow dissipation mechanism inherent in the small-scale three-dimensional
-
Implicit reduced Vlasov–Fokker–Planck–Maxwell model based on high-order mixed elements J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Jan Nikl; Ilja Göthel; Milan Kuchařík; Stefan Weber; Michael Bussmann
Detailed description of the transport processes in plasma is crucial for many disciplines. When the mean-free-path of the electrons is comparable or exceeds a characteristic length scale of the plasma profile, non-local behavior can be observed. Predictions of the diffusion theory are not valid and non-local electric and magnetic fields are generated. Kinetic modeling of these phenomena on time scales
-
Calibration of projection-based reduced-order models for unsteady compressible flows J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-17 Victor Zucatti; William Wolf; Michel Bergmann
An analysis of calibration for reduced-order models (ROMs) is presented in this work. The Galerkin and least-squares Petrov-Galerkin (LSPG) methods are tested on compressible flows involving a disparity of temporal scales. A novel calibration strategy is proposed for the LSPG method and two test cases are analyzed. The first consists of a subsonic airfoil flow where boundary layer instabilities are
-
Development of an algebraic fractional step scheme for the primitive formulation of the compressible Navier-Stokes equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-11 Samuel Parada; Ramon Codina; Joan Baiges
In this work we address the compressible Navier-Stokes equations written in the so-called primitive formulation. The proposed methodology is a finite-element solver based on a fractional step scheme in time, which allows to uncouple the calculation of the problem unknowns providing important savings in computational cost. In addition, we include a stabilization technique within the Variational Multi-Scale
-
Nonlinear weighting process in ghost-cell immersed boundary methods for compressible flow J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-17 Hanahchim Choung; Vignesh Saravanan; Soogab Lee; Haeseong Cho
Computational challenges arise for the immersed boundary method (IBM) when dealing with compressible flow, where discontinuous and smoothly varying flow regions appear near the immersed boundary. The conventional ghost-cell IBM provides inaccurate results for smoothly varying regions when a low-order interpolation is used, or it suffers from a numerical instability for the discontinuous flow when a
-
GINNs: Graph-Informed Neural Networks for multiscale physics J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-13 Eric J. Hall; Søren Taverniers; Markos A. Katsoulakis; Daniel M. Tartakovsky
We introduce the concept of a Graph-Informed Neural Network (GINN), a hybrid approach combining deep learning with probabilistic graphical models (PGMs) that acts as a surrogate for physics-based representations of multiscale and multiphysics systems. GINNs address the twin challenges of removing intrinsic computational bottlenecks in physics-based models and generating large data sets for estimating
-
Principles of a network element method J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-15 Julien Coatléven
We introduce a new variational numerical method which does not require any background mesh to compute the scheme coefficients. Replacing the mesh by a point cloud endowed with connectivity, which we call a discretization network, and following the virtual element framework we derive a consistent, coercive and stable numerical scheme. We illustrate the good behavior of the method on several numerical
-
Efficient mesh deformation using radial basis functions with a grouping-circular-based greedy algorithm J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-17 Hong Fang; He Zhang; Fanli Shan; Ming Tie; Xing Zhang; Jinghua Sun
A grouping-circular-based (GCB) greedy algorithm is proposed to improve the efficiency of mesh deformation. By incorporating the multigrid concept that the computational errors on the fine mesh can be approximated with those on the coarse mesh, this algorithm stochastically divides all boundary nodes into m groups and uses the locally maximum radial basis functions (RBF) interpolation error of the
-
A structure preserving numerical scheme for Fokker-Planck equations of neuron networks: Numerical analysis and exploration J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-15 Jingwei Hu; Jian-Guo Liu; Yantong Xie; Zhennan Zhou
In this work, we are concerned with the Fokker-Planck equations associated with the Nonlinear Noisy Leaky Integrate-and-Fire model for neuron networks. Due to the jump mechanism at the microscopic level, such Fokker-Planck equations are endowed with an unconventional structure: transporting the boundary flux to a specific interior point. While the equations exhibit diversified solutions from various
-
An adaptive high-order piecewise polynomial based sparse grid collocation method with applications J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-05 Zhanjing Tao; Yan Jiang; Yingda Cheng
This paper constructs adaptive sparse grid collocation method onto arbitrary order piecewise polynomial space. The sparse grid method is a popular technique for high dimensional problems, and the associated collocation method has been well studied in the literature. The contribution of this work is the introduction of a systematic framework for collocation onto high-order piecewise polynomial space
-
Third order positivity-preserving direct discontinuous Galerkin method with interface correction for chemotaxis Keller-Segel equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-17 Changxin Qiu; Qingyuan Liu; Jue Yan
In this article, we apply direct discontinuous Galerkin method with interface correction (DDGIC) to solve chemotaxis Keller-Segel equation and prove the quadratic polynomial solution satisfying positivity-preserving with third order of accuracy. We show DDGIC method can obtain optimal convergence order for the cell density variable without introducing extra auxiliary variables to approximate the gradient
-
Enhanced weakly-compressible MPS method for violent free-surface flows: Role of particle regularization techniques J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-19 Mojtaba Jandaghian; Abdelkader Krimi; Amir Reza Zarrati; Ahmad Shakibaeinia
-
A consistent and parallelized height function based scheme for applying contact angle to 3D volume-of-fluid simulations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-11 Tian-Yang Han; Jie Zhang; Hua Tan; Ming-Jiu Ni
A height function based numerical scheme has been developed to apply contact angle in the parallel computational framework of 3D volume-of-fluid method. Given that most of the reported height function methods are implemented in 2D, the differences between the implementations of 2D and 3D height function methods are first focused. Then, a novel scheme is proposed to preserve the accuracy and consistency
-
Optimal preconditioners on solving the Poisson equation with Neumann boundary conditions J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-10 Byungjoon Lee; Chohong Min
MILU preconditioner is well known [16], [3] to be the optimal choice among all the ILU-type preconditioners in solving the Poisson equation with Dirichlet boundary conditions. However, it is less known which is an optimal preconditioner in solving the Poisson equation with Neumann boundary conditions. The condition number of an unpreconditioned matrix is as large as O(h−2), where h is the step size
-
Solving inverse problems using conditional invertible neural networks J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-11 Govinda Anantha Padmanabha; Nicholas Zabaras
Inverse modeling for computing a high-dimensional spatially-varying property field from indirect sparse and noisy observations is a challenging problem. This is due to the complex physical system of interest often expressed in the form of multiscale PDEs, the high-dimensionality of the spatial property of interest, and the incomplete and noisy nature of observations. To address these challenges, we
-
MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-15 Florian Holderied; Stefan Possanner; Xin Wang
We present a STRUcture-Preserving HYbrid code - STRUPHY - for the simulation of magneto-hydrodynamic (MHD) waves interacting with a population of energetic particles far from thermal equilibrium (kinetic species). The implemented model features linear, ideal MHD equations in curved, three-dimensional space, coupled nonlinearly to the full-orbit Vlasov equations via a current coupling scheme. The algorithm
-
Wall-modeled large eddy simulation of turbulent heat transfer by the lattice Boltzmann method J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-11 Y. Kuwata; K. Suga
A novel implementation route of the wall-function method to the lattice Boltzmann method (LBM) is proposed to extend the applicability of the LBM for high Reynolds number turbulent heat transfer in complex geometries. The proposed immersed virtual wall method assumes the virtual wall layer beneath the wall which satisfies the slip wall conditions allowing the subsurface heat and fluid flows within
-
Efficient and accurate gradient flow methods for computing ground states of spinor Bose-Einstein condensates J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-05 Yongyong Cai; Wei Liu
An accurate and efficient normalized gradient flow (NGF) method is developed to compute the ground state of a spinor Bose-Einstein condensate (BEC) with an arbitrary positive integer spin, which is described as a vector wave function that minimizes the Gross-Pitaevskii energy functional under the two constraints of the total mass and magnetization. Similar to the NGF for computing the ground states
-
A high-order moment limiter for the discontinuous Galerkin method on triangular meshes J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-10 Krishna Dutt; Lilia Krivodonova
We propose a moment limiter of arbitrary high order for the discontinuous Galerkin method on unstructured triangular meshes. The limiter works by hierarchically limiting solution coefficients (moments) by comparing them to reconstructed directional derivatives along specific directions. Limiting along these directions is performed using one-dimensional minmod slope limiters. The stencil used to reconstruct
-
An arbitrary-order Cell Method with block-diagonal mass-matrices for the time-dependent 2D Maxwell equations J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-08 Bernard Kapidani; Lorenzo Codecasa; Joachim Schöberl
We introduce a new numerical method for the time-dependent Maxwell equations on unstructured meshes in two space dimensions. This relies on the introduction of a new mesh, which is the barycentric-dual cellular complex of the starting simplicial mesh, and on approximating two unknown fields with integral quantities on geometric entities of the two dual complexes. A careful choice of basis functions
-
An alternative extended linear system for boundary value problems on locally perturbed geometries J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-04 Yabin Zhang; Adrianna Gillman
This manuscript presents a new extended linear system for integral equation based techniques for solving boundary value problems on locally perturbed geometries. The new extended linear system is similar to a previously presented technique for which the authors have constructed a fast direct solver. The key features of the work presented in this paper are that the fast direct solver is more efficient
-
A friction interface model for multi-material interactions in a Eulerian framework J. Comput. Phys. (IF 2.985) Pub Date : 2020-12-10 Wanli Wang; Cheng Wang; Tonghui Yang; Dongping Chen
A friction model for hyper-elastic solid materials was proposed for the simulation of complex shear impacts in a Eulerian framework. The interfacial status of multi-material interactions was obtained in Harten, Lax, and van Leer discontinuity (HLLD) Riemann solver. The inverse deformation gradient tensor was utilized in governing equations to describe the shape change of hyper-elastic solid materials
-
A variational interface-preserving and conservative phase-field method for the surface tension effect in two-phase flows J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-04 Xiaoyu Mao; Vaibhav Joshi; Rajeev Jaiman
We present a finite element based variational interface-preserving and conservative phase-field formulation for the modeling of incompressible two-phase flows with surface tension dynamics. The preservation of the hyperbolic tangent interface profile of the convective Allen-Cahn phase-field formulation relies on a novel time-dependent mobility model. The mobility coefficient is adjusted adaptively
-
Space-time multilevel Monte Carlo methods and their application to cardiac electrophysiology J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-05 S. Ben Bader; P. Benedusi; A. Quaglino; P. Zulian; R. Krause
We present a novel approach aimed at high-performance uncertainty quantification for time-dependent problems governed by partial differential equations. In particular, we consider input uncertainties described by a Karhunen-Loève expansion and compute statistics of high-dimensional quantities-of-interest, such as the cardiac activation potential. Our methodology relies on a close integration of multilevel
-
A hybrid mimetic spectral element method for three-dimensional linear elasticity problems J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-04 Yi Zhang; Joël Fisser; Marc Gerritsma
We introduce a domain decomposition structure-preserving method based on a hybrid mimetic spectral element method for three-dimensional linear elasticity problems in curvilinear conforming structured meshes. The method is an equilibrium method which satisfies pointwise equilibrium of forces. The domain decomposition is established through hybridization which first allows for an inter-element normal
-
Multiscale and monolithic arbitrary Lagrangian–Eulerian finite element method for a hemodynamic fluid-structure interaction problem involving aneurysms J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-05 Wenrui Hao; Pengtao Sun; Jinchao Xu; Lian Zhang
In this paper, a multiscale and monolithic arbitrary Lagrangian–Eulerian finite element method (ALE-FEM) is developed for a multiscale hemodynamic fluid-structure interaction (FSI) problem involving an aortic aneurysm growth to quantitatively predict the long-term aneurysm risk in the cardiovascular environment, where the blood fluid profile, the hyperelastic arterial wall, and the aneurysm pathophysiology
-
Adaptive deformation of 3D unstructured meshes with curved body fitted boundaries with application to unsteady compressible flows J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-08 Luca Cirrottola; Mario Ricchiuto; Algiane Froehly; Barbara Re; Alberto Guardone; Giuseppe Quaranta
We present an adaptive moving mesh method for unstructured meshes which is a three-dimensional extension of the previous works of Ceniceros et al. [10], Tang et al. [40] and Chen et al. [11]. The iterative solution of a variable diffusion Laplacian model on the reference domain is used to adapt the mesh to moving sharp solution fronts while imposing slip conditions for the displacements on curved boundary
-
Stable Perfectly Matched Layers with Lorentz transformation for the convected Helmholtz equation J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-05 Philippe Marchner; Hadrien Bériot; Xavier Antoine; Christophe Geuzaine
Perfectly Matched Layers (PMLs) appear as a popular alternative to non-reflecting boundary conditions for wave-type problems. The core idea is to extend the computational domain by a fictitious layer with specific absorption properties such that the wave amplitude decays significantly and does not produce back reflections. In the context of convected acoustics, it is well-known that PMLs are exposed
-
A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-08 Xiaoxi Zhao; Birendra Jha
Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated
-
Customized data-driven RANS closures for bi-fidelity LES–RANS optimization J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-26 Yu Zhang; Richard P. Dwight; Martin Schmelzer; Javier F. Gómez; Zhong-hua Han; Stefan Hickel
Multi-fidelity optimization methods promise a high-fidelity optimum at a cost only slightly greater than a low-fidelity optimization. This promise is seldom achieved in practice, due to the requirement that low- and high-fidelity models correlate well. In this article, we propose an efficient bi-fidelity shape optimization method for turbulent fluid-flow applications with Large-Eddy Simulation (LES)
-
A comparative study of machine learning models for predicting the state of reactive mixing J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-19 B. Ahmmed; M.K. Mudunuru; S. Karra; S.C. James; V.V. Vesselinov
Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass growth. To investigate the evolution of mixing dynamics under different scenarios (e.g
-
Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-26 Wai Hong Ronald Chan; Michael S. Dodd; Perry L. Johnson; Parviz Moin
-
LANS-α turbulence modeling for coastal sea: An application to Alderney Race J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-26 A.-C. Bennis; F. Adong; M. Boutet; F. Dumas
The Lagrangian-Averaged Navier-Stokes-α (LANS-α) turbulence model was implemented for the first time in a coastal hydrodynamic model. We present in this paper the details of the implementation, as well as the difficulties encountered. To overcome the difficulties, a convolution filter was used instead of the Helmoltz operator, and incompressibility was imposed in both rough and smooth velocities. The
-
A stable and accurate scheme for solving the Stefan problem coupled with natural convection using the Immersed Boundary Smooth Extension method J. Comput. Phys. (IF 2.985) Pub Date : 2021-02-01 Jinzi Mac Huang; Michael J. Shelley; David B. Stein
The dissolution of solids has created spectacular geomorphologies ranging from centimeter-scale cave scallops to the kilometer-scale “stone forests” of China and Madagascar. Mathematically, dissolution processes are modeled by a Stefan problem, which describes how the motion of a phase-separating interface depends on local concentration gradients, coupled to a fluid flow. Simulating these problems
-
An adaptive Hessian approximated stochastic gradient MCMC method J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-21 Yating Wang; Wei Deng; Guang Lin
Bayesian approaches have been successfully integrated into training deep neural networks. One popular family is stochastic gradient Markov chain Monte Carlo methods (SG-MCMC), which have gained increasing interest due to their ability to handle large datasets and the potential to avoid overfitting. Although standard SG-MCMC methods have shown great performance in a variety of problems, they may be
-
An analysis of three formulations of the tensor artificial viscosity in two-dimensional Cartesian geometry J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-27 Zhiwei Lin; Shaoen Jiang; Lu Zhang; Longyu Kuang
This paper analyzes three formulations of the tensor artificial viscosity in two-dimensional Cartesian geometry, namely, Campbell-Shashkov viscosity, Lipnikov-Shashkov viscosity and Wendroff viscosity. We first present a general derivation of these three ones. The discrepancy between Campbell-Shashkov viscosity and Lipnikov-Shashkov viscosity is then provided, followed by the non-triviality proof of
-
Learning and correcting non-Gaussian model errors J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-29 Danny Smyl; Tyler N. Tallman; Jonathan A. Black; Andreas Hauptmann; Dong Liu
All discretized numerical models contain modeling errors – this reality is amplified when reduced-order models are used. The ability to accurately approximate modeling errors informs statistics on model confidence and improves quantitative results from frameworks using numerical models in prediction, tomography, and signal processing. Further to this, the compensation of highly nonlinear and non-Gaussian
-
Application of discrete mechanics model to jump conditions in two-phase flows J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-21 Jean-Paul Caltagirone
Discrete mechanics is presented as an alternative to the equations of fluid mechanics, in particular to the Navier-Stokes equation. The derivation of the discrete equation of motion is built from the intuitions of Galileo, the principles of Galilean equivalence and relativity. Other more recent concepts such as the equivalence between mass and energy and the Helmholtz-Hodge decomposition complete the
-
Phase-change modeling based on a novel conservative phase-field method J. Comput. Phys. (IF 2.985) Pub Date : 2021-01-13 Reza Haghani-Hassan-Abadi; Abbas Fakhari; Mohammad-Hassan Rahimian
Based on the Allen-Cahn equation, we propose a phase-field model for liquid-vapor phase-change phenomena. We first extend the conservative form of the Allen-Cahn equation to include phase-change effects, and then develop a lattice-Boltzmann model to numerically solve the governing equations for interface motion and hydrodynamics, both of which include mass transfer effects. The net heat flux at the
Contents have been reproduced by permission of the publishers.