Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonequilibrium physics”

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 127 records · Page 7

Comparisons of thermochemical nonequilibrium viscous flowfield predictions for AFE vehicle

Numerical solutions for the Aeroassist Flight Experiment vehicle were obtained from three methods at a trajectory point corresponding to the maximum aerodynamic heating. The flow regime and vehicle's speed require a viscous model of laminar flow and finite-rate chemical and thermal modeling. The computational domain covers both forebody and base such that the shock layer and near wake flowfield are included. Because of differences in computational grids, methods of solution, and models of rate equations, the results are generally in poor agreement. Temperature and species concentrations are strongly affected by the physical model equations and associated parameters. The chemistry model based on 11 species is found to yield lower translational temperature profile near stagnation than those from a seven-species model. The vibrational temperature varies according to the modeling details. All solutions indicate that strong neutral and/or molecular dissociation and weak ionization take place at the forebody and vibrational freezing is present in the afterbody expansion region where the vibrational temperature is higher than the translational temperature. Some forms of shear layer emanating from the aerobrake skirt coalesce in the region of reversed flow behind the vehicle.

Tam, Luen T.↗

Elements of radiative interactions in gaseous systems

Basic formulations, analyses, and numerical procedures are presented to study radiative interactions in gray as well as nongray gases under different physical and flow conditions. After preliminary fluid-dynamical considerations, essential governing equations for radiative transport are presented that are applicable under local and nonlocal thermodynamic equilibrium conditions. Auxiliary relations for relaxation times and spectral absorption model are also provided. For specific applications, several simple gaseous systems are analyzed. The first system considered consists of a gas bounded by two parallel plates having the same temperature. For this system, both vibrational nonequilibrium effects and radiation conduction interactions are studied. The second system consists of fully developed laminar flow and heat transfer in a parallel plate duct under the boundary condition of a uniform surface heat flux. For this system, effects of gray surface emittance are studied. With the single exception of a circular geometry, the third system is identical to the second system. Here, the influence of nongray walls is also studied, and a correlation between the parallel plates and circular tube results is presented. The particular gases selected are CO, CO2, H2O, CH4, N2O, NH3, OH, and NO. The temperature and pressure range considered are 300 to 2000 K, and 0.1 to 100 atmosphere, respectively. Illustrative results obtained for different cases are discussed and some specific conclusions are provided.

Tiwari, Surendra N.↗

CME Initiation

One of the most challenging problems in solar physics is understanding the processes responsible for giant magnetic disruptions such as the event of July 14, 2000, which consisted of a massive filament ejection, a fast coronal mass ejection (CME), prolonged X-class flaring, and an intense particle storm. These major events are of critical importance because they drive the most destructive forms of space weather and they provide a unique opportunity to study, in revealing detail, MHD instability and nonequilibrium -- processes that are at the heart of plasma astrophysics. It is now widely accepted that CMEs/eruptive flares represent the explosive release of magnetic energy stored in the corona. Therefore, in order to understand the phenomenon, we must answer the following questions: What is the field structure responsible for the disruption and why is the energy released explosively? In this talk we address these two questions using the latest theories and numerical models for CMEs/eruptive flares.

Antiochos, Spiro K.↗

The measurement of boundary layers on a compressor blade in cascade. Volume 1: Experimental technique, analysis and results

Measurements were made of the boundary layers and wakes about a highly loaded, double-circular-arc compressor blade in cascade. These laser Doppler velocimetry measurements have yielded a very detailed and precise data base with which to test the application of viscous computational codes to turbomachinery. In order to test the computational codes at off-design conditions, the data were acquired at a chord Reynolds number of 500,000 and at three incidence angles. Moreover, these measurements have supplied some physical insight into these very complex flows. Although some natural transition is evident, laminar boundary layers usually detach and subsequently reattach as either fully or intermittently turbulent boundary layers. These transitional separation bubbles play an important role in the development of most of the boundary layers and wakes measured in this cascade and the modeling or computing of these bubbles should prove to be the key aspect in computing the entire cascade flow field. In addition, the nonequilibrium turbulent boundary layers on these highly loaded blades always have some region of separation near the trailing edge of the suction surface. These separated flows, as well as the subsequent near wakes, show no similarity and should prove to be a challenging test for the viscous computational codes.

Zierke, William C.↗

Computational Aerothermodynamic Simulation Issues on Unstructured Grids

