Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equilibrium programming”

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 271 records · Page 15

coh3

CoH3 (CoH ver.3) is an optical model, exciton pre-equilibrium, and Hauser-Feshbach statistical model code, which calculates nuclear reaction cross sections for medium to heavy targets in the keV to MeV energy region. This program is written in standard C++, divided into approximately 200 source and header files. CoH solves the Schroedinger equation for optical potentials defined in the code, and calculates differential elastic scattering, reaction, and total cross sections, for neutron, proton, deuteron, triton, 3He, and alpha-particle. Deformed optical potentials are solved with the coupled-channels method, in which the ground state rotational band members, or vibrational phonon states are coupled. The optical model gives particle transmission coefficients that are fed into the statistical model calculations. CoH includes the pre-equilibrium model (exciton model), the direct/semidirect capture model, and the multi-stage Hauser-Feshbach statistical decay with width fluctuation correction based on the Gaussian orthogonal ensemble. For weakly coupled levels, the DWBA (distorted wave Born approximation) method is used to calculate the direct inelastic scattering process to the excited states.

Kawano, Toshihiko↗

CHARMM-GUI Nanomaterial Modeler for Modeling and Simulation of Nanomaterial Systems

Molecular modeling and simulation are invaluable tools for nanoscience that predict mechanical, physicochemical, and thermodynamic properties of nanomaterials and provide molecular-level insight into underlying mechanisms. However, building nanomaterial-containing systems remains challenging due to the lack of reliable and integrated cyberinfrastructures. Here we present Nanomaterial Modeler in CHARMM-GUI, a web-based cyberinfrastructure that provides an automated process to generate various nanomaterial models, associated topologies, and configuration files to perform state-of- the-art molecular dynamics simulations using most simulation packages. The nanomaterial models are based on the interface force field, one of the most reliable force fields (FFs). The transferability of nanomaterial models among the simulation programs was assessed by single-point energy calculations, which yielded 0.01% relative absolute energy differences for various surface models and equilibrium nanoparticle shapes. Three widely used Lennard-Jones (LJ) cutoff methods are employed to evaluate the compatibility of nanomaterial models with respect to conventional biomolecular FFs: simple truncation at r = 12 Å (12 cutoff), force-based switching over 10 to 12 Å (10-12 fsw), and LJ particle mesh Ewald with no cutoff (LJPME). The FF parameters with these LJ cutoff methods are extensively validated by reproducing structural, interfacial, and mechanical properties. As a result, we find that the computed density and surface energies are in good agreement with reported experimental results, although the simulation results increase in the following order: 10-12 fsw <12 cutoff < LJPME. Nanomaterials in which LJ interactions are a major component show relatively higher deviations (up to 4% in density and 8% in surface energy differences) compared with the experiment. Nanomaterial Modeler's capability is also demonstrated by generating complex systems of nanomaterial-biomolecule and nanomaterial-polymer interfaces with a combination of existing CHARMM-GUI modules. We hope that Nanomaterial Modeler can be used to carry out innovative nanomaterial modeling and simulations to acquire insight into the structure, dynamics, and underlying mechanisms of complex nanomaterial-containing systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of two-temperature kinetic model for dissociating and weakly-ionizing nitrogen

The validity of the author's two-temperature, chemical/kinetic model which the author has recently improved is assessed by comparing the calculated results with the existing experimental data for nitrogen in the dissociating and weakly ionizing regime produced behind a normal shock wave. The computer program Shock Tube Radiation Program (STRAP) based on the two-temperature model is used in calculating the flow properties behind the shock wave and the Nonequilibrium Air Radiation (NEQAIR) program, in determining the radiative characteristics of the flow. Both programs were developed earlier. Comparison is made between the calculated and the existing shock tube data on (1) spectra in the equilibrium region, (2) rotational temperature of the N2(+) B state, (3) vibrational temperature of the N2(+) B state, (4) electronic excitation temperature of the N2 B state, (5) the shape of time-variation of radiation intensities, (6) the times to reach the peak in radiation intensity and equilibrium, and (7) the ratio of nonequilibrium to equilibrium radiative heat fluxes. Good agreement is seen between the experimental data and the present calculation except for the vibrational temperature. A possible reason for the discrepancy is given.

Park, C.↗

The chemical effects of the Martian environment on power system component materials: A theoretical approach

