A roe-average algorithm for a granular-gas model with non-conservative terms

A Roe-average algorithm has been derived for a granular-gas model, proposed by Goldshtein and Shapiro [Goldshtein, Shapiro, Mechanics of collisional motion of granular materials: Part 1. General hydrodynamic equations, J. Fluid Mech. 282 (1995) 75–114], which contains non-conservative terms in the Euler-like hyperbolic governing equations apart from sink terms, which arise from inelastic collision of granules and are present only in the energy equation. The non-conservative terms introduce non-isentropic effects in acoustic-wave propagation within granular media and they also contribute to the Rankine–Hugoniot relations across a discontinuity. A Roe-average algorithm, based on the same granular-gas model, was derived in the literature [V. Kamenetsky, A. Goldshtein, M. Shapiro, D. Degani, Evolution of a shock wave in a granular gas, Phys. Fluids, 12 (2000) 3036–3049] which then required the implementation of a shock-fitting technique at a discontinuity. In the present work, Roe-averaged variables have been obtained from the Rankine–Hugoniot jump relations and the non-conservative terms have been incorporated in the numerical flux formula consistent with upwind principles associated with the granular speed of sound. Results for unsteady one-dimensional granular flows, colliding with a wall, demonstrate the capability of the proposed algorithm to capture strong shocks in addition to flow features not found in molecular gases, such as a fluidized region downstream of the shock and a compacted solid-block region adjacent to the wall.     

A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry

The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes.     

A conservative discrete compatibility-constraint low-Mach pressure-correction algorithm for time-accurate simulations of variable density flows

In this paper, we develop the discrete compatibility-constraint pressure-correction algorithm for transient simulations of variable density flows at low-Mach numbers. The constraint for the velocity field is constructed from a combination of the discrete equations of continuity and scalar (e.g. energy) transport, imposing that the newly predicted state must be compatible, in agreement with the equation of state. This way, mass and scalar conservation are guaranteed and the equation of state is exactly fulfilled at every time step. For comparison reasons, two other types of well-known pressure-correction algorithms are also used. The first class, denoted as continuity-constraint pressure-correction, is based on a constraint for the velocity field that is derived solely from the continuity equation. The second class, denoted as analytical compatibility-constraint pressure-correction, constructs the constraint from an analytical combination of the material derivative of the equation of state and the continuity and scalar equations. The algorithms are tested for three example fluid configurations: a single-fluid ideal gas, a two-fluid inert mixture and a two-fluid reacting mixture. The latter is special in the sense that the equation of state is non-linear and not everywhere differentiable. The continuity-constraint pressure-correction algorithm yields unstable solutions if density ratios are high. The analytical compatibility-constraint pressure-correction algorithm yields stable results, but the predicted states do not correspond to the equation of state. The discrete compatibility-constraint pressure-correction algorithm performs well on all test cases: the simulation results are stable and exactly match the equation of state.     

An immersed boundary–lattice-Boltzmann method for the simulation of the flow past an impulsively started cylinder

We present a lattice-Boltzmann method coupled with an immersed boundary technique for the simulation of bluff body flows. The lattice-Boltzmann method for the modeling of the Navier–Stokes equations, is enhanced by a forcing term to account for the no-slip boundary condition on a non-grid conforming boundary. We investigate two alternatives of coupling the boundary forcing term with the grid nodes, namely the direct and the interpolated forcing techniques. The present LB–IB methods are validated in simulations of the incompressible flow past an impulsively started cylinder at low and moderate Reynolds numbers. We present diagnostics such as the near wall vorticity field and the drag coefficient and comparisons with previous computational and experimental works and assess the advantages and drawbacks of the two techniques.     

A high order kinetic flux-vector splitting method for the reduced five-equation model of compressible two-fluid flows

We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) [21]. The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two space dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge–Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor [14]. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.     

Splitting multisymplectic integrators for Maxwell’s equations

