Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric 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 811 records · Page 45

Attributes of GRB Pulses: Analysis of BATSE TTE Data

Our extensive program of modeling GRB profiles is aimed at elucidating the physical processes responsible for the burst itself, as well as possible extrinsic phenomena (e.g. time dilation) as discussed in another paper in this Symposium (Norris et al., "GRB PROFILES AS COSMIC PROBES"). We have developed special methods to extract the wealth of short time-scale information contained in the BATSE time-tag event (TTE) data. Our algorithm yields a piecewise-constant representation of the light curve -- using only the raw photon arrival times, and based on Bayesian change-point methods. This representation in effect lets the data determine the bin size and location, and avoids unwanted effects due to arbitrary choices of the bin parameters. We have determined widths, separations, and amplitudes of pulses contained in the bursts, without invoking a specific pulse model. The effect of cosmic time dilation can be easily seen in a direct plot of amplitude vs. time scale for individual pulses, without the need to lump the data into a small number of brightness classes. We are also performing noise equalization on these data (to reduce a well-known bias of pulse width as a function of signal-to-noise ratio), as well as fits of parametric pulse-shape models -- including explicit energy dependence of the pulse parameters. Such refinements are expected to improve the quality and physical significance of these results.

Norris, J. P.↗

BISON capability and validation for U3Si2, Cr2O3-doped UO2, FeCrAl, and Cr-coated Zircaloy ATF concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the BISON fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on, the focus was on U3Si2 fuel and FeCrAl cladding as part of a high-impact problem through the Nuclear Energy and Advanced Model- ing Simulation program [1]. Then, developments for Cr2O3- doped UO2 fuel, as well as SiC-SiC and Cr-coated Zircaloy cladding began based upon industry interests. Near the end of fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors took over further Accident Tolerant Fuel (ATF) work in BISON in support of Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC clarified their list of priority fuel and cladding concepts, which included Cr2O3-doped UO2 and U3Si2 fuels, as well as Cr-coated Zircaloy and FeCrAl claddings. Further, reports similar in form to NUREG/CR-7024 [2] that was developed for traditional light water reactor (LWR) materials, UO2 fuel with Zircaloy claddings were created for the priority fuel and cladding concepts. In this paper select details from the reports entitled "ATF material model development and validation for priority fuel concepts" [3] and "ATF material model development and validation for priority cladding concepts" [4] are presented. First, a brief overview of the available models and their associated uncertainty are provided for the priority fuel and cladding concepts. Then, for the concepts with experimental data, comparisons between simulation and experiment are shown taking into account forward uncertainty quantification (UQ) techniques. For the concepts without publicly available experimental data, parametric studies incorporating forward UQ are completed to verify claims made about the materials’ behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parametric Analysis of Life Support Systems for Future Space Exploration Missions

The National Aeronautics and Space Administration is in a process of evaluating future targets for space exploration. In order to maintain the welfare of a crew during future missions, a suite of life support technology is responsible for oxygen and water generation, carbon dioxide control, the removal of trace concentrations of organic contaminants, processing and recovery of water, and the storage and reclamation of solid waste. For each particular life support subsystem, a variety competing technologies either exist or are under aggressive development efforts. Each individual technology has strengths and weaknesses with regard to launch mass, power and cooling requirements, volume of hardware and consumables, and crew time requirements for operation. However, from a system level perspective, the favorability of each life support architecture is better assessed when the sub-system technologies are analyzed in aggregate. In order to evaluate each specific life support system architecture, the measure of equivalent system mass (ESM) was employed to benchmark system favorability. Moreover, the results discussed herein will be from the context of loop-closure with respect to the air, water, and waste sub-systems. Specifically, closure relates to the amount of consumables mass that crosses the boundary of the vehicle over the lifetime of a mission. As will be demonstrated in this manuscript, the optimal level of loop closure is heavily dependent upon mission requirements such as duration and the level of extra-vehicular activity (EVA) performed. Sub-system level trades were also considered as a function of mission duration to assess when increased loop closure is practical. Although many additional factors will likely merit consideration in designing life support systems for future missions, the ESM results described herein provide a context for future architecture design decisions toward a flexible path program.

