SEARCH · Engineering Papers

Results for “Building thermal dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

helios: An R package to process heating and cooling degrees for GCAM

helios is an open-source R package that estimates population-weighted heating and cooling degree-hours (HDH and CDH) and degree-days (HDD and CDD) at various temporal (e.g., energy dispatch segments, monthly, yearly) and spatial scales (e.g., U.S. states, global political regions, countries). The degree hour and degree day outputs from helios are used to inform electricity demand load in the Global Change Analysis Model (GCAM) as well as in GCAM-USA (which is the version of GCAM with U.S. state-level details). helios uses a workflow with four steps: processing raw data; calculating heating and cooling degrees; visualizing performance diagnostics; and outputing results in various formats. There are two sources of widely-used climate data compatible with helios: (1) hourly climate data with 12-km resolution that are dynamically downscaled with the Weather Research and Forecasting (WRF) model and projected using a thermal global warming (TGW) approach; and (2) daily climate data with 0.5-degree resolution from the Coupled Model Intercomparison Project (CMIP) that is bias-adjusted and statistical downscaled by the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP). In summary, helios is a model that standardizes methodology of heating and cooling degrees-hours and degree-days using publicly available data and advance the understanding of the impact of spatial and temporal temperature variability on building energy services.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Multi-objective sizing and dispatch for building thermal and battery storage towards economic and environmental synergy

The role of building thermal and battery storage is pivotal in advancing smart cities and achieving sustainability goals through effective energy management. Despite their significance, there are several limitations in the sizing approach and value stream analysis with various objectives for their widespread adoption in buildings. This work proposes a flexible and scalable multi-objective optimization framework for optimal sizing and dispatch of building thermal and battery storage, addressing conflicting objectives simultaneously using mixed-integer linear programming. The weighted-sum method is adapted, combining multiple objectives into a single function. The two-stage procedure iterates over different weights, generating optimal solutions and forming the Pareto front. Case studies are performed to assess the energy, economic, and environmental benefits of building energy storage systems for a large office building in three climate locations. The results demonstrate that the proposed framework efficiently determines optimal sizing and dispatch strategies, addressing the balance between economic viability and emission reduction. The dynamic relationship between time-of-use energy charges and emission factors leads to significantly different strategies based on whether economic or environmental concerns are prioritized. This research enhances our understanding of the benefits of TES and BES systems in buildings, providing valuable guidance to stakeholders.

25 ENERGY STORAGE↗

DOE OSTI2025

Experimental test of model predictive control in a variable air volume system

Model predictive control (MPC) has been widely studied as a promising approach for improving energy efficiency and operational flexibility in buildings, yet its real-world performance for commercial variable air volume (VAV) systems remains insufficiently characterized. In particular, the impacts of model mismatch on control robustness, real-time computational burden, and device-level operation are rarely evaluated using long-term field data. Here, this study presents a comprehensive experimental evaluation of MPC applied to a full-scale VAV system in Oak Ridge National Laboratory’s Flexible Research Platform-2 building with constant cooling/heating temperature setpoints and no occupancy. The study offers three key advantages over existing work: (1) it uses a representative building in a full-scale experimental test, capturing realistic system dynamics and complexity; (2) it evaluates a relatively sophisticated MPC formulation using two different optimization solvers (Gurobi and PSO), fully accounting for computational complexity and methodological diversity; and (3) it systematically assesses potential negative impacts on various building devices, benchmark against a well-established baseline, ASHRAE Guideline 36 (G36). To isolate zone- and air-handling-unit–level supervisory control effects, the supply fan was operated with a fixed static pressure setpoint under all strategies, and the trim-and-response static pressure reset in G36 was not enabled. Results show that MPC maintained thermal comfort while improving energy efficiency. Abrupt solar radiation variations degraded performance. Computation times ranged from ∼1 s (Gurobi) to ∼ 70 s (PSO). Compared with G36, MPC achieves 33% energy savings and reduces median reheat coil output by approximately a factor of 5–10 for a representative cooling day under matched weather conditions. However, it increases the maximum discomfort deviation from 0.5 to 1°C and results in a 32% increase in staging frequency. In addition, PSO-based MPC introduced damper oscillations, also affecting actuator longevity.

