Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonequilibrium dynamics”

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 163 records · Page 9

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

Improved Heat Transfer Prediction for High-Speed Flows over Blunt Bodies using Adaptive Mixed-Element Unstructured Grids

Aerothermodynamic prediction is important for the design and analysis of many aerospace vehicles. Computational Fluid Dynamics (CFD) tools used for these predictions commonly rely on structured shock-aligned grids due to the strong shocks exhibited at high speeds. This work details the incorporation of an improved HLLE++ scheme for both perfect gas and thermochemical nonequilibrium flows which robustly captures shock waves on both non-shock-aligned and shock-aligned grids. The HLLE++ scheme reduces the impact of eigenvalue limiting on flow solutions which is required for Roe's scheme for high Mach number flows. An adapted grid approach is also demonstrated using mixed-element unstructured grids to improve heat prediction. Thin prismatic boundary layers are used with tetrahedra in the farfield and for capturing shock waves. Various results are presented for hypersonic flows over blunt bodies including cylinders, spheres, and capsules.

CFD↗

Improved Heat Transfer Prediction for High-Speed Flows over Blunt Bodies using Adaptive Mixed-Element Unstructured Grids

Aerothermodynamic prediction is important for the design and analysis of many aerospace vehicles. Computational Fluid Dynamics (CFD) tools used for these predictions commonly rely on structured shock-aligned grids due to the strong shocks exhibited at high speeds. This work details the incorporation of an improved HLLE++ scheme for both perfect gas and thermochemical nonequilibrium flows which robustly captures shock waves on both non-shock-aligned and shock-aligned grids. The HLLE++ scheme reduces the impact of eigenvalue limiting on flow solutions which is required for Roe's scheme for high Mach number flows. An adapted grid approach is also demonstrated using mixed-element unstructured grids to improve heat prediction. Thin prismatic boundary layers are used with tetrahedra in the farfield and for capturing shock waves. Various results are presented for hypersonic flows over blunt bodies including cylinders, spheres, and capsules.

CFD↗

Laser Rayleigh-Scattering During Space Shuttle Entry

Report presents detailed study of capabilities and requirements for equipment of proposed ultraviolet Rayleigh-scattering instrument to be carried aboard Space Shuttle. Density of atmosphere around flightpath measured. With accuracy and resolution, measurements adequate for detection of small-scale meteorological structure affecting analysis of flight dynamic data of reentering spacecraft. Also discusses extensions of concept to measurements of location of, and density as function of position in, shock wave of Space Shuttle. Measurements provide baseline data for verification of computer models of high-enthalpy hypersonic, nonequilibrium, and viscous conditions.

Mckenzie, Robert L.↗

A linear quadratic Gaussian with loop transfer recovery proximity operations autopilot for spacecraft

An automatic control scheme for spacecraft proximity operations is presented. The controller is capable of holding the vehicle at a prescribed location relative to a target, or maneuvering it to a different relative position using straight line-of-sight translations. The autopilot uses a feedforward loop to initiate and terminate maneuvers, and for operations at nonequilibrium set-points. A multivariate feedback loop facilitates precise position and velocity control in the presence of sensor noise. The feedback loop is formulated using the Linear Quadratic Gaussian (LQG) with Loop Transfer Recovery (LTR) design procedure. Linear models of spacecraft dynamics, adapted from Clohessey-Wiltshire Equations, are augmented and loop shaping techniques are applied to design a target feedback loop. The loop transfer recovery procedure is used to recover the frequency domain properties of the target feedback loop. The resulting compensator is integrated into an autopilot which is tested in a high fidelity Space Shuttle Simulator. The autopilot performance is evaluated for a variety of proximity operations tasks envisioned for future Shuttle flights.

Chen, George T.↗

Hypersonic CFD applications for the National Aero-Space Plane

Design and analysis of the NASP depends heavily upon developing the critical technology areas that cover the entire engineering design of the vehicle. These areas include materials, structures, propulsion systems, propellants, integration of airframe and propulsion systems, controls, subsystems, and aerodynamics areas. Currently, verification of many of the classical engineering tools relies heavily on computational fluid dynamics. Advances are being made in the development of CFD codes to accomplish nose-to-tail analyses for hypersonic aircraft. Additional details involving the partial development, analysis, verification, and application of the CFL3D code and the SPARK combustor code are discussed. A nonequilibrium version of CFL3D that is presently being developed and tested is also described. Examples are given of portion calculations for research hypersonic aircraft geometries and comparisons with experiment data show good agreement.

Richardson, Pamela F.↗

Real-Gas Aerothermodynamics Test Facilities

