Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Simulation”

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,459 records · Page 81

A numerical and experimental investigation of electrothermal aircraft deicing

Experimental data were obtained for the deicing characteristics of a stationary UH-1H helicopter blade which had been fitted with an electrothermal deicer assembly. The tests were run in the NASA Lewis Icing Research Tunnel, and yielded transient temperature responses for the substrate, heater and abrasion shield at selected positions around the blade. The data at the abrasion shield-ice interface clearly documented when melting, shedding or refreezing occurred. Comparisons were made between the experimental data and a one-dimensional numerical model. The agreement was generally very good, with the simulations being shown to be capable of predicting the transient temperature responses along with phase change and ice shedding. At many blade positions, the model was capable of accurately simulating the thermal response of the electrothermal deicer assembly.

Leffel, K.↗

JPS heater and sensor lightning qualification

Simulated lightning strike testing of the Redesigned Solid Rocket Motor (RSRM) field joint protection system heater assembly was performed at Thiokol Corp., Wendover Lightning Facility. Testing consisted of subjecting the lightning evaluation test article to simulated lightning strikes and evaluating the effects of heater cable transients on cables within the systems tunnel. The maximum short circuit current coupled onto a United Space Boosters, Inc. operational flight cable within the systems tunnel, induced by transients from all cables external to the systems tunnel, was 92 amperes. The maximum open-circuit voltage coupled was 316 volts. The maximum short circuit current coupled onto a United Space Boosters, Inc. operational flight cable within the systems tunnel, induced by heater power cable transients only, was 2.7 amperes; the maximum open-circuit voltage coupled was 39 volts. All heater power cable induced coupling was due to simulated lightning discharges only, no heater operating power was applied during the test. The results showed that, for a worst-case lightning discharge, the heater power cable is responsible for a 3.9 decibel increase in voltage coupling to operational flight cables within the systems tunnel. Testing also showed that current and voltage levels coupled onto cables within the systems tunnel are partially dependant on the relative locations of the cables within the systems tunnel.

Cook, M.↗

Global Positioning System Simulator Field Operational Procedures

Global Positioning System (GPS) simulation is an important activity in the development or qualification of GPS signal receivers for space flight. Because a GPS simulator is a critical resource it is highly desirable to develop a set of field operational procedures to supplement the basic procedures provided by most simulator vendors. Validated field procedures allow better utilization of the GPS simulator in the development of new test scenarios and simulation operations. These procedures expedite simulation scenario development while resulting in scenarios that are more representative of the true design, as well as enabling construction of more complex simulations than previously possible, for example, spacecraft maneuvers. One difficulty in the development of a simulation scenario is specifying various modes of test vehicle motion and associated maneuvers requiring that a user specify some (but not all) of a few closely related simulation parameters. Currently this can only be done by trial and error. A stand-alone procedure that implements the simulator maneuver motion equations and solves for the motion profile transient times, jerk and acceleration would be of considerable value. Another procedure would permit the specification of some configuration parameters that would determine the simulated GPS signal composition. The resulting signal navigation message, for example, would force the receiver under test to use only the intended C-code component of the simulated GPS signal. A representative class of GPS simulation-related field operational procedures is described in this paper. These procedures were developed and used in support of GPS integration and testing for many successful spacecraft missions such as SAC-A, EO-1, AMSAT, VCL, SeaStar, sounding rockets, and by using the industry standard Spirent Global Simulation Systems Incorporated (GSSI) STR series simulators.

Kizhner, Semion↗

Computer program for analysis of high speed, single row, angular contact, spherical roller bearing, SASHBEAN. Volume 1: User's guide

