The aim of this paper is to demonstrate how the variational equations from can be formulated and solved in some abstract Banach spaces without any a priori construction of special linearization schemes. This should be useful e.g. in the analysis of heat conduction problems and modelling of flow in porous media.

For $n=2m\ge 4$, let $\Omega \in {ℝ}^n$ be a bounded smooth domain and $𝒩\subset {ℝ}^L$ a compact smooth Riemannian manifold without boundary. Suppose that $\left\{u_k\right\}\in W^{m,2}(\Omega ,𝒩)$ is a sequence of weak solutions in the critical dimension to the perturbed $m$-polyharmonic maps $$\frac{\mathrm{d}}{\mathrm{d}t}{|}_{t=0}{E}_m\left(\Pi (u+t\xi )\right)=0$$ with ${\Phi }_k\to 0$ in ${\left(W^{m,2}(\Omega ,𝒩)\right)}^{*}$ and $u_k\rightharpoonup u$ weakly in $W^{m,2}(\Omega ,𝒩)$. Then $u$ is an $m$-polyharmonic map. In particular, the space of $m$-polyharmonic maps is sequentially compact for the weak-$W^{m,2}$ topology.

In this paper, we study a model for the magnetization in thin ferromagnetic films. It comes as a variational problem for $S^1$-valued maps $m^{\text{'}}$ (the magnetization) of two variables $x^{\text{'}}$: $E_{\epsilon }\left(m^{\text{'}}\right)=\epsilon \int {|{\nabla }^{\text{'}}·m^{\text{'}}|}^{2}d{x}^{\text{'}}+\frac{1}{2}\int {\left||{\nabla }^{\text{'}}{|}^{-1/2}{\nabla }^{\text{'}}·m^{\text{'}}\right|}^{2}d{x}^{\text{'}}$. We are interested in the behavior of minimizers as $\epsilon \to 0$. They are expected to be $S^1$-valued maps $m^{\text{'}}$ of vanishing distributional divergence ${\nabla }^{\text{'}}·m^{\text{'}}=0$, so that appropriate boundary conditions enforce line discontinuities. For finite $\epsilon >0$, these line discontinuities are approximated by smooth transition layers, the so-called Néel walls. Néel...

Domain decomposition techniques provide a powerful tool for the numerical approximation of partial differential equations. We focus on mortar finite element methods on non-matching triangulations. In particular, we discuss and analyze dual Lagrange multiplier spaces for lowest order finite elements. These non standard Lagrange multiplier spaces yield optimal discretization schemes and a locally supported basis for the associated constrained mortar spaces. As a consequence, standard efficient iterative...

Domain decomposition techniques provide a powerful tool for the numerical approximation of partial differential equations. We focus on mortar finite element methods on non-matching triangulations. In particular, we discuss and analyze dual Lagrange multiplier spaces for lowest order finite elements. These non standard Lagrange multiplier spaces yield optimal discretization schemes and a locally supported basis for the associated constrained mortar spaces. As a consequence, standard efficient iterative...

In this paper we study a unidirectional and elastic fiber composite. We use the homogenization method to obtain numerical results of the plane strain bulk modulus and the transverse shear modulus. The results are compared with the Hashin-Shtrikman bounds and are found to be close to the lower bounds in both cases. This indicates that the lower bounds might be used as a first approximation of the plane strain bulk modulus and the transverse shear modulus. We also point out the connection with the...

A lot of papers and books analyze analytical a posteriori error estimates from the point of view of robustness, guaranteed upper bounds, global efficiency, etc. At the same time, adaptive finite element methods have acquired the principal position among algorithms for solving differential problems in many physical and technical applications. In this survey contribution, we present and compare, from the viewpoint of adaptive computation, several recently published error estimation procedures for...

The aim of this paper is to compare and realize three efficient iterative methods, which have mesh independent convergence, and to propose some improvements for them. We look for the numerical solution of a nonlinear model problem using FEM discretization with gradient and Newton type methods. Three numerical methods have been carried out, namely, the gradient, Newton and quasi-Newton methods. We have solved the model problem with these methods, we have investigated the differences between them...

The development of iterative methods for solving linear algebraic equations has brought the question of when the employment of these methods is more advantageous than the use of the direct ones. In the paper, a comparison of the direct and iterative methods is attempted. The methods are applied to solving a certain class of boundary-value problems for elliptic partial differential equations which are used for the numerical modeling of electromagnetic fields in geophysics. The numerical experiments...

We compare numerical experiments from the String Gradient Weighted Moving Finite Element method and a Parabolic Moving Mesh Partial Differential Equation method, applied to three benchmark problems based on two different partial differential equations. Both methods are described in detail and we highlight some strengths and weaknesses of each method via the numerical comparisons. The two equations used in the benchmark problems are the viscous Burgers’ equation and the porous medium equation, both...

We prove a strong comparison principle for the solution of the Levi equation $L\left(u\right)=\sum {}_{i=1}{}^n\left(\right(1+u{}_t{}^2\left)\right(u{}_{x_i{x}_i}+u{}_{y_i{y}_i})+(u{}_{x_i}{}^2+u{}_{y_i}{}^2)u{}_{tt}+2(u{}_{y_i}-u{}_{x_i}u{}_t)u{}_{x_{i}t}-2(u{}_{x_i}+u{}_{y_i}u{}_t)u{}_{y_{i}t}+k(x,y,t\left)\right(1+\left|Du\right|{}^2){}^{3/2}=0$, applying Bony Propagation Principle.

The theory of compensated compactness of Murat and Tartar links the algebraic condition of rank-r convexity with the analytic condition of weak lower semicontinuity. The former is an algebraic condition and therefore it is, in principle, very easy to use. However, in applications of this theory, the need for an efficient classification of rank-r convex forms arises. In the present paper, we define the concept of extremal 2-forms  and characterize them in the rotationally invariant jointly...

In the present paper, we motivate and describe a numerical approach in order to detect the creation of fractures in a facet of a crystal evolving by anisotropic mean curvature. The result appears to be in accordance with the known examples of facet-breaking. Graphical simulations are included.