Swickrath, Michael J.↗

Rapid Airplane Parametric Input Design (RAPID)

RAPID is a methodology and software system to define a class of airplane configurations and directly evaluate surface grids, volume grids, and grid sensitivity on and about the configurations. A distinguishing characteristic which separates RAPID from other airplane surface modellers is that the output grids and grid sensitivity are directly applicable in CFD analysis. A small set of design parameters and grid control parameters govern the process which is incorporated into interactive software for 'real time' visual analysis and into batch software for the application of optimization technology. The computed surface grids and volume grids are suitable for a wide range of Computational Fluid Dynamics (CFD) simulation. The general airplane configuration has wing, fuselage, horizontal tail, and vertical tail components. The double-delta wing and tail components are manifested by solving a fourth order partial differential equation (PDE) subject to Dirichlet and Neumann boundary conditions. The design parameters are incorporated into the boundary conditions and therefore govern the shapes of the surfaces. The PDE solution yields a smooth transition between boundaries. Surface grids suitable for CFD calculation are created by establishing an H-type topology about the configuration and incorporating grid spacing functions in the PDE equation for the lifting components and the fuselage definition equations. User specified grid parameters govern the location and degree of grid concentration. A two-block volume grid about a configuration is calculated using the Control Point Form (CPF) technique. The interactive software, which runs on Silicon Graphics IRIS workstations, allows design parameters to be continuously varied and the resulting surface grid to be observed in real time. The batch software computes both the surface and volume grids and also computes the sensitivity of the output grid with respect to the input design parameters by applying the precompiler tool ADIFOR to the grid generation program. The output of ADIFOR is a new source code containing the old code plus expressions for derivatives of specified dependent variables (grid coordinates) with respect to specified independent variables (design parameters). The RAPID methodology and software provide a means of rapidly defining numerical prototypes, grids, and grid sensitivity of a class of airplane configurations. This technology and software is highly useful for CFD research for preliminary design and optimization processes.

Smith, Robert E.↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Parametric Analysis of Life Support Systems for Future Space Exploration Missions

Having adopted a flexible path approach to space exploration, the National Aeronautics and Space Administration is in a process of evaluating future targets for space exploration. In order to maintain the welfare of a crew during future missions, a suite of life support technology is responsible for oxygen and water generation, carbon dioxide control, the removal of trace concentrations of organic contaminants, processing and recovery of water, and the storage and reclamation of solid waste. For each particular life support subsystem, a variety competing technologies either exist or are under aggressive development efforts. Each individual technology has strengths and weaknesses with regard to launch mass, power and cooling requirements, volume of hardware and consumables, and crew time requirements for operation. However, from a system level perspective, the favorability of each life support architecture is better assessed when the sub-system technologies are analyzed in aggregate. In order to evaluate each specific life support system architecture, the measure of equivalent system mass (ESM) was employed to benchmark system favorability. Moreover, the results discussed herein will be from the context of loop-closure with respect to the air, water, and waste sub-systems. Specifically, closure relates to the amount of consumables mass that crosses the boundary of the vehicle over the lifetime of a mission. As will be demonstrated in this manuscript, the optimal level of loop closure is heavily dependent upon mission requirements such as duration and the level of extra- vehicular activity (EVA) performed. Sub-system level trades were also considered as a function of mission duration to assess when increased loop closure is practical. Although many additional factors will likely merit consideration in designing life support systems for future missions, the ESM results described herein provide a context for future architecture design decisions toward a flexible path program.

Swickrath, Michael J.↗

Design of Mini-Plate-1 Irradiation Test for Qualification of High-Density, Low-Enriched U-10Mo Monolithic Fuel

