Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Is Water a Universal Solvent for Life?

There are strong reasons to believe that the laws, principles and constraints of physics and chemistry are universal. It is much less clear how this universality translates into our understanding of the origins of life. Conventionally, discussions of this topic focus on chemistry that must be sufficiently rich to seed life. Although this is clearly a prerequisite for the emergence of living systems, I propose to focus instead on self-organization of matter into functional structures capable of reproduction, evolution and responding to environmental changes. In biology, most essential functions are largely mediated by noncovalent interactions (interactions that do not involve making or breaking chemical bonds). Forming chemical bonds is only a small part of what living systems do. There are specific implications of this point of view for universality. I will concentrate on one of these implications. Strength of non-covalent interactions must be properly tuned. If they were too weak, the system would exhibit undesired, uncontrolled response to natural fluctuations of physical and chemical parameters. If they were too strong kinetics of biological processes would be slow and energetics costly. This balance, however, is not a natural property of complex chemical systems. Instead, it has to be achieved with the aid of an appropriate solvent for life. In particular, potential solvents for life must be characterized by a high dielectric constant to ensure solubility of polar species and sufficient flexibility of biological structures stabilized by electrostatic interactions. Among these solvents, water exhibits a remarkable trait that it also promotes solvophobic (hydrophobic) interactions between non-polar species, typically manifested by a tendency of these species to aggregate and minimize their contacts with the aqueous solvent. Hydrophobic interactions are responsible, at least in part, for many self-organization phenomena in biological systems, such as the formation of cellular boundary structures or protein folding. Strengths of electrostatic and hydrophobic interactions are similar and can be balanced over a wide range of temperatures, which considerably increases the repertoire of interactions that can be used to modulate biological functions. Some properties of water, e.g. its chemical activity against polymerization reactions, are considered as unfavorable to life. In actuality, this might be a favorable trait because life requires a balance between constructive and destructive processes. For example, molecules synthesized in response to specific conditions must be degraded once these conditions change. Otherwise regulation of biological processes would be virtually impossible. Water might not be the only liquid with favorable properties for supporting life. It has been proposed that formamide, which might be present elsewhere in the universe in sufficient quantities to warrant interest, could be a potential alternative to water for the origin of life. However, this will remain highly hypothetical until it is demonstrated in further studies on its physical, chemical and biological properties it is capable of mediating self-organization of matter and providing proper balance between different types of non-covalent interactions.

Pohorill, Andrew↗
DOE OSTI2026

Mechanistic Insights for Plasma-Catalytic CO 2 Reduction over TiO 2 in a Dielectric Barrier Discharge Reactor

Reaction kinetics experiments coupled with phenomenological kinetic modeling and parameter estimation are used to elicit insights into the mechanism and active sites for the plasma-catalytic dissociation of CO 2 on TiO 2 . Experimental and model insights showed that gas-phase reactions contribute at least two-thirds of the overall product formation at explored conditions; weak temperature dependence, strong sensitivity to specific energy input (SEI), apparent first order in CO 2 , and positive influence of cofed argon (Ar) and oxygen (O 2 ) for the gas-phase contributions all suggest that expected plasma reaction steps such as electron-impact and high-energy collisions are the dominant modes for CO 2 dissociation. The Arrhenius-like expression for gas contributions resulted in a preexponential of 4.40 × 10 –3 s –1 , an E SEI,g of 7.90 × 10 –4 mol/kJ, and an E a,g of 1.00 × 10 –3 J/mol. For surface contributions, the small apparent barrier of 16.3 kJ/mol, relatively weaker dependence on SEI, first-order dependence on CO 2 , and insensitivity to cofed Ar and O 2 all point to CO 2 dissociation on TiO 2 surface facets without vacancies and aided by plasma (leading to vibrationally excited CO 2 and/or a reactive surface with significant surface charge accumulation). The Arrhenius-like expression resulted in a preexponential of 7.81 × 10 –2 s –1 , an E SEI,s of 1.90 × 10 –3 mol/kJ, and an E a,s of 1.63 × 10 4 J/mol. The derived kinetic model further enabled a systematic evaluation of the effect of inputs (plasma power, flow rate, CO 2 inlet concentration, and temperature) to identify process trends and optimal operating conditions.

catalyst↗
DOE OSTI2023

Comprehensive Kinetics on the C 7 H 7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C 7 H 7 potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Near-critical point phenomena in fluids (19-IML-1)

