Central Difference Scheme Research Articles

Monotonic unbounded schemes transformer (MUST): an approach to remove undershoots and overshoots in family of unbounded schemes using finite-volume method

PurposeThis paper presents a Monotonic Unbounded Schemes Transformer (MUST) approach to bound/monotonize (remove undershoots and overshoots) unbounded spatial differencing schemes automatically, and naturally. Automatically means the approach (1) captures the critical cell Peclet number when an unbounded scheme starts to produce physically unrealistic solution automatically, and (2) removes the undershoots and overshoots as part of the formulation without requiring human interventions. Naturally implies, all the terms in the discretization equation of the unbounded spatial differencing scheme are retained.Design/methodology/approachThe authors do not formulate new higher-order scheme. MUST transforms an unbounded higher-order scheme into a bounded higher-order scheme.FindingsThe solutions obtained with MUST are identical to those without MUST when the cell Peclet number is smaller than the critical cell Peclet number. For cell Peclet numbers larger than the critical cell Peclet numbers, MUST sets the nodal values to the limiter value which can be derived for the problem at-hand. The authors propose a way to derive the limiter value. The authors tested MUST on the central differencing scheme, the second-order upwind differencing scheme and the QUICK differencing scheme. In all cases tested, MUST is able to (1) capture the critical cell Peclet numbers; the exact locations when overshoots and undershoots occur, and (2) limit the nodal value to the value of the limiter values. These are achieved by retaining all the discretization terms of the respective differencing schemes naturally and accomplished automatically as part of the discretization process. The authors demonstrated MUST using one-dimensional problems. Results for a two-dimensional convection–diffusion problem are shown in Appendix to show generality of MUST.Originality/valueThe authors present an original approach to convert any unbounded scheme to bounded scheme while retaining all the terms in the original discretization equation.

A comparative study on numerical methods for Fredholm integro-differential equations of convection-diffusion problem with integral boundary conditions

This paper numerically solves Fredholm integro-differential equations with small parameters and integral boundary conditions. The solution of these equations has a boundary layer at the right boundary. A central difference scheme approximates the second-order derivative, a backward difference (upwind scheme) approximates the first-order derivative, and the trapezoidal rule is used for the integral term with a Shishkin mesh. It is shown that theoretically, the proposed scheme is uniformly convergent with almost first-order convergence. Further to improve the order of convergence from first order to second order, we use the post-processing and the hybrid scheme. Two numerical examples are computed to support the theoretical results.

A novel framework of the lattice Boltzmann model for multilayer shallow water systems

This study proposes a novel framework of the lattice Boltzmann model for multilayer shallow water equations, considering the mass and momentum exchanges between layers (LABMSWE+). Compared with the original LABMSWE model consisting of N two-dimensional lattice Boltzmann method for shallow water equation (LABSWE) models, the new model includes 1+N LABSWE models. The singular LABSWE model with unit relaxation time is introduced to update the total water depth, and thus, the layer water depths can be obtained explicitly through the fixed layer ratios. The N-layer LABSWE models with the multiple-relaxation-time operator evolve the layer velocities. These two modules are coupled by the total water depth and depth-averaged velocities. The constructed model avoids the freely variable layer thicknesses, which is considered as the main source of the instability. In addition, the mass exchanges enable this model to simulate vertical circulation flows, which are beyond the application of the LABMSWE model. Several numerical tests are then conducted to validate the proposed model. The results show that it exactly satisfies the C-property. In addition, the central difference scheme is more stable and accurate than the upwind and nonequilibrium schemes in the computing of the mass exchanges. The numerical results have an excellent agreement with analytical solutions and reference data, while some unstable and nonphysical results are obtained by the original LABMSWE model. Moreover, the computational time is about 40%–60% of that for the MIKE3, a finite volume solver for the three-dimensional shallow water equations by the Danish Hydraulic Institute.

Global Stability and Forced Response Analysis of Swirling Flows In Aviation Combustors