The United States High Performance Research Reactor project is tasked with fuel development and qualification leading to conversion of higher power research and test reactors in the US from high-enriched uranium (HEU) to low-enriched uranium (LEU) fuels. Here, this manuscript identifies the functional and operational design requirements of the first miniature test plate (mini-plate [MP]) irradiation campaign (MP-1) of commercially fabricated LEU U-10Mo monolithic plate-type fuel and is the precursor to a large parametric mini-plate test (MP-2) aimed at producing the data to support regulatory qualification of the LEU U-10Mo monolithic fuel. The manuscript (1) provides a general description of the selected U-10Mo LEU fuel and (2) defines the overall experiment design and functional requirements to accomplish the specific test objective of MP-1, which is to confirm that the commercially manufactured LEU U-10Mo monolithic fuel meets the established requirements of geometric stability, mechanical integrity and stable and predictable behavior. The fuel testing parameters are established by the need to bound performance behavior within the operational envelope of the reactors being converted.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Automated detector simulation and reconstruction parametrization using machine learning

Rapidly applying the effects of detector response to physics objects (e.g. electrons, muons, showers of particles) is essential in high energy physics. Presently available tools for the transformation from truth-level physics objects to reconstructed detector-level physics objects involve manually defining resolution functions. These resolution functions are typically derived in bins of variables that are correlated with the resolution (e.g. pseudorapidity and transverse momentum). This process is time consuming, requires manual updates when detector conditions change, and can miss important correlations. Machine learning offers a way to automate the process of building these truth-to-reconstructed object transformations and can capture complex correlation for any given set of input variables. Such machine learning algorithms, with sufficient optimization, could have a wide range of applications: improving phenomenological studies by using a better detector representation, allowing for more efficient production of Geant4 simulation by only simulating events within an interesting part of phase space, and studies on future experimental sensitivity to new physics.

47 OTHER INSTRUMENTATION↗

Differential absorption lidar (DIAL) for water vapor and aerosol profiling from airborne and space-based platforms

Many current challenges in weather and climate prediction involve cloud and aerosol processes, which are intimately linked to processes that initiate or take place within the planetary boundary layer (PBL). The current space-based observing system provides a wealth of information on the spatial distribution of key atmospheric variables such as profiles of water vapor (WV) and temperature, but lacks the sensitivity, accuracy, and vertical resolution in the PBL to overcome the observational gaps that inhibit the parametrization and advancement of current climate and weather models. To overcome some of these observational gaps and improve understanding and monitoring of these key processes, mature remote sensing techniques like Differential Absorption Lidar (DIAL) can be employed from various platforms to capture processes on a wide range of spatial scales. This presentation will provide an overview of NASA’s water vapor DIAL capabilities, present the strengths, weaknesses, and synergies of the DIAL technique relative to traditional methods, and provide representative data examples from the HALO airborne water vapor DIAL which was designed as a technology testbed for future space-based DIAL missions as well as a workhorse instrument in support of suborbital process studies. HALO is an airborne multifunction lidar system, simultaneously profiling WV with the DIAL technique and aerosol and cloud optical properties and PBL heights via the high spectral resolution lidar (HSRL) technique. HALO has a compact design that allows it to fly on most research aircraft, enabling nadir profiling from the lower stratosphere down through the PBL, which is critical for studying processes such as convection and cloud formation that can link the PBL and free troposphere. We will end with an outlook on the benefits and challenges of transitioning the mature DIAL technique to space to compliment other active and passive sounders and enable new PBL observations for the future Earth Observing System. The Atmospheric Boundary layer Lidar pathfindEr (ABLE), a program aimed at advancing technologies to enable the first demonstration of DIAL from space with cross-cutting observables of PBL heights, aerosol and cloud profiling, and methane columns from space will also be presented.

Brian J. Carroll↗

Comparative Analysis of Confidence Metrics for Nuclear Criticality Safety

