Lie group analysis, numerical and non-traveling wave solutions for the (2+1)-dimensional diffusion–advection equation with variable coefficients

In this paper, the variable-coefficient diffusion-advection （DA） equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by determining the complete sets of point symmetries of this equation, and then exact and numerical solutions are reported for the reduced second-order nonlinear ordinary differential equations. Further, an extended （Gl/G）-expansion method is applied to the DA equation to construct some new non-traveling wave solutions.

Stability Analysis of a Numerical Integrator for Solving First Order Ordinary Differential Equation

In this paper, we used an interpolation function to derive a Numerical Integrator that can be used for solving first order Initial Value Problems in Ordinary Differential Equation. The numerical quality of the Integrator has been analyzed to authenticate the reliability of the new method. The numerical test showed that the finite difference methods developed possess the same monotonic properties with the analytic solution of the sampled Initial Value Problems.

A Comparative Study of Adomian Decomposition Method with Variational Iteration Method for Solving Linear and Nonlinear Differential Equations

This study compares the Adomian Decomposition Method (ADM) and the Variational Iteration Method (VIM) for solving nonlinear differential equations in engineering. Differential equations are essential for modeling dynamic systems in various disciplines, including biological processes, heat transfer, and control systems. This study addresses first, second, and third-order nonlinear differential equations using Mathematica for data generation and graphing. The ADM, developed by George Adomian, uses Adomian polynomials to handle nonlinear terms, which can be computationally intensive. In contrast, VIM, developed by He, directly iterates the correction functional, providing a more straightforward and efficient approach. This study highlights VIM’s rapid convergence and effectiveness of VIM, particularly for nonlinear problems, where it simplifies calculations and offers direct solutions without polynomial derivation. The results demonstrate VIM’s superior efficiency and rapid convergence of VIM compared with ADM. The VIM’s minimal computational requirements make it practical for real-time applications and complex system modeling. Our findings align with those of previous research, confirming VIM’s efficiency of VIM in various engineering applications. This study emphasizes the importance of selecting appropriate methods based on specific problem requirements. While ADM is valuable for certain nonlinearities, VIM’s approach is ideal for many engineering scenarios. Future research should explore broader applications and hybrid methods to enhance the solution’s accuracy and efficiency. This comprehensive comparison provides valuable guidance for selecting effective numerical methods for differential equations in engineering.

Numerical Analysis of a Spiral-groove Dry-gas Seal Considering Micro-scale Effects

A dry-gas seal system is a non-contact seal technology that is widely used in different industrial applications.Spiral-groove dry-gas seal utilizes fluid dynamic pressure effects to realize the seal and lubrication processes,while forming a high pressure gas film between two sealing faces due to the deceleration of the gas pumped in or out.There is little research into the effects and the influence on seal performance,if the grooves and the gas film are at the micro-scale.This paper investigates the micro-scale effects on spiral-groove dry-gas seal performance in a numerical solution of a corrected Reynolds equation.The Reynolds equation is discretized by means of the finite difference method with the second order scheme and solved by the successive-over-relaxation（SOR） iterative method.The Knudsen number of the flow in the sealing gas film is changed from 0.005 to 0.120 with a variation of film depth and sealing pressure.The numerical results show that the average pressure in the gas film and the sealed gas leakage increase due to micro-scale effects.The open force is enlarged,while the gas film stiffness is significantly decreased due to micro-scale effects.The friction torque and power consumption remain constant,even in low sealing pressure and spin speed conditions.In this paper,the seal performance at different rotor face spin speeds is also described.The proposed research clarifies the micro-scale effects in a spiral-groove dry-gas seal and their influence on seal performance,which is expected to be useful for the improvement of the design of dry-gas seal systems operating in the slip flow regime.