Understanding the effects of gravity is essential if the behavior of fluids is to be predicted in spacecraft and orbital stations, and, more generally, to give a better understanding of the hydrodynamics in these systems. An understanding is sought of the behavior of fluids in space. What should emerge from the International Microgravity Lab (IML-1) mission is a better understanding of the kinetics of growth in off-critical conditions, in both liquid mixtures and pure fluids. This complex phenomenon is the object of intensive study in physics and materials sciences area. It is also expected that the IML-1 flight will procure key results to provide a better understanding of how a pure fluid can be homogenized without gravity induced convections, and to what extent the 'Piston Effect' is effective in thermalizing the compressible fluids.

Beysens, D.↗
DOE OSTI2022

Thermodynamically Driven Formation of Intercalated Cu Carpets from Supported Cu Pyramids on MoS 2

Thermodynamic and kinetic analyses based on our first-principles density functional theory calculations are used to interpret the experimentally observed formation of Cu carpets intercalated under the top layer of a 2H-MoS 2 substrate. Spontaneous Cu transport from Cu pyramids on top of the MoS 2 substrate through surface point defects to the growing Cu carpet is shown to be driven by a slightly lower chemical potential for the Cu carpet. We demonstrate that the competition between a preference for a thicker Cu carpet and the cost of elastic stretching of the top MoS 2 layer results in a selected Cu carpet thickness. Furthermore, we also propose that Cu transport occurs primarily via vacancy-mediated diffusion through constricting point defect portals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Establishing structure-sensitivity of ceria reducibility: real-time observations of surface-hydrogen interactions

The first layer of atoms on an oxide catalyst provides the first sites for adsorption of reactants and the last sites before products or oxygen are desorbed. We employ a unique combination of morphological, structural, and chemical analyses of a model ceria catalyst with different surface terminations under an H 2 environment to unequivocally establish the effect of the last layer of atoms on surface reduction. (111) and (100) terminated epitaxial islands of ceria are simultaneously studied in situ allowing for a direct investigation of the structure–reducibility relationship under identical conditions. Kinetic rate constants of Ce 4+ to Ce 3+ transformation and equilibrium concentrations are extracted for both surface terminations. Unlike the kinetic rate constants, which are practically the same for both types of islands, more pronounced oxygen release, and overall higher reducibility were observed for (100) islands compared to (111) ones. The findings are in agreement with coordination-limited oxygen vacancy formation energies calculated by density functional theory. The results point out the important aspect of surface terminations in redox processes, with particular impact on the catalytic reactions of a variety of catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Overall Traveling-Wave-Tube Efficiency Improved By Optimized Multistage Depressed Collector Design