Nuclear criticality safety standards provide guidance on the requirements and recommendations to establish confidence in computerized model results used to support operation with fissionable materials. By design, the guidance is not prescriptive, leaving the analysts free to determine how various sources of uncertainties are to be statistically aggregated. This report compares the analyses and key assumptions behind four notable methodologies documented in the nuclear criticality safety literature: the parametric, nonparametric, Whisper, and TSURFER methodologies. Because of the involved use of statistics entangled with heuristic recipes, the results of these methodologies are often difficult to interpret. Also, they are augmented by additional large administrative margins, eliminating the incentive to understand their differences. With the new resurgent wave of advanced nuclear systems focused on economizing operation—including advanced reactors, fuel cycles, and fuel concepts—there is a strong need to develop a clear understanding of uncertainties and their fusion methodologies to reduce uncertainties in a scientifically defensible manner. This report offers a deep dive into the various assumptions of the four noted methodologies, their adequacy, and their limitations, to provide guidance on developing confidence for the emergent nuclear systems. These systems are expected to be challenged by the scarcity of experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Imaging spectrophotometry of the nuclear outflow of NGC 1068

This observational program (in conjunction with R. B. Tully (IfA, Honolulu), and J. Bland (Rice U., Houston)) aims to constrain the kinematic organization and dominant excitation mechanisms of ionized gas in active galaxies. More generally, researchers are interested in the dynamics of radiative, supersonic flows in the Interstellar Medium (ISM). Imaging Fabry-Perot interferometers and low-noise Charge Coupled Devices (CCDs) are used for complete spatial coverage of the complex gas distribution in circumnuclear narrow-line regions (NLRs). Extranuclear emission line widths in NLRs can exceed 3000 km s(-1), so to avoid inter-order confusion researchers use an etalon of 4000 km s(-1) free spectral range to map (N II) lambda lambda 6548, 6583 and H alpha. To maximize spatial resolution, researchers select nearby active systems independent of luminosity but known to possess interesting morphologies and/or high-velocity extranuclear ionized gas. Monochromatic images Full Width Half Maximum (FWHM) approx. 65 km s(-1) have thus far been obtained in 1 second or better seeing at the U. Hawaii 2.2m, CFH 3.6m, and CTIO 4.0m telescopes. These are stacked into grids of line profiles, of spectrophotometric quality, at sub-arcsecond increments across a 3 second field. To handle the approx. 20,000 to 300,000 useful spectra that arise from a typical night's work, researchers have developed a complete analysis and reduction package for VAX and Sun image workstations. Reduction involves parametrization of approx. 10 to the 8th raw data points to a few maps (e.g., velocities of each kinematic subsystem, continuum-free line fluxes) containing approx. 10 to the 5th pixels. Researchers identify kinematic and structural symmetries by examining these maps and the point to point variations of the synthesized line profiles. The combination of monochromatic images and full spatial sampling of line profiles has allowed them to isolate such symmetries and has led to reliable kinematic deprojections.

Cecil, Gerald↗

Lead-Lag Control for Helicopter Vibration and Noise Reduction

