Search results

This paper aims to tackle the problem of thermo‐solutal convection and macrosegregation during ingot solidification of metal alloys. Complex flow structures associated with the development of channels segregate and sharp gradients in the solutal field call for the implementation of accurate methods for numerical modeling of alloy solidification. In particular, the solute transport equation is convection dominated and requires special non‐oscillarity type high‐order schemes to handle the regions of channels segregates.

In the present study, a time‐splitting approach has been adopted to separately handle solute advection and diffusion. This splitting technique allows the application of accurate total variation dimensioning (TVD) schemes for solution of solute advection. Applications of second‐order Lax‐Wendroff TVD SUPERBEE and fifth‐order weighted essentially non‐oscillatory (WENO) schemes are described in the present article. Classical numerical solution of solute transport using hybrid and central‐difference schemes are also employed for the purpose of comparisons. Numerical simulations for solidification of Pb‐18%Sn in a two‐dimensional rectangular cavity have been carried out using different numerical schemes.

Numerical results show the difficulty of obtaining grid‐independent solutions with respect to local details in the region of channels. Grid convergence patterns and numerical uncertainty are found to be dependent on the applied scheme. In general, the first‐order hybrid scheme is diffusive and under predicts the formation of channels. The second‐order central‐difference scheme brings about oscillations with possible non‐physical extremes of solute composition in the region of channel segregates due to sharp gradients in the solutal field. The results obtained using TVD and WENO schemes contain no oscillations and show an excellent capture of channels formation and resolution of the interface between solute‐rich and depleted bands. Different stages of channels formation are followed by analyzing thermo‐solutal convection and macrosegregation at different times during solidification.

Accurate prediction of local variation in the solutal and flow fields in the channels regions requires grid refinement up to scales in the order of microscopic dendrite arm spacing. This imposes limitations in terms of large computational time and applicability of available macroscopic models based on classical volume‐averaging techniques.

The present study is very useful for numerical simulation of macrosegregation during ingot casting of metal alloys.

The paper provides the methodology and application of TVD schemes to predict channel segregates during columnar solidification of metal alloys. It also demonstrates the limitations of classical schemes for simulation of alloy solidification.

The performance of a spinning secant‐ogive‐cylinder‐boattail projectile in the transonic regime in terms of aerodynamic drag has been analyzed numerically in this study. To obtain an accurate prediction of the spinning effect on individual drag components and total drag of a projectile for the shell design, the implicit, diagonalied, symmetric Total Variation Diminishing (TVD) scheme, accompanied by a suitable grid, is employed to solve the thin‐layer axisymmetric Navier‐Stokes equations associated with the Baldwin‐Lomax turbulence model. The computed results show that, in comparison with the non‐spinning case, to increase the spin rate can result in increases in viscous drag and nose pressure drag, but can cause decreases in boattail drag and base drag. The variations of these drag components result in only a small (less than 5%) increase in total drag; thus the performance of the transonic projectiles is found to be insensitive to the spin rate.

Proposes a measure to stabilize the fourth(fifth)‐order high resolution schemes for the compressible Navier‐Stokes equations. Solves the N‐S equations of the volume fluxes and the low‐Reynolds number k‐ε turbulence model in general curvilinear co‐ordinates by the delta‐form implicit finite difference methods. Notes that, in order to simulate the flow containing weak discontinuities accurately, it is very effective to use some higher‐order TVD upstream‐difference schemes in the right‐hand side of the equations of these methods; however, the higher‐order correction terms of such schemes in general amplify the numerical disturbances. Therefore, restricts these terms here by operating the minmod functions to the curvatures so as to suppress the occurrence of new inflection points. Computes an unsteady transonic turbine cascade flow where vortex streets occur from the trailing edge of blades and interact with shock waves. Finds that the stabilization measure improves not only the computational results but also the convergency for such a complicated flow problem.

The purpose of this paper is to present an efficient improved version of Implicit Pressure-Explicit Saturation (IMPES) method for the solution of incompressible two-phase flow model based on the discontinuous Galerkin (DG) numerical scheme.

The governing equations, based on the wetting-phase pressure-saturation formulation, are discretized using various primal DG schemes. The authors use H(div) velocity reconstruction in Raviart-Thomas space (RT_0 and RT_1), the weighted average formulation, and the scaled penalties to improve the spatial discretization. It uses a new improved IMPES approach, by using the second-order explicit Total Variation Diminishing Runge-Kutta (TVD-RK) as temporal discretization of the saturation equation. The main purpose of this time stepping technique is to speed up computation without losing accuracy, thus to increase the efficiency of the method.

Utilizing pressure internal interpolation technique in the improved IMPES scheme can reduce CPU time. Combining the TVD property with a strong multi-dimensional slope limiter namely, modified Chavent-Jaffre leads to a non-oscillatory scheme even in coarse grids and highly heterogeneous porous media.

The presented locally conservative scheme can be applied only in 2D incompressible two-phase flow modeling in non-deformable porous media. In addition, the capillary pressure discontinuity between two adjacent rock types assumed to be negligible.

The proposed numerical scheme can be efficiently used to model the incompressible two-phase flow in secondary recovery of petroleum reservoirs and tracing immiscible contamination in aquifers.

The paper describes a novel version of the DG two-phase flow which illustrates the effects of improvements in special discretization. Also the new improved IMPES approach used reduces the computation time. The non-oscillatory scheme is an efficient algorithm as it maintains accuracy and saves computation time.

The purpose of this study is to propose a novel relaxation modified factor to accelerate the numerical solution of the radiative transfer equation (RTE) with several high-resolution total variation diminishing schemes. The methodology proposed is denoted as the X-factor method.

