SEARCH · Search NASA

Results for “capability curve”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2020

UFCC Calibration Measurements for UO2 Samples

The uranium UF 6 coincidence counter (UFCC) was designed to measure 5A/B cylinders containing UF6 in the passive neutron self-interrogation mode. The purpose of the present report is to extend the UFCC capability to other sample types such as oxides and liquids that might be found during verification activities. A 252 Cf source was introduced into the detector to make the active mode measurements possible, because only the UF 6 samples have adequate neutron emission for self-interrogation. A customized HDPE insert was designed to position the 252 Cf source under the UO 2 sample cans. The sample enrichments varied from depleted to 93% to establish the shape of the calibration curve and the measurement uncertainties.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Transient Storage Model Parameter Optimization Using the Simulated Annealing Method

Hyporheic exchange in streams is critical to ecosystem functions such as nutrient cycling along river corridors, especially for slowly moving or small stream systems. The transient storage model (TSM) has been widely used for modeling of hyporheic exchange. TSM calibration, for hyporheic exchange, is typically used to estimate four parameters, including the mass exchange rate coefficient, the dispersion coefficient, stream cross-sectional area, and hyporheic zone cross-sectional area. Prior studies have raised concerns regarding the non-uniqueness of the inverse problem for the TSM, that is, the occurrence of different parameter vectors resulting in TSM solution that reproduces the observed in-stream tracer break through curve (BTC) with the same error. This leads to practical non-identifiability in determining the unknown parameter vector values even when global-optimal values exist, and the parameter optimization becomes practically non-unique. To address this problem, we applied the simulated annealing method to calibrate the TSM to BTCs, because it is less susceptible to local minima-induced non-identifiability. A hypothetical (or synthetic) tracer test data set with known parameters was developed to demonstrate the capability of the simulated annealing method to find the global minimum parameter vector, and it identified the “hypothetically-true” global minimum parameter vector even with input data that were modified with up to 10% noise without increasing the number of iterations required for convergence. The simulated annealing TSM was then calibrated using two in-stream tracer tests conducted in East Fork Poplar Creek, Tennessee. In conclusion, simulated annealing was determined to be appropriate for quantifying the TSM parameter vector because of its search capability for the global minimum parameter vector.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1

Abstract. Groundwater plays a key role in meeting water demands, supplying over 40 % of irrigation water globally, with this role likely to grow as water demands and surface water variability increase. A better understanding of the future role of groundwater in meeting sectoral demands requires an integrated hydro-economic evaluation of its cost and availability. Yet substantial gaps remain in our knowledge and modeling capabilities related to groundwater availability, recharge, feasible locations for extraction, extractable volumes, and associated extraction costs, which are essential for large-scale analyses of integrated human–water system scenarios, particularly at the global scale. To address these needs, we developed Superwell, a physics-based groundwater extraction and cost accounting model that operates at sub-annual temporal and at the coarsest 0.5° (≈50 km × 50 km) gridded spatial resolution with global coverage. The model produces location-specific groundwater supply–cost curves that provide the levelized cost to access different quantities of available groundwater. The inputs to Superwell include recent high-resolution hydrogeologic datasets of permeability, porosity, aquifer thickness, depth to water table, recharge, and hydrogeological complexity zones. It also accounts for well capital and maintenance costs, as well as the energy costs required to lift water to the surface. The model employs a Theis-based scheme coupled with an amortization-based cost accounting formulation to simulate groundwater extraction and quantify the cost of groundwater pumping. The result is a spatiotemporally flexible, physically realistic, economics-based model that produces groundwater supply–cost curves. We show examples of these supply–cost curves and the insights that can be derived from them across a set of scenarios designed to explore model outcomes. The supply–cost curves produced by the model show that most (90 %) nonrenewable groundwater in storage globally is extractable at costs lower than USD 0.57 m−3, while half of the volume remains extractable at under USD 0.108 m−3. The global unit cost is estimated to range from a minimum of USD 0.004 m−3 to a maximum of USD 3.971 m−3. We also demonstrate and discuss examples of how these cost curves could be used by linking Superwell's outputs with other models to explore coupled human–environmental system challenges, such as water resources planning and management, or broader analyses of multisectoral feedbacks.

Global Change Analysis Model (GCAM)↗

DOE OSTI2023

The interaction of wholesale electricity market structures under futures with decarbonization policy goals: A complexity conundrum

