SEARCH · Engineering Papers

Results for “Transient Modeling”

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 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Improvements to the Modeling of the TREAT Reactor and Experiments

This paper summarizes the latest improvements and lessons learned from the modeling and simulation of the transient test reactor at Idaho National Laboratory using the MAM- MOTH reactor physics application. MAMMOTH is a MOOSE-based, Finite Element Method application that specializes in the analysis of the spatial dynamics behavior of nuclear reactors. Since early 2018 several transient tests have been conducted at TREAT, thus providing the opportunity to apply and benchmark modern modeling and simulation tools. MAMMOTH was used to provide predictions of the power coupling factor between the core and the experiment for various experiments. Even though the power coupling factor predictions agree very well with the experimental data, within the bounds of the experimental uncertainty, one shortcoming was the underprediction of the total energy deposited in the core and in the experiment. Determination of the sources for this discrepancy are ongoing, but several key problems have been identified and resolved, thus providing valuable insights for future research. This paper discusses several of these lessons learned. First, the heat capacity data for the TREAT fuel has some significant problems due to limitations of the measurement techniques used circa 1960s. The sensitivity of the peak power and the total energy deposition to various representations of the heat capacity is approximately 5%. Second, the effects of the biological shield and thermal column on the modeling of the core are non-negligible, since they affect the mean generation time and the effective reflection of neutrons back into the core, which is suspected to be important during the core heat up. Matching the reactor period resolves the fact that the reduced spatial domain used in the MAMMOTH model underpredicts the mean generation time. The neutron reflection from these regions is marginally improved with the use of an albedo boundary condition. Third, modeling of the control rod movement with a multi-scheme method is introduced and its current limitations are exposed. Fourth, we explore the effects of using a homogenized model with Superhomogenization equivalence and how that differs from fully heterogeneous simulations. Finally, the energy condensation effects for this graphite core are significant. Solutions with 10 and 26 energy groups show the benefits of using a finer coarse group structure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Studying the Impact of IBR Modeling on the Commonly Applied Transmission Line Protective Elements

It is a well-accepted fact that detailed inverter-based resource (IBR) models are needed for electromagnetic transient (EMT) studies, especially for protection studies and protection designs. However, it is not clear which modeling aspects should be included in the EMT model to accurately represent the IBR fault responses. Therefore, this paper addresses this gap by evaluating the impact of key modeling aspects (DC source, inverter model, power control, current/voltage control, and current limiting) on the response of protective relay elements both qualitatively and quantitatively. Both grid-following (GFL) and grid-forming (GFM) IBR models are comprehensively evaluated through a simplified real-world transmission network by varying the IBR modeling aspects, fault locations, and types. The insights and findings from this paper will benefit the power system protection community because this work provides modeling requirements and guidelines for IBR modeling so that correct IBR models will be used for protection scheme settings design.

grid-following inverter↗

NASA NTRSDate not supplied

Dependence of Lightning Cloud-Top Optical Emissions on Channel Characteristics Using Boltzmann Radiative Transfer Modeling

With an estimated global flash rate of 46 lightning per second, a significant amount of lightning nitrogen oxides (LNOx) are produced within the climate system. In turn, LNOx affect the concentrations of ozone (O 3 , an important greenhouse gas) as well as the very reactive hydroxyl radicals (OH, that affect all greenhouse gases). Thus, space-based estimation of LNOx is highly sought [Koshak 2017; Koshak et al. (2015)]. Total lightning channel length, L , is a proxy for the total flash energy (and hence LNOx production). Geostationary lightning mappers provide cloud-top 2D views of the diffuse (i.e., blurry) optical images of lightning, thereby making it difficult to determine L (and other channel features). Transient Boltzmann radiative transfer modeling is applied here to investigate the dependence of the spatio-temporal pattern of cloud-top lightning optical emissions on L and other channel characteristics (e.g., orientation, propagation direction, amplitude, location, ...).

William J Koshak↗

NASA NTRS1992

Coupled loads analysis for Space Shuttle payloads

Described here is a method for determining the transient response of, and the resultant loads in, a system exposed to predicted external forces. In this case, the system consists of four racks mounted on the inside of a space station resource node module (SSRNMO) which is mounted in the payload bay of the space shuttle. The predicted external forces are forcing functions which envelope worst case forces applied to the shuttle during liftoff and landing. This analysis, called a coupled loads analysis, is used to couple the payload and shuttle models together, determine the transient response of the system, and then recover payload loads, payload accelerations, and payload to shuttle interface forces.

Eldridge, J.↗

DOE OSTI2022

Validation of URANS and STRUCT-ε turbulence models for stratified sodium flow

