Search results

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.

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 main purpose of this paper is to introduce the gradient discretisation method (GDM) to a system of reaction diffusion equations subject to non-homogeneous Dirichlet boundary conditions. Then, the authors show that the GDM provides a comprehensive convergence analysis of several numerical methods for the considered model. The convergence is established without non-physical regularity assumptions on the solutions.

In this paper, the authors use the GDM to discretise a system of reaction diffusion equations with non-homogeneous Dirichlet boundary conditions.

The authors provide a generic convergence analysis of a system of reaction diffusion equations. The authors introduce a specific example of numerical scheme that fits in the gradient discretisation method. The authors conduct a numerical test to measure the efficiency of the proposed method.

This work provides a unified convergence analysis of several numerical methods for a system of reaction diffusion equations. The generic convergence is proved under the classical assumptions on the solutions.

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.

This article offered two meshfree algorithms, namely the local radial basis functions-finite difference (LRBF-FD) approximation and local radial basis functions-differential quadrature method (LRBF-DQM) to simulate the multidimensional hyperbolic wave models and work is an extension of Jiwari (2015).

In the evolvement of the first algorithm, the time derivative is discretized by the forward FD scheme and the Crank-Nicolson scheme is used for the rest of the terms. After that, the LRBF-FD approximation is used for spatial discretization and quasi-linearization process for linearization of the problem. Finally, the obtained linear system is solved by the LU decomposition method. In the development of the second algorithm, semi-discretization in space is done via LRBF-DQM and then an explicit RK4 is used for fully discretization in time.

For simulation purposes, some 1D and 2D wave models are pondered to instigate the chastity and competence of the developed algorithms.

The developed algorithms are novel for the multidimensional hyperbolic wave models. Also, the stability analysis of the second algorithm is a new work for these types of model.

This paper aims to describe a proposal for an innovative method of normal shock wave–turbulent boundary layer interaction (SBLI) and shock-induced separation control.

The concept is based on the introduction of a tangentially moving wall upstream of the shock wave and in the interaction region. The SBLI control mechanism may be implemented as a closed belt floating on an air cushion, sliding over two cylinders and forming the outer skin of the suction side of the airfoil. The presented exploratory numerical study is conducted with SPARC solver (steady 2D RANS). The effect of the moving wall is presented for the NACA 0012 airfoil operating in transonic conditions.

To assess the accuracy of obtained solutions, validation of the computational model is demonstrated against the experimental data of Harris, Ladson & Hill and Mineck & Hartwich (NASA Langley). The comparison is conducted not only for the reference (impermeable) but also for the perforated (permeable) surface NACA 0012 airfoils. Subsequent numerical analysis of SBLI control by moving wall confirms that for the selected velocity ratios, the method is able to improve the shock-upstream boundary layer and counteract flow separation, significantly increasing the airfoil aerodynamic performance.

The moving wall concept as a means of normal shock wave–turbulent boundary layer interaction and shock-induced separation control has been investigated in detail for the first time. The study quantified the necessary operational requirements of such a system and practicable aerodynamic efficiency gains and simultaneously revealed the considerable potential of this promising idea, stimulating a new direction for future investigations regarding SBLI control.

The purpose of this study is to investigate the interaction between magnetotactic bacteria and Fe₃O₄–water nanofluid (NF) in a wavy enclosure in the presence of 2D natural convection flow.

Uniform magnetic field (MF), Brownian and thermophoresis effects are also contemplated. The dimensionless, time-dependent equations are governed by stream function, vorticity, energy, nanoparticle concentration and number of bacteria. Radial basis function-based finite difference method for the space derivatives and the second-order backward differentiation formula for the time derivatives are performed. Numerical outputs in view of isolines as well as average Nusselt number, average Sherwood number and flux density of microorganisms are presented.

