Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal response functions”

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 613 records · Page 34

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗

MEMS Reaction Control and Maneuvering for Picosat Beyond LEO

The MEMS Reaction Control and Maneuvering for Picosat Beyond LEO project will further develop a multi-functional small satellite technology for low-power attitude control, or orientation, of picosatellites beyond low Earth orbit (LEO). The Film-Evaporation MEMS Tunable Array (FEMTA) concept initially developed in 2013, is a thermal valving system which utilizes capillary forces in a microchannel to offset internal pressures in a bulk fluid. The local vapor pressure is increased by resistive film heating until it exceeds meniscus strength in a nozzle which induces vacuum boiling and provides a stagnation pressure equal to vapor pressure at that point which is used for propulsion. Interplanetary CubeSats can utilize FEMTA for high slew rate attitude corrections in addition to desaturating reaction wheels. The FEMTA in cooling mode can be used for thermal control during high-power communication events, which are likely to accompany the attitude correction. Current small satellite propulsion options are limited to orbit correction whereas picosatellites are lacking attitude control thrusters. The available attitude control systems are either quickly saturated reaction wheels or movable high drag surfaces with long response times.

Low-power attitude control↗

Thermal Model Performance for the James Webb Space Telescope OTIS Cryo-Vacuum Test

The James Webb Space Telescope (JWST), set to launch in mid-2020, is currently undergoing a series of system-level environmental tests to verify its workmanship and end-to-end functionality. As part of this series, the Optical Telescope Element and Integrated Science Instrument Module (OTIS) Cryo-Vacuum (CV) test, the most complex cryogenic test executed to date by NASA, has recently been completed at the Johnson Space Center's Chamber A facility. The OTIS CV test was intended as a comprehensive test of the integrated instrument and telescope systems to fully understand its optical, structural, and thermal performance within its intended flight environment. Due to its complexity, extensive pre-test planning was required to ensure payload safety and compliance with all limits and constraints. A system-level pre-test thermal model was constructed which fully captured the behavior of the payload, ground support equipment, and surrounding test chamber. This thermal model simulated both the transient cooldown to and warmup from a 20 K flight-like environment, as well as predicted the payload performance at cryo-stable conditions. The current work is an assessment of thermal model pre-test prediction performance against actual payload response during the OTIS CV test. Overall, the thermal model performed exceedingly well at predicting schedule and payload response. Looking in depth, this work examines both the benefits and shortcomings of assumptions made pre-test to simplify model execution when compared against test data. It explores in detail the role of temperature-dependent emissivities during transition to cryogenic temperatures, as well as the impact that model geometry simplifications have on tracking of critical hardware limits and constraints. This work concludes with a list of recommendations to improve the accuracy of thermal modeling for future large cryogenic tests. The insight gained from the OTIS CV test thermal modeling will benefit planning and execution for upcoming cryogenic missions.

Testing↗

Intermetallics Based on Sodium Chalcogenides Promote Stable Electrodeposition–Electrodissolution of Sodium Metal Anodes

Sodiophilic micro-composite films of sodium-chalcogenide intermetallics (Na 2 Te and Na 2 S) and Cu particles are fabricated onto commercial copper foam current collectors (Na 2 Te@CF and Na 2 S@CF). For the first time a controllable capacity thermal infusion process is demonstrated. Enhanced wetting by the metal electrodeposition leads to state-of-the-art electrochemical performance. For example, Na 2 Te@CF-based half-cells demonstrate stable cycling at 6 mA cm -2 and 6 mAh cm -2 , corresponding to 54 µm of Na electrodeposited/electrodissolved by geometric area. Sodium metal batteries with Na 3 V 2 (PO 4 ) 3 cathodes are stable at 30C (7 mA cm -2 ) and for 10 000 cycles at 5C and 10C. Cross-sectional cryogenic focused ion beam (cryo-FIB) microscopy details deposited and remnant dissolved microstructures. Sodium metal electrodeposition onto Na 2 Te@CF is dense, smooth, and free of dendrites or pores. On unmodified copper foam, sodium grows in a filament-like manner, not requiring cycling to achieve this geometry. Substrate–metal interaction critically affects the metal–electrolyte interface, namely the thickness and morphology of the solid electrolyte interphase. Density functional theory and mesoscale simulations provide insight into support-adatom energetics, nucleation response, and early-stage morphological evolution. On Na 2 Te sodium atomic dispersion is thermodynamically more stable than isolated clusters, leading to conformal adatom coverage of the surface.

