Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local dynamic operator”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Thermal Hydraulics Safety Analyses of TRIGA Reactor Upon Fuel Reconfiguration

Safety analyses of the Training, Research, Isotope, General Atomics (TRIGA) reactor at Kansas State University(KSU) is performed using thermal hydraulic simulations based on changes to reactivity or the volumetric heat generation rate due to changes in the positioning of the enriched fuel elements in the reactor core. If the reactivity or neutronics parameters are known, then they can be used to simulate the thermal energy state of the material which is crucial to assess the safety of the nuclear fuel inside a nuclear reactor. Several thermal hydraulics models and codes have been developed in the past to assess reactor safety. RELAP (Reactor Excursion and Leak Analysis Program) is one such tool widely used throughout the world for safety analyses. Although RELAP is considered a reliable tool by regulators, vendors, and operators, it lacks the numerical fidelity provided by 3-D computational fluid dynamics (CFD) models. This work involves the critical task of assessing safety of a research reactor upon fuel reconfiguration based on models developed in RELAP and comparing them with a commercially available 3-D CFD program. The TRIGA Mark-II reactor, designed by General Atomics and located at Kansas State University, is a class of small pool type inherently safe research reactor licensed to operate at steady state power of up to 1.25 MW. It has a cylindrical core measuring 0.45 meters diametrically and 0.38 meters axially positioned in a tank measuring 1.98 meters diametrically and 6.25 meters axially. It has uranium zirconium hydride (UZrH) fuel rods housed in stainless steel cladding material arranged in five concentric rings (B,C,D,E and F rings, in that order, moving outwards from the center). TRIGA fuel is nominally 20% enriched, but the reactivity of the fuel can be changed by using burnable poisons (not used at KSU), or by varying the total weight percent of uranium in the UZrH fuel meat. The TRIGA reactor is designed to be an inherently safe reactor and is designed to be cooled by natural convection. Thus it is of paramount importance to avoid the film boiling regime; failing to do so will subject the fuel to severe conditions, possibly resulting in its failure due to fuel/clad interactions. A high DNBR ratio is a primary consideration in the reactor core design and thus any modification in the core should ensure that this criteria is always met. In order to fully utilize the capabilities of KSU TRIGA reactor facility, it should be operated at its full capacity, i.e. 1.25 MW, which requires increasing the fuel loading in some of the E- and F-ring elements from 8.5% U (nominal) to 12.0% U (nominal). Only four elements are to be switched to 12% loading, causing local power peaking in these elements. The degree of power peaking to be expected in these 12%-loaded fuel elements was calculated based on a detailed MCNP (Monte Carlo N-Particle code) model of the reactor core. Some preliminary results of the calculations are presented in this paper based on changes in peaking factors and power levels. To predict reactor thermal behavior, axial changes in the peaking factors obtained from MCNP models were used to calculate volumetric energy density functions for a single channel one dimensional RELAP model and a three dimensional COMSOL model of the TRIGA MARK-II reactor. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

What is the importance of climate model bias when projecting the impacts of climate change on land surface processes?

Regional climate change impact (CCI) studies have widely involved downscaling and bias-correcting (BC) Global Climate Model (GCM)-projected climate for driving land surface models. However, BC may cause uncertainties in projecting hydrologic and biogeochemical responses to future climate due to the impaired spatiotemporal covariance of climate variables and a breakdown of physical conservation principles. Here we quantify the impact of BC on simulated climate-driven changes in water variables(evapotranspiration, ET; runoff; snow water equivalent, SWE; and water demand for irrigation), crop yield, biogenic volatile organic compounds (BVOC), nitric oxide (NO) emissions, and dissolved inorganic nitrogen (DIN) export over the Pacific Northwest (PNW) Region. We also quantify the impacts on net primary production (NPP) over a small watershed in the region (HJ Andrews). Simulation results from the coupled ECHAM5/MPI-OM model with A1B emission scenario were firstly dynamically downscaled to 12 km resolutions with WRF model. Then a quantile mapping based statistical downscaling model was used to downscale them into 1/16th degree resolution daily climate data over historical and future periods. Two series climate data were generated according to the option of bias-correction (i.e. with bias-correction (BC) and without bias-correction, NBC). Impact models were then applied to estimate hydrologic and biogeochemical responses to both BC and NBC meteorological datasets. These im20 pact models include a macro-scale hydrologic model (VIC), a coupled cropping system model (VIC-CropSyst), an ecohydrologic model (RHESSys), a biogenic emissions model (MEGAN), and a nutrient export model (Global-NEWS). Results demonstrate that the BC and NBC climate data provide consistent estimates of the climate-driven changes in water fluxes (ET, runoff, and water demand), VOCs (isoprene and monoterpenes) and NO emissions, mean crop yield, and river DIN export over the PNW domain. However, significant differences rise from projected SWE, crop yield from dry lands, and HJ Andrews’s ET between BC and NBC data. Even though BC post-processing has no significant impacts on most of the studied variables when taking PNW as a whole, their effects have large spatial variations and some local areas are substantially influenced. In addition, there are months during which BC and NBC post-processing produces significant differences in projected changes, such as summer runoff. Factor-controlled simulations indicate that BC post-processing of precipitation and temperature both substantially contribute to these differences at region scales. We conclude that there are trade-offs between using BC climate data for offline CCI studies vs. direct modeled climate data. These trade-offs should be considered when designing integrated modeling frameworks for specific applications; e.g., BC may be more important when considering impacts on reservoir operations in mountainous watersheds than when investigating impacts on biogenic emissions and air quality (where VOCs are a primary indicator).