The computer program SASHBEAN (Sikorsky Aircraft Spherical Roller High Speed Bearing Analysis) analyzes and predicts the operating characteristics of a Single Row, Angular Contact, Spherical Roller Bearing (SRACSRB). The program runs on an IBM or IBM compatible personal computer, and for a given set of input data analyzes the bearing design for it's ring deflections (axial and radial), roller deflections, contact areas and stresses, induced axial thrust, rolling element and cage rotation speeds, lubrication parameters, fatigue lives, and amount of heat generated in the bearing. The dynamic loading of rollers due to centrifugal forces and gyroscopic moments, which becomes quite significant at high speeds, is fully considered in this analysis. For a known application and it's parameters, the program is also capable of performing steady-state and time-transient thermal analyses of the bearing system. The steady-state analysis capability allows the user to estimate the expected steady-state temperature map in and around the bearing under normal operating conditions. On the other hand, the transient analysis feature provides the user a means to simulate the 'lost lubricant' condition and predict a time-temperature history of various critical points in the system. The bearing's 'time-to-failure' estimate may also be made from this (transient) analysis by considering the bearing as failed when a certain temperature limit is reached in the bearing components. The program is fully interactive and allows the user to get started and access most of its features with a minimal of training. For the most part, the program is menu driven, and adequate help messages were provided to guide a new user through various menu options and data input screens. All input data, both for mechanical and thermal analyses, are read through graphical input screens, thereby eliminating any need of a separate text editor/word processor to edit/create data files. Provision is also available to select and view the contents of output files on the monitor screen if no paper printouts are required. A separate volume (Volume-2) of this documentation describes, in detail, the underlying mathematical formulations, assumptions, and solution algorithms of this program.

Aggarwal, Arun K.↗

Evaluation of Erosion of the Dummy “EE” Plate 19 in YA Type ATR Fuel Element During Reactor PALM Cycles

Advanced Test Reactor (ATR) Cycle 153B-1 was a 14-day, high-power, powered axial locator mechanism (PALM) operating cycle that completed on April 12, 2013. Cycle 153B-1 was a typical operating cycle for the ATR, and did not result in any unusual plant transients. ATR was started up and shut down as scheduled. The PALM drive physically moves the selected experiments into and out of the core to simulate reactor startup and heat up, and shutdown and cooldown transients, while the reactor remains in steady-state conditions. However, after the cycle was over, when the fuel elements were removed from the core and inspected, several thousand flow-assisted erosion pits and “horseshoeing” defects were readily observed on the surface of the several YA-type fuel elements (these are aluminum “dummy” plates that contain no fuel). In order to understand these erosion phenomena, a thermal-hydraulic model of coolant channel 20 on a YA-M fuel element was generated. The boundaries of the model were the aluminum EE plate of a YA-M fuel element and a beryllium reflector block with 13 horizontal saw cuts which represented regions of zero flow. The heat generated in fuel plates 1 through 18 was modeled to be passing through the aluminum EE plate. The coolant channel 20 width was set at 0.058 in. (58 mils). It was established that the horizontal saw cuts had a significant effect on the temperature of the coolant. The flow, which was expected to vary linearly with gradual heating of the coolant as it passed through the channel, was extremely turbulent. The temperature rise, which was expected to be a smooth “S” curve, was represented by a series temperature rise “humps,” which occurred at each horizontal saw cut in the beryllium reflector block. Each of the 13 saw cuts had a chamfered edge which resulted in the coolant flow being re-directed as a jet across the coolant channel into the surface of the EE plate, which explained the temperature rise and the observed scalloping and pitting degradation on the YA-M fuel elements. In the case of scalloping (horseshoeing) a surprising similarity of that defect to those appearing on aluminum plate rolled in over-lubrication conditions, were established. In turn, this made us think that the principal feature responsible for the appearance of these defects, was horizontal cuts in the beryllium reflector block created to arrest the propagation of large vertical crack(s) in Be in PALM cycles with higher overall fluence. This assumption was fully confirmed by the results of thermo-hydraulic simulations. The neutronics data for these modeling experiments were provided using advanced irradiation simulations (MCNP, HELIOS). In the case of pitting erosion the following corrective measures were proposed based upon the results of JMatPro v.8.2 modeling (TTT- and CCT-diagrams): change the fabrication process by adding blister anneal before program anneal, immediately after cold rolling of AA6061plate. This step will allow achieving complete recrystallization, eliminating of strengthening due to metastable precipitates, and reduce the possibility of forming sharp microstructural features upon the surface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region, India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.

