Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

Crank-Nicolson Quasi-Compact Scheme for the Nonlinear Two-Sided Spatial Fractional Advection-Diffusion Equations

The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference operators and combining the compact technique, in the time direction is discretized by the Crank-Nicolson method. Through the energy method, the stability and convergence of the numerical scheme in the sense of L2-norm are proved, and the convergence order is . Some examples are given to show that our numerical scheme is effective.

Alternating Group Explicit Iterative Methods for One-Dimensional Advection-Diffusion Equation

The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.

Discontinuous Legendre Wavelet Galerkin Method for One-Dimensional Advection-Diffusion Equation

This paper presents discontinuous Legendre wavelet Galerkin (DLWG) approach for solving one-dimensional advection-diffusion equation (ADE). Variational formulation of this type equation and corresponding numerical fluxes are devised by utilizing the advantages of both the Legendre wavelet bases and discontinuous Galerkin (DG) method. The distinctive features of the proposed method are its simple applicability for a variety of boundary conditions and able to effectively approximate the solution of PDEs with less storage space and execution. The results of a numerical experiment are provided to verify the efficiency of the designed new technique.

RETROSPECTIVE TIME INTEGRAL SCHEME AND ITS APPLICATIONS TO THE ADVECTION EQUATION

To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.

STABILITY OF VISCOUS SHOCK WAVES FOR THE ONE-DIMENSIONAL COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH DENSITY-DEPENDENT VISCOSITY

We study the large-time behavior toward viscous shock waves to the Cauchy problem of the one-dimensional compressible isentropic Navier-Stokes equations with density- dependent viscosity. The nonlinear stability of the viscous shock waves is shown for certain class of large initial perturbation with integral zero which can allow the initial density to have large oscillation. Our analysis relies upon the technique developed by Kanel~ and the continuation argument.

One-dimensional dynamic equations of a piezoelectric semiconductor beam with a rectangular cross section and their application in static and dynamic characteristic analysis

Within the framework of continuum mechanics, the double power series ex- pansion technique is proposed, and a series of reduced one-dimensional （1D） equations for a piezoelectric semiconductor beam are obtained. These derived equations are universal, in which extension, flexure, and shear deformations are all included, and can be degen- erated to a number of special cases, e.g., extensional motion, coupled extensional and flexural motion with shear deformations, and elementary flexural motion without shear deformations. As a typical application, the extensional motion of a ZnO beam is analyzed sequentially. It is revealed that semi-conduction has a great effect on the performance of the piezoelectric semiconductor beam, including static deformations and dynamic be- haviors. A larger initial carrier density will evidently lead to a lower resonant frequency and a smaller displacement response, which is a little similar to the dissipative effect. Both the derived approximate equations and the corresponding qualitative analysis are general and widely applicable, which can clearly interpret the inner physical mechanism of the semiconductor in the piezoelectrics and provide theoretical guidance for further experimental design.

NONLINEAR STABILITY OF VISCOUS SHOCK WAVES FOR ONE-DIMENSIONAL NONISENTROPIC COMPRESSIBLE NAVIER–STOKES EQUATIONS WITH A CLASS OF LARGE INITIAL PERTURBATION

We study the nonlinear stability of viscous shock waves for the Cauchy problem of one-dimensional nonisentropic compressible Navier–Stokes equations for a viscous and heat conducting ideal polytropic gas. The viscous shock waves are shown to be time asymptotically stable under large initial perturbation with no restriction on the range of the adiabatic exponent provided that the strengths of the viscous shock waves are assumed to be sufficiently small.The proofs are based on the nonlinear energy estimates and the crucial step is to obtain the positive lower and upper bounds of the density and the temperature which are uniformly in time and space.

Enriched reproducing kernel particle method for fractional advection–diffusion equation

The reproducing kernel particle method （RKPM） has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modeled by a fractional advectiondiffusion equation （FADE）, which exhibits a boundary layer with low regularity. We formulate this method on a moving least-square approach. Via the enrichment of fractional-order power functions to the traditional integer-order basis for RKPM, leading terms of the solution to the FADE can be exactly reproduced, which guarantees a good approximation to the boundary layer. Numerical tests are performed to verify the proposed approach.

Chebyshev Pseudo-Spectral Method for Solving Fractional Advection-Dispersion Equation

Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.

Comparison between Non-Gaussian Puff Model and a Model Based on a Time-Dependent Solution of Advection-Diffusion Equation

A comparison between a non-Gaussian puff model and an advanced time-dependent model to simulate the pollutant dispersion in the Planetary Boundary Layer is presented. The puff model is based on a general technique for solving the K-equation, using the truncated Gram-Charlier expansion (type A) of the concentration field and finite set equations for the corresponding moments. The other model (named ADMM: Analytical Dispersion Multilayers Model) is an semi- analytical solution to the time-dependent two-dimensional advection-diffusion equation based on a discretization of the PBL in N sub-layers;in each sub-layers the advection-diffusion equation is solved by the Laplace transform technique, considering an average value for eddy diffusivity and the wind speed. A preliminary performance evaluation is shown in the case of continuous emission from an elevated source in a variable boundary layer. Both models were able to correctly reproduce the concentration field measured and so to be used as operative air pollution models.

