 Award ID(s):
 2208159
 NSFPAR ID:
 10407649
 Date Published:
 Journal Name:
 Acta Mathematica Vietnamica
 ISSN:
 02514184
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

We propose a new iterative scheme to compute the numerical solution to an overdetermined boundary value problem for a general quasilinear elliptic PDE. The main idea is to repeatedly solve its linearization by using the quasireversibility method with a suitable Carleman weight function. The presence of the Carleman weight function allows us to employ a Carleman estimate to prove the convergence of the sequence generated by the iterative scheme above to the desired solution. The convergence of the iteration is fast at an exponential rate without the need of an initial good guess. We apply this method to compute solutions to some general quasilinear elliptic equations and a large class of firstorder HamiltonJacobi equations. Numerical results are presented.more » « less

This paper aims to reconstruct the initial condition of a hyperbolic equation with an unknown damping coeﬃcient. Our approach involves approximating the hyperbolic equation’s solution by its truncated Fourier expansion in the time domain and using the recently developed polynomialexponential basis. This truncation process facilitates the elimination of the time variable, consequently, yielding a system of quasilinear elliptic equations. To globally solve the system without needing an accurate initial guess, we employ the Carleman contraction principle. We provide several numerical examples to illustrate the eﬃcacy of our method. The method not only delivers precise solutions but also showcases remarkable computational eﬃciency.more » « less

Abstract Fixedpoint fast sweeping WENO methods are a class of efficient highorder numerical methods to solve steadystate solutions of hyperbolic partial differential equations (PDEs). The GaussSeidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steadystate solutions. A nice property of fixedpoint fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems. Hence, they are easy to be applied to a general hyperbolic system. To deal with the difficulties associated with numerical boundary treatment when highorder finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains, inverse LaxWendroff (ILW) procedures were developed as a very effective approach in the literature. In this paper, we combine a fifthorder fixedpoint fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions. Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids. Numerical results show highorder accuracy and good performance of the method. Furthermore, the method is compared with the popular thirdorder total variation diminishing RungeKutta (TVDRK3) timemarching method for steadystate computations. Numerical examples show that for most of examples, the fixedpoint fast sweeping method saves more than half CPU time costs than TVDRK3 to converge to steadystate solutions.

null (Ed.)This paper is concerned with the inverse scattering problem which aims to determine the spatially distributed dielectric constant coefficient of the 2D Helmholtz equation from multifrequency backscatter data associated with a single direction of the incident plane wave. We propose a globally convergent convexification numerical algorithm to solve this nonlinear and illposed inverse problem. The key advantage of our method over conventional optimization approaches is that it does not require a good first guess about the solution. First, we eliminate the coefficient from the Helmholtz equation using a change of variables. Next, using a truncated expansion with respect to a special Fourier basis, we approximately reformulate the inverse problem as a system of quasilinear elliptic PDEs, which can be numerically solved by a weighted quasireversibility approach. The cost functional for the weighted quasireversibility method is constructed as a Tikhonovlike functional that involves a Carleman Weight Function. Our numerical study shows that, using a version of the gradient descent method, one can find the minimizer of this Tikhonovlike functional without any advanced a priori knowledge about it.more » « less

In this paper, the convergence of a homotopy method (1.1) for solving the steady state problem of Burgers’ equation is considered. When ν is fixed, we prove that the solution of (1.1) converges to the unique steady state solution as epsilon → 0, which is independent of the initial conditions. Numerical examples are presented to confirm this conclusion by using the continuous finite element method. In contrast, when ν = epsilon → 0, numerically we show that steady state solutions obtained by (1.1) indeed depend on initial conditions.more » « less