Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid water path”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 451 records · Page 25

Elastic wave analyzer for icy sub-surfaces (EWAIS) in the solar system

The sub-surface structures of water-bearing icy worlds (e.g., Europa, Enceladus, and Titan) and the structures of icy regions on Mars contain critical information about their origin and evolution. These bodies have been identified as high-priority targets in the NRC Planetary Science Decadal Survey, which noted a lack of technology readiness as one impediment to related missions. An Elastic Wave Analyzer for Icy Sub-surfaces (EWAIS) instrument is being developed to acquire data that can provide answers to fundamental science questions related to these bodies. The EWAIS instrument generates and receives elastic waves using an array of piezoelectric transducers capable of operating at temperatures as low as 30K. A developed breadboard acquires reflected signals from discontinuities in the traveled materials’ impedance (the product of wave velocity and density) along the wave path. Analysis of the acquired data can be used to determine the elastic properties, presence of cracks, locations of cavities and other discontinuities, as well as the thicknesses of ice and liquid layers. The EWAIS instrument is being developed with a novel dual-frequency transmitter array that will enable reaching distances of kilometers through ice, as well as resolving smaller acoustic reflectors nearer the EWAIS instrument. The selected dual frequencies correspond to two depth ranges that can be analyzed, with resolutions of tens of centimeters in the region of tens to hundreds of meters and tens of meters on the kilometers scale, respectively. The EWAIS will be applicable for operation from any in-situ platform, including surface assets (e.g., lander and rover) and sub-surface assets (e.g., melting probes).

Bar-Cohen, Y.↗

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗

Assessing parallel path cooling tower performance via artificial neural networks

Real-time monitoring of a research nuclear reactor, a system in which all generated power is dissipated to the environment, can be performed via analysis of the heat rejection from the cooling system. Given an inlet water temperature and flow rate, the reactor power can be well-approximated from the outlet water temperature; however, the instrumentation to measure outlet conditions may not be robust or accurate. If we know how a cooling tower performs from historical data, but cannot measure the outlet temperature, a mathematical representation of the system can be inverted to obtain the outlet water temperature that describes the cooling capacity. Unfortunately, model inversion processes are computationally expensive. To address this, an artificial neural network (ANN) is implemented to assess the performance of a multi-cell cooling tower for a nuclear reactor. This approach leverages the Merkel model to obtain an extensive data set describing performance of the cooling tower cells throughout a wide array of potential operating conditions. The Merkel model is expressed as a function of four parameters: the inlet and outlet water temperatures, inlet air wet bulb temperature, and ratio of liquid-to-gas mass flow rates (L/G), which together provide a non-dimensional number indicative of cooling tower performance, called the Merkel integral. Computing a 4-dimensional data structure that describes finite combinations of the Merkel integral, an inverse model is then generated using an ANN to determine the cell outlet water temperature from the other three model parameters along with the computed Merkel integral. Compared to traditional model inversion methods, the ANN reduces the computational time by approximately 4 orders of magnitude, with effectively no sacrifice to solution accuracy, and could be applied for different cooling towers in the event the performance curve is known. Finally, three use cases of the ANN are then reviewed: (1) determining the cell outlet water temperatures when gas flow at rated conditions (GFRC) is known, (2) performing the prior case without knowledge of the GRFC, and (3) assessing performance differences between the individual tower cells.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Beyond Kepler: Direct Imaging of Earth-like Planets

