Stability of solutions to differential equations with several variable delays. I

We consider a class of scalar linear differential equations with several variable delays and constant coefficients. A family of equations of the class is defined by coefficients and maximum admissible values of delays. We obtain conditions that are necessary and sufficient for the stability of solutions to all equations of the family. It is ascertained that the conditions are determined entirely by properties of the solution to the initial problem for an autonomous equation that belongs to the family. Some alternatives of required conditions are obtained in the form of estimates for solutions to autonomous equations in a finite interval.     

Plane wave approximation in linear elasticity

We consider the approximation of solutions of the time-harmonic linear elastic wave equation by linear combinations of plane waves. We prove algebraic orders of convergence both with respect to the dimension of the approximating space and to the diameter of the domain. The error is measured in Sobolev norms and the constants in the estimates explicitly depend on the problem wavenumber. The obtained estimates can be used in the h- and p-convergence analysis of wave-based finite element schemes.     

On finite time BV blow-up for the p-system

Abstract

The paper studies the possible blowup of the total variation for entropy weak solutions of the p-system, modeling isentropic gas dynamics. It is assumed that the density remains uniformly positive, while the initial data can have arbitrarily large total variation (measured in terms of Riemann invariants). Two main results are proved. (I) If the total variation blows up in finite time, then the solution must contain an infinite number of large shocks in a neighborhood of some point in the t-x plane. (II) Piecewise smooth approximate solutions can be constructed whose total variation blows up in finite time. For these solutions the strength of waves emerging from each interaction is exact, while rarefaction waves satisfy the natural decay estimates stemming from the assumption of genuine nonlinearity.     

一个刚性守恒律方程组的全隐式差分方法

１．引言本文研究如下非线性刚性守恒律方程组的全隐式差分逼近． 方程（1．１）中的源项ｇ（ｕ,ｖ）定义为 ｇ（ｕ,ｖ）＝ｖ－（１－μ）ｆ（ｕ）,（１．２）其中ｆ是ｕ的一个给定函数,δ是一个小正参数,称为松弛时间,μ是参数．方程组（１．１）频繁出现于粘弹性力学中． 在零松弛时间限（δ→０）下,从（１．１）可得到如下方程组该方程组通常称为“平衡”模型,而方程组（１．１）称为“非平衡”模型． 文中将假设μ满足 ０＜ μ＜ １,（１．４）以便保证拟稳定性条件［１９,２０］和次特征条件［１１,２,３］： λ１≤λ＊…     

A Study in the BV Space of a Denoising—Deblurring Variational Problem

In this paper we study, in the framework of functions of bounded variation, a general variational problem arising in image recovery, introduced in [3]. We prove the existence and the uniqueness of a solution using lower semicontinuity results for convex functionals of measures. We also give a new and fine characterization of the subdifferential of the functional, together with optimality conditions on the solution, using duality techniques of Temam for the theory of time-dependent minimal surfaces. We study the associated evolution equation in the context of nonlinear semigroup theory and we give an approximation result in continuous variables, using Γ -convergence. Finally, we discretize the problems by finite differences schemes and we present several numerical results for signal and image reconstruction. Accepted 5 March 2001. Online publication 20 June 2001.     

FINITE DIFFERENCE FRACTIONAL STEP METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

Characteristic finite difference fractional step schemes are put forward. The electric potential equation is described by a seven-point finite difference scheme, and the electron and hole concentration equations are treated by a kind of characteristic finite difference fractional step methods. The temperature equation is described by a fractional step method. Thick and thin grids are made use of to form a complete set. Piecewise threefold quadratic interpolation, symmetrical extension, calculus of variations, commutativity of operator product, decomposition of high order difference operators and prior estimates are also made use of. Optimal order estimates in l2 norm are derived to determine the error of the approximate solution. The well-known problem is thorongley and completely solred.     

ERROR ESTIMATES FOR SPARSE OPTIMAL CONTROL PROBLEMS BY PIECEWISE LINEAR FINITE ELEMENT APPROXIMATION

