Search

In the paper (HA14), unfortunately, a computational error occurred in one estimate. Although the wrong estimate does not affect the main results, we want to present the necessary corrections. Essentially, Lemma 5.2 has to be corrected and, since it is used in the proof of Theorem 5.1, the proof of this theorem also requires an adaptation. (i) The corrected formulation of Lemma 5.2 is as follows. Lemma 5.2 For Lagrange finite elements with a shape-regular family of affine meshes { T n h } h>0 there is a constant C > 0 independent of q and h such that for all w ∈ Wh and q = 2m, m ∈N: CΛq−2 p (∇w,∇Ip h (wq−1))T ∫ T ‖∇w‖l2‖w‖ q−2 0,∞,T dx, ∀T ∈ T n h , (5.1) where Λp = ‖ ∑ndof i=1 |φi|‖0,∞,T is the Lebesgue constant.

n/a

Differential-algebraic equations (DAEs) are widely used for modeling of dynamical systems. In numerical analysis of DAEs, consistent initialization and index reduction are important preprocessing prior to numerical integration. Existing DAE solvers commonly adopt structural preprocessing methods based on combinatorial optimization. Unfortunately, the structural methods fail if the DAE has numerical or symbolic cancellations. For such DAEs, methods have been proposed to modify them to other DAEs to which the structural methods are applicable, based on the combinatorial relaxation technique. Existing modification methods, however, work only for a class of DAEs that are linear or close to linear. This paper presents two new modification methods for nonlinear DAEs: the substitution method and the augmentation method. Both methods are based on the combinatorial relaxation approach and are applicable to a large class of nonlinear DAEs. The substitution method symbolically solves equations for some derivatives based on the implicit function theorem and substitutes the solution back into the system. Instead of solving equations, the augmentation method modifies DAEs by appending new variables and equations. The augmentation method has advantages that the equation solving is not needed and the sparsity of DAEs is retained. It is shown in numerical experiments that both methods, especially the augmentation method, successfully modify high-index DAEs that the DAE solver in MATLAB cannot handle., Comment: A preliminary version of this paper is to appear in Proceedings of the 44th International Symposium on Symbolic and Algebraic Computation (ISSAC 2019), Beijing, China, July 2019

Olkin & Vaeth (1981) consider a two-way classification model with k rows and p columns in which the residuals (ei1, ..., eip) for row i are independently distributed with a multivariate normal distribution with zero means and common covariance matrix ((ij)). They concentrate their attention particularly on the situation where the row classification corresponds to k subjects, each treated in one of two or more ways; the column classification corresponds to p repetitions, perhaps in time, from each subject, probably resulting in correlated errors. Replication of subjects within treatments enables the elements of the covariance matrix to be estimated from the data, which in turn enables maximum likelihood estimates of the row and column parameters to be computed and likelihood ratio tests to be carried out without having to make assumptions about the form of the covariance matrix. The careful analysis of data of this kind has been the subject of numerous papers, Wishart (1938) perhaps being the first to concentrate on this topic. Despite this, however, the use of erroneous methods in published material is widespread and this provided the stimulus for the publication of our earlier paper (Rowell & Walters, 1976), and also for the present communication. We note here that there is confusion in Olkin & Vaeth's row/column terminology in their worked example: the rows in their tables are referred to as columns in the text, and vice versa. In what follows, when referring to their numerical example, the rows classification will refer to the high/low factor, and the columns classification to the time factor.

In 1993, Shub and Smale posed the problem of finding a sequence of univariate polynomials of degree $N$ with condition number bounded above by $N$. In a previous paper by C. Belt\'an, U. Etayo, J. Marzo and J. Ortega-Cerd\`a, it was proved that the optimal value of the condition number is of the form $O(\sqrt{N})$, and the sequence demanded by Shub and Smale was described by a closed formula (for large enough $N\geqslant N_0$ with $N_0$ unknown) and by a search algorithm for the rest of the cases. In this paper we find concrete estimates for the constant hidden in the $O(\sqrt{N})$ term and we describe a simple formula for a sequence of polynomials whose condition number is at most $N$, valid for all $N=4M^2$, with $M$ a positive integer., Comment: 21 pages

