Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic processes”

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 577 records · Page 32

Central American mountains inhibit eastern North Pacific seasonal tropical cyclone activity

Abstract The eastern North Pacific (ENP) has the highest density of tropical cyclones (TCs) on earth, and yet the controls on TCs, from individual events to seasonal totals, remain poorly understood. One effect that has not been fully considered is the unique geography of the Central American mountains. Although observational studies suggest these mountains can readily fuel individual TCs through dynamical processes, here we show that these mountains indeed play the opposite role on the seasonal timescale, hindering seasonal ENP TC activity by up to 35%. We found that these mountains significantly interrupt the abundant moisture transport from the Caribbean Sea to the ENP, limiting deep convection over the open ocean area where TCs preferentially occur. This study advances our fundamental understanding of ENP TC genesis mechanisms across the weather-to-climate timescales, and also highlights the importance of topography representation in improving the ENP regional climate simulations, as well as TC seasonal predictions and future projections.

54 ENVIRONMENTAL SCIENCES↗

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗

Phenomenological models of two-particle correlation distributions on transverse momentum in relativistic heavy-ion collisions

Two-particle, pair-number correlation distributions on two-dimensional transverse momentum (p t1 , p t2 ) constructed from the particle production in relativistic heavy-ion collisions allow access to dynamical processes in these systems beyond what can be studied with angular correlations alone. Only a few measurements of this type have been reported in the literature, and phenomenological models, which facilitate physical interpretation of the correlation structures, are nonexistent. Ongoing efforts at the Relativistic Heavy-Ion Collider (RHIC) will provide a significant volume of these correlation measurements in the future. In anticipation of these new data, two phenomenological models are developed which describe two-dimensional (2D) correlation distributions on transverse momentum. One model is based on a collision event-by-event fluctuating blast wave. The other is based on event-by-event fluctuations in fragmenting color-flux tubes and in jets. Both models are shown to be capable of accurately describing the measured single-particle p t distributions for minimum-bias Au+Au collisions at √s NN = 200 GeV. Both models are then applied to preliminary, charged-particle correlation measurements on 2D transverse momentum. Here, the capabilities of the two models for describing the overall structure of these correlations, the stability of the fitting results with respect to collision centrality, and the resulting trends of the dynamical fluctuations are evaluated. In general, both phenomenological models are capable of qualitatively describing the major correlation structures on transverse momentum and can be used to establish the required magnitudes and centrality trends of the fluctuations. Both models will be useful for interpreting the forthcoming correlation data from the RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The MERiT Onboard the CeREs: A Novel Instrument to Study Energetic Particles in the Earth's Radiation Belts

The Miniaturized Electron pRoton Telescope, MERiT, is a low‐mass, low‐power, compact instrument using an innovative combination of particle detectors, sensor electronics, and onboard processing. MERiT is flying on the Compact Radiation belt Explorer, CeREs, a 3U CubeSat launched into a low earth orbit of 500‐km altitude and inclination of 85° on 16 December 2018. The primary and secondary science goals of CeREs are to investigate electron microbursts and to study solar particles. MERiT comprises a stack of solid state detectors (SSD) behind space facing avalanche photo diodes (APDs) surrounded by W‐Al shielding to reduce side‐penetrating particle background. The APD‐SSD combination enables measurement of electrons from 5 to 200 keV and 1 to 8 MeV; protons from 200–400 keV and 7–100 MeV in differential channels with energy resolution ΔE/E≈30% for both electrons and protons. MERiT measures microbursts with a high time resolution ranging from 4 to 16 ms and solar particles with a cadence of 1 s. MERiT energy channels and cadences are software configurable via algorithms and lookup tables residing on a field‐programmable gate array. The lookup tables can be changed via ground commands. MERiT geometry factor is 31 sq.cm‐sr and optimized to measure microbursts with the instrument viewing the local zenith in orbit. MERiT enables investigation of dynamical processes of radiation belt electron energization and loss, solar electron and proton transport, and their access to the Earth's polar caps. We describe the MERiT sensor design, calibration, operational modes, data products, and science goals.

Kanekal, S. G.↗

Effect of the abolition of intersubunit salt bridges on allosteric protein structural dynamics