Simulations of a transient stratified sodium experiment are carried out using a classic unsteady RANS model and the second-generation URANS model, STRUCT-ε. Turbulence modeling challenges and their implications to stratified flow prediction are discussed in the context of other sources of error. Input errors are discussed and addressed; discretization error is calculated to be less than 5% of the inlet velocity, for 80% of the domain; and remaining errors in temperature distributions are attributed to the turbulence model. Qualitative flow features from the simulations are presented and discussed. Compared to the experiment, the STRUCT-ε turbulence model provides a more physically accurate prediction of temperature and momentum mixing in key regions of the domain. Quantitative measures such as the L 2 norm of the temperature discrepancy demonstrate the improved performance of the STRUCT-ε approach. In conclusion, the magnitude of the temperature fluctuations is very well-predicted by the STRUCT-ε, while URANS overpredicts them by approximately 50%.

42 ENGINEERING↗

DOE OSTI2026

Isopentane Disproportionation in Lewis Acidic Chloroaluminate Ionic Liquid

Lewis acidic chloroaluminate ionic liquid diluted in dichloromethane catalyzes the disproportionation of alkanes, such as isopentane, via carbenium ions in a reaction readily initiated by carbenium ion precursors such as tert-butyl chloride (TBC). The carbenium ion-AlCl 4 – ion-pairs stabilized by the ionic liquid-dichloromethane solution are the key intermediates in two distinctive kinetic regimes, i.e., a transient regime (0–5 min) and a steady-state regime (after 5 min). The transient regime constitutes the majority of isopentane conversion and is governed by the initial carbenium ion concentration. In the steady-state regime, disproportionation occurs at a considerably lower rate, affected by the carbenium ion concentration, the concentration of the ionic liquid, and the reaction temperature. The formation of alkenes observed in the 1 H NMR spectra of the reacting substrates, along with the DFT calculations, suggests that deprotonation of carbenium ion-pairs reduces their concentration, decreasing, in turn, the reaction rate. Kinetic modeling indicates that the transient regime is significantly controlled by the hydride transfer (k HT ) and the deprotonation rate constants (k DP ), while the steady-state regime is additionally influenced by the alkene protonation rate constant (k –DP ). The overall activation energy of the reaction at the steady state, expressed as E a,steady-state regime = E a,HT – E a,DP + E a,–DP , was 54 kJ/mol. The reaction mechanism and the kinetics highlight the potential of Lewis acid-catalyzed conversions of hydrocarbons under remarkably mild conditions.

Deprotonation↗

DOE OSTI2024

In-cylinder spray evolution in a motored central-injection gasoline engine: Imaging and simulating the effects of flash-boiling and intake crossflow

Accurate predictions of fuel spray behavior and mixture formation in simulations of direct-injection spark-ignition (DISI) engines are fundamental to ensure proper description of all subsequent processes including ignition, combustion, and emissions. In this work, the spray evolution in a single-cylinder optical DISI engine was studied experimentally and numerically with the goal of enabling predictive computational fluid dynamics (CFD) modeling of in-cylinder sprays. The authors explored a wide range of operating conditions characterized by several fuel injection temperatures and engine speeds, using a well-characterized nine-component gasoline surrogate known as PACE-20. The effect of flash boiling and intake crossflow on the spray is discussed, with a focus on evaluating the ability of the spray models to capture highly transient spray behavior. In the experiments, the fuel temperature was varied between 20°C and 80°C, allowing for non-flash- to flash-boiling transition to emerge with enhanced flashing intensity at the highest temperatures. Spray collapse resulted in vapor-rich regions, owing to the locally lower inertia of the fluid. Varying the engine speed from 650 to 1950 rpm promoted increasingly more turbulent in-cylinder crossflow which interacted with the spray during the injection event and resulted in enhanced spray dispersion. The CFD model was able to capture the spray morphology transition at different fuel temperatures and engine speeds adequately. Further, it is shown that the spray breakup model could capture the transitional spray behavior induced by flash boiling atomization and intake flow via proper initialization of the spray cone angle and calibration of the spray models’ constants.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗

DOE OSTI2023

Dynamic systems modeling of the spallation neutron source cryogenic moderator system to optimize transient control and prepare for power upgrades

Through support of the US Department of Energy's Office of Basic Energy Sciences, Oak Ridge National Laboratory has begun applying machine learning methods to improve accelerator and target performance of the Spallation Neutron Source (SNS). These methods are being applied to the control optimization and power upgrade of the SNS Cryogenic Moderator System (CMS). A numerical model of the CMS has been developed to study these optimizations and system modifications using EcosimPro. This paper compares steady-state and transient numerical results with experimental data. Control optimization studies focused on dampening mass flow, temperature, and pressure fluctuations during sudden losses of accelerator beam power. This analysis was conducted by adjusting five proportional-integral-derivative controllers connected to four flow control valves and one heater. Future efforts include power uprate studies focused on increasing the CMS cooling capacity. The current accelerator power is 1.4 MW; the first target station is being upgraded to 2.0 MW as part of the Proton Power Upgrade effort. The CMS cooling capacity is sufficient for 2.0 MW operation.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Signature Analysis Utilizing a Dynamic Molten Salt Reactor Model for MC&A

