Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat equation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

Time-dependent Photoionization Modeling of Warm Absorbers: High-resolution Spectra and Response to Flaring Light Curves

Time-dependent photoionization modeling of warm absorber (WA) outflows in active galactic nuclei can play an important role in understanding the interaction between WAs and the central black hole. The WA may be out of the equilibrium state because of the variable nature of the central continuum. In this paper, with the help of time dependent photoionization modeling, we study how the WA gas properties change with time and how it reacts to changing radiation fields. Incorporating a flaring incident light curve, we investigate the behavior of WAs using a photoionization code that simultaneously and consistently solves the time-dependent equations of level population, heating and cooling, and radiative transfer. We simulate the physical processes in the gas clouds, such as ionization, recombination, heating, cooling, and the transfer of ionizing radiation through the cloud and present high-resolution time-resolved absorption spectra. We demonstrate that time-dependent radiative transfer is important and the calculations that omit this effect quantitatively and systematically underestimate the absorption. Time-dependent photoionization models provide crucial insights into the characteristics of WAs and can be used to constrain their density and spatial distribution.

Astrophysicists↗

Time-Dependent Photoionization Modeling of Warm Absorbers: High-Resolution Spectra and Response to Flaring Light Curves

Time-dependent photoionization modeling of warm absorber (WA) outflows in active galactic nuclei can play an important role in understanding the interaction between WAs and the central black hole. The WA may be out of the equilibrium state because of the variable nature of the central continuum. In this paper, with the help of time dependent photoionization modeling, we study how the WA gas properties change with time and how it reacts to changing radiation fields. Incorporating a flaring incident light curve, we investigate the behavior of WAs using a photoionization code that simultaneously and consistently solves the time-dependent equations of level population, heating and cooling, and radiative transfer. We simulate the physical processes in the gas clouds, such as ionization, recombination, heating, cooling, and the transfer of ionizing radiation through the cloud and present high-resolution time-resolved absorption spectra. We demonstrate that time-dependent radiative transfer is important and the calculations that omit this effect quantitatively and systematically underestimate the absorption. Time-dependent photoionization models provide crucial insights into the characteristics of WAs and can be used to constrain their density and spatial distribution.

Dev R. Sadaula↗

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗

Flow in a model turbine stator

In view of the complex nature of the flowfield in the hot section of gas turbine engines, the need to predict heat transfer and flow losses, the possible appearance of separation and strong secondary flows, etc., the present effort is focusing upon a Navier-Stokes approach to the three dimensional turbine stator problem. The advantages of a full Navier-Stokes approach are clear since when combined with a suitable turbulence model these equations represent the flow and heat transfer physics. In particular, the Navier-Stokes equations accurately represent possible separated regions and regions of significant secondary flow. In addition, the Navier-Stokes approach allows representation of the entire flow field by a single set of equations, thus avoiding problems associated with representing different regions of the flow by different equations and then matching flow regions.

Buggeln, R. C.↗

Overview of Icing Physics Relevant to Scaling

An understanding of icing physics is required for the development of both scaling methods and ice-accretion prediction codes. This paper gives an overview of our present understanding of the important physical processes and the associated similarity parameters that determine the shape of Appendix C ice accretions. For many years it has been recognized that ice accretion processes depend on flow effects over the model, on droplet trajectories, on the rate of water collection and time of exposure, and, for glaze ice, on a heat balance. For scaling applications, equations describing these events have been based on analyses at the stagnation line of the model and have resulted in the identification of several non-dimensional similarity parameters. The parameters include the modified inertia parameter of the water drop, the accumulation parameter and the freezing fraction. Other parameters dealing with the leading edge heat balance have also been used for convenience. By equating scale expressions for these parameters to the values to be simulated a set of equations is produced which can be solved for the scale test conditions. Studies in the past few years have shown that at least one parameter in addition to those mentioned above is needed to describe surface-water effects, and some of the traditional parameters may not be as significant as once thought. Insight into the importance of each parameter, and the physical processes it represents, can be made by viewing whether ice shapes change, and the extent of the change, when each parameter is varied. Experimental evidence is presented to establish the importance of each of the traditionally used parameters and to identify the possible form of a new similarity parameter to be used for scaling.

Anderson, David N.↗

A Simplified, Closed-Form Method for Screening Spacecraft Orbital Heating Variations