58 GEOSCIENCES↗

Pressure transient analysis during CO 2 push-pull tests into faults for EGS characterization

Here, with the goal of detecting and characterizing faults and fractures in enhanced geothermal systems (EGS), a new technology involving CO 2 push-pull testing, active-source geophysical imaging, and well logging has recently been proposed. This technique takes advantage of (1) the contrasting properties of supercritical CO 2 and water which cause CO 2 to appear distinct from surrounding brine in seismic and other geophysical logging approaches, (2) the non-wetting nature of CO 2 which keeps it localized to the faults and fractures to create contrast potentially sufficient for active seismic and well-logging approaches to image faults and fracture zones at EGS sites. In this study, we evaluate the feasibility of using pressure transient monitoring during CO 2 push-pull tests to complement active seismic and wireline well logging for EGS characterization. For this purpose, we developed a 2D model of a prototypical geothermal site (Desert Peak, NV) that includes a single fault. The fault zone consists of a slip plane, fault gouge, and damage zone, and is bounded by the surrounding matrix of the country rock. Through numerical simulation using iTOUGH2, we found that the pressure transient at the monitoring wells in the fault gouge shows unique traits due to the multiphase flow conditions developed by CO 2 injection, and varies sensitively on the arrival of the CO 2 plume and the degree of CO 2 saturation. A sensitivity analysis shows the pressure transient is most sensitive to the fault gouge permeability, but also depends on multiphase flow parameters and the boundary conditions of the fault. Lastly, an inversion study reveals that the fault gouge permeability can be best estimated with the pressure transient data, whereas additional CO 2 saturation data do not improve the accuracy of the inversion significantly.

58 GEOSCIENCES↗

SIMULATIONS OF VISCOUS ACCRETION FLOW AROUND BLACK HOLES IN A TWO-DIMENSIONAL CYLINDRICAL GEOMETRY

We simulate shock-free and shocked viscous accretion flows onto a black hole in a two-dimensional cylindrical geometry, where initial conditions were chosen from analytical solutions. The simulation code used the Lagrangian total variation diminishing plus remap routine, which enabled us to attain high accuracy in capturing shocks and to handle the angular momentum distribution correctly. The inviscid shock-free accretion disk solution produced a thick disk structure, while the viscous shock-free solution attained a Bondi-like structure, but in either case, no jet activity nor any quasi-periodic oscillation (QPO)-like activity developed. The steady-state shocked solution in the inviscid as well as in the viscous regime matched theoretical predictions well. However, increasing viscosity renders the accretion shock unstable. Large-amplitude shock oscillation is accompanied by intermittent, transient inner multiple shocks. This oscillation of the inner part of the disk is interpreted as the source of QPO in hard X-rays observed in micro-quasars. Strong shock oscillation induces strong episodic jet emission. The jets also show the existence of shocks, which are produced as one shell hits the preceding one. The periodicities of the jets and shock oscillation are similar; the jets for the higher viscosity parameter appear to be stronger and faster.

79 ASTRONOMY AND ASTROPHYSICS↗

Open Source Modeling for Power Systems Dynamics

Time-domain simulations for power systems (both phasor and electromagnetic transients) are mostly dominated by commercial tools, such as PSS/E or PSCAD. However, with the increasing penetration of inverter-based resources and changes in model paradigm, new open-source tools are looking into become a possible alternative for a range of stakeholders. In this presentation we will focus on how to improve communication with industry and other stakeholders to properly showcase the value of open source simulators. Participants will identify current barriers facing the adoption of open-source time-domain tools, and will discuss how to improve the bench-marking processes across different tools, while maintaining software modularity to include new features and models as new challenges continue to happen in energy systems.