As a helicopter transitions from hover to forward flight, the main rotor blades experience an asymmetry in flow field around the azimuth, with the blade section tangential velocities increasing on the advancing side and decreasing on the retreating side. To compensate for the reduced dynamic pressure on the retreating side, the blade pitch angles over this part of the rotor disk are increased. Eventually, a high enough forward speed is attained to produce compressibility effects on the advancing side of the rotor disk and stall on the retreating side. The onset of these two phenomena drastically increases the rotor vibratory loads and power requirements, thereby effectively establishing a limit on the maximum achievable forward speed. The alleviation of compressibility and stall (and the associated decrease in vibratory loads and power) would potentially result in an increased maximum forward speed. In the past, several methods have been examined and implemented to reduce the vibratory hub loads. Some of these methods are aimed specifically at alleviating vibration at very high flight speeds and increasing the maximum flight speed, while others focus on vibration reduction within the conventional flight envelope. Among the later are several types passive as well as active schemes. Passive schemes include a variety of vibration absorbers such as mechanical springs, pendulums, and bifilar absorbers. These mechanism are easy to design and maintain, but incur significant weight and drag penalties. Among the popular active control schemes in consideration are Higher Harmonic Control (HHC) and Individual Blade Control (IBC). HHC uses a conventional swash plate to generate a multi-cyclic pitch input to the blade. This requires actuators capable of sufficiently high power and bandwidth, increasing the cost and weight of the aircraft. IBC places actuators in the rotating reference frame, requiring the use of slip rings capable of transferring enough power to the actuators. Both schemes cause an increase in pitch link loads. Trailing Edge Flap (TEF) deployment can also used to generate unsteady aerodynamic forces and moments that counter the original vibratory loads, and thereby reduce rotor vibrations. While the vibrations absorbers, HHC, IBC, and TEF concepts discussed above attempt to reduce the vibratory loads, they do not specifically address the phenomena causing the vibrations at high advance ratios. One passive method that attempts to directly alleviate compressibility and stall, instead of reducing the ensuing vibrations, is the use of advanced tip designs. Taper, sweep, anhedral, and the manipulation of other geometric properties of the blade tips can reduce the severity of stall and compressibility effects , as well as reduce rotor power. A completely different approach to solve these problems is the tiltrotor configuration. As the forward velocity of the aircraft increases, the rotors, in this case, are tilted forward until they are perpendicular to the flow and act as propellers. This eliminates the edgewise flow encountered by conventional rotors and circumvents all the problems associated with flow asymmetry. However, the success involves a tremendous increase in cost and complexity of the aircraft. Another possible approach that has been proposed for the alleviation of vibratory loads at high forward flight speeds involves the use of controlled lead-lag motions to reduce the asymmetry in flow. A correctly phased 1/rev controlled lag motion could be introduced such that it produces a backward velocity on the advancing side and a forward velocity on the retreating side, to delay compressibility effects and stall to a higher advance ratio. Using a large enough lead-lag amplitude, the tip velocities could be reduced to levels encountered in hover. This concept was examined by two groups in the 1950's and early 1960's. In the United States, the Research Labs Division of United Aircraft developed a large lead-lag motion rotor, meant to achieve lag motion amplitudes up to 45 degrees. In order to reduce the required actuation force, the blade hinges were moved to 40% of the blade radius to increase the rotating lag frequency to approximately 1/rev. The blade hinges were redesigned to produce a flap-lag coupling so the large flapwise aerodynamic loads could be exploited to actuate the blades in the lag direction. A wind tunnel test of this rotor concept revealed actuation and blade motion scheduling problems. The project was eventually discontinued due to these problems and high blade stresses. Around the same time, at Boelkow in Germany, a similar lead-lag rotor program was conducted under the leadership of Hans Derschmidt. Here, too, the blade hinges were moved outboard to 34% radius to reduce the actuation loads. The main difference between this and the United Aircraft program was the use of a mechanical actuation scheme with maximum lead-lag motions of 400. This program was also discontinued for unclear reasons. The present study is directed toward conducting a comprehensive analytical examination to evaluate the effectiveness of controlled lead-lag motions in reducing vibratory hub loads and increasing maximum flight speed. Since both previous studies on this subject were purely experimental, only a limited data set and physical understanding of the problem was obtained. With the currently available analytical models and computational resources, the present effort is geared toward developing an in-depth physical understanding of the precise underlying mechanisms by which vibration reduction may be achieved. Additionally, in recognition of the fact that large amplitude lead-lag motions would - (i) be difficult to implement, and (ii) produce very large blade stresses; the present study examines the potential of only moderate-to-small lead-lag motions for reduction of vibratory hub loads. Using such an approach, the emphasis is not on eliminating the periodic variations in tangential velocity at the blade tip, but at best reducing these variations slightly so that compressibility and stall are delayed to slightly higher advance ratios. This study was conducted in two steps. In the first step, a hingeless helicopter rotor was modeled using rigid blades undergoing flap-lag-torsion rotations about spring restrained hinges and bearings. This model was then modified by separating the lead-lag degree of freedom into two components, a free and a prescribed motion. Using this model, a parametric study of the effect of phase and amplitude of a prescribed lead-lag motion on hub vibration was conducted. The data gathered was analyzed to obtain an understanding of the basic physics of the problem and show the capability of this method to reduce vibration and expand the flight envelope. In the second half of the study, the similar analysis was conducted using an elastic blade model to confirm the effects predicted by the simpler model.