Depressed Collector Design The microwave traveling wave tube (TWT) is used widely for space communications and high-power airborne transmitting sources. One of the most important features in designing a TWT is overall efficiency. Yet, overall TWT efficiency is strongly dependent on the efficiency of the electron beam collector, particularly for high values of collector efficiency. For these reasons, the NASA Glenn Research Center developed an optimization algorithm based on simulated annealing to quickly design highly efficient multistage depressed collectors (MDC's). Simulated annealing is a strategy for solving highly nonlinear combinatorial optimization problems. Its major advantage over other methods is its ability to avoid becoming trapped in local minima. Simulated annealing is based on an analogy to statistical thermodynamics, specifically the physical process of annealing: heating a material to a temperature that permits many atomic rearrangements and then cooling it carefully and slowly, until it freezes into a strong, minimum-energy crystalline structure. This minimum energy crystal corresponds to the optimal solution of a mathematical optimization problem. The TWT used as a baseline for optimization was the 32-GHz, 10-W, helical TWT developed for the Cassini mission to Saturn. The method of collector analysis and design used was a 2-1/2-dimensional computational procedure that employs two types of codes, a large signal analysis code and an electron trajectory code. The large signal analysis code produces the spatial, energetic, and temporal distributions of the spent beam entering the MDC. An electron trajectory code uses the resultant data to perform the actual collector analysis. The MDC was optimized for maximum MDC efficiency and minimum final kinetic energy of all collected electrons (to reduce heat transfer). The preceding figure shows the geometric and electrical configuration of an optimized collector with an efficiency of 93.8 percent. The results show the improvement in collector efficiency from 89.7 to 93.8 percent, resulting in an increase of three overall efficiency points. In addition, the time to design a highly efficient MDC was reduced from a month to a few days. All work was done in-house at Glenn for the High Rate Data Delivery Program. Future plans include optimizing the MDC and TWT interaction circuit in tandem to further improve overall TWT efficiency.

Vaden, Karl R.↗
NASA NTRS2008

Vertical Jump Height is more Strongly Associated with Velocity and Work Performed Prior to Take-off

Vertical jump (VJ) height is commonly used as a measure of athletic capability in strength and power sports. Although VJ has been shown to be a predictor of athletic performance, it is not clear which kinetic ground reaction force (GRF) variables, such as peak force (PF), peak power (PP), peak velocity (PV), total work (TW) or impulse (Imp) are the best correlates. To determine which kinetic variables (PF, PP, PV, TW, and Imp) best correlate with VJ height. Twenty subjects (14 males, 6 females) performed three maximal countermovement VJs on a force platform (Advanced Mechanical Technology, Inc., Watertown, MA, USA). VJ jump height was calculated as the difference between standing reach and the highest reach point measured using a Vertec. PF, PP, PV, TW, and Imp were calculated using the vertical GRF data sampled at 1000 Hz from the lowest point in the countermovement through the concentric portion until take-off. GRF data were normalized to body mass measured using a standard scale (Detecto, Webb City, MO, USA). Correlation coefficients were computed between each GRF variable and VJ height using a Pearson correlation. VJ height (43.4 plus or minus 9.1 cm) was significantly correlated (p less than 0.001) with PF (998 plus or minus 321 N; r=0.51), PP (1997 plus or minus 772 W; r=0.69), PV (2.66 plus or minus 0.40 m (raised dot) s(sup -1); r=0.85), TW (259 plus or minus 93.0 kJ; r=0.82), and Imp (204 plus or minus 51.1 N(raised dot)s; r=0.67). Although all variables were correlated to VJ height, PV and TW were more strongly correlated to VJ height than PF, PP, and Imp. Therefore, since TW is equal to force times displacement, the relative displacement of the center of mass along with the forces applied during the upward movement of the jump are critical determinants of VJ height. PV and TW are key determinants of VJ height, and therefore successful training programs to increase VJ height should focus on rapid movement (PV) and TW by increasing power over time rather than focusing on PF alone.

Bentley, J. R.↗
NASA NTRS2008

Topological Evolution of a Fast Magnetic Breakout CME in 3-Dimensions

W present the extension of the magnetic breakout model for CME initiation to a fully 3-dimensional, spherical geometry. Given the increased complexity of the dynamic magnetic field interactions in 3-dimensions, we first present a summary of the well known axisymmetric breakout scenario in terms of the topological evolution associated with the various phases of the eruptive process. In this context, we discuss the completely analogous topological evolution during the magnetic breakout CME initiation process in the simplest 3-dimensional multipolar system. We show that an extended bipolar active region embedded in an oppositely directed background dipole field has all the necessary topological features required for magnetic breakout, i.e. a fan separatrix surface between the two distinct flux systems, a pair of spine fieldlines, and a true 3-dimensional coronal null point at their intersection. We then present the results of a numerical MHD simulation of this 3-dimensional system where boundary shearing flows introduce free magnetic energy, eventually leading to a fast magnetic breakout CME. The eruptive flare reconnection facilitates the rapid conversion of this stored free magnetic energy into kinetic energy and the associated acceleration causes the erupting field and plasma structure to reach an asymptotic eruption velocity of greater than or approx. equal to 1100 km/s over an approx.15 minute time period. The simulation results are discussed using the topological insight developed to interpret the various phases of the eruption and the complex, dynamic, and interacting magnetic field structures.

Lynch, B. J.↗
DOE OSTI2024

Broadband Rotational Spectroscopy in Uniform Supersonic Flows: Chirped Pulse/Uniform Flow for Reaction Dynamics and Low Temperature Kinetics

ConspectusThe study of gas-phase chemical reactions at very low temperatures first became possible with the development and implementation of the CRESU (French acronym for Reaction Kinetics in Uniform Supersonic Flows) technique. CRESU relies on a uniform supersonic flow produced by expansion of a gas through a Laval (convergent-divergent) nozzle to produce a wall-less reactor at temperatures from 10 to 200 K and densities of 1016-1018 cm-3 for the study of low temperature kinetics, with particular application to astrochemistry. In recent years, we have combined uniform flows with revolutionary advances in broadband rotational spectroscopy to yield an instrument that affords near-universal detection for novel applications in photodissociation, reaction dynamics, and kinetics. This combination of uniform supersonic flows with chirped-pulse Fourier-transform microwave spectroscopy (Chirped-Pulse/Uniform Flow, CPUF) permits detection of any species with a modest dipole moment, thermalized to the uniform temperature of the gas flow, with isomer, conformer, and vibrational state specificity. In addition, the use of broadband, high-resolution, and time-dependent (microsecond time scale) micro- and mm-wave spectroscopy makes it an ideal tool for characterizing both transient and stable molecules, as well as studying their spectroscopy and dynamics.In this Account, we review recent advances made using the CPUF technique, including studies of photodissociation, radical-radical reaction dynamics, and low temperature kinetics. These studies highlight both the strength of universal and multiplexed detection and the challenges of coupling it to a high-density collisional environment. Product branching and product evolution as a function of time have been measured for astrochemically relevant systems, relying on the detailed characterization of these flow conditions via experiments and fluid dynamics simulations. In the photodissociation of isoxazole, an unusual heterocyclic molecule with a very low-energy conical intersection, we have identified 7 products in 5 reaction channels and determined the product branching, pointing to both direct and indirect pathways. We have also approached the same system from separated NO and C3H3 reactants to explore a broader range of the potential energy surface, demonstrating the power of multichannel branching measurements for complex radical-radical reactions. We determined the product branching in the C3H2 isomers in the photodissociation of the propargyl radical and identified the importance of a hydrogen atom catalyzed isomerization to the lowest energy cyclic form. This then motivated a study of direct D-H exchange reaction in radicals, in which we demonstrate that it is an important and overlooked pathway for deuterium fractionation in astrochemical environments. Recently, we have shown the measurement of low temperature kinetics inside an extended Laval nozzle, after which a shock-free secondary expansion to low temperature and density affords an ideal environment for detection by rotational spectroscopy. These results highlight the power and potential of the CPUF approach, and future prospects will also be discussed in light of these developments.

Dias, Nureshan↗
DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Cleavage at the nsp5–nsp6 site of SARS-CoV-2 main protease intermediate precursor is faster from a monomer than a dimer form

Our previous studies of severe acute respiratory syndrome coronavirus 2 main protease (MPro) precursor monomer indicate that the initial N-terminal nonstructural protein (nsp)4/nsp5 cleavage occurs intramolecularly, with a small fraction of the active site loop equilibrium being in the active state. To understand the influence of dimer formation of MPro upon N-terminal cleavage on the subsequent C-terminal nsp5/nsp6 intermolecular cleavage kinetics, the stepwise processing of a monomeric, inactive precursor containing the native terminal cleavage sites of MPro (MBP- (−6) MPro C145A(+3) -GB1-6H, 86.2 kDa) by mature WT MPro (MPro WT ) was investigated. Differential scanning fluorimetry and analytical ultracentrifugation measurements of various MPro constructs suggest that the C145A mutation decreases the dimer dissociation constant (K dimer ) by ∼26-fold, relative to WT C145 and H41A. The monomeric precursor’s nsp4–nsp5 site appears to saturate MProWT’s active sites and cleave faster, followed by a slower first-order cleavage at the C-terminal site. No detectable product resulting from the C-terminal cleavage is observed until most of the N-terminal cleavage is complete. The initial intermediate product (termed MPro C145A-IP ) is a homodimer with an estimated K dimer of <0.05 μM. In contrast, the first-order kinetics observed for the cleavage of the monomeric form of the intermediate product is at least 300 times faster than that of the dimer form. Room-temperature X-ray structure of the MPro C145A-IP –ensitrelvir complex is like that of the MPro WT –ensitrelvir complex and reveals a dynamic C-terminal region including MPro residues 302 to 306. These results are interpreted from the point of view of a mechanism in which nsp5–nsp6 cleavage may occur from a monomeric intermediate, and dimer formation restricts this cleavage.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Resolving Discrepancies between State-of-the-Art Theory and Experiment for HO 2 + HO 2 via Multiscale Informatics

Recent high-level theoretical calculations predict a mild temperature dependence for HO 2 + HO 2 inconsistent with state-of-the-art experimental determinations that upheld the stronger temperature dependence observed in early experiments. Via MultiScale Informatics analysis of the theoretical and experimental data, we identified an alternative interpretation of the raw experimental data that uses HO 2 + HO 2 rate constants nearly identical to theoretical predictions---implying that the theoretical and experimental data are actually consistent, at least when considering the raw data from experimental studies. Here, similar analyses of typical signals from low-temperature experiments indicate that an HOOOOH intermediate---identified by recent theory but absent from earlier interpretations---yields modest effects that are smaller than, but may have contributed to, the scatter in data among different experiments. More generally, the findings demonstrate that modern chemical theories and experiments have progressed to a point where meaningful comparison requires joint consideration of their data simultaneously.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

SAXS analysis of the intrinsic tenase complex bound to a lipid nanodisc highlights intermolecular contacts between factors VIIIa/IXa

Abstract The intrinsic tenase (Xase) complex, formed by factors (f) VIIIa and fIXa, forms on activated platelet surfaces and catalyzes the activation of factor X to Xa, stimulating thrombin production in the blood coagulation cascade. The structural organization of the membrane-bound Xase complex remains largely unknown, hindering our understanding of the structural underpinnings that guide Xase complex assembly. Here, we aimed to characterize the Xase complex bound to a lipid nanodisc with biolayer interferometry (BLI), Michaelis–Menten kinetics, and small-angle X-ray scattering (SAXS). Using immobilized lipid nanodiscs, we measured binding rates and nanomolar affinities for fVIIIa, fIXa, and the Xase complex. Enzyme kinetic measurements demonstrated the assembly of an active enzyme complex in the presence of lipid nanodiscs. An ab initio molecular envelope of the nanodisc-bound Xase complex allowed us to computationally model fVIIIa and fIXa docked onto a flexible lipid membrane and identify protein–protein interactions. Our results highlight multiple points of contact between fVIIIa and fIXa, including a novel interaction with fIXa at the fVIIIa A1–A3 domain interface. Lastly, we identified hemophilia A/B-related mutations with varying severities at the fVIIIa/fIXa interface that may regulate Xase complex assembly. Together, our results support the use of SAXS as an emergent tool to investigate the membrane-bound Xase complex and illustrate how mutations at the fVIIIa/fIXa dimer interface may disrupt or stabilize the activated enzyme complex.

Hematology↗
NASA NTRS1992

Navier-Stokes solutions with surface catalysis for Martian atmospheric entry

In this study numerical solutions have been obtained for two-dimensional axisymmetric hypersonic nonequilibrium CO2 flow over a high angle blunt cone with appropriate surface boundary conditions to account for energy and mass conservation at the body surface. The flowfield is described by the Navier-Stokes equations and multicomponent conservation laws which account for both translational and internal vibrational nonequilibrium effects. Complete forebody solutions have been obtained for the peak heating point of the Mars entry trajectory specified in the proposed NASA MESUR (Mars Environmental Survey) project. In these solutions, radiative equilibrium wall temperature and surface heating distributions are determined over the MESUR aeroshell forebody for entry velocity equal to 7 km/sec with varying degrees of surface catalysis. The effects of gas kinetics, surface catalysis, transport properties, and vibrational relaxation times on the surface heating are examined. The results identify some important issues in the prediction of surface heating for flows in thermochemical nonequilibrium and show that the Navier-Stokes code used herein is effective for thermal protection system design and materials selection.

Chen, Y.-K.↗
NASA NTRS2000

Preliminary Design Optimization For A Supersonic Turbine For Rocket Propulsion

In this study, we present a method for optimizing, at the preliminary design level, a supersonic turbine for rocket propulsion system application. Single-, two- and three-stage turbines are considered with the number of design variables increasing from 6 to 11 then to 15, in accordance with the number of stages. Due to its global nature and flexibility in handling different types of information, the response surface methodology (RSM) is applied in the present study. A major goal of the present Optimization effort is to balance the desire of maximizing aerodynamic performance and minimizing weight. To ascertain required predictive capability of the RSM, a two-level domain refinement approach has been adopted. The accuracy of the predicted optimal design points based on this strategy is shown to he satisfactory. Our investigation indicates that the efficiency rises quickly from single stage to 2 stages but that the increase is much less pronounced with 3 stages. A 1-stage turbine performs poorly under the engine balance boundary condition. A portion of fluid kinetic energy is lost at the turbine discharge of the 1-stage design due to high stage pressure ratio and high-energy content, mostly hydrogen, of the working fluid. Regarding the optimization technique, issues related to the design of experiments (DOE) has also been investigated. It is demonstrated that the criteria for selecting the data base exhibit significant impact on the efficiency and effectiveness of the construction of the response surface.

Papila, Nilay↗