One-Dimensional Explicit Tolesa Numerical Scheme for Solving First Order Hyperbolic Equations and Its Application to Macroscopic Traffic Flow Model

In this paper, a new numerical scheme for solving first-order hyperbolic partial differential equations is proposed and is implemented in the simulation study of macroscopic traffic flow model with constant velocity and linear velocity-density relationship. Macroscopic traffic flow model is first developed by Lighthill Whitham and Richards (LWR) and used to study traffic flow by collective variables such as flow rate, velocity and density. The LWR model is treated as an initial value problem and its numerical simulations are presented using numerical schemes. A variety of numerical schemes are available in literature to solve first order hyperbolic equations. Of these the well-known ones include one-dimensional explicit: Upwind, Downwind, FTCS, and Lax-Friedrichs schemes. Having been studied carefully the space and time mesh sizes, and the patterns of all these schemes, a new scheme has been developed and named as one-dimensional explicit Tolesa numerical scheme. Tolesa numerical scheme is one of the conditionally stable and highest rates of convergence schemes. All the said numerical schemes are applied to solve advection equation pertaining traffic flows. Also the one-dimensional explicit Tolesa numerical scheme is another alternative numerical scheme to solve advection equation and apply to traffic flows model like other well-known one-dimensional explicit schemes. The effect of density of cars on the overall interactions of the vehicles along a given length of the highway and time are investigated. Graphical representations of density profile, velocity profile, flux profile, and in general the fundamental diagrams of vehicles on the highway with different time levels are illustrated. These concepts and results have been arranged systematically in this paper.

Air Pollution Steady-State Advection-Diffusion Equation: The General Three-Dimensional Solution

Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Many models simulating air pollution dispersion are based upon the solution (numerical or analytical) of the advection-diffusion equation as- suming turbulence parameterization for realistic physical scenarios. We present the general steady three-dimensional solution of the advection-diffusion equation considering a vertically inhomogeneous atmospheric boundary layer for arbitrary vertical profiles of wind and eddy-diffusion coefficients. Numerical results and comparison with experimental data are shown.

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation

The interpolation method in a semi-Lagrangian scheme is decisive to its performance. Given the number of grid points one is considering to use for the interpolation, it does not necessarily follow that maximum formal accuracy should give the best results. For the advection equation, the driving force of this method is the method of the characteristics, which accounts for the flow of information in the model equation. This leads naturally to an interpolation problem since the foot point is not in general located on a grid point. We use another interpolation scheme that will allow achieving the high order for the box initial condition.

Advection Dispersion Equation and BMO Space

In this paper, we provide a new way of characterizing the upper and lower bound for the concentration and the gradient of concentration in advection dispersion equation under the condition that source term, concentration and stirring term belong to BMO space.

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.

Potential Symmetries, One-Dimensional Optimal System and Invariant Solutions of the Coupled Burgers’ Equations

In this paper, we discuss one-dimensional optimal system and the invariant solutions of Coupled Burgers’ equations. By using Wu-differential characteristic set algorithm with the aid of Mathematica software, the classical symmetries of the Coupled Burgers’ equations are calculated, and the one-dimensional optimal system of Lie algebra is constructed. And we obtain the invariant solution of the Coupled Burgers’ equations corresponding to one element in one dimensional optimal system by using the invariant method. The results generalize the exact solutions of the Coupled Burgers’ equations.

LOCAL ONE-DIMENSIONAL ASE-I SCHEME FOR 2D DIFFUSION EQUATION

A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some numerical experiments show the method is not only simple but also more accurate.

Dynamical Soliton Wave Structures of One-Dimensional Lie Subalgebras via Group-Invariant Solutions of a Higher-Dimensional Soliton Equation with Various Applications in Ocean Physics and Mechatronics Engineering

Having realized various significant roles that higher-dimensional nonlinear partial differ-ential equations(NLPDEs)play in engineering,we analytically investigate in this paper,a higher-dimensional soliton equation,with applications particularly in ocean physics and mechatronics(electrical electronics and mechanical)engineering.Infinitesimal generators of Lie point symmetries of the equation are computed using Lie group analysis of differen-tial equations.In addition,we construct commutation as well as Lie adjoint representation tables for the nine-dimensional Lie algebra achieved.Further,a one-dimensional optimal system of Lie subalgebras is also presented for the soliton equation.This consequently enables us to generate abundant group-invariant solutions through the reduction of the understudy equation into various ordinary differential equations(ODEs).On solving the achieved nonlinear differential equations,we secure various solitonic solutions.In conse-quence,these solutions containing diverse mathematical functions furnish copious shapes of dynamical wave structures,ranging from periodic,kink and kink-shaped nanopteron,soliton(bright and dark)to breather waves with extensive wave collisions depicted.We physically interpreted the resulting soliton solutions by imploring graphical depictions in three dimensions,two dimensions and density plots.Moreover,the gained group-invariant solutions involved several arbitrary functions,thus exhibiting rich physical structures.We also implore the power series technique to solve part of the complicated differential equa-tions and give valid comments on their results.Later,we outline some applications of our results in ocean physics and mechatronics engineering.