Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDiscrete transparent boundary conditions for the two-dimensional leap-frog scheme
80
0
0.0
(
0
)
Added by
Christophe Besse
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Authors
Christophe Besse
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
In this paper, we consider artificial boundary conditions for the linearized mixed Korteweg-de Vries (KDV) Benjamin-Bona-Mahoney (BBM) equation which models water waves in the small amplitude, large wavelength regime. Continuous (respectively discrete) artificial boundary conditions involve non local operators in time which in turn requires to compute time convolutions and invert the Laplace transform of an analytic function (respectively the Z-transform of an holomor-phic function). In this paper, we propose a new, stable and fairly general strategy to carry out this crucial step in the design of transparent boundary conditions. For large time simulations, we also introduce a methodology based on the asymptotic expansion of coefficients involved in exact direct transparent boundary conditions. We illustrate the accuracy of our methods for Gaussian and wave packets initial data.
We show that the Strang splitting method applied to a diffusion-reaction equation with inhomogeneous general oblique boundary conditions is of order two when the diffusion equation is solved with the Crank-Nicolson method, while order reduction occurs in general if using other Runge-Kutta schemes or even the exact flow itself for the diffusion part. We prove these results when the source term only depends on the space variable, an assumption which makes the splitting scheme equivalent to the Crank-Nicolson method itself applied to the whole problem. Numerical experiments suggest that the second order convergence persists with general nonlinearities.
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].
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.
In this paper, we propose a strategy to determine the Dirichlet-to-Neumann (DtN) operator for infinite, lossy and locally perturbed hexagonal periodic media. We obtain a factorization of this operator involving two non local operators. The first one is a DtN type operator and corresponds to a half-space problem. The second one is a Dirichlet-to-Dirichlet (DtD) type operator related to the symmetry properties of the problem. The half-space DtN operator is characterized via Floquet-Bloch transform, a family of elementary strip problems and a family of stationary Riccati equations. The DtD operator is the solution of an affine operator valued equation which can be reformulated as a non standard integral equation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
