Search results

This paper aims to clarify some aspects of the application of the Godunov method for the Baer–Nunziato equations solution on the example of the problem of shock wave – dense particles cloud interaction.

The statement of the problem corresponds to the natural experiment. Mathematical model is based on the Baer–Nunziato system of equations with algebraic right-hand side source terms that takes into account the interphase friction force. Two numerical approaches are used: Harten-Lax-van Leer method and Godunov method.

For the robust simulation using Godunov method, the application of the pressure relaxation procedure is proposed. The comparative analysis of the simulation results using two methods is carried out. The Godunov method provides significantly smaller numerical diffusion of the solid phase volume fraction in the cloud that leads to the much better agreement of the pressure curves on transducers and the dynamics of the cloud motion with the experimental data.

Godunov method for the Baer–Nunziato equations is applied for the simulation of the natural experiment on the shock wave particles cloud interaction. Up to now, the examples of the application of the Godunov method for the Baer–Nunziato equations to the investigation of the practical problems have been limited by the works of the authors of the method and the field of detonation in the heterogeneous explosives. For the robust simulations in the presence of interphase boundaries, it is proposed to use the Godunov method together with the pressure relaxation procedure.

This paper aims to perform a comparative study between capabilities of two numerical schemes from two main branches of numerical methods for solving hyperbolic conservation equations.

The accuracy and performance of a newly developed high‐resolution central scheme vs a higher‐order Godunov‐based method are evaluated in the context of black‐oil reservoir simulations. Both methods are modified enabling study of applications that are not strictly hyperbolic and exhibit local linear degeneracies in their wave structure.

The numerical computations show that while both schemes produce results with virtually the same accuracy, the Godunov method reproduces slightly more accurate results at the expense of calculation of eigen‐structures.

The paper investigates only one dimensional problems, but the idea can be easily extended to multi‐dimensional problems.

The paper shows the strengths and weaknesses of two practical numerical methods.

Such comparative study has not been published elsewhere and in particular, the performance of numerical methods on simulating hysteresis effect in hydrocarbon reservoirs has not been investigated in detail before.

The purpose of this paper is twofold. Firstly, to present a detailed account of the generalized Lagrangian formulation of Hui and Zhao, in which the stream function ζ and Lagrangian distance λ, are used as independent variables, and secondly to assess and compare the performance of various flux limiters in this formulation with their corresponding performance in the Eulerian formulation. The generalized Lagrangian formulation is obtained by a transformation from the cartesian co‐ordinates (x, y) to the Lagrangian co‐ordinates (λ, ζ). In this manner, the number of independent variables for steady, 3‐D flow is reduced from four to three, placing this formulation on the same footing as the Eulerian formulation even for steady flows (as opposed to the conventional Lagrangian formulation which apparently still requires four independent variables even for steady flows). The generalized Lagrangian formulation with the Godunov scheme (using flux limiters) appears to have distinct advantages over the corresponding Eulerian formulation, particularly with respect to accuracy. Furthermore, the method requires no grid generation.

The purpose of this paper is to investigate the performance of a specific class of high‐resolution central schemes in conjunction with the black oil models for hydrocarbon reservoir simulation.

A generalized black oil model is adopted, in which the solubility of gas in both oil and water and evaporation of oil are considered, leading to a system of equations prone to degeneracy. A computer code is generated and three test cases are solved to evaluate the performance of various schemes in terms of accuracy and discontinuity handling.

It is shown that, although some of the central schemes are highly sensitive to the choice of Courant‐Friedrich‐Levy (CFL) number and produce overly diffusive results, a certain type of this class is insensitive to the CFL number and can conveniently handle degenerate equations appearing in the reservoir simulation. The obtained results are compared with those available in the literature, showing merits of this class of schemes in complex reservoir simulation models.

This paper gives the one‐dimensional implementation of the above‐mentioned schemes. Extension to higher dimensional black oil model is currently under development by the authors.

The specific class of high‐resolution central schemes investigated here presents the same level of accuracy as more complicated numerical methods, yet keeping it much more simple, by avoiding Riemann solvers.

The high‐resolution central scheme used in this work has been newly developed and applied to simple scalar hyperbolic equations. It has been adopted for the black oil for the first time.

This paper aims to improve the reliability of numerical methods for predicting the transient heat transfers in combustion chambers heated internally by moving heat sources.

A two-phase fluid dynamic model was used to govern the non-uniformly distributed moving heat sources. A Riemann-problem-based numerical scheme was provided to update the fluid field and provide convective boundary conditions for the heat transfer. The heat conduction in the solids was investigated by using a thermo-mechanical coupled model to obtain a reliable expanding velocity of the heat sources. The coupling between the combustion and the heat transfer is realized based on user subroutines VDFLUX and VUAMP in the commercial software ABAQUS.

The capability of the numerical scheme in capturing discontinuities in initial conditions and source terms was validated by comparing the predicted results of commonly used verification cases with the corresponding analytical solutions. The coupled model and the numerical methods are capable of investigating heat transfer problems accompanied by extreme conditions such as transient effects, high-temperature and high-pressure working conditions.

The work provides a reliable numerical method to obtain boundary conditions for predicting the heat transfers in solids heated by expanding multiphase reactive flows.

