Search results for " Numerical analysis"
showing 10 items of 106 documents
On the a posteriori error analysis for linear Fokker-Planck models in convection-dominated diffusion problems
2018
This work is aimed at the derivation of reliable and efficient a posteriori error estimates for convection-dominated diffusion problems motivated by a linear Fokker-Planck problem appearing in computational neuroscience. We obtain computable error bounds of the functional type for the static and time-dependent case and for different boundary conditions (mixed and pure Neumann boundary conditions). Finally, we present a set of various numerical examples including discussions on mesh adaptivity and space-time discretisation. The numerical results confirm the reliability and efficiency of the error estimates derived.
Spectral approach to D-bar problems
2017
We present the first numerical approach to D-bar problems having spectral convergence for real analytic, rapidly decreasing potentials. The proposed method starts from a formulation of the problem in terms of an integral equation that is numerically solved with Fourier techniques. The singular integrand is regularized analytically. The resulting integral equation is approximated via a discrete system that is solved with Krylov methods. As an example, the D-bar problem for the Davey-Stewartson II equations is considered. The result is used to test direct numerical solutions of the PDE.© 2017 Wiley Periodicals, Inc.
Reliable numerical solution of a class of nonlinear elliptic problems generated by the Poisson-Boltzmann equation
2020
We consider a class of nonlinear elliptic problems associated with models in biophysics, which are described by the Poisson-Boltzmann equation (PBE). We prove mathematical correctness of the problem, study a suitable class of approximations, and deduce guaranteed and fully computable bounds of approximation errors. The latter goal is achieved by means of the approach suggested in [S. Repin, A posteriori error estimation for variational problems with uniformly convex functionals. Math. Comp., 69:481-500, 2000] for convex variational problems. Moreover, we establish the error identity, which defines the error measure natural for the considered class of problems and show that it yields computa…
Optimal recovery of a radiating source with multiple frequencies along one line
2020
We study an inverse problem where an unknown radiating source is observed with collimated detectors along a single line and the medium has a known attenuation. The research is motivated by applications in SPECT and beam hardening. If measurements are carried out with frequencies ranging in an open set, we show that the source density is uniquely determined by these measurements up to averaging over levelsets of the integrated attenuation. This leads to a generalized Laplace transform. We also discuss some numerical approaches and demonstrate the results with several examples.
Experimental analysis of "bovedas tabicadas"
2012
A particular type of vaulted structures, known as bóvedas tabicadas (Catalan vaults), made with alternate layers of bricks and mortar and characterized by a very low thickness in comparison to the other two dimensions, are studied. The research has been faced under the historical aspect as well as under the mechanical one, developing experimental and numerical analysis regarding real existing structures. In particular, the experimental analysis, performed by effecting compression and bending tests on samples taken from a real structure, has been devoted to the constitutive material identification and the characterization is made by considering the material as an ideal homogeneous; further e…
Energy-stable linear schemes for polymer-solvent phase field models
2017
We present new linear energy-stable numerical schemes for numerical simulation of complex polymer-solvent mixtures. The mathematical model proposed by Zhou, Zhang and E (Physical Review E 73, 2006) consists of the Cahn-Hilliard equation which describes dynamics of the interface that separates polymer and solvent and the Oldroyd-B equations for the hydrodynamics of polymeric mixtures. The model is thermodynamically consistent and dissipates free energy. Our main goal in this paper is to derive numerical schemes for the polymer-solvent mixture model that are energy dissipative and efficient in time. To this end we will propose several problem-suited time discretizations yielding linear scheme…
New degrees of freedom for differential forms on cubical meshes
2022
We consider new degrees of freedom for higher order differential forms on cubical meshes. The approach is inspired by the idea of Rapetti and Bossavit to define higher order Whitney forms and their degrees of freedom using small simplices. We show that higher order differential forms on cubical meshes can be defined analogously using small cubes and prove that these small cubes yield unisolvent degrees of freedom. Significantly, this approach is compatible with discrete exterior calculus and expands the framework to cover higher order methods on cubical meshes, complementing the earlier strategy based on simplices.
Convergence of dynamic programming principles for the $p$-Laplacian
2018
We provide a unified strategy to show that solutions of dynamic programming principles associated to the $p$-Laplacian converge to the solution of the corresponding Dirichlet problem. Our approach includes all previously known cases for continuous and discrete dynamic programming principles, provides new results, and gives a convergence proof free of probability arguments.
Inverse problems and invisibility cloaking for FEM models and resistor networks
2013
In this paper we consider inverse problems for resistor networks and for models obtained via the finite element method (FEM) for the conductivity equation. These correspond to discrete versions of the inverse conductivity problem of Calderón. We characterize FEM models corresponding to a given triangulation of the domain that are equivalent to certain resistor networks, and apply the results to study nonuniqueness of the discrete inverse problem. It turns out that the degree of nonuniqueness for the discrete problem is larger than the one for the partial differential equation. We also study invisibility cloaking for FEM models, and show how an arbitrary body can be surrounded with a layer …
A fast Fourier transform based direct solver for the Helmholtz problem
2018
This article is devoted to the efficient numerical solution of the Helmholtz equation in a two‐ or three‐dimensional (2D or 3D) rectangular domain with an absorbing boundary condition (ABC). The Helmholtz problem is discretized by standard bilinear and trilinear finite elements on an orthogonal mesh yielding a separable system of linear equations. The main key to high performance is to employ the fast Fourier transform (FFT) within a fast direct solver to solve the large separable systems. The computational complexity of the proposed FFT‐based direct solver is O(N log N) operations. Numerical results for both 2D and 3D problems are presented confirming the efficiency of the method discussed…