High-Order Bound-Preserving Finite Difference Methods for Multispecies and Multireaction Detonations

In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.

Stability and accuracy of central difference method for real-time dynamic substructure testing considering mass participation coefficient

For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.

Analysis of Extended Fisher-Kolmogorov Equation in 2D Utilizing the Generalized Finite Difference Method with Supplementary Nodes

In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.

Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

Physicochemical properties and biological activities of polysaccharides from the peel of Dioscorea opposita Thunb.extracted by four different methods

Polysaccharides are the important biologically active components found in the peel of Dioscorea opposita Thunb.(DTTP).The influences of 4 extraction methods,namely hot water extraction(W),acidic extraction(HA),hot-compressed water extraction(HCW)and enzyme-assisted extraction(EAE),on the yields,physicochemical properties,hypoglycemic and antioxidant activities of polysaccharides from DTTP were studied and compared.Among these DTTP polysaccharides,DTTP-HA was outstanding in several respects.DTTP-HA was the most water soluble;it had the highest total carbohydrate content(85.08%),the highest uronic acid(13.20%)and the highest thermal stability.DTTP-HA and DTTP-W possessed a triple-helix structure.These 4 kinds of polysaccharides have the same types of monosaccharides,but in different molar percentages.Extraction method had a significant impact on the microstructures of the extracted polysaccharides.DTT-HA exhibited irregular structure with many holes.Among the 4 extracted methods,the DTTP-HA and DTTP-W initially exhibited higher hypoglycemic and antioxidant activities.The better bioactivities of DTTP-HA may be related to the above factors.The findings indicated that acid extraction is an effective method to extract polysaccharides with high biological activities from DTTP.

Asymptotic Behavior of the Finite Difference and the Finite Element Methods for Parabolic Equations

The asymptotic convergence of the solution of the parabolic equation is proved. By the eigenvalues estimation, we obtain that the approximate solutions by the finite difference method and the finite element method are asymptotically convergent. Both methods are considered in continnous time.

Finite Difference Methods for the Time Fractional Advection-diffusion Equation

In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.

Solution of a One-Dimension Heat Equation Using Higher-Order Finite Difference Methods and Their Stability

One-dimensional heat equation was solved for different higher-order finite difference schemes, namely, forward time and fourth-order centered space explicit method, backward time and fourth-order centered space implicit method, and fourth-order implicit Crank-Nicolson finite difference method. Higher-order schemes have complexity in computing values at the neighboring points to the boundaries. It is required there a specification of the values of field variables at some points exterior to the domain. The complexity was incorporated using Hicks approximation. The convergence and stability analysis was also computed for those higher-order finite difference explicit and implicit methods in case of solving a one dimensional heat equation. The obtained numerical results were compared with exact solutions. It is found that backward time and fourth-order centered space implicit scheme along with Hicks approximation performed well over the other mentioned higher-order approaches.

Efficient Finite Difference/Spectral Method for the Time Fractional Ito Equation Using Fast Fourier Transform Technic

A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.

Numerical Study by Imposing the Finite Difference Method for Unsteady Casson Fluid Flow with Heat Flux

This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.

An Effective Meshless Approach for Inverse Cauchy Problems in 2D and 3D Electroelastic Piezoelectric Structures

In the past decade,notable progress has been achieved in the development of the generalized finite difference method(GFDM).The underlying principle of GFDM involves dividing the domain into multiple sub-domains.Within each sub-domain,explicit formulas for the necessary partial derivatives of the partial differential equations(PDEs)can be obtained through the application of Taylor series expansion and moving-least square approximation methods.Consequently,the method generates a sparse coefficient matrix,exhibiting a banded structure,making it highly advantageous for large-scale engineering computations.In this study,we present the application of the GFDM to numerically solve inverse Cauchy problems in two-and three-dimensional piezoelectric structures.Through our preliminary numerical experiments,we demonstrate that the proposed GFDMapproach shows great promise for accurately simulating coupled electroelastic equations in inverse problems,even with 3%errors added to the input data.

ONE-PARAMETER FINITE DIFFERENCE METHODS AND THEIR ACCELERATED SCHEMES FOR SPACE-FRACTIONAL SINE-GORDON EQUATIONS WITH DISTRIBUTED DELAY

This paper deals with numerical methods for solving one-dimensional(1D)and twodimensional(2D)initial-boundary value problems(IBVPs)of space-fractional sine-Gordon equations(SGEs)with distributed delay.For 1D problems,we construct a kind of oneparameter finite difference(OPFD)method.It is shown that,under a suitable condition,the proposed method is convergent with second order accuracy both in time and space.In implementation,the preconditioned conjugate gradient(PCG)method with the Strang circulant preconditioner is carried out to improve the computational efficiency of the OPFD method.For 2D problems,we develop another kind of OPFD method.For such a method,two classes of accelerated schemes are suggested,one is alternative direction implicit(ADI)scheme and the other is ADI-PCG scheme.In particular,we prove that ADI scheme can arrive at second-order accuracy in time and space.With some numerical experiments,the computational effectiveness and accuracy of the methods are further verified.Moreover,for the suggested methods,a numerical comparison in computational efficiency is presented.

Radiation Effect on Heat Transfer Analysis of MHD Flow of Upper Convected Maxwell Fluid between a Porous and a Moving Plate

