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.

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.

Interaction between Topographic Conditions and Entrainment Rate in Numerical Simulations of Debris Flow

Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship between the topographic conditions and debris flow initiation is not well understood. This study compared the use of several entrainment rate equations in numerical simulations of debris flows to examine the effect of topographic conditions on the flow. One-dimensional numerical simulations were performed based on the shallow water equations and three entrainment rate equations were tested. These entrainment rate equations were based on the same idea that erosion and the deposition of debris flows occur via the difference between the equilibrium and current conditions of debris flows, while they differed in the expression of the concentration, channel angle, and sediment amount. The comparison was performed using a straight channel with various channel angles and a channel with a periodically undulating surface. The three entrainment rate equations gave different amounts of channel bed degradation and hydrographs for a straight channel with a channel angle greater than 21° when water was supplied from upstream at a steady rate. The difference was caused by the expression of the entrainment rate equations. For channels with little undulation, the numerical simulations gave results almost identical to those for straight channels with the same channel angle. However, for channels with large undulations, the hydrographs differed from those for straight channels with the same channel angle when the channel angle was less than 21°. Rapid erosion occurred and the hydrograph showed a significant peak, especially in cases using the entrainment equation expressed by channel angle. This was caused by the effects of the steep undulating sections, since the effect increased with the magnitude of the undulation, suggesting that a debris flow in an upstream area develops differently according to the topographic conditions. These results also inferred that numerical simulations of debris flow can differ depending on the spatial resolution of the simulation domain, as the resolution determines the reproducibility of the undulations.

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.

Well interference evaluation considering complex fracture networks through pressure and rate transient analysis in unconventional reservoirs

Severe well interference through complex fracture networks(CFNs)can be observed among multi-well pads in low permeability reservoirs.The well interference analysis between multi-fractured horizontal wells(MFHWs)is vitally important for reservoir effective development.Well interference has been historically investigated by pressure transient analysis,while it has shown that rate transient analysis has great potential in well interference diagnosis.However,the impact of complex fracture networks(CFNs)on rate transient behavior of parent well and child well in unconventional reservoirs is still not clear.To further investigate,this paper develops an integrated approach combining pressure and rate transient analysis for well interference diagnosis considering CFNs.To perform multi-well simulation considering CFNs,non-intrusive embedded discrete fracture model approach was applied for coupling fracture with reservoir models.The impact of CFN including natural fractures and frac-hits on pressure and rate transient behavior in multi-well system was investigated.On a logelog plot,interference flow and compound linear flow are two new flow regimes caused by nearby producers.When both NFs and frac-hits are present in the reservoir,frac-hits have a greater impact on well#1 which contains frac-hits,and NFs have greater impact on well#3 which does not have frac-hits.For all well producing circumstances,it might be challenging to see divergence during pseudosteady state flow brought on by frac-hits on the logelog plot.Besides,when NFs occur,reservoir depletion becomes noticeable in comparison to frac-hits in pressure distribution.Application of this integrated approach demonstrates that it works well to characterize the well interference among different multi-fractured horizontal wells in a well pad.Better reservoir evaluation can be acquired based on the new features observed in the novel model,demonstrating the practicability of the proposed approach.The findings of this study can help for better evaluating well interference degree in multi-well systems combing PTA and RTA,which can reduce the uncertainty and improve the accuracy of the well interference analysis based on both field pressure and rate data.

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.

Numerical analysis of plasma arc behaviors

Completely understanding the physical mechanisms of the plasma arc is critical to its application in welding of medium thickness plates. In this study, a mathematical model is developed to analyze the temperature, fluid flow, electromagnetic fields and pressure distribution in plasma arc welding. The correlations between the torch structure （ nozzle diameter） and the plasma are properties are analyzed qualitatively. The influence of the plasma gas flow rate on the plasma arc behavior is also simulated numerically. The temperature distribution and current density of the plasma are change greatly with a little variation of the nozzle diameter and^or the plasma gas flow rate. Compared to the tungsten-inert-gas arc with almost same conditions, the heat intensity, fluid velocity and pressure at the anode suoCace rise by one order of magnitude for a plasma arc. The analysis results lay solid foundation for effective usage ofplnsma arc welding.

Numerical Algorithms for Solving One Type of Singular Integro-Differential Equation Containing Derivatives of the Time Delay States

This study presents numerical algorithms for solving a class of equations that partly consists of derivatives of the unknown state at previous certain times, as well as an integro-differential term containing a weakly singular kernel. These equations are types of integro-differential equation of the second kind and were originally obtained from an aeroelasticity problem. One of the main contributions of this study is to propose numerical algorithms that do not involve transforming the original equation into the corresponding Volterra equation, but still enable the numerical solution of the original equation to be determined. The feasibility of the proposed numerical algorithm is demonstrated by applying examples in measuring the maximum errors with exact solutions at every computed nodes and calculating the corresponding numerical rates of convergence thereafter.

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.

OPTIMAL CONVERGENCE RATE OF THE LANDAU EQUATION WITH FRICTIONAL FORCE

The Cauchy problem of the Landau equation with frictional force is investigated. Based on Fourier analysis and nonlinear energy estimates, the optimal convergence rate to the steady state is obtained under some conditions on initial data.

Numerical Solution of Advection Diffusion Equation Using Semi-Discretization Scheme

Numerical diffusion and oscillatory behavior characteristics are averted applying numerical solutions of advection-diffusion equation are themselves immensely sophisticated. In this paper, two numerical methods have been used to solve the advection diffusion equation. We use an explicit finite difference scheme for the advection diffusion equation and semi-discretization on the spatial variable for advection-diffusion equation yields a system of ordinary differential equations solved by Euler’s method. Numerical assessment has been executed with specified initial and boundary conditions, for which the exact solution is known. We compare the solutions of the advection diffusion equation as well as error analysis for both schemes.

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.

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.

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.

基于数值分析法的水肥一体机设计研究

以进一步提升水肥一体机的智能性、数字性作业为目标,针对整机的灌施精准性要求展开设计。充分考虑水肥一体机的运用特点和结构组成,结合水肥的参数融合关系,建立用于水肥一体机准确控制作业的状态方程模型,设计以数值分析处理函数为核心的软件控制模块,搭建同步动作实施的硬件执行平台。展开多作物的灌施作业试验,结果表明:基于数值分析方法的水肥一体机系统架构设计合理,数值分析算法融入有效,整机试验的数值计算准确率可达98.00%以上;水肥混合均匀有度,整体管路灌施顺畅,灌施指令准确率可达99.00%以上,整机灌施效率较高,充分验证了数值分析方法应用的正确性与优越性,可以促进类似农机装备与高等数学多维度融合。

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.

Sensitivity analysis of influencing parameters in cavern stability

In order to analyze the stability of the underground rock structures,knowing the sensitivity of geomechanical parameters is important.To investigate the priority of these geomechanical properties in the stability of cavern,a sensitivity analysis has been performed on a single cavern in various rock mass qualities according to RMR using Phase 2.The stability of cavern has been studied by investigating the side wall deformation.Results showed that most sensitive properties are coefficient of lateral stress and modulus of deformation.Also parameters of Hoek-Brown criterion and r c have no sensitivity when cavern is in a perfect elastic state.But in an elasto-plastic state,parameters of Hoek-Brown criterion and r c affect the deformability;such effect becomes more remarkable with increasing plastic area.Other parameters have different sensitivities concerning rock mass quality(RMR).Results have been used to propose the best set of parameters for study on prediction of sidewall displacement.