Search NASA⌕ Search

SEARCH · Search NASA

Results for “High velocity change”

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.

323 records · Page 18

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Analysis of beam ion driven Alfvén eigenmode stability induced by Tungsten contamination in EAST

Alfvén eigenmodes (AE) activity is observed in the EAST high β N and low B T discharge 93910, operation scenario dedicated to explore the ITER baseline scenario. AEs are triggered after the plasma is contaminated by Tungsten that causes an abrupt variation of the thermal plasma and energetic particles (EPs) profiles. The aim of the present study is to analyze the AE stability in the 93910 discharge using the gyro-fluid code FAR3d, identifying the AE stability trends by comparing the plasma before and after the Tungsten contamination. Tungsten contamination causes the destabilization of Toroidal AEs (TAE) and Energetic particle modes (EPMs) in the same frequency range and radial location with respect to the experimental observation and M3D-K/GTAW code results. Next, a set of parametric studies are performed to analyze the effect of the thermal plasma and EP parameters on the AE stability. The analysis indicates a lower EP β threshold for the AEs destabilization if the EP energy increases, an improved AE stability of on-axis NBI configurations due to the stronger continuum damping in the inner plasma region as well as a large enhancement of the EP drive as the thermal ion density increases due to a higher ratio of the EP and Alfven velocities. Consequently, the simulations indicate the increment of the thermal ion density after the Tungsten contamination could be the main cause of the AE/EPM destabilization.

AE↗

Laser Plasma Interaction in the MG Magnetic Field (Final Report)

Strong magnetic fields change properties of plasma. Employment of the magnetic fields can improve conditions for inertial confinement fusion. Compression and heating of the magnetized plasma is a base of the MagLIF pulsed power approach to the controlled fusion. Study of expansion and heating of plasma in the external magnetic fields, development of plasma instabilities, parametric and other effects are important for basic physics and applications. The main goal of the research for the grant was experimental investigation of laser produced plasma in the well characterized and controlled MG magnetic fields. We developed a platform for investigation of high intensity laser-plasma interaction (LPI) in the MG fields using a university-scale pulsed power machine. A Zebra machine at UNR generated transverse magnetic fields of 3 MG on the surface of the rod load and longitudinal fields of 1.4 MG in the coil loads. The magnetic fields were characterized with two-color Faraday rotation diagnostics. The pulsed power machine provided almost static magnetic fields for LPI compared to the much faster laser driven magnetic fields. We found that eddy currents in the rising magnetic field generate plasma on the metal targets. Eddy currents did not impact Si and CH dielectric targets. A Leopard laser with intensity in the focus up to 10 19 W/cm 2 produced plasma in the MG magnetic fields. The modified reflected Faraday diagnostics was used to study laser driven B-field at the MTW laser at Laboratory for Laser Energetics, UR. We demonstrated a slow dependence of the magnetic field in the coil on the applied laser intensity. The magnetic field increased by a factor of 2 if laser intensity increased by a factor of 30. This scaling is important for comparison of experiments at different laser facilities. Expansion of the plasma in the 2-3 MG external B-field was studied. Laser-produced plasma in the azimuthal magnetic field took the unique form of a thin disc expanding radially with a velocity of 250 km/s and confined in the vertical direction. A HYDRA MHD program at the University of Rochester was applied for simulations of LPI in the magnetic field B >3 MG. Simulations for the plasma disc parameters and the expansion dynamics were in a good agreement with experiments. Generation of narrow dense plasma jets in the longitudinal magnetic field of 0.6-0.8 MG was studied. Narrow plasma jets reached a length of 3-4 mm with the electron plasma density of (0.2-1.2)x10 20 cm -3 . A jet tip propagated with the velocity of 160-200 km/s. MHD simulations showed a good agreement of the dynamics of the formation of plasma jets with experiments. These jets are relevant to astrophysical jets. Simulations of K-shell spectra of Si plasma was performed for LPI in the B-field. The laser prepulse generated plasma near the laser target. A PrismSPECT model with a MeV electron beam produced by the laser was used. The spectral modelling showed the increased plasma density by a factor by 2-3 in the magnetic field. This was in agreement with the observed dynamics of the plasma plume confined by the axial magnetic field. The two-plasmon decay was studied in the 2.5-3 MG transversal magnetic fields. Strong 2-3.5 nm spectral widening and a 2-4 nm shifts of “red” and “blue” 3/2ω 0 spectral components were observed. The large shift and widening exceed the expected temperature and magnetic shift. PIC simulations are performing to clarify physics of this effect. The experimental research program for the grant is completed. New plasma effects in the MG fields were studied at the university-scale pulsed power generator. Two graduate and two undergraduate UNR and UR students carried experiments at the Zebra generator, supported and developed plasma diagnostics, provided data processing, and performed MHD simulations of LPI in the strong magnetic field. Research results are published in 9 referred papers and presented in 12 conferences. Program manager: Kramer U. Akli; Collaborators: R. C. Mancini and H. Sawada, University of Nevada, Reno (UNR); R. Betti and A. V. Maximov, University of Rochester (UR).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Near real-time streaming analysis of big fusion data