The study in this manuscript aims to analyse the impact of thermal radiation on the two-dimensional magnetohydrodynamic flow of upper convected Maxwell(UCM)fluid between parallel plates.The lower plate is porous and stationary,while the top plate is impermeable and moving.The equations that describe the flow are transformed into non-linear ordinary differential equations with boundary conditions by employing similarity transformations.The Homotopy Perturbation Method(HPM)is then employed to approach the obtained non-linear ordinary differential equations and get an approximate analytical solution.The analysis includes plotting the velocity profile for different Reynolds number values and temperature distribution curves for distinct physical parameters such as Reynolds number,Deborah number,magnetic parameter,porosity parameter,radiation parameter,and Prandtl number.In the case of injection,the temporal profile declines with an increase in radiation parameter as the plates move away from each other,and an opposite trend is observed as plates move towards each other.Furthermore,the skin friction coefficient and heat transfer rate are analysed for the impact of these parameters using HPM.The numerical values obtained using HPM are compared using the classical finite difference method.The results show good agreement between the semi-analytical and numerical solutions.

Steady Natural Convection from a Vertical Hot Plate with Variable Radiation

The natural convection from a vertical hot plate with radiation and constant flux is studied numerically to know the velocity and temperature distribution characteristics over a vertical hot plate.The governing equations of the natural convection in two-dimension are solved with the implicit finite difference method,whereas the discretized equations are solved with the iterative relaxation method.The results show that the velocity and the temperature increase along the vertical wall.The influence of the radiation parameter in the boundary layer is significant in increasing the velocity and temperature profiles.The velocity profiles increase with the increase of the radiation parameter.The temperature profiles near the wall plate parallel each other due to the constant heat flux applied to the wall.The influence of the radiation parameter is significant either in velocity or temperature characteristics.At the same time,the effect of the Prandtl number greater than 0.71 is not sensitive to the velocity and temperature variations elsewhere.

Magnetohydrodynamic Conjugate Free Convective Heat Transfer Analysis of an Isothermal Horizontal Circular Cylinder with Temperature Dependent Viscosity

Heat transfers due to MHD-conjugate free convection from the isothermal horizontal circular cylinder while viscosity is a function of temperature is investigated. The governing equations of the flow and connected boundary conditions are made dimensionless using a set of non-dimensional parameters. The governing equations are solved numerically using the finite difference method. Numerical results are obtained for various values of viscosity variation parameter, Prandtl number, magnetic parameter, and conjugate conduction parameter for the velocity and the temperature within the boundary layer as well as the skin friction coefficients and heat transfer rate along the surface.

Zeno and the Wrong Understanding of Motion—A Philosophical-Mathematical Inquiry into the Concept of Finitude as a Peculiarity of Infinity

In contrast to the solutions of applied mathematics to Zeno’s paradoxes, I focus on the concept of motion and show that, by distinguishing two different forms of motion, Zeno’s apparent paradoxes are not paradoxical at all. Zeno’s paradoxes indirectly prove that distances are not composed of extensionless points and, in general, that a higher dimension cannot be completely composed of lower ones. Conversely, lower dimensions can be understood as special cases of higher dimensions. To illustrate this approach, I consider Cantor’s only apparent proof that the real numbers are uncountable. However, his widely accepted indirect proof has the disadvantage that it depends on whether there is another way to make the real numbers countable. Cantor rightly assumes that there can be no smallest number between 0 and 1, and therefore no beginning of counting. For this reason he arbitrarily lists the real numbers in order to show with his diagonal method that this list can never be complete. The situation is different if we start with the largest number between 0 and 1 (0.999…) and use the method of an inverted triangle, which can be understood as a special fractal form. Here we can construct a vertical and a horizontal stratification with which it is actually possible to construct all real numbers between 0 and 1 without exception. Each column is infinite, and each number in that column is the starting point of a new triangle, while each row is finite. Even in a simple sine curve, we experience finiteness with respect to the y-axis and infinity with respect to the x-axis. The first parts of this article show that Zeno’s assumptions contradict the concept of motion as such, so it is not surprising that this misconstruction leads to contradictions. In the last part, I discuss Cantor’s diagonal method and explain the method of an inverted triangle that is internally structured like a fractal by repeating this inverted triangle at each column. The consequence is that we encounter two very different methods of counting. Vertically it is continuous, horizontally it is discrete. While Frege, Tarski, Cantor, Gödel and the Vienna Circle tried to derive the higher dimension from the lower, a procedure that always leads to new contradictions and antinomies (Tarski, Russell), I take the opposite approach here, in which I derive the lower dimension from the higher. This perspective seems to fail because Tarski, Russell, Wittgenstein, and especially the Vienna Circle have shown that the completeness of the absolute itself is logically contradictory. For this reason, we agree with Hegel in assuming that we can never fully comprehend the Absolute, but only its particular manifestations—otherwise we would be putting ourselves in the place of the Absolute, or even God. Nevertheless, we can understand the Absolute in its particular expressions, as I will show with the modest example of the triangle proof of the combined horizontal and vertical countability of the real numbers, which I developed in rejection of Cantor’s diagonal proof. .

A Hybrid ESA-CCD Method for Variable-Order Time-Fractional Diffusion Equations

In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.

Improved finite difference method for pressure distribution of aerostatic bearing

An improved finite difference method （FDM）is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aerostatic bearings. A detailed theoretical analysis of the pressure distribution of the orifice-compensated aerostatic journal bearing is presented. The nonlinear dimensionless Reynolds equation of the aerostatic journal bearing is solved by the finite difference method. Based on the principle of flow equilibrium, a new iterative algorithm named the variable step size successive approximation method is presented to adjust the pressure at the orifice in the iterative process and enhance the efficiency and convergence performance of the algorithm. A general program is developed to analyze the pressure distribution of the aerostatic journal bearing by Matlab tool. The results show that the improved finite difference method is highly effective, reliable, stable, and convergent. Even when very thin gas film thicknesses （less than 2 Win）are considered, the improved calculation method still yields a result and converges fast.