ASHRAE guideline 36↗

NASA NTRS2013

Human and Robotic Mission to Small Bodies: Mapping, Planning and Exploration

This study investigates the requirements, performs a gap analysis and makes a set of recommendations for mapping products and exploration tools required to support operations and scientific discovery for near- term and future NASA missions to small bodies. The mapping products and their requirements are based on the analysis of current mission scenarios (rendezvous, docking, and sample return) and recommendations made by the NEA Users Team (NUT) in the framework of human exploration. The mapping products that sat- isfy operational, scienti c, and public outreach goals include topography, images, albedo, gravity, mass, density, subsurface radar, mineralogical and thermal maps. The gap analysis points to a need for incremental generation of mapping products from low (flyby) to high-resolution data needed for anchoring and docking, real-time spatial data processing for hazard avoidance and astronaut or robot localization in low gravity, high dynamic environments, and motivates a standard for coordinate reference systems capable of describing irregular body shapes. Another aspect investigated in this study is the set of requirements and the gap analysis for exploration tools that support visualization and simulation of operational conditions including soil interactions, environment dynamics, and communications coverage. Building robust, usable data sets and visualisation/simulation tools is the best way for mission designers and simulators to make correct decisions for future missions. In the near term, it is the most useful way to begin building capabilities for small body exploration without needing to commit to specific mission architectures.

Neffian, Ara V.↗

DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗

DOE OSTI2021

Heat and mass transfer during a sudden loss of vacuum in a liquid helium cooled tube - Part III: Heat deposition in He II

A sudden loss of vacuum can be catastrophic for superconducting particle accelerators. In such an event, air leaks into the liquid-helium-cooled accelerator beamline tube and condenses on its inner surface, causing rapid boiling of the helium and dangerous pressure build-up. Understanding the coupled heat and mass transfer processes is important for the design of the beamline cryogenic system. Our past experimental study on nitrogen gas propagating in a copper tube cooled by normal liquid helium (He I) has revealed a nearly exponential slowing down of the gas front. A theoretical model that accounts for the interplay of the gas dynamics and the condensation was developed, which successfully reproduced various key observations. However, since many accelerator beamlines are actually cooled by superfluid helium (He II) in which the heat transfer is via a non-classical thermal counterflow mode, we need to extend our work to the He II cooled tube. This paper reports our systematic measurements using He II and the numerical simulations based on a modified model that accounts for the He II heat-transfer characteristics. By tuning the He II peak heat-flux parameter in our model, we have reproduced the observed gas dynamics in all experimental runs. The fine-tuned model is then utilized to reliably evaluate the heat deposition in He II. Finallly, this work not only advances our understanding of condensing gas dynamics but also has practical implications to the design codes for beamline safety.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Site demonstration and performance evaluation of MPC for a large chiller plant with TES for renewable energy integration and grid decarbonization

Thermal energy storage (TES) for a cooling plant is a crucial resource for load flexibility. Traditionally, simple, heuristic control approaches, such as the storage priority control which charges TES during the nighttime and discharges during the daytime, have been widely used in practice, and shown reasonable performance in the past benefiting both the grid and the end-users such as buildings and district energy systems. However, the increasing penetration of renewables changes the situation, exposing the grid to a growing duck curve, which encourages the consumption of more energy in the daytime, and volatile renewable generation which requires dynamic planning. The growing pressure of diminishing greenhouse gas emissions also increases the complexity of cooling TES plant operations as different control strategies may apply to optimize operations for energy cost or carbon emissions. This paper presents a model predictive control (MPC), site demonstration and evaluation results of optimal operation of a chiller plant, TES and behind-meter photovoltaics for a campus-level district cooling system. The MPC was formulated as a mixed-integer linear program for better numerical and control properties. Compared with baseline rule-based controls, the MPC results show reductions of the excess PV power by around 25%, of the greenhouse gas emission by 10%, and of peak electricity demand by 10%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRSDate not supplied

Americium-Fuelled Radioisotope Stirling Generator Elegant Breadboard