Gandhi, Farhan↗

Cryogenic Propellant Feed System Analytical Tool Development

The Propulsion Systems Branch at NASA s Lyndon B. Johnson Space Center (JSC) has developed a parametric analytical tool to address the need to rapidly predict heat leak into propellant distribution lines based on insulation type, installation technique, line supports, penetrations, and instrumentation. The Propellant Feed System Analytical Tool (PFSAT) will also determine the optimum orifice diameter for an optional thermodynamic vent system (TVS) to counteract heat leak into the feed line and ensure temperature constraints at the end of the feed line are met. PFSAT was developed primarily using Fortran 90 code because of its number crunching power and the capability to directly access real fluid property subroutines in the Reference Fluid Thermodynamic and Transport Properties (REFPROP) Database developed by NIST. A Microsoft Excel front end user interface was implemented to provide convenient portability of PFSAT among a wide variety of potential users and its ability to utilize a user-friendly graphical user interface (GUI) developed in Visual Basic for Applications (VBA). The focus of PFSAT is on-orbit reaction control systems and orbital maneuvering systems, but it may be used to predict heat leak into ground-based transfer lines as well. PFSAT is expected to be used for rapid initial design of cryogenic propellant distribution lines and thermodynamic vent systems. Once validated, PFSAT will support concept trades for a variety of cryogenic fluid transfer systems on spacecraft, including planetary landers, transfer vehicles, and propellant depots, as well as surface-based transfer systems. The details of the development of PFSAT, its user interface, and the program structure will be presented.

Lusby, Brian S.↗

Transient responses of phosphoric acid fuel cell power plant system

An analytical and computerized study of the steady state and transient response of a phosphoric acid fuel cell (PAFC) system was completed. Parametric studies and sensitivity analyses of the PAFC system's operation were accomplished. Four non-linear dynamic models of the fuel cell stack, reformer, shift converters, and heat exchangers were developed based on nonhomogeneous non-linear partial differential equations, which include the material, component, energy balance, and electrochemical kinetic features. Due to a lack of experimental data for the dynamic response of the components only the steady state results were compared with data from other sources, indicating reasonably good agreement. A steady state simulation of the entire system was developed using, nonlinear ordinary differential equations. The finite difference method and trial-and-error procedures were used to obtain a solution. Using the model, a PAFC system, that was developed under NASA Grant, NCC3-17, was improved through the optimization of the heat exchanger network. Three types of cooling configurations for cell plates were evaluated to obtain the best current density and temperature distributions. The steady state solutions were used as the initial conditions in the dynamic model. The transient response of a simplified PAFC system, which included all of the major components, subjected to a load change was obtained. Due to the length of the computation time for the transient response calculations, analysis on a real-time computer was not possible. A simulation of the real-time calculations was developed on a batch type computer. The transient response characteristics are needed for the optimization of the design and control of the whole PAFC system. All of the models, procedures and simulations were programmed in Fortran and run on IBM 370 computers at Cleveland State University and the NASA Lewis Research Center.

Lu, Cheng-Yi↗

Design Study of Propulsion and Drive Systems for the Large Civil TiltRotor (LCTR2) Rotorcraft