Abstract Swirling jets are canonical flow fields used to stabilize flames in many gas turbine combustion systems. These flows amplify background acoustic disturbances and perturb the flame, producing combustion noise and potentially combustion instabilities. Hydrodynamic flow response to acoustic excitations in reacting flows is an important modeling task for understanding combustion noise and combustion instabilities. This study utilizes a global hydrodynamic stability analysis in a Bi-Global and Tri-Global framework to model the vortical hydrodynamic modes of a reacting, swirling LES mean flow based on a commercial nozzle. After conducting an unforced, natural Bi-Global stability analysis to study the unstable eigenmodes of the flow, linearized equations of motion are used to study flow response of the flow to forcing. An empirical velocity transfer function is constructed by determining the flow response to a harmonically varying 150-1500 Hz u' velocity disturbance at the inflow boundary. These disturbances are strongly amplified, with amplification factors ranging from about 20 to 50 and peaking at about 1050 Hz, or st = 0.375. A comparison study with a cartesian, three-dimensional reacting base flow is also performed where the base-flow varies spatially in three-dimensions. The study utilizes a stable, high accuracy, but computationally effective centered finite difference scheme based on a three-dimensional structured mesh. Illustrative results present the qualitative modes and quantitative transfer functions for Bi-Global and Tri-Global stability analyses.

PARAMETER-UNIFORM HYBRID NUMERICAL SCHEME FOR SINGULARLY PERTURBED PARABOLIC CONVECTION-DIFFUSION PROBLEMS WITH DISCONTINUOUS INITIAL CONDITIONS

This article presents a parameter-uniform hybrid numerical technique for singularly perturbed parabolic convection-diffusion problems (SPPCDP) with discontinuous initial conditions (DIC). It utilizes the classical backward-Euler technique for time discretization and a hybrid finite difference scheme ( which is a proper combination of the midpoint upwind scheme in the outer regions and the classical central difference scheme in the interior layer regions(generated by the DIC)) for spatial discretization. The scheme produces parameter-uniform numerical approximations on a piecewise- uniform Shishkin mesh. When the perturbation parameter ε(0 < ε ≪ 1) is small, it becomes difficult to solve these problems using the classical numerical methods (standard central difference or a standard upwind scheme) on uniform meshes because discontinuous initial conditions frequently appear in the solutions of this class of problems. The method is shown to converge uniformly in the discrete supremum with nearly second-order spatial accuracy. The suggested method is subjected to a stability study, and parameter-uniform error estimates are generated. In order to support the theoretical findings, numerical results are presented.

Evaluating the Electromagnetic Field Problem Generated by the Power Transmission Line Using High Order Compact Method

This paper compares a finite difference method-based computational scheme for evaluating electromagnetic field problems in power transmission lines in cases of on surface and in Ground. The scheme uses Maxwell's partial differential equations to represent electric and magnetic field components and approximate boundary conditions. The High Order Compact Method (HOC) is applied to estimate the electric field (Ez )in ground and compares it with the standard central difference scheme. The HOC produces more accurate results than the traditional central approximation at 99.7% for the electric field Ez. It also calculates both of Electric and magnetic intensity to evaluate the effect of Electric and magnetic intensities in ground and on surface respect to distance.

CFD study of heat transfer in power‐law fluids over multiple corrugated circular cylinders in a heat exchanger

AbstractThe heat transfer in power‐law fluids across three corrugated circular cylinders placed in a triangular pitch arrangement is studied computationally in a confined channel. Continuity, momentum, and energy balance equations were solved using ANSYS FLUENT (Version 18.0). The flow is assumed to be steady, incompressible, two‐dimensional, and laminar. A square domain of side 300Dh is selected after a detailed domain study. An optimized grid with 98,187 cells is used in the study. The convergence criteria of 10−7 for the continuity, x‐momentum, and y‐momentum balances and 10−12 for the energy equation were used. Constant density and non‐Newtonian power‐law viscosity modules were used. The diffusive term is discretized using a central difference scheme. Convective terms are discretized using the Second‐Order Upwind scheme. Pressure–velocity coupling between continuity and momentum equations was implemented using the semi‐implicit method for pressure‐linked equation scheme. Streamlines show wake development behind the cylinders, which is very dominant at large ReN and n. Isotherm contours are cramped at higher values of ReN and PrN, implying higher heat transfer. Global parameters, like, Cd and Nu, are computed for the wide ranges of controlling dimensionless parameters, such as power‐law index (0.3 ≤ n ≤ 1.5), Reynolds (0.1 ≤ ReN ≤ 40), and Prandtl (0.72 ≤ PrN ≤ 500) numbers. The NuLocal plot attains a pitch near the corrugation of the surface due to abrupt changes in velocity and temperature gradients. Nu increases with ReN and/or PrN and decreases with n under ot herwise identical situations. Nu is correlated with pertinent parameters, namely, ReN, PrN, and n.

A linear second-order maximum bound principle preserving finite difference scheme for the generalized Allen–Cahn equation