A salt bridge, one of the representative structural factors established by non-covalent interactions, plays a crucial role in stabilizing the structure and regulating the protein function, but its role in dynamic processes has been elusive. Here, to scrutinize the structural and functional roles of the salt bridge in the process of performing the protein function, we investigated the effects of salt bridges on the allosteric structural transition of homodimeric hemoglobin (HbI) by applying time-resolved X-ray solution scattering (TRXSS) to the K30D mutant, in which the interfacial salt bridges of the wild type (WT) are abolished. The TRXSS data of K30D are consistent with the kinetic model that requires one monomer intermediate in addition to three structurally distinct dimer intermediates (I 1 , I 2 , and I 3 ) observed in WT and other mutants. The kinetic and structural analyses show that K30D has an accelerated biphasic transition from I 2 to I 3 by more than nine times compared to WT and lacks significant structural changes in the transition from R-like I 2 to T-like I 3 observed in WT, unveiling that the loss of the salt bridges interrupts the R–T allosteric transition of HbI. Besides, the correlation between the bimolecular CO recombination rates in K30D, WT, and other mutants reveals that the bimolecular CO recombination is abnormally decelerated in K30D, indicating that the salt bridges also affect the cooperative ligand binding in HbI. These comparisons of the structural dynamics and kinetics of K30D and WT show that the interfacial salt bridges not only assist the physical connection of two subunits but also play a critical role in the global structural signal transduction of one subunit to the other subunit via a series of well-organized structural transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Milky Way globular clusters in γ-rays: analysing the dynamical formation of millisecond pulsars

ABSTRACT Globular clusters (GCs) are evolved stellar systems containing entire populations of millisecond pulsars (MSPs), which are efficient γ-ray emitters. Observations of this emission can be used as a powerful tool to explore the dynamical processes leading to binary system formation in GCs. In this work, 9 yr of Fermi Large-Area Telescope data were used to investigate the γ-ray emission from all GCs in the Milky Way. Twenty-three clusters were found as γ-ray bright, with two of them never having been reported before. It was also found that magnetic braking probably has a smaller impact on the formation rate of binary systems in metal-rich GCs than previously suggested, while a large value for the two-body encounter rate seems to be a necessary condition. The influence of the encounter rate per formed binary was for the first time explored in conjunction with γ-ray data, giving evidence that if this quantity is very high, binary systems will get destroyed before having time to evolve into MSPs, thus decreasing the total number of MSPs in a GC. No extended emission was found even for clusters whose optical extent is ≈0.5°; all of them are point-like sources spatially in agreement with the optical cores of the GCs, supporting previous X-ray results of heavier objects sinking into the clusters’ cores via dynamical friction. The possibility of extrapolating these results to ultra-compact dwarf galaxies is discussed, as these systems are believed to be the intermediate case between GCs and dwarf galaxies.

de Menezes, Raniere↗

Grand challenges in the digitalisation of wind energy

The availability of large amounts of data is starting to impact how the wind energy community works. From turbine design to plant layout, construction, commissioning, and maintenance and operations, new processes and business models are springing up. This is the process of digitalisation, and it promises improved efficiency and greater insight, ultimately leading to increased energy capture and significant savings for wind plant operators, thus reducing the levelised cost of energy. Digitalisation is also impacting research, where it is both easing and speeding up collaboration, as well as making research results more accessible. This is the basis for innovations that can be taken up by end users. But digitalisation faces barriers. This paper uses a literature survey and the results from an expert elicitation to identify three common industry-wide barriers to the digitalisation of wind energy. Comparison with other networked industries and past and ongoing initiatives to foster digitalisation show that these barriers can only be overcome by wide-reaching strategic efforts, and so we see these as “grand challenges” in the digitalisation of wind energy. They are, first, creating FAIR data frameworks; secondly, connecting people and data to foster innovation; and finally, enabling collaboration and competition between organisations. The grand challenges in the digitalisation of wind energy thus include a mix of technical, cultural, and business aspects that will need collaboration between businesses, academia, and government to solve. Working to mitigate them is the beginning of a dynamic process that will position wind energy as an essential part of a global clean energy future.

17 WIND ENERGY↗

Microwave catalytic conversion of acetylene for co-production of hydrogen and carbon nanotubes