Competitive wholesale electricity markets can help facilitate energy system decarbonization by incentivizing investments in clean energy technologies that meet evolving system needs. We explore market structure impacts on generator operations and deployment by risk-averse, heterogeneous investor firms using the Electricity Markets and Investment Suite - Agent-based Simulation (EMIS-AS) model. Here we apply clean energy targets of 45%-100% by 2035 considering energy, ancillary services, capacity, and clean energy credit products and pricing and eligibility rules. Results highlight a complexity conundrum, whereby finding the "right" market design to achieve decarbonization goals and avoid unintended consequences can be a highly-nuanced, non-incremental challenge. Carefully designed energy-only markets can achieve the same clean energy targets as capacity market structures but with different revenue and profitability outcomes. Operating reserve demand curve-based scarcity pricing can substitute capacity markets for similar deployment outcomes. Carbon pricing alone is most effective at achieving decarbonization levels at low clean energy targets, and clean energy credit markets and carbon pricing are substitutionary at high clean energy targets. Restricting technology participation in capacity and operating reserve markets can impact deployment and operations, even for nonrestricted technologies. Adding an inertia product with fast frequency response yields insufficient provision at high clean energy targets, but work is needed to understand frequency requirements and capabilities.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

A robust three-parameter reference curve for condensed phase materials

Many isothermal equations of state lead to numerical problems under tension conditions due to negative bulk modulus. We summarize the properties of four existing models and devise a continuous non-piecewise reference curve to address this issue. The final functional form contains three fitting parameters that are directly related to the bulk modulus and its pressure derivative at the reference state and at infinite pressures. We also include bounds that guarantee the monotonicity of its pressure derivative in all volume space. To demonstrate the qualitative features of each isotherm, we conduct calibrations against experimental data for cyclotetramethylene tetranitramine, magnesium oxide, gold, and platinum. Furthermore, pressure residuals and root mean square deviations show that the proposed functional form is capable of adequately fitting the data compared to other three-parameter models.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Relativistic two-component projection-based quantum embedding for open-shell systems

In this work, we present a relativistic quantum embedding formalism capable of variationally treating relativistic effects, including scalar-relativity and spin–orbit coupling. We extend density functional theory (DFT)-in-DFT projection-based quantum embedding to a relativistic two-component formalism, where the full spin magnetization vector form is retained throughout the embedding treatment. To benchmark various relativistic embedding schemes, spin–orbit splitting of the nominally t 2g valence manifold of W(CO) 6 , exchange coupling of [(H 3 N) 4 Cr(OH) 2 Cr(NH 3 ) 4 ] 4+ , and the dissociation potential curve of WF 6 are investigated. The relativistic embedding formalism introduced in this work is well suited for efficient modeling of open-shell systems containing late transition metal, lanthanide, and actinide molecular complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Kinetics of amorphous defect phases measured through ultrafast nanocalorimetry

Recognition of the role of extended defects on local phase transitions has led to the conceptualization of the defect phase, localized thermodynamically stable interfacial states that have since been applied in a myriad of material systems to realize significant enhancements in material properties. Here, we explore the kinetics of grain boundary confined amorphous defect phases, utilizing the high temperature and scanning rates afforded by ultrafast differential scanning calorimetry to apply targeted annealing/quenching treatments at high rates capable of capturing the kinetic behavior. Four Al-based nanocrystalline alloys, including two binary systems, Al–Ni and Al–Y, and two ternary systems, Al–Mg–Y and Al–Ni–Y, are selected to probe the materials design space (enthalpy of mixing, enthalpy of segregation, chemical complexity) for amorphous defect phase formation and stability, with correlative transmission electron microscopy applied to link phase evolution and grain stability to nanocalorimetry signatures. A series of targeted isothermal annealing heat treatments is utilized to construct a Time–Temperature-Transformation curve for the Al–Ni system, from which a critical cooling rate of 2400 °C/s was determined for the grain boundary confined disordered-to-ordered transition. Finally, a thermal profile consisting of 1000 repeated annealing sequences was created to quantify the recovery of the amorphous defect phase following sequential annealing treatments, with results indicating remarkable microstructural stability after annealing at temperatures above 90% of the melting temperature. This work contributes to a deeper understanding of grain boundary localized thermodynamics and kinetics, with potential implications for the design and optimization of advanced materials with enhanced stability and performance.

Amorphous defect phases↗

DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases: Preprint