Liu, M. L.↗
DOE OSTI2019

DFT-FE – A massively parallel adaptive finite-element code for large-scale density functional theory calculations

We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching ~100,000 electrons) using Kohn-Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn-Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn-Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes---both in CPU-times and wall-times, and on both non-periodic and periodic systems---at systems sizes beyond a few thousand electrons, with over 5-10 fold speedups in systems with more than 10,000 electrons. In conclusion, the benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated on up to 192,000 MPI tasks.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2018

Modeling of Orbital and Attitude Dynamics of Nanosatellites Controlled via Active Electrostatic Charging

The large-scale exploration of airless bodies, such as asteroids and moons, is gaining interest, however it is limited by mobility issues: the lack of atmosphere, low gravity, and unknown soil properties pose difficult challenges for many forms of traditional locomotion. The environment in proximity of these bodies is also electrically charged due to interactions with solar wind and UV radiation. The EGlider (Electrostatic Glider) concept would be able to overcome these mobility issues by leveraging the natural environment, allowing operations in close proximity of the surface, while enabling long duration missions by minimizing propellant consumption. The EGlider is an advanced concept for small satellite mobility and propulsion, which relies on the electric fields naturally present around airless bodies in order to generate forces and torques useful for maneuvering. It does so by extending electrically charged appendages, which enable it to electrostatically soar above the surface. By differentially charging its electrodes it can also produce torques to control its attitude. The charges are maintained by continuous active ion or electron emission from the spacecraft, in order to cancel out the neutralizing influx of charges from surrounding plasma. An investigation of the spacecraft-plasma interaction was carried out. This included studying the effect of electrode geometry and calculating the charge-to-mass ratios required to enable several mission scenarios. Long, thin wire electrodes were identified to be the most power-efficient and would allow power-to-weight ratios achievable with current nanosatellite technologies. High electrode potential represents the main limiting factor for the system design. In order to test the feasibility of active control by means of differential charging, a simple 2D interaction model was developed, and a feedback controller to stabilize the vehicle was tested in a simulation environment. The results confirmed that good performance could be obtained for both position and attitude control. Finally, a dedicated software was developed for future simulation and testing of control strategies for the EGlider. This software allows to study the trajectory and attitude of an arbitrarily configured spacecraft in the proximity of an arbitrarily defined main airless body. The spacecraft can be assembled from basic parts, each with specified electrical, mass and optical properties. Efficient models allow to calculate gravitational and electrical interactions with the rotating main body and the local plasma field, as well as solar radiation pressure effects. Control models can be implemented as simple plug-in functions and easily tested. The preliminary validation campaign showed good matching with the reference cases that have been analyzed.

Corpino, Sabrina↗
DOE OSTI2015

Alcator C-Mod: research in support of ITER and steps beyond

This paper presents an overview of recent highlights from research on Alcator C-Mod. Significant progress has been made across all research areas over the last two years, with particular emphasis on divertor physics and power handling, plasma–material interaction studies, edge localized mode-suppressed pedestal dynamics, core transport and turbulence, and RF heating and current drive utilizing ion cyclotron and lower hybrid tools. Specific results of particular relevance to ITER include: inner wall SOL transport studies that have led, together with results from other experiments, to the change of the detailed shape of the inner wall in ITER; runaway electron studies showing that the critical electric field required for runaway generation is much higher than predicted from collisional theory; core tungsten impurity transport studies reveal that tungsten accumulation is naturally avoided in typical C-Mod conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

TRACER-Coastal Urban Boundary-Layer Interactions with Convection (TRACER-CUBIC) Field Campaign Report