Approximate Solutions of Primary Resonance for Forced Duffing Equation by Means of the Homotopy Analysis Method

Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approximate methods.Classical perturbation methods such as LP method,KBM method,multi-scale method and the averaging method on weakly nonlinear vibration system is effective,while the strongly nonlinear system is difficult to apply.Approximate solutions of primary resonance for forced Duffing equation is investigated by means of homotopy analysis method （HAM）.Different from other approximate computational method,the HAM is totally independent of small physical parameters,and thus is suitable for most nonlinear problems.The HAM provides a great freedom to choose base functions of solution series,so that a nonlinear problem may be approximated more effectively.The HAM provides us a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter and the auxiliary function.Therefore,HAM not only may solve the weakly non-linear problems but also may be suitable for the strong non-linear problem.Through the approximate solution of forced Duffing equation with cubic non-linearity,the HAM and fourth order Runge-Kutta method of numerical solution were compared,the results show that the HAM not only can solve the steady state solution,but also can calculate the unsteady state solution,and has the good computational accuracy.

On Trigonometric Numerical Integrator for Solving First Order Ordinary Differential Equation

In this paper, we used an interpolation function with strong trigonometric components to derive a numerical integrator that can be used for solving first order initial value problems in ordinary differential equation. This numerical integrator has been tested for desirable qualities like stability, convergence and consistency. The discrete models have been used for a numerical experiment which makes us conclude that the schemes are suitable for the solution of first order ordinary differential equation.

A ONE-STEP EXPLICIT FORMULA FOR THE NUMERICAL SOLUTION OF STIFF ORDINARY DIFFERENTIAL EQUATION

In this paper, a new one-step explicit method of fourth order is derived. The new method is proved to be A-stable and L-stable, and it gives exact results when applied to the test equation y’=λy with Re(λ)【0, Also several numerical examples are included.

EXPONENTIAL DICHOTOMIES OF NONLINEAR DISCRETE SYSTEMS AND ITS APPLICATION TO NUMERICAL ANALYSIS AND COMPUTATION

In this paper we consider the upwind difference scheme of a kind of boundary value problems for nonlinear, second order, ordinary differential equations. Singular perturbation method is applied to construct the asymptotic approximation of the solution to the upwind difference equation. Using the theory of exponential dichotomies we show that the solution of an order-reduced equation is a good approximation of the solution to the upwind difference equation except near boundaries. We construct correctors which yield asymptotic approximations by adding them to the solution of the order-reduced equation. Finally, some numerical examples are illustrated.

A Third-Order Scheme for Numerical Fluxes to Guarantee Non-Negative Coefficients for Advection-Diffusion Equations

According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. In case of advection-diffusion equations, so far there have been not found stable schemes with positive difference coefficients in a family of numerical schemes exceeding the second-order accuracy. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter by using the same stencil number as convemtional third-order shemes such as KAWAMURA and UTOPIA schemes. We extend the present method into multi-dimensional equations. Numerical experiments for linear and nonlinear advection-diffusion equations were performed and the present scheme’s applicability to nonlinear Burger’s equation was confirmed.

Numerical Analysis of the Magnetization Behavior in Magnetic Resonance Imaging in the Presence of Multiple Chemical Exchange Pools

The purpose of this study was to demonstrate a simple and fast method for solving the time-dependent Bloch-McConnell equations describing the behavior of magnetization in magnetic resonance imaging (MRI) in the presence of multiple chemical exchange pools. First, the time-dependent Bloch- McConnell equations were reduced to a homogeneous linear differential equation, and then a simple equation was derived to solve it using a matrix operation and Kronecker tensor product. From these solutions, the longitudinal relaxation rate (R1ρ) and transverse relaxation rate in the rotating frame (R2ρ) and Z-spectra were obtained. As illustrative examples, the numerical solutions for linear and star-type three-pool chemical exchange models and linear, star- type, and kite-type four-pool chemical exchange models were presented. The effects of saturation time (ST) and radiofrequency irradiation power (ω1) on the chemical exchange saturation transfer (CEST) effect in these models were also investigated. Although R1ρ and R2ρ were not affected by the ST, the CEST effect observed in the Z-spectra increased and saturated with increasing ST. When ω1 was varied, the CEST effect increased with increasing ω1 in R1ρ, R2ρ, and Z-spectra. When ω1 was large, however, the spillover effect due to the direct saturation of bulk water protons also increased, suggesting that these parameters must be determined in consideration of both the CEST and spillover effects. Our method will be useful for analyzing the complex CEST contrast mechanism and for investigating the optimal conditions for CEST MRI in the presence of multiple chemical exchange pools.