49 EE - Wind and Water Power Program - Wind (EE-4W↗

DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Constraints on the Physical Properties of GW190814 through Simulations Based on DECam Follow-up Observations by the Dark Energy Survey

On 2019 August 14, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star–black hole merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90% confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline, and we systematically reduced the candidate counterparts through catalog matching, light-curve properties, host-galaxy photometric redshifts, Southern Astrophysical Research spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light-curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Because the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally disrupted neutron star (NS) ejecta properties to quantify our detection efficiency of a counterpart if the merger included an NS. We find that, if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than 10 -8.56 , and has a velocity between 0.18c and 0.21c are disfavored at the 2σ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the 4σ level in 95% of future follow-up observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Phase Transition and Criticality of Methane Confined in Nanopores

For the first time, the phase transition and criticality of methane confined in nanoporous media are measured. The measurement is performed by establishing an experimental setup utilizing a differential scanning calorimeter (DSC) capable of operating under very low temperatures as well as high pressures to detect the capillary phase transition of methane inside nanopores. By performing experiments along isochoric cooling paths, both the capillary condensation and the bulk condensation of methane are detected. Here, the pore critical point (PCP) of nanoconfined methane is also determined and then used to derive the parameters of a previously developed self-consistent equation of state based on the Generalized van der Waals (vdW) partition function. Using this parameter, the equation of state can predict the capillary-condensation curves that agree well with the experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

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↗

DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗

DOE OSTI2024

Validation of a Tapered Impact Damper for Traffic Signal Structure Vibrations Using the Method of Harmonic Balance

Significant vibration amplitudes and cycles can be produced when traffic signal structures with low inherent damping are excited near one of their natural frequencies. For the mitigation of wind-induced vibrations, dynamic vibration absorbers coupled to the structure are often used. Here, this research investigates the performance of a tapered impact damper, consisting of a hanging spring-mass oscillator inside a housing capable of reducing vibration amplitude over a broader frequency range than the conventional tuned mass damper. A nonlinear, two degree-of-freedom model is developed with coordinates representing the traffic structure and the tapered impact damper. This research focuses on the application of the harmonic balance method to approximate the periodic solutions of the nonlinear equations to compute the nonlinear dynamics of the damped traffic signal structure. After designing and manufacturing a tapered impact damper, the traffic signal structure is tested with and without the damper using free vibration snapback tests. The experimental frequency and damping backbone curves are used to validate the analytical model, and the effectiveness of the damper is discussed.

42 ENGINEERING↗

DOE OSTI2025

Angularly resolved spectral reconstruction of x rays via filter pack attenuation

We have designed a new filter pack array to measure angular variations in x-ray spectra during a single shot. The filter pack was composed of repeating identical columns of aluminum and copper filters of varying thicknesses. These columns were located at different positions to measure the spectrum at each corresponding angle. This array was utilized in an experiment to measure the energy evolution of betatron x rays in a laser wakefield accelerator by curving the wakefield with a transverse density gradient, streaking the x rays across the array in front of an x-ray charge-coupled device (CCD) camera. After subtracting the background and “flattening” the image to remove spatial nonuniformities, a critical energy was calculated for each position that produced the best agreement with the measured signal. There was a clear change in critical energy with angle, shedding light on the dynamics of the electrons that traveled through the accelerator. Furthermore, these angles correspond to distinct emission times, covering a timescale of tens of picoseconds. The filter pack was capable of recovering these angular details without the impact of errors introduced by shot-to-shot variability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Novel Calibration Approach for Monitoring Aerosol Hydrogen Isotopes Using Laser-Induced Breakdown Spectroscopy for Molten Salt Reactor Off-Gas Streams

Online monitoring is a key challenge for the continued development of molten salt reactor (MSR) technology. Laser-induced breakdown spectroscopy (LIBS) has previously been demonstrated to be a viable tool for monitoring aerosolized species and noble gases in real time, but the ability to discern varying isotopes in these streams has not yet been investigated for MSR applications. Tritium will form in MSRs from ternary fission and from (n,α)-reactions occurring in lithium-containing salts. This study compares three spectrometers of varying resolutions and types for measuring hydrogen isotope shifts in LIBS spectra of wetted filter paper. For each spectrometer, multivariate models were built (i.e., principal component regression, partial least squares regression, and multivariate curve resolution) to quantify the isotope ratio. The top models were then modified and corrected to apply the models to aerosol samples with varying isotope ratios. This novel calibration strategy offers an 82% reduction in volume of the calibration samples needed and is a more viable pathway for calibrating deployable LIBS systems. Lastly, this calibration model was compared with an all-aerosol trained model for monitoring hydrogen isotopes during a real-time test where the protium/deuterium ratio, along with representative salt species (i.e., lithium, sodium, and potassium) were adjusted dynamically. Results of this test validated the predictive capabilities of the transferred model and highlighted the capabilities of LIBS for real-time monitoring of MSR effluent streams.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.4 is a program that obtains a steady-state flow and temperature solution, in geometrical extent, for a nuclear reactor core, or a single fuel assembly, including the effect of manufacturing and modelling uncertainties by the means of hot channel factors. It is based on an evolutionary sequence of codes originally used for calculating plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over the past 30 years. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset of nucleate boiling ratio (ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability for select design and safety parameters is available to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