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.

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.

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.

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.

A Numerical Study of One-Dimensional Hyperbolic Telegraph Equation

In this paper, an approximate solution for the one-dimensional hyperbolic telegraph equation by using the q-homotopy analysis method （q-HAM） is proposed.The results shows that the convergence of the q- homotopy analysis method is more accurate than the convergence of the homotopy analysis method （HAM）.

Nonhomogeneous(H,Q)-Process:The Backward and Forward Equations

As for the backward and forward equation of nonhomogeneous(H, Q) -processes,we proof them in a new way. On the base of that, this paper gives the direct computational formalfor one dimensional distribution of the nonhomogeneous(H, Q) -process.

TIME-HARMONIC DYNAMIC GREEN'S FUNCTIONS FOR ONE-DIMENSIONAL HEXAGONAL QUASICRYSTALS

Quasicrystals have additional phason degrees of freedom not found in conventional crystals. In this paper, we present an exact solution for time-harmonic dynamic Green＇s function of one-dimensional hexagonal quasicrystals with the Laue classes 6/mh and 6/mhmm. Through the introduction of two new functions φ and ψ, the original problem is reduced to the determination of Green＇s functions for two independent Helmholtz equations. The explicit expressions of displacement and stress fields are presented and their asymptotic behaviors are discussed. The static Green＇s function can be obtained by letting the circular frequency approach zero.

Solvability of Non-homogeneous Mixed Type Multi-point BVPs for Second Order Differential Equations with p-Laplacian

A class of multi-point boundary value problems are studied.Easily verified suffcient conditions to guarantee the existence of at least one solutions of above mentioned BVPs are established.The examples are presented to illustrate the main results.

Estimation of vertical diffusion coefficient based on a onedimensional temperature diffusion equation with an inverse method

Diapycnal mixing is important in oceanic circulation. An inverse method in which a semi-explicit scheme is applied to discretize the one-dimensional temperature diffusion equation is established to estimate the vertical temperature diffusion coefficient based on the observed temperature profiles. The sensitivity of the inverse model in the idealized and actual conditions is tested in detail. It can be found that this inverse model has high feasibility under multiple situations ensuring the stability of the inverse model, and can be considered as an efficient way to estimate the temperature diffusion coefficient in the weak current regions of the ocean. Here, the hydrographic profiles from Argo floats are used to estimate the temporal and spatial distribution of the vertical mixing in the north central Pacific based on this inverse method. It is further found that the vertical mixing in the upper ocean displays a distinct seasonal variation with the amplitude decreasing with depth, and the vertical mixing over rough topography is stronger than that over smooth topography It is suggested that the high-resolution profiles from Argo floats and a more reasonable design of the inverse scheme will serve to understand mixing processes.

One-dimensional Mathematical Model of Numerical Simulation of Urban Heat Island

With the development of society and the acceleration of urbanization, urban heat island phenomenon is becoming increasingly prominent. In this paper, enrichment capacity of plant light energy is supplemented based on OKe model according to law of energy conservation. By contrasting each component of energy balance equation between the cities and rural areas, the causes for urban heat island are explained.

Hyperbolic Monge-Ampère Equation

In this paper, based on the Lie symmetry method, the symmetry group of a hyperbolic Monge-Ampère equation is obtained first, then the one-dimensional optimal system of the obtained symmetries is given, and finally the group-invariant solutions are investigated.

One-dimensional monotone nonautonomous dynamical systems

This work is devoted to the study of the dynamics of one-dimensional monotone nonautonomous(cocycle) dynamical systems. A description of the structures of their invariant sets, omega limit sets,Bohr/Levitan almost periodic and almost automorphic motions, global attractors, and pinched and minimalsets is given. An application of our general results is given to scalar differential and difference equations.

A New Directional Algebraic Fast Multipole Method Based Iterative Solver for the Lippmann-Schwinger Equation Accelerated with HODLR Preconditioner