Is there another Earth out there? Is there life on it? People have been asking these questions for over two thousand years, and we finally stand on the verge of answering them. The Kepler space telescope is NASA's first mission designed to study Earthlike exoplanets (exo-Earths), and it will soon tell us how often exo-Earths occur in the habitable zones of their stars. The next natural step after Kepler is spectroscopic characterization of exo-Earths, which would tell us whether they possess an atmosphere, oxygen, liquid water, as well as other biomarkers. In order to do this, directly imaging an exo-Earth may be necessary (at least for Sun-like stars). Directly imaging an exo-Earth is challenging and likely requires a flagship-size optical space telescope with an unprecedented imaging system capable of achieving contrasts of 1(exp 10) very close to the diffraction limit. Several coronagraphs and external occulters have been proposed to meet this challenge and are in development. After first overviewing the history and current state of the field, my talk will focus on the work proceeding at the Ames Coronagraph Experiment (ACE) at the NASA Ames Research Center, where we are developing the Phase Induced Amplitude Apodization (PIAA) coronagraph in a collaboration with JPL. PIAA is a powerful technique with demonstrated aggressive performance that defines the state of the art at small inner working angles. At ACE, we have achieved contrasts of 2(exp -8) with an inner working angle of 2 lambda/D and 1(exp -6) at 1.4 lambda/D. On the path to exo-Earth imaging, we are also pursuing a smaller telescope concept called EXCEDE (EXoplanetary Circumstellar Environments and Disk Explorer), which was recently selected for technology development (Category III) by NASA's Explorer program. EXCEDE will do fundamental science on debris disks as well as serve as a technological and scientific pathfinder for an exo-Earth imaging mission.

Belikov, Ruslan↗

Development and experimental verification of an active noise cancellation (ANC) method for non-invasive sensing

Ultrasonic wave propagation may be used in non-invasive sensing to identify the properties of a fluid in a container. An ultrasonic excitation can be coupled into the liquid or the gas through the walls of the container and detected by a receiver on the other side. However, in addition to this fluid path, another signal will also propagate through the structure and is referred as noise in this study. For most materials, this structural noise will be of higher amplitude compared to the fluid signal due to the large acoustic impedance difference between the solid container and the fluid. In this paper, an active noise cancellation (ANC) method is developed and demonstrated on an experimental setup. It was shown that this method cancels more than 90% of the noise signal in a given time window of interest. A sequential version of this cancellation method is also proposed, and it was shown that its performance improves as the number of canceling excitations is increased. The performance was found to be independent of the location of canceling transducers and position of the silence window but dependent on the width of the silence window. As an example, the ANC method was applied to noise cancellation of a signal travelling through water in a stainless-steel container. The signal-to-noise ratio (SNR) of the water signal in the window of interest was increased by 148% from 13.98 dB to 34.68 dB by using the proposed ANC method.

42 ENGINEERING↗

Venus as a Laboratory for Exoplanetary Science

The current goals of the astrobiology community are focused on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments (surface liquid water) and detectable signatures thereof is the study of where the boundaries of such environments can occur. Such studies provide the basis for understanding how a once inhabitable planet might come to be uninhabitable. The archetype of such a planet is arguably Earth’s sibling planet, Venus. Given the need to define the conditions that can rule out bio-related signatures of exoplanets, Venus provides a unique opportunity to explore the processes that led to a completely uninhabitable environment by our current definition of the term. Here we review the current state of knowledge regarding Venus, particularly in the context of remote-sensing techniques that are being or will be employed in the search for and characterization of exoplanets. We discuss candidate Venus analogs identified by the Kepler and TESS exoplanet missions and provide an update to exoplanet demographics that can be placed in the potential runaway greenhouse regime where Venus analogs are thought to reside. We list several major outstanding questions regarding the Venus environment and the relevance of those questions to understanding the atmospheres and interior structure of exoplanets. Finally, we outline the path towards a deeper analysis of our sibling planet and the synergy to exoplanetary science.

Venus↗

Drop Interactions with the Conical Shock Structure Generated by a Mach 4.5 Projectile

This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of $192$ $\mu$m. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier–Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant $C$ 2 = $2$ compared to the accepted value $C$ 2 = $2/3$. In conclusion, results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Geophysical Characterization and Monitoring of the 200 Area Treated Effluent Disposal Facility to Support Permit Renewal; TEDF Geophysical Monitoring