In this paper, we consider the numerical method for generalized Allen–Cahn equation with nonlinear mobility and convection term. We propose a linear second-order finite difference scheme which preserves the discrete maximum bound principle (MBP). The scheme is discretized by stabilized Crank–Nicolson in time, upwind scheme for convection term and central-difference scheme for diffusion term. We show that the proposed scheme preserves the discrete MBP under some constraints on temporal step size and stabilizing parameter. Optimal L∞-error estimate is obtained for our scheme. Several numerical experiments are performed to validate our theoretical results.

A Posteriori Error Analysis of Defect Correction Method for Singular Perturbation Problems With Discontinuous Coefficient and Point Source

Abstract The paper presents a defect correction method to solve singularly perturbed problems with discontinuous coefficient and point source. The method combines an inexpensive, lower-order stable, upwind difference scheme and a higher-order, less stable central difference scheme over a layer-adapted mesh. The mesh is designed so that most mesh points remain in the regions with rapid transitions. A posteriori error analysis is presented. The proposed numerical method is analyzed for consistency, stability, and convergence. The error estimates of the proposed numerical method satisfy parameter-uniform second-order convergence on the layer-adapted grid. The convergence obtained is optimal because it is free from any logarithmic term. The numerical analysis confirms the theoretical error analysis.

Time-Stepping Error Estimates of Linearized Grünwald–Letnikov Difference Schemes for Strongly Nonlinear Time-Fractional Parabolic Problems

A fully discrete scheme is proposed for numerically solving the strongly nonlinear time-fractional parabolic problems. Time discretization is achieved by using the Grünwald–Letnikov (G–L) method and some linearized techniques, and spatial discretization is achieved by using the standard second-order central difference scheme. Through a Grönwall-type inequality and some complementary discrete kernels, the optimal time-stepping error estimates of the proposed scheme are obtained. Finally, several numerical examples are given to confirm the theoretical results.

Numerical approximation of parabolic singularly perturbed problems with large spatial delay and turning point

PurposeSingular perturbation turning point problems (SP-TPPs) involving parabolic convection–diffusion Partial Differential Equations (PDEs) with large spatial delay are studied in this paper. These type of equations are important in various fields of mathematics and sciences such as computational neuroscience and require specialized techniques for their numerical analysis.Design/methodology/approachWe design a numerical method comprising a hybrid finite difference scheme on a layer-adapted mesh for the spatial discretization and an implicit-Euler scheme on a uniform mesh in the temporal variable. A combination of the central difference scheme and the simple upwind scheme is used as the hybrid scheme.FindingsConsistency, stability and convergence are investigated for the proposed scheme. It is established that the present approach has parameter-uniform convergence of OΔτ+K−2(lnK)2, where Δτ and K denote the step size in the time direction and number of mesh-intervals in the space direction.Originality/valueParabolic SP-TPPs exhibiting twin boundary layers with large spatial delay have not been studied earlier in the literature. The presence of delay portrays an interior layer in the considered problem’s solution in addition to twin boundary layers. Numerical illustrations are provided to demonstrate the theoretical estimates.

Analysis of the dynamic response for Kirchhoff plates by the element-free Galerkin method

Dynamic response analysis involves examining structural behavior under dynamic loading conditions. Transient dynamic analysis emerges as a method employed to assess the responses of deformable bodies. Due to the depth analysis of transient responses for thin plates, the associated error analysis work becomes particularly important. In this work, we conduct stability and convergence analysis of the element-free Galerkin (EFG) method for the dynamic response of Kirchhoff plate model. We start by discretizing the thin plate dynamics using the EFG method for the spatial domain and the central difference scheme for the temporal domain. Stability and convergence analysis for transient responses, including deflection and moment responses, are derived similarly to those in two-dimensional transient structural dynamic analysis. The key to the analysis is transforming the governing equations and clamped boundary conditions into a weak formulation with penalty terms using the penalty method. The main results demonstrate a significant correlation between the error estimate and both nodal spacing and time step size. Additionally, the continuity of the approximation functions and the selection of penalty factors significantly affect numerical accuracy. To validate the theoretical results, we conduct numerical experiments involving clamped square and L-shaped plates with uniform and nonuniform node distributions, as well as a perforated plate model. Furthermore, we investigate the impact of node influence domains and penalty factors on numerical accuracy, facilitating parameter selection for practical engineering applications.

