Search results for "Numerical Analysis"
showing 10 items of 883 documents
On the use of a meshless solver for PDEs governing electromagnetic transients
2009
In this paper some key elements of the Smoothed Particle Hydrodynamics methodology suitably reformulated for analyzing electromagnetic transients are investigated. The attention is focused on the interpolating smoothing kernel function which strongly influences the computational results. Some issues are provided by adopting the polynomial reproducing conditions. Validation tests involving Gaussian and cubic B-spline smoothing kernel functions in one and two dimensions are reported.
Global sensitivity analysis in wastewater treatment modelling
2019
Global sensitivity analysis (GSA) is a valuable tool to support the use of mathematical models. GSA allows the identifcation of the effect of model and input factor uncertainty on the model response, also considering the effect due to the interactions among factors. During recent years, the wastewater modelling feld has embraced the use of GSA. Wastewater modellers have tried to transfer the knowledge and experience from other disciplines and other water modelling felds.
Representation of capacity drop at a road merge via point constraints in a first order traffic model
2018
We reproduce the capacity drop phenomenon at a road merge by implementing a non-local point constraint at the junction in a first order traffic model. We call capacity drop the situation in which the outflow through the junction is lower than the receiving capacity of the outgoing road, as too many vehicles trying to access the junction from the incoming roads hinder each other. In this paper, we first construct an enhanced version of the locally constrained model introduced by Haut et al. (Proceedings 16th IFAC World Congress. Prague, Czech Republic 229 (2005) TuM01TP/3), then we propose its counterpart featuring a non-local constraint and finally we compare numerically the two models by c…
Robust and Efficient IMEX Schemes for Option Pricing under Jump-Diffusion Models
2013
We propose families of IMEX time discretization schemes for the partial integro-differential equation derived for the pricing of options under a jump diffusion process. The schemes include the families of IMEX-midpoint, IMEXCNAB and IMEX-BDF2 schemes. Each family is defined by a convex parameter c ∈ [0, 1], which divides the zeroth-order term due to the jumps between the implicit and explicit part in the time discretization. These IMEX schemes lead to tridiagonal systems, which can be solved extremely efficiently. The schemes are studied through Fourier stability analysis and numerical experiments. It is found that, under suitable assumptions and time step restrictions, the IMEX-midpoint fa…
Solving the pentahedron problem
2015
Nowadays, all geometric modelers provide some tools for specifying geometric constraints. The 3D pentahedron problem is an example of a 3D Geometric Constraint Solving Problem (GCSP), composed of six vertices, nine edges, five faces (two triangles and three quadrilaterals), and defined by the lengths of its edges and the planarity of its quadrilateral faces. This problem seems to be the simplest non-trivial problem, as the methods used to solve the Stewart platform or octahedron problem fail to solve it. The naive algebraic formulation of the pentahedron yields an under-constrained system of twelve equations in eighteen unknowns. Even if the use of placement rules transforms the pentahedron…
An IMEX-Scheme for Pricing Options under Stochastic Volatility Models with Jumps
2014
Partial integro-differential equation (PIDE) formulations are often preferable for pricing options under models with stochastic volatility and jumps, especially for American-style option contracts. We consider the pricing of options under such models, namely the Bates model and the so-called stochastic volatility with contemporaneous jumps (SVCJ) model. The nonlocality of the jump terms in these models leads to matrices with full matrix blocks. Standard discretization methods are not viable directly since they would require the inversion of such a matrix. Instead, we adopt a two-step implicit-explicit (IMEX) time discretization scheme, the IMEX-CNAB scheme, where the jump term is treated ex…
A rigidity problem on the round sphere
2015
We consider a class of overdetermined problems in rotationally symmetric spaces, which reduce to the classical Serrin's overdetermined problem in the case of the Euclidean space. We prove some general integral identities for rotationally symmetric spaces which imply a rigidity result in the case of the round sphere.
The Calderón problem with partial data on manifolds and applications
2013
We consider Calderon's inverse problem with partial data in dimensions $n \geq 3$. If the inaccessible part of the boundary satisfies a (conformal) flatness condition in one direction, we show that this problem reduces to the invertibility of a broken geodesic ray transform. In Euclidean space, sets satisfying the flatness condition include parts of cylindrical sets, conical sets, and surfaces of revolution. We prove local uniqueness in the Calderon problem with partial data in admissible geometries, and global uniqueness under an additional concavity assumption. This work unifies two earlier approaches to this problem (\cite{KSU} and \cite{I}) and extends both. The proofs are based on impr…
Partial data inverse problems for the Hodge Laplacian
2017
We prove uniqueness results for a Calderon type inverse problem for the Hodge Laplacian acting on graded forms on certain manifolds in three dimensions. In particular, we show that partial measurements of the relative-to-absolute or absolute-to-relative boundary value maps uniquely determine a zeroth order potential. The method is based on Carleman estimates for the Hodge Laplacian with relative or absolute boundary conditions, and on the construction of complex geometric optics solutions which reduce the Calderon type problem to a tensor tomography problem for 2-tensors. The arguments in this paper allow to establish partial data results for elliptic systems that generalize the scalar resu…
Nonlinear hyperbolic equations in surface theory: integrable discretizations and approximation results
2006
A numerical scheme is developed for solution of the Goursat problem for a class of nonlinear hyperbolic systems with an arbitrary number of independent variables. Convergence results are proved for this difference scheme. These results are applied to hyperbolic systems of differential-geometric origin, like the sine-Gordon equation describing the surfaces of the constant negative Gaussian curvature (K-surfaces). In particular, we prove the convergence of discrete K--surfaces and their Backlund transformations to their continuous counterparts. This puts on a firm basis the generally accepted belief (which however remained unproved untill this work) that the classical differential geometry of…