Search results for " numerical analysis."
showing 10 items of 103 documents
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…
A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries
2021
Author(s): Gulizzi, Vincenzo; Almgren, Ann S; Bell, John B | Abstract: We present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an $hp$ adaptive mesh refinement ($hp$-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The $hp$-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a…
On a nonlinear Schrödinger equation for nucleons in one space dimension
2021
We study a 1D nonlinear Schrödinger equation appearing in the description of a particle inside an atomic nucleus. For various nonlinearities, the ground states are discussed and given in explicit form. Their stability is studied numerically via the time evolution of perturbed ground states. In the time evolution of general localized initial data, they are shown to appear in the long time behaviour of certain cases.
Functional a posteriori error estimates for boundary element methods
2019
Functional error estimates are well-established tools for a posteriori error estimation and related adaptive mesh-refinement for the finite element method (FEM). The present work proposes a first functional error estimate for the boundary element method (BEM). One key feature is that the derived error estimates are independent of the BEM discretization and provide guaranteed lower and upper bounds for the unknown error. In particular, our analysis covers Galerkin BEM and the collocation method, what makes the approach of particular interest for scientific computations and engineering applications. Numerical experiments for the Laplace problem confirm the theoretical results.