Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

基于蒙特卡洛法的Euler-Bernoulli梁基频和振型求解方法

将Rayleigh法和蒙特卡洛法相结合,在Euler-Bernoulli梁理论假设下求解了均匀梁、变截面梁和附带集中质量的变截面梁自由振动问题.对原本连续的梁结构模型进行离散化处理,利用蒙特卡洛法给出梁结构的假设振型.将假设得到的梁结构振型函数代入Rayleigh法,多次计算过程中,将历次基频所得值与计算所得最小值进行比较,根据其相对误差判断是否满足收敛条件,进而求得基频及对应的振型.结果表明,不同计算模型中基频最大误差不超过10%,能够满足工程需求,且精度和时间的控制参数调整灵活,使用者可根据自身需要自行调节.该方法理论简明,适用范围广泛,能够快速准确地求解诸多类型的梁结构基频和振型.

On the Approximation of Fractal-Fractional Differential Equations Using Numerical Inverse Laplace Transform Methods

Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.

扩散系数Holder连续的随机微分方程的截断Euler-Maruyama方法

本文研究漂移系数超线性增长和扩散系数Holder连续的随机微分方程的截断Euler-Maruyama方法的强收敛性.研究结果显示强收敛率依赖于Holder指数.本文给出一个例子验证所得的结果.

随机年龄结构固定资产系统倒向Euler法的p阶矩耗散性

在单边Lipschitz条件下,研究了一类随机年龄结构固定资产系统倒向Euler法数值解的p阶矩耗散性.当0<p<1时,步长满足一定条件可以得到系统的p阶矩耗散性;而p=2时,在对步长没有任何限制条件的情况下,得到了系统的均方耗散性.最后,通过数值例子验证了理论结果的可行性和有效性.

Analysis of regular and chaotic dynamics of the Euler-Bernoulli beams using finite difference and finite element methods

Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.

能量泛函及Euler-Lagrange方程在图像降噪中的应用研究

研究了自适应分数阶偏微分方程修正模型的能量泛函及Euler-Lagrange方程。首先,定义了自适应分数阶偏微分方程修正模型的能量泛函,其中包含未知函数和拉格朗日乘子的集合。然后,通过求解能量泛函的极值方程,推导出了Euler-Lagrange方程。最后,讨论了Euler-Lagrange方程在自适应分数阶偏微分方程修正模型中的应用。

The Semi-implicit Euler Method for Stochastic Pantograph Equations with Jumps

在这份报纸，我们在场半含蓄的 Euler (SIE ) 随机的缩放仪方程与的数字答案跳并且证明 SIE 近似答案在本地 Lipschitz 条件下面在吝啬平方的意义收敛到准确答案。

High Order Semi-implicit Multistep Methods for Time-Dependent Partial Differential Equations

We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not possible.As shown in Boscarino et al.(J.Sci.Comput.68:975-1001,2016)for Runge-Kutta methods,these semi-implicit techniques give a great flexibility,and allow,in many cases,the construction of simple linearly implicit schemes with no need of iterative solvers.In this work,we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype lineal'advection-diffusion equation and in the setting of strong stability preserving(SSP)methods.Our findings are demonstrated on several examples,including nonlinear reaction-diffusion and convection-diffusion problems.

A Comparative Study on Numerical Solutions of Initial Value Problems (IVP) for Ordinary Differential Equations (ODE) with Euler and Runge Kutta Methods

This paper mainly presents Euler method and fourth-order Runge Kutta Method (RK4) for solving initial value problems (IVP) for ordinary differential equations (ODE). The two proposed methods are quite efficient and practically well suited for solving these problems. In order to verify the ac-curacy, we compare numerical solutions with the exact solutions. The numerical solutions are in good agreement with the exact solutions. Numerical comparisons between Euler method and Runge Kutta method have been presented. Also we compare the performance and the computational effort of such methods. In order to achieve higher accuracy in the solution, the step size needs to be very small. Finally we investigate and compute the errors of the two proposed methods for different step sizes to examine superiority. Several numerical examples are given to demonstrate the reliability and efficiency.

Comparative Study on Results of Euler,Improved Euler and Run­ge-Kutta Methods for Solving the Engineering Unknown Problems

The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods are quite efficient and practically well suited for solving the unknown engineering problems.This paper aims to enhance the teaching and learning quality of teachers and students for various levels.At each point of the interval,the value of y is calculated and compared with its exact value at that point.The next interesting point is the observation of error from those methods.Error in the value of y is the difference between calculated and exact value.A mathematical equation which relates various functions with its derivatives is known as a differential equation.It is a popular field of mathematics because of its application to real-world problems.To calculate the exact values,the approximate values and the errors,the numerical tool such as MATLAB is appropriate for observing the results.This paper mainly concentrates on identifying the method which provides more accurate results.Then the analytical results and calculates their corresponding error were compared in details.The minimum error directly reflected to realize the best method from different numerical methods.According to the analyses from those three approaches,we observed that only the error is nominal for the fourth-order Runge-Kutta method.

