اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدThe Numerical Simulation of Radiative Shocks I: The elimination of numerical shock instabilities using a localized oscillation filter
43
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We address a numerical instability that arises in the directionally split computation of hydrodynamic flows when shock fronts are parallel to a grid plane. Transverse oscillations in pressure, density and temperature are produced that are exacerbated by thermal instability when cooling is present, forming post--shock `stripes. These are orthogonal to the classic post--shock ringing fluctuations. The resulting post--shock `striping substantially modifies the flow. We discuss three different methods to resolve this problem. These include (1) a method based on artificial viscosity; (2) grid--jittering and (3) a new localized oscillation filter that acts on specific grid cells in the shock front. These methods are tested using a radiative wall shock problem with an embedded shear layer. The artificial viscosity method is unsatisfactory since, while it does reduce post--shock ringing, it does not eliminate the stripes and the excessive shock broadening renders the calculation of cooling inaccurate, resulting in an incorrect shock location. Grid--jittering effectively counteracts striping. However, elsewhere on the grid, the shear layer is unphysically diffused and this is highlighted in an extreme case. The oscillation filter method removes stripes and permits other high velocity gradient regions of the flow to evolve in a physically acceptable manner. It also has the advantage of only acting on a small fraction of the cells in a two or three dimensional simulation and does not significantly impair performance.
قيم البحث
اقرأ أيضاً
We present 3D hydrodynamic adiabatic simulations of a shock interacting with a dense, elongated cloud. We compare how the nature of the interaction changes with the filaments length and its orientation to the shock, and with the shock Mach number and
the density contrast of the filament. We then examine the differences with respect to 3D spherical-cloud calculations. We find significant differences in the morphology of the interaction when M=10 and chi=100: in many cases 3 parallel rolls are formed, and spread further apart with time, and periodic vortex shedding can occur off the ends of oblique filaments. Sideways-on filaments are accelerated more quickly, and initially lose mass more quickly than spherical clouds due to their greater surface area to volume ratio. However, at late stages they lose mass more slowly, due to the reduced relative speed between the filament and the postshock flow. The acceleration and mixing timescales can vary by a factor of 2 as the filament orientation changes. Oblique filaments can achieve transverse velocities up to 10% of the shock speed. Some aspects of our simulations are compared against experimental and numerical work on rigid cylinders.
We propose a Hermite spectral method for the spatially inhomogeneous Boltzmann equation. For the inverse-power-law model, we generalize an approximate quadratic collision operator defined in the normalized and dimensionless setting to an operator for
arbitrary distribution functions. An efficient algorithm with a fast transform is introduced to discretize this new collision operator. The method is tested for one-dimensional benchmark microflow problems.
We present a Fortran 95 code for simulating the evolution of astrophysical systems using particles to represent the underlying fluid flow. The code is designed to be versatile, flexible and extensible, with modular options that can be selected either
at compile time or at run time. We include a number of general purpose modules describing a variety of physical processes commonly required in the astrophysical community. The code can be used as an N-body code to evolve a set of particles in two or three dimensions using either a Leapfrog or Runge-Kutta-Fehlberg integrator, with or without individual timesteps for each particle. Particles may interact gravitationally as $N$-body particles, and all or any subset may also interact hydrodynamically, using the Smoothed Particle Hydrodynamic (SPH) method. Massive point particles (`stars) which may accrete nearby SPH or $N$-body particles may also be included. The default free boundary conditions can be replaced by a module to include periodic boundaries. Cosmological expansion may also be included. An interface with special purpose `GRAPE hardware may also be selected. If available, forces obtained from the GRAPE coprocessors may be transparently substituted for those obtained from the default tree based calculation. The code may be run without modification on single processors or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. In comparison to the Gadget-2 code of Springel 2005, the gravitational force calculation is $approx 3.5 - 4.8$ times faster with VINE when run on 8 Itanium~2 processors in an SGI Altix, while producing nearly identical outcomes in our test problems. We present simulations of several test problems, including a merger simulation of two elliptical galaxies with 800000 particles.
A simple model of random Brownian walk of a spherical mesoscopic particle in viscous liquids is proposed. The model can be both solved analytically and simulated numerically. The analytic solution gives the known Eistein-Smoluchowski diffusion law $<
r^2> = Dt$ where the diffusion constant $D$ is expressed by the mass and geometry of a particle, the viscosity of a liquid and the average effective time between consecutive collisions of the tracked particle with liquid molecules. The latter allows to make a simulation of the Perrin experiment and verify in detailed study the influence of the statistics on the expected theoretical results. To avoid the problem of small statistics causing departures from the diffusion law we introduce in the second part of the paper the idea of so called Artificially Increased Statistics (AIS) and prove that within this method of experimental data analysis one can confirm the diffusion law and get a good prediction for the diffusion constant even if trajectories of just few particles immersed in a liquid are considered.
In the paper, the pricing of Quanto options is studied, where the underlying foreign asset and the exchange rate are correlated with each other. Firstly, we adopt Bayesian methods to estimate unknown parameters entering the pricing formula of Quanto
options, including the volatility of stock, the volatility of exchange rate and the correlation. Secondly, we compute and predict prices of different four types of Quanto options based on Bayesian posterior prediction techniques and Monte Carlo methods. Finally, we provide numerical simulations to demonstrate the advantage of Bayesian method used in this paper comparing with some other existing methods. This paper is a new application of the Bayesian methods in the pricing of multi-asset options.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