This chapter provides an overview of the current ground-based aerothermodynamic testing capabilities in Western Europe and the United States. The focus is on facilities capable of producing real-gas effects (dissociation, ionization, and thermochemical nonequilibrium) pertinent to the study of atmospheric flight in the Mach number range of 5 < M < 50. Perceived mission needs of interest to the Americans and Western Europeans are described where such real-gas flows are important. The role of Computational Fluid Dynamics (CFD) in modern ground testing is discussed, and the capabilities of selected American and European real-gas facilities are described. An update on the current instrumentation in aerothermodynamic testing is also outlined. Comments are made regarding the use of new facilities which have been brought on line during the past 3-5 years. Finally, future needs for aerothermodynamic testing, including instrumentation, are discussed and recommendations for implementation are reported.

James 0. Arnold↗

The Effects of B/L-Dependent Heating on the Formation and Evolution of a Multi-Threaded Prominence

We have developed a comprehensive, multi-threaded, three-dimensional model of the plasma dynamics and energetics of a prominence and its overlying arcade (Luna et al. 2012). In this model, the basic magnetic structure is that of two interacting sheared arcades, while the cool condensations composing the prominence are formed by the well-studied thermal nonequilibrium mechanism. In a given filament-channel flux tube, the mass is evaporated from the chromosphere by heating localized near the footpoints, and condenses in the form of transient blobs or a persistent thread. Our previous studies of thermal nonequilibrium used steady or impulsive heating functions with no dependence on local physical. However, parametric active-region models with steady heating proportional to B/L, where B is the flux-tube magnetic field strength at each footpoint and L is the flux-tube length, yield the best agreement with observations (e.g., Schrijver et al. 2008). We have determined the effects of this active-region heating function on our model for the formation and evolution of prominence mass. We have also expanded the range of our computational domain to include more ofthe overlying arcade (the so-called "cavity") than in Luna et al. (2012), and have increased the number of selected flux tubes from 125 to 533. We will illustrate the time-dependent plasma behavior produced by the B/L heating function with synthetic images in several ALA passbands, and compare the resulting prominence properties with those predicted by our model with steady heating.

Karpen, Judith T.↗

Magnetohydrodynamic equilibrium. III - Helically symmetric fields. IV - Nonequilibrium of nonsymmetric hydrodynamic topologies

It is pointed out that plasma confinement in stable equilibrium states constitutes a fundamental and still unresolved question in plasma astrophysics and thermonuclear fusion research. The problem has two parts related to the equilibrium states themselves and their mechanical stability. The question of the existence of general solutions of the field and fluid equations for the steady dynamical interaction of inviscid compressible fluids of high electrical conductivity with magnetic and gravity fields is considered. In the absence of fluid motions, the presented equations become the familiar equations of magnetostatics supplemented by an equation of state. Starting from this simplest case of magnetostatic equilibrium, the investigation proceeds to the more complex case of magnetohydrodynamic equilibrium. Examples of helically symmetric fields are presented to illustrate the use of the formulation for treating the dynamics of helically symmetric hydromagnetic flows.

Tsinganos, K. C.↗

A study of numerical methods for hyperbolic conservation laws with stiff source terms

In the present study of the behavior of typical numerical methods in the case of a model advection equation having a parameter-dependent source term, two approaches to the incorporation of the source terms are used: MacCormack-type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. The latter are found to perform slightly better. The model scalar equation is used to show that the incorrectness of the propagation speeds of discontinuities observed in the stiff case is due to the introduction of nonequilibrium values through numerical dissipation in the advection step.

Leveque, R. J.↗

Issues and approach to develop validated analysis tools for hypersonic flows: One perspective

Critical issues concerning the modeling of low density hypervelocity flows where thermochemical nonequilibrium effects are pronounced are discussed. Emphasis is on the development of validated analysis tools, and the activity in the NASA Ames Research Center's Aerothermodynamics Branch is described. Inherent in the process is a strong synergism between ground test and real gas computational fluid dynamics (CFD). Approaches to develop and/or enhance phenomenological models and incorporate them into computational flowfield simulation codes are discussed. These models were partially validated with experimental data for flows where the gas temperature is raised (compressive flows). Expanding flows, where temperatures drop, however, exhibit somewhat different behavior. Experimental data for these expanding flow conditions is sparse and reliance must be made on intuition and guidance from computational chemistry to model transport processes under these conditions. Ground based experimental studies used to provide necessary data for model development and validation are described. Included are the performance characteristics of high enthalpy flow facilities, such as shock tubes and ballistic ranges.

Deiwert, George S.↗

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗

Issues and approach to develop validated analysis tools for hypersonic flows: One perspective

