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.

A Modified Lagrange Method for Solving Convex Quadratic Optimization Problems

In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.

Dispersive propagation of optical solitions and solitary wave solutions of Kundu-Eckhaus dynamical equation via modified mathematical method

In this research work,we constructed the optical soliton solutions of nonlinear complex Kundu-Eckhaus(KE)equation with the help of modified mathematical method.We obtained the solutions in the form of dark solitons,bright solitons and combined dark-bright solitons,travelling wave and periodic wave solutions with general coefficients.In our knowledge earlier reported results of the KE equation with specific coefficients.These obtained solutions are more useful in the development of optical fibers,dynamics of solitons,dynamics of adiabatic parameters,dynamics of fluid,problems of biomedical,industrial phenomena and many other branches.All calculations show that this technique is more powerful,effective,straightforward,and fruitfulness to study analytically other higher-order nonlinear complex PDEs involves in mathematical physics,quantum physics,Geo physics,fluid mechanics,hydrodynamics,mathematical biology,field of engineering and many other physical sciences.

Solving Different Types of Differential Equations Using Modified and New Modified Adomian Decomposition Methods

The Modified Adomian Decomposition Method (MADM) is presented. A number of problems are solved to show the efficiency of the method. Further, a new solution scheme for solving boundary value problems with Neumann conditions is proposed. The scheme is based on the modified Adomian decomposition method and the inverse linear operator theorem. Several differential equations with Neumann boundary conditions are solved to demonstrate the high accuracy and efficiency of the proposed scheme.

A modified single edge V-notched beam method for evaluating surface fracture toughness of thermal barrier coatings

The surface fracture toughness is an important mechanical parameter for studying the failure behavior of air plasma sprayed(APS)thermal barrier coatings(TBCs).As APS TBCs are typical multilayer porous ceramic materials,the direct applications of the traditional single edge notched beam(SENB)method that ignores those typical structural characters may cause errors.To measure the surface fracture toughness more accurately,the effects of multilayer and porous characters on the fracture toughness of APS TBCs should be considered.In this paper,a modified single edge V-notched beam(MSEVNB)method with typical structural characters is developed.According to the finite element analysis(FEA),the geometry factor of the multilayer structure is recalculated.Owing to the narrower V-notches,a more accurate critical fracture stress is obtained.Based on the Griffith energy balance,the reduction of the crack surface caused by micro-defects is corrected.The MSEVNB method can measure the surface fracture toughness more accurately than the SENB method.

Application of the Modified Adomian Decomposition Method on a Mathematical Model of COVID-19

In this study, we constructed and analysed a mathematical model of COVID-19 in order to comprehend the transmission dynamics of the disease. The reproduction number (RC) was calculated via the next generation matrix method. We also used the Lyaponuv method to show the global stability of both the disease free and endemic equilibrium points. The results showed that the disease-free equilibrium point is globally asymptotically stable if RC RC > 1. We further used the Adomian decomposition method and the modified Adomian decomposition method to obtain the solutions of the model. Numerical analysis of the model was done using Sagemath 9.0 software.