Finite Difference Schemes for Time-Space Fractional Diffusion Equations in One-and Two-Dimensions

In this paper,finite difference schemes for solving time-space fractional diffusion equations in one dimension and two dimensions are proposed.The temporal derivative is in the Caputo-Hadamard sense for both cases.The spatial derivative for the one-dimensional equation is of Riesz definition and the two-dimensional spatial derivative is given by the fractional Laplacian.The schemes are proved to be unconditionally stable and convergent.The numerical results are in line with the theoretical analysis.

Numerical analysis of the return flow solar air heater(RF-SAH)with assimilation of V-type artificial roughness

A novel design of Return Flow Solar Air Heater(RFSAH)with different arrangements of baffles especially V-Type Artificial roughness is simulated and numerically analyzed with energy balance equations.To enhance the effectiveness of baffles,numerous studies have been conducted.The performance of the RFSAH is studied in terms of thermal efficiency,thermo-hydraulic efficiency,and optimization of baffle parameters.Maximum Thermal efficiency and thermo-hydraulic efficiency are found in RFSAH with baffle on both sides of the absorber plate and mass flow rate above 0.2kg/s.Sensitivity analysis of the influencing parameters is carried out and reported the best performance of the system on selective geometrical parameters(ψ=0.7,β=20%,e/H=1,p/e=0.8,α=60°).The results obtained from the present model are validated with the published experimental results and have been found in quite reasonable agreement with an average error of 16.45%.Thermal and Thermohydraulic efficiency of RFSAH with a baffle on both sides of the absorber plate is maximum among baffles below,above,and on both sides of the absorber plate.It is observed that the thermal efficiency of RFSAH is greater than SF-SAH.The proposed optimum baffles roughness is suggested to increase the air upholding time period for more efficient output.

A NUMERICAL AND ANALYTICAL STUDY ON A TAIL-FLAPPING MODEL FOR FISH FAST C-START

The force production physics and the flow control mechanism of fish fast C-start are studied numerically and theoretically by using a tail-flapping model.The problem is simplified to a 2-D foil that rotates rapidly to and fro on one side about its fixed leading edge in water medium.The study involves the simulation of the flow by solving the two-dimensional unsteady incompressible Navier- Stokes equations and employing a theoretical analytic modeling approach.Firstly,reasonable thrust magnitude and its time history are obtained and checked by fitting predicted results coming from these two approaches.Next,the flow fields and vortex structures are given,and the propulsive mechanism is interpreted.The results show that the induction of vortex distributions near the trailing edge of the tail are important in the time-averaged thrust generation,though the added inertial effect plays an important role in producing an instant large thrust especially in the first stage.Furthermore,dynamic and energetic effects of some kinematic controlling factors are discussed.For enhancing the time- averaged thrust but keeping a favorable ratio of it to time-averaged input power within the limitations of muscle ability,it is recommended to have a larger deflection amplitude in a limited time interval and with no time delay between the to-and-fro strokes.

Numerical Analysis of a High-Order Scheme for Nonlinear Fractional Differential Equations with Uniform Accuracy

We introduce a high-order numerical scheme for fractional ordinary differential equations with the Caputo derivative.The method is developed by dividing the domain into a number of subintervals,and applying the quadratic interpolation on each subinterval.The method is shown to be unconditionally stable,and for general nonlinear equations,the uniform sharp numerical order 3−νcan be rigorously proven for sufficiently smooth solutions at all time steps.The proof provides a gen-eral guide for proving the sharp order for higher-order schemes in the nonlinear case.Some numerical examples are given to validate our theoretical results.

An Effective Numerical Method for the Solution of a Stochastic Coronavirus(2019-nCovid)Pandemic Model

