skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An endeavor from the Glowinski-Le Tallec splitting for approximating the solution of Kawarada equation
This paper studies an extended application of the Glowinski-Le Tallec splitting for approximating solutions of linear and nonlinear partial differential equations. It is shown that the three-level, six-component operator decomposition, originally designed for Lagrangian optimizations, provides a stable second-order operator splitting approximation for the solutions of evolutional partial differential equations. It is also found that the Glowinski-Le Tallec formula not only provides an effective enhancement to conventional two-level, four-component ADI and LOD methods, but also introduces a flexible way for constructing multi-parameter operator splitting strategies in respective spaces where broad spectrums of mathematical models may exist for important natural phenomena and applications. The extended operator splitting is utilized for solving a singular and nonlinear Kawarada problem satisfactorily. Multiple simulation results are presented.  more » « less
Award ID(s):
2318032
PAR ID:
10505723
Author(s) / Creator(s):
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Journal of Mathematical Analysis and Applications
Volume:
534
Issue:
1
ISSN:
0022-247X
Page Range / eLocation ID:
128051
Subject(s) / Keyword(s):
Operator splitting, Approximation order, Kawarada equations, Quenching blow-up, Finite difference methods, Stability
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Lump solutions are analytical rational function solutions localized in all directions in space. We analyze a class of lump solutions, generated from quadratic functions, to nonlinear partial differential equations. The basis of success is the Hirota bilinear formulation and the primary object is the class of positive multivariate quadratic functions. A complete determination of quadratic functions positive in space and time is given, and positive quadratic functions are characterized as sums of squares of linear functions. Necessary and sufficient conditions for positive quadratic functions to solve Hirota bilinear equations are presented, and such polynomial solutions yield lump solutions to nonlinear partial differential equations under the dependent variable logarithmic derivative transformations. Applications are made for a few generalized KP and BKP equations. 
    more » « less
  2. null (Ed.)
    Abstract Combining the classical theory of optimal transport with modern operator splitting techniques, we develop a new numerical method for nonlinear, nonlocal partial differential equations, arising in models of porous media, materials science, and biological swarming. Our method proceeds as follows: first, we discretize in time, either via the classical JKO scheme or via a novel Crank–Nicolson-type method we introduce. Next, we use the Benamou–Brenier dynamical characterization of the Wasserstein distance to reduce computing the solution of the discrete time equations to solving fully discrete minimization problems, with strictly convex objective functions and linear constraints. Third, we compute the minimizers by applying a recently introduced, provably convergent primal dual splitting scheme for three operators (Yan in J Sci Comput 1–20, 2018). By leveraging the PDEs’ underlying variational structure, our method overcomes stability issues present in previous numerical work built on explicit time discretizations, which suffer due to the equations’ strong nonlinearities and degeneracies. Our method is also naturally positivity and mass preserving and, in the case of the JKO scheme, energy decreasing. We prove that minimizers of the fully discrete problem converge to minimizers of the spatially continuous, discrete time problem as the spatial discretization is refined. We conclude with simulations of nonlinear PDEs and Wasserstein geodesics in one and two dimensions that illustrate the key properties of our approach, including higher-order convergence our novel Crank–Nicolson-type method, when compared to the classical JKO method. 
    more » « less
  3. A new bifurcation from simple eigenvalue theorem is proved for general nonlinear functional equations. It is shown that in this bifurcation scenario, the bifurcating solutions are on a curve which is tangent to the line of trivial solutions, while in typical bifurcations the curve of bifurcating solutions is transversal to the line of trivial ones. The stability of bifurcating solutions can be determined, and examples from partial differential equations are shown to demonstrate such bifurcations. 
    more » « less
  4. null (Ed.)
    In this paper, we propose a local discontinuous Galerkin (LDG) method for nonlinear and possibly degenerate parabolic stochastic partial differential equations, which is a high-order numerical scheme. It extends the discontinuous Galerkin (DG) method for purely hyperbolic equations to parabolic equations and shares with the DG method its advantage and flexibility. We prove the L 2 -stability of the numerical scheme for fully nonlinear equations. Optimal error estimates ( O ( h (k+1) )) for smooth solutions of semi-linear stochastic equations is shown if polynomials of degree k are used. We use an explicit derivative-free order 1.5 time discretization scheme to solve the matrix-valued stochastic ordinary differential equations derived from the spatial discretization. Numerical examples are given to display the performance of the LDG method. 
    more » « less
  5. This paper is concerned with the problem of counting solutions of stationary nonlinear Partial Differential Equations (PDEs) when the PDE is known to admit more than one solution. We suggest tackling the problem via a sampling-based approach. The method allows one to find solutions that are stable, in the sense that they are stable equilibria of the associated time-dependent PDE. We test our proposed methodology on the McKean–Vlasov PDE, more precisely on the problem of determining the number of stationary solutions of the McKean–Vlasov equation. This article is part of the theme issue ‘Partial differential equations in data science’. 
    more » « less