High power solutions for spacecraft are becoming increasingly sought after in the new era of space exploration, with energy intensive in-situ lunar surface exploration becoming a key part of the roadmap for major space agencies, as well as commercial companies. Nuclear power provides the most feasible solution to achieve continued operation through the long lunar night, and in extended periods in Permanently Shadowed Regions (PSRs). Americium-241 (241Am) is used as the fuel of choice in European Radioisotope Power Systems (RPS). The University of Leicester have developed a number of systems using this fuel type, including a 200 Wth heat source. Building on existing international partnerships, the University of Leicester is continuing to work with NASA Glenn Research Center (GRC) to maximise the electrical output that can be generated from the heat source by progressing the radioisotope-Stirling generator concept to elegant breadboard level, building on the long heritage of dynamic power conversion developments at NASA Glenn. Presented in the following conference paper is the next iteration of the Americium-fuelled Radioisotope Stirling Generator (Am-RSG) Elegant Breadboard, following on from a successful technology demonstration design conducted at NASA Glenn in early 2025. The design laid-out in this paper is a more flight-like version of the technology demonstrator, and is intended to pair with the next generation of Stirling convertor technology, the SunPower Robust Stirling Converter (SRSC). The SRSC is able to operate more efficiently at lower thermal power levels than the previous Advanced Stirling Convertor (ASC), making it more suitable for the 241Am system.

Lunar↗

DOE OSTI2021

Fluid dynamic simulation and analysis of water-cooling systems for the Electron-Ion Collider

The Electron-Ion Collider is the newest large-scale project at Brookhaven National Laboratory. The collider’s purpose is to provide further advancements in the knowledge of the universe’s origin by accelerating particles near the speed of light. Our project for this 3.8 km ring was to create a thermal hydraulic steady-state simulation design of the water-cooling system to be cost-effective and energy efficient, as envisioned by Charlie Foltz, the EIC Infrastructure Division Director. The system would include a supply and return header, which cools several thousand components of the ring. The water would then be returned and cooled down using a system of cooling towers and plate and frame heat exchangers. Due to the size of the system and the complexity of the network analysis, a fluid dynamic simulation software, AFT Fathom, was used. Since previous methods of maintaining systems relied on building upon smaller real-life models and implementing empirical data, this flow model was unique and first of a kind in the domain of accelerator design, construction and operation. Therefore, our hydraulic team piloted a new method to perform network analysis on a large scale cooling system. We successfully created several test scenarios for system behavior in a shorter time compared to the method of performing hand calculations. Cooling specifications for heat rejection, pressure drop, flow rate, and pipe sizing were changed based on the individual systems of the vacuum, radio frequency (RF), magnet and power supply, and cryogenics sections. Finally, we used DOE guidelines to perform life-cycle cost analysis with net present value and carbon saving analysis on the systems where pipe size could be optimized.

43 PARTICLE ACCELERATORS↗

NASA NTRSDate not supplied

High Efficiency Megawatt Motor Stator Thermal Performance

Enabling single aisle electric aircraft propulsion requires power dense, megawatt scale, high efficiency electric machines. To that end NASA has been developing the High Efficiency Megawatt Motor (HEMM) 1.46 MW (16kW/kg), 98% efficient electric machine. The success of this effort is highly dependent on the stator’s thermal (cooling) design. The design to date has been based upon computation fluid dynamic and finite element analysis models that have been validated with testing of stator sub sections that replicate expected thermal conditions. However, the test to validate the model can only be so accurate without fully representing the full geometry and components of the stator. A full HEMM stator has been fabricated, potted, and installed in a housing complete with vacuum tube. This apparatus was used to validate the full thermal environment of HEMM and act as a final validation of the stator design before fully fabricating the machine. This paper discusses the modeling and test results from this penultimate HEMM stator build.

electric machine↗

NASA NTRS2001

EUNIS: An Extreme-Ultraviolet Normal-Incidence Spectrometer