The X-factor method was compared with the technique deferred-correction (DC) for the calculations of a two-dimensional cavity with absorting-emiting-scatteting gray media using the discrete ordinates method. Four parameters were considered to evaluate: the absorption coefficient, the emissivity of boundary surface, the scattering albedo and under-relaxation factor.

The results showed the central processing unit (CPU) time of X-factor method was lower than DC. The reductions of CPU time with the X-factor method were observed from 0.6 to 75.4%.

The superiority of the X-factor method over DC was showed with the reduction of CPU time of the numerical solution of RTE for evaluated cases.

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 evaluate the performance of a numerical method for the solution to shallow-water equations on a staggered grid, in simulations for shear instabilities at two convective Froude numbers.

The simulations start from a small perturbation to a base flow with a hyperbolic-tangent velocity profile. The subsequent development of the shear instabilities is studied from the simulations using a number of flux-limiting schemes, including the second-order MINMOD, the third-order ULTRA-QUICK and the fifth-order WENO schemes for the spatial interpolation of the nonlinear fluxes. The fourth-order Runge-Kutta method advances the simulation in time.

The simulations determine two parameters: the fractional growth rate of the linear instabilities; and the vorticity thickness of the first nonlinear peak. Grid refinement using 32, 64, 128, 256 and 512 nodes over one wave length determines the exact values by extrapolation and the computational error for the parameters. It also determines the overall order of convergence for each of the flux-limiting schemes used in the numerical simulations.

The four-digit accuracy of the numerical simulations presented in this paper are comparable to analytical solutions. The development of this reliable numerical simulation method has paved the way for further study of the instabilities in shear flows that radiate waves.

The purpose of this paper is to develop and apply a high-order numerical method based on finite volume approximation for quadratic integrate-and-fire (QIF) neuron model with the help of population density approach.

The authors apply the population density approach for the QIF neuron model to derive the governing equation. The resulting mathematical model cannot be solved with existing analytical or numerical techniques owing to the presence of delay and advance. The numerical scheme is based along the lines of approximation: spatial discretization is performed by weighted essentially non-oscillatory (WENO) finite volume approximation (FVM) and temporal discretization are performed by strong stability-preserving explicit Runge–Kutta (SSPERK) method. Compared with existing schemes of orders 2 and 3 from the literature, the proposed scheme is found to be more efficient and it produces accurate solutions with few grid cells. In addition to this, discontinuity is added in the application of the model equation to illustrate the high performance of the proposed scheme.

The developed scheme works nicely for the simulation of the resulting model equation. The authors discussed the role of inhibitory and excitatory parts in variation of neuronal firing. The validation of the designed scheme is measured by its comparison with existing schemes in the literature. The efficiency of the designed scheme is demonstrated via numerical simulations.

It is expected that the present study will be a useful tool to tackle the complex neuron model and related studies.

The novel aspect of this paper is the application of the numerical methods to study the modified version of leaky integrate-and-fire neuron based on a QIF neuron. The model of the current study is inspired from the base model given in Stein (1965) and modified version in Kadalbajoo and Sharma (2005) and Wang and Zhang (2014). The applicability was confirmed by taking some numerical examples.

The purpose of this paper is to provide an improved Mie-Grüneisen mixture model to simulate underwater explosion (UNDEX).

By using Mie-Grüneisen equations of state (EOS) to model explosive charge, liquid water and solid structure, the whole fluid field is considered as a multi-phases mixture under Mie-Grüneisen EOS. Then by introducing auxiliary variables in Eulerian model and using mass fraction to establish a diffusion balance, a new improved Mie-Grüneisen mixture model is presented here. For the new reconstructed mixture model, a second order MUSCL scheme with TVD limiter is employed to solve the multi-phase Riemann problem.

Numerical examples show that the results obtained by Mie-Grüneisen mixture model are quite closed to theoretical and empirical data. The model can be also used in 2-D fluid-structure problem of UNDEX effectively and it is proved that the deformation of structure can be clearly described by mass fraction.

The FVM model based on mass fraction can only describe the motion of compressible material under impact. Material failure or large deformation needs a modification about the EOS or implementations of other models (i.e. FEM model).

An improved non-oscillation Mie-Grüneisen mixture model, which based on mass fraction, is given in the present paper. The present Mie-Grüneisen mixture model provides a simplified and efficient way to simulate UNDEX. The feasibility of this model to simulate the detonation impacts on different mediums, including water and other metal mediums, is tested and verified here. Then the model is applied to the simulation of underwater contact explosion problem. In the simulation, deformation of structure under explosion loads, as well as second shock wave, are studied here.

A consistent implementation of the general computational framework of unified second-order time accurate integrators via the well-known GSSSS framework in conjunction with the traditional Finite Difference Method is presented to improve the numerical simulations of reactive two-phase flows.

In the present paper, the phase interaction evaluation in the present implementation of the reactive two-phase flows has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process for solving the two phase flow dynamics with reactions.

Numerical examples, including the classical Sod shock tube problem and a reactive two-phase flow problem, are exploited to validate the proposed time integration framework and families of algorithms consistently to second order in time accuracy; this is in contrast to the traditional practices which only seem to obtain first-order time accuracy because of the inconsistent time level implementation with respect to the interaction of two phases. The comparisons with the traditional implementation and the advantages of the proposed implementation are given in terms of the improved numerical accuracy in time. The proposed approaches provide a correct numerical simulation implementation to the reactive two-phase flows and can obtain better numerical stability and computational features.

The new algorithmic framework and the consistent time level evaluation extended with the GS4 family encompasses a multitude of past and new schemes and offers a general purpose and unified implementation for fluid dynamics.