In the paper, we describe a novel kind of multisymplectic method for three-dimensional (3-D) Maxwell’s equations. Splitting the 3-D Maxwell’s equations into three local one-dimensional (LOD) equations, then applying a pair of symplectic Runge–Kutta methods to discretize each resulting LOD equation, it leads to splitting multisymplectic integrators. We say this kind of schemes to be LOD multisymplectic scheme (LOD-MS). The discrete conservation laws, convergence, dispersive relation, dissipation and stability are investigated for the schemes. Theoretical analysis shows that the schemes are unconditionally stable, non-dissipative, and of first order accuracy in time and second order accuracy in space. As a reduction, we also consider the application of LOD-MS to 2-D Maxwell’s equations. Numerical experiments match the theoretical results well. They illustrate that LOD-MS is not only efficient and simple in coding, but also has almost all the nature of multisymplectic integrators.     

A numerical algorithm for viscous incompressible interfacial flows

We present a new algorithm to numerically simulate two-dimensional viscous incompressible flows with moving interfaces. The motion is updated in time by using the backward difference formula through an iterative procedure. At each iteration, the pseudo-spectral technique is applied in the horizontal direction. The resulting semi-discretized equations constitute a boundary value problem in the vertical coordinate which is solved by decoupling growing and decaying solutions. Numerical tests justify that this method achieves fully second-order accuracy in both the temporal variable and vertical coordinate. As an application of this algorithm, we study the motion of Stokes waves in the presence of viscosity. Our numerical results are consistent with the recently published asymptotic solution for Stokes waves in slightly viscous fluids.     

A second order numerical scheme for the improved Boussinesq equation

A finite-difference scheme arising from the use of rational approximants to the matrix-exponential term in a three-time level recurrence relation is used for the numerical solution of the improved Boussinesq equation (IBq). The resulting linear scheme, which is analyzed for local truncation error and stability, is tested numerically and conclusions with corresponding results known in the bibliography are derived.     

A hybrid scheme for computing incompressible two-phase flows

We propose a hybrid scheme for computations of incompressible two-phase flows. The incompressible constraint has been replaced by a pressure Poisson-like equation and then the pressure is updated by the modified marker and cell method. Meanwhile, the moment equations in the incompressible Navier-Stokes equations are solved by our semidiscrete Hermite central-upwind scheme, and the interface between the two fluids is considered to be continuous and is described implicitly as the 0.5 level set of a smooth function being a smeared out Heaviside function. It is here named the hybrid scheme. Some numerical experiments are successfully carried out, which verify the desired efficiency and accuracy of our hybrid scheme.     

A remark on compact matrix quantum groups

A simpler set of axioms of the theory of compact matrix quantum groups (pseudogroups) is found.Supported by Japan Society for Promoting Science. On leave from the Department of Mathematical Methods in Physics, Faculty of Physics, University of Warsaw, Hoa 74, 00-682 Warsaw, Poland.     

A finite element variational multiscale method for incompressible flows based on two local gauss integrations

In this article, we present a finite element variational multiscale (VMS) method for incompressible flows based on two local Gauss integrations, and compare it with common VMS method which is defined by a low order finite element space L_h$L_h$ on the same grid as X_h$X_h$ for the velocity deformation tensor and a stabilization parameter α$\alpha$. The best algorithmic feature of our method is using two local Gauss integrations to replace projection operator. We theoretically discuss the relationship between our method and common VMS method for the Taylor–Hood elements, and show that the nonlinear system derived from our method by finite element discretization is much smaller than that of common VMS method computationally.     

A bubble-stabilized least-squares finite element method for steady MHD duct flow problems at high Hartmann numbers

