Search results for "Numerical"
showing 10 items of 2002 documents
Numerical study of the Kadomtsev–Petviashvili equation and dispersive shock waves
2018
A detailed numerical study of the long time behaviour of dispersive shock waves in solutions to the Kadomtsev-Petviashvili (KP) I equation is presented. It is shown that modulated lump solutions emerge from the dispersive shock waves. For the description of dispersive shock waves, Whitham modulation equations for KP are obtained. It is shown that the modulation equations near the soliton line are hyperbolic for the KPII equation while they are elliptic for the KPI equation leading to a focusing effect and the formation of lumps. Such a behaviour is similar to the appearance of breathers for the focusing nonlinear Schrodinger equation in the semiclassical limit.
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) …
Accretion shock on CTTSs and its X-ray emission
2009
High spectral resolution X-ray observations of classical T Tauri stars (CTTSs) demonstrate the presence of plasma at T~2-3×10^6 K and ne~10^11-10^13 cm-3. Stationary models suggest that this emission is due to shock-heated accreting material. We address this issue by a 1-D hydrodynamic model of the impact of the accretion flow onto a chromosphere of a CTTS with the aim of investigating the stability of accretion shock and the role of the chromosphere. Our simulations include the effects of gravity, radiative losses from optically thin plasma, the thermal conduction and a detailed modeling of the stellar chromosphere. Here we present the results of a simulation based on the parameters of the…
Recollimation shocks in relativistic jets
2018
Recollimation shocks (RS) appear associated with relativistic flows propagating through pressure mismatched atmospheres. Astrophysical scenarios invoking the presence of such shocks include jets from AGNs and X-ray binaries and GRBs. We shall start reviewing the theoretical background behind the structure of RS in overpressured jets. Next, basing on numerical simulations, we will focus on the properties of RS in relativistic steady jets threaded by helical magnetic fields depending on the dominant type of energy. Synthetic radio maps from the simulation of the synchrotron emission for a selection of models in the context of parsec-scale extragalactic jets will also be discussed.© 2018 World…
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.