Search results for " numerical analysis."
showing 10 items of 103 documents
Monotone cubic spline interpolation for functions with a strong gradient
2021
Abstract Spline interpolation has been used in several applications due to its favorable properties regarding smoothness and accuracy of the interpolant. However, when there exists a discontinuity or a steep gradient in the data, some artifacts can appear due to the Gibbs phenomenon. Also, preservation of data monotonicity is a requirement in some applications, and that property is not automatically verified by the interpolator. Hence, some additional techniques have to be incorporated so as to ensure monotonicity. The final interpolator is not actually a spline as C 2 regularity and monotonicity are not ensured at the same time. In this paper, we study sufficient conditions to obtain monot…
Numerical Modelling for Assessing the Structural Efficiency of CAM® Reinforcement System for Masonry Walls
2022
A large portion of the Italian building heritage is made of masonry construc-tions, which were erected in the first decades of the last century and were conceived to support gravitational loads only. Many research activities have been carried out in order to propose retrofitting techniques aimed at improving seismic behaviour of existing ma-sonry buildings. The present paper focuses on the CAM® reinforcing system. Such a retrofitting technique, which is widely used in Italy, consists in the application of pre-tensioned stainless steel ribbons on existing masonry walls, conferring to them addi-tional strength, ductility and beneficial confinement effects. The preliminary outcomes of a common…
Wave Propagation in a 3-D Optical Waveguide
2003
In this paper we study the problem of wave propagation in a 3-D optical fiber. The goal is to obtain a solution for the time-harmonic field caused by a source in a cylindrically symmetric waveguide. The geometry of the problem, corresponding to an open waveguide, makes the problem challenging. To solve it, we construct a transform theory which is a nontrivial generalization of a method for solving a 2-D version of this problem given by Magnanini and Santosa.\cite{MS} The extension to 3-D is made complicated by the fact that the resulting eigenvalue problem defining the transform kernel is singular both at the origin and at infinity. The singularities require the investigation of the behavio…
Guaranteed lower bounds for cost functionals of time-periodic parabolic optimization problems
2019
In this paper, a new technique is shown for deriving computable, guaranteed lower bounds of functional type (minorants) for two different cost functionals subject to a parabolic time-periodic boundary value problem. Together with previous results on upper bounds (majorants) for one of the cost functionals, both minorants and majorants lead to two-sided estimates of functional type for the optimal control problem. Both upper and lower bounds are derived for the second new cost functional subject to the same parabolic PDE-constraints, but where the target is a desired gradient. The time-periodic optimal control problems are discretized by the multiharmonic finite element method leading to lar…
Numerical Analysis of Bearing Capacity of a Ring Footing on Geogrid Reinforced Sand
2021
A ring footing is found to be of practical importance in supporting symmetrical constructions for example silos, oil storage container etc. In the present paper, numerical analysis was carried out with explicit code FLAC3D 7.0 to investigate bearing capacity of a ring footing on geogrid reinforced sand. Effects of the ratio n of its inner/outer diameter (Di/D) of a ring footing, an optimum depth to lay the geogrid layer were examined. It was found that an intersection zone was developed in soil under inner-side (aisle) of ring footing, contributing to its bearing capacity. Substantial increase of bearing capacities could be realized if ratio n of a ring footing was around 0.6. Numerical res…
On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method
2017
The Scheduled Relaxation Jacobi (SRJ) method is an extension of the classical Jacobi iterative method to solve linear systems of equations ($Au=b$) associated with elliptic problems. It inherits its robustness and accelerates its convergence rate computing a set of $P$ relaxation factors that result from a minimization problem. In a typical SRJ scheme, the former set of factors is employed in cycles of $M$ consecutive iterations until a prescribed tolerance is reached. We present the analytic form for the optimal set of relaxation factors for the case in which all of them are different, and find that the resulting algorithm is equivalent to a non-stationary generalized Richardson's method. …
Multi-domain spectral approach with Sommerfeld condition for the Maxwell equations
2021
We present a multidomain spectral approach with an exterior compactified domain for the Maxwell equations for monochromatic fields. The Sommerfeld radiation condition is imposed exactly at infinity being a finite point on the numerical grid. As an example, axisymmetric situations in spherical and prolate spheroidal coordinates are discussed.
Scheduled Relaxation Jacobi method: improvements and applications
2016
Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficien…
A partially reflecting random walk on spheres algorithm for electrical impedance tomography
2015
In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias…
Efficient numerical integration of neutrino oscillations in matter
2016
A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a general numerical integrator, a fact that could smooth the way for massive numerical integration concomitant with experimental data analyses. Detailed illustrations about numerical procedure and computer time costs are provided.