An optimized multigrid method (NSFLEX‐MG) for the NSFLEX‐code (Navier‐Stokes solver using characteristic flux extrapolation) of MBB (Messerschmitt Bolkow Blohm GmbH) is described. The method is based on a correction scheme and implicit relaxation procedures and is applied to two‐dimensional test cases. The principal feature of the flow solver is a Godunov‐type averaging procedure based on the eigenvalues analysis of the Euler equations by means of which the inviscid fluxes are evaluated at the finite volume faces. Viscous fluxes are centrally differenced at each cell face. The performance of NSFLEX‐MG is demonstrated for a large range of Mach numbers for compressible inviscid and viscous (laminar and turbulent) flows over a RAE‐2822 airfoil and over a NACA‐0012 airfoil.

The purpose of this work is to present the implementation of weighted essentially non-oscillatory (WENO) wavelet methods for solving multiphase flow problems. The particular interest is gas–liquid two-phase mixture with velocity non-equilibrium. Numerical simulations are carried out on different scenarios of one-dimensional Riemann problems for gas–liquid flows. Results are validated and qualitatively compared with solutions provided by other standard numerical methods.

This paper extends the framework of WENO wavelet adaptive method to a fully hyperbolic two-phase flow model in a conservative form. The grid adaptivity in each time step is provided by the application of a thresholded interpolating wavelet transform. This facilitates the construction of a small yet effective sparse point representation of the solution. The method of Lax–Friedrich flux splitting is used to resolve the spatial operator in which the flux derivatives are approximated by the WENO scheme.

Hyperbolic models of two-phase flow in conservative form are efficiently solved, as shocks and rarefaction waves are precisely captured by the chosen methodology. Substantial computational gains are obtained through the grid reduction feature while maintaining the quality of the solutions. The results indicate that WENO wavelet methods are robust and sufficient to accurately simulate gas–liquid mixtures.

Resolution of two-phase flows is rarely studied using WENO wavelet methods. It is the first time such a study on the relative velocity is reported in two-phase flows using such methods.

The purpose of this paper is to present a new hybrid Euler flux fonction for use in a finite-volume Euler/Navier-Stokes code and adapted to compressible flow problems.

The proposed scheme, called AUFSRR can be devised by combining the AUFS solver and the Roe solver, based on a rotated Riemann solver approach (Sun and Takayama, 2003; Ren, 2003). The upwind direction is determined by the velocity-difference vector and idea is to apply the AUFS solver in the direction normal to shocks to suppress carbuncle and the Roe solver across shear layers to avoid an excessive amount of dissipation. The resulting flux functions can be implemented in a very simple manner, in the form of the Roe solver with modified wave speeds, so that converting an existing AUFS flux code into the new fluxes is an extremely simple task.

The proposed flux functions require about 18 per cent more CPU time than the Roe flux. Accuracy, efficiency and other essential features of AUFSRR scheme are evaluated by analyzing shock propagation behaviours for both the steady and unsteady compressible flows. This is demonstrated by several test cases (1D and 2D) with standard finite-volume Euler code, by comparing results with existing methods.

The hybrid Euler flux function is used in a finite-volume Euler/Navier-Stokes code and adapted to compressible flow problems.

The AUFSRR scheme is devised by combining the AUFS solver and the Roe solver, based on a rotated Riemann solver approach.

The aim of the study is to present a piecewise parabolic method (PPM) for numerical simulation of barotropic and nonbarotropic two‐fluid flows in more than one space dimension.

In transition layers of two components, a fluid mixture model system is introduced. Besides, conserving the mass, momentum and energy for the mixture, the model is supplemented with an advection equation for the volume fraction of one of the two fluid components to recover the pressure and track interfaces. The Tait and stiffened gas equations of state are used to describe thermodynamic properties of the barotropic and nonbarotropic components, respectively. To close the model system, a mixture equation of state is derived. The classical third‐order PPM is extended to the two‐fluid case and used to solve the model system.

The feasibility of this method has been demonstrated by good results of sample applications. Each of the material interfaces is resolved with two grid cells and there is no any pressure oscillation on the interfaces.

With the mixture model system, there may be energy gain or loss for the nonbarotropic component on the material interfaces.

The method can be applied to a wide range of practical problems.

The method is simple. It not only has the advantage of Lagrangian‐type schemes but also keeps the robustness of Eulerian schemes.

Partially‐decoupled upwind‐based total‐variation‐diminishing (TVD) finite‐difference schemes for the solution of the conservation laws governing two‐dimensional non‐equilibrium vibrationally relaxing and chemically reacting flows of thermally‐perfect gaseous mixtures are presented. In these methods, a novel partially‐decoupled flux‐difference splitting approach is adopted. The fluid conservation laws and species concentration and vibrational energy equations are decoupled by means of a frozen flow approximation. The resulting partially‐decoupled gas‐dynamic and thermodynamic subsystems are then solved alternately in a lagged manner within a time marching procedure, thereby providing explicit coupling between the two equation sets. Both time‐split semi‐implicit and factored implicit flux‐limited TVD upwind schemes are described. The semi‐implicit formulation is more appropriate for unsteady applications whereas the factored implicit form is useful for obtaining steady‐state solutions. Extensions of Roe's approximate Riemann solvers, giving the eigenvalues and eigenvectors of the fully coupled systems, are used to evaluate the numerical flux functions. Additional modifications to the Riemann solutions are also described which ensure that the approximate solutions are not aphysical. The proposed partially‐decoupled methods are shown to have several computational advantages over chemistry‐split and fully coupled techniques. Furthermore, numerical results for single, complex, and double Mach reflection flows, as well as corner‐expansion and blunt‐body flows, using a five‐species four‐temperature model for air demonstrate the capabilities of the methods.