The newly developed unifying discontinuous formulation named the correction procedure via reconstruction (CPR) for conservation laws is extended to solve the Navier-Stokes equations for 3D mixed grids. In the current development, tetrahedrons and triangular prisms are considered. The CPR method can unify several popular high order methods including the discontinuous Galerkin and the spectral volume methods into a more efficient differential form....

The goal of this paper is to obtain a well-balanced, stable, fast, and robust HLLC-type approximate Riemann solver for a hyperbolic nonconservative PDE system arising in a turbidity current model. The main difficulties come from the nonconservative nature of the system. A general strategy to derive simple approximate Riemann solvers for nonconservative systems is introduced, which is applied to the turbidity current model to obtain two different HLLC solvers. Some results concerning the non-negativity...

The goal of this paper is to obtain a well-balanced, stable, fast, and robust HLLC-type approximate Riemann solver for a hyperbolic nonconservative PDE system arising in a turbidity current model. The main difficulties come from the nonconservative nature of the system. A general strategy to derive simple approximate Riemann solvers for nonconservative systems is introduced, which is applied to the turbidity current model to obtain two different...

The present paper is devoted to the computation of single phase or two phase flows using the single-fluid approach. Governing equations rely on Euler equations which may be supplemented by conservation laws for mass species. Emphasis is given on numerical modelling with help of Godunov scheme or an approximate form of Godunov scheme called VFRoe-ncv based on velocity and pressure variables. Three distinct classes of closure laws to express the internal energy in terms of pressure, density and additional...

The present paper is devoted to the computation of single phase or two phase flows using the single-fluid approach. Governing equations rely on Euler equations which may be supplemented by conservation laws for mass species. Emphasis is given on numerical modelling with help of Godunov scheme or an approximate form of Godunov scheme called VFRoe-ncv based on velocity and pressure variables. Three distinct classes of closure laws to express the internal energy in terms of pressure, density...

In this paper we present a numerical study of the hyperbolic model for convection-diffusion transport problems that has been recently proposed by the authors. This model avoids the infinite speed paradox, inherent to the standard parabolic model and introduces a new parameter called relaxation time. This parameter plays the role of an “inertia” for the movement of the pollutant. The analysis presented herein is twofold: first, we perform an accurate study of the 1D steady-state equations and its...

In this paper we deal with a semilinear hyperbolic chemotaxis model in one space dimension evolving on a network, with suitable transmission conditions at nodes. This framework is motivated by tissue-engineering scaffolds used for improving wound healing. We introduce a numerical scheme, which guarantees global mass densities conservation. Moreover our scheme is able to yield a correct approximation of the effects of the source term at equilibrium. Several numerical tests are presented to show the...

This paper examines a linear extrapolation time-discretization of a $2D$ phase relaxation model with temperature dependent convection and reaction. The model consists of a diffusion-advection PDE for temperature and an ODE with double obstacle $\pm 1$ for phase variable. Under a stability constraint, this scheme is shown to converge with optimal orders $\mathcal{O}\left(\tau {\left|\log \tau \right|}^{1/2}\right)$ for temperature and enthalpy, and $\mathcal{O}\left({\tau }^{1/2}{\left|\log \tau \right|}^{1/2}\right)$ for heat flux as time-step $\tau \downarrow 0$.

In this work we study a fully discrete mixed scheme, based on continuous finite elements in space and a linear semi-implicit first-order integration in time, approximating an Ericksen–Leslie nematic liquid crystal model by means of a Ginzburg–Landau penalized problem. Conditional stability of this scheme is proved via a discrete version of the energy law satisfied by the continuous problem, and conditional convergence towards generalized Young measure-valued solutions to the Ericksen–Leslie problem...

This paper proposes a linear discrete-time scheme for general nonlinear cross-diffusion systems. The scheme can be regarded as an extension of a linear scheme based on the nonlinear Chernoff formula for the degenerate parabolic equations, which proposed by Berger et al. [RAIRO Anal. Numer.13 (1979) 297–312]. We analyze stability and convergence of the linear scheme. To this end, we apply the theory of reaction-diffusion system approximation. After discretizing the scheme in space, we obtain a versatile,...

This paper proposes a linear discrete-time scheme for general nonlinear cross-diffusion systems. The scheme can be regarded as an extension of a linear scheme based on the nonlinear Chernoff formula for the degenerate parabolic equations, which proposed by Berger et al. [RAIRO Anal. Numer.13 (1979) 297–312]. We analyze stability and convergence of the linear scheme. To this end, we apply the theory of reaction-diffusion system approximation. After discretizing the scheme in space, we obtain a versatile,...

An extension of the local projection stabilization (LPS) finite element method for convection-diffusion-reaction equations is presented and analyzed, both in the steady-state and the transient setting. In addition to the standard LPS method, a nonlinear crosswind diffusion term is introduced that accounts for the reduction of spurious oscillations. The existence of a solution can be proved and, depending on the choice of the stabilization parameter, also its uniqueness. Error estimates are derived...

There is evidence that cancer develops when cells acquire a sequence of mutations that alter normal cell characteristics. This sequence determines a hierarchy among the cells, based on how many more mutations they need to accumulate in order to become cancerous. When cells divide, they exhibit telomere loss and differentiate, which defines another cell hierarchy, on top of which is the stem cell. We propose a mutation-generation model, which combines...

In this paper we construct analytic-numerical solutions for initial-boundary value systems related to the equation $u_t-A{u}_{xx}-Bu=0$, where $B$ is an arbitrary square complex matrix and $A$ ia s matrix such that the real part of the eigenvalues of the matrix $\frac{1}{2}(A+A^H)$ is positive. Given an admissible error $\epsilon$ and a finite domain $G$, and analytic-numerical solution whose error is uniformly upper bounded by $\epsilon$ in $G$, is constructed.