To better understand the complicated web of processes governing convective cloud life cycle and aerosol-convection interactions, the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility supported deployment of a variety of advanced atmospheric measurement systems to the greater Houston, Texas, area from 1 October 2021 to 30 September 2022 as part of the Tracking Aerosol Convection Interactions Experiment (TRACER). Houston was selected as a study area because isolated convection and a variety of aerosol conditions are common in this region. This one-year ARM Mobile Facility (AMF) deployment featured a four-month intensive operational period (IOP) during summer 2022 (1 June–30 September). The ARM instrumentation was deployed at three sites along an east-west transect from La Porte, Texas to an ancillary site in a less-polluted rural region southwest of downtown Houston (Figure 1). At the La Porte Site, which is located near the Houston ship channel in an area that experiences significant pollution, the first ARM Mobile Facility (AMF1) was deployed. During the IOP, the ARM tethered balloon system (TBS) operated at the ancillary site. The second-generation C-Band Scanning ARM Precipitation Radar (CSAPR) operated near Pearland, Texas, roughly halfway between the Laporte and ancillary sites. As part of the TRACER- Coastal Urban Boundary-Layer Interactions with Convection (CUBIC) project, three boundary-layer profiling systems) were deployed along a north-south transect spanning from the University of Houston Coastal Center to the Aldine site north of downtown Houston (also blue dot in Figure 1) during the TRACER IOP. These systems included the National Oceanic and Atmospheric Administration (NOAA) National Severe Storms Laboratory CLAMPS2 (C2), which was deployed at the UHCC, the University of Wisconsin SPARC, which was deployed at the ARM CSAPR site near Pearland (orange diamond in middle of map in Figure 1), and the University of Oklahoma CLAMPS1 (C1), which was deployed at Aldine. These three systems have been successfully operated in various field campaigns, providing data sets that collectively offer new insights into atmospheric-boundary-layer (ABL) processes, sea-breeze (SB) circulations, and convection initiation (CI). For the TRACER IOP window, these systems ran continuously between 1 June and 26 September, 2022. Due to commitments to NOAA projects, the Doppler lidar at the UHCC site was not available until 24 June 2022. The CLAMPS and SPARC profiling systems are self-contained platforms that have benefited from several years of development and deployment. Instruments and data processing were maintained remotely, which made the 4-month deployment for the TRACER-CUBIC IOP period possible. The same basic instrument configuration comprises each system: a scanning Doppler wind lidar for flow characterization and passive profiler(s) for characterizing planetary-boundary-layer (PBL) thermodynamic properties. Each platform includes a Halo Streamline Doppler wind lidar, an Atmospheric Emitted Radiance Interferometer (AERI), and a surface meteorology station. CLAMPS1 and CLAMPS2 each also include a microwave radiometer (MWR, Figure 1d). The TRACER-CUBIC hypotheses included (i) Interactions of SB and urban circulations and how they affect the PBL structure in the Houston environment, causing spatially (horizontally and vertically) and temporally highly variable flow patterns, (ii) heat, moisture, and aerosol transport and mixing depend on these flow dynamics, and (iii) an improved understanding of the flow patterns and PBL structure are critical for investigating the processes leading to CI. To test these hypotheses, the project aimed at (i) characterizing SB circulations and their impacts on the diurnal evolution of the structure of the ABL, (ii) studying the evolution of Houston’s complex urban boundary layer, and (iii) identifying effects of urban-induced circulations on pre-convective environments. TRACER-CUBIC observations generally provide good coverage during the summer IOP. Initial screening of the data indicates a good number of cases with bay-breeze (BB) and/or SB signatures, local CI, and interesting boundary-layer features such as strong nocturnal low-level jets (LLJs, Table 1). The numbers listed in rows 3-5 in this table will be further updated as part of ongoing in-depth analyses and systematic identification of local circulations and CI events. More detailed information about the data availability and quality for each instrument is provided in the “readme” files that were submitted to the ARM Data Center along with each archived data sets. These “readme” files also provide instrument descriptions, information about the data collection and processing procedures, data formats, and any additional information relevant for further data analysis.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Compositional Impact of Io Volcanic Emissions on Jupiter's Magnetosphere and the Icy Galilean Moons

The magnetospheric ion population of Jupiter is dominated by the 1000 kg/s of iogenic material constantly ejected by IO volcanism as neutral gas (approx. 1 kg/s goes out as high speed dust grains), subsequent atmospheric losses to the IO torus, and radial transport of torus ions throughout the magnetosphere. As that magnetosphere is greatly distended in radial size by the iogenic plasma loading, so are surfaces of the other Galilean moons also significantly, and perhaps even dominantly, affected by iogenic plasma bombardment, e.g. at the level up to 0.2 kg/s heavy ions (mostly O and S) onto Europa as per local plasma ion measurements. In comparison, cometary impacts onto IO deliver about 0.02 kg/s of impact ejecta to Europa via ballistic transfer through the Jupiter system. The magnetosphere of this system operates as a powerful engine to produce and transport ions from the IO source to the surfaces of these other moons, and any future orbiter missions to these moons must account for surface distributions of the iogenic material and its chemical effects before real assessments can be made of sensible chemical materials otherwise arising from primordial formation and subsequent evolution of these moons. This is a fundamental problem of space weathering that must be addressed for all planetary bodies with thin atmospheres and direct surface exposure to their space plasma environments. Long-standing debates from Galileo Orbiter measurements about the origins of hydrate sulfates at Europa present examples of this problem, as to whether the sulfates arise from oceanic minerals or from iogenic sulfur chemistry. Any orbiter or landed mission to Europa for astrobiological investigations would further need to separate the potential chemical biosignatures of life or its precursors from the highly abundant background of iogenic material. Although no single ion carries a tag identifying it as of iogenic or other origin, the elemental abundance distributions of ions to be measured throughout the jovian magnetosphere and in the local moon environments can act as tracers if we know from direct measurements and models the distributions at the mostly likely sources, i.e. at IO. However, our knowledge of these abundances are very limited from earlier in-situ and remote measurements, mainly confined to major (S, O) and some minor (Na, K, Cl) species with abundances at or above a few percent relative to O. Future in-situ plasma measurements by the planned Jupiter Europa Orbiter and Jupiter Ganymede Orbiter missions should extend the abundance coverage to minor and even trace elemental species. For Europa astrobiological investigations it is also important to specify iogenic inputs and surface processing of isotopic species. We discuss the range of abundance distributions arising from models for IO hot volcanic emissions, and from the subsequent dynamics of ion injection, magnetospheric transport, and icy moon surface bombardment.

Cooper, John↗
DOE OSTI2026

Prototype Modeling for a Light-Trapping Planar-Cavity Enclosed Particle Solar Receiver

