This paper analyzes an abstract two-level algorithm for hybridizable discontinuous Galerkin (HDG) methods in a unified fashion. We use an extended version of the Xu-Zikatanov (X-Z) identity to derive a sharp estimate of the convergence rate of the algorithm, and show that the theoretical results also are applied to weak Galerkin (WG) methods. The main features of our analysis are twofold: one is that we only need the minimal regularity of the model problem; the other is that we do not require the triangulations to be quasi-uniform. Numerical experiments are provided to confirm the theoretical results.

This article studies propagating traveling waves in a class of reaction-diffusion systems which model isothermal autocatalytic chemical reactions as well as microbial growth and competition in a flow reactor. In the context of isothermal autocatalytic systems, two different cases will be studied. The first is autocatalytic chemical reaction of order m without decay. The second is chemical reaction of order m with a decay of order n, where m and n are positive integers and m>n≥1. A typical system in autocatalysis is A+2B→3B and B→C involving two chemical species, a reactant A and an auto-catalyst B and C an inert chemical species.

The numerical computation gives more accurate estimates on minimum speed of traveling waves for autocatalytic reaction without decay, providing useful insight in the study of stability of traveling waves.

For autocatalytic reaction of order m = 2 with linear decay n = 1, which has a particular important role in chemical waves, it is shown numerically that there exist multiple traveling waves with 1, 2 and 3 peaks with certain choices of parameters.

In the first part, we study the convergence of discrete solutions to splitting schemes for two-phase flow with different mass densities suggested in [Guillen-Gonzalez, Tierra, J.Comput.Math. (6)2014]. They have been formulated for the diffuse interface model in [Abels, Garcke, Grün, M3AS, 2012, DOI:10.1142/S0218202511500138] which is consistent with thermodynamics. Our technique covers various discretization methods for phase-field energies, ranging from convex-concave splitting to difference quotient approaches for the double-well potential. In the second part of the paper, numerical experiments are presented in two space dimensions to identify discretizations of Cahn-Hilliard energies which are ϕ-stable and which do not reduce the acceleration of falling droplets. Finally, 3d simulations in axial symmetric geometries are shown to underline even more the full practicality of the approach.

The miscible displacement of one incompressible fluid by another in a porous medium is governed by a system of two equations. One is elliptic form equation for the pressure and the other is parabolic form equation for the concentration of one of the fluids. Since only the velocity and not the pressure appears explicitly in the concentration equation, we use a mixed finite element method for the approximation of the pressure equation and mixed finite element method with characteristics for the concentration equation. To linearize the mixed-method equations, we use a two-grid algorithm based on the Newton iteration method for this full discrete scheme problems. First, we solve the original nonlinear equations on the coarse grid, then, we solve the linearized problem on the fine grid used Newton iteration once. It is shown that the coarse grid can be much coarser than the fine grid and achieve asymptotically optimal approximation as long as the mesh sizes satisfy h = H2 in this paper. Finally, numerical experiment indicates that two-grid algorithm is very effective.

The present work demonstrates the extension of a thermal-hydraulic model, THRUST, with an objective to simulate the fast transient flow dynamics in a supercritical water channel of circular cross section. THRUST is a 1-D model which solves the nonlinearly coupled mass, axial momentum and energy conservation equations in time domain based on a characteristics-dependent fully implicit finite difference scheme using an Eulerian approach. The model developed accounts for the compressibility of the supercritical flow by considering the finite value of acoustic speed in the solution algorithm and treats the boundary conditions naturally. A supercritical water channel of circular cross section, for which the experimental data is available at steady state operating conditions, is chosen for the transient simulations to start with. Two different case studies are undertaken with a purpose to assess the capability of the model to analyze the fast transient processes caused by the large reduction in system pressure. The first transient case study is where the initial exit pressure is reduced by 1 MPa exponentially in a time span of 5 s. In the second case study, the transient is initiated with a sudden step decrease in the exit pressure by the same amount. Results obtained for both the case studies show the desired performance from the model developed.

It is known that large time-stepping method are useful for simulating phase field models. In this work, an adaptive time-stepping strategy is proposed based on numerical energy stability and equi-distribution principle. The main idea is to use the energy variation as an indicator to update the time step, so that the resulting algorithm is free of user-defined parameters, which is different from several existing approaches. Some numerical experiments are presented to illustrate the effectiveness of the algorithms.