Application of Extended Eddy-viscosity and Elliptic-Relaxation Approaches to Turbulent Convective Flow in a Partially Divided Cavity

The paper reports on the numerical turbulence model in predicting mass, momentum and heat transfer in a partially divided cavity heated from the side using buoyancy-extended eddy-viscosity and elliptic relaxation approach with the algebraic expressions for the Reynold stress tensor and turbulent heat flux vector. The CDS (central differencing scheme) and LUDS (linear upwind differencing scheme) were used as the discretization method and the governing equations were solved using the finite volume method and Navier-Stokes solver. Validation of the model has been carried out by experimental data of convective flow in the cavity as well as by numerical data DNS (direct numerical simulation). The model agrees very well with the experiment and DNS and it is also able to demonstrate the performance which is comparable to that of the previous advanced second-moment closure model (SMC) in the literature. The results show that the model is suitable for use in simulations of the turbulent convective flow in a cavity with partition and it has the potential to be applied to more complex cavities and a wide range of turbulence levels.

A fast solvable operator-splitting scheme for time-dependent advection diffusion equation

It is well known that the implicit central difference discretization for unsteady advection diffusion equation (ADE) suffers from being time-consuming to solve when the advection term dominates. In this paper, we propose an operator-splitting scheme for the unsteady ADE, in which the ADE is firstly discretized by Crank-Nicolson (CN) scheme in time and central difference scheme in space; and then the discrete advection-diffusion problem is split as advection sub-problem and diffusion sub-problem at each time-level. The significance of the new scheme is that these sub-problems can be fast and directly solved within a linearithmic complexity (a linear-times-logarithm complexity) by means of fast sine transforms (FSTs). In particular, the complexity is independent of the dominance of the advection term. Theoretically, we show that proposed scheme is unconditionally stable and of second-order convergence in time and space. Numerical results are reported to show the efficiency of the proposed scheme.

Unveiling the Power of Implicit Six-Point Block Scheme: Advancing numerical approximation of two-dimensional PDEs in physical systems.

In the era of computational advancements, harnessing computer algorithms for approximating solutions to differential equations has become indispensable for its unparalleled productivity. The numerical approximation of partial differential equation (PDE) models holds crucial significance in modelling physical systems, driving the necessity for robust methodologies. In this article, we introduce the Implicit Six-Point Block Scheme (ISBS), employing a collocation approach for second-order numerical approximations of ordinary differential equations (ODEs) derived from one or two-dimensional physical systems. The methodology involves transforming the governing PDEs into a fully-fledged system of algebraic ordinary differential equations by employing ISBS to replace spatial derivatives while utilizing a central difference scheme for temporal or y-derivatives. In this report, the convergence properties of ISBS, aligning with the principles of multi-step methods, are rigorously analyzed. The numerical results obtained through ISBS demonstrate excellent agreement with theoretical solutions. Additionally, we compute absolute errors across various problem instances, showcasing the robustness and efficacy of ISBS in practical applications. Furthermore, we present a comprehensive comparative analysis with existing methodologies from recent literature, highlighting the superior performance of ISBS. Our findings are substantiated through illustrative tables and figures, underscoring the transformative potential of ISBS in advancing the numerical approximation of two-dimensional PDEs in physical systems.

Laminar convective heat transfer combining buoyancy and thermal radiation of Carreau fluid within a concentric and eccentric annulus

The present study concerns a steady-state laminar combined thermal radiation and buoyancy-driven conjugate magnetohydrodynamic (MHD) natural convective flow and heat transfer inside a fully porous medium eccentric circular annulus bounded by two infinite horizontal cylinders occupied by a Carreau fluid permeated by a uniform transverse magnetic field (B0 ) and inclined by an angle (α). The convective and thermal radiation are approximated by the Boussinesq and Roseland models, respectively. The numerical integration was performed using the implicit finite volume method. The convective terms in the momentum and energy equations were Discretized using the second-order upwind method. On the other hand, the diffusive terms were Discretized using the central difference scheme. The pressure-velocity coupling was solved using the SIMPLE approach. The findings suggest a clear trend of enhancement in the heat transfer rate when the aspect ratio and annulus eccentricity increase. This pattern is more pronounced in the negative vertical eccentricities. It is worth noting that the eccentric annulus, which has an annulus aspect ratio of 1.6 and a vertical eccentricity of −0.8, offers the best performance regarding the heat transfer mechanism inside the annulus gap. The heat transfer mechanism turns out to be strongly affected by thermal radiation, Hartmann number, and magnetic field inclination angle. The increased thermal radiation parameter and the decreased flow behavior index improve the heat transfer mechanism between the Carreau fluid and the inner cylinder in the annular gap. Furthermore, the average Nusselt number increases considerably as the induced magnetic field inclines.

