Higher-order energy consistent time integrators for nonlinear thermoviscoelastodynamics

An advantage of the temporal fe method is that higher-order accurate time integrators can be constructed easily. A further important advantage is the inherent energy consistency if applied to equations of motion. The temporal fe method is therefore used to construct higher-order energy-momentum conserving time integrators for nonlinear elastodynamics (see Ref. [1]). Considering finite motions of a flexible solid body with internal dissipation, an energy consistent time integration is also of great advantage (see the references [2, 3]). In this paper, we show that an energy consistent time integration is also advantageous for dynamics with dissipation arising from conduction of heat as well as from a viscous material. The energy consistency is preserved by using a new enhanced hybrid Galerkin (ehG) method. The obtained numerical schemes satisfy the energy balance exactly, independent of their accuracy and the used time step size. This guarantees numerical stability. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

高阶波动方程的时空估计与低能量散射

本文研究了高阶波动方程的低能量散射理论，基本工具是高阶线性波动方程解的时空估计．与经典的二阶波动方程解的时空估计证明不同，我们采用泛函分析的方法与待定指标技巧，首次给出了高阶线性波动方程的时空估计，藉此与非线性函数在齐次Ｓｏｂｏｌｅｖ空间中的估计，获得了高阶波动方程的低能量散射结论．与此同时，也得到了具临界增长的高阶波动方程的柯西问题在低能量条件下的整体存在唯一性．     

A higher-order moment method of the lattice Boltzmann model for the conservation law equation

In this paper, we proposed a higher-order moment method in the lattice Boltzmann model for the conservation law equation. In contrast to the lattice Bhatnagar–Gross–Krook (BGK) model, the higher-order moment method has a wide flexibility to select equilibrium distribution function. This method is based on so-called a series of partial differential equations obtained by using multi-scale technique and Chapman–Enskog expansion. According to Hirt’s heuristic stability theory, the stability of the scheme can be controlled by modulating some special moments to design the third-order dispersion term and the fourth-order dissipation term. As results, the conservation law equation is recovered with higher-order truncation error. The numerical examples show the higher-order moment method can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the conservation law equation.     

Uniqueness and symmetry of ground states for higher-order equations

We establish uniqueness and radial symmetry of ground states for higher-order nonlinear Schrödinger and Hartree equations whose higher-order differentials have small coefficients. As an application, we obtain error estimates for higher-order approximations to the pseudo-relativistic ground state. Our proof adapts the strategy of Lenzmann (Anal PDE 2:1–27, 2009) using local uniqueness near the limit of ground states in a variational problem. However, in order to bypass difficulties from lack of symmetrization tools for higher-order differential operators, we employ the contraction mapping argument in our earlier work (Choi et al. 2017. arXiv:1705.09068) to construct radially symmetric real-valued solutions, as well as improving local uniqueness near the limit.     

Higher-order total variation bounds for expectations of periodic functions and simple integer recourse approximations

We derive bounds on the expectation of a class of periodic functions using the total variations of higher-order derivatives of the underlying probability density function. These bounds are a strict improvement over those of Romeijnders et al. (Math Program 157:3–46, 2016b), and we use them to derive error bounds for convex approximations of simple integer recourse models. In fact, we obtain a hierarchy of error bounds that become tighter if the total variations of additional higher-order derivatives are taken into account. Moreover, each error bound decreases if these total variations become smaller. The improved bounds may be used to derive tighter error bounds for convex approximations of more general recourse models involving integer decision variables.     

Iterative operator-splitting methods with higher-order time integration methods and applications for parabolic partial differential equations

In this paper we design higher-order time integrators for systems of stiff ordinary differential equations. We combine implicit Runge–Kutta and BDF methods with iterative operator-splitting methods to obtain higher-order methods. The idea of decoupling each complicated operator in simpler operators with an adapted time scale allows to solve the problems more efficiently. We compare our new methods with the higher-order fractional-stepping Runge–Kutta methods, developed for stiff ordinary differential equations. The benefit is the individual handling of each operator with adapted standard higher-order time integrators. The methods are applied to equations for convection–diffusion reactions and we obtain higher-order results. Finally we discuss the applications of the iterative operator-splitting methods to multi-dimensional and multi-physical problems.     