Critical issues concerning the modeling of low-density hypervelocity flows where thermochemical nonequilibrium effects are pronounced are discussed. Emphasis is on the development of validated analysis tools. A description of the activity in the Ames Research Center's Aerothermodynamics Branch is also given. Inherent in the process is a strong synergism between ground test and real-gas computational fluid dynamics (CFD). Approaches to develop and/or enhance phenomenological models and incorporate them into computational flow-field simulation codes are discussed. These models have been partially validated with experimental data for flows where the gas temperature is raised (compressive flows). Expanding flows, where temperatures drop, however, exhibit somewhat different behavior. Experimental data for these expanding flow conditions are sparse; reliance must be made on intuition and guidance from computational chemistry to model transport processes under these conditions. Ground-based experimental studies used to provide necessary data for model development and validation are described. Included are the performance characteristics of high-enthalpy flow facilities, such as shock tubes and ballistic ranges.

Deiwert, George S.↗

Soft X-ray spectroscopy of solar flares - An overview

An overview of the current status of high spectral resolution soft X-ray observations of solar flares is given. The review concentrates primarily on recent results and interpretations of results obtained from orbiting Bragg crystal spectrometers flow during the last solar maximum on the US DoD P78-1 spacecraft, the NASA SMM, and the ISAS Hinotori spacecraft. Results and several key issues regarding interpretation of the spectra are presented. Specifically, the dynamics of coronal flare plasmas as revealed by X-ray line profiles and wavelength shifts are discussed. Recent results concerning the theory of chromospheric evaporation are given. The temperature of coronal flare plasma is discussed within the context of a differential mission measure. Results concerning electron density measurements, nonequilibrium processes, and relative element abundances are also reviewed.

Doschek, G. A.↗

Interstellar cloud evolution and the abundance of formaldehyde

The time scale for essentially complete conversion of C(+) to CO in interstellar clouds can be comparable to, or greater than, dynamical time scales for evolution, therefore suggesting steady state time independent abundances to be inappropriate. The solutions for the time-dependent carbon chemistry in dense clouds, with density not less than 500/cu cm, indicate that significant amounts of neutral carbon will be present throughout a cloud's lifetime. These nonequilibrium values of C I can explain the relatively large abundances observed for formaldehyde, isotopes of carbon monoxide, and other trace molecules

Langer, W. D.↗

Computational flow predictions for hypersonic drag devices

The effectiveness of two types of hypersonic decelerators is examined: mechanically deployable flares and inflatable ballutes. Computational fluid dynamics (CFD) is used to predict the flowfield around a solid rocket motor (SRM) with a deployed decelerator. The computations are performed with an ideal gas solver using an effective specific heat ratio of 1.15. The results from the ideal gas solver are compared to computational results from a thermochemical nonequilibrium solver. The surface pressure coefficient, the drag, and the extend of the compression corner separation zone predicted by the ideal gas solver compare well with those predicted by the nonequilibrium solver. The ideal gas solver is computationally inexpensive and is shown to be well suited for preliminary design studies. The computed solutions are used to determine the size and shape of the decelerator that are required to achieve a drag coefficient of 5. Heat transfer rates to the SRM and the decelerators are predicted to estimate the amount of thermal protection required.

Tokarcik, Susan A.↗

What Dominates the Coronal Emission Spectrum During the Cycle of Impulsive Heating and Cooling?

The smoking gun of small-scale, impulsive events heating the solar corona is expected to be the presence of a hot ( > 5 MK) plasma component. Evidence for this has been scarce, but has gradually begun to accumulate due to recent studies designed to constrain the high temperature part of the emission measure distribution. However, the detected hot component is often weaker than models predict and this is due in part to the common modeling assumption that the ionization balance remains in equilibrium. The launch of the latest generation of space-based observing instrumentation aboard Hinode and the Solar Dynamics Observatory (SDO) has brought the matter of the ionization state of the plasma firmly to the forefront. It is timely to consider exactly what emission current instruments would detect when observing a corona heated impulsively on small-scales by nanoflares. Only after we understand the full effects of nonequilibrium ionization can we draw meaningful conclusions about the plasma that is (or is not) present. We have therefore performed a series of hydrodynamic simulations for a variety of different nanoflare properties and initial conditions. Our study has led to several key conclusions. 1. Deviations from equilibrium are greatest for short-duration nanoflares at low initial coronal densities. 2. Hot emission lines are the most affected and are suppressed sometimes to the point of being invisible. 3. The emission detected in all of the SDO-AIA channels is generally dominated by warm, over-dense, cooling plasma. 4. It is difficult not to create coronal loops that emit strongly at 1.5 MK and in the range 2 to 5 MK, which are the most commonly observed kind, for a broad range of nanoflare scenarios. 5. The Fe XV (284.16 ) emission in most of our models is about 10 times brighter than the Ca XVII (192.82 ) emission, consistent with observations. Our overarching conclusion is that small-scale, impulsive heating inducing a nonequilibrium ionization state leads to predictions for observable quantities that are entirely consistent with what is actually observed.

Bradshaw, Stephen J.↗