A closed-form analytical technique has been developed to screen orbital average heating variations as a function of beta angle, altitude, surface area, and surface optical properties. Using planetary view factor equations for surfaces parallel-to and normal-to the local vertical, a cylindrical umbral shadow approximation, and a simplified albedo flux model, heating rate equations are formulated and then integrated to obtain orbital average heating. The results are compared to detailed analytical predictions using Monte Carlo integration and an assessment of error is presented.

Rickman, S. L.↗

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics of sub-micron cosmic dust particles

Laboratory tests with simulated micrometeoroids to measure the heat transfer coefficient are discussed. Equations for ablation path length for electrically accelerated micrometeoroids entering a gas target are developed which yield guidelines for the laboratory measurement of the heat transfer coefficient. Test results are presented for lanthanum hexaboride (LaB sub 6) microparticles in air, argon, and oxygen targets. The tests indicate the heat transfer coefficient has a value of approximately 0.9 at 30 km/sec, and that it increases to approximately unity at 50 km/sec and above. Test results extend to over 100 km/sec. Results are also given for two types of small particle detectors. A solid state capacitor type detector was tested from 0.61 km/sec to 50 km/sec. An impact ionization type detector was tested from 1.0 to 150 km/sec using LaB sub 6 microparticles.

Roy, N. L.↗

Computing Isentropic Flow Properties of Air/R-134a Mixtures

MACHRK is a computer program that calculates isentropic flow properties of mixtures of air and refrigerant R-134a (tetrafluoroethane), which are used in transonic aerodynamic testing in a wind tunnel at Langley Research Center. Given the total temperature, total pressure, static pressure, and mole fraction of R-134a in a mixture, MACHRK calculates the Mach number and the following associated flow properties: dynamic pressure, velocity, density, static temperature, speed of sound, viscosity, ratio of specific heats, Reynolds number, and Prandtl number. Real-gas effects are taken into account by treating the gases comprising the mixture as both thermally and calorically imperfect. The Redlich-Kwong equation of state for mixtures and the constant-pressure ideal heat-capacity equation for the mixture are used in combination with the departure- function approach of thermodynamics to obtain the equations for computing the flow properties. In addition to the aforementioned calculations for air/R-134a mixtures, a research version of MACHRK can perform the corresponding calculations for mixtures of air and R-12 (dichlorodifluoromethane) and for air/SF6 mixtures. [R-12 was replaced by R-134a because of environmental concerns. SF6 has been considered for use in increasing the Reynolds-number range.]

Kvaternik, Ray↗

GP-BayesOpInf

SAND2025-01851O GP-BayesOpInf is a software tool that uses algorithms to combine Gaussian process regression, principal component analysis, and linear Bayesian inference to produce a probabilistic reduced-order model for time-dependent systems. Numerical examples include the compressible Euler equations for an ideal gas, a heat diffusion process with a nonlinear reaction term, and a set of ordinary differential equations describing a compartmental model in epidemiology. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft

This report compiles the various research activities conducted under the auspices of the NASA Grant NAG3-1026, "Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft" during the period of April 1989 to April 1994. The effort involved the development of multigrid based algorithms and computer programs for the calculation of the flow and temperature fields generated by Short Take-off and Vertical Landing (STOVL) aircraft, while hovering in ground proximity. Of particular importance has been the interaction of the exhaust jets with the head wind which gives rise to the hot gas ingestion process. The objective of new STOVL designs to reduce the temperature of the gases ingested into the engine. The present work describes a solution algorithm for the multi-dimensional elliptic partial-differential equations governing fluid flow and heat transfer in general curvilinear coordinates. The solution algorithm is based on the multigrid technique which obtains rapid convergence of the iterative numerical procedure for the discrete equations. Initial efforts were concerned with the solution of the Cartesian form of the equations. This algorithm was applied to a simulated STOVL configuration in rectangular coordinates. In the next phase of the work, a computer code for general curvilinear coordinates was constructed. This was applied to model STOVL geometries on curvilinear grids. The code was also validated in model problems. In all these efforts, the standard k-Epsilon model was used.

Vanka, S. P.↗

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗

Coupled flow, thermal and structural analysis of aerodynamically heated panels

A finite element approach to coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.

Thornton, Earl A.↗

Coupled flow, thermal and structural analysis of aerodynamically heated panels

A finite element approach for coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.

Thornton, Earl A.↗