We prove that multilevel Picard approximations are capable of approximating solutions of semilinear heat equations in $L^{p}$-sense, ${p}\in [2,\infty)$, in the case of gradient-dependent, Lipschitz-continuous nonlinearities, in the sense that the computational effort of the multilevel Picard approximations grow at most polynomially in both the dimension $d$ and the reciprocal $1/\epsilon$ of the prescribed accuracy $\epsilon$.
相關內容
We consider the discretization of a class of nonlinear parabolic equations by discontinuous Galerkin time-stepping methods and establish a priori as well as conditional a posteriori error estimates. Our approach is motivated by the error analysis in [9] for Runge-Kutta methods for nonlinear parabolic equations; in analogy to [9], the proofs are based on maximal regularity properties of discontinuous Galerkin methods for non-autonomous linear parabolic equations.
This work considers the Galerkin approximation and analysis for a hyperbolic integrodifferential equation, where the non-positive variable-sign kernel and nonlinear-nonlocal damping with both the weak and viscous damping effects are involved. We derive the long-time stability of the solution and its finite-time uniqueness. For the semi-discrete-in-space Galerkin scheme, we derive the long-time stability of the semi-discrete numerical solution and its finite-time error estimate by technical splitting of intricate terms. Then we further apply the centering difference method and the interpolating quadrature to construct a fully discrete Galerkin scheme and prove the long-time stability of the numerical solution and its finite-time error estimate by designing a new semi-norm. Numerical experiments are performed to verify the theoretical findings.
In this contribution, we address the numerical solutions of high-order asymptotic equivalent partial differential equations with the results of a lattice Boltzmann scheme for an inhomogeneous advection problem in one spatial dimension. We first derive a family of equivalent partial differential equations at various orders, and we compare the lattice Boltzmann experimental results with a spectral approximation of the differential equations. For an unsteady situation, we show that the initialization scheme at a sufficiently high order of the microscopic moments plays a crucial role to observe an asymptotic error consistent with the order of approximation. For a stationary long-time limit, we observe that the measured asymptotic error converges with a reduced order of precision compared to the one suggested by asymptotic analysis.
Finite element spaces by Whitney $k$-forms on cubical meshes in $\mathbb{R}^n$ are presented. Based on the spaces, compatible discretizations to $H\Lambda^k$ problems are provided, and discrete de Rham complexes and commutative diagrams are constructed.
A finite element method is introduced to track interface evolution governed by the level set equation. The method solves for the level set indicator function in a narrow band around the interface. An extension procedure, which is essential for a narrow band level set method, is introduced based on a finite element $L^2$- or $H^1$-projection combined with the ghost-penalty method. This procedure is formulated as a linear variational problem in a narrow band around the surface, making it computationally efficient and suitable for rigorous error analysis. The extension method is combined with a discontinuous Galerkin space discretization and a BDF time-stepping scheme. The paper analyzes the stability and accuracy of the extension procedure and evaluates the performance of the resulting narrow band finite element method for the level set equation through numerical experiments.
Stable distributions are a celebrated class of probability laws used in various fields. The $\alpha$-stable process, and its exponentially tempered counterpart, the Classical Tempered Stable (CTS) process, are also prominent examples of L\'evy processes. Simulating these processes is critical for many applications, yet it remains computationally challenging, due to their infinite jump activity. This survey provides an overview of the key properties of these objects offering a roadmap for practitioners. The first part is a review of the stability property, sampling algorithms are provided along with numerical illustrations. Then CTS processes are presented, with the Baeumer-Meerschaert algorithm for increment simulation, and a computational analysis is provided with numerical illustrations across different time scales.
We consider linear second order differential equation y''= f with zero Dirichlet boundary conditions. At the continuous level this problem is solvable using the Green function, and this technique has a counterpart on the discrete level. The discrete solution is represented via an application of a matrix -- the Green matrix -- to the discretised right-hand side, and we propose an algorithm for fast construction of the Green matrix. In particular, we discretise the original problem using the spectral collocation method based on the Chebyshev--Gauss--Lobatto points, and using the discrete cosine transformation we show that the corresponding Green matrix is fast to construct even for large number of collocation points/high polynomial degree. Furthermore, we show that the action of the discrete solution operator (Green matrix) to the corresponding right-hand side can be implemented in a matrix-free fashion.
We develop and analyze stochastic inexact Gauss-Newton methods for nonlinear least-squares problems and for nonlinear systems ofequations. Random models are formed using suitable sampling strategies for the matrices involved in the deterministic models. The analysis of the expected number of iterations needed in the worst case to achieve a desired level of accuracy in the first-order optimality condition provides guidelines for applying sampling and enforcing, with \minor{a} fixed probability, a suitable accuracy in the random approximations. Results of the numerical validation of the algorithms are presented.
We consider the problem of estimating the error when solving a system of differential algebraic equations. Richardson extrapolation is a classical technique that can be used to judge when computational errors are irrelevant and estimate the discretization error. We have simulated molecular dynamics with constraints using the GROMACS library and found that the output is not always amenable to Richardson extrapolation. We derive and illustrate Richardson extrapolation using a variety of numerical experiments. We identify two necessary conditions that are not always satisfied by the GROMACS library.
We study a cost-aware programming language for higher-order recursion dubbed $\textbf{PCF}_\mathsf{cost}$ in the setting of synthetic domain theory (SDT). Our main contribution relates the denotational cost semantics of $\textbf{PCF}_\mathsf{cost}$ to its computational cost semantics, a new kind of dynamic semantics for program execution that serves as a mathematically natural alternative to operational semantics in SDT. In particular we prove an internal, cost-sensitive version of Plotkin's computational adequacy theorem, giving a precise correspondence between the denotational and computational semantics for complete programs at base type. The constructions and proofs of this paper take place in the internal dependent type theory of an SDT topos extended by a phase distinction in the sense of Sterling and Harper. By controlling the interpretation of cost structure via the phase distinction in the denotational semantics, we show that $\textbf{PCF}_\mathsf{cost}$ programs also evince a noninterference property of cost and behavior. We verify the axioms of the type theory by means of a model construction based on relative sheaf models of SDT.
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191