GLOBAL EXISTENCE AND OPTIMAL CONVERGENCE RATES OF SOLUTIONS FOR THREE-DIMENSIONAL ELECTROMAGNETIC FLUID SYSTEM

In this article, we study the electromagnetic fluid system in three-dimensional whole space R~3. Under assumption of small initial data, we establish the unique global solution by energy method. Moreover, we obtain the time decay rates of the higher-order spatial derivatives of the solution by combining the L~p-L~q estimates for the linearized equations and an elaborate energy method when the L~1-norm of the perturbation is bounded.     

An exact local error representation of exponential operator splitting methods for evolutionary problems and applications to linear Schrödinger equations in the semi-classical regime

In this paper, we are concerned with the derivation of a local error representation for exponential operator splitting methods when applied to evolutionary problems that involve critical parameters. Employing an abstract formulation of differential equations on function spaces, our framework includes Schrödinger equations in the semi-classical regime as well as parabolic initial-boundary value problems with high spatial gradients. We illustrate the general mechanism on the basis of the first-order Lie splitting and the second-order Strang splitting method. Further, we specify the local error representation for a fourth-order splitting scheme by Yoshida. From the given error estimate it is concluded that higher-order exponential operator splitting methods are favourable for the time-integration of linear Schrödinger equations in the semi-classical regime with critical parameter 0<ε?1, provided that the time stepsize h is sufficiently smaller than \(\sqrt[p]{\varepsilon}\), where p denotes the order of the splitting method.     

Smooth Frame-Path Termination for Higher Order Sigma-Delta Quantization

We study the performance of finite frames for the encoding of vectors by applying standard higher-order sigma-delta quantization to the frame coefficients. Our results are valid for any quantizer with accuracy ε > 0 operating in the no-overload regime. The frames under consideration are obtained from regular sampling of a path in a Hilbert space. In order to achieve error bounds that are comparable to results on higher-order sigma-delta for the quantization of oversampled bandlimited functions, we construct frame paths that terminate smoothly in the zero vector, that is, with an appropriate number of vanishing derivatives at the endpoint.     

New construction of higher-order local continuous platforms for Error Correction Methods

Error correction method (ECM)~\cite{kim2011a,kim2011b} which has been recently developed, is based on the construction of a local approximation to the solution on each time step, and has the excellent convergence order $O(h^{2p+2})$, provided the local approximation has a local residual error $O(h^p)$. In this paper, we construct a higher-order continuous local platform to develop higher-order semi-explicit one-step ECM for solving initial value time dependent differential equations. It is shown that special choices of parameters for the local platform can lead to the improvement of the well-known explicit fourth and fifth order Runge-Kutta methods. Numerical experiments demonstrate the theoretical results     

Spurious modes in finite-element discretizations of the wave equation may not be all that bad

If higher-order finite elements are used to discretize the wave equation, spurious modes may occur. These modes are classified as unphysical and supposedly make elements of high order useless for accurate computations. This is in conflict with numerical experiments which appear to provide good results. Here Fourier analysis is used to investigate the behaviour of the numerical error for a number of higher-order one-dimensional finite elements. It is shown that the spurious modes have a contribution to the numerical error that behaves in a reasonable manner, and that higher-order elements can be more accurate than lower-order elements. Lumped elements with Gauss–Lobatto nodes appear to be the best choice.     

Geometric Integrators for Higher-Order Variational Systems and Their Application to Optimal Control

Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for higher-order Lagrangian systems. Given a regular higher-order Lagrangian \(L:T^{(k)}Q\rightarrow {\mathbb {R}}\) with \(k\ge 1\), the resulting discrete equations define a generally implicit numerical integrator algorithm on \(T^{(k-1)}Q\times T^{(k-1)}Q\) that approximates the flow of the higher-order Euler–Lagrange equations for L. The algorithm equations are called higher-order discrete Euler–Lagrange equations and constitute a variational integrator for higher-order mechanical systems. The general idea for those variational integrators is to directly discretize Hamilton’s principle rather than the equations of motion in a way that preserves the invariants of the original system, notably the symplectic form and, via a discrete version of Noether’s theorem, the momentum map. We construct an exact discrete Lagrangian \(L_d^e\) using the locally unique solution of the higher-order Euler–Lagrange equations for L with boundary conditions. By taking the discrete Lagrangian as an approximation of \(L_d^e\), we obtain variational integrators for higher-order mechanical systems. We apply our techniques to optimal control problems since, given a cost function, the optimal control problem is understood as a second-order variational problem.     