Numerical study of shock wave boundary layer interaction flows using a novel turbulence model

This study proposes a novel partial average fluctuation velocity (PAFV) turbulence based on the PAFV and the average strain rate tensor, which is used to calculate the shock wave boundary layer interaction flows, including cascade flow and transonic compressor flow. In order to illustrate the adequacy, the calculation results of several other turbulence models were also compared with PAFV, including k-ω, shear stress transport, SA-neg (negative Spalart-Allmaras), and stress-BSL (ω-based Reynolds stress model). The suitability of the PAFV turbulence model was initially tested by calculating the transonic flow in the L030-4 cascade. Subsequently, numerical computations were performed for the 3D (three-dimensional) National Aeronautics and Space Administration Rotor 67 transonic compressor rotor. In the numerical simulation, the flow control equations were transformed in a general curvilinear coordinate system for precision enhancement, and the spatial discretization of control equations was executed using the finite difference method. The convection term was processed using the Steger–Warming flux vector splitting method, the diffusion term was discretized using a second-order central difference scheme, and the time derivative term was discretized using the third-order Runge–Kutta method. A parallel computing strategy with multi-block partitioning was employed to speed up the calculations and facilitate faster convergence alongside a local time step method. All numerical computation procedures were written in Python, revealing the following: (1) The PAFV turbulence model aptly simulates the transonic flow within the cascade channel. The shock wave pattern in the cascade channel aligns well with the experimental schlieren images. (2) During the Rotor 67 transonic compressor rotor calculation, the flow–pressure ratio, flow efficiency, and inlet and outlet total temperatures and pressures align well with the experimental data. Furthermore, the resolution of the flow field structure related to the interaction between the tip leakage vortex and channel shock wave is high. (3) The turbulence model uses PAFV as the turbulence velocity scale, thereby suppressing the fluctuation velocity and turbulence viscosity near shock wave areas, preventing excessive turbulence viscosity. Moreover, PAFV inherently exhibits transport properties and anisotropic characteristics, making it well-suited for handling transonic compressor flows.

On the Stability of a Central Difference Scheme with a Stabilizing Correction for the 3D Transport Equation

On the Stability of a Central Difference Scheme with a Stabilizing Correction for the 3D Transport Equation

Unstructured Conservative Level-Set (UCLS) method for reactive mass transfer in bubble swarms at high density ratio

This research presents a parallel Unstructured Conservative Level-Set (UCLS) method for reactive mass transfer in bubbles with a high-density ratio. This method uses the multiple-marker UCLS approach to circumvent the numerical coalescence of bubbles, combined with the finite-volume method to discretize transport equations on 3D collocated unstructured meshes. The fractional-step projection method solves the pressure-velocity coupling. The central difference scheme discretizes the diffusive term, whereas convective term within the momentum transport equation, level-set advection equations, and mass transfer equation are discretized by unstructured flux-limiters schemes. Such a combination of numerical techniques preserves the numerical stability in bubbly flows with a high Reynolds number and high-density ratio. Numerical and physical findings on the effect of physical properties ratios on the reactive mass transfer are reported, including their effect on the Reynolds number, interfacial area and Sherwood number.

Numerical simulation of coupled Klein–Gordon–Schrödinger equations: RBF partition of unity

The coupled Klein–Gordon–Schrödinger equations have significant implications in quantum field theory, particle physics, cosmology, and nonlinear dynamics. In this study, we propose an efficient method for numerically simulating this system. The proposed approach involves employing the radial basis function partition of unity for spatial discretization. This method utilizes scaled Lagrange basis functions with polyharmonic spline kernels, taking advantage of the scalability property of the polyharmonic kernel to ensure stability in the approximation process. The resulting spatially discretized system yields a nonlinear time-dependent set of differential equations. To solve this system, we combine an implicit central finite difference scheme with a predictor–corrector procedure to overcome the nonlinearity. In the numerical results section, we assess the efficiency, accuracy, and versatility of the proposed method by conducting several simulations in large domains over extended periods. We present and compare these results with those obtained using alternative methods, demonstrating the effectiveness of the proposed approach in accurately solving the coupled Klein–Gordon–Schrödinger equations.