Experiments on fusion plasmas produce high-dimensional data time series with ever-increasing magnitude and velocity, but turn-around times for analysis of this data have not kept up. For example, many data analysis tasks are often performed in a manual, ad-hoc manner some time after an experiment. In this article, we introduce the Delta framework that facilitates near real-time streaming analysis of big and fast fusion data. By streaming measurement data from fusion experiments to a high-performance compute center, Delta allows computationally expensive data analysis tasks to be performed in between plasma pulses. This article describes the modular and expandable software architecture of Delta and presents performance benchmarks of individual components as well as of an example workflow. Focusing on a streaming analysis workflow where electron cyclotron emission imaging (ECEi) data is measured at KSTAR on the National Energy Research Scientific Computing Center's (NERSC's) supercomputer we routinely observe data transfer rates of about 4 Gigabit per second. In NERSC, a demanding turbulence analysis workflow effectively utilizes multiple nodes and graphical processing units and executes them in under 5 min. We further discuss how Delta uses modern database systems and container orchestration services to provide web-based real-time data visualization. For the case of ECEi data we demonstrate how data visualizations can be augmented with outputs from machine learning models. Here, by providing session leaders and physics operators, results of higher-order data analysis using live visualizations may make more informed decisions on how to configure the machine for the next shot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-Powered Autonomous Sensing System for Arctic Ocean using a Frequency-multiplied Cylindrical Triboelectric Nanogenerator

An autonomous sensing system for collecting environmental observations in the Arctic region is critical for the estimation and prediction of climate change. Ocean waves are a great source of energy for these sensing systems but there has been limited research done on small-scale energy harvesting applications in the Arctic Ocean (subsea or surface). The available wave energy in the Arctic Ocean is lower than the typical ocean wave energy due to the low wave frequency, height, and operating months. Although the available wave energy depends on the specific location in the Arctic Ocean, the wave height decreases everywhere during the winter (Jan-March) [1]. Our target location in this work is the Beaufort and Chukchi Seas, which are ice covered for about 195 days per year leaving only a few months (June-November) for wave energy harvesting [2]. During these few months, the average wave frequency is 0.2 Hz and the most common wave frequency is around 0.15 Hz. Many energy harvesting methods are unsuitable for use in the Arctic Ocean because of to the cold temperature, low wave frequency, and low wave height. Triboelectric nanogenerators (TENG) are one of the few energy harvesting methods that excel in these conditions. Zhong et al. [3] designed a stacked pendulum-structured TENG for low-frequency ocean wave energy and reported the device generated a peak power density of 11.2 W/m3 under the low wave frequency of 0.2 Hz. The autonomous sensing system we developed for the Arctic Ocean, the Arctic-TENG, is based on a frequency-multiplied cylindrical TENG (FMC-TENG) which is an optimized TENG configuration for Arctic Ocean conditions due to the high power density under low-frequency wave conditions [4]. Figure 1 shows an FMC-TENG with multiple pairs of free-standing triboelectric-layer mode materials (Aluminum and FEP). The mass and magnet attached to the rotor store gravitational potential energy which is released as kinetic energy when the potential energy overcomes the repulsive magnetic force generated by the opposing magnet attached to the stator. This force unbalance triggers a sudden rotation and swinging motion of the mass which increases the angular velocity of the system and therefore enhances the output power of the TENG system. The main components of the Arctic-TENG are a 3D-printed rotor and stator, electric and dielectric material adhesive tape, and bearings. All these materials were cold-soaked and tested at -40 °C in a chest freezer. Both the 3D-printed parts and adhesives were confirmed to have a minimal effect from the cold temperature. Multiple bearings were tested and the starting torque of each one was compared both at room temperature and -40 °C. The bearing with the lowest starting torque at -40 °C was selected for use in the Arctic-TENG. The Arctic-TENG was tested using an out-of-water motor-driven wave simulator (Figure 2). The wave simulator allows for controlled testing at a wave height of 0.2 m and frequencies between 0.1 Hz to 0.5 Hz. This wave simulator was used for both room temperature and -40 °C testing. The Arctic-TENG generated significantly more power at -40 °C compared to room temperature for each frequency tested. The system stored energy in a supercapacitor via a power management circuit, and the amount of energy stored per day was calculated at different wave frequencies. Based on the conditions at the proposed deployment location (days of non-ice-covered ocean and wave frequency), the amount of stored energy per year from Arctic-TENG was calculated to be enough energy for two transmissions every day. Durability testing was completed on the Arctic-TENG at room temperature to determine the lifetime. The rotor of the Arctic-TENG was attached to a DC motor and spun continuously for several millions of cycles without degradation of the electrical output, demonstrating the feasibility for long-term operation.

