Finite-element approximation for singularly perturbed nonlinear elliptic boundary-value problems

Numerical solution of a two-dimensional nonlinear singularly perturbed elliptic partial differential equation ∈ Δu = f(x, u), 0 < x, y < 1, with Dirichlet boundary condition is discussed here. The modified Newton method of third-order convergence is employed to linearize the nonlinear problem in place of the standard Newton method. The finite-element method is used to find the solution of the nonlinear differential equation. Numerical results are provided to demonstrate the usefulness of the method.     

Two-Scale Convergence Method and Scattering Problems

In this paper, we deal with the steady-state acoustic wave equation in the space ℝ³ diffracted by an obstacle made by an inhomogeneous medium and located in a bounded domain. The inhomogeneity of the medium depends on a parameter ε > 0. If the solution u _ε converges to a solution u ₀ of the limit problem as ε → 0, as in the homogenization process, then we can use the two-scale convergence method to study the convergence of the gradient.__________Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 10, Suzdal Conference-4, 2003.     

On the structure of biharmonic functions satisfying the clamped plate conditions on a right angle

Let u(r,θ) be biharmonic and bounded in the circular sector ¦θ¦ < π/4, 0 < r < ρ (ρ > 1) and vanish together with δu/δθ when ¦θ¦ = π/4. We consider the transform û(p,θ) = ∝₀¹r^{p − 1}u(r,θ)dr. We show that for any fixed θ₀ u(p,θ₀) is meromorphic with no real poles and cannot be entire unless u(r, θ₀) ≡ 0. It follows then from a theorem of Doetsch that u(r, θ₀) either vanishes identically or oscillates as r → 0.     

Asymptotic behaviour of nonlinear parabolic equations with critical exponents. A dynamical systems approach

We investigate the large-time behaviour of solutions to the nonlinear heat-conduction equation with absorption u_t = Δ(u^{σ + 1}) − u^β in Q = R^N × (0, ∞) (E) with N 1, σ > 0 and critical absorption exponent β = σ + 1 + 2/N; the initial function u(x, 0) = 0 is assumed to be integrable, nonnegative and compactly supported. We prove that u converges as t → ∞ to a unique self-similar function which is a contracted version of one of the asymptotic profiles of the nonabsorptive problem u_t = Δ(u^{σ + 1}), the same for any initial data. The cornerstone of the proof is a result about ω-limits of (infinite-dimensional) asymptotical dynamical systems. Combining this result with an asymptotic evaluation of the mass function as well as typical PDE estimates gives the behaviour of (E) for large times.Similar unusual asymptotic behaviour is obtained for the equation u_t = div(¦Du¦^σ Du) − u^β with same conditions on σ and u(x, 0) and critical value for β = σ + 1 + (σ + 2)/N.     

A Parameter-Uniform B-Spline Collocation Method for Singularly Perturbed Semilinear Reaction-Diffusion Problems

We consider a Dirichlet boundary value problem for a class of singularly perturbed semilinear reaction-diffusion equations. A  B-spline collocation method on a piecewise-uniform Shishkin mesh is developed to solve such problems numerically. The convergence analysis is given and the method is shown to be almost second-order convergent, uniformly with respect to the perturbation parameter ε in the maximum norm. Numerical results are presented to validate the theoretical results.     

Moment methods in two point Padé approximation

In the separable Hilbert space (H, ·, ·) the following “operator moment problem” is solved: given a complex sequence (c_k)_{k ε Z} generated by a meromorphic function f, find T ε B(H) and u₀ ε H such that T^ku₀, u₀ = c_k (k ε Z). If the sequence (c_k)_{k ε Z} is “normal,” an adapted form of Vorobyev's method of moments yields a sequence of two point Padé approximants to f. A sufficient condition for convergence of this sequence of approximants is given.     

Every T-space is equivalent to a T-space of continuous functions

A T-space U of degree k is a (k + 1)-dimensional vector space over (the real line) of real-valued functions defined on a linearly ordered set, satisfying the condition: for every nonzero u ε U, Z(u), the number of distinct zeros of u and ^-(u), the number of alternations in sign of u(t) with increasing t, each do not exceed k. It is demonstrated that given a T-space U of degree k > 0 on an arbitrary linearly ordered set T, there is a subset T′ of the real line and a nonsingular linear map L:U→ C(T′), the set of continuous functions on T′, such that the following hold: L(U) is a T-space of degree k; for u ε U, Z(u) = Z(L(u)), S^−(u) = S^−(L(u); and for some order-preserving bijection Θ:T → T′, u(t) = O if and only if L(u)(Θ(t) = 0. It is also shown that a T-space on a subset T can be extended to a T-space on the closure of T in ]inf T, sup T], provided that there are no “interval gaps” in T. Examples show that, in general, a T-space cannot be extended across an “interval gap” in its domain, and cannot be extended to both the infimum and supremum of its domain. Conditions for a T-space to be Markov, and to admit an adjoined function are derived.     