ENERGY PLANNING, POLICY, AND ECONOMY,MATHEMATICS A↗

Ring Model Development and Validation for Prismatic HTGR Core Thermal-hydraulics and Safety Analysis

Because of the complex core geometry, prismatic high temperature gas-cooled reactors (prismatic HTGRs) often exhibit complex thermal fluid behaviors during both normal operating and transient conditions. Most HTGR designs rely on passive safety system for decay heat removal, such as the reactor cavity cooling system (RCCS). During postulated accidents like Pressurized Conduction Cooldown (PCC) event, the decay heat is first radially transferred from the core region to the reactor vessel outer surface, then to the RCCS cooling panels. The peak fuel temperature is controlled by heat transfer mechanisms with two distinctive characteristic length scales, i.e., the core-wise effective heat conduction and local heat conduction in the fuel pellet scale. As both length scales are essential to determine the fuel temperature, from the modeling perspective, computer codes must be able to capture heat transfer in both scales. This is challenging for both computational fluid dynamics (CFD) tools because of extremely large amount of computation resources required, and for system analysis codes because of the challenge to model the complex core geometry. Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, efforts have been pursued to support HTGR technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior during safety-related transients. In our previous studies, the ring model was adopted in SAM and further developed to simulate the normal operating condition and a PCC event using the MHTGR-350 design of General Atomics as the reference design. This current work represents a continuation of these previous efforts, and the focus is to critically review and examine simplifications and assumptions made to develop the ring model, and to perform code validation using experimental data from an integral-effect test facility, the High Temperature Test Facility (HTTF) at the Oregon State University. In this study, the test PG-27 from the HTTF test suite was selected for code validation purpose. The test PG-27 is a transient test designed to simulate the PCC event of the MHTGR design. Very good agreements between SAM prediction and experimental measurements were found in both coolant and solid structure temperatures during the transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Uncertainty and Sensitivity Analysis of TREAT Transient Test no. 15

The U.S. Department of Energy (DOE) is preparing for the resumption of transient testing using the Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory (INL). In order to support the analysis of TREAT operation and experiment planning, the MAMMOTH reactor physics team has recently begun development of uncertainty and sensitivity analysis for the TREAT core simulated via MAMMOTH. MAMMOTH is a high-fidelity three-dimensional transient multi-physics code developed at INL under the DOE-NE Advanced Modeling and Simulation program (NEAMS). This research is part of the development presenting some initial uncertainty and sensitivity analysis of TREAT Transient Test no. 15. In this summary, we employ the Sobol sampling strategy implemented in RAVEN to efficiently perform the sensitivity and uncertainty analysis. The idea is using Sobol functional decomposition with generalized polynomial chaos (gPC) expansions to estimate the responses of interest. The gPC coefficients are then computed via the sparse grid spectral projection method. In addition, this strategy includes evaluation of the second or higher order sensitivity indices, allowing the study of interactions between several input variables. The principal focus of this analysis is to quantify the uncertainties in the total power and maximum fuel temperatures with respect to uncertainties in the initial transient conditions, such as rated power, fuel temperature, graphite fuel density and control rod boron concentration. In addition, both the first order and second order time-dependent sensitivity indices are computed and analyzed. The results show that the interactions between initial control rod boron concentration and the rated power will contribute at most 1% to the total uncertainty of the power peak.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Numerical simulations of microflare evolution in the solar transition region and corona

Several observers report transient ultraviolet brightenings, often referred to as microflares, in the solar atmosphere. In this paper, the results are presented of a series of one-dimensional numerical simulations examining possible relationships between microflares and the generation of dynamical chromospheric and transition region features. Low-energy and medium-energy microflares eject long-lived cool, dense gas plugs into the corona, with the gas plug traversing the loop apex in the medium energy case. In the case of high-energy microflares, the gas plug is rapidly heated to the temperature of the surrounding corona, and the results resemble the dynamics occurring in standard solar flare thick-target electron beam models.