Jung, Hyunjun↗

INVESTIGATION OF TURBULENT JETS, SPRAYS, AND SUPERCRITICAL MIXING BY TIME RESOLVED RAINBOW SCHLIEREN DEFLECTOMETRY

As technological advancements have enabled higher efficiency compression ignition (CI) engines (diesel engines) by operating at higher pressures and temperatures, the mechanisms of fuel-air mixing are postulated to change. Thus, there is a greater need for experimental data to develop and verify new phenomena. However, the fast time scales and harsh environments associated with CI engines present many experimental difficulties and requires high quality non-intrusive diagnostics to obtain useful data. The diagnostics itself must be validated prior to its implementation in an engine environment. The purpose of this research is to advance rainbow schlieren deflectometry (RSD) to study turbulent fuel sprays including the effects of supercritical mixing in CI engine environments. This objective is met by developing and applying quantitative RSD, for the first time, to three separate aspects of fuel sprays: 1) turbulent mixing, 2) phase boundaries, and 3) non-ideal gas mixing. First, quantitative RSD is demonstrated for turbulent mixing in a canonical helium jet. Rayleigh scattering data in the literature is used to validate the concentration measurements obtained by RSD. The results also provide explanations for phenomena occurring in the helium jet and provide data not yet present in the literature. Second, RSD is applied to high pressure multi-phase fuel sprays and new methodologies are developed to distinguish the liquid region from the vapor region and to quantify the in-between region on a probabilistic basis. Isolating liquid and vapor regions generally requires two separate diagnostics, e.g., Mie scattering (liquid region) and Rayleigh scattering (vapor region). However, the results of this study demonstrate that RSD can be used as a single diagnostic to measure both liquid and vapor boundaries. Third, the optical-to-thermodynamic relations associated with refrac-tometry are investigated under non-ideal gas mixing conditions to develop a generalized relationship valid for both ideal or non-ideal gas mixing. The last part of this research develops the theoretical tools to understand and experimentally realize supercritical mixing at diesel conditions. Guided by the theoretical analysis, RSD is applied to compare supercritical versus transcritical mixing at CI conditions using a constant pressure test rig. It is found that while fully supercritical mixing ers benefits in diesel engines, i.e, faster mixing times (50% higher velocities), the absence of droplets which cause high levels of particulate emissions, etc., it would be difficult to achieve in practice because of the requirements of fuels with low critical temperature and/or substantial fuel preheating.

02 PETROLEUM↗

Influence of Vertical Hydrologic Exchange Flow, Channel Flow, and Biogeochemical Kinetics on CH 4 Emissions From Rivers

