Search results

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 determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of microfluidic droplet ensembles, water flow in unsaturated flows, infiltration, etc., as functions of the powers of the convection and diffusion fluxes and upstream boundary condition; to study numerically the evolution of the wave for two different initial conditions; and to assess the accuracy of several finite difference methods for the solution of the degenerate, nonlinear, advection--diffusion equation that governs the model.

The theory of ordinary differential equations and several explicit, finite difference methods that use first- and second-order, accurate upwind, central and compact discretizations for the convection terms are used to determine the analytical solution for steadily propagating waves and the evolution of the wave fronts from hyperbolic tangent and piecewise linear initial conditions to steadily propagating waves, respectively. The amplitude and phase errors of the semi-discrete schemes are determined analytically and the accuracy of the discrete methods is assessed.

For non-zero upstream boundary conditions, it has been found both analytically and numerically that the shock wave is smooth and its steepness increases as the power of the diffusion term is increased and as the upstream boundary value is decreased. For zero upstream boundary conditions, smooth, cusped and sharp shock waves may be encountered depending on the powers of the convection and diffusion terms. For a linear diffusion flux, the shock wave is smooth, whereas, for a quadratic diffusion flux, the wave exhibits a cusped front whose left spatial derivative decreases as the power of the convection term is increased. For higher nonlinear diffusion fluxes, a sharp shock wave is observed. The wave speed decreases as the powers of both the convection and the diffusion terms are increased. The evolution of the solution from hyperbolic tangent and piecewise linear initial conditions shows that the wave back adapts rapidly to its final steady value, whereas the wave front takes much longer, especially for piecewise linear initial conditions, but the steady wave profile and speed are independent of the initial conditions. It is also shown that discretization of the nonlinear diffusion flux plays a more important role in the accuracy of first- and second-order upwind discretizations of the convection term than either a conservative or a non-conservative discretization of the latter. Second-order upwind and compact discretizations of the convection terms are shown to exhibit oscillations at the foot of the wave’s front where the solution is nil but its left spatial derivative is largest. The results obtained with a conservative, centered second--order accurate finite difference method are found to be in good agreement with those of the second-order accurate, central-upwind Kurganov--Tadmor method which is a non-oscillatory high-resolution shock-capturing procedure, but differ greatly from those obtained with a non-conservative, centered, second-order accurate scheme, where the gradients are largest.

A new, one-dimensional model for microfluidic droplet transport, water flow in unsaturated flows, infiltration, etc., that includes high-order convection fluxes and degenerate diffusion, is proposed and studied both analytically and numerically. Its smooth, cusped and sharp shock wave solutions have been determined analytically as functions of the powers of the nonlinear convection and diffusion fluxes and the boundary conditions. These solutions are used to assess the accuracy of several finite difference methods that use different orders of accuracy in space, and different discretizations of the convection and diffusion fluxes, and can be used to assess the accuracy of other numerical procedures for one-dimensional, degenerate, convection--diffusion equations.

The purpose of this paper is to investigate the two-phase flow problems involving gas–liquid mixture.

The governed equations consist of a range of conservation laws modeling a classification of two-phase flow phenomena subjected to a velocity nonequilibrium for the gas–liquid mixture. Effects of the relative velocity are accounted for in the present model by a kinetic constitutive relation coupled to a collection of specific equations governing mass and volume fractions for the gas phase. Unlike many two-phase models, the considered system is fully hyperbolic and fully conservative. The suggested relaxation approach switches a nonlinear hyperbolic system into a semilinear model that includes a source relaxation term and characteristic linear properties. Notably, this model can be solved numerically without the use of Riemann solvers or linear iterations. For accurate time integration, a high-resolution spatial reconstruction and a Runge–Kutta scheme with decreasing total variation are used to discretize the relaxation system.

The method is used in addressing various nonequilibrium two-phase flow problems, accompanied by a comparative study of different reconstructions. The numerical results demonstrate the suggested relaxation method’s high-resolution capabilities, affirming its proficiency in delivering accurate simulations for flow regimes characterized by strong shocks.

While relaxation methods exhibit notable performance and competitive features, as far as we are aware, there has been no endeavor to address nonequilibrium two-phase flow problems using these methods.

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.

A new approach for solving steady incompressible Navier‐Stokes equations is presented in this paper. This method extends the upwind Riemann‐problem‐based techniques to viscous flows, which is obtained by applying modified artificial compressibility Navier‐Stokes equations and fully discrete high‐order numerical schemes for systems of advection‐diffusion equations. In this approach, utilizing the local Riemann solutions the steady incompressible viscous flows can be solved in a similar way to that of inviscid hyperbolic conservation laws. Numerical experiments on the driven cavity problem indicate that this approach can give satisfactory solutions.

This paper aims to construct a sixth-order weighted essentially nonoscillatory scheme for simulating the nonlinear degenerate parabolic equations in a finite difference framework.

To design this scheme, we approximate the second derivative in these equations in a different way, which of course is still in a conservative form. In this way, unlike the common practice of reconstruction, the approximation of the derivatives of odd order is needed to develop the numerical flux.

The results obtained by the new scheme produce less error compared to the results of other schemes in the literature that are recently developed for the nonlinear degenerate parabolic equations while requiring less computational times.