Concentrating solar thermal (CST) systems present a promising avenue for affordable and reliable energy production. Solar receivers are key components that determine the efficiency and longevity of these systems. Particle-based solar receivers have emerged as a compelling alternative to traditional technologies, offering several advantages that address limitations in current CST systems. This is especially true as next-generation CST technologies target applications including electricity generation, thermochemical processes, and industrial process heat, many of which necessitate higher operating temperatures than current commercial molten salt systems. Molten-salt thermal energy storage (TES) systems, commonly used in CSP, face challenges related to freezing and corrosion. Particle-based TES systems, in contrast, do not experience these issues, as particles are stable at high temperatures, exceeding 1000 degrees Celsius. This capability allows for a wider range of applications, including those requiring higher temperatures for industrial processes and efficient electricity generation. A novel innovation in particle-based solar receiver technology is the light-trapping planar cavity receiver (LTPCR) configuration developed by NREL. The LTPCR design consists of small cavity-like structures using opaque planar surfaces, enabling efficient capture and absorption of solar energy. A high incident flux concentration at the cavity aperture is absorbed on the receiver walls, and subsequently transferred to particles on the inside of cavities. The particles flow through the system, forming a fluidized bed inside of the receiver panels, effectively capturing the absorbed solar heat. Air is used as a fluidizing medium in this process to enhance particle heat transfer and mixing. The effectiveness of this design lies in its ability to manage solar flux conditions and ensure high solar-to-thermal receiver efficiency. A 100-kW prototype is currently being tested at the King Saud University in Saudi Arabia to assess the receiver performance. A range of modeling analyses for the optical, thermal, and mechanical effects were conducted to assess the performance of the receiver under on-sun conditions. The solar flux resulting from the KSU heliostat field was modeled using NREL SolTrace software and produced up to 600 kW/m2 at the receiver aperture. The solar flux absorbed on the receiver walls was then used within a computational fluid dynamics (CFD) model to predict wall temperature distributions along with radiation and convection loss. A two-phase CFD model was developed for the fluidized bed of silica sand inside the receiver panels to predict local wall-to-particle heat transfer coefficients, particle temperature distributions, and outlet temperature of the particles. We have also conducted analyses to understand the thermomechanical behavior of these innovative enclosed light-trapping solar receivers optimized for particle heating. We used finite element analysis (FEA) to predict the receiver's performance using temperature distributions obtained from CFD and based on the resulting stress profiles, evaluated creep-fatigue damage with a goal of achieving a 30-year service life. Analysis showed a significant impact of the particle-to-wall heat transfer coefficients (HTCs) on receiver performance, with higher HTCs resulting in reduced stress and increased lifespan. For instance, when using Inconel 740H, increasing the HTC from 800 W/m2 K to 1400 W/m2 K increased the creep life from 4,000 hours to over 100,000 hours. This highlights the importance of understanding and optimizing heat transfer in the design of high-efficiency receivers.

14 SOLAR ENERGY↗
DOE OSTI2025

Grid Resiliency with a 100% Renewable Microgrid

San Diego Gas & Electric Company (SDG&E) installed America’s first and largest utility-scale microgrid in Borrego Springs in 2013. The first generation Borrego Springs Microgrid utilized diesel generators to form and stabilize the microgrid island, with support from grid-scale batteries and local solar photovoltaic (PV) generation. In this project, SDG&E in partnership with National Renewable Energy Laboratory (NREL) demonstrated through modeling, simulation and utility field testing that blackstart and islanding of the microgrid can be led with 100% renewable, inverter based resources (IBRs), to help reduce community reliance on conventional generation resources. Through equipment upgrades, grid-forming island leader capability was transitioned to a battery IBR instead of the Borrego Springs Microgrid diesel generators. A new microgrid controller was integrated to the microgrid and programmed to control and manage multiple energy storage systems. Synchrophasor and other power quality data verified autonomous, high-speed response of the IBRs through blackstart, islanding, and load step testing. Results of project field evaluations provide distribution systems operators (DSO) with increased confidence that renewable, IBR can replace traditional generators to blackstart and island microgrids and rapidly establish stable island frequency with rapid changes in peak power demand. Importantly, the project validated the integration feasibility of a distributed energy resource management system (DERMS) controller that manages multiple grid-forming and grid-following IBRs, establishing a standard design interface to reduce the complexity of integrating new DERs in the future and supporting replication by the industry. As a result of learnings in this project, SDG&E has implemented the microgrid controller strategy at multiple other microgrid sites, thereby validating the replicability of the solution. Hardware-in-the-loop (HIL) simulations including power and controller HIL hardware — along with electromagnetic transient (EMT) simulations of Borrego Springs Microgrid —informed adjustments to inverter parameters and were important to characterize the performance of the IBRs in relevant operating conditions before deployment. The EMT and HIL simulations of islanding the entire community are important contributions in providing confidence in IBR performance prior to future islanding of the community in the field. High-fidelity EMT and/or HIL simulation of IBRs can de-risk field operations, and its relevance and importance as a tool is increasing as distribution grids and microgrids become more complex and dynamic with an increasing proportion of renewable generation, distributed energy storage, and two-way power and energy flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Grid-Connected Modular Soft-Switching Solid State Transformers (M-S4T)