Moving from traditional fixed fuel nuclear reactor systems to a mobile, dynamic fuel system that has dissolved special nuclear material in a molten salt is a paradigm shift in several respects. One major consideration is how to develop effective nuclear material controls and accounting for these novel reactor systems. The Molten Salt Reactor Experiment is one example of a critical molten salt system and provided a significant reference library based on the documented effort. But it was low thermal power that was not intended to reflect a commercial scale electricity production design. Therefore, to facilitate and assist vendors with domestic licensing considerations, research is underway to identify methods that could be used for domestic safeguards approaches for these novel reactor systems. This research presents the results of a signature analysis of data generated from three simulated scenarios using the molten salt demonstration reactor model defined in the Transient Simulation Framework of Reconfigurable Modules. Each scenario provides 1 hour isotope inventories over a 180 day period. The scenarios investigated provide test cases to examine if direct gamma-ray spectroscopy of the fuel can be used to identify changes comparing a base case (no reactivity control and fixed fission contribution) to an insertion of reactivity (10 pcm no change in fission composition) and a change in fission composition. The analysis demonstrates that monitoring the total count rate in a highly collimated high-resolution photon energy spectrum is sensitive to perturbations imposed into the reactor model. The total photon count rate changes ≈2% for the fission composition change and ≈4.5% for the reactivity insertion compared to the base case. However, both scenarios show an increase in the total photon count rate. The total photon count rate can be used to identify changes due to power (number of fissions) but not due to a change in the material undergoing fission. To distinguish between the two cases of increased power, the photon spectrum would require an intensive analysis technique. A photon peak count rate ratio analysis could be used to identify changes in the fissile material fission generation in the core through identification of a static peak that shows little variation to the source of fission and a highly varying peak. The photon peak strength will ultimately be determined by the isotope’s fission yield. A preliminary analysis investigating the coupling of the isotope’s fission yield and its concentration in the fuel salt derived from the modeling has been performed. A ratio analysis of the photon peak count rates of 140 La to 99 Mo demonstrated that the reactivity insertion creates a distinct difference in the ratio compared to the fission composition change scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1973

Rheological effects on friction in elastohydrodynamic lubrication

An analytical and experimental investigation is presented of the friction in a rolling and sliding elastohydrodynamic lubricated contact. The rheological behavior of the lubricant is described in terms of two viscoelastic models. These models represent the separate effects of non-Newtonian behavior and the transient response of the fluid. A unified description of the non-Newtonian shear rate dependence of the viscosity is presented as a new hyperbolic liquid model. The transient response of viscosity, following the rapid pressure rise encountered in the contact, is described by a compressional viscoelastic model of the volume response of a liquid to an applied pressure step. The resulting momentum and energy equations are solved by an iterative numerical technique, and a friction coefficient is calculated. The experimental study was performed, with two synthetic paraffinic lubricants, to verify the friction predictions of the analysis. The values of friction coefficient from theory and experiment are in close agreement.

Trachman, E. G.↗

DOE OSTI2022

Implementation and testing of physics-based pulverization model in BISON

This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRSDate not supplied

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗

NASA NTRS1981

Closed-form solutions of performability

Methods which yield closed form performability solutions for continuous valued variables are developed. The models are similar to those employed in performance modeling (i.e., Markovian queueing models) but are extended so as to account for variations in structure due to faults. In particular, the modeling of a degradable buffer/multiprocessor system is considered whose performance Y is the (normalized) average throughput rate realized during a bounded interval of time. To avoid known difficulties associated with exact transient solutions, an approximate decomposition of the model is employed permitting certain submodels to be solved in equilibrium. These solutions are then incorporated in a model with fewer transient states and by solving the latter, a closed form solution of the system's performability is obtained. In conclusion, some applications of this solution are discussed and illustrated, including an example of design optimization.

Meyer, J. F.↗

NASA NTRS2002

Optimal Charging of Nickel-Hydrogen Batteries for Life Extension

