Search NASA⌕ Search

SEARCH · Search NASA

Results for “fast time simulation”

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 649 records · Page 36

Simulated Encounters of the Parker Solar Probe with a Coronal-Hole Jet

Solar coronal jets are small, transient, collimated ejections most easily observed in coronal holes (CHs). The upcoming Parker Solar Probe (PSP) mission provides the first opportunity to encounter CH jets in situ near the Sun and examine their internal structure and dynamics. Using projected mission orbital parameters, we have simulated PSP encounters with a fully three-dimensional magnetohydrodynamic (MHD) model of a CH jet. We find that three internal jet regions, featuring different wave modes and levels of compressibility, have distinct identifying signatures detectable by PSP. The leading Alfvén wave front and its immediate wake are characterized by transAlfvénic plasma flows with mild density enhancements. This front exhibits characteristics of a fast switch-on MHD shock, whose arrival is signaled by the sudden onset of large-amplitude transverse velocity and magnetic-field oscillations highly correlated in space and time. The trailing portion is characterized by supersonic but subAlfvénicout flows of dense plasma with uncorrelated velocity and magnetic-field oscillations. This compressible region contains most of the jet's mass. The volume between the immediate wake and dense jet, the remote wake,mixes and transitions the characteristics of the two other regions. In addition to probing each region separately, we also simulate a corotational PSP-jet encounter. In this scenario, the simulated spacecraft hovers over the jet producing CH, as may occur during the mission's corotational phases, sampling each jet region in turn. We estimate that PSP will encounter numerous CH jets over the lifetime of the mission.

Roberts, Merrill A.↗

Machine Protection Studies of Fast Instability-Driven Beam Losses in the EIC Electron Storage Ring

Fast transverse beam instabilities represent one of the most challenging protection scenarios for the Electron-Ion Collider (EIC) electron storage ring (ESR). Coherent betatron oscillations driven by collective effects can grow rapidly and lead to substantial beam losses within only a few tens of turns, posing a threat to collimators, superconducting (SC) magnets, cryogenic systems, and detector components. A dedicated simulation framework based on distributed transverse excitations was developed to evaluate instability-driven beam losses in the ESR. A conservative instability rise time on the order of 10 turns was used to define a bounding protection scenario. Simulations show that the original two-collimator configuration is insufficient to localize losses under all operating conditions, whereas an optimized four-collimator configuration reduces loss leakage into protected regions by more than two orders of magnitude. Nevertheless, the most severe instability scenarios produce catastrophic beam impacts, with nearly half of the stored beam being intercepted by the collimators within a single turn, exceeding the preliminary estimated survivability limits of the baseline collimator design. The results demonstrate that collimation alone is insufficient to guarantee safe ESR operation and must be complemented by a fast machine-protection system. Dedicated turn-by-turn beam-position monitors and fast beam-loss monitors provide sufficient advance warning to detect instability growth and initiate a beam abort before equipment-damaging loss levels are reached. These findings establish quantitative requirements for the ESR machine-protection architecture and support the implementation of a multilayer protection strategy for both accelerator and detector systems.

43 PARTICLE ACCELERATORS↗

Micro–macro finite element modeling method for rub response in abradable coating materials

Gas turbine engines experience “rub” when the rotating blades come in contact with a static abradable coating. This results in extreme strain rates and dynamics inside a high-temperature/high-pressure environment. Current rub models are phenomenological and do not reflect the underlying microstructures, thus limiting their prediction accuracy. In this work, a microstructure-informed, reduced order modeling framework is introduced for simulating abradable coating “rub" behavior. This framework comprises a microscale model constructed based on digitized abradable microstructure and explicitly simulates the mechanical behavior of each constituent phases and their interactions. After calibration and validation with experiment data, the calibrated microscale model is used to generate data across a vast range of applied strain rates and temperature with various load paths. Then, the virtually generated data are used to fit the macroscopic-reduced order model, which enables fast component scale rub simulation without compromising the integrity of the complex material behavior. In conclusion, the proposed effort will address the technical challenge of predicting abradable material behavior during rub through the application of multiscale modeling from microstructure to engines behavior, effectively reducing the development costs and time of new abradable material for better “rub” properties.