This research develops a new weighted essentially nonoscillatory scheme for solving the nonlinear degenerate parabolic equations in multidimensional space. Besides, any selection of the constants (sum equals one is the only requirement for them), named the linear weights, will obtain the desired accuracy.

The purpose of this paper is to present a fully conservative numerical algorithm for solving the coupled shallow water hydro-sediment-morphodynamic equations governing fluvial processes, and also to clarify the performance of a conventional algorithm, which redistributes the variable water-sediment mixture density to the source terms of the governing equations and accordingly the hyperbolic operator is rendered similar to that of the conventional shallow water equations for clear water flows.

The coupled shallow water hydro-sediment-morphodynamic equations governing fluvial processes are arranged in full conservation form, and solved by a well-balanced weighted surface depth-gradient method along with a slope-limited centred scheme. The present algorithm is verified for a spectrum of test cases, which involve complex flows with shock waves and sediment transport processes with contact discontinuities over irregular topographies. The computational results of the conventional algorithm are compared with those of the present algorithm and evaluated by available referenced data.

The fully conservative numerical algorithm performs satisfactorily over the spectrum of test cases, and the conventional algorithm is confirmed to work similarly well.

A fully conservative numerical algorithm, without redistributing the water-sediment mixture density, is proposed for solving the coupled shallow water hydro-sediment-morphodynamic equations. It is clarified that the conventional algorithm, involving redistribution of the water-sediment mixture density, performs similarly well. Both algorithms are equally applicable to problems encountered in computational river modelling.

For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the restrictive time step introduced by this higher order term. The purpose of this paper is to employ a fractional splitting method to isolate the convective, the nonlinear second-order and the fourth-order differential terms.

The full equation is then solved by consistent schemes for each differential term independently. In addition to validating the second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.

The scheme is second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.

The authors believe that this is the first time the equation is handled numerically using the fractional step method. Apart from the fact that the fractional step method substantially reduces computational time, it has the advantage of simplifying a complex process efficiently. This method permits the treatment of each segment of the original equation separately and piece them together, in a way that will be explained shortly, without destroying the properties of the equation.

This paper aims to numerically study the compositional flow of two- and three-phase fluids in one-dimensional porous media and to make a comparison between several upwind and central numerical schemes.

Implicit pressure explicit composition (IMPEC) procedure is used for discretization of governing equations. The pressure equation is solved implicitly, whereas the mass conservation equations are solved explicitly using different upwind (UPW) and central (CEN) numerical schemes. These include classical upwind (UPW-CLS), flux-based decomposition upwind (UPW-FLX), variable-based decomposition upwind (UPW-VAR), Roe’s upwind (UPW-ROE), local Lax–Friedrichs (CEN-LLF), dominant wave (CEN-DW), Harten–Lax–van Leer (HLL) and newly proposed modified dominant wave (CEN-MDW) schemes. To achieve higher resolution, high-order data generated by either monotone upstream-centered schemes for conservation laws (MUSCL) or weighted essentially non-oscillatory (WENO) reconstructions are used.

It was found that the new CEN-MDW scheme can accurately solve multiphase compositional flow equations. This scheme uses most of the information in flux function while it has a moderate computational cost as a consequence of using simple algebraic formula for the wave speed approximation. Moreover, numerically calculated wave structure is shown to be used as a tool for a priori estimation of problematic regions, i.e. degenerate, umbilic and elliptic points, which require applying correction procedures to produce physically acceptable (entropy) solutions.

This paper is concerned with one-dimensional study of compositional two- and three-phase flows in porous media. Temperature is assumed constant and the physical model accounts for miscibility and compressibility of fluids, whereas gravity and capillary effects are neglected.

The proposed numerical scheme can be efficiently used for solving two- and three-phase compositional flows in porous media with a low computational cost which is especially useful when the number of chemical species increases.

A new central scheme is proposed that leads to improved accuracy and computational efficiency. Moreover, to the best of authors knowledge, this is the first time that the wave structure of compositional model is investigated numerically to determine the problematic situations during numerical solution and adopt appropriate correction techniques.

The discontinuous Galerkin finite element method (DGFEM) is very suited for realizing high order resolution approximations on unstructured grids for calculating the hyperbolic conservation law. However, it requires a significant amount of computing resources. Therefore, this paper aims to investigate how to solve the Euler equations in parallel systems and improve the parallel performance.

Discontinuous Galerkin discretization is used for the compressible inviscid Euler equations. The multi-level domain decomposition strategy was used to deal with the computational grids and ensure the calculation load balancing. The total variation diminishing (TVD) Runge–Kutta (RK) scheme coupled with the multigrid strategy was employed to further improve parallel efficiency. Moreover, the Newton Block Gauss–Seidel (GS) method was adopted to accelerate convergence and improve the iteration efficiency.

Numerical experiments were implemented for the compressible inviscid flow problems around NACA0012 airfoil, over M6 wing and DLR-F6 configuration. The parallel acceleration is near to a linear convergence. The results indicate that the present parallel algorithm can reduce computational time significantly and allocate memory reasonably, which has high parallel efficiency and speedup, and it is well-suited to large-scale scientific computational problems on multiple instruction stream multiple data stream model.

The parallel DGFEM coupled with TVD RK and the Newton Block GS methods was presented for hyperbolic conservation law on unstructured meshes.