Optimization problems with L¹-control cost functional subject to an elliptic partial differential equation(PDE)are considered.However,different from the finite dimensiona l¹-regularization optimization,the resulting discretized L¹norm does not have a decoupled form when the standard piecewise linear finite element is employed to discretize the continuous problem.A common approach to overcome this difficulty is employing a nodal quadrature formula to approximately discretize the L¹-norm.In this paper,a new discretized scheme for the L¹-norm is presented.Compared to the new discretized scheme for L¹-norm with the nodal quadrature formula,the advantages of our new discretized scheme can be demonstrated in terms of the order of approximation.Moreover,finite element error estimates results for the primal problem with the new discretized scheme for the L¹-norm are provided,which confirms that this approximation scheme will not change the order of error estimates.To solve the new discretized problem,a symmetric Gauss-Seidel based majorized accelerated block coordinate descent(sGS-mABCD)method is introduced to solve it via its dual.The proposed sGS-mABCD algorithm is illustrated at two numerical examples.Numerical results not only confirm the finite element error estimates,but also show that our proposed algorithm is efficient.     

Approximation of the solution of a fourth order boundary value problem with nonsmooth coefficient

Summary A class of generalized finite element methods for the approximate solution of fourth order two point boundary value problem with nonsmooth coefficient is presented. The methods are based on the use of problem dependentL-splines incorporating the nonsmoothness of the coefficient. Stability is proved and optimal error estimates in theH ² norm are derived for the solution and postprocessed solution, under the assumption that the coefficient is of bounded variation. The relation of these methods to mixed methods is discussed.This research was sponsored by the Senate Research Committee of Syracuse University, Syracuse, NY 13210     

An adaptive time stepping method with efficient error control for second-order evolution problems

This work is concerned with time stepping fnite element methods for abstract second order evolution problems.We derive optimal order a posteriori error estimates and a posteriori nodal superconvergence error estimates using the energy approach and the duality argument.With the help of the a posteriori error estimator developed in this work,we will further propose an adaptive time stepping strategy.A number of numerical experiments are performed to illustrate the reliability and efciency of the a posteriori error estimates and to assess the efectiveness of the proposed adaptive time stepping method.     

Arch beam models: finite element analysis and superconvergence

Summary This paper studies finite element methods for a class of arch beam models. For both standard and mixed methods, existence and uniqueness results are proved, optimal rates of convergence are obtained and the superconvergence property is established. Reduced integration is shown to be an efficient method for arch beam problems and selected reduced integration is found to be identical to the mixed method. The significance of the analysis is threefold. The mixed method and the reduced integration methods converge uniformly at the optimal rate with respect to the arch thickness parameter, so they are locking free. Second, mixed method and reduced integration keep the superconvergence properties of the standard method. Finally, this is the first attempt to investigate the superconvergence of finite element methods for arch beam problems. We set up two types of superconvergence results: displacement at the nodal points and gradient at the Gauss points.This work was partially supported by the National Science Fundation grant CCR-88-20279     

一类非线性扩散方程解的blow up速率估计

利用F riedm an-M cleod方法和变动尺度方法研究了一类具有非线性边界条件的非线性扩散方程解的b low up问题,证明了解在有限时间b low up,并且得到了b low up速率估计.     

Numerical Solutions for second-kind Volterra integral equations by Galerkin methods

In this paper, we study the global convergence for the numerical solutions of nonlinear Volterra integral equations of the second kind by means of Galerkin finite element methods. Global superconvergence properties are discussed by iterated finite element methods and interpolated finite element methods. Local superconvergence and iterative correction schemes are also considered by iterated finite element methods. We improve the corresponding results obtained by collocation methods in the recent papers [6] and [9] by H. Brunner, Q. Lin and N. Yan. Moreover, using an interpolation post-processing technique, we obtain a global superconvergence of the O(h ^2r )-convergence rate in the piecewise-polynomial space of degree not exceeding (r–1). As a by-product of our results, all these higher order numerical methods can also provide an a posteriori error estimator, which gives critical and useful information in the code development.     

不可压缩流动的并行数值方法

不可压缩流动的数值模拟是计算流体力学的重要组成部分. 基于有限元离散方法, 本文设计了不可压缩Navier-Stokes (N-S)方程支配流的若干并行数值算法. 这些并行算法可归为两大类: 一类是基于两重网格离散方法, 首先在粗网格上求解非线性的N-S方程, 然后在细网格的子区域上并行求解线性化的残差方程, 以校正粗网格的解; 另一类是基于新型完全重叠型区域分解技巧, 每台处理器用一局部加密的全局多尺度网格计算所负责子区域的局部有限元解. 这些并行算法实现简单, 通信需求少, 具有良好的并行性能, 能获得与标准有限元方法相同收敛阶的有限元解. 理论分析和数值试验验证了并行算法的高效性     

Stabilized nodal‐based finite element methods for the magnetic induction problem

