Search results for "Numerical"
showing 10 items of 2002 documents
Investigation of an entropic stabilizer for the lattice-Boltzmann method
2015
The lattice-Boltzmann (LB) method is commonly used for the simulation of fluid flows at the hydrodynamic level of description. Due to its kinetic theory origins, the standard LB schemes carry more degrees of freedom than strictly needed, e.g., for the approximation of solutions to the Navier-stokes equation. In particular, there is freedom in the details of the so-called collision operator. This aspect was recently utilized when an entropic stabilizer, based on the principle of maximizing local entropy, was proposed for the LB method [I. V. Karlin, F. Bosch, and S. S. Chikatamarla, ¨ Phys. Rev. E 90, 031302(R) (2014)]. The proposed stabilizer can be considered as an add-on or extension to b…
A numerical study of postshock oscillations in slowly moving shock waves
2003
Abstract Godunov-type methods and other shock capturing schemes can display pathological behavior in certain flow situations. This paper discusses the numerical anomaly associated to slowly moving shocks. We present a series of numerical experiments that illustrate the formation and propagation of this pathology, and allows us to establish some conclusions and question some previous conjectures for the source of the numerical noise. A simple diagnosis on an explicit Steger-Warming scheme shows that some intermediate states in the first time steps deviate from the true direction and contaminate the flow structure. A remedy is presented in the form of a new flux split method with an entropy i…
Hydrodynamic modelling of ejecta shrapnel in the Vela supernova remnant
2013
Many supernova remnants (SNRs) are characterized by a knotty ejecta structure. The Vela SNR is an excellent example of remnant in which detached clumps of ejecta are visible as X-ray emitting bullets that have been observed and studied in great detail. We aim at modelling the evolution of ejecta shrapnel in the Vela SNR, investigating the role of their initial parameters (position and density) and addressing the effects of thermal conduction and radiative losses. We performed a set of 2-D hydrodynamic simulations describing the evolution of a density inhomogeneity in the ejecta profile. We explored different initial setups. We found that the final position of the shrapnel is very sensitive …
Probing the Internal Structure of Magnetized, Relativistic Jets with Numerical Simulations
2016
From an observational point of view, unveiling the physical processes behind the nature of the jets emanating from radio-loud AGN demands the resolution of the structure across the jet with the highest angular resolutions. Relying on a magneto-fluid dynamical description, numerical simulations can help to characterize the internal structure of jets (transversal structure, magnetic field structure, internal shocks, etc.). In the first part of the paper, we shall discuss equilibrium models of magnetized, relativistic, infinite, axisymmetric jets with rotation propagating through a homogeneous, static, unmagnetized ambient medium. Then, these transversal equilibrium profiles will be used to bu…
Millimeter-wave and microwave signal generation by low-bandwidth electro-optic phase modulation
2009
We propose, analyze and numerically illustrate a photonic-based technique for waveform generation of electrical signals approaching the 50 GHz bandwidth with time apertures as large as a few nanoseconds, by low-frequency, up to 2 GHz, electro-optic phase modulation of time-stretched optical pulses. Synthesis of the electrical waveform relies on phase-to-amplitude conversion of the modulated signal by a group delay dispersion circuit designed to behave as a transversal filter with N taps. Although arbitrary waveform generation capabilities are limited, a wide variety of user-defined signals are numerically demonstrated by appropriately designing the low-frequency signal driving the electro-o…
Numerical study of the primitive equations in the small viscosity regime
2018
In this paper we study the flow dynamics governed by the primitive equations in the small viscosity regime. We consider an initial setup consisting on two dipolar structures interacting with a no slip boundary at the bottom of the domain. The generated boundary layer is analyzed in terms of the complex singularities of the horizontal pressure gradient and of the vorticity generated at the boundary. The presence of complex singularities is correlated with the appearance of secondary recirculation regions. Two viscosity regimes, with different qualitative properties, can be distinguished in the flow dynamics.
Analysis of Cylindrical Dielectric Resonators in Rectangular Cavities Using a State-Space Integral-Equation Method
2006
In this letter, a state-space integral-equation method in the s-domain has been employed for the accurate analysis of rectangular cavities loaded with cylindrical dielectric resonators. The dielectric obstacles have been treated in terms of their polarization equivalent charge and current densities. The dielectric resonator can be placed at any arbitrary position inside the cavity. The presented technique allows to calculate in a very efficient way a large number of solenoidal modes. The resonant frequencies of dielectric-loaded cavities are calculated and compared with data from literature and a commercial finite element method software, showing a good agreement
A Coupled Solid-Fluid Method for Modeling Subduction
2007
International audience; We present a novel dynamic approach for solid/fluid coupling by joining two different numerical methods: Boundary Element Method (BEM) and Finite Element Method (FEM). FEM results describe the thermo-mechanical evolution of the solid while the fluid is solved with the BEM. The bidirectional feedback between the two domains evolves along a Lagrangian interface where the FEM domain is embedded inside the BEM domain. The feedback between the two codes is based on the calculation of a specific drag tensor for each boundary/finite element. The approach is presented here to solve the complex problem of the descent of a cold subducting oceanic plate into a hot fluid like ma…
Numerical evolution of matter in dynamical axisymmetric black hole spacetimes
2000
We have developed a numerical code to study the evolution of self-gravitating matter in dynamic black hole axisymmetric spacetimes in general relativity. The matter fields are evolved with a high-resolution shock-capturing scheme that uses the characteristic information of the general relativistic hydrodynamic equations to build up a linearized Riemann solver. The spacetime is evolved with an axisymmetric ADM code designed to evolve a wormhole in full general relativity. We discuss the numerical and algorithmic issues related to the effective coupling of the hydrodynamical and spacetime pieces of the code, as well as the numerical methods and gauge conditions we use to evolve such spacetime…
Generalized Langevin dynamics: construction and numerical integration of non-Markovian particle-based models.
2018
We propose a generalized Langevin dynamics (GLD) technique to construct non-Markovian particle-based coarse-grained models from fine-grained reference simulations and to efficiently integrate them. The proposed GLD model has the form of a discretized generalized Langevin equation with distance-dependent two-particle contributions to the self- and pair-memory kernels. The memory kernels are iteratively reconstructed from the dynamical correlation functions of an underlying fine-grained system. We develop a simulation algorithm for this class of non-Markovian models that scales linearly with the number of coarse-grained particles. Our GLD method is suitable for coarse-grained studies of syste…