GSFC is in the process of assembling an Extreme-Ultraviolet Normal-Incidence Spectrometer called EUNIS, to be flown as a sounding rocket payload. This instrument builds on the many technical innovations pioneered by our highly successful SERTS experiment over its past ten flights. The new design will have somewhat improved spatial and spectral resolutions, as well as 100 times greater sensitivity, permitting EUV spectroscopy with a temporal resolution near 1-second for the first time ever. To achieve such high time cadence, a novel Active-Pixel-Sensor detector is being developed as a key component of our design. The high sensitivity of EUNIS will allow entirely new studies of transient coronal phenomena, such as the rapid loop dynamics seen by TRACE, and searches for non-thermal motions indicative of magnetic reconnection or wave heating. The increased sensitivity will also permit useful EUV spectra at heights of 2-3-R$ \odot$ above the limb, where the transition between the static corona and the solar wind might occur. In addition, the new design features two independent optical systems, more than doubling the spectral bandwidth covered on each flight. Its 300-370\AA\ bandpass includes He-II 304\AA\ and strong lines from Fe-XI-XVI, extending the current SERTS range of 300-355\AA\ to further improve our ongoing series of calibration under-flights for SOHO/CDS and EIT. The second bandpass of 170-230\AA\ has a sequence of very strong Fe-IX-XIV lines, and will allow under-flight support for two more channels on SOHO/EIT, two channels on TRACE, one on Solar-B/EIS, and all four channels on the STEREO/EUVI instrument. First flight of the new EUNIS payload is scheduled for 2002 October.

Thomas, Roger J.↗

DOE OSTI2021

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2005

Evaluation of an Active Clearance Control System Concept

Reducing blade tip clearances through active tip clearance control in the high pressure turbine can lead to significant reductions in emissions and specific fuel consumption as well as dramatic improvements in operating efficiency and increased service life. Current engines employ scheduled cooling of the outer case flanges to reduce high pressure turbine tip clearances during cruise conditions. These systems have relatively slow response and do not use clearance measurement, thereby forcing cold build clearances to set the minimum clearances at extreme operating conditions (e.g., takeoff, reburst) and not allowing cruise clearances to be minimized due to the possibility of throttle transients (e.g., step change in altitude). In an effort to improve upon current thermal methods, a first generation mechanically-actuated active clearance control (ACC) system has been designed and fabricated. The system utilizes independent actuators, a segmented shroud structure, and clearance measurement feedback to provide fast and precise active clearance control throughout engine operation. Ambient temperature performance tests of this first generation ACC system assessed individual seal component leakage rates and both static and dynamic overall system leakage rates. The ability of the nine electric stepper motors to control the position of the seal carriers in both open- and closed-loop control modes for single and multiple cycles was investigated. The ability of the system to follow simulated engine clearance transients in closed-loop mode showed the system was able to track clearances to within a tight tolerance (0.001 in. error).

Steinetz, Bruce M.↗

NASA NTRS2005

Evaluation of an Active Clearance Control System Concept

Steinetz, Bruce M.↗

NASA NTRS2016

Predictive Evaluations of Oxygen-Rich Hydrocarbon Combustion Gas-Centered Swirl Coaxial Injectors using a Flamelet-Based 3-D CFD Simulation Approach

NASA Marshall Space Flight Center (MSFC) has embarked upon a joint project with the Air Force to improve the state-of-the-art of space application combustion device design and operational understanding. One goal of the project is to design, build and hot-fire test a 40,000 pound-thrust Oxygen/Rocket Propellant-2 (RP-2) Oxygen-Rich staged engine at MSFC. The overall project goals afford the opportunity to test multiple different injector designs and experimentally evaluate the any effect on the engine performance and combustion dynamics. To maximize the available test resources and benefits, pre-test, combusting flow, Computational Fluid Dynamics (CFD) analysis was performed on the individual injectors to guide the design. The results of the CFD analysis were used to design the injectors for specific, targeted fluid dynamic features and the analysis results also provided some predictive input for acoustic and thermal analysis of the main Thrust Chamber Assembly (TCA). MSFC has developed and demonstrated the ability to utilize a computationally efficient, flamelet-based combustion model to guide the pre-test design of single-element Gas Centered Swirl Coaxial (GCSC) injectors. Previous, Oxygen/RP-2 simulation models utilizing the Loci-STREAM flow solver, were validated using single injector test data from the EC-1 Air Force test facility. The simulation effort herein is an extension of the validated, CFD driven, single-injector design approach applied to single injectors which will be part of a larger engine array. Time-accurate, Three-Dimensional, CFD simulations were performed for five different classes of injector geometries. Simulations were performed to guide the design of the injector to achieve a variety of intended performance goals. For example, two GCSC injectors were designed to achieve stable hydrodynamic behavior of the propellant circuits while providing the largest thermal margin possible within the design envelope. While another injector was designed to purposefully create a hydrodynamic instability in the fuel supply circuit as predicted by the CFD analysis. Future multi-injector analysis and testing will indicate what if any changes occur in the predicted behavior for the single-element injector when the same injector geometry is placed in a multi-element array.

