Elliptic Differential Equations Research Articles

Global bifurcation for monotone fronts of elliptic equations

We present two results on global continuation of monotone front-type solutions to elliptic PDEs posed on infinite cylinders. This is done under quite general assumptions, and in particular applies even to fully nonlinear equations as well as quasilinear problems with transmission boundary conditions. Our approach is rooted in the analytic global bifurcation theory of Dancer (1973) and Buffoni–Toland (2003), but extending it to unbounded domains requires contending with new potential limiting behavior relating to loss of compactness. We obtain an exhaustive set of alternatives for the global behavior of the solution curve that is sharp, with each possibility having a direct analogue in the bifurcation theory of second-order ODEs.As a major application of the general theory, we construct global families of internal hydrodynamic bores. These are traveling front solutions of the full two-phase Euler equation in two dimensions. The fluids are confined to a channel that is bounded above and below by rigid walls, with incompressible and irrotational flow in each layer. Small-amplitude fronts for this system have been obtained by several authors. We give the first large-amplitude result in the form of continuous curves of elevation and depression bores. Following the elevation curve to its extreme, we find waves whose interfaces either overturn (develop a vertical tangent) or become exceptionally singular in that the flow in both layers degenerates at a single point on the boundary. For the curve of depression waves, we prove that either the interface overturns or it comes into contact with the upper wall.

A closure theorem for $\Gamma$-convergence and $H$-convergence with applications to non-periodic homogenization

In this work we examine the stability of some classes of integrals, in particular with respect to homogenization. The prototypical case is the homogenization of quadratic energies with periodic coefficients perturbed by a term vanishing at infinity, which has recently been examined in the framework of elliptic PDEs. We use localization techniques and higher-integrability Meyers-type results to provide a closure theorem by \Gamma -convergence within a large class of integral functionals. From such a result we derive stability theorems in homogenization which comprise the case of perturbations with zero average on the whole space. The results are also extended to the stochastic case, and specialized to the G -convergence of operators corresponding to quadratic forms. A corresponding analysis is also carried out for non-symmetric operators using the localization properties of H -convergence. Finally, we treat the case of perforated domains with Neumann boundary condition, and their stability.

Accelerating the convergence of Newton’s method for nonlinear elliptic PDEs using Fourier neural operators

It is well known that Newton’s method can have trouble converging if the initial guess is too far from the solution. Such a problem particularly occurs when this method is used to solve nonlinear elliptic partial differential equations (PDEs) discretized via finite differences. This work focuses on accelerating Newton’s method convergence in this context. We seek to construct a mapping from the parameters of the nonlinear PDE to an approximation of its discrete solution, independently of the mesh resolution. This approximation is then used as an initial guess for Newton’s method. To achieve these objectives, we elect to use a Fourier neural operator (FNO). The loss function is the sum of a data term (i.e., the comparison between known solutions and outputs of the FNO) and a physical term (i.e., the residual of the PDE discretization). Numerical results, in one and two dimensions, show that the proposed initial guess accelerates the convergence of Newton’s method by a large margin compared to a naive initial guess, especially for highly nonlinear and anisotropic problems, with larger gains on coarse grids.

An extension of the approximate component mode synthesis method to the heterogeneous Helmholtz equation

Abstract In this work, we propose and analyze an extension of the approximate component mode synthesis (ACMS) method to the two-dimensional heterogeneous Helmholtz equation. The ACMS method has originally been introduced by Hetmaniuk and Lehoucq as a multiscale method to solve elliptic partial differential equations. The ACMS method uses a domain decomposition to separate the numerical approximation by splitting the variational problem into two independent parts: local Helmholtz problems and a global interface problem. While the former are naturally local and decoupled such that they can be easily solved in parallel, the latter requires the construction of suitable local basis functions relying on local eigenmodes and suitable extensions. We carry out a full error analysis of this approach focusing on the case where the domain decomposition is kept fixed, but the number of eigenfunctions is increased. The theoretical results in this work are supported by numerical experiments verifying algebraic convergence for the method. In certain, practically relevant cases, even super-algebraic convergence for the local Helmholtz problems can be achieved without oversampling.