Exact Traveling Wave Solutions to Phi-4 Equation and Joseph-Egri (TRLW) Equation and Calogro-Degasperis (CD) Equation by Modified (G'/G2)-Expansion Method

In this study, we will introduce the modified (G'/G2)-expansion method to explore some of the exact traveling wave solutions of some nonlinear partial differential equations namely, Phi-4 equation, Joseph-Egri (TRLW) equation, and Calogro-Degasperis (CD) equation. As a result, we have obtained solutions for the equations expressed in terms of trigonometric, hyperbolic and rational functions. Moreover, some selected solutions are plotted using some specific values for the parameters.

Propagation of traveling wave solutions to the Vakhnenko-Parkes dynamical equation via modified mathematical methods

In this paper, we investigate some new traveling wave solutions to Vakhnenko-Parkes equation via three modified mathematical methods. The derived solutions have been obtained including periodic and solitons solutions in the form of trigonometric, hyperbolic, and rational function solutions. The graphical representations of some solutions by assigning particular values to the parameters under prescribed conditions in each solutions and comparing of solutions with those gained by other authors indicate that these employed techniques are more effective, efficient and applicable mathematical tools for solving nonlinear problems in applied science.

Density of Ni-Al Alloys in Liquid and Solid-Liquid Coexistence State Measured by a Modified Pycnometric Method

The density of Ni-AI alloys in both liquid state and solid-liquid coexistence state was measured with a modified pycnometric method. It was found that the density of Ni-AI alloys decreases with increasing temperature and Al concentration in the alloys. The molar volume of liquid Ni-AI binary alloys increases with the increase of temperature and Al concentration. The partial molar volume of Al in Ni-AI binary alloy was calculated approximately. The molar volume of liquid Ni-AI alloy determined in the present work shows a negative deviation from the ideal linear molar volume.

Grouping response method for equivalent static wind loads based on a modified LRC method

Wind loading is one of the most important loads for controlling the design of large-span roof structures. Equivalent static wind loads, which can generally aim at determining a specific response, are widely used by structural designers. A method for equivalent static wind loads applicable to multi-responses is proposed in this paper. A modified load- response-correlation （LRC） method corresponding to a particular peak response is presented, and the similarity algorithm implemented for the group response is described. The main idea of the algorithm is that two responses can be put into one group if the value of one response is close to that of the other response, when the structure is subjected to equivalent static wind loads aiming at the other response. Based on the modified LRC, the grouping response method is put forward to construct equivalent static wind loading. This technique can simultaneously reproduce peak responses for some grouped responses. To verify its computational accuracy, the method is applied to an actual large-span roof structure. Calculation results show that when the similarity of responses in the same group is high, equivalent static wind loads with high accuracy and reasonable magnitude of equivalent static wind distribution can be achieved.

Density of Liquid Ni-Mo Alloys Measured by a Modified Sessile Drop Method

The density of liquid binary Ni-Mo alloys with molybdenum concentration from 0 to 20% (mass fraction) wasmeasured by a modified sessile drop method. It has been found that the density of the liquid Ni-Mo alloys decreaseswith increasing temperature, but increases with the increase of molybdenum concentration in the alloys. The molarvolume of liquid Ni-Mo binary alloys increases with the increase of temperature and molybdenum concentration. Thepartial molar volume of molybdenum in Ni-Mo binary alloy has been approximately calculated as [13.18-2.65×10^(-3)T+(-47.94+3.10×10^(-2)T)×10^(-2)X_(Mo)]×10^(-6)m^3·mol^(-1). The molar volume of Ni-Mo alloy determined inthe present work shows a negative deviation from the ideal linear mixing molar volume.

A modified domain reduction method for numerical simulation of wave propagation in localized regions

A modified domain reduction method(MDRM) that introduces damping terms to the original DRM is presented in this paper. To verify the proposed MDRM and compare the computational accuracy of these two methods, a numerical test is designed. The numerical results of the MDRM and DRM are compared using an extended meshed model. The results show that the MDRM significantly improved the computational accuracy of the DRM. Then, the MDRM is compared with two existing conventional methods, namely Liao's transmitting boundary and viscous-spring boundary with Liu's method. The MDRM shows its great advancement in computational accuracy, stability and range of applications. This paper also discusses the influence of boundary location on computational accuracy. It can be concluded that smaller models tend to have larger errors. By introducing two dimensionless parameters, φ_1 and φ_2, the rational distance between the observation point and the MDRM boundary is suggested. When φ_1 >2 or φ_2>13, the relative PGA error can be limited to 5%. In practice, the appropriate model size can be chosen based on these two parameters to achieve desired computational accuracy.

Lateral Bearing Capacity of Modified Suction Caissons Determined by Using the Limit Equilibrium Method

The modified suction caisson(MSC) adds a short-skirted structure around the regular suction caissons to increase the lateral bearing capacity and limit the deflection. The MSC is suitable for acting as the offshore wind turbine foundation subjected to larger lateral loads compared with the imposed vertical loads. Determination of the lateral bearing capacity is a key issue for the MSC design. The formula estimating the lateral bearing capacity of the MSC was proposed in terms of the limit equilibrium method and was verified by the test results. Parametric studies on the lateral bearing capacity were also carried out. It was found that the lateral bearing capacity of the MSC increases with the increasing length and radius of the external skirt, and the lateral bearing capacity increases linearly with the increasing coefficient of subgrade reaction. The maximum lateral bearing capacity of the MSC is attained when the ratio of the radii of the internal compartment to the external skirt equals 0.82 and the ratio of the lengths of the external skirt to the internal compartment equals 0.48, provided that the steel usage of the MSC is kept constant.

Solving mKdV-sinh-Gordon equation by a modified variable separated ordinary differential equation method

By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact solutions including some new formal solutions are successfully picked up for the mKdV-sinh-Gordon equation by this approach.

Calculation of circulation in the South China Sea during summer of 2000 by the modified inverse method

On the basis of hydrographic data obtained in August 2000 cruise, the circulation in the South China Sea (SCS) is computed by the modified inverse method in combination with SSH data from TOPEX/ERS-2 analysis. For study of the dynamical mechanism, which causes the pattern of summer circulation in the SCS, the diagnostic model (Yuan et al. 1982. Acta Oceanologica Sinica,4(1):1-11; Yuan and Su. 1992. Numerical Computation of Physical Oceanography.474-542) is used to simulate numerically the summer circulation in the SCS. The following results have been obtained. (1) The central and southwestern SCSs are dominated mainly by anticy-clonic circulation systems. They are mainly as follows. 1) There is strong anticyclonic eddy southeast of Vietnam (W1). Its horizontal scale is about 300 km, and it extends vertically from the surface to the about 1 000 m level. 2) There are a warm eddy W2 southeast of Zhongsha Islands and the anticyclonic circulation system W3 west off the Luzon Island. 3) There is a stronger cyclonic eddy C1 between the anticyclonic eddies W1 and W2.4) A strong northward coastal jet is present near the coast of Vietnam, and separates from the coast of Vietnam at about 12° N to the northeast.(2)The northern SCS is dominated mainly by a cyclonic circulation system. There is a cyclonic circulation system near and north of Section N2. (3) The southeastern SCS is dominated mainly by the cyclonic circulation system. (4) Comparing the results of circulation in the SCS during the summer of 2000 with those during the summer of 1998, it is found that they agree qualitatively, but there is the some difference between them in quantity.This shows that the circulation in the SCS has obviously seasonal feature. (5) The dynamical mechanism which products the basic pattern of summer circulation is because the following two reasons: 1) the joint effect of the baroclinity and relief (JEBAR) is essential dynamical cause; and 2) it is next important dynamical cause that the interaction between the wind stress and bottom topography under the southerly monsoon. (6) Comparing the hydrographic structure and distribution of stream functions with the SSH data from TOPEX/ERS-2 analysis in the SCS during August of 2000, they agree qualitatively.