The Treated Effluent Disposal Facility (TEDF), located in the 200 East Area of the Hanford Site, is a site where non-hazardous and non-radioactive liquid wastes are disposed into two state-permitted infiltration basins. In 2016, the Washington State Department of Ecology denied a permit renewal request for TEDF due to the inability to adequately assess the impact of TEDF discharge water on the underlying groundwater quality. TEDF overlies the relatively impermeable Ringold Lower Mud (RLM) unit, whose upper contact lies in the vadose zone approximately 30 m below ground surface and approximately 10 m above the water table boundary. The RLM is assumed to isolate TEDF discharge water from the natural groundwater aquifer, which is monitored using wellbores screened below the RLM. Therefore, samples collected from monitoring wells near TEDF are not considered representative of TEDF discharge water. Rather, TEDF discharge water is assumed to mound on top of the RLM to form a perched aquifer. To support permitting of the TEDF, a new monitoring well is required that can be used to sample water from the presumed perched water aquifer above the RLM. Ideally, the screened section of the well would be located at the peak of the perched water mound(s), which presumably occurs at the point of maximum vertical flux from TEDF to the RLM, or equivalently where the dominant infiltration flow paths reach the RLM. This report describes how time-lapse 3D electrical resistivity tomography (ERT) was used in conjunction with nominal TEDF discharge operations to image the dominant flow paths from each pond to the RLM. Results are summarized in Figure ES.1. Figure ES.1A shows a satellite image of the TEDF overlain by an array of surface ERT electrodes. The solid and dashed black circles denote the zones of maximum vertical flux at the RLM within the south and north ponds, respectively, and presumably the regions where perched water peaks during discharge. Figure ES.1B shows time-lapse difference imaging results approximately 15 days after switching discharge from the south pond to the north pond. Blue iso-surfaces beneath the north pond denote regions of increasing bulk electrical conductivity caused by increasing saturation due to infiltrating water. Red iso-surfaces beneath the south pond denote regions that were previously saturated during the south pond discharge and are now de-saturating, causing a decrease in bulk electrical conductivity. In both cases, the zones of maximum change mark the dominant flow paths to the RLM. Figure ES.1C shows time-lapse imaging results approximately 15 days after switching discharge from the north pond to the south pond. In this case, blue iso-surfaces mark the dominant flow paths to the RLM from the south pond. The red iso-surfaces mark the dominant flow paths that existed during discharge to the north pond. The regions of maximum change in bulk conductivity (due to saturation or desaturation) that mark the primary flow paths are equivalent in both cases and denoted by the dashed solid and black circles. If a perched water zone forms on the RLM, it is likely to mound within or near the dashed circle during discharge to the north pond, and in the solid circle during discharge to the south pond. In other words, if perched water mounding occurs, the circles mark the optimum locations suggested by the ERT imaging for monitoring boreholes to be placed, enabling samples to be collected that are representative of TEDF discharge impacts on groundwater quality.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Numerical modeling and performance analysis of a membrane-based air dehumidifier using ionic liquid desiccant

Membrane-based Liquid Desiccant Air Dehumidifier (MLDAD) has the potential to overcome the shortcomings of the conventional open-tower liquid desiccant air dehumidifier. The MLDAD uses membranes with high water vapor permeability to separate air stream from the liquid desiccant and thus eliminates carry-over and the resulting corrosion issues. Recently, the ionic liquid desiccant becomes a promising alternative to conventional liquid desiccant. The ionic liquid desiccant has a large potential of dehumidification, which is also non-corrosive to metals and non-crystallizable. A two-dimensional numerical heat and mass transfer model of the MLDAD using a recently identified ionic liquid desiccant, [EMIM]OAc, is presented here. This model can simulate the performance of the MLDAD with various designs, including different selections of liquid desiccant and membrane materials, dimensions of the MLDAD, and flow patterns. For porous membranes, this model accounts for several micro-scale mass transfer mechanisms of vapor transportation across the membrane. It can also model the mass transfer performance of dehumidifiers using non-porous membranes based on the experimentally measured membrane permeance. The numerical model was validated against the performance data available from literature and experimental tests. The maximum discrepancy of the latent effectiveness between the measured and the model-predicted results was about 6% when porous membranes were used. A parametric study was conducted with the numerical model. Finally, the results indicated that the membrane permeability, air path height, and solution to air mass ratio were the most critical parameters, determining the dehumidification performance of the MLDAD.