Sobolev regularity of bivariate isogeometric finite element spaces in case of a geometry map with degenerate corner

We investigate Sobolev regularity of bivariate functions obtained in Isogeometric Analysis when using geometry maps that are degenerate in the sense that the first partial derivatives vanish at isolated points. In particular, we show how the known C1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$C^1$$\\end{document}-conditions for D-patches have to be tightened to guarantee square integrability of second partial derivatives, as required when computing finite element approximations of elliptic fourth order PDEs like the biharmonic equation.

Nonexistence of the compressible Euler equations with space-dependent damping in high dimensions

Abstract Compressible Euler equations with space-dependent damping in high dimensions R n ( n = 2 , 3 ) {{\bf{R}}}^{n}\hspace{0.33em}\hspace{0.33em}\left(n=2,3) are considered in this article. Assuming that the small initial velocity and small perturbation of the initial density have compact support, we establish finite-time blow-up results for the Euler system, by combining energy estimate and new test functions constructed by the solutions of the following linear elliptic partial differential equations system: − G 1 ( x ) + ∇ ⋅ G 2 → ( x ) = 0 , − G 2 → ( x ) + ∇ G 1 ( x ) = μ G 2 → ( x ) ( 1 + ∣ x ∣ ) λ . \left\{\begin{array}{l}-{G}_{1}\left(x)+\nabla \cdot \overrightarrow{{G}_{2}}\left(x)=0,\\ -\overrightarrow{{G}_{2}}\left(x)+\nabla {G}_{1}\left(x)=\frac{\mu \overrightarrow{{G}_{2}}\left(x)}{{(1+| x| )}^{\lambda }}.\end{array}\right. This result generalizes the one in the literature from 1 − D 1-D to high dimension R n ( n = 2 , 3 ) {{\bf{R}}}^{n}\hspace{0.33em}\hspace{0.33em}\left(n=2,3) .

Convergence of an adaptive C0-interior penalty Galerkin method for the biharmonic problem

Abstract We develop a basic convergence analysis for an adaptive ${C}^{0}\mathsf{IPG}$ method for the biharmonic problem that provides convergence without rates for all practically relevant marking strategies and all penalty parameters assuring coercivity of the method. The analysis hinges on embedding properties of (broken) Sobolev and BV spaces, and the construction of a suitable limit space. In contrast to the previous convergence result of adaptive discontinuous Galerkin methods for elliptic PDEs, here we have to deal with the fact that the Lagrange finite element spaces may possibly contain no proper $C^{1}$-conforming subspace. This prevents from a straight forward generalization and requires the development of some new key technical tools.

Overcoming the curse of dimensionality in the numerical approximation of high-dimensional semilinear elliptic partial differential equations

Recently, so-called full-history recursive multilevel Picard (MLP) approximation schemes have been introduced and shown to overcome the curse of dimensionality in the numerical approximation of semilinear parabolic partial differential equations (PDEs) with Lipschitz nonlinearities. The key contribution of this article is to introduce and analyze a new variant of MLP approximation schemes for certain semilinear elliptic PDEs with Lipschitz nonlinearities and to prove that the proposed approximation schemes overcome the curse of dimensionality in the numerical approximation of such semilinear elliptic PDEs.

Error analysis of kernel/GP methods for nonlinear and parametric PDEs

We introduce a priori Sobolev-space error estimates for the solution of arbitrary nonlinear, and possibly parametric, PDEs that are defined in the strong sense, using Gaussian process and kernel based methods. The primary assumptions are: (1) a continuous embedding of the reproducing kernel Hilbert space of the kernel into a Sobolev space of sufficient regularity; and (2) the stability of the differential operator and the solution map of the PDE between corresponding Sobolev spaces. The proof is articulated around Sobolev norm error estimates for kernel interpolants and relies on the minimizing norm property of the solution. The error estimates demonstrate dimension-benign convergence rates if the solution space of the PDE is smooth enough. We illustrate these points with applications to high-dimensional nonlinear elliptic PDEs and parametric PDEs. Although some recent machine learning methods have been presented as breaking the curse of dimensionality in solving high-dimensional PDEs, our analysis suggests a more nuanced picture: there is a trade-off between the regularity of the solution and the presence of the curse of dimensionality. Therefore, our results are in line with the understanding that the curse is absent when the solution is regular enough.

The Extended Weierstrass Transformation Method for the Biswas–Arshed Equation with Beta Time Derivative

In this article, exact solutions of the Biswas–Arshed equation are obtained using the extended Weierstrass transformation method (EWTM). This method is widely used in solid-state physics, electrodynamics, and mathematical physics, and it yields exact solution functions involving trigonometric, rational trigonometric, Weierstrass elliptic, wave, and rational functions. The process involves expanding the solution functions of an elliptic differential equation into finite series by transforming them into Weierstrass functions. Furthermore, it generates parametric solutions for nonlinear algebraic equation systems, which are particularly useful in mathematical physics. These solutions are derived using the Mathematica package program. To analyze the behavior of these determined solution functions, the article employs separate two- and three-dimensional graphs showing the real and imaginary components, along with contour and density graphs. These visuals aid in comprehending the physical characteristics exhibited by these solution functions.

Qualitative analysis of optimisation problems with respect to non-constant Robin coefficients

Following recent interest in the qualitative analysis of some optimal control and shape optimisation problems, we provide in this article a detailed study of the optimisation of Robin boundary conditions in PDE constrained calculus of variations. Our main model consists of an elliptic PDE of the form −∆u β = f (x, u β) endowed with the Robin boundary conditions ∂ν u β +β(x)u β = 0. The optimisation variable is the function β, which is assumed to take values between 0 and 1 and to have a fixed integral. Two types of criteria are under consideration: the first one is non-energetic criteria. In other words, we aim at optimising functionals of the form J (β) =\int_{Ω or ∂Ω} j(u β). We prove that, depending on the monotonicity of the function j, the optimisers may be of bang-bang type (in other words, the optimisers write 1Γ for some measurable subset Γ of ∂Ω) or, on the contrary, that they may only take values strictly between 0 and 1. This has consequence for a related shape optimisation problem, in which one tries to find where on the boundary Neumann (∂ν u = 0) and constant Robin conditions (∂ν u + u = 0) should be placed in order to optimise criteria. The proofs for this first case rely on new fine oscillatory techniques, used in combination with optimality conditions. We then investigate the case of compliance-type functionals. For such energetic functionals, we give an in-depth analysis and even some explicit characterisation of optimal β * .

Box Constraints and Weighted Sparsity Regularization for Identifying Sources in Elliptic PDEs

We explore the possibility for using boundary data to identify sources in elliptic PDEs. Even though the associated forward operator has a large null space, it turns out that box constraints, combined with weighted sparsity regularization, can enable rather accurate recovery of sources with constant magnitude/strength. In addition, for sources with varying strength, the support of the inverse solution will be a subset of the support of the true source. We present both an analysis of the problem and a series of numerical experiments. Our work only addresses discretized problems. The reason for introducing the weighting procedure is that standard (unweighted) sparsity regularization fails to provide adequate results for the source identification task considered in this paper. This investigation is also motivated by applications, e.g. recovering mass distributions from measurements of gravitational fields and inverse scattering. We develop the methodology and the analysis in terms of Euclidean spaces, and our results can therefore be applied to many problems. For example, the results are equally applicable to models involving the screened Poisson equation as to models using the Helmholtz equation, with both large and small wave numbers.

The operational matrices for Elliptic Partial Differential Equations with mixed boundary conditions

Abstract The purpose of this research is to implement the orthogonal polynomials associated with operational matrices to get the approximate solutions for solving two-dimensional elliptic partial differential equations (E-PDEs) with mixed boundary conditions. The orthogonal polynomials are based on the Standard polynomial (xi ), Legendre, Chebyshev, Bernoulli, Boubaker, and Genocchi polynomials. This study focuses on constructing quick and precise analytic approximations using a simple, elegant, and potent technique based on an orthogonal polynomial representation of the solution as a double power series. Consequently, a linear partial differential equation is transformed into a linear algebraic system which is solved by the Mathematica®12. Approximate solutions can be found if the answers are polynomials in and of itself. Three applications involving well-known linear problems Laplace, Poisson, and Helmholtz equations have been solved by using the proposed methods, and a comparison of the approaches has been provided. Furthermore, the computation of the error norm L∞ has been done to show the accuracy of the suggested approaches. The results clearly demonstrate how precise, efficient, and dependable the proposed methods are in obtaining rough solutions to the problem. Bernoulli was one of the best methods in most examples.

Застосування методу двобічних наближень до знаходження додатних аксіально-симетричних розв’язків крайових задач з монотонними нелінійностями

This paper considers the application of the method of two-sided approximations for finding positive solutions to boundary value problems for nonlinear elliptic differential equations with axial symmetry. The case of Dirichlet boundary conditions (first kind) is considered, with a power-type monotone nonlinearity, where the exponent ranges from 0 to 1. By transitioning to polar coordinates in the boundary value problem for the elliptic equation, due to the axial symmetry of the solution, the problem is reduced to a boundary value problem for an ordinary differential equation on an interval with respect to a function that depends only on the polar radius, thus eliminating the dependence on the angle. The pole of the polar coordinate system becomes a singular point, where a boundedness condition must be imposed. For the boundary value problem, a Green's function is constructed to further reduce the problem to a Hammerstein integral equation. The integral equation is treated as a nonlinear operator equation in a Banach space of continuous functions on the interval, partially ordered by a cone of non-negative continuous functions on this interval. The operator is examined for properties such as monotonicity (isotonicity), positivity, boundedness, and concavity. Next, the initial approximation is found as the endpoints of the invariant conical segment for the isotone operator in such a way as to ensure the highest convergence rate of the iterative process. The following iterative sequences of two-sided approximations are constructed: the first sequence, which is non-decreasing with respect to the cone, and the second sequence, which is non-increasing with respect to the cone. At each iteration, the arithmetic mean of the upper and lower approximations is taken as the next approximation. The iterative process is terminated when the error estimate of the solution satisfies the specified accuracy. The theoretical results obtained in this work were validated through a computational experiment. The dependence of the solution on the parameters in the right-hand side was analyzed and illustrated with corresponding graphs

Adaptive reduced basis trust region methods for parameter identification problems

In this contribution, we are concerned with model order reduction in the context of iterative regularization methods for the solution of inverse problems arising from parameter identification in elliptic partial differential equations. Such methods typically require a large number of forward solutions, which makes the use of the reduced basis method attractive to reduce computational complexity. However, the considered inverse problems are typically ill-posed due to their infinite-dimensional parameter space. Moreover, the infinite-dimensional parameter space makes it impossible to build and certify classical reduced-order models efficiently in a so-called “offline phase”. We thus propose a new algorithm that adaptively builds a reduced parameter space in the online phase. The enrichment of the reduced parameter space is naturally inherited from the Tikhonov regularization within an iteratively regularized Gauß-Newton method. Finally, the adaptive parameter space reduction is combined with a certified reduced basis state space reduction within an adaptive error-aware trust region framework. Numerical experiments are presented to show the efficiency of the combined parameter and state space reduction for inverse parameter identification problems with distributed reaction or diffusion coefficients.

Variational methods for some singular stochastic elliptic PDEs

Variational methods for some singular stochastic elliptic PDEs

Derivation of analytic formulae for several resonance frequencies of the SparkJet actuator

Building on previous research that emphasized the importance of focusing on resonance frequencies for a fundamental understanding of the thrust and complex internal flow phenomena of the SparkJet actuators, this study theoretically derives analytic formulae for several important resonance frequencies. The research addresses the typical configuration of the SparkJet actuators, which consists of a cylindrical cavity and orifice connected by a conical converging nozzle. While the resonance frequencies of the SparkJet actuators can be obtained by solving the eigenvalue problem presented in previous studies, this eigenvalue problem, despite being a typical boundary value problem in the form of an elliptic partial differential equation, is challenging to solve using conventional numerical methods such as iterative methods, because the eigenvalue is included as an unknown in the operator. Consequently, Boundary Element Method (BEM) or methods using the Green functions have been proposed to obtain numerical solutions, but these still require handling large matrix data, resulting in significant computational costs and memory consumption. To overcome this, the current study first simplifies the eigenvalue problem using the conformal mapping presented in previous research. Then, referencing prior studies that claim that three types of resonance frequencies (Helmholtz resonance frequency and resonance frequencies representing streamwise or radial directional oscillation) significantly affect the thrust of the SparkJet actuators, characteristics of these frequencies are mathematically defined. Using these mathematical characteristics, the study derives the asymptotic and approximate solutions of the simplified eigenvalue problem, from which the resonance frequencies are obtained. The analytic formulae proposed in this study theoretically explain the already known geometric tendencies of the resonance frequencies and reveal new geometric factors influencing resonance frequencies, which were previously unknown. This approach is expected to facilitate obtaining several important resonance frequencies of the SparkJet actuator more promptly and accurately and to provide a deeper understanding of the nature of the complex oscillation phenomena inside the actuator.

Homogenisation of nonlinear Dirichlet problems in randomly perforated domains under minimal assumptions on the size of perforations

AbstractIn this paper we study the convergence of nonlinear Dirichlet problems for systems of variational elliptic PDEs defined on randomly perforated domains of $$\mathbb {R}^n$$ R n . Under the assumption that the perforations are small balls whose centres and radii are generated by a stationary short-range marked point process, we obtain in the critical-scaling limit an averaged nonlinear analogue of the extra term obtained in the classical work of Cioranescu and Murat (Res Notes Math III, 1982). In analogy to the random setting recently introduced by Giunti, Höfer and Velázquez (Commun Part Differ Equ 43(9):1377–1412, 2018) to study the Poisson equation, we only require that the random radii have finite $$(n-q)$$ ( n - q ) -moment, where $$1<q<n$$ 1 < q < n is the growth-exponent of the associated energy functionals. This assumption on the one hand ensures that the expectation of the nonlinear q-capacity of the spherical holes is finite, and hence that the limit problem is well defined. On the other hand, it does not exclude the presence of balls with large radii, that can cluster up. We show however that the critical rescaling of the perforations is sufficient to ensure that no percolating-like structures appear in the limit.

Convexity of solutions to elliptic PDE's

This article concerns the convexity or concavity of solutions to special second order elliptic partial differential equations in convex domains. We concentrate our investigation to boundary blow up solutions as well as to solutions of particular singular equations. Following a method due to Korevaar and Kennington, we find a new sufficient condition for proving convexity or concavity. This sufficient condition is useful when the semilinear component of the equation is the sum of two or more terms. For more information see https://ejde.math.txstate.edu/Volumes/2024/51/abstr.html

On the regularity of optimal potentials in control problems governed by elliptic equations

Abstract In this paper we consider optimal control problems where the control variable is a potential and the state equation is an elliptic partial differential equation of Schrödinger type, governed by the Laplace operator. The cost functional involves the solution of the state equation and a penalization term for the control variable. While the existence of an optimal solution simply follows by the direct methods of the calculus of variations, the regularity of the optimal potential is a difficult question and under the general assumptions we consider, no better regularity than the BV \mathrm{BV} one can be expected. This happens in particular for the cases in which a bang-bang solution occurs, where optimal potentials are characteristic functions of a domain. We prove the BV \mathrm{BV} regularity of optimal solutions through a regularity result for PDEs. Some numerical simulations show the behavior of optimal potentials in some particular cases.