36 MATERIALS SCIENCE↗

99 TcO 4 - removal from legacy defense nuclear waste by an alkaline-stable 2D cationic metal organic framework

Removal of 99 TcO 4 - from legacy defense nuclear tank waste at Savannah River Site is highly desirable for the purpose of nuclear safety and environmental protection, but currently not achievable given the extreme conditions including high alkalinity, high ionic strength, and strong radiation field. Herein, we present a potential solution to this long-term issue by developing a two-dimensional cationic metal organic framework SCU-103, showing ultrahigh stability in alkaline aqueous media and great resistance to both β and γ radiation. More importantly, it is very effective for 99 TcO 4 - separation from aqueous media as demonstrated by fast exchange kinetics, high sorption capacity, and superior selectivity, leading to the successful removal of 99 TcO 4 - from actual Savannah River Site high level tank waste for the first time, to the best of our knowledge. In addition, the uptake mechanism is comprehensively elucidated by molecular dynamics simulation and density functional theory calculation, showing a unique chemical recognition of anions with low charge density.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Individualized empirical baselines for evaluating the energy performance of existing buildings

The evaluation of building energy performance requires a baseline for comparison. Common empirical baselines are usually used for existing buildings since they are fast and convenient. However, the same type of building at the same location will receive the same baseline despite their difference in usage. Individualized baselines by creating building energy models are possible solutions, but it is labor intensive and time-consuming. To fill the gap, this study is to develop individualized empirical baselines for existing buildings in a fast way. First, common empirical baselines are created based on survey data. Then, to get training samples, building energy models for large-scale existing buildings are created and simulated. So finally, based on simulation results, mathematical models to get individualized empirical baselines in a fast way are created. U.S. medium office buildings were used as an example to demonstrate the method. We developed 30 mathematical models for medium office buildings in two vintages (constructed before 1980 and after 1980) and 15 climate zones. The mean absolute percentage errors (MAPE) between the individualized empirical baselines and the modeled baselines for those 30 mathematical models are all lower than 5.5%. An engineer can obtain the individualized empirical baseline for an existing building in a few seconds by using the open-source tool we developed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗

Modeling of flash sintering of ionic ceramics

A fundamental understanding of the influence of defects in ionic ceramics at the atomic, microstructural, and macroscopic levels, before, during, and after the flash sintering event is key to the development of ceramic processing operations that lead to fast, low cost, and environmentally safe fabrication of materials. The observed phenomenology of the flash process encompasses multiple time and length scales and have resulted on a wide variety of what sometimes appears to be contradictory explanations. Here, this article summarizes the latest developments on the modeling and simulation of flash sintering, specifically those related to the understanding of the equilibrium and kinetic properties and the corresponding microstructural evolution of ionic ceramics. Challenges and opportunities in the development of theoretical analyses that include unidentified multiphysical effects are discussed, as they pertain to the processing of technologically relevant ceramic materials for advanced structures and devices.

36 MATERIALS SCIENCE↗

Characterization of DESI fiber assignment incompleteness effect on 2-point clustering and mitigation methods for DR1 analysis

