Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA new formulation for the numerical proof of the existence of solutions to elliptic problems
84
0
0.0
(
0
)
Added by
Kouta Sekine
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the infinite Newton-type fixed point equation $w = - {mathcal L}^{-1} {mathcal F}(hat{u}) + {mathcal L}^{-1} {mathcal G}(w)$, where ${mathcal L}$ is a linearized operator, ${mathcal F}(hat{u})$ is a residual, and ${mathcal G}(w)$ is a local Lipschitz term. Therefore, the estimations of $| {mathcal L}^{-1} {mathcal F}(hat{u}) |$ and $| {mathcal L}^{-1}{mathcal G}(w) |$ play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement for matrix problems, we represent the inverse operator ${mathcal L}^{-1}$ as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as ${mathcal L}^{-1}$ are presented in the appendix.
rate research
Read More
The existence and uniqueness of weak solutions to dynamical low-rank evolution problems for parabolic partial differential equations in two spatial dimensions is shown, covering also non-diagonal diffusion in the elliptic part. The proof is based on a variational time-stepping scheme on the low-rank manifold. Moreover, this scheme is shown to be closely related to practical methods for computing such low-rank evolutions.
In this paper, we present a numerical verification method of solutions for nonlinear parabolic initial boundary value problems. Decomposing the problem into a nonlinear part and an initial value part, we apply Nakaos projection method, which is based on the full-discrete finite element method with constructive error estimates, to the nonlinear part and use the theoretical analysis for the heat equation to the initial value part, respectively. We show some verified examples for solutions of nonlinear problems from initial value to the neighborhood of the stationary solutions, which confirm us the actual effectiveness of our method.
This paper addresses the question whether there are numerical schemes for constant-coefficient advection problems that can yield convergent solutions for an infinite time horizon. The motivation is that such methods may serve as building blocks for long-time accurate solutions in more complex advection-dominated problems. After establishing a new notion of convergence in an infinite time limit of numerical methods, we first show that linear methods cannot meet this convergence criterion. Then we present a new numerical methodology, based on a nonlinear jet scheme framework. We show that these methods do satisfy the new convergence criterion, thus establishing that numerical methods exist that converge on an infinite time horizon, and demonstrate the long-time accuracy gains incurred by this property.
A variational formulation is introduced for the Oseen equations written in terms of vor-ti-city and Bernoulli pressure. The velocity is fully decoupled using the momentum balance equation, and it is later recovered by a post-process. A finite element method is also proposed, consisting in equal-order Nedelec finite elements and piecewise continuous polynomials for the vorticity and the Bernoulli pressure, respectively. The {it a priori} error analysis is carried out in the $mathrm{L}^2$-norm for vorticity, pressure, and velocity; under a smallness assumption either on the convecting velocity, or on the mesh parameter. Furthermore, an {it a posteriori} error estimator is designed and its robustness and efficiency are studied using weighted norms. Finally, a set of numerical examples in 2D and 3D is given, where the error indicator serves to guide adaptive mesh refinement. These tests illustrate the behaviour of the new formulation in typical flow conditions, and they also confirm the theoretical findings.
In this paper, an important discovery has been found for nonconforming immersed finite element (IFE) methods using integral-value degrees of freedom for solving elliptic interface problems. We show that those IFE methods can only achieve suboptimal convergence rates (i.e., $O(h^{1/2})$ in the $H^1$ norm and $O(h)$ in the $L^2$ norm) if the tangential derivative of the exact solution and the jump of the coefficient are not zero on the interface. A nontrivial counter example is also provided to support our theoretical analysis. To recover the optimal convergence rates, we develop a new nonconforming IFE method with additional terms locally on interface edges. The unisolvence of IFE basis functions is proved on arbitrary triangles. Furthermore, we derive the optimal approximation capabilities of both the Crouzeix-Raviart and the rotated-$Q_1$ IFE spaces for interface problems with variable coefficients via a unified approach different from multipoint Taylor expansions. Finally, optimal error estimates in both $H^1$- and $L^2$- norms are proved and confirmed with numerical experiments.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