In the foreseeable future, an expedition may be undertaken to explore the planet Mars. Some of the power source options being considered for such a mission are photovoltaics, regenerative fuel cells and nuclear reactors. In addition to electrical power requirements, environmental conditions en route to Mars, in the planetary orbit and on the Martian surface must be simulated and studied in order to anticipate and solve potential problems. Space power systems components such as photovoltaic arrays, radiators, and solar concentrators may be vulnerable to degradation in the Martian environment. Natural characteristics of Mars which may pose a threat to surface power systems include high velocity winds, dust, ultraviolet radiation, large daily variation in temperature, reaction to components of the soil, atmosphere and atmospheric condensates as well as synergistic combinations. Most of the current knowledge of the characteristics of the Martian atmosphere and soil composition was obtained from the Viking 1 and 2 missions in 1976. A theoretical study is presented which was used to assess the effects of the Martian atmospheric conditions on the power systems components. A computer program written at NASA-Lewis for combustion research that uses a free energy minimization technique was used to calculate chemical equilibrium for assigned thermodynamic states of temperature and pressure. The power system component materials selected for this study include: silicon dioxide, silicon, carbon, copper, and titanium. Combinations of environments and materials considered include: (1) Mars atmosphere with power surface material, (2) Mars atmosphere and dust component with power surface material, and (3) Mars atmosphere and hydrogen peroxide or superoxide or superoxide with power system material. The chemical equilibrium calculations were performed at a composition ratio (oxidant to reactant) of 100. The temperature for the silicon dioxide material and silicon, which simulate photovoltaic cells, were 300 and 400 K; for carbon, copper and titanium, which simulate radiator surfaces, 300, 500, and 1000 K. All of the systems were evaluated at pressures of 700, 800, and 900 Pa, which stimulate the Martian atmosphere.

Perez-Davis, Marla E.↗

Gas-phase measurements of combustion interaction with materials for radiation-cooled chambers

Foil samples of Ir and Pt are exposed to combustion products in a controlled premixed environment at atmospheric pressure. Electrical heating of the foil samples is used to control the surface temperature and to elevate it above the radiative equilibrium temperature within the test apparatus. Profiles of temperature and OH concentration in the boundary layer adjacent to the specimen surface are measured by laser-induced fluorescence. Measured OH concentrations are significantly higher than equilibrium concentrations calculated for the known mixture ratio and the measured temperature profiles. This result indicates that superequilibrium concentrations of H-atoms and O-atoms are also present in the boundary layer, due to partial equilibrium of the rapid binary reactions of the H2/O2 chemical kinetic system. These experiments are conducted as part of a research program to investigate fundamental aspects of the interaction of combustion gases with advanced high-temperature materials for radiation-cooled thrusters.

Barlow, R. S.↗

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The basic objective of this research is to extend the capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. In the efforts related to LES, we were primarily involved with assessing the performance of the various modern methods based on the Probability Density Function (PDF) methods for providing closures for treating the subgrid fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we concentrated on understanding some of the relevant physics of compressible reacting flows by means of statistical analysis of the data generated by DNS of such flows. In the research conducted in the second year of this program, our efforts focused on the modeling of homogeneous compressible turbulent flows by PDF methods, and on DNS of non-equilibrium reacting high speed mixing layers. Some preliminary work is also in progress on PDF modeling of shear flows, and also on LES of such flows.

Givi, P.↗

COSMIC monthly progress report

Activities of the Computer Software Management and Information Center (COSMIC) are summarized for the month of May 1994. Tables showing the current inventory of programs available from COSMIC are presented and program processing and evaluation activities are summarized. Nine articles were prepared for publication in the NASA Tech Brief Journal. These articles (included in this report) describe the following software items: (1) WFI - Windowing System for Test and Simulation; (2) HZETRN - A Free Space Radiation Transport and Shielding Program; (3) COMGEN-BEM - Composite Model Generation-Boundary Element Method; (4) IDDS - Interactive Data Display System; (5) CET93/PC - Chemical Equilibrium with Transport Properties, 1993; (6) SDVIC - Sub-pixel Digital Video Image Correlation; (7) TRASYS - Thermal Radiation Analyzer System (HP9000 Series 700/800 Version without NASADIG); (8) NASADIG - NASA Device Independent Graphics Library, Version 6.0 (VAX VMS Version); and (9) NASADIG - NASA Device Independent Graphics Library, Version 6.0 (UNIX Version). Activities in the areas of marketing, customer service, benefits identification, maintenance and support, and dissemination are also described along with a budget summary.

Source record↗