42 ENGINEERING↗

Numerical modeling and performance analysis of a membrane-based air dehumidifier using ionic liquid desiccant

Membrane-based Liquid Desiccant Air Dehumidifier (MLDAD) has the potential to overcome the shortcomings of the conventional open-tower liquid desiccant air dehumidifier. The MLDAD uses membranes with high water vapor permeability to separate air stream from the liquid desiccant and thus eliminates carry-over and the resulting corrosion issues. Recently, the ionic liquid desiccant becomes a promising alternative to conventional liquid desiccant. The ionic liquid desiccant has a large potential of dehumidification, which is also non-corrosive to metals and non-crystallizable. A two-dimensional numerical heat and mass transfer model of the MLDAD using a recently identified ionic liquid desiccant, [EMIM]OAc, is presented in this paper. We report this model can simulate the performance of the MLDAD with various designs, including different selections of liquid desiccant and membrane materials, dimensions of the MLDAD, and flow patterns. For porous membranes, this model accounts for several micro-scale mass transfer mechanisms of vapor transportation across the membrane. It can also model the mass transfer performance of dehumidifiers using non-porous membranes based on the experimentally measured membrane permeance. The numerical model was validated against the performance data available from literature and experimental tests. The maximum discrepancy of the latent effectiveness between the measured and the model-predicted results was about 6% when porous membranes were used. A parametric study was conducted with the numerical model. The results indicated that the membrane permeability, air path height, and solution to air mass ratio were the most critical parameters, determining the dehumidification performance of the MLDAD.

42 ENGINEERING↗

Uncertainty of Passive Imager Cloud Optical Property Retrievals to Instrument Radiometry and Model Assumptions: Examples from MODIS

The optical and microphysical structure of clouds is of fundamental importance for understanding a variety of cloud radiation and precipitation processes. With the advent of MODIS on the NASA EOS Terra and Aqua platforms, simultaneous global-daily 1 km retrievals of cloud optical thickness (COT) and effective particle radius (CER) are provided, as well as the derived water path (WP). The cloud product (MOD06/MYD06 for MODIS Terra and Aqua, respectively) provides separate retrieval datasets for various two-channel retrievals, typically a VISNIR channel paired with a 1.6, 2.1, and 3.7 m spectral channel. The MOD06 forward model is derived from on a homogeneous plane-parallel cloud. In Collection 5 processing (completed in 2007 with a modified Collection 5.1 completed in 2010), pixel-level retrieval uncertainties were calculated for the following non-3-D error sources: radiometry, surface spectral albedo, and atmospheric corrections associated with model analysis uncertainties (water vapor only). The latter error source includes error correlation across the retrieval spectral channels. Estimates of uncertainty in 1 aggregated (Level-3) means were also provided assuming unity correlation between error sources for all pixels in a grid for a single day, and zero correlation of error sources from one day to the next. I n Collection 6 (expected to begin in late summer 2013) we expanded the uncertainty analysis to include: (a) scene-dependent calibration uncertainty that depends on new band and detector-specific Level 1B uncertainties, (b) new model error sources derived from the look-up tables which includes sensitivities associated with wind direction over the ocean and uncertainties in liquid water and ice effective variance, (c) thermal emission uncertainties in the 3.7 m band associated with cloud and surface temperatures that are needed to extract reflected solar radiation from the total radiance signal, (d) uncertainty in the solar spectral irradiance at 3.7 m, and (e) addition of stratospheric ozone uncertainty in visible atmospheric corrections. A summary of the approach and example Collection 6 results will be shown.

Cloud Properties Observations↗

Spectroscopic features of dissolved iodine in pristine and gamma-irradiated nitric acid solutions

