Search results for "Numerical Analysis"
showing 10 items of 883 documents
Numerical Algorithms Based on Characteristic Domain Decomposition for Obstacle Problems
1997
A new numerical solution algorithm for obstacle problems is proposed, where the characteristic domain decomposition into active and inactive subdomains separated by the free boundary is approximated by a Schwarz method. Such an approach gives an opportunity to apply fast linear system solvers to genuinely non-linear obstacle problems. Other solution algorithms, like projected relaxation methods and active set strategies, are compared to the new solution algorithm. Numerical experiments related to the elastoplastic torsion problem are included showing the efficiency of the new approach.
Computation of a few smallest eigenvalues of elliptic operators using fast elliptic solvers
2001
The computation of a few smallest eigenvalues of generalized algebraic eigenvalue problems is studied. The considered problems are obtained by discretizing self-adjoint second-order elliptic partial differential eigenvalue problems in two- or three-dimensional domains. The standard Lanczos algorithm with the complete orthogonalization is used to compute some eigenvalues of the inverted eigenvalue problem. Under suitable assumptions, the number of Lanczos iterations is shown to be independent of the problem size. The arising linear problems are solved using some standard fast elliptic solver. Numerical experiments demonstrate that the inverted problem is much easier to solve with the Lanczos…
The Calderón problem for the fractional Schrödinger equation
2020
We show global uniqueness in an inverse problem for the fractional Schr\"odinger equation: an unknown potential in a bounded domain is uniquely determined by exterior measurements of solutions. We also show global uniqueness in the partial data problem where the measurements are taken in arbitrary open, possibly disjoint, subsets of the exterior. The results apply in any dimension $\geq 2$ and are based on a strong approximation property of the fractional equation that extends earlier work. This special feature of the nonlocal equation renders the analysis of related inverse problems radically different from the traditional Calder\'on problem.
Comparative Study of the a Posteriori Error Estimators for the Stokes Problem
2007
The research presented is focused on a comparative study of a posteriori error estimation methods to various approximations of the Stokes problem. Mainly, we are interested in the performance of functional type a posterior error estimates and their comparison with other methods. We show that functional type a posteriori error estimators are applicable to various types of approximations (including non-Galerkin ones) and robust with respect to the mesh structure, type of the finite element and computational procedure used. This allows the construction of effective mesh adaptation procedures in all cases considered. Numerical tests justify the approach suggested.
The Exponential Dichotomy under Discretization on General Approximation Scheme
2011
This paper is devoted to the numerical analysis of abstract parabolic problem 𝑢 ( 𝑡 ) = 𝐴 𝑢 ( 𝑡 ) ; 𝑢 ( 0 ) = 𝑢 0 , with hyperbolic generator 𝐴 . We are developing a general approach to establish a discrete dichotomy in a very general setting in case of discrete approximation in space and time. It is a well-known fact that the phase space in the neighborhood of the hyperbolic equilibrium can be split in a such way that the original initial value problem is reduced to initial value problems with exponential decaying solutions in opposite time direction. We use the theory of compact approximation principle and collectively condensing approximation to show that such a decomposition o…
Spritz: General relativistic magnetohydrodynamics with neutrinos
2020
We here present a new version of the publicly available general relativistic magnetohydrodynamic (GRMHD) code $\texttt{Spritz}$, which now includes an approximate neutrino leakage scheme able to handle neutrino cooling and heating. The leakage scheme is based on the publicly available $\texttt{ZelmaniLeak}$ code, with a few modifications in order to properly work with $\texttt{Spritz}$. We discuss the involved equations, physical assumptions, and implemented numerical methods, along with a large battery of general relativistic tests performed with and without magnetic fields. Our tests demonstrate the correct implementation of the neutrino leakage scheme, paving the way for further improvem…
Asymptotic stability of solutions to Volterra-renewal integral equations with space maps
2012
Abstract In this paper we consider linear Volterra-renewal integral equations (VIEs) whose solutions depend on a space variable, via a map transformation. We investigate the asymptotic properties of the solutions, and study the asymptotic stability of a numerical method based on direct quadrature in time and interpolation in space. We show its properties through test examples.
A Dirichlet problem for the Laplace operator in a domain with a small hole close to the boundary
2016
We study the Dirichlet problem in a domain with a small hole close to the boundary. To do so, for each pair $\boldsymbol\varepsilon = (\varepsilon_1, \varepsilon_2 )$ of positive parameters, we consider a perforated domain $\Omega_{\boldsymbol\varepsilon}$ obtained by making a small hole of size $\varepsilon_1 \varepsilon_2 $ in an open regular subset $\Omega$ of $\mathbb{R}^n$ at distance $\varepsilon_1$ from the boundary $\partial\Omega$. As $\varepsilon_1 \to 0$, the perforation shrinks to a point and, at the same time, approaches the boundary. When $\boldsymbol\varepsilon \to (0,0)$, the size of the hole shrinks at a faster rate than its approach to the boundary. We denote by $u_{\bolds…
The Effects of Orography on the Extratropical Transition of Tropical Cyclones: A Case Study of Typhoon Sinlaku (2008)
2018
Abstract Extratropical transition (ET) can cause high-impact weather in midlatitude regions and therefore constitutes an ongoing threat at the end of a tropical cyclone’s (TC) life cycle. Most of the ET events occur over the ocean, but some TCs recurve and undergo ET along coastal regions; however, the latter category is less investigated. Typhoon Sinlaku (2008), for example, underwent ET along the southern coast of Japan. It was one of the typhoons that occurred during the T-PARC field campaign, providing unprecedented high-resolution observational data. Sinlaku is therefore an excellent case to investigate the impact of a coastal region, and in particular orography, on the evolution of ET…
Numerical simulation of internal boundary-layer development and comparison with atmospheric data
2006
A finite-volume numerical model is employed to investigate the adaptation of the atmospheric boundary layer to a change in the underlying surface roughness, such as that existing in the transition from land to the free surface of a water body. Numerical results are validated by comparison with neutral stratification atmospheric data and compared with the internal boundary-layer (IBL) heights computed using a number of existing empirical formulae. The numerical analysis allows an extension of the fetch range in which the existing formulae, calibrated only by comparison with short fetch data, may be applied. An argument is offered that the spatial variability of the water surface roughness sh…