We consider the time‐dependent magnetic induction model as a step towards the resistive magnetohydrodynamics model in incompressible media. Conforming nodal‐based finite element approximations of the induction model with inf‐sup stable finite elements for the magnetic field and the magnetic pseudo‐pressure are investigated. Based on a residual‐based stabilization technique proposed by Badia and Codina, SIAM J. Numer. Anal. 50 (2012), pp. 398–417, we consider a stabilized nodal‐based finite element method for the numerical solution. Error estimates are given for the semi‐discrete model in space. Finally, we present some examples, for example, for the magnetic flux expulsion problem, Shercliff's test case and singular solutions of the Maxwell problem. Copyright © 2015 John Wiley & Sons, Ltd.     

A constructive method for deriving finite elements of nodal type

Summary In this paper, we propose an algorithm to derive nodal methods corresponding to various two and three-dimensional nonconforming and mixed finite elements. We show that this algorithm can be used to obtain several classical schemes as well as some more recently developed schemes, and that it leads to a simple proof of unisolvence for these methods. Finally we use our method to obtain a three dimensional nodal scheme of BDM type.     

A spatial shear deformable beam finite element based on the absolute nodal coordinate formulation

In the present paper a three-dimensional beam finite element undergoing large deformations is proposed. Since the definition of the proposed finite element is based on the absolute nodal coordinate formulation (ANCF), no rotational coordinates occur in the formulation. In the current approach, the orientation of the cross section is parameterized by means of slope vectors. Since those are no unit vectors, the cross-section can deform, similar to existing thick beam and shell elements. The nodal displacements and the directional derivatives of the displacements are chosen as nodal coordinates, but in contrast to standard ANCF elements, the proposed formulation is based on the two transversal slope vectors per node only. Different approaches for the virtual work of elastic forces are presented: a continuum mechanics based formulation, as well as a structural mechanics based formulation, which is in accordance with classical nonlinear beam finite elements. Since different interpolation functions as in standard ANCF elements are used, a much better convergence rate (up to order four) can be obtained. Therefore, the present element has high potential for application in geometrically nonlinear problems. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of consistent and inconsistent finite difference schemes and an acceleration technique

This paper states and generalizes in part some recent results on finite difference methods for Dirichlet problems in a bounded domain Ω which the author has obtained by himself or with coworkers. After stating a superconvergence property of finite difference solution for the case where the exact solution u belongs to , it is remarked that such a property does not hold in general if . Next, a convergence theorem is given for inconsistent schemes under some assumptions. Furthermore, it is shown that the accuracy of the approximate solution can be improved by a coordinate transformation. Numerical examples are also given.     

Interior maximum norm estimates for finite element discretizations of the Stokes equations

Interior estimates are proved in the L ^∞ norm for stable finite element discretizations of the Stokes equations on translation invariant meshes. These estimates yield information about the quality of the finite element solution in subdomains a positive distance from the boundary. While they have been established for second-order elliptic problems, these interior, or local, maximum norm estimates for the Stokes equations are new. By applying finite differenciation methods on a translation invariant mesh, we obtain optimal convergence rates in the mesh size h in the maximum norm. These results can be used for analyzing superconvergence in finite element methods for the Stokes equations.     

Variational methods for a singular SPDE yielding the universality of the magnetization ripple

The magnetization ripple is a microstructure formed in thin ferromagnetic films. It can be described by minimizers of a nonconvex energy functional leading to a nonlocal and nonlinear elliptic SPDE in two dimensions driven by white noise, which is singular. We address the universal character of the magnetization ripple using variational methods based on Γ-convergence. Due to the infinite energy of the system, the (random) energy functional has to be renormalized. Using the topology of Γ-convergence, we give a sense to the law of the renormalized functional that is independent of the way white noise is approximated. More precisely, this universality holds in the class of (not necessarily Gaussian) approximations to white noise satisfying the spectral gap inequality, which allows us to obtain sharp stochastic estimates. As a corollary, we obtain the existence of minimizers with optimal regularity.     

On the convergence of finite state mean-field games through Γ-convergence

In this study, we consider the long-term convergence (trend toward an equilibrium) of finite state mean-field games using Γ-convergence. Our techniques are based on the observation that an important class of mean-field games can be viewed as the Euler–Lagrange equation of a suitable functional. Therefore, using a scaling argument, one can convert a long-term convergence problem into a Γ-convergence problem. Our results generalize previous results related to long-term convergence for finite state problems.