The synthesis of physical models for gas chemistry and turbulence from the structured grid codes LAURA and VULCAN into the unstructured grid code FUN3D is described. A directionally Symmetric, Total Variation Diminishing (STVD) algorithm and an entropy fix (eigenvalue limiter) keyed to local cell Reynolds number are introduced to improve solution quality for hypersonic aeroheating applications. A simple grid-adaptation procedure is incorporated within the flow solver. Simulations of flow over an ellipsoid (perfect gas, inviscid), Shuttle Orbiter (viscous, chemical nonequilibrium) and comparisons to the structured grid solvers LAURA (cylinder, Shuttle Orbiter) and VULCAN (flat plate) are presented to show current capabilities. The quality of heating in 3D stagnation regions is very sensitive to algorithm options in general, high aspect ratio tetrahedral elements complicate the simulation of high Reynolds number, viscous flow as compared to locally structured meshes aligned with the flow.

Gnoffo, Peter A.↗

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon similar correlation techniques that have been is use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. Recent entry flight testing performed with the Orbiter Discovery now provides a means to validate this engineering correlation approach to higher confidence. These results only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, and to discuss the characteristics which are desirable in a next generation advancement. The details of the paper will review the experimental facilities and techniques that were utilized to perform the implementation of the Orbiter RTF BLT Vsn 2 prediction capability. Statistically significant results for multiple engineering correlations from a ground testing campaign will be reviewed in order to describe why only certain correlations were selected for complete implementation to support the Shuttle Program. Historical Orbiter flight data on early boundary layer transition due to protruding gap fillers will be described in relation to the selected empirical correlations. In addition, Orbiter entry flight testing results from the BLT Flight Experiment will be discussed in relation to these correlations. Applicability of such correlations to the entry design problem will be reviewed, and finally a perspective on the desirable characteristics for a next generation capability based on high fidelity physical models will be provided.

Campbell, Charles H.↗

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.

Maccormack, Robert W.↗

The propagation of torsion along flux tubes subject to dynamical nonequilibrium

It is noted that the dynamical nonequilibrium of close-packed flux tubes is driven by the torsion in the individual tubes. Because of this, whenever tubes with the same sense of twisting come into contact, there is reconnection of their azimuthal field components. The reconnection consumes the local torsion, and this causes the propagation of torsional Alfven waves into the region from elsewhere along the tubes. The formal problem of the propagation of the torsion along twisted flux tubes is presented, along with some of the basic physical properties worked out in the limit of small torsion. It is noted that in tubes with finite twisting the propagation of torsional Alfven waves can be a more complicated phenomenon. Application to the sun suggests that the propagation of torsion from below the visible surface up into the corona is an important energy supply to the corona for a period of perhaps 10-20 hours after the emergence of the flux tubes through the surface of the sun, bringing up torsion from depths of 10,000 km or more. Torsion is of course continually furnished by the manipulation and shuffling of the field by the convection.

Parker, E. N.↗

Radiative energy transfer in molecular gases

Basic formulations, analyses, and numerical procedures are presented to study radiative interactions in gray as well as nongray gases under different physical and flow conditions. After preliminary fluid-dynamical considerations, essential governing equations for radiative transport are presented that are applicable under local and nonlocal thermodynamic equilibrium conditions. Auxiliary relations for relaxation times and spectral absorption models are also provided. For specific applications, several simple gaseous systems are analyzed. The first system considered consists of a gas bounded by two parallel plates having the same temperature. Within the gas there is a uniform heat source per unit volume. For this system, both vibrational nonequilibrium effects and radiation conduction interactions are studied. The second system consists of fully developed laminar flow and heat transfer in a parallel plate duct under the boundary condition of a uniform surface heat flux. For this system, effects of gray surface emittance are studied. With the single exception of a circular geometry, the third system is considered identical to the second system. Here, the influence of nongray walls is also studied.

Tiwari, Surendra N.↗

Physics and design of advanced IR bolometers and photoconductors

The state-of-the-art in photoconductors and bolometers is reviewed, centering on the materials development and concepts of extrinsic Si and Ge devices. It is suggested that in the field of 2-D detector arrays the hybrid circuits containing a photoconductor or a photodiode array, bonded with In solder to a switched MOS readout device with up to 60 x 60 elements, will produce superior performance. Other systems discussed include charge-coupled devices, developed for visible-light imaging and adapted to IR imaging, and charge-injection devices, a monolithic combination of a standard photoconductor and a metal-insulator-semiconductor capacitor. Finally, it is pointed out that a nonequilibrium theory of bolometer noise demonstrates the possibility of reducing the traditionally accepted noise limits by 60 percent for Johnson noise and 30 percent for the noise from phonons in the thermal link of the bolometer.