Analytical simulation of meander morphology from equilibrium to long-term evolution: Impacts of channel geometry and vegetation-induced coarsening

This study introduces an innovative approach to modeling meandering river morphology, integrating and investigating the effects of geometric characteristics and vegetation-induced channel coarsening. The developed comprehensive framework combines several advanced techniques: Genetic Programming for refining the scour factor of transverse bed slope, a Leaf Area Index (LAI)-enhanced analytical model for quantifying vegetative flow resistance, and an upstream-weighted moving average method for efficient approximation of the convolution integral in meander migration calculations. The model is validated against both an idealized Kinoshita meander and a natural bend of the Tumen River (China) in equilibrium, demonstrating its robustness across diverse scales and conditions. The model's ability to simulate the long-term evolution, including cutoff events, provides valuable insight for river management strategies. The current findings demonstrate that channel geometry, particularly width-to-depth ratio, plays a dominant role in meander evolution, with wider channels prone to more complex and rapid morphological changes. Vegetation effects are most pronounced in channels with moderate width-to-depth ratios, where they can significantly influence migration rates and bed topography. A combination of channel widening and deepening, coupled with strategic vegetation management, can effectively enhance navigability while maintaining channel stability in the studied Tumen River reach. Sensitivity analyses highlight the complex interplay between hydraulic conditions, sediment characteristics, and vegetation in shaping river morphology. This research advances understanding of the multifaceted nature of meandering river systems and offers practical tools for informed decision-making in river engineering and environmental management, particularly in the context of climate change and increasing anthropogenic pressures on fluvial ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Purification and Concentration of 99 Mo from a UREX raffinate

The medical isotope 99m Tc ( t ½ = 6.0 h) accounts for over 80% of isotopes used in diagnostic nuclear medicine today and is normally derived from its transient equilibrium parent 99 Mo ( t ½ = 66 h). To reduce U.S. dependence on 99 Mo derived from HEU by foreign suppliers, the NNSA/M3 program—under the American Medical Isotope Production Act—has been tasked with facilitating the work of domestic 99 Mo suppliers that do not utilize HEU. Superconducting electron linear accelerators with high-Z converter targets can generate bremsstrahlung photons and neutron fluxes that can induce photonuclear reactions and LEU fission. After sufficient production intervals, targets can be rotated out and processed while another batch is irradiated. The process flow chemistry is shown in Figure 1. The irradiated U 3 O 8 targets are retrieved and dissolved in HNO 3 ; the volatile fission products are expelled and captured during this step. The HNO 3 liquor bearing the uranyl, fission and activation products is injected into the UREX suite of liquid-liquid extraction banks. The uranyl is partitioned in typical fashion using tri-$\textit{n}$-butyl phosphate (TBP) in a hydrocarbon diluent. The raffinate—comprising mainly 99 Mo and fission products—is then injected into the MoLLE (Molybdenum Liquid-Liquid Extraction) flowsheet, where the Mo is selectively extracted by an organophosphorous acid extractant such as di-(2-ethylhexyl phosphoric acid) (HDEHP). Trace amounts of Nb, Zr, Np, I, and Te are co-extracted. The Mo is stripped using acetohydroxamic acid (AHA) and fed onto an anion exchange column. Following a series of hydroxide, HCl, and oxalic acid wash steps, the Mo can be recovered in NaOH/NaCl. A manuscript describing this process chemistry in more detail was recently published.

07 ISOTOPE AND RADIATION SOURCES↗

Stagnation region heating of a phenolic-nylon ablator during return from planetary missions

The results of a stagnation-line analysis of the radiative heating of a phenolic-nylon ablator are reported. The analysis includes flowfield coupling with the ablator surface, equilibrium chemistry, multicomponent or binary diffusion and a coupled line and continuum radiation calculations. Users manual and operating instructions for the computer programs listed in the report are included.

Esch, D. D.↗

Program and charts for determining shock tube, and expansion tunnel flow quantities for real air

A computer program in FORTRAN 4 language was written to determine shock tube, expansion tube, and expansion tunnel flow quantities for real-air test gas. This program permits, as input data, a number of possible combinations of flow quantities generally measured during a test. The versatility of the program is enhanced by the inclusion of such effects as a standing or totally reflected shock at the secondary diaphragm, thermochemical-equilibrium flow expansion and frozen flow expansion for the expansion tube and expansion tunnel, attenuation of the flow in traversing the acceleration section of the expansion tube, real air as the acceleration gas, and the effect of wall boundary layer on the acceleration section air flow. Charts which provide a rapid estimation of expansion tube performance prior to a test are included.

