Search results for " Numerical method"
showing 10 items of 42 documents
A shallow water SPH model with PML boundaries
2015
Abstract We focus on the study and implementation of Smoothed Particle Hydrodynamics (SPH) numerical code to deal with non-reflecting boundary conditions, starting from the Perfect Matched Layer (PML) approach. Basically, the method exploits the concept of a physical damping which acts on a fictitious layer added to the edges of computational domain. In this paper, we develop the study of time dependent shallow waves propagating on a finite 2D-XY plane domain and their behavior in the presence of circular and, more generic, rectangular boundary absorbing layers. In particular, an analysis of variation of the layer׳s thickness versus the absorbing efficiency is conducted. In our model, the m…
Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations
2004
We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments such as shock tubes and the relativistic flat-faced step test are performed. As an astrophysical application the article includes two-dimensional simulations of the propagation of relativistic jets using both Cartesian and cylindrical coordinates. The simulations reported clearly show the capabili…
Feedback Classification and Optimal Control with Applications to the Controlled Lotka-Volterra Model
2023
Let M be a σ-compact C^∞ manifold of dimension n ≥ 2 and consider a single-input control system: ẋ(t) = X (x(t)) + u(t) Y (x(t)), where X , Y are C^∞ vector fields on M. We prove that there exist an open set of pairs (X , Y ) for the C^∞ –Whitney topology such that they admit singular abnormal rays so that the spectrum of the projective singular Hamiltonian dynamics is feedback invariant. It is applied to controlled Lotka–Volterra dynamics where such rays are related to shifted equilibria of the free dynamics.
Size-intensive decomposition of orbital energy denominators
2000
We introduce an alternative to Almlöf and Häser’s Laplace transform decomposition of orbital energy denominators used in obtaining reduced scaling algorithms in perturbation theory based methods. The new decomposition is based on the Cholesky decomposition of positive semidefinite matrices. We show that orbital denominators have a particular short and size-intensive Cholesky decomposition. The main advantage in using the Cholesky decomposition, besides the shorter expansion, is the systematic improvement of the results without the penalties encountered in the Laplace transform decomposition when changing the number of integration points in order to control the convergence. Applications will…
Modeling Atmospheric Turbulence via Rapid Distortion Theory: Spectral Tensor of Velocity and Buoyancy
2017
Abstract A spectral tensor model is presented for turbulent fluctuations of wind velocity components and temperature, assuming uniform vertical gradients in mean temperature and mean wind speed. The model is built upon rapid distortion theory (RDT) following studies by Mann and by Hanazaki and Hunt, using the eddy lifetime parameterization of Mann to make the model stationary. The buoyant spectral tensor model is driven via five parameters: the viscous dissipation rate ε, length scale of energy-containing eddies L, a turbulence anisotropy parameter , gradient Richardson number (Ri) representing the local atmospheric stability, and the rate of destruction of temperature variance . Model outp…
Relativistic MHD simulations of extragalactic jets
2005
We have performed a comprehensive parameter study of the morphology and dynamics of axisymmetric, magnetized, relativistic jets by means of numerical simulations. The simulations have been performed with an upgraded version of the GENESIS code which is based on a second-order accurate finite volume method involving an approximate Riemann solver suitable for relativistic ideal magnetohydrodynamic flows, and a method of lines. Starting from pure hydrodynamic models we consider the effect of a magnetic field of increasing strength (up to β ≡ |b|2/2p ≈ 3.3 times the equipartition value) and different topology (purely toroidal or poloidal). We computed several series of models investigating the …
Numerical study for a new methodology of flaws detection in train axles
2013
Train loads and travel speeds have increased over time, requiring more efficient non-destructive inspection methods. Railway axles are critical elements; despite being designed to last more than 20 years several cases of premature failure have been recorded. Train axles are inspected regularly, but the limits associated to the traditional inspection technologies create a growing interest towards new solutions. Here a novel non-destructive inspection method of in-service axles based on non-contact data collection is presented. The propagation of surface waves, generated by a thermo-elastic laser source, is investigated using a finite element method based on dynamic explicit integration. Coup…
Electrical conductance of carbon nanotubes with misaligned ends
2013
During a manufacturing process, when a straight carbon nanotube is placed on a substrate, e.g., production of transistors, its two ends are often misaligned. In this study, we investigate the effects of multiwall carbon nanotubes’ (MWCNTs) outer diameter and chirality on the change in conductance due to misalignment of the two ends. The length of the studied MWCNTs was 120 nm, while the diameters ranged between 4 and 7 nm. A mixed finite element-tight-binding approach was carefully designed to realize reduction in computational time by orders of magnitude in calculating the deformation-induced changes in the electrical transport properties of the nanotubes. Numerical results suggest that ar…
A multi-sphere particle numerical model for non-invasive investigations of neuronal human brain activity
2013
In this paper, a multi-sphere particle method is built- up in order to estimate the solution of the Poisson's equation with Neumann boundary conditions describing the neuronal human brain activity. The partial difierential equations governing the relationships between neural current sources and the data produced by neuroimaging technique, are able to compute the scalp potential and magnetic fleld distributions generated by the neural activity. A numerical approach is proposed with current dipoles as current sources and going on in the computation by avoiding the mesh construction. The current dipoles are into an homogeneous spherical domain modeling the head and the computational approach i…
Monotonic solution of flow and transport problems in heterogeneous media using Delaunay unstructured triangular meshes
2013
Transport problems occurring in porous media and including convection, diffusion and chemical reactions, can be well represented by systems of Partial Differential Equations. In this paper, a numerical procedure is proposed for the fast and robust solution of flow and transport problems in 2D heterogeneous saturated media. The governing equations are spatially discretized with unstructured triangular meshes that must satisfy the Delaunay condition. The solution of the flow problem is split from the solution of the transport problem and it is obtained with an approach similar to the Mixed Hybrid Finite Elements method, that always guarantees the M-property of the resulting linear system. The…