Reduction of the Gibbs phenomenon for smooth functions with jumps by the -algorithm

Recently, Brezinski has proposed to use Wynn's ε-algorithm in order to reduce the Gibbs phenomenon for partial Fourier sums of smooth functions with jumps, by displaying very convincing numerical experiments. In the present paper we derive analytic estimates for the error corresponding to a particular class of hypergeometric functions, and obtain the rate of column convergence for such functions, possibly perturbed by another sufficiently differentiable function. We also analyze the connection to Padé–Fourier and Padé–Chebyshev approximants, including those recently studied by Kaber and Maday.     

Characterization of generalized convex functions by their best approximation in sign-monotone norms

Let {u₀, u₁,… u_{n − 1}} and {u₀, u₁,…, u_n} be Tchebycheff-systems of continuous functions on [a, b] and let f ε C[a, b] be generalized convex with respect to {u₀, u₁,…, u_{n − 1}}. In a series of papers ([1], [2], [3]) D. Amir and Z. Ziegler discuss some properties of elements of best approximation to f from the linear spans of {u₀, u₁,…, u_{n − 1}} and {u₀, u₁,…, u_n} in the L_p-norms, 1 p ∞, and show (under different conditions for different values of p) that these properties, when valid for all subintervals of [a, b], can characterize generalized convex functions. Their methods of proof rely on characterizations of elements of best approximation in the L_p-norms, specific for each value of p. This work extends the above results to approximation in a wider class of norms, called “sign-monotone,” [6], which can be defined by the property: ¦ f(x)¦ ¦ g(x)¦,f(x)g(x) 0, a x b, imply f g . For sign-monotone norms in general, there is neither uniqueness of an element of best approximation, nor theorems characterizing it. Nevertheless, it is possible to derive many common properties of best approximants to generalized convex functions in these norms, by means of the necessary condition proved in [6]. For {u₀, u₁,…, u_n} an Extended-Complete Tchebycheff-system and f ε C⁽ⁿ⁾[a, b] it is shown that the validity of any of these properties on all subintervals of [a, b], implies that f is generalized convex. In the special case of f monotone with respect to a positive function u₀(x), a converse theorem is proved under less restrictive assumptions.     

Interpretation of the Lavrentiev phenomenon by relaxation

We consider functionals of the calculus of variations of the form F(u)= ∝₀¹ f(x, u, u′) dx defined for u ε W^1,∞(0, 1), and we show that the relaxed functional with respect to weak W^1,1(0, 1) convergence can be written as

Full-size image

, where the additional term L(u), called the Lavrentiev term, is explicitly identified in terms of F.     

The Complexity of Indefinite Elliptic Problems with Noisy Data

We study the complexity of second-order indefinite elliptic problems −div(au) +bu=f(with homogeneous Dirichlet boundary conditions) over ad-dimensional domain Ω, the error being measured in theH¹(Ω)-norm. The problem elementsfbelong to the unit ball ofW^{r, p}, (Ω), wherep [2, ∞] andr>d/p. Information consists of (possibly adaptive) noisy evaluations off,a, orb(or their derivatives). The absolute error in each noisy evaluation is at most δ. We find that thenth minimal radius for this problem is proportional ton^−r/d+ δ and that a noisy finite element method with quadrature (FEMQ), which uses only function values, and not derivatives, is a minimal error algorithm. This noisy FEMQ can be efficiently implemented using multigrid techniques. Using these results, we find tight bounds on the -complexity (minimal cost of calculating an -approximation) for this problem, said bounds depending on the costc(δ) of calculating a δ-noisy information value. As an example, if the cost of a δ-noisy evaluation isc(δ) = δ^−s(fors> 0), then the complexity is proportional to (1/)^{d/r + s}.     

Existence of three solutions for a class of elliptic eigenvalue problems

In this paper, we consider a problem of the type −Δu = λ(f(u) + μg(u)) in Ω, u¦∂Ω = 0, where Ω Rⁿ is an open-bounded set, f, g are continuous real functions on R, and λ, μ ε R. As an application of a new approach to nonlinear eigenvalues problems, we prove that, under suitable hypotheses, if ¦μ¦ is small enough, then there is some λ > 0 such that the above problem has at least three distinct weak solutions in W₀^1,2(Ω).     