The objective of this project is to develop and verify the concept of a flexible and modular soft-switching solid-state transformer (M-S4T) for direct grid-connected applications. The ability to directly connect power electronics converters to the medium voltage grid (4 kV – 13 kV), and to potentially replace the passive and bulky, but ubiquitous 60 hertz service transformer in the 25 kVA to 100 kVA range, with a more flexible and controllable device, has been regarded as the ‘holy grail’ in grid control. However, this has proven to be extremely difficult. This project has developed the solutions to several key challenges of the direct grid-connected power electronics and realized a 7.2 kV M-S4T prototype. First, a protection method to protect the M-S4T from the high voltages (110 kV for the 13 kV system) that occur on the grid due to transients and lightning strikes have been developed and experimentally verified. Second, the realization and the operation of the M-S4T based on high-voltage SiC devices (>3.3 kV) and a medium-frequency medium-voltage low-leakage transformer in a single-stage solid-state transformer with zero-voltage switching, low dv/dt, and low electromagnetic interference has been successfully demonstrated up to 7.5 kV peak. Third, an oil-cooling system and stable communication and distributed control system for converter module voltage sharing have been developed and experimentally verified. The developed M-S4T has realized a modular universal high-performance power conversion system. This conversion system is scalable to different voltage and power levels and adaptable to four-quadrant bidirectional operation. Moreover, the use of passive cooling techniques meets the equipment life requirements, and the lightning protection scheme fulfills the basic insulation level specifications for direct grid connection. Such power conversion system opens up near-term opportunities, including energy storage, solar PV, or electric vehicle charging with significant cost and footprint savings. In the longer term, the possibility of replacing the utility distribution transformer with an M-S4T will be transformative for future distribution grids with a compact footprint and full controllability to enable high renewable energy and storage penetration. In addition to the main project, this report expands on the Plus-Up projected including as part of the main award. This project developed and demonstrated the technology for autonomous collaborative inverters that can be connected in an ad hoc manner to the grid. The aim of the project was to: (1) evaluate the existing techniques for grid-connected inverters and find their limitations; (2) develop detailed requirements for grid-connected inverters in the modern grid with millions of active nodes; (3) design a unified control strategy that brings more autonomy and intelligence to grid-connected inverters, and addresses parts of the issues with the existing techniques. The proposed technique, called UniCon, enables inverters to 1) connect/disconnect to/from the grid in an ad hoc manner; (2) work based on local sensing. Slow communication could be used for a more optimized behavior; (3) work automatically in both grid-forming/grid-following mode; (4) handle large disturbances, e.g., big load step and fault, in an oscillation-free manner; (5) work collaboratively with other inverters in steady-state and during transients. UniCon can be implemented in the middle-level control; hence it is agnostic to the vendor and to the implementation of the inner voltage/current and protection loops. Furthermore, a new synchronization scheme, based on deep learning, was developed that can extract the grid voltage phase and amplitude in a stable manner. The method is cheap to implement can improve the dynamic performance of the grid-connected inverters during fast transients, e.g., fault. The proposed control scheme was validated by (1) MATLAB/Simulink; (2) hardware-in-the-loop results, and; (3) experimental results using three inverters that form a microgrid in a down-scaled feeder. Lastly, both the M-S4T and UniCon have achieved promising tangible paths to markets. In the case of the M-S4T, the underlying technology — the Soft Switching Solid State Transformer (S4T) developed at the Georgia Tech Center for Distributed Energy (GT-CDE) has been licensed by GridBlock from the Georgia Tech Research Corporation, and GridBlock has been working with manufacturing partner Jabil (one of the largest US-based contract manufacturers) and system integrator Power Secure (largest deployer of microgrids in the US with 4.7 GW under management), to meet the strong initial demand. Similarly, GridBlock has an exclusive license to the UniCon technology, developed under this award by GT-CDE. The UniCon provides an intermediate control layer that enables the implementation of the higher-level ‘transactive’ control commands for the system. The architecture of the system - slow communications with the cloud for system optimization and setpoints, and the use of locally measured quantities for real-time control, provide a very robust and secure way of implementing a real-time must-run grid that is also secure and stable. This is a brand-new functionality that is critical for the future grid and key to GridBlock’s business model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1996

Visualization of High Latitude Ion Upflow in Support of the Image Mission