Natural gas conversion to hydrogen and solid carbon can drastically reduce the carbon footprint. Microwave plasma pyrolysis is an emerging process for chemical industries to directly convert methane (CH 4 ), the major component of natural gas, to hydrogen and carbon, which offers benefits such as fast process dynamics, flexibility, and high product yield. However, the formation of unwanted by-products, like acetylene (C 2 H 2 ), will require extra cost for gas separation. Plasma pyrolysis of CH 4 produces large amount of acetylene at downstream, which is currently used mainly for welding. Here, to avoid gas separation and upgrade acetylene from the downstream of methane plasma reactor, a novel approach toward its transformation to carbon nanotubes (CNT) and pure hydrogen (H 2 ) over Ni-based bimetallic catalyst driven by microwave irradiation has been reported in this work. The dominant gas product was hydrogen. Low concentrations of methane, ethane, and ethylene were observed in the product stream. Catalytic acetylene decomposition was carried out at 400 °C over Al 2 O 3 supported Ni catalyst. The results showed that, at 400 °C, acetylene was dehydrogenated to CNT and hydrogen was the dominant gas in product stream over 10Ni-1Pd-Al 2 O 3 catalyst. Characterizations of spent catalysts was conducted using Raman spectroscopy and transmission electron microscopy (TEM) to investigate the properties of carbon deposited over the catalyst during the catalytic acetylene decomposition. The results highlight that methane can be efficiently converted to hydrogen and CNT through 2-step process, microwave plasma and microwave catalytic conversion of intermediate acetylene, operated in a single reactor system driven by microwave.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Meridional Variations of C2H2 and C2H6 in Jupiter's Atmosphere from Cassini CIRS Infrared Spectra

The abundances of hydrocarbons such as acetylene (C2H2) and ethane (C2H6) in Jupiter's atmosphere are important physical quantities, constraining our models of the chemical and dynamical processes. However, our knowledge of these quantities and their vertical and latitudinal variations has remained sparse. The flyby of the Cassini spacecraft with Jupiter at the end of 2000 provided an excellent opportunity to observe the infrared spectrum with the Composite Infrared Spectrometer (CIRS) instrument, mapping the spatial variation of emissions from 10-1400 cm-1. CIRS spectra taken at the highest resolution (0.5 cm-1) in early December 2000 have been analysed to infer atmospheric temperatures in the stratosphere at 0.5-20 mbar via the v4 of CH4, and in the troposphere at 100-400 mbar, via the hydrogen collision-induced continuum absorption at 600-800 cm. Simultaneously, we have searched for meridional abundance variations in C2H2 and C2H6 via the v5 and vg bands respectively. Tropospheric absorption and stratospheric emission are highly anti-correlated at the CIM resolution, introducing a non-uniqueness into the retrievals, which means that vertical gradient and column abundance cannot be simultaneously found without additional constraints. If we assume the profile shapes from photochemical model calculations, we show that the column abundance of C2H2 must decrease sharply towards the poles, while C2H6 is constant or slightly increasing. The relevance of these results to current photochemical and dynamical knowledge of Jupiter's atmosphere is discussed.

Nixon, C. A.↗

Meridional Variations of C2H2 and C2H6 in Jupiter's Atmosphere from Cassini CIRS Infrared Spectra

Hydrocarbons such as acetylene (C2H2) and ethane (C2H6) are important tracers in Jupiter's atmosphere, constraining our models of the chemical and dynamical processes. However, our knowledge of the vertical and meridional variations of their abundances has remained sparse. During the flyby of the Cassini spacecraft in December 2000, the Composite Infrared Spectrometer (CIRS) instrument was used to map the spatial variation of emissions from 10-1400 cm(sup -1) (1000-7 microns). In this paper we analyze a zonally-averaged set of CIRS spectra taken at the highest (0.5 cm(sup -1)) resolution, to infer atmospheric temperatures in the stratosphere at 0.5-20 mbar via the v4 band of CH4, and in the troposphere at 150-400 mbar, via the H2 absorption at 600-800 cm(sup -1). Simultaneously, we retrieve the abundances of C2H2 and C2H6 via the v5 and vg bands respectively. Tropospheric absorption and stratospheric emission are highly anti-correlated at the CIRS resolution, introducing a non-uniqueness into the retrievals, such that vertical gradient and column abundance cannot both be found without additional constraints. Assuming profile gradients from photochemical calculations, we show that the column abundance of C2H2 decreases sharply towards the poles by a factor approximately 4, while C2H6 is unchanged in the north and increasing in the south, by a factor approximately 1.8. An explanation for the meridional trends is proposed in terms of a combination of photochemistry and dynamics. Poleward, the decreasing UV flux is predicted to decrease the abundances of C2H2 and C2H6 by factors 2.7 and 3.5 respectively at a latitude 70 deg. However, the lifetime of C2H6 in the stratosphere (5 x 10(exp 9)) is much longer than the dynamical timescale for meridional motions inferred from SL-9 debris (5 x 10(exp 8 s)), and therefore the constant or rising abundance towards high latitudes likely indicates that meridional mixing dominates over photochemical effects. For C2H2, the opposite occurs, with the relatively short photochemical lifetime (3 x 10(exp 7 s)), compared to meridional mixing times, ensuring that the expected photochemical trends are visible.