Euler函数方程φ(xy)=28(φ(x)+φ(y))的正整数解

利用初等方法研究了不定方程φ(xy)=28(φ(x)+φ(y))的可解性问题,并给出了该方程的全部正整数解,其中φ(n)是Euler函数。

Moving-Particle Semi-Implicit Method for Vortex Patterns and Roll Damping of 2D Ship Sections

A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.

The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate

In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.

Geophysical Study: Estimation of Deposit Depth Using Gravimetric Data and Euler Method (Jalalabad Iron Mine, Kerman Province of IRAN)

Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the drilling boreholes. The purpose of this study is to use field geophysics to calculate the depth of mineral reserve. The study area is located 38 km from Zarand city called Jalalabad iron mine. In this study, gravimetric data were measured and mineral depth was calculated using the Euler method. 1314 readings have been performed in this area. The rocks of the region include volcanic and sedimentary. The source of the mineralization in the area is hydrothermal processes. After gravity measuring in the region, the data were corrected, then various methods such as anomalous map remaining in levels one and two, upward expansion, first and second-degree vertical derivatives, analytical method, and analytical signal were drawn, and finally, the depth of the deposit was estimated by Euler method. As a result, the depth of the mineral deposit was calculated to be between 20 and 30 meters on average.

Consistency and Validity of the Mathematical Models and the Solution Methods for BVPs and IVPs Based on Energy Methods and Principle of Virtual Work for Homogeneous Isotropic and Non-Homogeneous Non-Isotropic Solid Continua

Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous, non-isotropic matter without using (or in the absence of) the mathematical models of the BVPs and the IVPs. These methods are also used for deriving mathematical models for BVPs and IVPs associated with isotropic, homogeneous as well as non-homogeneous, non-isotropic continuous matter. In energy methods when applied to IVPs, one constructs energy functional (I) consisting of kinetic energy, strain energy and the potential energy of loads. The first variation of this energy functional (δI) set to zero is a necessary condition for an extremum of I. In this approach one could use δI = 0 directly in constructing computational processes such as the finite element method or could derive Euler’s equations (differential or partial differential equations) from δI = 0, which is also satisfied by a solution obtained from δI = 0. The Euler’s equations obtained from δI = 0 indeed are the mathematical model associated with the energy functional I. In case of BVPs we follow the same approach except in this case, the energy functional I consists of strain energy and the potential energy of loads. In using the principle of virtual work for BVPs and the IVPs, we can also accomplish the same as described above using energy methods. In this paper we investigate consistency and validity of the mathematical models for isotropic, homogeneous and non-isotropic, non-homogeneous continuous matter for BVPs that are derived using energy functional consisting of strain energy and the potential energy of loads. Similar investigation is also presented for IVPs using energy functional consisting of kinetic energy, strain energy and the potential energy of loads. The computational approaches for BVPs and the IVPs designed using energy functional and principle of virtual work, their consistency and validity are also investigated. Classical continuum mechanics (CCM) principles i.e. conservation and balance laws of CCM with consistent constitutive theories and the elements of calculus of variations are employed in the investigations presented in this paper.

THE FINITE DIFFERENCE STREAMLINE DIFFUSION METHODS FOR TIME-DEPENDENT CONVECTION-DIFFUSION EQUATIONS

In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.

Euler-Bernoulli海洋立管涡致强迫振动响应研究

针对海洋立管(Pipe-in-pipe,PIP)系统在海水作用下发生的振动问题,开展了对PIP系统在涡致强迫振动下的动力学响应研究,分析了在涡致强迫振动下海洋立管外管直径、轴向拉力、外激力频率对海洋立管位移响应的影响规律。基于Euler-Bernoulli双梁模型,采用Lamb-Oseen涡模型,建立了动力学模型,利用格林函数法求得该强迫振动的稳态响应。结果表明,随着管道直径增加,外激力增加,产生最大力幅值的位置离管道越远;轴向拉力对外部管道的影响较大,对内部管道的影响较小;无因次频率取0.4时,外部管道位移超出允许变形极限,内外管壁发生周期碰撞,易对海洋立管造成损伤。

A Preconditioned Gridless Method for Solving Euler Equations at Low Mach Numbers

A preconditioned gridless method is developed for solving the Euler equations at low Mach numbers.The preconditioned system in a conservation form is obtained by multiplying apreconditioning matrix of the type of Weiss and Smith to the time derivative of the Euler equations,which are discretized using agridless technique wherein the physical domain is distributed by clouds of points.The implementation of the preconditioned gridless method is mainly based on the frame of the traditional gridless method without preconditioning,which may fail to converge for low Mach number simulations.Therefore,the modifications corresponding to the affected terms of preconditioning are mainly addressed.The numerical results show that the preconditioned gridless method still functions for compressible transonic flow simulations and additionally,for nearly incompressible flow simulations at low Mach numbers as well.The paper ends with the nearly incompressible flow over a multi-element airfoil,which demonstrates the ability of the method presented for treating flows over complicated geometries.