While iodine speciation is important for a wide range of nuclear activities, understanding the mechanisms of the transformations of iodine between chemical forms and the sensitivity of these transitions to solution conditions and exposure to radiation remains an active area of research. This work curates spectroscopic data from several experimental techniques and establishes their sensitivity and limitations in detecting changes in iodine speciation in both neutral and acidic regimes. The techniques include Raman spectroscopy, Fourier Transform Infrared (FTIR) spectroscopy, 127 I NMR spectroscopy, and ultraviolet-visible (UV-Vis) spectroscopy. Analysis of these data indicates that these commonly accessible spectroscopies often have dynamic ranges of measurable concentrations that do not always overlap between all techniques. The experimental techniques are disparately sensitive to iodide (I - ), molecular iodine (I 2 ), iodate (IO 3 - ), and periodate (IO 4 - ) species. Raman, FTIR, and NMR spectra were subsequentially analyzed using two-dimensional correlation analyses to generate high-resolution autocorrelation spectra. Here, the use of these spectroscopies is then extended to tracking acidification-induced and gamma irradiation-induced transformations of dissolved sodium iodate in deionized water and concentrated nitric acid. Both dissolution into nitric acid and irradiation with a gamma source are demonstrated to perturb the iodine speciation promoting their assembly into molecular iodine (I 2 ) and/or triiodide (I 3 - ). While I 2 and I 3 - species are undetectable with FTIR spectroscopy and 127 I NMR spectroscopy, the species can be detected with UV-Vis spectroscopy, and in some instances, I 3 - can be detected with Raman spectroscopy in the low wavenumber region. Ultimately, the results of this work provide a path to designing optimal combinations of techniques to detect forms of iodine across a wide range of concentrations and conditions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

High-Speed Transport of Fluid Drops and Solid Particles via Surface Acoustic Waves

A compact sampling tool mechanism that can operate at various temperatures, and transport and sieve particle sizes of powdered cuttings and soil grains with no moving parts, has been created using traveling surface acoustic waves (SAWs) that are emitted by an inter-digital transducer (IDT). The generated waves are driven at about 10 MHz, and it causes powder to move towards the IDT at high speed with different speeds for different sizes of particles, which enables these particles to be sieved. This design is based on the use of SAWs and their propelling effect on powder particles and fluids along the path of the waves. Generally, SAWs are elastic waves propagating in a shallow layer of about one wavelength beneath the surface of a solid substrate. To generate SAWs, a piezoelectric plate is used that is made of LiNbO3 crystal cut along the x-axis with rotation of 127.8 along the y-axis. On this plate are printed pairs of fingerlike electrodes in the form of a grating that are activated by subjecting the gap between the electrodes to electric field. This configuration of a surface wave transmitter is called IDT. The IDT that was used consists of 20 pairs of fingers with 0.4-mm spacing, a total length of 12.5 mm. The surface wave is produced by the nature of piezoelectric material to contract or expand when subjected to an electric field. Driving the IDT to generate wave at high amplitudes provides an actuation mechanism where the surface particles move elliptically, pulling powder particles on the surface toward the wavesource and pushing liquids in the opposite direction. This behavior allows the innovation to separate large particles and fluids that are mixed. Fluids are removed at speed (7.5 to 15 cm/s), enabling this innovation of acting as a bladeless wiper for raindrops. For the windshield design, the electrodes could be made transparent so that they do not disturb the driver or pilot. Multiple IDTs can be synchronized to transport water or powder over larger distances. To demonstrate the transporting action, a video camera was used to record the movement. The speed of particles was measured from the video images.

Bar-Cohen, Yoseph↗

Plasma breakdown in bubbles passing between two pin electrodes