Nixon, C. A.↗

Processes Controlling Water Vapor in the Winter Arctic Tropopause Region

This work describes transport and thermodynamic processes that control water vapor near the tropopause during the SAGE III-Ozone Loss and Validation Experiment (SOLVE), held during the Arctic 1999/2000 winter season. Aircraft-based water vapor, carbon monoxide, and ozone measurements were analyzed so as to establish how deeply tropospheric air mixes into the Arctic lowermost stratosphere and what the implications are for cloud formation and water vapor removal in this region of the atmosphere. There are three major findings. First, troposphere-to-stratosphere exchange extends into the Arctic stratosphere to about 13 km. Penetration is to similar levels throughout the winter, however, because ozone increases with altitude most rapidly in the early spring, tropospheric air mixes with the highest values of ozone in that season. The effect of this upward mixing is to elevate water vapor mixing ratios significantly above their prevailing stratospheric values of above 5ppmv. Second, the potential for cloud formation in the stratosphere is highest during early spring, with about 20% of the parcels which have ozone values of 300-350 ppbv experiencing ice saturation in a given 10 day period. Third, during early spring, temperatures at the troposphere are cold enough so that 5-10% of parcels experience relative humidities above 100%, even if the water content is as low as 5 ppmv. The implication is that during this period, dynamical processes near the Arctic tropopause can dehydrate air and keep the Arctic tropopause region very dry during early spring.

Pfister, Leonhard↗

A comparative study of rapidly and slowly rotating dynamical regimes in a terrestrial general circulation model

As a preliminary step in the development of a general circulation model for general planetary use, a simplified vesion of thef GISS Model I GCM has been run at various rotation periods to investigate differences between the dynamical regimes of rapidly and slowly rotating planets. To isolate the dynamical processes, the hydrologic cycle is suppressed and the atmosphere is forced with perpetual annual mean solar heating. All other parameters except the rotation period remain fixed at their terrestrial values. Experiments were conducted for rotation periods of 2/3, 1, 2, 4, 8, 16, 64 and 256 days. The results are in qualitative agreement with similar experiments carried out previously with other GCMs and with certain aspects of one Venus GCM simulation. As rotation rate decreases, the energetics shifts from baroclinc to quasi-barotropic when the Rossby radius of deformation reaches planetary scale. The Hadley cell expands poleward and replaces eddies as the primary mode of large-scale heat transport. Associated with this is a poleward shift of the baroclinic zone and jet stream and a reduction of the equator-pole temperature contrast. Midlatitude jet strength peaks at 8 days period, as does the weak positive equatorial zonal wind which occurs at upper levels at all rotation periods. Eddy momentum transport switches from poleward to equatorward at the same period. Tropospheric mean static stability generally increases in the tropics and decreases in midlatitudes as rotation rate decreases, but the global mean static stability is independent of rotation rate. The peak in the eddy kinetic energy spectrum shifts toward lower wavenumbers, reaching wavenumber 1 at a period of 8 days. Implications of these results for the dynamics of Venus and Titan are discussed. Specifically, it is suggested that the extent of low-level convection determines whether the Gierasch mechanism contributes significantly to equatorial superrotation on these planets.

Del Genio, Anthony D.↗

Dynamics and X-ray emission of a galactic superwind interacting with disk and halo gas

