2,414 research outputs found
On stability in the maximum norm of difference scheme for nonlinear parabolic equation with nonlocal condition
We construct and analyze the backward Euler method for one nonlinear one-dimensional parabolic equation with nonlocal boundary condition. The main objective of this article is to investigate the stability and convergence of the difference scheme in the maximum norm. For this purpose, we use the M-matrices theory. We describe some new approach for the estimation of the error of solution and construct the majorant for it. Some conclusions and discussion of our approach are presented
Probabilistic foundation of nonlocal diffusion and formulation and analysis for elliptic problems on uncertain domains
2011 Summer.Includes bibliographical references.In the first part of this dissertation, we study the nonlocal diffusion equation with so-called Lévy measure ν as the master equation for a pure-jump Lévy process. In the case ν ∈ L1(R), a relationship to fractional diffusion is established in a limit of vanishing nonlocality, which implies the convergence of a compound Poisson process to a stable process. In the case ν ∉ L1(R), the smoothing of the nonlocal operator is shown to correspond precisely to the activity of the underlying Lévy process and the variation of its sample paths. We introduce volume-constrained nonlocal diffusion equations and demonstrate that they are the master equations for Lévy processes restricted to a bounded domain. The ensuing variational formulation and conforming finite element method provide a powerful tool for studying both Lévy processes and fractional diffusion on bounded, non-simple geometries with volume constraints. In the second part of this dissertation, we consider the problem of estimating the distribution of a quantity of interest computed from the solution of an elliptic partial differential equation posed on a domain Ω(θ) ⊂ R2 with a randomly perturbed boundary, where (θ) is a random vector with given probability structure. We construct a piecewise smooth transformation from a partition of Ω(θ) to a reference domain Ω in order to avoid the complications associated with solving the problems on Ω(θ). The domain decomposition formulation is exploited by localizing the effect of the randomness to boundary elements in order to achieve a computationally efficient Monte Carlo sampling procedure. An a posteriori error analysis for the approximate distribution, which includes a deterministic error for each sample and a stochastic error from the effect of sampling, is also presented. We thus provide an efficient means to estimate the distribution of a quantity of interest via a Monte Carlo sampling procedure while also providing a posteriori error estimates for each sample
Bayesian Inverse Quantum Theory
A Bayesian approach is developed to determine quantum mechanical potentials
from empirical data. Bayesian methods, combining empirical measurements and "a
priori" information, provide flexible tools for such empirical learning
problems. The paper presents the basic theory, concentrating in particular on
measurements of particle coordinates in quantum mechanical systems at finite
temperature. The computational feasibility of the approach is demonstrated by
numerical case studies. Finally, it is shown how the approach can be
generalized to such many-body and few-body systems for which a mean field
description is appropriate. This is done by means of a Bayesian inverse
Hartree-Fock approximation.Comment: LaTex, 32 pages, 19 figure
Non-local control in the conduction coefficients: well posedness and convergence to the local limit
We consider a problem of optimal distribution of conductivities in a system
governed by a non-local diffusion law. The problem stems from applications in
optimal design and more specifically topology optimization. We propose a novel
parametrization of non-local material properties. With this parametrization the
non-local diffusion law in the limit of vanishing non-local interaction
horizons converges to the famous and ubiquitously used generalized Laplacian
with SIMP (Solid Isotropic Material with Penalization) material model. The
optimal control problem for the limiting local model is typically ill-posed and
does not attain its infimum without additional regularization. Surprisingly,
its non-local counterpart attains its global minima in many practical
situations, as we demonstrate in this work. In spite of this qualitatively
different behaviour, we are able to partially characterize the relationship
between the non-local and the local optimal control problems. We also
complement our theoretical findings with numerical examples, which illustrate
the viability of our approach to optimal design practitioners
- …