The study of the magnetosphere is a 400 year old science that began with the publication by Gilbert, in 1600, of his hypotheses that the Earth was a giant magnet. Since then we have learned many things about the magnetosphere, particularly in the last 40 years of the space age, but we still have many unanswered questions. In spite of the many thousands of observations of this system we still lack a global understanding of how it works. This is due to its large size and tenuous nature that mean that any measurement made of the fields or particles involved only give one a knowledge of the local conditions at a given time. To gain a global perspective through such observations would require the simultaneous operation of thousands of satellites spread throughout the magnetospheric system in addition to observations made on the ground. Such a program would be impractical at least from financial considerations. What is needed for the advancement of magnetospheric physics is to develop the same capabilities that astrophysicists, solar physicists and meteorologists have been using for years --- the ability to stand back from the object under study and see it in its entirety. The challenge for doing this for the magnetosphere is that the particle densities are very low and the material is, for the most part, not luminous. In the last 25 years several ideas have been proposed that would allow at least the imaging of certain portions of the magnetosphere. These include imaging of the plasmasphere through the resonant scattering of solar 304 A from He+ ions, imaging of various hot plasma populations (i.e. the ring current, plasmasheet, upflowing ionospheric ions, etc.) from the neutral atoms that result when ions of these populations charge exchange with the hydrogen geocorona, and imaging the aurora at various wavelengths in the far ultraviolet. In addition, a novel technique for probing various boundaries in the magnetosphere by bouncing low frequency radio waves off of them has been extensively studied. Such a technique is analogous to the way the under water world can be probed with sonar. About five years ago NASA convened a science working group to study the possibility of flying a magnetospheric imaging mission. This resulted in a number of proposals for such a mission, one of which was selected to be the first MIDEX mission, to be launched in early 2000. The mission is called IMAGE (Imager for Magnetopause to Aurora Global Exploration) and its P.I. is J. Burch at SwRI. The IMAGE spacecraft will carry imagers to view the plasmasphere, aurora, ring current, inner plasmasheet, and upflowing ionospheric ions as well as a radio sounder to probe the location, shape and dynamics of the magnetopause, plasmapause, etc. Between its selection last April and the non advocacy mission review, which takes place next spring, the IMAGE teams needs to further refine the design of the mission and its instruments. The theory and modeling (T&M) subgroup of this team has the task of demonstrating what kind of images the instruments on IMAGE will see as well as showing that useful scientific information can be extracted from such images. As a central element to the efforts of the T&M subgroup we have decided to simulate and create synthetic images for the magnetic cloud event of October, 1995. In this event a large cloud, with high plasma densities and strong magnetic fields, ejected from the sun collided with the earth's magnetosphere triggering a three day period of intense magnetic storms and substorms. This event was observed from a number of different spacecraft and on the ground so we have a good data set to work with. In our work we will place the IMAGE spacecraft in the magnetosphere on its proposed orbit, with its proposed instruments, to see what it would see had it been there. Existing models of the plasmasphere, ring current and magnetopause will be run for this event to give the structures for the imaging instruments. There are several models which are lacking and which need to be developed. These include a model for the cusp, the inner plasmasheet and the upflowing ions. My task this summer was to develop the upflowing ion model and use it to create synthetic images.

Wilson, Gordon R.↗
DOE OSTI2015

Ultrafast dynamics of liquid water: Energy relaxation and transfer processes of the OH stretch and the HOH bend

The vibrational energy relaxation and transfer processes of the OH stretching and the HOH bending vibrations in liquid water are investigated via the theoretical calculation of the pump-probe spectra obtained from non-equilibrium molecular dynamics simulations with the TTM3-F interaction potential. The excitation of the OH stretch induces an instantaneous response of the high frequency librational motions in the 600-1000 cm-1 range. In addition, the excess energy of the OH stretch of a water molecule quickly transfers to the OH stretches of molecules in its first hydration shell with a time constant of ~50 fs, followed by relaxation to the HOH bends of the surrounding molecules with a time constant of 230 fs. The excitation of the HOH bend also results in the ultrafast excitation of the high frequency librational motions. The energy of the excited HOH bend of a water molecule decays, with a time constant of 200 fs, mainly to the relaxation of the HOH bends of its surrounding molecules. The energies of the HOH bends were found to transfer quickly to the intermolecular motions via the coupling with the high frequency librational motions. The excess energy of the OH stretch or the HOH bend relaxes to the high frequency intermolecular librational motions and eventually to the hot ground state with a time scale of ~1 ps via the coupling with the librational and translational motions. The energy relaxation and transfer processes were found to depend on the local hydrogen bonding network; the relaxations of the excess energy of the OH stretch and the HOH bend of four- and five-coordinated molecules are faster than those of a three-coordinated molecule due to the delocalization of the vibrational motions of the former (four- and five-coordinated molecules) compared to those of the later (three-coordinated molecules). The present results highlight the importance of the high frequency intermolecular librational modes in facilitating the ultrafast energy relaxation process in liquid water via their strong nonlinear couplings with the intramolecular OH stretching and HOH bending vibrations. S.S.X. acknowledges the support of the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. The calculation was carried out using the computing resources at the Research Center for Computational Science in Okazaki, Japan.

Imoto, Sho↗
DOE OSTI2015

Temporal Characterization of Hydrates System Dynamics beneath Seafloor Mounds. Integrating Time-Lapse Electrical Resistivity Methods and In Situ Observations of Multiple Oceanographic Parameters

