96 research outputs found
A Hybrid Godunov Method for Radiation Hydrodynamics
From a mathematical perspective, radiation hydrodynamics can be thought of as
a system of hyperbolic balance laws with dual multiscale behavior (multiscale
behavior associated with the hyperbolic wave speeds as well as multiscale
behavior associated with source term relaxation). With this outlook in mind,
this paper presents a hybrid Godunov method for one-dimensional radiation
hydrodynamics that is uniformly well behaved from the photon free streaming
(hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and
to the strong equilibrium diffusion (hyperbolic) limit. Moreover, one finds
that the technique preserves certain asymptotic limits. The method incorporates
a backward Euler upwinding scheme for the radiation energy density and flux as
well as a modified Godunov scheme for the material density, momentum density,
and energy density. The backward Euler upwinding scheme is first-order accurate
and uses an implicit HLLE flux function to temporally advance the radiation
components according to the material flow scale. The modified Godunov scheme is
second-order accurate and directly couples stiff source term effects to the
hyperbolic structure of the system of balance laws. This Godunov technique is
composed of a predictor step that is based on Duhamel's principle and a
corrector step that is based on Picard iteration. The Godunov scheme is
explicit on the material flow scale but is unsplit and fully couples matter and
radiation without invoking a diffusion-type approximation for radiation
hydrodynamics. This technique derives from earlier work by Miniati & Colella
2007. Numerical tests demonstrate that the method is stable, robust, and
accurate across various parameter regimes.Comment: accepted for publication in Journal of Computational Physics; 61
pages, 15 figures, 11 table
Proceedings for the ICASE Workshop on Heterogeneous Boundary Conditions
Domain Decomposition is a complex problem with many interesting aspects. The choice of decomposition can be made based on many different criteria, and the choice of interface of internal boundary conditions are numerous. The various regions under study may have different dynamical balances, indicating that different physical processes are dominating the flow in these regions. This conference was called in recognition of the need to more clearly define the nature of these complex problems. This proceedings is a collection of the presentations and the discussion groups
Mini-Workshop: Innovative Trends in the Numerical Analysis and Simulation of Kinetic Equations
In multiscale modeling hierarchy, kinetic theory plays a vital role in connecting microscopic Newtonian mechanics and macroscopic continuum mechanics. As computing power grows, numerical simulation of kinetic equations has become possible and undergone rapid development over the past decade. Yet the unique challenges arising in these equations, such as highdimensionality, multiple scales, random inputs, positivity, entropy dissipation, etc., call for new advances of numerical methods. This mini-workshop brought together both senior and junior researchers working on various fastpaced growing numerical aspects of kinetic equations. The topics include, but were not limited to, uncertainty quantification, structure-preserving methods, phase transitions, asymptotic-preserving schemes, and fast methods for kinetic equations
Solving the system of radiation magnetohydrodynamics for solar physical simulations in 3d
In this study we present a finite-volumen scheme for solving the equations of radiation magnetohydrodynamics in two and three space dimensions. Among other applications this system is used to model the plasma in the solar convection zone and in the solar photosphere. It is a non--linear system of balance laws derived from the Euler equations of gas dynamics and the Maxwell equations; the energy transport through radiation is also included in the model. The starting point of our presentation is a standard explicit first and second order finite-volume scheme on both structured and unstructured grids. We first study the convergence of a finite-volume scheme applied to a scalar model problem for the full system of radiation magnetohydrodynamics. We then present modifications of the base scheme. These make it possible to approximate the system of magnetohydrodynamics with an arbitrary equation of state; they reduce errors due to a violation of the divergence constraint on the magnetic field, and they lead to an improved accuracy in the approximation of solution near an equilibrium state. These modifications significantly increase the robustness of the scheme and are essential for an accurate simulation of processes in the solar atmosphere ...thesi
Advanced Computational Fluid Dynamics for Emerging Engineering Processes
As researchers deal with processes and phenomena that are geometrically complex and phenomenologically coupled the demand for high-performance computational fluid dynamics (CFD) increases continuously. The intrinsic nature of coupled irreversibility requires computational tools that can provide physically meaningful results within a reasonable time. This book collects the state-of-the-art CFD research activities and future R&D directions of advanced fluid dynamics. Topics covered include in-depth fundamentals of the Navier-Stokes equation, advanced multi-phase fluid flow, and coupling algorithms of computational fluid and particle dynamics. In the near future, true multi-physics and multi-scale simulation tools must be developed by combining micro-hydrodynamics, fluid dynamics, and chemical reactions within an umbrella of irreversible statistical physics
Towards a solution of the closure problem for convective atmospheric boundary-layer turbulence
We consider the closure problem for turbulence in the dry convective atmospheric boundary
layer (CBL). Transport in the CBL is carried by small scale eddies near the surface and large
plumes in the well mixed middle part up to the inversion that separates the CBL from the
stably stratified air above. An analytically tractable model based on a multivariate Delta-PDF
approach is developed. It is an extension of the model of Gryanik and Hartmann [1] (GH02)
that additionally includes a term for background turbulence. Thus an exact solution is derived
and all higher order moments (HOMs) are explained by second order moments, correlation
coefficients and the skewness. The solution provides a proof of the extended universality
hypothesis of GH02 which is the refinement of the Millionshchikov hypothesis (quasi-
normality of FOM). This refined hypothesis states that CBL turbulence can be considered as
result of a linear interpolation between the Gaussian and the very skewed turbulence regimes.
Although the extended universality hypothesis was confirmed by results of field
measurements, LES and DNS simulations (see e.g. [2-4]), several questions remained
unexplained. These are now answered by the new model including the reasons of the
universality of the functional form of the HOMs, the significant scatter of the values of the
coefficients and the source of the magic of the linear interpolation. Finally, the closures
61
predicted by the model are tested against measurements and LES data. Some of the other
issues of CBL turbulence, e.g. familiar kurtosis-skewness relationships and relation of area
coverage parameters of plumes (so called filling factors) with HOM will be discussed also
- …