A functional random-walk model of the many- particle system
Generalized random walk method presented as master equation in nonequilibrium statistical mechanics
SEARCH · Engineering Papers
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.
Generalized random walk method presented as master equation in nonequilibrium statistical mechanics
Correlations of orbiter windward surface entry heating data from the first five flights are presented with emphasis on boundary layer transition and the effects of catalytic recombination. Results show that a single roughness boundary layer transition correlation developed for spherical element trips works well for the orbiter tile system. Also, an engineering approach for predicting heating in nonequilibrium flow conditions shows good agreement with the flight test data in the time period of significant heating. The results of these correlations, when used to predict orbiter heating for a high cross mission, indicate that the thermal protection system on the windward surface will perform successfully in such a mission.
Radiation intensity profiles behind incident normal shock waves in pure CO2 have been measured spectroscopically in the Langley Arc-Driven Shock Tube. These profiles, which were obtained for shock velocities between 9 and 13 km/sec and ambient densities corresponding to Venus altitudes between 100 and 80 km, were measured in the vacuum ultraviolet regime. Wavelengths of 127.7, 158.0, 177.5, and 195.0 nm were monitored simultaneously using a four-channel vacuum spectrograph equipped with sodium-salicylate-coated photomultipliers, thereby including the CO(4+) band system which is the most prominent radiator. Measured nonequilibrium overshoots are modeled to provide a means of estimating the effect of nonequilibrium radiation heating to the stagnation region of proposed aero-shells for Venusian entry. These results indicate a significant increase in radiative heating due to nonequilibrium effects. The measurements are believed to represent the most accurate data available on the effect of nonequilibrium radiative heat transfer for Venus entry. This accuracy is primarily due to improved spectrographic instrumentation, which is discussed in some detail regarding its application in related studies.
The kinematic behavior of a dilute, two-dimensional 'gas' of parallel twisted magnetic flux tubes in a highly conducting fluid is considered. Collisions between flux tubes with the same sense of twist are inelastic due to the reconnection of their oppositely directed azimuthal magnetic field components at impact. In some collisions, the tension in the reconnected magnetic flux is sufficient to bind the two flux tubes together. Collisions between flux tubes with the opposite sense of twist are elastic. It is shown that a tenuous gas of individual twisted flux tubes condenses into two large-scale regions of magnetic field with opposite twists. This calculation illustrates the results of Montgomery and co-workers on the migration of magnetic energy towards small wavenumbers in two-dimensional magnetohydrodynamic (MHD) turbulence. In particular, this problem illustrates the dynamical nonequilibrium of the relaxation phase in turbulent MHD systems.
A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Pathfinder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.
A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Path finder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.
A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Pathfinder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.
Thermochemical nonequilibrium flow simulation capabilities have been previously implemented, verified, and validated for central processing unit (CPU) systems in NASA’s unstructured-grid computational fluid dynamics solver FUN3D. Many exascale-class high-performance computing systems will rely on graphics processing unit (GPU) architectures for high throughput and energy efficiency; thus, CPU-based scientific computing software unable to effectively utilize these systems must be updated. In this work, we present a CUDA C++ implementation of FUN3D’s thermochemical nonequilibrium flow simulation capabilities targeting NVIDIA Tesla GPUs. An overview of the porting and optimization strategy is described and performance comparisons with other recent architectures are presented. Scaling to thousands of GPUs is demonstrated, yielding computational performance equivalent to that of several million CPU cores. The implementation enables efficient, high-fidelity, scale-resolving simulations of thermochemical nonequilibrium flows for many applications including atmospheric entry, hypersonics, and combustion.
In recent years there has been considerable interest in intersubband-transition-based infrared semiconductor quantum well (QW) lasers because of their potential applications. In the mid-infrared range, both electrically-injected quantum cascade lasers [1] and optically-pumped multiple QW lasers [2] have been experimentally realized. In these studies, optical gain is due to population inversion between the lasing subbands. It was also proposed that stimulated Raman scattering in QW systems can produce net infrared optical gain [3j. In such a nonlinear optical scheme, the appearance of optical gain that may lead to intersubband Raman lasers does not rely on the population inversion. Since, in tile resonant Raman process (Raman gain is the largest in this case), the pump field induces population redistribution among subbands in the QW s ystem, it seems that a realistic estimate of the optical gain has to include this effect. Perturbative calculations used in the previous work [3] may overestimate the Raman gain. In this paper we present a nonperturbative calculation of terahertz gain of optically-pumped semiconductor step quantum wells. Limiting optical transitions within the conduction band of QW, we solve the pump-field-induced nonequilibrium distribution function for each subband of the QW system from a set of coupled rate equations. Both intrasubband and intersubband relaxation processes in the quantum well system are included. Taking into account the coherent interactions between pump and THz (signal) waves, we we derive the susceptibility of the QW system for the THz field. For a GaAs/AlGaAs step QW, we calculate the Thz gain spectrum for different pump frequencies and intensities. Under moderately strong pumping (approximately 0.3 MW/sq cm), a significant THz gain (approximately 300/m) is predicted. It is also shown that the coherent wave interactions (resonant stimulated Raman processes) contribute significantly to the THz gain.
Terahertz optical gain due to intersubband transitions in optically-pumped semiconductor quantum wells (QW's) is calculated nonperturbatively. We solve the pump- field-induced nonequilibrium distribution function for each subband of the QW system from a set of rate equations that include both intrasubband and intersubband relaxation processes. The gain arising from population inversion and stimulated Raman processes is calculated in a unified manner. We show that the coherent pump and signal wave interactions contribute significantly to the THz gain. Because of the optical Stark effect and pump-induced population redistribution, optical gain saturation at larger pump intensities is predicted.
An upwind-biased point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional viscous hypersonic flows in chemical and thermal nonequilibrium is described. Details of the algorithm development, in the context of an 11-species two-temperature reacting gas model, are emphasized. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large block-tridiagonal systems. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km) over the Aeroassist Flight Experiment configuration, obtained on a multidomain grid, are discussed.
This study deals with a systematic investigation of nonequilibrium processes in supersonic combustion. The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with nonequilibrium chemistry and thermodynamics, coupled with radiation, for hydrogen-air systems. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. For a basic understanding of the flow physics, premixed flows undergoing finite rate chemical reactions are investigated. Results obtained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species, and that this can have a noticeable influence on the flowfield. The second part of this study deals with premixed reacting flows under thermal nonequilibrium conditions. Here, the critical problem is coupling of the vibrational relaxation process with the radiative heat transfer. The specific problem considered is a premixed expanding flow in a supersonic nozzle. Results indicate the presence of nonequilibrium conditions in the expansion region of the nozzle. This results in reduction of the radiative interactions in the flowfield. Next, the present study focuses on investigation of non-premixed flows under chemical nonequilibrium conditions. In this case, the main problem is the coupled turbulence-chemistry interaction. The resulting formulation is validated by comparison with experimental data on reacting supersonic coflowing jets. Results indicate that the effect of heat release is to lower the turbulent shear stress and the mean density. The last part of this study proposes a new theoretical formulation for the coupled turbulence-radiation interactions. Results obtained for the coflowing jets experiment indicate that the effect of turbulence is to enhance the radiative interactions.
The interaction between the solar wind and the interstellar medium is modeled self-consistently using numerical solutions of the time-dependent gasdynamic equations in spherical and cylindrical coordinates. For the results presented here it is assumed that the solar system moves through the surrounding medium with a supersonic velocity. After an initial (nonequilibrium) state has been specified, the numerical solution follows the evolution in time until the interaction relaxes to a dynamic equilibrium. As would be expected, the solutions show the formation of a bow shock upstream of the traveling solar system to deflect the interstellar plasma around the cavity created by the solar wind. A terminiation shock also forms to slow and compress the solar wind plasma. For the simulation in spherical coordinates, the downstream portion of the termination shock reaches equilibrium more than three times further from the Sun than the equilibrium distance to the termination shock on the upstream side.
It has been suggested previously that the performance of scramjet propulsion system may be improved by the use of magnetohydrodynamic (MHD) energy bypass: an MHD generator could be made to decelerate the flow entering the combustor, thereby improving combustion efficiency, and the electrical power generated could be made to accelerate the flow exiting from the combustor prior to expanding through the nozzle. In one of such proposed schemes, the MHD generator is proposed to be operated at a low temperature and ionization is to be achieved under nonequilibrium by the application of an external power. In the present work, the required power of such an external source is calculated assuming a 100%-efficient nonequilibrium ionization scheme. The power required is that needed to prevent the degree of ionization from reaching equilibrium with the low gas temperature. The flow is seeded with potassium or cesium. Specific impulse is calculated with and without turbulent friction. The results show that, for typical intended flight conditions, the specific impulse obtained is substantially higher than that of a typical scramjet, but the required external-power is several times that of the power generated in the MHD generator.
The extension of Van Leer's preconditioning procedure to generalized finite-rate chemistry is discussed. Application to viscous flow is begun with the proper preconditioning matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic flow speeds typical of complete propulsion-system simulations. The extended preconditioning matrix necessarily accounts for both thermal and chemical nonequilibrium. Numerical analysis reveals the possible theoretical improvements from using a preconditioner for all Mach number regimes. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number areas. Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning technique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are scaled so that all waves propagate at the same rate, an essential property to eliminate inherent eigenvalue stiffness. This session discusses a synthesis of the thermochemical nonequilibrium flux-splitting developed by Grossman and Cinnella and the characteristic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and three-dimensional flows with generalized finite-rate chemistry. For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, the preconditioning matrix extended to generalized finite-rate chemistry must accommodate a flexible system of moving waves. Fortunately, no new kind of wave appears in the system. The only existing waves are entropy and vorticity waves, which move with the fluid, and acoustic waves, which propagate in Mach number dependent directions. The nonequilibrium vibrational energies and species densities in the unknown state vector act strictly as convective waves. The essential concept for extending the preconditioning to generalized chemistry models is determining the differential variables which symmetrize the flux Jacobians. The extension is then straight-forward. This algorithm research effort will be released in a future version of the production level computational code coined the General Aerodynamic Simulation Program (GASP), developed by Walters, Slack, and McGrory.
An examination is conducted of the basic assumptions of conventional wave-CISK theory, against observational data on mesoscale-organized convective systems. In response to the inadequacy of the quasi-equilibrium assumption, a nonequilibrium model of clouds actively responding to mesoscale low level mass fluxes in a preexisting and conditionally unstable atmosphere is formulated. The structures obtained for the propagating and growing waves are noted to be comparable with observations of midlatitude mesoscale convective systems. An examination is made of the energetics of this modified wave-CISK model, in order to deepen understanding of its instability and scale selection.
In semiconductor quantum well structures, the intersubband energy separation can be adjusted to the terahertz (THz) frequency range by changing the well width and material combinations. The electronic and optical properties of these nanostructures can also be controlled by an applied dc electric field. These unique features lead to a large frequency tunability of the quantum well devices. In the on-going project of modeling of the THz lasers, we investigate the possibility of using optical pumping to generate THz radiation based on intersubband transitions in semiconductor quantum wells. We choose the optical pumping because in the electric current injection it is difficult to realize population inversion in the THz frequency range due to the small intersubband separation (4-40 meV). We considered both small conduction band offset (GaAs/AlGaAs) and large band offset (InGaAs/AlAsSb) quantum well structures. For GaAs/AlGaAs quantum wells, mid-infrared C02 lasers are used as pumping sources. For InGaAs/AlAsSb quantum wells, the resonant intersubband transitions can be excited by the near-infrared diode lasers. For three- and four-subband quantum wells, we solve the pumpfield-induced nonequilibrium distribution function for each subband of the quantum well system from a set of rate equations that include both intrasubband and intersubband relaxation processes. Taking into account the coherent interactions between pump and THz (signal) waves, we calculate the optical gain for the THz field. The gain arising from population inversion and stimulated Raman processes is calculated in a unified manner. A graph shows the calculated THz gain spectra for three-subband GaAs/AlGaAs quantum wells. We see that the coherent pump and signal wave interactions contribute significantly to the gain. The pump intensity dependence of the THz gain is also studied. The calculated results are shown. Because of the optical Stark effect and pump-induced population redistribution, the maximum THz gain saturates at larger pump intensities.
The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with nonequilibrium chemistry and thermodynamics, coupled with radiation, for hydrogen-air systems. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed, expanding flow in a supersonic nozzle. The reacting flow consists of seven species, one of which is the inert N2 molecule. The thermal state of the gas is modeled with one translational-rotational temperature and five vibrational temperatures. The harmonic oscillator model is used in the formulation for vibrational relaxation. The tangent slab approximation is used in the radiative flux formulation. A pseudo-gray model is used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate the presence of nonequilibrium in the expansion region. This reduces the radiative interactions and can have a significant influence on the flowfield.