The process of grinding and polishing optical surfaces using a Computer Numerically Controlled machine produces a machine material removal profile. The profiles achievable by the machine depend on the nature of the tool used in the process, and the tool center motions. In this paper machine material removal profiles are developed as mollifications of given workpiece profiles for a variety of tool configurations. The form of the mollification, in effect, defines the tool center motion. Convergence of the machine's material removal profile to the given workpiece profile as the support of the tool goes to zero is established under mild assumptions. Numerical examples are included

The paper focuses on the analysis of the Euler projection Galerkin finite element method (FEM) for the dynamics of magnetization in ferromagnetic materials, described by the Landau–Lifshitz equation with the point-wise constraint $|{\textbf{m}}|=1$. The method is based on a simple sphere projection that projects the numerical solution onto a unit sphere at each time step, and the method has been used in many areas in the past several decades. However, error analysis for the commonly used method has not been done since the classical energy approach cannot be applied directly. In this paper we present an optimal $\textbf{L}^2$ error analysis of the backward Euler sphere projection method by using quadratic or higher order finite elements under a time step condition $\tau =O(\epsilon _0 h)$ with some small $\epsilon _0>0$. The analysis is based on more precise estimates of the extra error caused by the sphere projection in both $\textbf{L}^2$ and $\textbf{H}^1$ norms, and the classical estimate of dual norm. Numerical experiment is provided to confirm our theoretical analysis.

This paper is concerned with optimal error estimates and recovery technique of a classical mixed finite element method for the thermistor problem, which is governed by a parabolic/elliptic system with strong nonlinearity and coupling. The method is based on a popular combination of the lowest-order Raviart–Thomas mixed approximation for the electric potential/field $(\phi , \boldsymbol{\theta })$ and the linear Lagrange approximation for the temperature $u$. A common question is how the first-order approximation influences the accuracy of the second-order approximation to the temperature in such a strongly coupled system, while previous work only showed the first-order accuracy $O(h)$ for all three components in a traditional way. In this paper, we prove that the method produces the optimal second-order accuracy $O(h^2)$ for $u$ in the spatial direction, although the accuracy for the potential/field is in the order of $O(h)$. And more importantly, we propose a simple one-step recovery technique to obtain a new numerical electric potential/field of second-order accuracy. The analysis presented in this paper relies on an $H^{-1}$-norm estimate of the mixed finite element methods and analysis on a nonclassical elliptic map. We provide numerical experiments in both two- and three-dimensional spaces to confirm our theoretical analyses.

In this paper we revisit the problem of performing a $QZ$ step with a so-called ‘perfect shift’, which is an ‘exact’ eigenvalue of a given regular pencil $\lambda B-A$ in unreduced Hessenberg triangular form. In exact arithmetic, the $QZ$ step moves that eigenvalue to the bottom of the pencil, while the rest of the pencil is maintained in Hessenberg triangular form, which then yields a deflation of the given eigenvalue. But in finite precision the $QZ$ step gets ‘blurred’ and precludes the deflation of the given eigenvalue. In this paper we show that when we first compute the corresponding eigenvector to sufficient accuracy, then the $QZ$ step can be constructed using this eigenvector, so that the deflation is also obtained in finite precision. An important application of this technique is the computation of the index of a system of differential algebraic equations, since an exact deflation of the infinite eigenvalues is needed to impose correctly the algebraic constraints of such differential equations.

We propose and analyse numerical algorithms based on weighted least squares for the approximation of a real-valued function on a general bounded domain $\Omega \subset \mathbb{R}^d$. Given any $n$-dimensional approximation space $V_n \subset L^2(\Omega)$, the analysis in [6] shows the existence of stable and optimally converging weighted least-squares estimators, using a number of function evaluations $m$ of the order $n \log n$. When an $L^2(\Omega)$-orthonormal basis of $V_n$ is available in analytic form, such estimators can be constructed using the algorithms described in [6,Section 5]. If the basis also has product form, then these algorithms have computational complexity linear in $d$ and $m$. In this paper we show that, when $\Omega$ is an irregular domain such that the analytic form of an $L^2(\Omega)$-orthonormal basis is not available, stable and quasi-optimally weighted least-squares estimators can still be constructed from $V_n$, again with $m$ of the order $n \log n$, but using a suitable surrogate basis of $V_n$ orthonormal in a discrete sense. The computational cost for the calculation of the surrogate basis depends on the Christoffel function of $\Omega$ and $V_n$. Numerical results validating our analysis are presented., Comment: Version of the paper accepted for publication