Abstract CH 4 emissions from inland water are highly uncertain in the current global CH 4 budget, especially for rivers and streams due to sparse measurements and the uncertainty of measurements caused by turbulent water flow. A previous study has revealed that vertical hydrologic exchange flow (VHEF) is the main regulator of CH 4 emissions from riverbed sediments. However, to what extent the understanding obtained from the plot‐scale can be extended to the reach scale and basin scale remains unknown. To address this challenge, we developed a process‐based model to estimate CH 4 flux at the air‐water interface using the attributes available in the national hydrography data set. It calculates the annual mean flux of VHEF, CH 4 production in sediments, and CH 4 transport in the river channel in a sequential manner. Model performance is evaluated by CH 4 efflux observed at the Hanford reach of the Columbia River. We show that reach‐wise sediment hydrologic and biogeochemical conditions estimated from the national hydrography data set could serve as a good indicator of CH 4 emissions from rivers. Aerobic methane oxidation and export to the downstream are the dominant ways of total CH 4 loss for the large lowland river. The hotspots of CH 4 emissions are likely to be at the reaches with fine sediments and slow channel velocity. This study demonstrates the possibility of quantifying CH 4 emissions at the reach scale and the modeling framework has the potential to be extended to the basin scale to improve estimates of CH 4 emissions from lotic inland water.

58 GEOSCIENCES↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗

Fast ion studies in the extended high-performance high β P plasma on EAST

Comprehending and optimizing fast ion behaviors is critical for the enhancement of performance in Experimental Advanced Superconducting Tokamak (EAST). This study explores the potential benefits of several factors that can improve the fast ion confinement. First, experiments show the change in the direction of the NBI2 from counter-I p to co-I p leads to a significant reduction in fast ion losses. TRANSP/NUBEAM simulation and tomography results based on fast-ion D-alpha measurements reveal that after the neutral beam injection (NBI) upgrade, the beam ion prompt loss is reduced by approximately 50%. Second, the upgraded ion cyclotron resonant frequency (ICRF) antenna at the N-port features twice the coupling resistance of the original antennas at EAST. This improved ICRF power coupling has enhanced the synergistic heating effect of NBI + ICRF, where the ICRF wave field accelerates beam ions at the harmonics. Experiments demonstrate that NBI + ICRF synergistic not only enhances plasma neutron yield and β P , but also accelerates beam ions to hundreds of keV. Further, the electron density and the neutral beam voltage have been optimized to reduce the fast ion slowing-down time and beam ion losses. Experimental and simulation results indicate that increasing the electron density reduces beam ion losses and enhances the bootstrap current fraction. While higher beam voltage results in a slight decrease in beam power absorption, it can increase the fraction of bootstrap current. With the understanding of these optimization of fast ion confinement, experiments have demonstrated fully non-inductive operation at high density (n e /n G ∼ 0.67, β P ∼ 3.1, β N ∼ 2.1, H 98,y2 ∼ 1.2) even without the support of co-I p beam NBI2. This investigation presents a potential regime to enhance fast ion confinement and extend performance in the high β P plasma for future experiments.

EAST tokamak↗

Constraining the Twomey effect from satellite observations: issues and perspectives

Abstract. The Twomey effect describes the radiative forcing associated with a change in cloud albedo due to an increase in anthropogenic aerosol emissions. It is driven by the perturbation in cloud droplet number concentration (ΔNd, ant) in liquid-water clouds and is currently understood to exert a cooling effect on climate. The Twomey effect is the key driver in the effective radiative forcing due to aerosol–cloud interactions, but rapid adjustments also contribute. These adjustments are essentially the responses of cloud fraction and liquid water path to ΔNd, ant and thus scale approximately with it. While the fundamental physics of the influence of added aerosol particles on the droplet concentration (Nd) is well described by established theory at the particle scale (micrometres), how this relationship is expressed at the large-scale (hundreds of kilometres) perturbation, ΔNd, ant, remains uncertain. The discrepancy between process understanding at particle scale and insufficient quantification at the climate-relevant large scale is caused by co-variability of aerosol particles and updraught velocity and by droplet sink processes. These operate at scales on the order of tens of metres at which only localised observations are available and at which no approach yet exists to quantify the anthropogenic perturbation. Different atmospheric models suggest diverse magnitudes of the Twomey effect even when applying the same anthropogenic aerosol emission perturbation. Thus, observational data are needed to quantify and constrain the Twomey effect. At the global scale, this means satellite data. There are four key uncertainties in determining ΔNd, ant, namely the quantification of (i) the cloud-active aerosol – the cloud condensation nuclei (CCN) concentrations at or above cloud base, (ii) Nd, (iii) the statistical approach for inferring the sensitivity of Nd to aerosol particles from the satellite data and (iv) uncertainty in the anthropogenic perturbation to CCN concentrations, which is not easily accessible from observational data. This review discusses deficiencies of current approaches for the different aspects of the problem and proposes several ways forward: in terms of CCN, retrievals of optical quantities such as aerosol optical depth suffer from a lack of vertical resolution, size and hygroscopicity information, non-direct relation to the concentration of aerosols, difficulty to quantify it within or below clouds, and the problem of insufficient sensitivity at low concentrations, in addition to retrieval errors. A future path forward can include utilising co-located polarimeter and lidar instruments, ideally including high-spectral-resolution lidar capability at two wavelengths to maximise vertically resolved size distribution information content. In terms of Nd, a key problem is the lack of operational retrievals of this quantity and the inaccuracy of the retrieval especially in broken-cloud regimes. As for the Nd-to-CCN sensitivity, key issues are the updraught distributions and the role of Nd sink processes, for which empirical assessments for specific cloud regimes are currently the best solutions. These considerations point to the conclusion that past studies using existing approaches have likely underestimated the true sensitivity and, thus, the radiative forcing due to the Twomey effect.