We present an in-depth analysis of the fiber assignment incompleteness in the Dark EnergySpectroscopic Instrument (DESI) Data Release 1 (DR1). This incompleteness is caused by therestricted mobility of the robotic fiber positioner in the DESI focal plane, which limits thenumber of galaxies that can be observed at the same time, especially at small angular separations.As a result, the observed clustering amplitude is suppressed in a scale-dependent manner, which,if not addressed, can severely impact the inference of cosmological parameters. We discuss themethods adopted for simulating fiber assignment on mocks and data. In particular, we introducethe fast fiber assignment (FFA) emulator, which was employed to obtain the power spectrumcovariance adopted for the DR1 full-shape analysis. We present the mitigation techniques,organised in two classes: measurement stage and model stage. We then use high fidelity mocks as areference to quantify both the accuracy of the FFA emulator and the effectiveness of the differentmeasurement-stage mitigation techniques. This complements the studies conducted in a parallelpaper for the model-stage techniques, namely the θ-cut approach. We find that pairwiseinverse probability (PIP) weights with angular upweighting recover the “true” clustering in allthe cases considered, in both Fourier and configuration space. Notably, we present the first everpower spectrum measurement with PIP weights from real data.

79 ASTRONOMY AND ASTROPHYSICS↗

Fast Whole-Engine Stirling Analysis

An experimentally validated approach is described for fast axisymmetric Stirling engine simulations. These simulations include the entire displacer interior and demonstrate it is possible to model a complete engine cycle in less than an hour. The focus of this effort was to demonstrate it is possible to produce useful Stirling engine performance results in a time-frame short enough to impact design decisions. The combination of utilizing the latest 64-bit Opteron computer processors, fiber-optical Myrinet communications, dynamic meshing, and across zone partitioning has enabled solution times at least 240 times faster than previous attempts at simulating the axisymmetric Stirling engine. A comparison of the multidimensional results, calibrated one-dimensional results, and known experimental results is shown. This preliminary comparison demonstrates that axisymmetric simulations can be very accurate, but more work remains to improve the simulations through such means as modifying the thermal equilibrium regenerator models, adding fluid-structure interactions, including radiation effects, and incorporating mechanodynamics.

Dyson, Rodger W.↗

Fast Whole-Engine Stirling Analysis

An experimentally validated approach is described for fast axisymmetric Stirling engine simulations. These simulations include the entire displacer interior and demonstrate it is possible to model a complete engine cycle in less than an hour. The focus of this effort was to demonstrate it is possible to produce useful Stirling engine performance results in a time-frame short enough to impact design decisions. The combination of utilizing the latest 64-bit Opteron computer processors, fiber-optical Myrinet communications, dynamic meshing, and across zone partitioning has enabled solution times at least 240 times faster than previous attempts at simulating the axisymmetric Stirling engine. A comparison of the multidimensional results, calibrated one-dimensional results, and known experimental results is shown. This preliminary comparison demonstrates that axisymmetric simulations can be very accurate, but more work remains to improve the simulations through such means as modifying the thermal equilibrium regenerator models, adding fluid-structure interactions, including radiation effects, and incorporating mechanodynamics.

Dyson, Rodger W.↗

A Framework to Analyze the Requirements of a Multiport Megawatt-Level Charging Station for Heavy-Duty Electric Vehicles

Widespread adoption of heavy-duty (HD) electric vehicles (EVs) will soon necessitate the use of megawatt (MW)-scale charging stations to charge high-capacity HD EV battery packs. Such a station design needs to anticipate possible station traffic, average and peak power demand, and charging/wait time targets to improve throughput and maximize revenue-generating operations. High-power direct current charging is an attractive candidate for MW-scale charging stations at the time of this study, but there are no precedents for such a station design for HD vehicles. We present a modeling and data analysis framework to elucidate the dependencies of a MW-scale station operation on vehicle traffic data and station design parameters and how that impacts vehicle electrification. This framework integrates an agent-based charging station model with vehicle schedules obtained through real-world vehicle telemetry data analysis to explore the station design and operation space. A case study applies this framework to a Class 8 vehicle telemetry dataset and uses Monte Carlo simulations to explore various design considerations for MW-scale charging stations and EV battery technologies. The results show a direct correlation between optimal charging station placement and major traffic corridors such as cities with ports, e.g., Los Angeles and Oakland. Corresponding parametric sweeps reveal that while good quality of service can be achieved with a mix of 1.2-megawatt and 100-kilowatt chargers, the resultant fast charging time of 35–40 min will need higher charging power to reach parity with refueling times.

