2,080 research outputs found
A unified gas-kinetic particle method for frequency-dependent radiative transfer equations with isotropic scattering process on unstructured mesh
In this paper, we extend the unified kinetic particle (UGKP) method to the
frequency-dependent radiative transfer equation with both absorption-emission
and scattering processes. The extended UGKP method could not only capture the
diffusion and free transport limit, but also provide a smooth transition in the
physical and frequency space in the regime between the above two limits. The
proposed scheme has the properties of asymptotic-preserving, regime-adaptive,
and entropy-preserving, which make it an accurate and efficient scheme in the
simulation of multiscale photon transport problems. The methodology of scheme
construction is a coupled evolution of macroscopic energy equation and the
microscopic radiant intensity equation, where the numerical flux in macroscopic
energy equation and the closure in microscopic radiant intensity equation are
constructed based on the integral solution. Both numerical dissipation and
computational complexity are well controlled especially in the optical thick
regime. A 2D multi-thread code on a general unstructured mesh has been
developed. Several numerical tests have been simulated to verify the numerical
scheme and code, covering a wide range of flow regimes. The numerical scheme
and code that we developed are highly demanded and widely applicable in the
high energy density engineering applications
Research and Technology
Langley Research Center is engaged in the basic an applied research necessary for the advancement of aeronautics and space flight, generating advanced concepts for the accomplishment of related national goals, and provding research advice, technological support, and assistance to other NASA installations, other government agencies, and industry. Highlights of major accomplishments and applications are presented
Best practices for modeling structural boundary conditions due to a localized fire
Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/154450/1/fam2774.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/154450/2/fam2774_am.pd
A High-Order Low-Order Algorithm with Exponentially-Convergent Monte Carlo for Thermal Radiative Transfer Problems
We have implemented a new high-order low-order (HOLO) algorithm for solving thermal radiative transfer (TRT) problems. Within each discrete time step, fixed-point iterations are performed between a high-order (HO) exponentially-convergent Monte Carlo (ECMC) solver and a low-order (LO) system of equations. The LO system is based on spatial and angular moments of the transport equation and a linear-discontinuous finite-element (LDFE) spatial representation, producing equations similar to the standard S2 equations. The LO solver is fully implicit in time and efficiently converges the non-linear temperature dependence with Newton's method. The HO solver provides a globally accurate solution for the angular intensity to a fixed-source, pure absorber transport problem. This global solution is used to compute consistency terms in the LO equations that require the HO and LO solutions to converge towards the same solution. The use of ECMC allows for the efficient reduction of statistical noise in the solution.
We investigated several extensions of this algorithm. A parametric closure of the LO system was used for the spatial variable, based on local relations computed with the HO solver. The spatial closure improves consistency between the two solvers compared to a standard LDFE spatial discretization of the LO system. The ECMC algorithm has been extended to integrate the angular intensity in time, with a consistent time closure of the LO radiation equations. The time closure increases accuracy in optically-thin problems compared to a backward Euler discretization. Finally, we have applied standard source iteration and Krylov procedures to iteratively solve the LO equations, with linear diffusion synthetic acceleration.
Herein, we present results for one-dimensional, gray test problems. Results demonstrate several desirable properties of this algorithm: the HOLO method preserves the equilibrium diffusion limit, prevents violation of the maximum principle, and can provide high-fidelity MC solutions to the TRT equations with minimal statistical noise. We have compared results with an implicit Monte Carlo (IMC) code and compared the efficiency of ECMC to standard Monte Carlo in this HOLO algorithm. Our HOLO algorithm is more accurate and more efficient than standard IMC. The extent to which this is so is problem-dependent
Perspectives on hypersonic viscous and nonequilibrium flow research
An attempt is made to reflect on current focuses in certain areas of hypersonic flow research by examining recent works and their issues. Aspects of viscous interaction, flow instability, and nonequilibrium aerothermodynamics pertaining to theoretical interest are focused upon. The field is a diverse one, and many exciting works may have either escaped the writer's notice or been abandoned for the sake of space. Students of hypersonic viscous flow must face the transition problems towards the two opposite ends of the Reynolds or Knudsen number range, which represents two regimes where unresolved fluid/gas dynamic problems abound. Central to the hypersonic flow studies is high-temperature physical gas dynamics; here, a number of issues on modelling the intermolecular potentials and inelastic collisions remain the obstacles to quantitative predictions. Research in combustion and scramjet propulsion will certainly be benefitted by advances in turbulent mixing and new computational fluid dynamics (CFD) strategies on multi-scaled complex reactions. Even for the sake of theoretical development, the lack of pertinent experimental data in the right energy and density ranges is believed to be among the major obstacles to progress in aerothermodynamic research for hypersonic flight. To enable laboratory simulation of nonequilibrium effects anticipated for transatmospheric flight, facilities capable of generating high enthalpy flow at density levels higher than in existing laboratories are needed (Hornung 1988). A new free-piston shock tunnel capable of realizing a test-section stagnation temperature of 10(exp 5) at Reynolds number 50 x 10(exp 6)/cm is being completed and preliminary tests has begun (H. Hornung et al. 1992). Another laboratory study worthy of note as well as theoretical support is the nonequilibrium flow experiment of iodine vapor which has low activation energies for vibrational excitation and dissociation, and can be studied in a laboratory with modest resources (Pham-Van-Diep et al. 1992)
Radiative transfer with scattering for domain-decomposed 3D MHD simulations of cool stellar atmospheres
We present the implementation of a radiative transfer solver with coherent
scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD)
simulations of stellar surface convection. The code is fully parallelized using
MPI domain decomposition, which allows for large grid sizes and improved
resolution of hydrodynamical structures. We apply the code to simulate the
surface granulation in a solar-type star, ignoring magnetic fields, and
investigate the importance of coherent scattering for the atmospheric
structure. A scattering term is added to the radiative transfer equation,
requiring an iterative computation of the radiation field. We use a
short-characteristics-based Gauss-Seidel acceleration scheme to compute
radiative flux divergences for the energy equation. The effects of coherent
scattering are tested by comparing the temperature stratification of three 3D
time-dependent hydrodynamical atmosphere models of a solar-type star: without
scattering, with continuum scattering only, and with both continuum and line
scattering. We show that continuum scattering does not have a significant
impact on the photospheric temperature structure for a star like the Sun.
Including scattering in line-blanketing, however, leads to a decrease of
temperatures by about 350\,K below log tau < -4. The effect is opposite to that
of 1D hydrostatic models in radiative equilibrium, where scattering reduces the
cooling effect of strong LTE lines in the higher layers of the photosphere.
Coherent line scattering also changes the temperature distribution in the high
atmosphere, where we observe stronger fluctuations compared to a treatment of
lines as true absorbers.Comment: A&A, in pres
Grand Challenges in Protoplanetary Disc Modelling
The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th
--11th March 2016 --- included several open sessions led by participants. This
paper reports on the discussions collectively concerned with the multiphysics
modelling of protoplanetary discs, including the self-consistent calculation of
gas and dust dynamics, radiative transfer and chemistry. After a short
introduction to each of these disciplines in isolation, we identify a series of
burning questions and grand challenges associated with their continuing
development and integration. We then discuss potential pathways towards solving
these challenges, grouped by strategical, technical and collaborative
developments. This paper is not intended to be a review, but rather to motivate
and direct future research and collaboration across typically distinct fields
based on \textit{community driven input}, to encourage further progress in our
understanding of circumstellar and protoplanetary discs
- …