A modified stiffness spreading method for layout optimization of truss structures

The stiffness spreading method (SSM) was initially proposed for layout optimization of truss structures,in which an artificial elastic material of low modulus is uniformly distributed in the design domain to create connections between discrete members.In this paper,a modified stiffness spreading method is proposed by replacing the artificial elastic material with auxiliary bars to connect real members of the truss structure.Since the background continuum mesh for the elastic material is no longer required,the computational cost is significantly reduced.Like SSM,the new method is advantageous in that an initial design may consist of disconnected bars allocated in the design domain,and mathematical programming methods can be applied for the efficient solution of the formulated optimization problem.A number of solution strategies are also developed to achieve more practical designs with lower computational cost.Numerical examples of both 2-D and 3-D truss structures are presented to demonstrate the feasibility,robustness and effectiveness of the proposed method.

Seismic passive earth resistance using modified pseudo-dynamic method

In earthquake prone areas, understanding of the seismic passive earth resistance is very important for the design of different geotechnical earth retaining structures. In this study, the limit equilibrium method is used for estimation of critical seismic passive earth resistance for an inclined wall supporting horizontal cohesionless backfill. A composite failure surface is considered in the present analysis. Seismic forces are computed assuming the backfill soil as a viscoelastic material overlying a rigid stratum and the rigid stratum is subjected to a harmonic shaking. The present method satisfies the boundary conditions. The amplification of acceleration depends on the properties of the backfill soil and on the characteristics of the input motion. The acceleration distribution along the depth of the backfill is found to be nonlinear in nature. The present study shows that the horizontal and vertical acceleration distribution in the backfill soil is not always in-phase for the critical value of the seismic passive earth pressure coefficient. The effect of different parameters on the seismic passive earth pressure is studied in detail. A comparison of the present method with other theories is also presented, which shows the merits of the present study.

High-pressure phonon dispersion of copper by using the modified analytic embedded atom method

By using the Born-von Kfirmfin theory of lattice dynamics and the modified analytic embedded atom method, we reproduce the experimental results of the phonon dispersion in fcc metal Cu at zero pressure along three high symmetry directions and four oft-symmetry directions, and then simulate the phonon dispersion curves of Cu at high pressures of 50, 100, and 150 GPa. The results show that the shapes of dispersion curves at high pressures are very similar to that at zero pressure. All the vibration frequencies of Cu in all vibration branches at high pressures are larger than the results at zero pressure, and increase correspondingly as pressure reaches 50, 100, and 150 GPa sequentially. Moreover, on the basis of phonon dispersion, we calculate the values of specific heat of Cu at different pressures. The prediction of thermodynamic quantities lays a significant foundation for guiding and judging experiments of thermodynamic properties of solids under high pressures.

A modified Tikhonov regularization method for a Cauchy problem of a time fractional diffusion equation

In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.

Applications of F-expansion method to the coupled KdV system

An extended F-expansion method for finding periodic wave solutions of nonlinear evolution equations in mathematical physics is presented, which can be thought of as a concentration of extended Jacobi elliptic function expansion method proposed more recently. By using the homogeneous balance principle and the extended F-expansion, more periodic wave solutions expressed by Jacobi elliptic functions for the coupled KdV equations are derived. In the limit cases, the solitary wave solutions and the other type of travelling wave solutions for the system are also obtained.