33 ADVANCED PROPULSION SYSTEMS↗

Parallel spatial direct numerical simulations on the Intel iPSC/860 hypercube

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube is documented. The direct numerical simulation approach is used to compute spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows. The feasibility of using the PSDNS on the hypercube to perform transition studies is examined. The results indicate that the direct numerical simulation approach can effectively be parallelized on a distributed-memory parallel machine. By increasing the number of processors nearly ideal linear speedups are achieved with nonoptimized routines; slower than linear speedups are achieved with optimized (machine dependent library) routines. This slower than linear speedup results because the Fast Fourier Transform (FFT) routine dominates the computational cost and because the routine indicates less than ideal speedups. However with the machine-dependent routines the total computational cost decreases by a factor of 4 to 5 compared with standard FORTRAN routines. The computational cost increases linearly with spanwise wall-normal and streamwise grid refinements. The hypercube with 32 processors was estimated to require approximately twice the amount of Cray supercomputer single processor time to complete a comparable simulation; however it is estimated that a subgrid-scale model which reduces the required number of grid points and becomes a large-eddy simulation (PSLES) would reduce the computational cost and memory requirements by a factor of 10 over the PSDNS. This PSLES implementation would enable transition simulations on the hypercube at a reasonable computational cost.

Joslin, Ronald D.↗

Transition in a Supersonic Boundary Layer Due to Acoustic Disturbances

The boundary layer receptivity process due to the interaction of three-dimensional slow and fast acoustic disturbances with a blunted flat plate is numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 10(exp 6)/inch. The computations are performed with and without two-dimensional isolated roughness element located near the leading edge. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The simulations showed that the linear instability waves are generated very close to the leading edge. The wavelength of the disturbances inside the boundary layer first increases gradually and becomes longer than the wavelength for the instability waves within a short distance from the leading edge. The wavelength then decreases gradually and merges with the wavelength for the Tollmien-Schlichting wave. The initial amplitudes of the instability waves near the neutral points, the receptivity coefficients, are about 1.20 and 0.07 times the amplitude of the free-stream disturbances for the slow and the fast waves respectively. It was also revealed that small isolated roughness element does not enhance the receptivity process for the given nose bluntness.

Balakumar, P.↗

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