Sterling, Alphonse C.↗

Heliosphere Responds to a Large Solar Wind Intensification: Decisive Observations from IBEX

Our heliosphere—the bubble in the local interstellar medium produced by the Sun's outflowing solar wind—has finally responded to a large increase in solar wind output and pressure in the second half of 2014. NASA's Interstellar Boundary Explorer ( IBEX ) mission remotely monitors the outer heliosphere by observing energetic neutral atoms (ENAs) returning from the heliosheath, the region between the termination shock and heliopause. IBEX observed a significant enhancement in higher energy ENAs starting in late 2016. While IBEX observations over the previous decade reflected a general reduction of ENA intensities, indicative of a deflating heliosphere, new observations show that the large (~50%), persistent increase in the solar wind dynamic pressure has modified the heliosheath, producing enhanced ENA emissions. The combination of these new observations with simulation results indicate that this pressure is re-expanding our heliosphere, with the termination shock and heliopause already driven outward in the locations closest to the Sun. The timing between the IBEX observations, a large transient pressure enhancement seen by Voyager 2 , and the simulations indicates that the pressure increase propagated through the heliosheath, reflected off the heliopause, and the enhanced density of the solar wind filled the heliosheath behind it before generating significantly enhanced ENA emissions. As a result, the coming years should see significant changes in anomalous cosmic rays, galactic cosmic radiation, and the filtration of interstellar neutral atoms into the inner heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiparticle Simulations of Quadrupolar Ion Detection in an Ion Cyclotron Resonance Cell with Four Narrow Aperture Detection Electrodes

The current paradigm in FT-ICR cell design is to approximate the ideal three-dimensional quadratic trapping potential as closely as possible to maintain ion cloud spatial coherence and achieve long transients, either with hyperbolically shaped electrodes, shimming electrodes, or by dynamic harmonization. In sharp contrast, the FT-ICR analyzer cell with four narrow aperture detection electrodes (NADEL) introduces significant anharmonic terms to the trapping potential. This analyzer cell is capable of quadrupolar detection by which one can measure a signal that is close to the unperturbed cyclotron frequency. This is far less sensitive to trapping potential and space charge shifts than the reduced cyclotron frequency measured in conventional ICR cells. The quadrupolar mode of ion detection in NADEL cells has been examined previously by SIMION simulations of ion clouds with up to 500 ions per simulation. Here, the behavior of the NADEL analyzer cell is examined through particle-in-cell (PIC) simulations, which allows us to examine the behavior of large populations (tens of thousands) of ions with space charge considerations, and to calculate the induced charge on the NADEL detection electrodes, and thus the transient signal. PIC simulations confirm a unique spatial distribution of the ions, with a coherent motion that results in long transient signals. Dependence of the ion cloud and image current signal on cell design, ion energy, and magnetron radius are examined. Coalescence effects are compared with those found in a dynamically harmonized cell. The insensitivity of the measured cyclotron frequency to space-charge is demonstrated both with simulations and experimentally. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Possible evidence of excitonic condensation in a topological insulator