Richardson, Brian R.↗

NASA NTRS2019

Maturation of Dynamic Power Convertors for Radioisotope Power Systems

Dynamic Radioisotope Power Systems (DRPS) are being developed by NASA's RPS Program in collaboration with the U.S. Department of Energy (DOE) for space science and exploration missions. A development effort is currently underway to mature dynamic power convertors for infusion into a potential future flight generator. This convertor maturation effort was formulated by the RPS Program at NASA Headquarters based on successful maturation models and utilizes expertise from agency technology and mission centers to support requirements development and technology assessments. The effort is being executed by Glenn Research Center's Dynamic Radioisotope Power Systems (DRPS) project and the Thermal Energy Conversion Branch. The Dynamic Power Convertor (DPC) contracts consist of three phases to enable design, fabrication, and independent assessment after delivery to the government. The contracts are intended to gather data on candidate dynamic conversion technologies to fill knowledge gaps, support assessments of dynamic conversion technologies, and elicit generator requirements. The SRG-110 (Stirling Radioisotope Generator) and ASRG (Advanced Stirling Radioisotope Generator) flight development projects provided Stirling convertor demonstration units and engineering models to verify and validate convertors against performance specifications and mission requirements. These units utilize temperature resistance materials and non-contacting bearings to demonstrate wear-free, long-life operation. This maturation effort builds on past lessons learned and new requirements focused on demonstrating convertor robustness to critical environments meant to stress key aspects of each convertor within the margins of the design. This effort seeks to realize the full potential of dynamic power conversion technologies for NASA's space science and exploration missions.

Wilson, Scott↗

NASA NTRS1989

The transient performance of a two-phase fluid reservoir

Thermal control of future large, high power spacecraft will require a two-phase fluid central bus. The two-phase fluid reservoir is a critical component in the two-phase fluid bus. It both controls the saturation temperature and provides a space for volumetric changes. A dynamic reservoir simulation model does not currently exist, but it is needed to expedite efforts and reduce risk. During 1989 an effort was made to develop a simulation model of the transient performance of a two-phase fluid reservoir. As a beginning, a preliminary model was developed. It is based upon component mathematical models in lumped parametric form and build upon five component mathematical models for calculating dynamic responses of two-phase fluid reservoirs, primary feedback elements, controller commands, heater actuators, and reservoir heaters. As much as possible, the model took advantage of the available SINDA'85/FLUINT thermal/fluid integrator. Additional calculation logic and computer subroutines were developed to complete implementation of the model. The model is capable of simulating dynamic response of an equilibrium two-phase fluid reservoir. Modification of the model to include the liquid/vapor nonequilibrium is required for applications of the model to simulate performance of reservoir in which the liquid and vapor phases of the reservoir fluid are not in equilibrium. In addition, the model in its present form, needs to be refined in several respects. More empirical data are needed to guide the model development. The model may then be used to conduct a full parametric study of two-phase fluid reservoirs. More complexities in two-phaes flow regions in laboratory and flight conditions may have to be considered eventually if empirical data cannot be simulated satisfactorily. System with other components arrangement also need to be simulated if optimization is ever to be attained. The present model does, however, preliminarily demonstrates that such analyses are quite possible and offers a far less expensive method to understand the transient of a two-phase fluid reservoir system than a totally headware approach.

Chi, Joseph↗