The Energy & Environmental Research Center (EERC), in partnership with the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the EERC’s State Energy Research Center (SERC), MPLX Operations LLC, and TC Energy Development Holdings Inc. (a subsidiary of TC Energy Corporation), studied the potential for subsurface hydrogen storage and recovery in the Williston Basin of western North Dakota. The project’s goal was to evaluate the feasibility of large-scale, secure geologic H 2 storage to support future hydrogen generation, storage, and use. This work included laboratory testing, H 2 –rock–fluid exposure experiments, literature reviews on H 2 embrittlement, and reservoir modeling and simulations. The study included an assessment of storage potential across three types of storage reservoirs using both reservoir simulation and DOE’s web-based tool SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration – Hydrogen Estimator for Logistical Planning), as well as investigation of potential H 2 production and markets for commercial-scale deployment. Building on prior EERC gas storage research, three storage options were selected for detailed evaluation: the Broom Creek Formation (a clastic saline reservoir), the Dickinson Lodgepole Mounds (DLM) complex (carbonate mud mound structures) of the Lodgepole Formation (an active oil and gas producing reservoir), and the Dunham Salt Interval of the Piper Formation (to be used for engineered salt cavern development). These targets were prioritized based on prior EERC research using datasets related to seal capacity, reservoir quality, mechanical integrity, and injectivity. Exposure tests on Broom Creek and DLM samples showed mineral dissolution and precipitation that increased brine salinity and altered reservoir rock surfaces. Although these results provide useful insight, they are limited by small sample sizes and short-term (30-day) exposure, requiring further study to assess long-term storage integrity. Salt formations were not tested because of their known nonreactivity and established mechanical stability. Results of reservoir simulations performed for a single site demonstrated that the Broom Creek Formation may be capable of receiving up to 42,000 tonnes of injected H 2 over 7 months via one well. H 2 recovery took place over 5 months, resulting in approximately 26,000 tonnes (~62% without cushion gas [CG]). This work suggests water production may be important and subsequent cycles of injection and production may perform more efficiently; however, significant site-specific work in the future is needed to assess actual reservoir performance of injection and withdrawal of H 2 storage. For oil reservoir potential, a multiple-well model was used to simulate injection of approximately 32,000 tonnes of H 2 into a single wellbore while simultaneously producing in place reservoir fluids from four offset wells to maintain reservoir pressure. The simulation results suggested a high recovery (~98%); in addition, cost advantages through existing infrastructure could be realized. Challenges in this reservoir include vi managing gas purity and leakage risks. In both scenarios, production of H 2 takes place in a single-well scenario with 10 cycles (7 months of injection and 5 months of production) over 10 years. Finally, the use of engineered caverns in the Dunhan Salt was evaluated, and the results suggest that while they have a smaller capacity (<1000 tonnes per cavern), they exhibit nearly complete gas recovery (>99%), fast response times, and low purity risk. While caverns in North Dakota may be smaller in capacity, fields can be developed in galleries to accommodate the volumetric needs and rapid turnaround times necessary to meet market demands. Geographic limitations and thin salt intervals in North Dakota may represent less total storage potential than salt domes elsewhere, but significant opportunities exist to expand this market for gas storage in North Dakota. A basinwide assessment was performed to estimate a first-of-its-kind value for H 2 storage on a large scale. DOE’s SHASTA-HELP, combined with EERC simulation work, was used to perform the assessment. Estimated H 2 storage potential varied widely for each formation type. The Broom Creek saline formation was estimated to have a storage potential of approximately 1.7–90.5 million tonnes (MMt). The DLM oil reservoirs were estimated to have 0.07–0.19 MMt of capacity. Notably, each of these estimates relies on significant assumptions regarding reservoir thickness, porosity, permeability, and CG needed for operation. Much research is needed to understand the true site-specific storage resource potential of each formation. Using the Dunham Salt Interval for cavern development may result in as much as 4.79 MMt (up to 2.87 MMt working gas) of H 2 storage potential. An important note for consideration is that the values presented here need significantly more geological characterization and engineering assessments prior to gaining confidence in performance. This will be a focal point for future research and development needs. The basinwide evaluation also indicated that North Dakota has significant H 2 generation potential, with estimates up to a possible 13 MMt annually, suggesting a substantial opportunity for H 2 market development and thus the need for commercial-scale H 2 storage to facilitate growth and resilience. Key Recommendations 1. Conduct detailed site characterization (3D seismic, well logs, core sampling) to reduce geologic uncertainty. 2. Perform techno-economic analyses incorporating market, regulatory, and incentive frameworks. 3. Investigate long-term interactions among H 2 , CGs, well materials, and formations to assess risks. 4. Develop pilot- and field-scale demonstrations to validate models and establish best practices.

03 NATURAL GAS↗

The dependence of confinement on the isotope mass in the core and the edge of AUG and JET-ILW H-mode plasmas

Experiments in ASDEX Upgrade (AUG) and JET with the ITER-like wall (JET-ILW) are performed to separate the pedestal and core contributions to confinement in H-modes with different main ion masses. A strong isotope mass dependence in the pedestal is found which is enhanced at high gas puffing. This is because the ELM type changes when going from D to H for matched engineering parameters, which is likely due to differences in the inter ELM transport with isotope mass. The pedestal can be matched in H and D plasmas by varying only the triangularity and keeping the engineering parameters relevant for core transport the same. With matched pedestals Astra/TGLF (Sat1geo) core transport simulations predict the experimental profiles equally well for H and D. These core transport simulations show a negligible mass dependence and no gyro-Bohm scaling is observed. However, to match the experimental observations at medium β it is required to take the fast-ion dilution and rotation into account. This is not enough for high β plasmas where for the first time a profile match between H and D plasmas was achieved experimentally. Under these conditions quasilinear modelling with TGLF over predicts the transport in the core of H and D plasmas alike.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Applications of a fast and efficient method for simulating orbital mission safety