We are exploring the possibility of extending the cycle life of battery systems by using a charging profile that minimizes cell damage. Only nickel-hydrogen cells are discussed at this time, but applications to lithium-ion cells are being considered. The process first requires the development of a fractional calculus based nonlinear dynamic model of the specific cells being used. The parameters of this model are determined from the cell transient responses. To extend cell cycle life, an instantaneous damage rate model is developed. The model is based on cycle life data and is highly dependent on cell voltage. Once both the cell dynamic model and the instantaneous damage rate model have been determined, the charging profile for a specific cell is determined by numerical optimization. Results concerning the percentage life extension for different charging strategies are presented. The overall procedure is readily adaptable to real-time implementations where the charging profile can maintain its minimum damage nature as the specific cell ages.

Hartley, Tom T.↗

DOE OSTI2023

Comprehensive Efficiency Analysis of Seven-Switch Current-Source Inverter Based on Voltage-Boost Function and Genetic Algorithm

This paper investigates the voltage-boost features of the emerging seven-switch current-source inverter (H7-CSI) to improve the efficiency of the CSI-based motor drive system. The switching transient and power loss model of the H7-CSI topology during voltage-boost operation are analyzed in detail. The effects of the boost function on the switching and conduction loss of the H7-CSI topology are investigated. A genetic algorithm has been introduced that makes it possible to optimize the H7-CSI’s modulation index to minimize its power loss. Experimental results for the H7-CSI unit show that this optimization can reduce the power loss by 10.6%.

42 ENGINEERING↗

NASA NTRS2003

Investigations of Mechanisms Associated with Nucleate Boiling Under Microgravity Conditions

In this work a building block type of approach is used so that a basic understanding of the processes that contribute to nucleate boiling heat fluxes under microgravity conditions can be developed. This understanding will lead to development of a mechanistic model for nucleate boiling heat transfer which could eventually be used as a design tool in space applications. Task Description Task 1: Fabrication of the Experimental Setup. Under this task, the test section and liquid holding and viewing chambers will be fabricated. Artificial cylinder cavities will be formed on silicon wafers. A single cavity and two or four cavities, with a prescribed spacing and size, will be formed. The desired nucleation wall superheat will be used to determine the size of the mouth of the cavities. Task 2: Experiments. The basic experiments for flow and temperature field around single and multiple (2 or 4 separated or merged bubbles growing on downward facing or inclined surfaces) will be carried out under normal gravity conditions. The experiments will be conducted at one atmosphere pressure, but liquid subcooling will be varied from 0 to 30C. Water and PF-5050 will be used as test liquids. Task 3: Analytical/Numerical Models. In this task, transient temperature and flow field in vapor and liquid will be determined during growth of a single bubble. Analysis will include the evolution of the vapor-liquid interface and development of microlayer underneath the bubbles. For merged bubbles, detailed calculations of flow and temperature field will be carried out for transient shapes of vapor stems supporting a large bubble and the corresponding evaporation rate. Flow and temperature field for a bubble sliding along a heated wall will also be determined. Microgravity conditions will be simulated and a framework of a numerical tool for prediction of nucleate boiling heat fluxes under microgravity conditions will be developed. Task 4: Experiments in a KC-135. To understand bubble growth and detachment behavior of single or large merged bubbles, boiling experiments will be conducted under low gravity (10-2 g) conditions of the aircraft. In these experiments, 'designed' surfaces will be used. Visual observations and heat transfer data will be taken, but holography will not be used. The apparatus used for laboratory experiments will also be employed for experiments in the aircraft. Task 5: Experiments in the Space Shuttle. Effort will be devoted for defining a boiling experiment to be conducted on a 'designed' surface. The experiment will provide microgravity data on bubble growth and departure. These data are needed for development of a credible model for nucleate boiling heat fluxes under microgravity conditions. The heat transfer data will also be obtained and will be used to validate the models.

Dhir, Vijay K.↗

NASA NTRS1977

Current status of droplet and liquid combustion

The present understanding of spray combustion in rocket engine, gas turbine, Diesel engine and industrial furnace applications is reviewed. In some cases, spray combustion can be modeled by ignoring the details of spray evaporation and treating the system in the same manner as a gaseous diffusion flame; however, in many circumstances, this type of simplification is not adequate and the turbulent two-phase flow must be considered. The behavior of individual droplets is a necessary component of two-phase models and recent work on transient droplet evaporation, ignition and combustion is considered, along with a discussion of important simplifying assumptions involved with modeling these processes. Methods of modeling spray evaporation and combustion processes are also discussed including: one-dimensional models for rocket engine and prevaporized combustion systems, lumped zone models (utilizing well-stirred reactor and plug flow regions) for gas turbine and furnace systems, locally homogeneous turbulent models, and two-phase models. The review highlights the need for improved injector characterization methods, more information of droplet transport characteristics in turbulent flow and continued development of more complete two-phase turbulent models.

Faeth, G. M.↗