Search results

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 study aims to propose a robust total variation diminishing (TVD) weighted average flux (WAF) finite volume scheme for investigating compressible gas–liquid mixture flows.

This study considers a two-phase flow composed of a liquid containing dispersed gas bubbles. To model this two-phase mixture, this paper uses a homogeneous equilibrium model (HEM) defined by two mass conservation laws for the two phases and a momentum conservation equation for the mixture. It is assumed that the velocity is the same for the two phases, and the density of phases is governed by barotropic laws. By applying the theory of hyperbolic equations, this study establishes an exact solution of the Riemann problem associated with the model equations, which allows to construct an exact Riemann solver within the first-order upwind Godunov scheme as well as a robust TVD WAF scheme.

The ability and robustness of the proposed TVD WAF scheme is validated by testing several two-phase flow problems involving different wave structures of the Riemann problem. Simulation results are compared against analytical solutions and other available numerical methods as well as experimental data in the literature. The proposed approach is much superior to other strategies in terms of the accuracy and ability of reconstruction.

The novelty of this work lies in its methodical extension of a TVD WAF scheme implementing an exact Riemann solver developed for compressible two-phase flows. Furthermore, other novelty lies on the quantitative calculation of different Riemann problem two-phase flows. Simulation results involve the verification of the constructed methods on the exact solutions of HEM without any restriction of variables.

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.

A higher-order implicit shock-capturing scheme is presented for the Euler equations based on time linearization of the implicit flux vector rather than the residual vector.

The flux vector is linearized through a truncated Taylor-series expansion whose leading-order implicit term is an inner product of the flux Jacobian and the vector of differences between the current and previous time step values of conserved variables. The implicit conserved-variable difference vector is evaluated at cell faces by using the reconstructed states at the left and right sides of a cell face and projecting the difference between the left and right states onto the right eigenvectors. Flux linearization also facilitates the construction of implicit schemes with higher-order spatial accuracy (up to third order in the present study). To enhance the diagonal dominance of the coefficient matrix and thereby increase the implicitness of the scheme, wave strengths at cell faces are expressed as the inner product of the inverse of the right eigenvector matrix and the difference in the right and left reconstructed states at a cell face.

The accuracy of the implicit algorithm at Courant–Friedrichs–Lewy (CFL) numbers greater than unity is demonstrated for a number of test cases comprising one-dimensional (1-D) Sod’s shock tube, quasi 1-D steady flow through a converging-diverging nozzle, and two-dimensional (2-D) supersonic flow over a compression corner and an expansion corner.

The algorithm has the advantage that it does not entail spatial derivatives of flux Jacobian so that the implicit flux can be readily evaluated using Roe’s approximate Jacobian. As a result, this approach readily facilitates the construction of implicit schemes with high-order spatial accuracy such as Roe-MUSCL.

A novel finite-volume-based higher-order implicit shock-capturing scheme was developed that uses time linearization of fluxes at cell interfaces.

The purpose of this paper is to present a novel numerical method to solve incompressible flows with natural and mixed convections using pseudo‐compressibility formulation.

The present method is derived using the framework of Harten Lax and van Leer with contact (HLLC) method of Toro, Spruce and Spears, that was originally developed for compressible gas dynamics equations. This work generalizes the algorithm described in the previous paper to the case where heat transfer is involved. Here, the solution of the Riemann problem is approximated by a three‐wave system.

A few test cases involving incompressible laminar flows inside 2D square cavity for various Rayleigh and Reynolds numbers are considered for validating the present method. The computed results from the present method are found to be quite promising.

Although pseudo‐compressibility formulation has been found to have superior performance and has the potential to have numerical treatments similar to compressible flow equations, only two numerical methods have been applied so far; namely Jameson method and Roes flux difference splitting method. A new sophisticated numerical method, following the framework of HLLC method, is derived and implemented for solving pseudo‐compressibility‐based incompressible flow equations with heat transfer.

Reynolds-averaged Navier–Stokes (RANS) models often perform poorly in shock/turbulence interaction regions, resulting in excessive wall heat load and incorrect representation of the separation length in shockwave/turbulent boundary layer interactions. The authors suggest that this can be traced back to inadequate numerical treatment of the inviscid fluxes. The purpose of this study is an extension to the well-known Harten, Lax, van Leer, Einfeldt (HLLE) Riemann solver to overcome this issue.

It explicitly takes into account the broadening of waves due to the averaging procedure, which adds numerical dissipation and reduces excessive turbulence production across shocks. The scheme is derived based on the HLLE equations, and it is tested against three numerical experiments.