A fast and efficient method for simulating potential interactions between space objects is presented. The specialized hazard analysis simulation is run at the Russian Space Surveillance Center. It simulates the motion of objects in space and determines close encounters between these objects. First, ephemerides are predicted for all space objects over the time interval of concern. Second, using the predicted ephemerides, pairs of satellites that would likely collide are grouped together. Throughout the process, the pairs initially chosen are refined until all pairs of satellites are found which come within a certain miss distance of each other. In addition to the miss distance, other characteristics of the encounter including angle of incidence, relative velocity, type of approaching object, and probability of collision come out of the simulation. There are many potential applications of this simulation and its results to missing planning and mission safety. Artificial objects can be introduced to the simulation to examine potential encounters between the new objects and current objects for mission design applications. The simulation can be run many times to create data for statistical analysis of encounter parameters. These as well as many other potential applications are presented and discussed.

Huth, Jeffrey↗

Isotopic mass effects of tritium-fueled high-performance TFTR supershots

An increase in total stored energy correlated with the addition of tritium fuel was observed in supershots during the TFTR DT campaign. This supershot regime had strikingly high, centrally peaked ion and electron temperatures, and the largest neutron emission rates observed in TFTR. This paper presents a study of the causes of this increase in stored energy in supershots. Twenty-six supershots have been recently reanalyzed with the TRANSP plasma analysis code. Early TRANSP simulations did not accurately match the measured magnitude and time evolutions of the neutron emission rates. This mismatch is attributed to neglecting apparent increases of trace amounts of heavy impurities during neutral beam injection. The new TRANSP runs were tuned to accommodate this and match the measured global neutron emission rates. These new runs also had improved fidelity in predicting the time histories and radial dependencies of measured DT neutron emission rates. That in turn adds confidence in the simulated thermal deuterium and tritium density profiles that are needed for calculating the average hydrogenic atomic mass profiles. Six subsets of these supershots had well matched toroidal field B tor , plasma current I p , flux geometry, and total injected neutral beam power. The mix of D and T beam ions was varied for different discharges. The magnitude of the increase of the thermal ion energy W i with added tritium was relatively small, and the total W tot increase is dominated by the increase in fast beam ions with T. Analyses at times before the occurrences of deleterious MHD instabilities yielded scaling of W tot with the volume-average isotopic mass consistent with previous publications. The relative fraction of fast energy ions is expected to be small in practical tokamak reactors. Thus the increase in stored energy W tot observed in TFTR supershots does not appear likely to be significantly helpful for producing useful fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mean Line Pump Flow Model in Rocket Engine System Simulation

A mean line pump flow modeling method has been developed to provide a fast capability for modeling turbopumps of rocket engines. Based on this method, a mean line pump flow code PUMPA has been written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The pump code can model axial flow inducers, mixed-flow and centrifugal pumps. The code can model multistage pumps in series. The code features rapid input setup and computer run time, and is an effective analysis and conceptual design tool. The map generation capability of the code provides the map information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of the code permit parametric design space exploration of candidate pump configurations and provide pump performance data for engine system evaluation. The PUMPA code has been integrated with the Numerical Propulsion System Simulation (NPSS) code and an expander rocket engine system has been simulated. The mean line pump flow code runs as an integral part of the NPSS rocket engine system simulation and provides key pump performance information directly to the system model at all operating conditions.

Veres, Joseph P.↗