Nonlinear stochastic modeling plays a significant role in disciplines such as psychology,finance,physical sciences,engineering,econometrics,and biological sciences.Dynamical consistency,positivity,and boundedness are fundamental properties of stochastic modeling.A stochastic coronavirus model is studied with techniques of transition probabilities and parametric perturbation.Well-known explicit methods such as Euler Maruyama,stochastic Euler,and stochastic Runge–Kutta are investigated for the stochastic model.Regrettably,the above essential properties are not restored by existing methods.Hence,there is a need to construct essential properties preserving the computational method.The non-standard approach of finite difference is examined to maintain the above basic features of the stochastic model.The comparison of the results of deterministic and stochastic models is also presented.Our proposed efficient computational method well preserves the essential properties of the model.Comparison and convergence analyses of the method are presented.

Maximum modulus principle estimates for one dimensional fractional diffusion equation

We present scheme I for solving one-dimensional fractional diffusion equation with variable coefficients based on the maximum modulus principle and two Grunwald approxima- tions. Scheme II is obtained by using classic Crank-Nicolson approximations in order to improve the time convergence. Schemes are proved to be unconditionally stable and second-order accuracy in spatial grid size for the problem with order of fractional derivative belonging to [（√17- 1）/2, 2] using the maximum modulus principle. A numerical example is given to confirm the theoretical analysis result.

Stochastic Approximate Solutions of Stochastic Differential Equations with Random Jump Magnitudes and Non-Lipschitz Coefficients

A class of stochastic differential equations with random jump magnitudes( SDEwRJMs) is investigated. Under nonLipschitz conditions,the convergence of semi-implicit Euler method for SDEwRJMs is studied. The main purpose is to prove that the semi-implicit Euler solutions converge to the true solutions in the mean-square sense. An example is given for illustration.

Solutions for Series of Exponential Equations in Terms of Lambert-W Function and Fundamental Constants

Series of exponential equations in the form of were solved graphically, numerically and analytically. The analytical solution was derived in terms of Lambert-W function. A general numerical solution for any y is found in terms of n or in base y. A solution is close to the fine structure constant. The equation which provided the solution as the fine structure constant was derived in terms of the fundamental constants.

Modelling of Energy Storage Photonic Medium by Wavelength-Based Multivariable Second-Order Differential Equation

Wavelength-dependent mathematical modelling of the differential energy change of a photon has been performed inside a proposed hypothetical optical medium.The existence of this medium demands certain mathematical constraints,which have been derived in detail.Using reverse modelling,a medium satisfying the derived conditions is proven to store energy as the photon propagates from the entry to exit point.A single photon with a given intensity is considered in the analysis and hypothesized to possess a definite non-zero probability of maintaining its energy and velocity functions analytic inside the proposed optical medium,despite scattering,absorption,fluorescence,heat generation,and other nonlinear mechanisms.The energy and velocity functions are thus singly and doubly differentiable with respect to wavelength.The solution of the resulting second-order differential equation in two variables proves that energy storage or energy flotation occurs inside a medium with a refractive index satisfying the described mathematical constraints.The minimum-value-normalized refractive index profiles of the modelled optical medium for transformed wavelengths both inside the medium and for vacuum have been derived.Mathematical proofs,design equations,and detailed numerical analyses are presented in the paper.

NUMERICAL ANALYSIS FOR A MEAN-FIELD EQUATION FOR THE ISING MODEL WITH GLAUBER DYNAMICS

In this paper, a mean-field equation of motion which is derived by Penrose (1991) for the dynamic Ising model with Glauber dynamics is considered. Various finite difference schemes such as explicit Euler scheme, predictor-corrector scheme and some implicit schemes are given and their convergence, stabilities and dynamical properties are discussed. Moreover, a Lyapunov functional for the discrete semigroup {S}(n>0) is constructed. Finally, numerical examples are computed and analyzed. it shows that the model is a better approximation to Cahn-Allen equation which is mentioned in Penrose (1991).