Structure of the blast wave with ionization
Plane, cylindrical and spherical blast waves structure with ionization at local thermodynamic equilibrium analyzed using successive approximation to nonsimilar solution
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.
Plane, cylindrical and spherical blast waves structure with ionization at local thermodynamic equilibrium analyzed using successive approximation to nonsimilar solution
An analysis for viscous radiating shock layer flow over a blunt body with massive blowing is presented. Detailed equilibrium computations for the thermodynamics and transport properties are included in the solution. Molecular band and atomic line, as well as continuum radiation are included. For a typical Jovian entry problem, a 40% reduction in total radiative flux at stagnation point, due to the injection of carbon gas, was found. The reduction increases to about 70% at about four body nose radii downstream. The increase of flux reduction indicates that most of the radiative energy is absorbed by the injected carbon gas and dumped into the wake. An examination of the spectral distribution of the surface radiative flux reveals that most of the flux reduction occurs in molecular bands of C2.
The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) has been modified to compute viscous equilibrium flow. Periodic calls to the thermodynamic and transport property curve-fits enable solutions to be computed for small percentage increase in computer time when compared with perfect gas times. The code is used to compute the hypersonic flow over slender and blunt cones, and solutions are compared with other computational techniques and flight data.
Examination of the protein crystallization process involves investigation of the liquid and solid state and a protein's properties in these states. Liquid state studies such as protein self association in solution by light scattering methods or other methods have been used to examine a protein Is properties and therefore its crystallization process and conditions. Likewise can osmotic pressure data be used to examine protein properties and various published osmotic pressure studies were examined by us to correlate osmotic pressure to protein solution properties. The solution behavior of serum albumin, alpha - chymotrypsin, beta - lactoglobulin and ovalbumin was examined over a range of temperatures, pH values and different salt types and concentrations. Using virial expansion and a local composition model the non ideal solution behavior in form of the activity coefficients (thermodynamic) was described for the systems. This protein activity coefficient data was related to a protein's solubility behavior and this process and the results will be presented.
The observation that arc magmas are the most oxidized on Earth have led petrologists to question whether the subduction process might cause oxidation of the sub-arc mantle source. A strong correlation between the input of slab-derived aqueous fluid and Fe3(+)/ΣFe in arc magmas has led to the hypothesis that slab fluids may facilitiate the transfer of redox potential from oxidized slab material to mantle wedge and subsequently to primary arc melts. Despite this intuitive link, identifying an efficient and ubiquitous chemical process to transfer oxidation state in slab fluids has been challenging. Pure H2O alone is an inefficient oxidizer in the mantle, necessitating another oxidant within the fluid. Commonly invoked components include S or C, which are likely heterogeneously distributed within subducting slabs, and direct transfer of Fe3(+), which is only possible in very solute-rich fluids. Here we present a new mechanism to explain the oxidation of the sub-arc mantle by slab-derived aqueous fluid, wherein dissolved H(+)(aq) is reduced to H2(aq). Redox equilibrium is satisfied by electron transfer from Fe2(+)(s), which is oxidized in the rock to Fe3(+)(s) in the reaction: 2Fe2(+)(s) + 2H(+)(aq) = 2Fe3(+)(s) + H2(aq). This simple reaction was previously overlooked due to the assumption that all H must be bound by H2O molecules (i.e., no dissolved H), applicable to a pure H2O fluid, but not one with a solute load as addressed here. The application of thermodynamic modeling tools from aqueous geochemistry (DEW, EQ3/6) to mantle petrology predicts this reaction between mantle rock and slab fluid with our own experimentally constrained chemistry and fO2 (QFM+2). Mass transfer model results show an increase in Fe3(+)/ΣFe in mantle rock from MORB- like (0.15) to arc-like (0.2–0.3) values at reasonable fluid fluxes, co-incident with increasing H2 activity in the fluid. This can explain the stark correlation between slab fluid input into the mantle wedge and the observed fO2 of arc magmas. Moreover, it does not require diffusion of hydrogen out of the system. Instead, H2 remains dissolved in fluid in equilibrium with oxidized mantle rock.
A viscous shock layer method has been developed which includes the effects of chemical and thermal nonequilibrium and is coupled with a radiation analysis which includes thermodynamic nonequilibrium effects. This code has been used to obtain solutions for a wide variety of nonequilibrium re-entry conditions in air. The results are tabulated and displayed graphically. Comparisons are made to similar results obtained for radiatively coupled equilibrium flow and conclusions drawn on the effect of nonequilibrium and in particular thermodynamic nonequilibrium on the radiative environment about re-entry vehicles.
This work studies pressurization and pressure control of a large liquid hydrogen storage tank. A finite element model is developed that couples a lumped thermodynamic formulation for the vapor region with a complete solution of the Navier-Stokes and energy equations for the flow and temperature fields in the liquid. Numerical results show that buoyancy effects are strong, even in microgravity, and can reposition a vapor bubble that is initially at the center of the tank to a region near the tank wall in a relatively short time. Long-term tank pressurization with the vapor bubble at the tank wall shows that after an initial transient lasting about a week, the final rate of pressure increase agrees with a purely thermodynamic analysis of the entire tank. However, the final pressure levels are quite different from thermodynamic predictions. Numerical results also show that there is significant thermal stratification in the liquid due to the effects of natural convection. A subcooled jet is used to provide simultaneous cooling and mixing in order to bring the tank pressure back down to its initial value. Three different jet speeds are examined. Although the lowest jet speed is ineffective at controlling the pressure because of insufficient penetration into the liquid region, the highest jet speed is shown to be quite effective at disrupting thermal stratification and reducing the tank pressure in reasonable time.
Coupling the Knudsen effusion method with mass spectrometry has proven to be one of the most useful experimental techniques for studying the equilibrium between condensed phases and complex vapors. The Knudsen effusion method involves placing a condensed sample in a Knudsen cell, a small "enclosure", that is uniformly heated and held until equilibrium is attained between the condensed and vapor phases. The vapor is continuously sampled by effusion through a small orifice in the cell. A molecular beam is formed from the effusing vapor and directed into a mass spectrometer for identification and pressure measurement of the species in the vapor phase. Knudsen cell mass spectrometry (KCMS) has been used for nearly fifty years now and continues to be a leading technique for obtaining thermodynamic data. Indeed, much of the well-established vapor specie data in the JANAF tables has been obtained from this technique. This is due to the extreme versatility of the technique. All classes of materials can be studied and all constituents of the vapor phase can be measured over a wide range of pressures (approximately 10(exp -4) to 10(exp -11) bar) and temperatures (500-2800 K). The ability to selectively measure different vapor species makes KCMS a very powerful tool for the measurement of component activities in metallic and ceramic solutions. Today several groups are applying KCMS to measure thermodynamic functions in multicomponent metallic and ceramic systems. Thermodynamic functions, especially component activities, are extremely important in the development of CALPHAD (Calculation of Phase Diagrams) type thermodynamic descriptions. These descriptions, in turn, are useful for modeling materials processing and predicting reactions such as oxide formation and fiber/matrix interactions. The leading experimental methods for measuring activities are the Galvanic cell or electro-motive force (EMF) technique and the KCMS technique. Each has specific advantages, depending on material and conditions. The EMF technique is suitable for lower temperature measurements, provided a suitable cell can be constructed. KCMS is useful for higher temperature measurements in a system with volatile components. In this paper, we briefly review the KCMS technique and identify the major experimental issues that must be addressed for precise measurements. These issues include temperature measurements, cell material and cell design and absolute pressure calibration. The resolution of these issues are discussed together with some recent examples of measured thermodynamic data.
Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.
Equilibrium thermodynamic and flow properties are presented in tabulated and graphical form for moving, standing, and reflected normal shock waves into helium-hydrogen mixtures representative of proposed outer planet atmospheres. The volumetric compositions of these mixtures are 0.35He-0.65H2, 0.20He-0.80H2, and 0.05He-0.95H2. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular-weight ratio, isentropic exponent, velocity, and species mole fractions. Incident (moving) shock velocities are varied from 4 to 70 km/sec for a range of initial pressure of 5 N/sq m to 100 kN/sq m. The present results are applicable to shock-tube flows and to free-flight conditions for a blunt body at high velocities. A working chart illustrating idealized shock-tube performance with a 0.20He-0.80H2 test gas and heated helium driver gas is also presented.
We measured the infrared (IR) nu(sub CH) absorption bands around 3.4 microns (2800 - 3000 cm(sup -1) in large laboratory-grown magnesium oxide (MgO) and natural olivine single crystals that crystallized from CO/CO2/H2O saturated melts. These bands are very similar to those from many astronomical sources, such as from dust in the diffuse interstellar medium (ISM), from the outflow of evolved stars, etc., and they are characteristic of aliphatic -CH2- and -CH3 entities. In our laboratory single crystals the VCH bands arise from C-H entities that were introduced by a solid solution process, and that are imbedded in the mineral matrix in form of polyatomic C(sub n) entities with C atoms bonded to O and to H. Heating breaks the C-H bonds, causing hydrogen to disperse in the mineral matrix. C-H bonds are re-established rapidly during annealing. We propose that dust grains probably contain the same type of internal C(sub n)-H entities in solid matrix rather than an external organic layer covering the grain surfaces. Thermodynamical arguments show that the concentration of organics in solid solution in small grains can be comparable to that found in astronomical environments.
The results of a theoretical investigation of thermodynamic effects on developed cavitation are presented. An approximate solution to the conservation equations for a two-phase laminar boundary layer is obtained. This analysis produces an expression for the temperature difference between the liquid and vapor phases which can be applied to developed cavity flows. Experimental data of cavity temperature depressions are correlated using this result. In addition, a theoretical estimate of the Nusselt number for the cavity is made using a turbulent boundary layer cavity model proposed by Brennen. The result agrees in part with empirically determined expressions for the cavity Nusselt number.
A numerical cloud model is presented which can describe the evolution of a cloud starting from moist aerosol-laden air through the diffusional growth regime. The model is designed for the direct support of cloud chamber laboratory experimentation, i.e., experiment preparation, real-time control and data analysis. In the model the thermodynamics is uncoupled from the droplet growth processes. Analytic solutions for the cloud droplet growth equations are developed which can be applied in most laboratory situations. The model is applied to a variety of representative experiments.
The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.
A statistical solution for the problem of non-LTE radiative transfer in the millions of atomic transitions responsible for line blanketing is proposed which takes into account the influence of scattering on the atomic populations. A multifrequency/multigray algorithm is used to analyze the thermal equilibrium. Although CO is shown to cool the upper atmosphere to a boundary temperature of 2640 K, the molecule has little influence in the region of the solar temperature minimum.
The durability of gas-turbine engine components can be significantly affected by the ingestion of siliceous particles, which can melt at high temperature and corrode protective coatings that are essential for long life requirements. The silicate debris consists mainly of CaO-MgO-Al2O3-SiO2 (CMAS) and is usually ingested by aircraft engines during and after take-off, sticking to their hot surfaces and resulting in the formation of calcium rare-earth silicate oxyapatites. The thermochemistry of coatings and their reaction products with molten silicate debris are crucial to understand in order to improve the durability of gas-turbine engines. Here we discuss results of high temperature drop solution calorimetry, drop-and-catch calorimetry (DnC) and differential thermal analysis (DTA) techniques for the thermodynamic properties of both thermal barrier coatings (TBCs) and environmental barrier coatings (EBCs) and their reaction with CMAS compositions. The enthalpies of solution of Y2Si2O7, Yb2Si2O7, 31YSZ, and 16RESZ based coatings and the oxyapatite are moderately positive. However, oxyapatite formation is only favorable over coating dissolution in terms of enthalpy for 7YSZ. The enthalpies of mixing between the coatings and the molten silicate are less exothermic for Yb2Si2O7 and CaYb4Si3O13 than for 7YSZ, indicating lower energetic stability of the latter against molten silicate corrosion. We also report for the first time the calorimetric measurements of the enthalpies of formation of rare-earth silicate based EBC coatings and oxyapatites (rare-earth, RE = Y, Yb, Gd, Dy, Er, Nd and Sm).
The International Space Station (ISS) employs an Internal Active Thermal Control System (IATCS) comprised of several single-phase water coolant loops. These coolant loops are distributed throughout the ISS pressurized elements. The primary element coolant loops (i.e., US Laboratory module) contain a fluid accumulator to accommodate thermal expansion of the system. Other element coolant loops are parasitic (i.e., Airlock), have no accumulator, and require an alternative approach to insure that the system Maximum Design Pressure (MDP) is not exceeded during the Launch to Activation phase. During this time the element loop is a stand alone closed individual system. The solution approach for accommodating thermal expansion was affected by interactions of system components and their particular limitations. The mathematical solution approach was challenged by the presence of certain unknown or not readily obtainable physical and thermodynamic characteristics of some system components and processes. The purpose of this paper is to provide a brief description of a few of the solutions that evolved over time, a novel mathematical solution to eliminate some of the unknowns or derive the unknowns experimentally, and the testing and methods undertaken.
The International Space Station (ISS) employs an Internal Active Thermal Control System (IATCS) comprised of several single-phase water coolant loops. These coolant loops are distributed throughout the ISS pressurized elements. The primary element coolant loops (i.e. U.S. Laboratory module) contain a fluid accumulator to accomodate thermal expansion of the system. Other element coolant loops are parasitic (i.e. Airlock), have no accumulator, and require an alternative approach to insure that the system maximum design pressure (MDP) is not exceeded during the Launch to Activation (LTA) phase. During this time the element loops is a stand alone closed system. The solution approach for accomodating thermal expansion was affected by interactions of system components and their particular limitations. The mathematical solution approach was challenged by the presence of certain unknown or not readily obtainable physical and thermodynamic characteristics of some system components and processes. The purpose of this paper is to provide a brief description of a few of the solutions that evolved over time, a novel mathematical solution to eliminate some of the unknowns or derive the unknowns experimentally, and the testing and methods undertaken.