In this paper we devise a stabilized least-squares finite element method using the residual-free bubbles for solving the governing equations of steady magnetohydrodynamic duct flow. We convert the original system of second-order partial differential equations into a first-order system formulation by introducing two additional variables. Then the least-squares finite element method using C⁰$C^0$ linear elements enriched with the residual-free bubble functions for all unknowns is applied to obtain approximations to the first-order system. The most advantageous features of this approach are that the resulting linear system is symmetric and positive definite, and it is capable of resolving high gradients near the layer regions without refining the mesh. Thus, this approach is possible to obtain approximations consistent with the physical configuration of the problem even for high values of the Hartmann number. Before incoorperating the bubble functions into the global problem, we apply the Galerkin least-squares method to approximate the bubble functions that are exact solutions of the corresponding local problems on elements. Therefore, we indeed introduce a two-level finite element method consisting of a mesh for discretization and a submesh for approximating the computations of the residual-free bubble functions. Numerical results confirming theoretical findings are presented for several examples including the Shercliff problem.     

A high-resolution mapped grid algorithm for compressible multiphase flow problems

We describe a simple mapped-grid approach for the efficient numerical simulation of compressible multiphase flow in general multi-dimensional geometries. The algorithm uses a curvilinear coordinate formulation of the equations that is derived for the Euler equations with the stiffened gas equation of state to ensure the correct fluid mixing when approximating the equations numerically with material interfaces. A γ-based and a α-based model have been described that is an easy extension of the Cartesian coordinates counterpart devised previously by the author [30]. A standard high-resolution mapped grid method in wave-propagation form is employed to solve the proposed multiphase models, giving the natural generalization of the previous one from single-phase to multiphase flow problems. We validate our algorithm by performing numerical tests in two and three dimensions that show second order accurate results for smooth flow problems and also free of spurious oscillations in the pressure for problems with interfaces. This includes also some tests where our quadrilateral-grid results in two dimensions are in direct comparisons with those obtained using a wave-propagation based Cartesian grid embedded boundary method.     

A finite volume method for approximating 3D diffusion operators on general meshes

A finite volume method is presented for discretizing 3D diffusion operators with variable full tensor coefficients. This method handles anisotropic, non-symmetric or discontinuous variable tensor coefficients while distorted, non-matching or non-convex n-faced polyhedron meshes can be used. For meshes of polyhedra whose faces have not more than four edges, the associated matrix is positive definite (and symmetric if the diffusion tensor is symmetric). A second-order (resp. first-order) accuracy is numerically observed for the solution (resp. gradient of the solution).     

A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part II: On an arbitrary collocated mesh

A conservative formulation of the Lorentz force is given here for magnetohydrodynamic (MHD) flows at a low magnetic Reynolds number with the current density calculated based on Ohm’s law and the electrical potential formula. This conservative formula shows that the total momentum contributed from the Lorentz force is conservative when the applied magnetic field is constant. For the case with a non-constant applied magnetic field, the Lorentz force has been divided into two parts: a strong globally conservative part and a weak locally conservative part.The conservative formula has been employed to develop a conservative scheme for the calculation of the Lorentz force on an unstructured collocated mesh. Only the current density fluxes on the cell faces, which are calculated using a consistent scheme with good conservation, are needed for the calculation of the Lorentz force. Meanwhile, a conservative interpolation technique is designed to get the current density at the cell center from the current density fluxes on the cell faces. This conservative interpolation can keep the current density at the cell center conservative, which can be used to calculate the Lorentz force at the cell center with good accuracy. The Lorentz force calculated from the conservative current at the cell center is equivalent to the Lorentz force from the conservative formula when the applied magnetic field is constant, which can conserve the total momentum. We will further prove that the simple interpolation scheme used in the Part I [M.-J. Ni, R. Munipalli, N.B. Morley, P.Y. Huang, M. Abdou, A current density conservative scheme for MHD flows at a low magnetic Reynolds number. Part I. On a rectangular collocated grid system, Journal of Computational Physics, in press, doi:10.1016/j.jcp.2007.07.025] of this series of papers is conservative on a rectangular grid and can keep the total momentum conservative in a rectangular grid.