Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOptimal artificial boundary conditions based on second-order correctors for three dimensional random elliptic media
81
0
0.0
(
0
)
Added by
Lihan Wang
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We are interested in numerical algorithms for computing the electrical field generated by a charge distribution localized on scale $ell$ in an infinite heterogeneous medium, in a situation where the medium is only known in a box of diameter $Lggell$ around the support of the charge. We propose a boundary condition that with overwhelming probability is (near) optimal with respect to scaling in terms of $ell$ and $L$, in the setting where the medium is a sample from a stationary ensemble with a finite range of dependence (set to be unity and with the assumption that $ell gg 1$). The boundary condition is motivated by quantitative stochastic homogenization that allows for a multipole expansion [BGO20]. This work extends [LO21] from two to three dimensions, and thus we need to take quadrupoles, next to dipoles, into account. This in turn relies on stochastic estimates of second-order, next to first-order, correctors. These estimates are provided for finite range ensembles under consideration, based on an extension of the semi-group approach of [GO15].
rate research
Read More
We consider a free boundary problem on three-dimensional cones depending on a parameter c and study when the free boundary is allowed to pass through the vertex of the cone. Combining analysis and computer-assisted proof, we show that when c is less than 0.43, the free boundary may pass through the vertex of the cone.
We develop a general strategy in order to implement (approximate) discrete transparent boundary conditions for finite difference approximations of the two-dimensional transport equation. The computational domain is a rectangle equipped with a Cartesian grid. For the two-dimensional leapfrog scheme, we explain why our strategy provides with explicit numerical boundary conditions on the four sides of the rectangle and why it does not require prescribing any condition at the four corners of the computational domain. The stability of the numerical boundary condition on each side of the rectangle is analyzed by means of the so-called normal mode analysis. Numerical investigations for the full problem on the rectangle show that strong instabilities may occur when coupling stable strategies on each side of the rectangle. Other coupling strategies yield promising results.
In this paper we develop an existence theory for the nonlinear initial-boundary value problem with singular diffusion $partial_t u = text{div}(k(x) abla G(u))$, $u|_{t=0}=u_0$ with Neumann boundary conditions $k(x) abla G(u)cdot u = 0$. Here $xin Bsubset mathbb{R}^d$, a bounded open set with locally Lipchitz boundary, and with $ u$ as the unit outer normal. The function $G$ is Lipschitz continuous and nondecreasing, while $k(x)$ is diagonal matrix. We show that any two weak entropy solutions $u$ and $v$ satisfy $Vert{u(t)-v(t)}Vert_{L^1(B)}le Vert{u|_{t=0}-v|_{t=0}}Vert_{L^1(B)}e^{Ct}$, for almost every $tge 0$, and a constant $C=C(k,G,B)$. If we restrict to the case when the entries $k_i$ of $k$ depend only on the corresponding component, $k_i=k_i(x_i)$, we show that there exists an entropy solution, thus establishing in this case that the problem is well-posed in the sense of Hadamard.
We propose and analyze an optimal mass transport method for a random genetic drift problem driven by a Moran process under weak-selection. The continuum limit, formulated as a reaction-advection-diffusion equation known as the Kimura equation, inherits degenerate diffusion from the discrete stochastic process that conveys to the blow-up into Dirac-delta singularities hence brings great challenges to both the analytical and numerical studies. The proposed numerical method can quantitatively capture to the fullest possible extent the development of Dirac-delta singularities for genetic segregation on one hand, and preserves several sets of biologically relevant and computationally favored properties of the random genetic drift on the other. Moreover, the numerical scheme exponentially converges to the unique numerical stationary state in time at a rate independent of the mesh size up to a mesh error. Numerical evidence is given to illustrate and support these properties, and to demonstrate the spatio-temporal dynamics of random generic drift.
We present an approach to handle Dirichlet type nonlocal boundary conditions for nonlocal diffusion models with a finite range of nonlocal interactions. Our approach utilizes a linear extrapolation of prescribed boundary data. A novelty is, instead of using local gradients of the boundary data that are not available a priori, we incorporate nonlocal gradient operators into the formulation to generalize the finite differences-based methods which are pervasive in literature; our particular choice of the nonlocal gradient operators is based on the interplay between a constant kernel function and the geometry of nonlocal interaction neighborhoods. Such an approach can be potentially useful to address similar issues in peridynamics, smoothed particle hydrodynamics and other nonlocal models. We first show the well-posedness of the newly formulated nonlocal problems and then analyze their asymptotic convergence to the local limit as the nonlocality parameter shrinks to zero. We justify the second order localization rate, which is the optimal order attainable in the absence of physical boundaries.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