In a series of papers published since 2017 the author introduced a simple alternative of the $n$-simplex type, to enhance the accuracy of approximations of second-order boundary value problems subject to Dirichlet boundary conditions, posed on smooth curved domains. This technique is based upon trial functions consisting of piecewise polynomials defined on straight-edged triangular or tetrahedral meshes, interpolating the Dirichlet boundary conditions at points of the true boundary. In contrast, the test functions are defined by the standard degrees of freedom associated with the underlying method for polytopic domains. While the mathematical analysis of the method for Lagrange and Hermite methods for two-dimensional second- and fourth-order problems was carried out in earlier paper by the author this paper is devoted to the study of the three-dimensional case. Well-posedness, uniform stability and optimal a priori error estimates in the energy norm are proved for a tetrahedron-based Lagrange family of finite elements. Novel error estimates in the $L^2$-norm, for the class of problems considered in this work, are also proved. A series of numerical examples illustrates the potential of the new technique. In particular, its superior accuracy at equivalent cost, as compared to the isoparametric technique, is highlighted.

Nonlinear elastic models are widely used to describe the elastic response of crystalline solids, for example, the well-known Cauchy–Born model. While the Cauchy–Born model only depends on the strain, effects of higher-order strain gradients are significant and higher-order continuum models are preferred in various applications such as defect dynamics and modeling of carbon nanotubes. In this paper we rigorously derive a higher-order nonlinear elasticity model for crystals from its atomistic description in one dimension. We show that, compared to the second-order accuracy of the Cauchy–Born model, the higher-order continuum model in this paper is of fourth-order accuracy with respect to the interatomic spacing in the thermal dynamic limit. In addition we discuss the key issues for the derivation of higher-order continuum models in more general cases. The theoretical convergence results are demonstrated by numerical experiments.

In this paper we analyze a class of trace finite element methods for the discretization of vector-Laplace equations. A key issue in the finite element discretization of such problems is the treatment of the constraint that the unknown vector field must be tangential to the surface (‘tangent condition’). We study three different natural techniques for treating the tangent condition, namely a consistent penalty method, a simpler inconsistent penalty method and a Lagrange multiplier method. The main goal of the paper is to present an analysis that reveals important properties of these three different techniques for treating the tangent constraint. A detailed error analysis is presented that takes the approximation of both the geometry of the surface and the solution of the partial differential equation into account. Error bounds in the energy norm are derived that show how the discretization error depends on relevant parameters such as the degree of the polynomials used for the approximation of the solution, the degree of the polynomials used for the approximation of the level set function that characterizes the surface, the penalty parameter and the degree of the polynomials used for the approximation of the Lagrange multiplier.

The Morley finite element method (FEM) is attractive for semilinear problems with the biharmonic operator as a leading term in the stream function vorticity formulation of two-dimensional Navier–Stokes problem and in the von Kármán equations. This paper establishes a best-approximation a priori error analysis and an a posteriori error analysis of discrete solutions close to an arbitrary regular solution on the continuous level to semilinear problems with a trilinear nonlinearity. The analysis avoids any smallness assumptions on the data, and so has to provide discrete stability by a perturbation analysis before the Newton–Kantorovich theorem can provide the existence of discrete solutions. An abstract framework for the stability analysis in terms of discrete operators from the medius analysis leads to new results on the nonconforming Crouzeix–Raviart FEM for second-order linear nonselfadjoint and indefinite elliptic problems with $L^\infty $ coefficients. The paper identifies six parameters and sufficient conditions for the local a priori and a posteriori error control of conforming and nonconforming discretizations of a class of semilinear elliptic problems first in an abstract framework and then in the two semilinear applications. This leads to new best-approximation error estimates and to a posteriori error estimates in terms of explicit residual-based error control for the conforming and Morley FEM.

Craig (1962) described an application of order statistics from a binomial distribution and gave an approximate expression for the mean value of the minimum of two binomial variables having the same probability of success and the same number of trials. In this note we suggest another approximation for the mean and also for the standard deviation of the minimum, using normal order statistics tabulated by Teichroew (1956).