Flow properties for the continuum viewpoint of a non-equilibrium gas-particle mixture.
Flow properties for nonequilibrium two phase flow of gas particle mixture from continuum mechanics approach
SEARCH · Search NASA
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.
Flow properties for nonequilibrium two phase flow of gas particle mixture from continuum mechanics approach
Silicon ingot growth technologies like the Ubiquitous Crystallization Process (UCP) are solidified within a shaping crucible. The rate at which heat can be lost from this crucible minus the rate at which heat is input from an external source determines the rate at which crystallization will occur. Occasionally, when the process parameters for solidification are exceeded, the normally large multi-centimeter grain size material assocated with the UCP will break down into regions containing extremely small, millimeter or less, grain size material. Accompanying this breakdown in grain growth is the development of so called sinuous grain boundaries. The breakdown in grain growth which results in this type of small grain structure with sinuous boundaries is usually associated with the rapid crystallization that would accompany a system failure. This suggests that there are limits to the growth velocity that one can obtain and still expect to produce material that would possess good photovoltaic properties. It is the purpose to determine the causes behind the breakdown of this material and what parameters will determine the best rates of solidification.
Transport properties for weakly ionized gases, binary Maxwell molecular gases, and fully ionized gases
HgCdTe is the current material of choice for high performance imagers operating at relatively high temperatures. Its lack of technological maturity compared with silicon and wide-band gap III-V compounds is more than offset by its outstanding IR sensitivity and by the relatively benign effect of its materials defects. This latter property has allowed non-equilibrium growth techniques, metal oxide chemical vapor deposition (MOCVD) and molecular beam epitaxy (MBE), to produce device quality long wavelength infrared (LWIR) HgCdTe even on common substrates like GaAs and GaAs/Si. Detector performance in these exotic materials structures is comparable in many ways with devices in equilibrium-grown material. Lifetimes are similar. RoA values at 77K as high as several hundred have been seen in HgCdTe/GaAs/Si with 9.5 micron cut-off wavelength. HgCdTe/GaAs layers with approx. 15 micron cut-off wavelengths have given average 77K RoAs of greater than 2. Hybrid focal plane arrays have been evaluated with excellent operability.
Selected mechanical properties of (plus or minus 45 degree sub 4s) graphite/epoxy composites were found to be affected by sub T sub g annealing. Postcured specimens of Thornel 300 graphite/Narmco 5208 epoxy were sub T sub G annealed at 413 K (140 C) for ca. 10 to the first through 10 to the fifth powers min., with a prior quenching from above T sub g. The ultimate tensile strength, strain-to-break, and toughness of the composite were found to decrease as a function of sub T sub g annealing time. The time-dependent change in properties can be explained on the basis of physical aging which is related to free volume changes in the non-equilibrium glassy state of network epoxies. The results imply possible changes in composite properties with service time.
Noble gas mixture large shock tunnel driven linear MHD generator, examining electron density, operating characteristics and electrical properties
One of the problems associated with the development of high-enthalpy hypersonic test facilities hos been the determination of nozzle flow properties. This note presents a technique for determining the flow properties of air in an equilibrium, non-equilibrium, or frozen state using two test section measurements. A knowledge of two stagnation properties is required before the technique can be applied. The stagnation pressure can be measured, and reference 1 offers a method for determining stagnation enthalpy. Results of the calculations using the method of reference 1 are presented in reference 2 in a chart (chart 21) which can be readily used to determine the stagnation enthalpy from measurements of stagnation pressure, mass flow of air through the tunnel, and throat diameter.
Many of the emerging applications of polymers on space vehicles require materials with outstanding thermal stability. These polymers must also be readily processable in order to facilitate their use. The syntheses and polymerization of a cardo dianhydride were investigated. This monomer was prepared via the reaction of N-methyl 4-nitrophthalimide with a cardo diol. Polyimides containing oxyalkylene linkages were studied. The effects of two additional structural modifications on the polymers' properties were investigated. The effects of carrying out the preparation of poly(amic acid)s under non-equilibrium conditions were examined. Approaches that were investigated included the in-situ neutralization of the generated amic acid and its in-situ esterification.
The analyses, designs, fabrication, and cold-flow acceptance testing of LOX/RP-1 preburner components required for a high-pressure staged-combustion rocket engine are discussed. Separate designs of injectors, combustion chambers, turbine simulators, and hot-gas mixing devices are provided for fuel-rich and oxidizer-rich operation. The fuel-rich design addresses the problem of non-equilibrium LOX/RP-1 combustion. The development and use of a pseudo-kinetic combustion model for predicting operating efficiency, physical properties of the combustion products, and the potential for generating solid carbon is presented. The oxygen-rich design addresses the design criteria for the prevention of metal ignition. This is accomplished by the selection of materials and the generation of well-mixed gases. The combining of unique propellant injector element designs with secondary mixing devices is predicted to be the best approach.
The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.
The magnetic field in the convection zone and photosphere of the sun exists mostly as concentrated tubes of magnetic flux. It is, therefore, necessary to study the basic properties of magnetic flux tubes to obtain a basis for understanding the behavior of the sun's magnetic field. The present investigation is concerned with the global equilibrium shape of a flux tube in the stratified solar atmosphere. A fundamental property of isolated flux tubes is magnetic buoyancy. Attention is given to flux tubes with external field, and twisted flux tubes. It is shown that the analysis of Parker (1975, 1979) and Spruit (1981) for calculating the equilibrium of a slender flux tube in a stratified atmosphere may be extended to more general situations. The slender tube approximation provides a method of solving the problem of modeling the overall curvature of flux tubes. It is found that for a twisted flux tube, there can be two possible equilibrium values of the height.
Hypersonic nonequilibrium flow and its underlying thermodynamic relations, including real gas effects
Kinetic effects, for example nucleation constraints and slow reactions, should have been important in nebular condensation. Consideration of these effects leads to the prediction of pressure-dependent compositions and physical properties of nebular condensates which is consistent with (1) the differences between different classes of chondritic meteorites, (2) some of the differences between planets, and (3) the presence of oxidized iron on the moon and in the eucrite parent body (presumably an asteroid) despite the low abundance of volatiles. Diffusion effects appear to be important for understanding oxygen isotope anomalies in refractory inclusions in Allende. The consideration of kinetic effects leads to more information concerning nebular processes than if equilibrium is assumed.
A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.
The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.
Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.
To select a target comet for a Comet Nucleus Sample Return Mission (CNSRM) it is necessary to have an experimental data base to evaluate the extent of diversity and similarity of comets. For example, the physical properties (e.g., low density) of chondritic porous (CP) interplanetary dust particles (IDPs) are believed to resemble these properties of cometary dust although it is yet to be demonstrated that the porous structure of CP IDPs is inherent to presolar dust particles stored in comet nuclei. Porous structures of IDPs could conceivably form during sublimation at the surface of active comet nuclei. Porous structures are also obtained during annealing of amorphous Mg-SiO smokes which initially forms porous aggregates of olivine + platey tridymite and which, upon continued annealing, react to fluffy enstatite aggregates. It is therefore uncertain that CP IDPs are entirely composed of unmetamorphosed presolar dust. Conceivably, new minerals and textures may form in situ in nuclei of active comets as a function of their individual thermal history. Unmetamorphosed comet dust is probably structurally amorphous. Thermal annealing of this dust can produce ultra fine-grained minerals and this ultrafine grain size of CP IDPs should be considered in assessments of aqueous alterations that could affect presolar dust in comet nuclei between 200 and 400 K. Devitrification and hydration may occur in situ in ice-dust mixtures and the mantle of active comet nuclei. Devitrification, or uncontrolled crystallization, of amorphous precursor dust can produce a range of chemical compositions of ultrafine-grained minerals and (non-equilibrium) mineral assemblages and textures in dust contained in comet nuclei as a function of period and trajectory of orbit and number of perihelion passages (not considering internal heating). Thus, experimental data on relevant processes and reaction rates between 200 and 400 K are needed in order to evaluate comet selection, penetration depth for sampling device and curation of samples for CNSRM.
Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.