Haller, E. E.↗

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗

Comparison of Coupled Radiative Flow Solutions with Project Fire 2 Flight Data

A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed for use in the design or thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing now and radiation solver and with the Project Fire 2 experimental data. Good agreement has been obtained over the entire Fire 2 trajectory with the experimentally determined values of the stagnation radiation intensity in the 0.2-6.2 eV range and with the total stagnation heating. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multitemperature thermal models, and finite rate chemistry. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.

Olynick, David R.↗

X-ray emission from supernova remnants MSH 14 - 63 and RCW 103

Electron-ion collisions, heavy element enrichment of emissions from reverse-shocked ejecta, abundances in the interstellar medium (ISM), the correspondence between the projected and measured age, ISM particle density, and distance, and classification factors for Type I or II supernovae are examined using a nonequilibrium ionization (NEI) model. The NEI model accounts for the absence of collisional ionization equilibrium and possibly thermal equilibrium between electrons and ions behind the shock fronts, and describes the emission spectra from ejecta and the shocked ISM. When applied to data from the SNR MSH 14 - 63 a non-Coulomb, ion-electron interaction is obtained that is consistent with a shocked ISM. For the SNR RCW 103, however, the processes behind the shock front are not discernable, although the elemental abundances will have solar values if the interactions are non-Coulomb. The age of RCW 103 must be approximated from physical data, as Chinese astronomers would not have been able to see the event from their position on earth two millenia ago.

Nugent, J. J.↗

Comparison of Coupled Radiative Navier-Stokes Flow Solutions with the Project Fire II Flight Data

A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed to use in the design of thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing flow and radiation solver and with the Project Fire II experimental data. Very good agreement has been obtained over the entire Fire II trajectory with the experimentally determined values of the stagnation radiation intensity in the .2 to 6.2 eV range and with the total stagnation heating. The agreement was significantly better than previous numerical predictions. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multi-temperature thermal models, finite-rate chemistry, and a quasi-steady-state or Boltzmann assumption for the calculation of the excited electronic states. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.

Olynick, David R.↗

Modeling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

The Interface Region Imaging Spectrograph has routinely observed the flaring Mg II near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg II in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg I in non-LTE (NLTE) results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg II need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.

Kerr, Graham S.↗

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to nonequilibrium reacting flows is presented. Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.

Suresh, Ambady↗

Nebula Models of Non-Equilibrium Mineralogy: Wark-Lovering Rims

Introduction: The meteorite record contains several examples of minerals that would not persist if allowed to come to equilibrium with a cooling gas of solar composition. This includes all minerals in CAIs and AOAs. Their survival is generally ascribed to physical removal of the object from the gas (isolation into a large parent object, or ejection by a stellar wind), but could also result from outward radial diffusion into cooler regions, which we discuss here. Accretion of CAIs into planetesimals has also been relied on to preserve them against loss into the sun. However, this suggestion faces several objections. Simple outward diffusion in turbulence has recently been modeled in some detail, and can preserve CAIs against loss into the sun [2]. Naturally, outward radial diffusion in turbulence is slower than immediate ejection by a stellar wind, which occurs on an orbital timescale. Here we ask whether these different transport mechanisms can be distinguished by nonequilibrium mineralogy, which provides a sort of clock. Our application here is to one aspect of CAI mineralogy - the Wark-Lovering rims (WLR); even more specifically, to alteration of one layer in the WLR sequence from melilite (Mel) to anorthite (An).

Cuzzi, J. N.↗

Nonlinear evolution of subsonic and supersonic disturbances on a compressible free shear layer

The effects of a nonlinear-nonequilibrium-viscous critical layer on the spatial evolution of subsonic and supersonic instability modes on a compressible free shear layer is considered. It is shown that the instability wave amplitude is governed by an integrodifferential equation with cubic-type nonlinearity. Numerical and asymptotic solutions to this equation show that the amplitude either ends in a singularity at a finite downstream distance or reaches an equilibrium value, depending on the Prandtl number, viscosity law, viscous parameter and a real parameter which is determined by the linear inviscid stability theory. A necessary condition for the existence of the equilibrium solution is derived, and whether or not this condition is met is determined numerically for a wide range of physical parameters including both subsonic and supersonic disturbances. it is found that no equilibrium solution exists for the subsonic modes unless the temperature ratio of the low-to-high-speed streams exceeds a critical value, while equilibrium solutions for the most rapidly growing supersonic mode exist over most of the parameter range examined.

Leib, S. J.↗