This study was designed to investigate temporal variations in hydrate system dynamics by measuring changes in volumes of hydrate beneath hydrate-bearing mounds on the continental slope of the northern Gulf of Mexico, the landward extreme of hydrate occurrence in this region. Direct Current Resistivity (DCR) measurements were made contemporaneously with measurements of oceanographic parameters at Woolsey Mound, a carbonate-hydrate complex on the mid-continental slope, where formation and dissociation of hydrates are most vulnerable to variations in oceanographic parameters affected by climate change, and where changes in hydrate stability can readily translate to loss of seafloor stability, impacts to benthic ecosystems, and venting of greenhouse gases to the water-column, and eventually, the atmosphere. We focused our study on hydrate within seafloor mounds because the structurally-focused methane flux at these sites likely causes hydrate formation and dissociation processes to occur at higher rates than at sites where the methane flux is less concentrated and we wanted to maximize our chances of witnessing association/dissociation of hydrates. We selected a particularly well-studied hydrate-bearing seafloor mound near the landward extent of the hydrate stability zone, Woolsey Mound (MC118). This mid-slope site has been studied extensively and the project was able to leverage considerable resources from the team’s research experience at MC118. The site exhibits seafloor features associated with gas expulsion, hydrates have been documented at the seafloor, and changes in the outcropping hydrates have been documented, photographically, to have occurred over a period of months. We conducted observatory-based, in situ measurements to 1) characterize, geophysically, the sub-bottom distribution of hydrate and its temporal variability, and 2) contemporaneously record relevant environmental parameters (temperature, pressure, salinity, turbidity, bottom currents) to detect short-term changes within the hydrates system, identify relationships/impacts of local oceanographic parameters on the hydrates system, and improve our understanding of how seafloor instability is affected by hydrates-driven changes. A 2009 DCR survey of MC118 demonstrated that we could image resistivity anomalies to a depth of 75m below the seafloor in water depths of 1km. We reconfigured this system to operate autonomously on the seafloor in a pre-programmed mode, for periods of months. We designed and built a novel seafloor lander and deployment capability that would allow us to investigate the seafloor at potential deployment sites and deploy instruments only when conditions met our criteria. This lander held the DCR system, controlling computers, and battery power supply, as well as instruments to record oceanographic parameters. During the first of two cruises to the study site, we conducted resistivity surveying, selected a monitoring site, and deployed the instrumented lander and DCR, centered on what appeared to be the most active locations within the site, programmed to collect a DCR profile, weekly. After a 4.5-month residence on the seafloor, the team recovered all equipment. Unfortunately, several equipment failures occurred prior to recovery of the instrument packages. Prior to the failures, however, two resistivity profiles were collected together with oceanographic data. Results show, unequivocally, that significant changes can occur in both hydrate volume and distribution during time periods as brief as one week. Occurrences appear to be controlled by both deep and near-surface structure. Results have been integrated with seismic data from the area and show correspondence in space of hydrate and structures, including faults and gas chimneys.

58 GEOSCIENCES↗
NASA NTRS2014

Robonaut 2 - IVA Experiments On-Board ISS and Development Towards EVA Capability

Robonaut 2 (R2) has completed its fixed base activities on-board the ISS and is scheduled to receive its climbing legs in early 2014. In its continuing line of firsts, the R2 torso finished up its on-orbit activities on its stanchion with the manipulation of space blanket materials and performed multiple tasks under teleoperation control by IVA astronauts. The successful completion of these two IVA experiments is a key step in Robonaut's progression towards an EVA capability. Integration with the legs and climbing inside the ISS will provide another important part of the experience that R2 will need prior to performing tasks on the outside of ISS. In support of these on-orbit activities, R2 has been traversing across handrails in simulated zero-g environments and working with EVA tools and equipment on the ground to determine manipulation strategies for an EVA Robonaut. R2 made significant advances in robotic manipulation of deformable materials in space while working with its softgoods task panel. This panel features quarter turn latches that secure a space blanket to the task panel structure. The space blanket covers two cloth cubes that are attached with Velcro to the structure. R2 was able to open and close the latches, pull back the blanket, and remove the cube underneath. R2 simulated cleaning up an EVA worksite as well, by replacing the cube and reattaching the blanket. In order to interact with the softgoods panel, R2 has both autonomously and with a human in the loop identified and localized these deformable objects. Using stereo color cameras, R2 identified characteristic elements on the softgoods panel then extracted the location and orientation of the object in its field of view using stereo disparity and kinematic transforms. R2 used both vision processing and supervisory control to successfully accomplish this important task. Teleoperation is a key capability for Robonaut's effectiveness as an EVA system. To build proficiency, crewmembers have attempted increasingly difficult tasks using R2 inside the Station. After donning motion capture equipment and a virtual reality visor, Expedition 34/35 flight engineer Tom Marshburn began operations with simple hand movements. Having gained confidence, Marshburn guided R2's arms in a leader-follower exercise with crewmate Chris Cassidy. He was also able to use the hand to grab a tumbling roll of tape, a task only demonstrable in microgravity. Later efforts saw Cassidy handle softgoods through shared control with ground operators, mimicking an activity previously achieved using only autonomy. Robotic climbing through the ISS on handrails requires both precision motion and compliant grasps in order to both position grippers on handrails/seat track and prevent large internal forces. R2 climbs using actively controlled compliance and torque limiting to meet both the precision and softness requirements. During a step, the attached leg is controlled to be strong and stiff in order to maintain precision trajectory tracking. The swing leg is controlled to be stiff but weak to minimize unintentional impact forces while maintaining precision. During a simulated dual limb grasp (as shown in Figure 1), the R2 controller maintains one limb rigid and one limb soft to prevent large internal forces from building up. R2's grippers also use a form of force control to limit grip force while not fully closed on either a handrail or seat track thus limiting unintentional forces on cables/objects that may be present in R2's translational path. The on-board torso R2 safety system relies on a single end-effector velocity limit to prevent potential impact forces from exceeding Station maximum load requirements. R2's mobile configuration required modifications to the velocity limiting safety function due to its large, dynamic inertia. R2's legs maneuver the robot's mass creating configuration dependent, joint-relative inertias. A single all-encompassing velocity limit to cover worst case inertia is prohibitively low. The upgraded R2 control and safety systems solve this problem using momentum limiting, momentum control, and kinetic energy minimization. Momentum and kinetic energy take the robot mass into account relieving low velocity restrictions on low inertia end-effectors while ensuring that the overall mass of R2 is limited from hazardous velocities. The momentum of R2's five safety nodes (each of the four end-effectors and the body) is monitored and compared to a single momentum limit. If any of the five nodes exceeds the safety limit, the motor power is removed and the robot comes to a stop. Momentum control/limiting also provides a simple, reliable method to integrate hand held tools into the safety system by providing the tool mass to the control system thus automatically reducing the allowable velocity of the end-effector with the tool. Work on the ground continues to build the skill set for an EVA Robonaut. Recent experiments (Figure 2) demonstrate how a teleoperator can use R2 to manipulate a tether hook, an important safety precaution on spacewalks. Another task displayed Robonaut's ability to pull back a protective jacket over a hose and search for damage, as well as inspect a quick-disconnect fitting for debris. Demonstrations such as these are indicative of EVA work done on ISS, specifically seen during a series of spacewalks over 2012 and 2013 where astronauts searched for an ammonia leak in one of the external cooling loops. Through experiments both on ISS and on the ground, R2 is evolving and providing the information needed to plan out the upgrades that will make an EVA Robonaut an effective tool. With the addition of legs, R2 will start climbing inside the space station and supply invaluable information on how the climbing strategies and task stabilization techniques must be refined. Ground R2 systems will continue to work with additional EVA tools and equipment in preparation for onboard IVA testing and future EVA applications.

