Search results for "Numerical Analysis"
showing 10 items of 883 documents
Capturing shock waves in inelastic granular gases
2005
Shock waves in granular gases generated by hitting an obstacle at rest are treated by means of a shock capturing scheme that approximates the Euler equations of granular gas dynamics with an equation of state (EOS), introduced by Goldshtein and Shapiro [J. Fluid Mech. 282 (1995) 75-114], that takes into account the inelastic collisions of granules. We include a sink term in the energy balance to account for dissipation of the granular motion by collisional inelasticity, proposed by Haff [J. Fluid Mech. 134 (1983) 401-430], and the gravity field added as source terms. We have computed the approximate solution to a one-dimensional granular gas falling on a plate under the acceleration of grav…
A flux-split algorithm applied to conservative models for multicomponent compressible flows
2003
In this paper we consider a conservative extension of the Euler equations for gas dynamics to describe a two-component compressible flow in Cartesian coordinates. It is well known that classical shock-capturing schemes applied to conservative models are oscillatory near the interface between the two gases. Several authors have addressed this problem proposing either a primitive consistent algorithm [J. Comput. Phys. 112 (1994) 31] or Lagrangian ingredients (Ghost Fluid Method by Fedkiw et al. [J. Comput. Phys. 152 (1999) 452] and [J. Comput. Phys. 169 (2001) 594]). We solve directly this conservative model by a flux-split algorithm, due to the first author (see [J. Comput. Phys. 125 (1996) …
Random wave run-up with a physically-based Lagrangian shoreline model
2014
Abstract In the present paper the run-up of random waves was calculated by means of a numerical method. In situ measurements based on a video imaging technique have been used for the validation of the present numerical model. The on-site run-up measurements have been carried out at Lido Signorino beach, near Marsala, Italy,along a transect, normal to the shore. A video camera and a linear array of rods have been used to obtain field data. Numerical simulations with a 1DH Boussinesq-type of model for breaking waves which takes into account the wave run-up by means of a Lagrangian shoreline model have been carried out. In such simulations random waves of given spectrum have been propagated in…
Iterative Reconstruction of Signals on Graph
2020
We propose an iterative algorithm to interpolate graph signals from only a partial set of samples. Our method is derived from the well known Papoulis-Gerchberg algorithm by considering the optimal value of a constant involved in the iteration step. Compared with existing graph signal reconstruction algorithms, the proposed method achieves similar or better performance both in terms of convergence rate and computational efficiency.
Numerical approach for signal delay in general distributed networks
2003
The authors consider a general network with telegraph equations modelling distributed elements and having, additionally, nonlinear capacitors. A global asymptotic exponential stability of the solution is given. A simple computable upper bound of the delay time is given. Numerical examples illustrate the usefulness of the results. >
Cell-average WENO with progressive order of accuracy close to discontinuities with applications to signal processing
2020
In this paper we translate to the cell-average setting the algorithm for the point-value discretization presented in S. Amat, J. Ruiz, C.-W. Shu, D. F. Y\'a\~nez, A new WENO-2r algorithm with progressive order of accuracy close to discontinuities, submitted to SIAM J. Numer. Anal.. This new strategy tries to improve the results of WENO-($2r-1$) algorithm close to the singularities, resulting in an optimal order of accuracy at these zones. The main idea is to modify the optimal weights so that they have a nonlinear expression that depends on the position of the discontinuities. In this paper we study the application of the new algorithm to signal processing using Harten's multiresolution. Se…
Cubic Local Splines on Non-uniform Grid
2015
In this chapter, two types of local cubic splines on non-uniform grids are described: 1. The simplest variation-diminishing splines and 2. The quasi-interpolating splines. The splines are computed by a simple fast computational algorithms that utilizes a relation between the splines and cubic interpolation polynomials. Those splines can serve as an efficient tool for real-time signal processing. As an input, they use either clean or noised arbitrarily-spaced samples. On the other hand, the capability to adapt the grid to the structure of an object and minimal requirements to the operating memory are great advantages for off-line processing of signals and multidimensional data arrays.
Spectrum cartography using adaptive radial basis functions: Experimental validation
2017
In this paper, we experimentally validate the functionality of a developed algorithm for spectrum cartography using adaptive Gaussian radial basis functions (RBF). The RBF are strategically centered around representative centroid locations in a machine learning context. We assume no prior knowledge about neither the power spectral densities (PSD) of the transmitters nor their locations. Instead, the received signal power at each location is estimated as a linear combination of different RBFs. The weights of the RBFs, their Gaussian decaying parameters and locations are jointly optimized using expectation maximization with a least squares loss function and a quadratic regularizer. The perfor…
Multidomain spectral method for the Gauss hypergeometric function
2018
International audience; We present a multidomain spectral approach for Fuchsian ordinary differential equations in the particular case of the hypergeometric equation. Our hybrid approach uses Frobenius’ method and Moebius transformations in the vicinity of each of the singular points of the hypergeometric equation, which leads to a natural decomposition of the real axis into domains. In each domain, solutions to the hypergeometric equation are constructed via the well-conditioned ultraspherical spectral method. The solutions are matched at the domain boundaries to lead to a solution which is analytic on the whole compactified real line R∪∞, except for the singular points and cuts of the Rie…
Numerical study of soliton stability, resolution and interactions in the 3D Zakharov–Kuznetsov equation
2021
International audience; We present a detailed numerical study of solutions to the Zakharov-Kuznetsov equation in three spatial dimensions. The equation is a three-dimensional generalization of the Korteweg-de Vries equation, though, not completely integrable. This equation is L-2-subcritical, and thus, solutions exist globally, for example, in the H-1 energy space.We first study stability of solitons with various perturbations in sizes and symmetry, and show asymptotic stability and formation of radiation, confirming the asymptotic stability result in Farah et al. (0000) for a larger class of initial data. We then investigate the solution behavior for different localizations and rates of de…