We present a fast iterative solver for scattering problems in 2D,where a penetrable object with compact support is considered.By representing the scattered field as a volume potential in terms of the Green’s function,we arrive at the Lippmann-Schwinger equation in integral form,which is then discretized using an appropriate quadrature technique.The discretized linear system is then solved using an iterative solver accelerated by Directional Algebraic Fast Multipole Method(DAFMM).The DAFMM presented here relies on the directional admissibility condition of the 2D Helmholtz kernel[1],and the construction of low-rank factorizations of the appropriate low-rank matrix sub-blocks is based on our new Nested Cross Approximation(NCA)[2].The advantage of the NCA described in[2]is that the search space of so-called far-field pivots is smaller than that of the existing NCAs[3,4].Another significant contribution of this work is the use of HODLR based direct solver[5]as a preconditioner to further accelerate the iterative solver.In one of our numerical experiments,the iterative solver does not converge without a preconditioner.We show that the HODLR preconditioner is capable of solving problems that the iterative solver can not.Another noteworthy contribution of this article is that we perform a comparative study of the HODLR based fast direct solver,DAFMMbased fast iterative solver,and HODLR preconditioned DAFMM based fast iterative solver for the discretized Lippmann-Schwinger problem.To the best of our knowledge,this work is one of the first to provide a systematic study and comparison of these different solvers for various problem sizes and contrast functions.In the spirit of reproducible computational science,the implementation of the algorithms developed in this article is made available at https://github.com/vaishna77/Lippmann_Schwinger_Solver.

Global Solution for Fully Nonlinear Parabolic Equations in One-Dimensional Space

We discuss the existence of global classical solution for the uniformly parabolic equation ■ut=a(x,t,u,ux,uxx)+b（x,t,u,ux）,（x,t）∈（-1,1）×(0,T], u（±1,t）=0,u（x,0）=■（x）, where a is strongly nonlinear with respect to uxxand ■ is not necessarily small.We also deal with nonuniform case.

One-Dimensional Symmetry and Liouville Type Results for the Fourth Order Allen-Cahn Equation in R^N

In this paper, the authors prove an analogue of Gibbons' conjecture for the extended fourth order Allen-Cahn equation in R^N, as well as Liouville type results for some solutions converging to the same value at infinity in a given direction. The authors also prove a priori bounds and further one-dimensional symmetry and rigidity results for semilinear fourth order elliptic equations with more general nonlinearities.

Study on the prediction and inverse prediction of detonation properties based on deep learning

The accurate and efficient prediction of explosive detonation properties has important engineering significance for weapon design.Traditional methods for predicting detonation performance include empirical formulas,equations of state,and quantum chemical calculation methods.In recent years,with the development of computer performance and deep learning methods,researchers have begun to apply deep learning methods to the prediction of explosive detonation performance.The deep learning method has the advantage of simple and rapid prediction of explosive detonation properties.However,some problems remain in the study of detonation properties based on deep learning.For example,there are few studies on the prediction of mixed explosives,on the prediction of the parameters of the equation of state of explosives,and on the application of explosive properties to predict the formulation of explosives.Based on an artificial neural network model and a one-dimensional convolutional neural network model,three improved deep learning models were established in this work with the aim of solving these problems.The training data for these models,called the detonation parameters prediction model,JWL equation of state(EOS)prediction model,and inverse prediction model,was obtained through the KHT thermochemical code.After training,the model was tested for overfitting using the validation-set test.Through the model-accuracy test,the prediction accuracy of the model for real explosive formulations was tested by comparing the predicted value with the reference value.The results show that the model errors were within 10%and 3%for the prediction of detonation pressure and detonation velocity,respectively.The accuracy refers to the prediction of tested explosive formulations which consist of TNT,RDX and HMX.For the prediction of the equation of state for explosives,the correlation coefficient between the prediction and the reference curves was above 0.99.For the prediction of the inverse prediction model,the prediction error of the explosive equation was within 9%.This indicates that the models have utility in engineering.