Boeing, Rolls Royce, and NASA have worked together to complete a parametric sizing study for NASA's Large Civil Tilt Rotor (LCTR2) concept 2nd iteration. Vehicle gross weight and fuel usage were evaluated as propulsion and drive system characteristics were varied to maximize the benefit of reduced rotor tip speed during cruise conditions. The study examined different combinations of engine and gearbox variability to achieve rotor cruise tip speed reductions down to 54% of the hover tip speed. Previous NASA studies identified that a 54% rotor speed reduction in cruise minimizes vehicle gross weight and fuel burn. The LCTR2 was the study baseline for initial sizing. This study included rotor tip speed ratios (cruise to hover) of 100%, 77% and 54% at different combinations of engine RPM and gearbox speed reductions, which were analyzed to achieve the lightest overall vehicle gross weight (GW) at the chosen rotor tip speed ratio. Different engine and gearbox technology levels are applied ranging from commercial off-the-shelf (COTS) engines and gearbox technology to entry-in-service (EIS) dates of 2025 and 2035 to assess the benefits of advanced technology on vehicle gross weight and fuel burn. Interim results were previously reported1. This technical paper extends that work and summarizes the final study results including additional engine and drive system study accomplishments. New vehicle sizing data is presented for engine performance at a single operating speed with a multispeed drive system. Modeling details for LCTR2 vehicle sizing and subject engine and drive sub-systems are presented as well. This study was conducted in support of NASA's Fundamental Aeronautics Program, Subsonic Rotary Wing Project.

Robuck, Mark↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Impact of Spontaneous Fission Neutron Emission Energy on Neutron Detector Response

A parametric study was simulated in MCNP to investigate the impact of spontaneous fission neutron energy on neutron detector response. A total of 880 configurations were simulated by varying 3 He pressure and polyethylene reflector thickness for 8 different isotopes. Results showed that deviations between the simulated isotopes become more prevalent at higher energies with 238 U showing the greatest deviation. It was also observed that configurations comparable to most detection setups, 10 atm and 6 - 8 cm polyethylene, were closer to the upper bound of deviations from a 252 Cf. 240 Pu, a common spontaneous fission source of interest, ranged from 2% - 3% within this space. This suggests that further investigation is required to determine exact correction factors when calibrating neutron detectors with a 252 Cf. Although demonstrated here for a 4π-detection system, future plans include simulating a more detailed model of the MC-15 detection setup to determine exact correction factors for future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Parametric Studies of Flow Separation using Air Injection

Boundary Layer separation causes the airfoil to stall and therefore imposes dramatic performance degradation on the airfoil. In recent years, flow separation control has been one of the active research areas in the field of aerodynamics due to its promising performance improvements on the lifting device. These active flow separation control techniques include steady and unsteady air injection as well as suction on the airfoil surface etc. This paper will be focusing on the steady and unsteady air injection on the airfoil. Although wind tunnel experiments revealed that the performance improvements on the airfoil using injection techniques, the details of how the key variables such as air injection slot geometry and air injection angle etc impact the effectiveness of flow separation control via air injection has not been studied. A parametric study of both steady and unsteady air injection active flow control will be the main objective for this summer. For steady injection, the key variables include the slot geometry, orientation, spacing, air injection velocity as well as the injection angle. For unsteady injection, the injection frequency will also be investigated. Key metrics such as lift coefficient, drag coefficient, total pressure loss and total injection mass will be used to measure the effectiveness of the control technique. A design of experiments using the Box-Behnken Design is set up in order to determine how each of the variables affects each of the key metrics. Design of experiment is used so that the number of experimental runs will be at minimum and still be able to predict which variables are the key contributors to the responses. The experiments will then be conducted in the 1ft by 1ft wind tunnel according to the design of experiment settings. The data obtained from the experiments will be imported into JMP, statistical software, to generate sets of response surface equations which represent the statistical empirical model for each of the metrics as a function of the key variables. Next, the variables such as the slot geometry can be optimized using the build-in optimizer within JMP. Finally, a wind tunnel testing will be conducted using the optimized slot geometry and other key variables to verify the empirical statistical model. The long term goal for this effort is to assess the impacts of active flow control using air injection at system level as one of the task plan included in the NASAs URETI program with Georgia Institute of Technology.

Zhang, Wei↗