The ignition of plasmas in liquids has applications from medical instrumentation to manipulation of liquid chemistry. Formation of plasmas directly in a liquid often requires prohibitively large voltages to initiate breakdown. Producing plasma streamers in bubbles submerged in a liquid with higher permittivity can significantly lower the voltage needed to initiate a discharge by reducing the electric field required to produce breakdown. The proximity of the bubble to the electrodes and the shape of the bubbles play critical roles in the manner in which the plasma is produced in, and propagates through, the bubble. In this paper, we discuss results from a three-dimensional direct numerical simulation (DNS) used to investigate the shapes of bubbles formed by injection of air into water. Comparisons are made to results from a companion experiment. A two-dimensional plasma hydrodynamics model was then used to capture the plasma streamer propagation in the bubble using a static bubble geometry generated by the DNS The simulations showed two different modes for streamer formation depending on the bubble shape. In an elliptical bubble, a short electron avalanche triggered a surface ionization wave (SIWs) resulting in plasma propagating along the surface of the bubble. In a circular bubble, an electron avalanche first traveled through the middle of the bubble before two SIWs began to propagate from the point closest to the grounded electrode where a volumetric streamer intersected the surface. In an elliptical bubble approaching a powered electrode in a pin-to-pin configuration, we experimentally observed streamer behavior that qualitatively corresponds with computational results. Optical emission captured over the lifetime of the streamer curve along the path of deformed bubbles, suggesting propagation of the streamer along the liquid/gas boundary interface. Plasma generation supported by the local field enhancement of the deformed bubble surface boundaries is a mechanism that is likely responsible for initiating streamer formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Small-Scale Drop Size Variability: Impact on Estimation of Cloud Optical Properties

Most cloud radiation models and conventional data processing techniques assume that the mean number of drops of a given radius is proportional to volume. The analysis of microphysical data on liquid water drop sizes shows that, for sufficiently small volumes, this proportionality breaks down; the number of cloud drops of a given radius is instead proportional to the volume raised to a drop size-dependent nonunit power. The coefficient of proportionality, a generalized drop concentration, is a function of the drop size. For abundant small drops the power is unity as assumed in the conventional approach. However, for rarer large drops, it falls increasingly below unity. This empirical fact leads to drop clustering, with the larger drops exhibiting a greater degree of clustering. The generalized drop concentration shows the mean number of drops per cluster, while the power characterizes the occurrence frequency of clusters. With a fixed total number of drops in a cloud, a decrease in frequency of clusters is accompanied by a corresponding increase in the generalized concentration. This initiates a competing process missed in the conventional models: an increase in the number of drops per cluster enhances the impact of rarer large drops on cloud radiation while a decrease in the frequency suppresses it. Because of the nonlinear relationship between the number of clustered drops and the volume, these two opposite tendencies do not necessarily compensate each other. The data analysis suggests that clustered drops likely have a stronger radiative impact compared to their unclustered counterpart; ignoring it results in underestimation of the contribution from large drops to cloud horizontal optical path.

Knyazikhin, Y.↗

Small Launch Vehicle Design Approaches: Clustered Cores Compared with Multi-Stage Inline Concepts

In an effort to better define small launch vehicle design options two approaches were investigated from the small launch vehicle trade space. The primary focus was to evaluate a clustered common core design against a purpose built inline vehicle. Both designs focused on liquid oxygen (LOX) and rocket propellant grade kerosene (RP-1) stages with the terminal stage later evaluated as a LOX/methane (CH4) stage. A series of performance optimization runs were done in order to minimize gross liftoff weight (GLOW) including alternative thrust levels, delivery altitude for payload, vehicle length to diameter ratio, alternative engine feed systems, re-evaluation of mass growth allowances, passive versus active guidance systems, and rail and tower launch methods. Additionally manufacturability, cost, and operations also play a large role in the benefits and detriments for each design. Presented here is the Advanced Concepts Office's Earth to Orbit Launch Team methodology and high level discussion of the performance trades and trends of both small launch vehicle solutions along with design philosophies that shaped both concepts. Without putting forth a decree stating one approach is better than the other; this discussion is meant to educate the community at large and let the reader determine which architecture is truly the most economical; since each path has such a unique set of limitations and potential payoffs.

Waters, Eric D.↗

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