The paper deals with the numerical resolution of the convection-diffusion system which arises when modeling combustion for turbulent flow. The considered model is of compressible turbulent reacting type where the turbulence-chemistry interactions are governed by additional balance equations. The system of PDE’s, that governs such a model, turns out to be in non-conservation form and usual numerical approaches grossly fail in the capture of viscous shock layers. Put in other words, classical finite...

The paper deals with the numerical resolution of the convection-diffusion system which arises when modeling combustion for turbulent flow. The considered model is of compressible turbulent reacting type where the turbulence-chemistry interactions are governed by additional balance equations. The system of PDE's, that governs such a model, turns out to be in non-conservation form and usual numerical approaches grossly fail in the capture of viscous shock layers. Put in other words, classical finite...

A new scheme for discretizing the P1 model on unstructured polygonal meshes is proposed. This scheme is designed such that its limit in the diffusion regime is the MPFA-O scheme which is proved to be a consistent variant of the Breil-Maire diffusion scheme. Numerical tests compare this scheme with a derived GLACE scheme for the P1 system.

We present an asynchronous multi-domain time integration algorithm with a dual domain decomposition method for the initial boundary-value problems for a parabolic equation. For efficient parallel computing, we apply the three-field domain decomposition method with local Lagrange multipliers to ensure the continuity of the primary unknowns at the interface between subdomains. The implicit method for time discretization and the multi-domain spatial decomposition enable us to use different time steps...

We present the current Reduced Basis framework for the efficient numerical approximation of parametrized steady Navier–Stokes equations. We have extended the existing setting developed in the last decade (see e.g. [S. Deparis, SIAM J. Numer. Anal. 46 (2008) 2039–2067; A. Quarteroni and G. Rozza, Numer. Methods Partial Differ. Equ. 23 (2007) 923–948; K. Veroy and A.T. Patera, Int. J. Numer. Methods Fluids 47 (2005) 773–788]) to more general affine and nonaffine parametrizations (such as volume-based...

One of the most challenging problems in the optimal control theory consists of solving the nonsmooth optimal control problems where several discontinuities may be present in the control variable and derivative of the state variable. Recently some extended spectral collocation methods have been introduced for solving such problems, and a matrix of differentiation is usually used to discretize and to approximate the derivative of the state variable in the particular collocation points. In such methods,...

This paper deals with nonlinear diffusion problems which include the Stefan problem, the porous medium equation and cross-diffusion systems. We provide a linear scheme for these nonlinear diffusion problems. The proposed numerical scheme has many advantages. Namely, the implementation is very easy and the ensuing linear algebraic systems are symmetric, which show low computational cost. Moreover, this scheme has the accuracy comparable to that of the wellstudied nonlinear schemes and make it possible...

An energy analysis is carried out for the usual semidiscrete Galerkin method for the semilinear equation in the region$\Omega$ (E) ${\left(t{u}_t\right)}_t={\sum }_{i,j=1}{\left(a_{ij}\left(x\right)u_{x_i}\right)}_{x_j}-a_0\left(x\right)u+f\left(u\right)$, subject to the initial and boundary conditions, $u=0$ on $\partial \Omega$ and $u(x,0)=u_0$. (E) is degenerate at $t=0$ and thus, even in the case $f\equiv 0$, time derivatives of $u$ will blow up as $t\to 0$. Also, in the case where $f$ is locally Lipschitz, solutions of (E) can blow up for $t>0$ in finite time. Stability and convergence of the scheme in $W^{2,1}$ is shown in the linear case without assuming $u_{tt}$ (which can blow up as $t\to 0$ is...

We consider the system of partial differential equations governing the one-dimensional flow of two superposed immiscible layers of shallow water. The difficulty in this system comes from the coupling terms involving some derivatives of the unknowns that make the system nonconservative, and eventually nonhyperbolic. Due to these terms, a numerical scheme obtained by performing an arbitrary scheme to each layer, and using time-splitting or other similar techniques leads to instabilities in...

We study a binary mixture of compressible viscous fluids modelled by the Navier-Stokes-Allen-Cahn system with isentropic or ideal gas law. We propose a finite volume method for the approximation of the system based on upwinding and artificial diffusion approaches. We prove the entropy stability of the numerical method and present several numerical experiments to support the theory.

We study error estimates and their convergence rates for approximate solutions of spectral Galerkin type for the equations for the motion of a viscous chemical active fluid in a bounded domain. We find error estimates that are uniform in time and also optimal in the L2-norm and H1-norm. New estimates in the H(-1)-norm are given.

We describe a numerical method for the equation $u_t+p{u}_x-\epsilon u_{xx}=f$ in $(0,1)\times (0,T)$ with Dirichlet boundary and initial conditions which is a combination of the method of characteristics and the finite-difference method. We prove both an a priori local error-estimate of a high order and stability. Example 3.3 indicates that our approximate solutions are disturbed only by a minimal amount of the artificial diffusion.

In this paper Rothe’s classical method is extended so that it can be used to solve some linear parabolic boundary value problems in non-cylindrical domains. The corresponding existence and uniqueness theorems are proved and some further results and generalizations are discussed and applied.