54 ENVIRONMENTAL SCIENCES↗

Uranium Oxide Elutriation at Higher Flow Rates and Longer Times

This work supports the recovery of high-assay low-enriched uranium (i.e., U enriched in 235 U to >5% and <20%) from used nuclear fuel. One step includes elutriating triuranium octoxide (U 3 O 8 ) from a fluidized bed of alumina (Al 2 O 3 ) that fills the main reaction chamber. The experiments outlined in this report assessed the effects of the time and flow rate on the elutriation of U 3 O 8 powder from the reaction vessel. Testing at Oak Ridge National Laboratory (ORNL) was performed in a one-fourth scale experimental apparatus based on a pilot-scale system, which is the Material Recovery Pilot Plant (MRPP) at Idaho National Laboratory (INL). The experimental system comprises a supply of pressurized nitrogen (N 2 ), a gas heater, the reaction vessel, a knock-out pot, and a water trap where elutriated particles are collected for sampling. An Al 2 O 3 fluidized bed is contained within the reaction vessel. A bed of Al 2 O 3 with a nominal size of 180 µm and U 3 O 8 with a nominal size of 7 µm was fluidized at 2 different flow rates so the U 3 O 8 physically separated from the bed. The apparatus was simple and had no valves or sharp bends downstream from the reactor. The flow exited the reactor through a side port at the top. The flow passed through a knock-out vessel before the flow was passed through a bubbler to capture the U 3 O 8 . Seven tests were completed with U 3 O 8 using an elutriating gas of N 2 heated to 200°C. Tests were performed with bed loadings that comprise 25.0 g of U 3 O 8 , or 3.6 wt% of the bed mass. Except for one test performed at 12 SLPM, the gas flow rate through the bed for the other experiments was set at 10 SLPM as measured by the mass flow meter. Linear velocities through the bed were calculated as 28 and 34 cm/s, or 0.9 and 1.1 ft·s -1 . The duration of the tests ranged from 20 to 180 min, and one test had a flow cycled between 0 and 5 SLPM every 5 min for 40 min. The testing successfully demonstrated elutriation of U 3 O 8 from the fluidized bed as up to a 100% of the U 3 O 8 was transferred out of the bed. The removal of U 3 O 8 had a logarithmic dependence on flow rate. The design changes made the apparatus easier to clean, and the mass balances were better than those used in previous studies. Some elutriated Al 2 O 3 collected in the tubing that connected the reactor to the knock-out pot, but the amount of Al 2 O 3 was minimal for 10 SLPM. The flow at 12 SLPM showed considerable transport of Al 2 O 3 bed material into the reaction-vessel headspace and the transfer line to the knock-out pot.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Improved high-gradient performance for medium-velocity superconducting half-wave resonators: Surface preparation and trapped flux mitigation