Existence of maximal solutions

Let Ω be a plane bounded region. Let U = {U_μ(P):μ(P)εL∞(Ω), u_μ ε H²_{2, 0}(Ω) and a(P, μ(P))u_μ,xx + 2b(P, μ(P))u_μ,xy + c(P, μ(P))u_μ,vv = ƒ(P) for P ε Ω; here we are given a(P, X), b(P, X), c(P, X) ε L_∞(Ω × E¹), ƒ(P) ε L_p(Ω) with p > 2, and our partial differential equation is uniformly elliptic. The functions μ(P) are called profiles. We establish sufficient conditions—which when they apply are constructive—that there exist a μ₀ ε L_∞(Ω) such that u_μ0 (P) u_μ(P) for all P ε Ω and for each μ ε L_∞(Ω). Similar results are obtained for a difference equation and convergence is proved.     

The wave equation for Dunkl operators

Let k = (k_α)_αε, be a positive-real valued multiplicity function related to a root system , and Δ_k be the Dunkl-Laplacian operator. For (x, t) ε ^N, × , denote by u_k(x, t) the solution to the deformed wave equation Δ_ku_k,(x, t) = δ_ttu_k(x, t), where the initial data belong to the Schwartz space on ^N. We prove that for k 0 and N l, the wave equation satisfies a weak Huygens' principle, while a strict Huygens' principle holds if and only if (N − 3)/2 + Σ_αε+k_α ε . Here ⁺ is a subsystem of positive roots. As a particular case, if the initial data are supported in a closed ball of radius R > 0 about the origin, the strict Huygens principle implies that the support of u_k(x, t) is contained in the conical shell {(x, t), ε ^N × | |t| − R x |t| + R}. Our approach uses the representation theory of the group SL(2, ), and Paley-Wiener theory for the Dunkl transform. Also, we show that the (t-independent) energy functional of u_k is, for large |t|, partitioned into equal potential and kinetic parts.     

Singularly perturbed Dirichlet boundary-value problem for a stationary system in the linear elasticity theory

We consider a singularly perturbed Dirichlet boundary-value problem for an elliptic operator of the linear elasticity theory in a bounded domain with a small cavity. The main result is the proof of the theorem about the convergence of eigenelements of the perturbed boundary-value problem to eigenelements of the corresponding limiting boundary-value problem, when the parameter ? which defines the diameter of the small cavity tends to zero.     

A Schrödinger singular perturbation problem

Consider the equation −ε²Δu_ε + q(x)u_ε = f(u_ε) in , u(∞) < ∞, ε = const > 0. Under what assumptions on q(x) and f(u) can one prove that the solution u_ε exists and lim_ε→0u_ε = u(x), where u(x) solves the limiting problem q(x)u = f(u)? These are the questions discussed in the paper.     

A Nonlinear Dirichlet Problem on the Unit Ball

In the present paper, we consider the nonlinear Dirichlet problem - Δu(x) u^β(x) = 0 is the unit ball and q is a continuous radially symmetric function on B which may be singular on ?B. We state some mild conditions for the function q so that the Dirichlet problem has a positive classical solution.     

Solution of the Dirichlet problem for a singularly perturbed reaction-diffusion equation in a square on a Bakhvalov grid

The Dirichlet problem for a singularly perturbed reaction-diffusion equation in a square is solved with the help of the classic five-point difference scheme and a grid that is the tensor product of 1D Bakhvalov grids. Without imposing additional matching conditions in the corners of the domain, it is shown that the grid solution to the problem has the accuracy O(N ⁻²) in the norm L _∞ ^h , where N is the number of grid nodes along each direction. The accuracy is uniform with respect to a small parameter. A simulation confirms the theoretical prediction.     

Existence and regularity properties of non-isotropic singular elliptic equations

We establish existence and sharp regularity results for solutions to singular elliptic equations of the order u ^−β , 0 < β < 1, with gradient dependence and involving a forcing term λ f(x, u). Our approach is based on a singularly perturbed technique. We show that if the forcing parameter λ > 0 is large enough, our solution is positive. For λ small solutions vanish on a nontrivial set and therefore they exhibit free boundaries. We also establish regularity results for the free boundary and study the asymptotic behavior of the problem as b\searrow 0{\beta\searrow 0} and b\nearrow 1{\beta\nearrow 1}. In the former, we show that our solutions u _β converge to a C ^1,1 function which is a solution to an obstacle type problem. When b\nearrow 1{\beta\nearrow 1} we recover the Alt-Caffarelli theory.