Convective mass transfer rises if any of Lewis number, Peclet number, Rayleigh number, bioconvection Rayleigh number and Brownian motion parameter increases, and the flux density of microorganisms is an increasing function of Rayleigh number, bioconvection Rayleigh number, Peclet number, Brownian and thermophoresis parameters. The rise in buoyancy ratio parameter between 0.1 and 1 and the rise in Hartmann number between 0 and 50 reduce all outputs average Nusselt, average Sherwood numbers and flux density of microorganisms.

This study implies the importance of the presence of magnetotactic bacteria and magnetite nanoparticles inside a host fluid in view of heat transfer and fluid flow. The limitation is to check the efficiency on numerical aspect. Experimental observations would be more effective.

In practical point of view, in a heat transfer and fluid flow system involving magnetite nanoparticles, the inclusion of magnetotactic bacteria and MF effect provide control over fluid flow and heat transfer.

This is a scientific study. However, this idea may be extended to sustainable energy or biofuel studies, too. This means that a better world may create better social environment between people.

The presence of magnetotactic bacteria inside a Fe₃O₄–water NF under the effect of a MF is a good controller on fluid flow and heat transfer. Since the magnetotactic bacteria is fed by nanoparticles Fe₃O₄ which has strong magnetic property, varying nanoparticle concentration and Brownian and thermophoresis effects are first considered.

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

This paper aims to report the study of a fluid buoy system that includes wave effects, with particular emphasis on validating the numerical results with experimental data.

A fluid–solid coupled algorithm is proposed to describe the motion of a rigid buoy under the effects of waves. The Navier–Stokes equations are solved with the open-source finite volume package Code Saturne, in which a free-surface capture technique and equations of motion for the solid are implemented. An ad hoc experiment on a laboratory scale is built. A buoy is placed into a tank partially filled with water; the tank is mounted into a shake table and subjected to controlled motion that promotes waves. The experiment allows for recording the evolution of the free surface at the control points using the ultrasonic sensors and the movement of the buoy by tracking the markers by postprocessing the recorded videos. The numerical results are validated by comparison with the experimental data.

The implemented free-surface technique, developed within the framework of the finite-volume method, is validated. The best-obtained agreement is for small amplitudes compatible with the waves evolving under deep-water conditions. Second, the algorithm proposed to describe rigid-body motion, including wave analysis, is validated. The numerical body motion and wave pattern satisfactorily matched the experimental data. The complete 3D proposed model can realistically describe buoy motions under the effects of stationary waves.

The novel aspects of this study encompass the implementation of a fluid–structure interaction strategy to describe rigid-body motion, including wave effects in a finite-volume context, and the reported free-surface and buoy position measurements from experiments. To the best of the authors’ knowledge, the numerical strategy, the validation of the computed results and the experimental data are all original contributions of this work.

High resolution simulations of body-internal electric field strengths induced by magneto-quasistatic fields from wireless power transfer systems are computationally expensive. The exposure simulation can be split into two separate simulation steps allowing the calculation of the magnetic flux density distribution, which serves as input into the second simulation step to calculate the body-internal electric fields. In this work, the magnetic flux density is interpolated from in situ measurements in combination with the scalar-potential finite difference scheme to calculate the resulting body-internal field. These calculations are supposed to take less than 5 s to achieve a near real-time visualization of these fields on mobile devices. The purpose of this work is to present an implementation of the simulation on graphics processing units (GPUs), allowing for the calculation of the body-internal field strength in about 3 s.

This work uses the co-simulation scalar-potential finite difference scheme to determine the body-internal electric field strength of human models with a voxel resolution of 2 × 2 × 2 mm³. The scheme is implemented on GPUs. This simulation scheme requires the magnetic flux density distribution as input, determined from radial basis functions.

Using NVIDIA A100 GPUs, the body-internal electric field strength with high-resolution models and 8.9 million degrees of freedom can be determined in about 2.3 s.

This paper describes in detail the used scheme and its implementation to make use of the computational performance of modern GPUs.