There is a general agreement that the conspicuous extranuclear X-ray, optical-line, and radio-contiuum emission of starbursts is associated with powerful galactic superwinds blowing from their centers. However, despite the significant advances in observational studies of superwinds, there is no consensus on the nature of the emitting material and even on the emission mechanisms themselves. This is to a great extent a consequence of a poor understanding of dynamical processes in the starburst superwind regions. To address this issue, we have conducted two-dimensional hydrodynamical simulations of galactic superwinds. While previous similar studies have used a single (disk) component to represent the ISM of the starburst galaxy, we analyze the interaction of the wind with a two-component disk-halo ambient interstellar medium and argue that this two-component representation is crucial for adequate modeling of starbursts. The emphasis of this study is on the geometry and structure of the wind region and the X-ray emission arising in the wind material and the shocked gas in the disk and the halo of the galaxy. The simulation results have shown that a clear-cut bipolar wind can easily develop under a range of very different conditions. On the other hand, a complex 'filamentary' structure associated with the entrained dense disk material is found to arise within the hot bubble blown out by the wind. The flow pattern within the bubble is dominated equally by the central biconic outflow and a system of whirling motions r elated to the origin and development of the 'filaments'. The filament parameters make them a good candidate for optical-emission-line filamentary gas observed in starburst halos. We find that the history of mass and energy deposition in the starburst region of the galaxy is crucial for wind dynamics. A 'mild' early wind, which arises as a result of the cumulative effect of stellar winds from massive stars, produces a bipolar vertical cavity in the disk and halo gas without strongly affecting the gaseous disk, thus creating conditions for virtually free vertical escape of the hot gas at the later, much more violent supernova-dominated phases of the starburst. We calculate the luminosity, mass, and effective temperature of the X-ray emitting gas in the 'soft' (0.1 to 0.7 keV, 0.7 to 2.2 keV, and 0.1 to 2.2 keV) and 'hard' (1.6 to 8.3 keV) energy bands and estimate the contribution of different gaseous components to the X-ray flux in these bands. Analysis of these parameters enables us to make conclusions regarding the nature of the X-ray-emitting material. We have inferred that the bulk of the soft thermal X-ray emission from starbursts arises in the wind-shocked material of the disk and halo gas rather than in the wind material itself. This enables us to predict that the integrated soft X-ray spectra of starbursts need not show an overabundance of heavy elements which are believed to be produced copiously in the centers of starbursts. Unlike soft X-ray emission, the hard component of thermal X-ray emission is found to originate in the wind material ejected from the starburst region. However, the derived ratio of hard-to-soft X-ray luminosities is too small compared to that observed in starbursts. We conclude therefore that the observed hard X-ray emission of starbursts is probably not associated with the thermal emission of hot wind or ambient shocked gas. Typical temperatures of the bulk of the soft X-ray-emitting material in our very different models have been found to agree well with the ones estimated on the basis of the ROSAT data for the soft component of X-ray emission of nearby starbursts. We predict that temperatures of the extranuclear soft X-ray-emitting gas in starburst galaxies with heavy element abundances near solar should be close to T(sub Xs = 2 to 5 x 10(exp 6)K.

Suchkov, Anatoly A.↗

Identification of major matrix metalloproteinase-20 proteolytic processing products of murine amelogenin and tyrosine-rich amelogenin peptide using a nuclear magnetic resonance spectroscopy based method

Objective: The aim of this study was to identify major matrix metalloproteinase-20 (MMP20) proteolytic processing products of amelogenin over time and determine if the tyrosine-rich amelogenin peptide (TRAP) was a substrate of MMP20. Design: Recombinant 15N-labeled murine amelogenin and 13C,15N-labeled TRAP were incubated with MMP20 under conditions where amelogenin self-assembles into nanospheres. Digestion products were fractionated by reverse-phase high-performance liquid chromatography at various time points. Product identification took advantage of the intrinsic disorder property of amelogenin that results in little change to its fingerprint 1H-15N heteronuclear single-quantum coherence nuclear magnetic resonance spectrum in 2% acetic acid upon removing parts of the protein, allowing cleavage site identification by observing which amide cross peaks disappear. Results: The primary product in five out of the six major reverse-phase high-performance liquid chromatography bands generated after a 24 hour incubation of murine amelogenin with MMP20 were: S55-L163, P2-L147, P2-E162, P2-A167, and P2-R176. After 72 hours these products were replaced with five major reverse-phase high-performance liquid chromatography bands containing: L46-A170, P2-S152, P2-F151, P2-W45, and short N-terminal peptides. TRAP was completely digested by MMP20 into multiple small peptides with the initial primary site of cleavage between S16 and Y17. Conclusions: Identification of the major MMP20 proteolytic products of amelogenin confirm a dynamic process, with sites towards the C-terminus more rapidly attacked than sites near the N-terminus. This observation is consistent with nanosphere models where the C-terminus is exposed and the N-terminus more protected. One previously reported end-product of the MMP20 proteolytic processing of amelogenin, TRAP, is shown to be an in vitro substrate for MMP20.

amelogenesis imperfecta, protease, amelogenin, ena↗

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗

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↗

The Urban Observatory: A Multi-Modal Imaging Platform for the Study of Dynamics in Complex Urban Systems

We describe an “Urban Observatory” facility designed for the study of complex urban systems via persistent, synoptic, and granular imaging of dynamical processes in cities. An initial deployment of the facility has been demonstrated in New York City and consists of a suite of imaging systems—both broadband and hyperspectral—sensitive to wavelengths from the visible (∼400 nm) to the infrared (∼13 micron) operating at cadences of ∼0.01–30 Hz (characteristically ∼0.1 Hz). Much like an astronomical survey, the facility generates a large imaging catalog from which we have extracted observables (e.g., time-dependent brightnesses, spectra, temperatures, chemical species, etc.), collecting them in a parallel source catalog. We have demonstrated that, in addition to the urban science of cities as systems, these data are applicable to a myriad of domain-specific scientific inquiries related to urban functioning including energy consumption and end use, environmental impacts of cities, and patterns of life and public health. We show that an Urban Observatory facility of this type has the potential to improve both a city’s operations and the quality of life of its inhabitants.

54 ENVIRONMENTAL SCIENCES↗

DOE NETL Gas Hydrate Project FWP 65213: Kinetic Parameters for the Exchange of Hydrate Formers (Final Report)

This report documents the research activities sponsored by the U.S. Department of Energy, National Energy Technology Laboratory (NETL), by the Pacific Northwest National Laboratory (PNNL) for the gas hydrate project entitled, “Kinetic Parameters for the Exchange of Hydrate Formers.” This project started June 1, 2013 and concluded December 31, 2018, and comprised five budget periods, under FWP-65213. The overarching goal of this project was to gain an improved understanding of the dynamic processes of gas hydrate accumulations in geologic media by combining laboratory studies, numerical simulation, and analysis of shipboard infrared imaging of hydrate core samples. This project comprised four principal components: 1) fundamental laboratory investigations, 2) numerical simulator development and verification, 3) hydrate core characterization and analysis, and 4) applied laboratory and numerical investigations. Laboratory and field experiments have demonstrated the four core methodologies for producing natural gas hydrates from geologic formations: 1) depressurization, 2) thermal stimulation, 3) inhibitor injection, and 4) guest molecule exchange. Converting these methodologies into viable commercial technologies remains challenging, either from an economic or technical perspective. This project was largely focused with the technical challenges associated with the guest molecule exchange method, which is often referred to as gas swapping. Previous research by PNNL had shown that gas swapping with pure CO 2 was not feasible, as first originally thought, principally due to the formation of secondary hydrates, which clogged or reduced the effective permeability of the hydrate-bearing formation. This research was directed at discovering the potential for using mixtures of gases to affect the gas swapping methodology at commercial scales. The Ignik Sikumi Gas Hydrate Field Trial demonstrated the injection of a N 2 -CO 2 gas mixture into a natural gas hydrate bearing layer, and exchange behavior with the CH 4 hydrate former or guest molecules. The original concept behind the use of a N 2 -CO 2 was to maintain gas conditions at the reservoir state. Pure CO 2 would have existed under supercritical conditions, with an associated density closer to a liquid than gas. One of the key outcomes from this research project has been the discovery, through laboratory experiments and numerical simulation, that gas mixtures, especially those with components having different gas hydrate equilibria conditions provide a control option for the guest molecule exchange methodology. For example, the dissociation or formation of hydrate could be controlled by adjusting the composition of the gas mixture. This research project was dynamic in that the discoveries during each budget period would set the course of research during successive budget periods. Therefore, this report is arranged in chronological order progressing through the five budget periods of the project.

03 NATURAL GAS↗