Diftler, Myron↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2026

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗
NASA NTRS2008

Analysis of a Dynamic Multi-Track Airway Concept for Air Traffic Management

The Dynamic Multi-track Airways (DMA) Concept for Air Traffic Management (ATM) proposes a network of high-altitude airways constructed of multiple, closely spaced, parallel tracks designed to increase en-route capacity in high-demand airspace corridors. Segregated from non-airway operations, these multi-track airways establish high-priority traffic flow corridors along optimal routes between major terminal areas throughout the National Airspace System (NAS). Air traffic controllers transition aircraft equipped for DMA operations to DMA entry points, the aircraft use autonomous control of airspeed to fly the continuous-airspace airway and achieve an economic benefit, and controllers then transition the aircraft from the DMA exit to the terminal area. Aircraft authority within the DMA includes responsibility for spacing and/or separation from other DMA aircraft. The DMA controller is responsible for coordinating the entry and exit of traffic to and from the DMA and for traffic flow management (TFM), including adjusting DMA routing on a daily basis to account for predicted weather and wind patterns and re-routing DMAs in real time to accommodate unpredicted weather changes. However, the DMA controller is not responsible for monitoring the DMA for traffic separation. This report defines the mature state concept, explores its feasibility and performance, and identifies potential benefits. The report also discusses (a) an analysis of a single DMA, which was modeled within the NAS to assess capacity and determine the impact of a single DMA on regional sector loads and conflict potential; (b) a demand analysis, which was conducted to determine likely city-pair candidates for a nationwide DMA network and to determine the expected demand fraction; (c) two track configurations, which were modeled and analyzed for their operational characteristic; (d) software-prototype airborne capabilities developed for DMA operations research; (e) a feasibility analysis of key attributes in the concept design; (f) a near-term, transitional application of the DMA concept as a proving ground for new airborne technologies; and (g) conclusions. The analysis indicates that the operational feasibility of a national DMA network faces significant challenges, especially for interactions between DMAs and between DMA and non-DMA traffic. Provided these issues are resolved, sectors near DMAs could experience significant local capacity benefits.

Wing, David J.↗
DOE OSTI2013

Comparisons of linear and nonlinear plasma response models for non-axisymmetric perturbations

With the installation of non-axisymmetric coil systems on major tokamaks for the purpose of studying the prospects of ELM-free operation, understanding the plasma response to the applied fields is a crucial issue. Application of different response models, using standard tools, to DIII-D discharges with applied non-axisymmetric fields from internal coils, is shown to yield qualitatively different results. The plasma response can be treated as an initial value problem, following the system dynamically from an initial unperturbed state, or from a nearby perturbed equilibrium approach, and using both linear and nonlinear models [A. D. Turnbull, Nucl. Fusion 52, 054016 (2012)]. Criteria are discussed under which each of the approaches can yield a valid response. In the DIII-D cases studied, these criteria show a breakdown in the linear theory despite the small 10{sup −3} relative magnitude of the applied magnetic field perturbations in this case. For nonlinear dynamical evolution simulations to reach a saturated nonlinear steady state, appropriate damping mechanisms need to be provided for each normal mode comprising the response. Other issues arise in the technical construction of perturbed flux surfaces from a displacement and from the presence of near nullspace normal modes. For the nearby equilibrium approach, in the absence of a full 3D equilibrium reconstruction with a controlled comparison, constraints relating the 2D system profiles to the final profiles in the 3D system also need to be imposed to assure accessibility. The magnetic helicity profile has been proposed as an appropriate input to a 3D equilibrium calculation and tests of this show the anticipated qualitative behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