Miller, C. G., III↗

Application of computational fluid dynamics to the design of the film cooled STME subscale nozzle for the National Launch System

The status of computational fluid dynamics (CFD) calculations for the Space Transportation Main Engine (STME) film/dump cooled nozzle design is presented, with an emphasis on the timely impact of CFD on the design of the sub-scale nozzle coolant system. The following aspects of the sub-scale coolant delivery system were analyzed with CFD: 1) a design trade study of a mechanical flow splitting device for uniform distribution of the subsonic cavity flow, 2) a design trade study of the subsonic cavity lip to achieve film integrity, and 3) an analysis of the primary flow interaction with the core/secondary coolant streams. All design calculations were performed with the Generalized Aerodynamic Simulation Program (GASP), a 3-D, multi-block, generalized Navier-Stokes code capable of solving with frozen, finite-rate or equilibrium chemical kinetics. The initial design of the subsonic cavity flow used square posts to distribute the sonic orifice jets into a uniform flow. Calculations for this design indicated that an unacceptable mal-distribution of film occurred. Design modifications involving curved and slotted posts were computed in an effort to uniformly distribute the secondary coolant flow. Analysis of these configurations showed that although the flowfield improved in uniformity, it was still unacceptable, especially at higher feed pressures. Results from these studies were then incorporated into a design that resulted in the insertion of a porous metal ring into the subsonic cavity. Subsequent water flow model studies showed that this concept was successful in uniformly distributing flow exiting the cavity. In addition to the design of the subsonic cavity, CFD was also used to analyze the secondary coolant lip and the primary flow interaction with the core/secondary coolant streams. A series of calculations were first performed to modify the subsonic cavity lip contour. The flow over the modified lip was then computed simultaneously with the primary injectors to determine the impact of the subsonic coolant stream on the primary slot jets. Pressure, temperature, velocity, and coolant mass fraction contours are presented for these configurations. Information is given in viewgraph form.

Garrett, Joseph L.↗

Thermodynamic and Transport Property Correlation Formulas for Equilibrium Air from 1,000 Degrees Kelvin to 15,000 Degrees Kelvin

The thermodynamic properties, density and temperature, as well as transport property parameters involving viscosity, Prandtl number (including diffusion effects), and gaseous radiation absorption coefficients have been correlated as a function of enthalpy at four pressure levels (10 (sup -1), 10 (sup 0), 10, and 10 (sup 2) atmospheres). The correlation formulas are written in a generalized form for which coefficients for a particular property and pressure level are tabulated. The correlation formulas are useful in digital computer programs for non-adiabatic viscous flow problems.

AIR↗

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗

A special purpose symbolic manipulation program for the derivation of the equations of motion for large flexible structures

This paper outlines a computer program especially tailored to the task of deriving explicit equations of motion for structures with point-connected substructures. The special purpose program is written in FORTRAN and is designed for performing the specific algebraic operations encountered in the derivation of explicit equations of motion. The derivation is by the Lagrangian approach. Using an orderly kinematical procedure and a discretization and/or truncation scheme, it is possible to write the kinetic and potential energy of each substructure in a compact vector-matrix form. Then, if each element of the matrices and vectors encountered in the kinetic and potential energy is a known algebraic expression, the computer program performs the necessary operations to evaluate the kinetic and potential energy of the system explicitly. Lagrange's equations for small motions about equilibrium can be deduced directly from the explicit form of the system kinetic and potential energy.

Hale, A. L.↗

Derivation of the equations of motion for complex structures by symbolic manipulation

This paper outlines a computer program especially tailored to the task of deriving explicit equations of motion for structures with point-connected substructures. The special purpose program is written in FORTRAN and is designed for performing the specific algebraic operations encountered in the derivation of explicit equations of motion. The derivation is by the Lagrangian approach. Using an orderly kinematical procedure and a discretization and/or truncation scheme, it is possible to write the kinetic and potential energy of each substructure in a compact vector-matrix form. Then, if each element of the matrices and vectors encountered in the kinetic and potential energy is a known algebraic expression, the computer program performs the necessary operations to evaluate the kinetic and potential energy of the system explicitly. Lagrange's equations for small motions about equilibrium can be deduced directly from the explicit form of the system kinetic and potential energy.

Hale, A. L.↗

Computing Viscoplastic Behavior Of A Material

Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.

Arya, V. K.↗