Sod’s shock tube case shows that the scheme succeeds in reducing turbulence amplification across shocks. A shock-free turbulent flat plate boundary layer indicates that smooth flow at moderate turbulence intensity is largely unaffected by the scheme. A shock/turbulent boundary layer interaction case with higher turbulence intensity shows that the added numerical dissipation can, however, impair the wall heat flux distribution.

The proposed scheme is motivated by implicit large eddy simulations that use numerical dissipation as subgrid-scale model. Introducing physical aspects of turbulence into the numerical treatment for RANS simulations is a novel approach.

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 study aims to evaluate blast loads on and the response of submerged structures.

An arbitrary Lagrangian–Eulerian method is developed to model fluid–structure interaction (FSI) problems of close-in underwater explosions (UNDEX). The “fluid” part provides the loads for the structure considers air, water and high explosive materials. The spatial discretization for the fluid domain is performed with a second-order vertex-based finite volume scheme with a tangent of hyperbola interface capturing technique. The temporal discretization is based on explicit Runge–Kutta methods. The structure is described by a large-deformation Lagrangian formulation and discretized via finite elements. First, one-dimensional test cases are given to show that the numerical method is free of mesh movement effects. Thereafter, three-dimensional FSI problems of close-in UNDEX are studied. Finally, the computation of UNDEX near a ship compartment is performed.

The difference in the flow mechanisms between rigid targets and deforming targets is quantified and evaluated.

Cavitation is modeled only approximately and may require further refinement/modeling.

The results demonstrate that the proposed numerical method is accurate, robust and versatile for practical use.

Better design of naval infrastructure [such as bridges, ports, etc.].

To the best of the authors’ knowledge, this study has been conducted for the first time.

The flux vector splitting (FVS) schemes are known for their higher resistance to shock instabilities and carbuncle phenomena in high-speed flow computations, which are generally accompanied by relatively large numerical diffusion. However, it is desirable to control the numerical diffusion of FVS schemes inside the boundary layer for improved accuracy in viscous flow computations. This study aims to develop a new methodology for controlling the numerical diffusion of FVS schemes for viscous flow computations with the help of a recently developed boundary layer sensor.

The governing equations are solved using a cell-centered finite volume approach and Euler time integration. The gradients in the viscous fluxes are evaluated by applying the Green’s theorem. For the inviscid fluxes, a new approach is introduced, where the original upwind formulation of an FVS scheme is first cast into an equivalent central discretization along with a numerical diffusion term. Subsequently, the numerical diffusion is scaled down by using a novel scaling function that operates based on a boundary layer sensor. The effectiveness of the approach is demonstrated by applying the same on van Leer’s FVS and AUSM schemes. The resulting schemes are named as Diffusion-Regulated van Leer’s FVS-Viscous (DRvLFV) and Diffusion-Regulated AUSM-Viscous (DRAUSMV) schemes.

The numerical tests show that the DRvLFV scheme shows significant improvement over its parent scheme in resolving the skin friction and wall heat flux profiles. The DRAUSMV scheme is also found marginally more accurate than its parent scheme. However, stability requirements limit the scaling down of only the numerical diffusion term corresponding to the acoustic part of the AUSM scheme.

To the best of the authors’ knowledge, this is the first successful attempt to regulate the numerical diffusion of FVS schemes inside boundary layers by applying a novel scaling function to their artificial viscosity forms. The new methodology can reduce the erroneous smearing of boundary layers by FVS schemes in high-speed flow applications.

The purpose of this paper is to present a new stabilised low-order finite element methodology for large strain fast dynamics.

The numerical technique describing the motion is formulated upon the mixed set of first-order hyperbolic conservation laws already presented by Lee et al. (2013) where the main variables are the linear momentum, the deformation gradient tensor and the total energy. The mixed formulation is discretised using the standard explicit two-step Taylor-Galerkin (2TG) approach, which has been successfully employed in computational fluid dynamics (CFD). Unfortunately, the results display non-physical spurious (or hourglassing) modes, leading to the breakdown of the numerical scheme. For this reason, the 2TG methodology is further improved by means of two ingredients, namely a curl-free projection of the deformation gradient tensor and the inclusion of an additional stiffness stabilisation term.

A series of numerical examples are carried out drawing key comparisons between the proposed formulation and some other recently published numerical techniques.

Both velocities (or displacements) and stresses display the same rate of convergence, which proves ideal in the case of industrial applications where low-order discretisations tend to be preferred. The enhancements introduced in this paper enable the use of linear triangular (or bilinear quadrilateral) elements in two dimensional nearly incompressible dynamics applications without locking difficulties. In addition, an artificial viscosity term has been added into the formulation to eliminate the appearance of spurious oscillations in the vicinity of sharp spatial gradients induced by shocks.