Search results for "Discretization"
showing 10 items of 237 documents
The FLO Diffusive 1D-2D Model for Simulation of River Flooding
2016
An integrated 1D-2D model for the solution of the diffusive approximation of the shallow water equations, named FLO, is proposed in the present paper. Governing equations are solved using the MArching in Space and Time (MAST) approach. The 2D floodplain domain is discretized using a triangular mesh, and standard river sections are used for modeling 1D flow inside the section width occurring with low or standard discharges. 1D elements, inside the 1D domain, are quadrilaterals bounded by the trace of two consecutive sections and by the sides connecting their extreme points. The water level is assumed to vary linearly inside each quadrilateral along the flow direction, but to remain constant …
An optimization-based approach for solving a time-harmonic multiphysical wave problem with higher-order schemes
2013
This study considers developing numerical solution techniques for the computer simulations of time-harmonic fluid-structure interaction between acoustic and elastic waves. The focus is on the efficiency of an iterative solution method based on a controllability approach and spectral elements. We concentrate on the model, in which the acoustic waves in the fluid domain are modeled by using the velocity potential and the elastic waves in the structure domain are modeled by using displacement.Traditionally, the complex-valued time-harmonic equations are used for solving the time-harmonic problems. Instead of that, we focus on finding periodic solutions without solving the time-harmonic problem…
Effects of the foot evolution on the behaviour of slow-moving landslides
2011
The paper presents a time-dependent 2D numerical model which has been developed with the purpose of highlighting the effects of the slope foot evolution on the behaviour of slow-moving landslides. The model allows to quantitatively analyse how foot mass variations can influence the stability and the movement rates of the landslide. The landslide body is modelled as composed of two rigid blocks sliding on two different planes and interacting through a common boundary, which position is assumed fixed during the analysis. A finite difference approach is used to discretize the time. For each time increment, changes in model parameters are allowed, including variations in shearing resistances, g…
Computation of vertically averaged velocities in irregular sections of straight channels
2015
Abstract. Two new methods for vertically averaged velocity computation are presented, validated and compared with other available formulas. The first method derives from the well-known Huthoff algorithm, which is first shown to be dependent on the way the river cross section is discretized into several subsections. The second method assumes the vertically averaged longitudinal velocity to be a function only of the friction factor and of the so-called "local hydraulic radius", computed as the ratio between the integral of the elementary areas around a given vertical and the integral of the elementary solid boundaries around the same vertical. Both integrals are weighted with a linear shape f…
CALIBRATION OF LÉVY PROCESSES USING OPTIMAL CONTROL OF KOLMOGOROV EQUATIONS WITH PERIODIC BOUNDARY CONDITIONS
2018
We present an optimal control approach to the problem of model calibration for L\'evy processes based on a non parametric estimation procedure. The calibration problem is of considerable interest in mathematical finance and beyond. Calibration of L\'evy processes is particularly challenging as the jump distribution is given by an arbitrary L\'evy measure, which form a infinite dimensional space. In this work, we follow an approach which is related to the maximum likelihood theory of sieves. The sampling of the L\'evy process is modelled as independent observations of the stochastic process at some terminal time $T$. We use a generic spline discretization of the L\'evy jump measure and selec…
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.
Generalized wave propagation problems and discrete exterior calculus
2018
We introduce a general class of second-order boundary value problems unifying application areas such as acoustics, electromagnetism, elastodynamics, quantum mechanics, and so on, into a single framework. This also enables us to solve wave propagation problems very efficiently with a single software system. The solution method precisely follows the conservation laws in finite-dimensional systems, whereas the constitutive relations are imposed approximately. We employ discrete exterior calculus for the spatial discretization, use natural crystal structures for three-dimensional meshing, and derive a “discrete Hodge” adapted to harmonic wave. The numerical experiments indicate that the cumulat…
Numerical methods for a nonlinear impact model: A comparative study with closed-form corrections
2011
A physically based impact model-already known and exploited in the field of sound synthesis-is studied using both analytical tools and numerical simulations. It is shown that the Hamiltonian of a physical system composed of a mass impacting on a wall can be expressed analytically as a function of the mass velocity during contact. Moreover, an efficient and accurate approximation for the mass outbound velocity is presented, which allows to estimate the Hamiltonian at the end of the contact. Analytical results are then compared to numerical simulations obtained by discretizing the system with several numerical methods. It is shown that, for some regions of the parameter space, the trajectorie…
A Domain Imbedding Method with Distributed Lagrange Multipliers for Acoustic Scattering Problems
2003
The numerical computation of acoustic scattering by bounded twodimensional obstacles is considered. A domain imbedding method with Lagrange multipliers is introduced for the solution of the Helmholtz equation with a second-order absorbing boundary condition. Distributed Lagrange multipliers are used to enforce the Dirichlet boundary condition on the scatterer. The saddle-point problem arising from the conforming finite element discretization is iteratively solved by the GMRES method with a block triangular preconditioner. Numerical experiments are performed with a disc and a semi-open cavity as scatterers.
On shape differentiation of discretized electric field integral equation
2013
Abstract This work presents shape derivatives of the system matrix representing electric field integral equation discretized with Raviart–Thomas basis functions. The arising integrals are easy to compute with similar methods as the entries of the original system matrix. The results are compared to derivatives computed with automatic differentiation technique and finite differences, and are found to be in an excellent agreement. Furthermore, the derived formulas are employed to analyze shape sensitivity of the input impedance of a planar inverted F-antenna, and the results are compared to those obtained using a finite difference approximation.