The transient excitonic condensate is a nonequilibrium electron-hole Bardeen-Cooper-Schrieffer state in a photoexcited semiconductor and semimetal, where electron-hole pairs undergo a phase transition and condense into a single coherent quantum state. Despite numerous experimental works to realize the predicted excitonic condensation phase, experimental evidence still remains elusive. This is largely due to the absence of direct measurements of a material's transient momentum-dependent electronic structure and the excitonic state in the condensation regime. Here, using time and angle-resolved photoemission spectroscopy, we find direct evidence of a transient excitonic condensate in the spin-polarized spatially indirect excitonic topological states in Bi2Te3. Accompanying the formation of the excitonic topological states by photoexcitation, we reveal a splitting of the hole's and electron's quasi-equilibrium chemical potential followed by the band flattening and backbending of the transient topological surface state. Moreover, within the same momentum range, we report a reshaping of the bulk valence band in the form of a Mexican-hat-like Bogoliubov dispersion-hallmarks of the excitonic condensation, followed by the opening of an energy gap at the Fermi level. The fluence and temperature dependence of these renormalization effects are reminiscent of excitonic condensation within Bardeen-Cooper-Schrieffer (BCS)-like behavior. These results, together with theoretical simulation, point to the possible formation of a transient excitonic condensate and provide opportunities to manipulate topologically protected Bose condensates with light.

ARPES↗

Evaluation of Test/Analysis Correlation Methods for Crash Applications

A project has been initiated to improve crash test and analysis correlation. The work in this paper concentrated on the test and simulation results for a fuselage section. Two drop tests of the section were conducted. The first test was designed to excite the linear structural response for comparison with finite element modal analysis results. The second test was designed to provide data for correlation with crash simulations. An MSC.Dytran model was developed to generate nonlinear transient dynamic results. Following minor modifications, the same model was executed in MSC.Nastran to generate modal analysis results. The results presented in this paper concentrate on evaluation of correlation methodologies for crash test data and finite element simulation results.

Lyle, Karen H.↗

GhoT of the GhoT/GhoS toxin/antitoxin system damages lipid membranes by forming transient pores

Highlights: • Determined that toxin GhoT disrupts membranes by forming transient pores. • Determined that membrane charge affects the activity of GhoT. • Determined through molecular dynamic simulations why GhoT F38R is inactive. GhoT is a bacterial toxin of the type V toxin/antitoxin system that allows Escherichia coli to reduce its metabolism in response to oxidative and bile stress. GhoT functions by increasing membrane permeability and reducing both ATP levels and the proton motive force. However, how GhoT damages the inner membrane has not been elucidated. Here we investigated how GhoT damages membranes by studying its interaction with lipid bilayers and determined that GhoT does not cause macroscopic disruption of the lipid bilayer to increase membrane permeability to the dye carboxyfluorescein. Using circular dichroism, we found that GhoT forms an alpha helical structure in lipid bilayers that agrees with the structure predicted by the I-TASSER protein structure prediction program. The structure generated using I-TASSER was used to conduct coarse-grained molecular dynamics simulations, which indicate that GhoT damages the cell membrane, as a multimer, by forming transient transmembrane pores.

60 APPLIED LIFE SCIENCES↗

Transient fuel performance analysis for the preliminary fuel concept of general atomics fast modular reactor

This study investigates the transient fuel performance of General Atomics Fast Modular Reactor (GA-FMR) during accident scenarios, focusing on the behavior of its innovative fuel system that combines high-assay low enriched uranium dioxide (HALEUO2) fuel with SiGA® ceramic matrix composite silicon carbide cladding. The preliminary fuel design’s response was analyzed during reactivity-initiated accidents (RIA) and loss of coolant accidents (LOCA) using BISON fuel performance analysis code, which included both the diffusion enhanced and BISON-FASTGRASS coupled UO 2 models. The RIA analysis demonstrated that effective reactivity control reduced fuel temperature, though with transient fission gas release resulting in additional tensile stress state on the cladding. LOCA simulations revealed differing predictions between the two models: the BISON UO 2 model showed more transient fission gas release but minimal pellet expansion, while the BISON-FASTGRASS UO 2 model predicted less pronounced fission gas release but more fuel swelling and thermal expansion, potentially leading to pellet-cladding mechanical interaction. Here, these findings highlight critical areas for fuel design optimization and identify knowledge gaps requiring further experimental and computational investigation to advance GA-FMR fuel development.

Lee, Soon K. [Argonne National Laboratory (ANL), A↗