36 MATERIALS SCIENCE↗

An Integrated Approach to Modeling Solar Electric Propulsion Vehicles During Long Duration, Near-Earth Orbit Transfers

Recent NASA interest in utilizing solar electronic propulsion (SEP) technology to transfer payloads, e.g. from low-Earth orbit (LEO) to higher energy geostationary-Earth orbit (GEO) or to Earth escape, has necessitated the development of high fidelity SEP vehicle models and simulations. These models and simulations need to be capable of capturing vehicle dynamics and sub-system interactions experienced during the transfer trajectories which are typically accomplished with continuous-burn (potentially interrupted by solar eclipse), long duration "spiral out" maneuvers taking several months or more to complete. This paper presents details of an integrated simulation approach achieved by combining a high fidelity vehicle simulation code with a detailed solar array model. The combined simulation tool gives researchers the functionality to study the integrated effects of various vehicle sub-systems (e.g. vehicle guidance, navigation and control (GN&C), electric propulsion system (EP)) with time varying power production. Results from a simulation model of a vehicle with a 50 kW class SEP system using the integrated tool are presented and compared to the results from another simulation model employing a 50 kW end-of-life (EOL) fixed power level assumption. These models simulate a vehicle under three degree of freedom dynamics (i.e. translational dynamics only) and include the effects of a targeting guidance algorithm (providing a "near optimal" transfer) during a LEO to near Earth escape (C (sub 3) = −2.0 km (sup 2) / sec (sup −2) spiral trajectory. The presented results include the impact of the fully integrated, time-varying solar array model (e.g. cumulative array degradation from traversing the Van Allen belts, impact of solar eclipses on the vehicle and the related temperature responses in the solar arrays due to operating in the Earth's thermal environment, high fidelity array power module, etc.); these are used to assess the impact on vehicle performance (i.e. propellant consumption) and transit times.

vehicle↗

Development of a reduced model for energetic particle transport by sawteeth in tokamaks

We report that the sawtooth instability is known for inducing transport and loss of energetic particles (EPs), and for generating seed magnetic islands that can trigger tearing modes. Both effects degrade the overall plasma performance. Several theories and numerical models have been previously developed to quantify the expected EP transport caused by sawteeth, with various degrees of sophistication to differentiate the response of EPs at different energies and on different orbits (e.g. passing vs. trapped), although the analysis is frequently limited to a single time slice during a tokamak discharge. This work describes the development and initial benchmark of a framework that enables a reduced model for EP transport by sawteeth retaining the full EP phase-space information. The model, implemented in the ORBIT hamiltonian particle-following code, can be used either as a standalone post-processor taking input data from codes such as TRANSP, or as a pre-processor to compute transport coefficients that can be fed back to TRANSP for time-dependent simulations including the effects of sawteeth on EPs. The advantage of the latter approach is that the evolution of the EP distribution can be simulated quantitatively for sawtoothing discharges, thus enabling a more accurate modeling of sources, sinks and overall transport properties of EP and thermal plasma species for comprehensive physics studies that require detailed information of the fast-ion distribution function and its evolution over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inflatable Softgoods Testing Capability Development for Habitation

The Habitation Systems Development Office (HP40) at NASA Marshall Space Flight Center supports systems engineering, integration, and project management for next generation space habitats. The office is responsible for formulation activities for Lunar Surface Habitat and Mars Transit Habitat. HP40 also manages commercial partnerships for exploration habitat development. The Lunar Surface Habitat is intended to initially support a minimum of two crew on the lunar surface for an approximately 30-day mission, with the ability to accommodate four crew for short durations. The Mars Transit Habitat will transport four crew to and from Mars orbit on a potentially 1,200 day mission. Both concepts are envisioned to leverage inflatable softgoods as the primary structural material for the habitation portion of the architecture. Inflatable softgoods are a complex, multi-material system consisting of multiple layers which perform specific functions. A bladder layer (polymer) contains the habitat’s atmosphere. The bladder is indexed to a woven restraint layer, usually Kevlar or Vectran, which serves as the structural layer. The remaining outer layers are for micro-meteoroid orbital debris (MMOD) shielding, protection from atomic Oxygen (if the material system is to be deployed in low earth orbit), and thermal insulation. The primary advantage of inflatable softgoods is that they can offer a larger volume per unit mass relative to a rigid metal structure. Inflatable softgoods have a strong development heritage within NASA, beginning with work in the 1960s in partnership with Goodyear Aerospace and continuing with Johnson Space Center’s TransHab project in the 1990s . Some development efforts subsequently transitioned into the private sector to mature the material technology for flight applications, including human habitation. In August 2022, NASA published its guidelines for certification of crewed inflatable softgoods structures (JSC-67721). This document outlines a series of tests which will be needed to certify inflatable softgoods material systems for use in crewed mission applications. It includes burst testing at the subscale and full-scale to determine the ultimate burst pressure of the softgoods architecture. Creep (long term deformation of the material under stress) is a key failure mode which must be characterized. Creep testing at the component and module levels is also recommended to generate time to failure curves and verify that the inflatable system has a predicted operational life which exceeds mission parameters. ET (Test Engineering) at MSFC and HP40 have worked together to establish softgoods testing capabilities for the center. This includes the ability to conduct module level burst tests and creep tests. The poster will highlight these capabilities and provide examples of testing which has been performed to date by ET. NASA MSFC anticipates that these test capabilities will continue to offer critical support for companies pursuing inflatable softgoods for exploration applications (lunar surface habitation, Mars transit) as well as commercial low earth orbit platforms.

habitation↗

Influence of Compositional Complexity on Amorphization Resistance of Swift Heavy Ion Irradiated Titanate Pyrochlores

Compositionally complex oxides have garnered attention recently for their potential technological applications in harsh environments such as thermal barrier coatings and nuclear waste forms. Therefore, their response to extreme conditions, including high temperature and intense irradiation fields, must be thoroughly investigated. Here, the structural evolution of two pyrochlore oxides with comparable cation size ratio, r A /r B , (Yb 0.2 Er 0.2 Dy 0.2 Tb 0.2 Gd 0.2 ) 2 Ti 2 O 7 and Ho 2 Ti 2 O 7 , was evaluated after irradiation with 946 MeV Au ions up to a fluence of 8 × 10 12 ions/cm2 using synchrotron X-ray diffraction, transmission electron microscopy, and Raman spectroscopy. The overall radiation response is comparable for both titanate oxides and is dominated by a loss of crystallinity. When compared to a series of conventional titanate pyrochlore compositions, the amorphous track diameter of (Yb 0.2 Er 0.2 Dy 0.2 Tb 0.2 Gd 0.2 ) 2 Ti 2 O 7 is slightly larger than that of Ho2Ti2O7 and more in line with the diameter of the endmember with the maximum A-site cation size (Gd 2 Ti 2 O 7 ). Density functional theory calculations suggest that this behavior may be linked to local lattice distortions and the associated energetics of cation antisite formation. TEM and Raman analyses show that a disordered, crystalline shell surrounds the amorphous ion tracks in (Yb 0.2 Er 0.2 Dy 0.2 Tb 0.2 Gd 0.2 ) 2 Ti 2 O 7 , and the corresponding short-range structure resembles a weberite-type atomic arrangement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-principles elucidation of defect-mediated Li transport in hexagonal boron nitride

Hexagonal boron nitride (hBN) is a promising candidate as a protective membrane or separator in Li-ion and Li–S batteries, given its excellent chemical stability, mechanical robustness, and high thermal conductivity. In addition, hBN can be functionalized by introducing defects and dopants, or be directly integrated into other active components of batteries, which further augments its appeal to the field. Here, we use first-principles simulations to evaluate the role of atomic defects in hBN in regulating the Li-ion diffusion mechanism and associated kinetics. Specifically, the following four distinct types of vacancy defects are considered: isolated single B and N vacancies, a B–N vacancy pair, and a B 3 N vacancy cluster. It is found that these defect sites generally favor Li intercalation and out-of-plane diffusion but slow down in-plane Li-ion diffusion due to a strong Li trapping effect at the defect sites. Such a trapping effect is, however, highly local such that it does not necessarily affect the overall Li-ion conductivity in defected hBN layers. The present systematic evaluation of the impact of atomic defects on Li ion migration and accompanied charge analysis of hBN lattice in response to Li-ion diffusion provide a mechanistic understanding of Li-ion transport behavior in defected hBN and highlight the potential of defect engineering to achieve optimal material performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive Forced Response Analysis of J2X Turbine Bladed-Discs with 360 Degree Variation in CFD Loading

The temporal frequency content of the dynamic pressure predicted by a 360 degree computational fluid dynamics (CFD) analysis of a turbine flow field provides indicators of forcing function excitation frequencies (e.g., multiples of blade pass frequency) for turbine components. For the Pratt and Whitney Rocketdyne J-2X engine turbopumps, Campbell diagrams generated using these forcing function frequencies and the results of NASTRAN modal analyses show a number of components with modes in the engine operating range. As a consequence, forced response and static analyses are required for the prediction of combined stress, high cycle fatigue safety factors (HCFSF). Cyclically symmetric structural models have been used to analyze turbine vane and blade rows, not only in modal analyses, but also in forced response and static analyses. Due to the tortuous flow pattern in the turbine, dynamic pressure loading is not cyclically symmetric. Furthermore, CFD analyses predict dynamic pressure waves caused by adjacent and non-adjacent blade/vane rows upstream and downstream of the row analyzed. A MATLAB script has been written to calculate displacements due to the complex cyclically asymmetric dynamic pressure components predicted by CFD analysis, for all grids in a blade/vane row, at a chosen turbopump running speed. The MATLAB displacements are then read into NASTRAN, and dynamic stresses are calculated, including an adjustment for possible mistuning. In a cyclically symmetric NASTRAN static analysis, static stresses due to centrifugal, thermal, and pressure loading at the mode running speed are calculated. MATLAB is used to generate the HCFSF at each grid in the blade/vane row. When compared to an approach assuming cyclic symmetry in the dynamic flow field, the current approach provides better assurance that the worst case safety factor has been identified. An extended example for a J-2X turbopump component is provided.

Elrod, David↗

High Accuracy In-Flight Wavelength Calibration of Imaging Spectrometry Data

Accurate wavelength calibration of imaging spectrometer data is essential if proper atmospheric transmission corrections are to be made to obtain apparent surface reflectance. Accuracies of 0.1 nm are necessary for a 10 nm-sampling instrument in order to match the slopes of the deep atmospheric water vapor features that dominate the 0.7-2.3 micrometer regions. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is calibrated in the laboratory to determine the wavelength position and full-width-half-maximum (FWHM) response for each of the 224 channels. The accuracies are limited by the quality of the monochromator used as a source. The accuracies vary from plus or minus to plus or minus 1.5 nm depending on the wavelength region, in general decreasing with increasing wavelength. Green et al. make corrections to the wavelength calibrations by using the known positions of 14 atmospheric absorption features throughout the 0.4-2.5 micrometer wavelength region. These features, having varying width and intensity, were matched to the MODTRAN model with a non-linear least squares fitting algorithm. A complete calibration was developed for all 224 channels by interpolation. Instrument calibration cannot be assumed to be stable to 0.1 nm over a flight season given the rigors of thermal cycling and launch and landing loads. The upcoming sensor HYDICE will require a means for in-flight spectral calibration of each scene because the calibration is both temperature and pressure sensitive. In addition, any sensor using a two-dimensional array has the potential for systematic wavelength shifts as a function of cross-track position, commonly called 'smile'. Therefore, a rapid means for calibrating complete images will be required. The following describes a method for determining instrument wavelength calibration using atmospheric absorption features that is efficient enough to be used for entire images on workstations. This study shows the effect of the surface reflectance on the calibration accuracy and the calibration history for the AVIRIS B spectrometer over the 1992 flight season.

Goetz, Alexander F. H.↗

Characteristic Ultraviolet/Optical Timescales in Active Galactic Nuclei

The UV/optical light curves of active galactic nuclei (AGNs) may be roughly characterized by a periodic fractional flux variations of approx. 10% on timescales of 1 month. The physical mechanism(s) responsible remain ill defined. We present a structure function analysis, i.e., measure the power distribution over a range of timescales tau, of 13 AGNs to constrain the origin of UV/optical emission. On timescales tau approx. 5- 60 days, the mean UV and optical power density spectra (PDS) are equivalent. This may suggest that the underlying energy generating mechanism is identical. The combined W/optical PDS is P(f) proportional to f(sup -alpha)- with alpha = 2.13(sup +0.22, sub -0.06). For sources with measured X-ray PDS indices, we find they are indistinguishable from their UV/optical counterparts. This supports scenarios whereby X-rays are generated via Compton upscattering of UV photons, to later radiatively drive optical variations. At the same time, we present evidence for characteristic variability timescales tau(sub char) of approx. 5-100 days in 10 sources. These variability timescales combined with reverberation based masses M suggest a M - tau(sub char) relationship; higher mass systems have larger characteristic timescales. The UV tau(sub char) may possibly reflect dynamical or accretion disk thermal timescales. We find suggestive evidence for a dichotomy, at tau approx. 30 days and M approx. l0(exp 7) solar mass, between short- and long-time scale optical variations. These optical variations may be attributable to dynamical and accretion disk thermal or starburst activity timescales, respectively.

Collier, Stefan↗

Pressure-Fed LOX/LCH4 Reaction Control System for Spacecraft: Transient Modeling and Thermal Vacuum Hotfire Test Results

An integrated cryogenic liquid oxygen, liquid methane (LOX/LCH4) reaction control system (RCS) was tested at NASA Glenn Research Center's Plum Brook Station in the Spacecraft Propulsion Research Facility (B-2) under vacuum and thermal vacuum conditions. The RCS is a subsystem of the Integrated Cryogenic Propulsion Test Article (ICPTA), a pressure-fed LOX/LCH4 propulsion system composed of a single 2,800 lbf main engine, two 28 lbf RCS engines, and two 7 lbf RCS engines. Propellants are stored in four 48 inch diameter 5083 aluminum tanks that feed both the main engine and RCS engines in parallel. Helium stored cryogenically in a composite overwrapped pressure vessel (COPV) flows through a heat exchanger on the main engine before being used to pressurize the propellant tanks to a design operating pressure of 325 psi. The ICPTA is capable of simultaneous main engine and RCS operation. The RCS engines utilize a coil-on-plug (COP) ignition system designed for operation in a vacuum environment, eliminating corona discharge issues associated with a high voltage lead. There are two RCS pods on the ICPTA, with two engines on each pod. One of these two engines is a heritage flight engine from Project Morpheus. Its sea level nozzle was removed and replaced by an 85:1 nozzle machined using Inconel 718, resulting in a maximum thrust of 28 lbf under altitude conditions. The other engine is a scaled down version of the 28 lbf engine, designed to match the core and overall mixture ratios as well as other injector characteristics. This engine can produce a maximum thrust of 7 lbf with an 85:1 nozzle that was additively manufactured using Inconel 718. Both engines are film-cooled and capable of limited duration gas-gas and gas-liquid operation, as well as steady-state liquid-liquid operation. Each pod contains one of each version, such that two engines of the same thrust level can be fired as a couple on opposite pods. The RCS feed system is composed of symmetrical 3/8 inch lines that tap off of the main propellant manifold to send LOX and LCH4 outboard to the RCS pods. A Thermodynamic Vent System (TVS) is used to condition propellants at each pod by venting through an orifice and then routing the cold expansion products back through tubing that is welded along a large portion of the main RCS feed lines. Prior to final installation on the ICPTA, the RCS engines were tested in a small vacuum chamber at the Johnson Space Center (JSC) Energy Systems Test Area (ESTA) to verify functionality of the new COP ignition system and check out operation of the vacuum nozzles. After engine-level testing, the RCS engines were installed on the vehicle and a series of integrated hot-fire tests were performed at JSC consisting of various pulsing and steady-state firings as well as integrated main engine/RCS operation. The ICPTA was then integrated into the Plum Brook B-2 facility for vacuum and thermal/vacuum testing. Testing in the B-2 facility was composed of multiple thermal and pressure environments. The first set of tests were performed under ambient temperature and altitude pressure conditions. These tests consisted of a range of minimum impulse bit (MIB) pulsing sequences with low duty cycle, analogous to a coast phase in which the RCS is primarily used for station keeping. The primary goal of this sequence is to understand how propellant conditions were effected without an active TVS. In this scenario, consistent gas-gas operation is desirable since it results in a smaller MIB and more efficient propellant consumption. Multiple skin thermocouples are mounted on the feedlines, in addition to a submerged thermocouple on each commodity, in order to gather thermal data on the system. Higher duty cycle pulsing tests were then performed, analogous to an ascent or landing mission phase. The primary goal of this sequence was to examine how well the engines self-conditioned without active TVS when starting from a quiescent state. The TVS was then activated during some tests to demonstrate the capability to quickly condition the engines for higher pulsing demand scenarios. A thermocouple at the TVS outlet allows for the calculation of energy absorbed by the vented propellant. Lastly, tests with longer pulses and multiple engines firing either in sequence or simultaneously were run in order to gather transient system response data on waterhammer. Six total high-speed pressure transducers are installed on the RCS system, one sensor at the end of each propellant manifold line on the pods, and one at the tap-off location for each commodity. This will allow for the accurate characterization of waterhammer in the system under various propellant conditions and firing sequences. Other instrumentation for this test series includes nozzle throat thermocouples, chamber pressure measurement, heat soakback measurement, and tank wall plume impingement temperature measurement. The next set of tests were performed to demonstrate simultaneous main engine and RCS operation. Data from this test will be used to examine if there is any change to nominal operation of the RCS as a result of feed system interaction or other phenomenon. Some of these tests began under high vacuum conditions (target ambient pressure less than 1x10(exp -3) torr) and others began at altitude conditions. The last set of tests were performed with the B-2 cold wall active. Under these tests, many of the same low duty cycle MIB tests were repeated in order to characterize how propellant conditions changed with the lower heat leak. In this scenario the RCS manifold experiences much less heat leak, resulting in a change to how well the engines self-condition. As a result, an increase in maximum waterhammer pressures and a change in natural frequency of the system was expected due to higher density propellants. The lower heat leak should also result in a change to the MIB pulse profile, and data will be examined to understand how MIB repeatability is affected in the different operating environments. Parallel to the test efforts, a set of transient model development efforts were made to predict RCS performance. The primary effort was aimed at producing a SINDA/FLUINT model to predict propellant conditioning up to the engine inlet as a function of different environmental and operating parameters, with the goal of predicting chamber pressure, TVS performance, and propellant consumption over time. Preliminary results for this effort will be presented in comparison with test data. Additional modeling efforts were made using SINDA/FLUINT to predict waterhammer in the system since the software is capable of handling multiphase transient fluid dynamics. These results will be compared with the high-speed pressure transducer test data for validation purposes.

Atwell, Matthew J.↗

Effect of Diphenyl Content on Viscoelasticity of Poly(dimethyl- co -diphenyl)siloxane Melt and Network

Polydimethylsiloxane (PDMS) is one of the most widely used polymeric materials for sealants, adhesives, lubricants, and thermal as well as electrical insulation. At low temperatures, however, PDMS is subject to crystallization that can cause deterioration in mechanical function. A common way to suppress such crystallization is through the incorporation of phenylsiloxane into the backbone of polysiloxane. Nevertheless, the introduction of phenyl components, even in small quantities, could potentially change the properties of the siloxane in a significant way. In this work, a series of mechanical tests and finite element simulations were performed to study the macroscale viscoelasticity of two poly(dimethyl-co-diphenyl)siloxane formulations in order to understand the effects of a few percent diphenyl contents on the viscoelasticity of the polysiloxane material. Here, we utilized the small-angle X-ray scattering to investigate the microscopic structures of the copolymers and broadband dielectric spectroscopy and rheology to probe the chain dynamics at the microscale. The results of these characterizations were used to inform the finite element simulations. We found that the degree of cross-linking does not significantly alter the microstructure but can profoundly affect the viscoelastic response of the copolymer networks. The corresponding hysteretic behavior is interpreted in terms of reptation-like motion and relaxation of the effective free chains in the cured polymer network. The relaxation of the copolymer chains is slowed significantly by even a small increase in the molar ratio of the diphenyl component.

36 MATERIALS SCIENCE↗

Obliquity variation in a Mars climate evolution model

The existence of layered terrain in both polar regions of Mars is strong evidence supporting a cyclic variation in climate. It has been suggested that periods of net deposition have alternated with periods of net erosion in creating the layered structure that is seen today. The cause for this cyclic climatic behavior is variation in the annually averaged latitudinal distribution of solar insolation in response to obliquity cycles. For Mars, obliquity variation leads to major climatological excursion due to the condensation and sublimation of the major atmospheric constituent, CO2. The atmosphere will collapse into the polar caps, or existing caps will rapidly sublimate into the atmosphere, dependent upon the polar surface heat balance and the direction of the change in obliquity. It has been argued that variations in the obliquity of Mars cause substantial departures from the current climatological values of the surface pressure and the amount of CO2 stored in both the planetary regolith and polar caps. In this new work we have modified the Haberle et al. model to incorporate variable obliquity by allowing the polar and equatorial insolation to become functions of obliquity, which we assume to vary sinusoidally in time. As obliquity varies in the model, there can be discontinuities in the time evolution of the model equilibrium values for surface pressure, regolith, and polar cap storage. The time constant, tau r, for the regolith to find equilibrium with the climate is estimated--depending on the depth, thermal conductivity, and porosity of the regolith--between 10(exp 4) and 10(exp 6) yr. Thus, using 2000-yr timesteps to move smoothly through the 0.1250 m.y. obliquity cycles, we have an atmosphere/regolith system that cannot be assumed in equilibrium. We have dealt with this problem by limiting the rate at which CO2, can move between the atmosphere and regolith, mimicking the diffusive nature and effects of the temperature and pressure waves, by setting the time rate of change of regolith storage proportional to the difference between equilibrium storage and current storage.

Tyler, D.↗

Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sierra/SD – Verification Test Manual – 5.22

Verification and validation (V&V) of scientific computing programs are important at Sandia National Labs due to the expanding role of computational simulation in managing the United States nuclear stockpile. The complexities of structural response calculations used to analyze physical problems, the varieties of codes applied to the calculations, and the importance of accurate predictions when assessing field conditions demand confidence in the consistency and accuracy of computer codes. Confidence in the accuracy of the predictions arising from computer simulations must ultimately be gained through verification and validation. The Sierra salinas structural dynamics analysis code, Sierra/SD, is used at the DOE Laboratories, and in several DOD projects. The roles of Sierra/SD in the qualification of weapon systems and components for normal and hostile environments throughout the Stockpile-to-Target Sequence include to, • Redesign weapon components. • Certify weapon components and systems for target environments such as hypersonic vehicles. • Certify that components will survive the thermal mechanical shock loads associated with hostile environments. • Evaluate current stockpile issues, including issues associated with uncertainty quantification. • Address many other problems that are encountered in stockpile management. The Sierra/SD verification plan is described, and an evolving set of key verification tests are described in detail. The verification tests ensure the correctness of the mathematics and numerical algorithms associated with functionality describing engineering phenomena. Development is in accordance with a set of tailored Software Quality Engineering (SQE) practices. SQE practices guide the overall verification and validation effort.

97 MATHEMATICS AND COMPUTING↗

Performance enhancement of aqueous ionic liquids with lower critical solution temperature (LCST) behavior through ternary mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid–liquid phase separation when mixed with water and heated above a lower critical solution temperature (LCST), resulting in a water-rich (WR) and an IL-rich (ILR) phase. These binary IL–water mixtures can be employed in a variety of thermodynamic processes such as forward osmosis (FO) desalination, for which two solution properties are desirable: low phase separation temperature and high osmotic strength (osmolality). However, these two properties are interlinked, with ILs that exhibit higher osmotic strengths typically requiring higher phase separation temperatures. This behavior tends to arise from the hydrophilicity of the IL cations, which enhances osmotic strength while also elevating the phase separation temperature. In this work, we highlight a pathway to overcome this tradeoff by developing ternary IL mixtures (two ILs with varying cation hydrophilicity mixed with water), which lowers the phase separation temperature while maintaining and even enhancing the osmotic strength of the solution. We characterize the mixing behavior (osmolality, phase separation temperature, WR phase purity, and WR to ILR phase mass ratio) of four ILs as a function of their concentration in solution. We find that an enhancement of up to 81.6% in the osmolality with a concomitant reduction of up to 15.4% in the phase separation temperature can be achieved using this approach. The ternary mixture is also shown to improve the phase separation kinetics by nearly 95% compared to the binary mixture. Overall, this work highlights a new pathway to improve the performance of LCST ILs for water and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