A development effort to improve the performance of superconducting radio-frequency half-wave resonators (SRF HWRs) is underway at the Facility for Rare Isotope Beams (FRIB), where 220 such resonators are in operation. Our goal was to achieve an intrinsic quality factor (𝑄 0 ) of ≥ 2 × 10 10 at an accelerating gradient (𝐸 acc ) of 12 MV/m. FRIB production resonators were prepared with buffered chemical polishing. First trials of electropolishing (EP) and post-EP low-temperature baking of FRIB HWRs allowed us to reach higher gradient (15 MV/m, limited by quench) with a higher quality factor at high gradient, but 𝑄 0 was still below our goal. Trapped magnetic flux during the Dewar test was found to be a source of 𝑄 0 reduction. Three strategies were used to reduce the trapped flux: (i) adding a local magnetic shield (LMGS) to supplement the “global” magnetic shield around the Dewar for reduction of the ambient magnetic field; (ii) performing a “uniform cooldown” (UC) to reduce the thermoelectric currents; and (iii) using a compensation coil to further reduce the ambient field with active field cancellation (AFC). The LMGS improved the 𝑄 0 , but not enough to reach our goal. With UC and AFC, we exceeded our goal, reaching 𝑄 0 = 2.8 ×10 10 at 𝐸 acc = 12 MV/m.

Cryogenics & vacuum technology↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

AmeriFlux US-GL1 Stannard Rock

This is the AmeriFlux version of the carbon flux data for the site US-GL1 Stannard Rock. Site Description - Stannard Rock is located 39 km from the nearest shore (Keweenaw Peninsula) in Lake Superior, 44 miles NNE of Marquette, Michigan, and 24 miles ESE of Manitou Island. The site is located on the historic Stannard Rock Lighthouse, which was completed in 1882. Eddy covariance instrumentation was installed in 2008 by a network of scientists from both US and Canada, eventually to be called the Great Lakes Evaporation Network (GLEN). The intent of GLEN has been to provide observations of over-lake meteorology and evaporation, improve forecasting of Great Lakes water levels, and support a wide variety of stakeholders, including the National Weather Service (NWS), Environment and Climate Change Canada, National Oceanic and Atmospheric Administration, U.S. Coast Guard, recreational boaters and commercial shipping, emergency management officials, and the Great Lakes research community. The eddy covariance station along with other ancillary meteorological instrumentation is located at an approximate elevation of 39.2 meters above mean lake water level. Meteorological data from the lighthouse are sent to the National Data Buoy Center, where they can be viewed in real-time at http://www.ndbc.noaa.gov/station_page.php?station=stdm4 Uncorrected half-hour fluxes were computed directly on the logger using a 30-minute block averaging period as high-frequency data were not available. The half hour flux measurements downloaded from the datalogger were post-processed with the following filters and corrections. Latent and sensible heat and carbon dioxide fluxes were corrected with 2-D coordinate rotation. Although the primary objective of data collection at Stannard Rock was to quantify the evaporative flux, additional preliminary measurements of the carbon dioxide concentration and flux from the LI-7500 are also included in this dataset but it is advised to use the carbon data with caution. Carbon data reported in this dataset includes turbulent fluxes of CO2 with no storage correction (FC, µmol m-2 s-1) and CO2 density in mole fraction of wet air (CO2), which was originally output on the datalogger as average CO2 density (mg m-2 s-1) and converted into µmol mol-1 using air temperature and pressure in post-processing. Webb, Pearman, and Leuning terms were applied to account for density fluctuations for water vapor and CO2. Sonic path length, high-frequency attenuation and sensor separation were accounted for according to Horst and Massman, and the oxygen absorption correction for the KH2O hygrometer was also applied. Latent and sensible heat fluxes were assumed to be unrealistic above an absolute value of 1000 W m-2 and were removed. Both carbon flux (FC) and carbon dioxide mole fraction in wet air (CO2) and were assumed to be unrealistic above 1000 µmol m-2 s-1 and 1000 µmol mol-1 respectively. Spikes in latent and sensible heat and carbon fluxes and densities (often due to periods of precipitation) were identified by computing the mean and standard deviation over a moving, overlapping window of 336 half-hours (7 days), similar to Shao et al., and were removed when the flux was more than 1.5 standard deviations from the moving window’s mean. While Vickers and Mahrt use a threshold of 3.5 standard deviations from the mean, a conservative value of 1.5 was chosen due to the noisy nature of over-lake data at this particular site. This process was repeated twice for latent and sensible heat, and carbon dioxide fluxes and densities and therefore it is possible that some real, realistic data was filtered out in this process. No detrending was performed. As per AmeriFlux standards, no friction velocity (USTAR, m s-1) filtering was performed.

Spence, Chris↗