Higher-order finite volume methods for elliptic boundary value problems

This paper studies higher-order finite volume methods for solving elliptic boundary value problems. We develop a general framework for construction and analysis of higher-order finite volume methods. Specifically, we establish the boundedness and uniform ellipticity of the bilinear forms for the methods, and show that they lead to an optimal error estimate of the methods. We prove that the uniform local-ellipticity of the family of the bilinear forms ensures its uniform ellipticity. We then establish necessary and sufficient conditions for the uniform local-ellipticity in terms of geometric requirements on the meshes of the domain of the differential equation, and provide a general way to investigate the mesh geometric requirements for arbitrary higher-order schemes. Several useful examples of higher-order finite volume methods are presented to illustrate the mesh geometric requirements.     

Global existence for a nonlinear system in thermoviscoelasticity with nonconvex energy

A three-dimensional thermoviscoelastic system derived from the balance laws of momentum and energy is considered. To describe structural phase transitions in solids, the stored energy function is not assumed to be convex as a function of the deformation gradient. A novel feature for multi-dimensional, nonconvex, and nonisothermal problems is that no regularizing higher-order terms are introduced. The mechanical dissipation is not linearized. We prove existence global in time. The approach is based on a fixed-point argument using an implicit time discretization and the theory of renormalized solutions for parabolic equations with L¹ data.     

Adaptive FEM with goal-oriented error estimation and an approximation of the dual problem for inelastic problems

An approximation of the dual problem for goal-oriented error estimation is derived, which preserves the numerical structure of the primal problem, usually found in standard finite element implementations for solid mechanics. For the error estimation we apply a higher-order recovery procedure. An application to elasto-plasticity is given. Numerical examples demonstrate the effectiveness of the procedure. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic behavior of solutions to the compressible nematic liquid crystal system in ℝ3

In this paper, we study a nematic liquid crystals system in three-dimensional whole space ?³ and obtain the time decay rates of the higher-order spatial derivatives of the solution by the method of spectral analysis and energy estimates if the initial data belongs to L¹?³ additionally.     

Compound Poisson Approximation to Convolutions of Compound Negative Binomial Variables

In this paper, the problem of compound Poisson approximation to the convolution of compound negative binomial distributions, under total variation distance, is considered. First, we obtain an error bound using the method of exponents and it is compared with existing ones. It is known that Kerstan’s method is more powerful in compound approximation problems. We employ Kerstan’s method to obtain better estimates, using higher-order approximations. These bounds are of higher-order accuracy and improve upon some of the known results in the literature. Finally, an interesting application to risk theory is discussed.     

ASYMPTOTIC BEHAVIOR OF SOLUTIONS TO THE COMPRESSIBLE NEMATIC LIQUID CRYSTAL SYSTEM IN R~3

In this paper,we study a nematic liquid crystals system in three-dimensional whole space R~3 and obtain the time decay rates of the higher-order spatial derivatives of the solution by the method of spectral analysis and energy estimates if the initial data belongs to L~1(R~3) additionally.     

A novel model of nonlocal thermoelasticity with time derivatives of higher order

The objective of this work is to introduce a new system of differential equations describing the nonlocal thermoelasticity theory with higher time derivatives and two-phase lags. In order to obtain this model, we used the nonlocal continuum theory proposed by Eringen and the methodology of the Taylor series expansion of higher-order time derivatives. Some generalized thermoelasticity theories follow as limited cases. This model is used to study the thermoelastic interaction in a nonlocal medium. The medium is exposed to an applied magnetic field and a periodic time heat source with a constant strength. Some comparisons have been displayed in figures to estimate the influences of the nonlocal parameter and magnetic field as well as the parameters of higher-order on all the field quantities.     

A necessary and sufficient condition for the existence of symmetric positive solutions of higher-order boundary value problems

A necessary and sufficient condition is obtained for the existence of symmetric positive solutions to higher-order nonlinear boundary value problems. The uniqueness, an iterative sequence and an error estimation for symmetric positive solutions are also discussed. Moreover, we give an example to illustrate the applicability of our results. Our analysis mainly relies on the monotone iterative technique.