Search NASA⌕ Search

SEARCH · Search NASA

Results for “column generation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Creation of Synthetic Electric Grids (SPP/MISO) Supporting PERFORM (Final Report)

Over the course of the project, two “realistic but not real” synthetic transmission-level grid models over the SPP-MISO and ERCOT footprints were created to provide more realistic data and increase the reliability and resiliency of the grids under a variety of scenarios. The synthetic ERCOT transmission grid is compatible with the distribution grid developed in collaboration with NREL. All generators are based on the EIA 860 data and a column with EIA plant code and Gen ID is added to generators of both grids so that they can be easily mapped. The improvements are also made to electric grids including N-1 contingencies with some remedial actions, improving the transmission lines to avoid lines in lakes, including an HVDC line to the SPP-MISO case, providing several generator parameters and their temporal constraints that were not included in EIA 860 form, generators’ cost curves, load offer curves, adding phase shifters and tap changers with impedance correction tables, adding reactive power control and partitioning the grids into active and reactive reserve zones and determine different types of the required reserve for each zone. Hourly load time series at the bus level were generated to create scenarios for solving power flow in different loading conditions. Weather measurement information and the models of renewable generators are used to directly include the impact of weather on the grids. Based on a variety of load and weather conditions the grids are improved to accommodate different conditions. The ERCOT 7k-bus grids were also modeled for the year 2030 with predicted improvements in renewable resources. The renewable generation model was also improved with historic weather data included. The impact of electric vehicles on the ERCOT grid is also modeled.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Recovery of high specific activity molybdenum-99 from accelerator-induced fission on low-enriched uranium for technetium-99m generators

Abstract A new process was developed to recover high specific activity (no carrier added) 99 Mo from electron-accelerator irradiated U 3 O 8 or uranyl sulfate targets. The process leverages a novel solvent extraction scheme to recover Mo using di(2-ethylhexyl) phosphoric acid following uranium and transuranics removal with tri- n -butyl phosphate. An anion-exchange concentration column step provides a final purification, generating pure 99 Mo intended for making 99 Mo/ 99m Tc generators. The process was demonstrated with irradiated uranium targets resulting in more than 95% 99 Mo recovery and without presence of fission products or actinides in the product.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗
DOE OSTI2022

The Earth Model Column Collaboratory (EMC 2 ) v1.1: an open-source ground-based lidar and radar instrument simulator and subcolumn generator for large-scale models

Abstract. Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

58 GEOSCIENCES↗
DOE OSTI2022

The influence of random packed column parameters on the liquid holdup and interfacial area

Abstract Carbon dioxide capture via solvent absorption in packed columns has emerged as a potential technology to mitigate coal‐fired power plant CO 2 emissions. Parameters, including packing types, solvent properties, and operating conditions, could potentially affect the packed column CO 2 capture efficiency. To understand the importance of those parameters and help packed column optimization, a design of experiments (DoEs) method was proposed to generate input parameter matrix. Combined with multiphase computational fluid dynamics (CFD), the random packed column parameter influence on the liquid holdup and interfacial area can be efficiently investigated. Surrogate‐based sensitivity analysis shows that the solvent flow rate and contact angle are key factors dictating liquid holdup and interfacial area. On the other hand, solvent viscosity has a marginal impact on the interfacial area. The sensitivity scores were calculated for each input parameter to guide the selection of dimensionless numbers for the liquid holdup and interfacial area correlation development.

42 ENGINEERING↗
DOE OSTI2020

Recycling Isotopically Enriched Mo from the RadioGenix Waste Stream

NorthStar’s RadioGenix generator utilizes an aqueous biphasic extraction chromatography column to separate Tc from low-specific-activity 99 Mo. During this separation process, a fraction of the source Mo material can be found in the waste. This is not of concern when natural Mo targets are used. However, the introduction of isotopically enriched 98 Mo for neutron capture or 100 Mo for accelerator-based production requires recovery of this valuable enriched Mo material. The generator waste is absorbed and solidified in a superabsorbent hydrogel polymer. In this contribution, we report results from our small-scale experimental study on development of a separation process for recovery of molybdenum from the polymeric gel waste using common chemicals and simple chemical techniques. Our investigation confirms that the acidic deswelling of the alkaline hydrogel with concentrated hydrochloric acid converts the occluded molybdate to a neutral molybdenyl dichloride species that is freely released from the gel. Mo can then be concentrated on a commercial tri-n-butyl phosphate resin, eluted with ammonium hydroxide, and crystallized as ammonium heptamolybdate for conversion back to Mo metal. Typical Mo recoveries are 90-95%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A comparison of the power potential for surface- and seabed-deployed tidal turbines in the San Juan Archipelago, Salish Sea, WA

The San Juan archipelago lies along the axis of tidal movement between Straits of Juan de Fuca and Strait of Georgia in the Salish Sea. The amplitude of the tidal exchange produces significant tidal currents between the islands, as well as in Rosario and Haro Straits. These are of interest as a future source of electrical power generation, given the islands’ dependence on electricity supply by a subsea cable from the mainland. Here, we evaluate the tidal current energy potential in this region through a re-analysis of measurements collected by the National Ocean Service (NOS) and a high-resolution numerical model. Given the considerable variations in water depth and vertical velocity profiles across candidate tidal energy sites, we consider the trade-offs between tidal turbines deployed from a floating platform and those anchored to the seabed. Measurement re-analysis indicates several locations that could support tidal current power generation by MW-scale turbines with an acceptable balance between turbine size, rated power, and capacity factor. Even for relatively large (30 m) turbine diameters, surface-deployed turbines would be expected to produce up to 30 % more electricity than the same turbine deployed near the seabed due to vertical shear, with this difference increasing for smaller diameter turbines. A benchmark comparison at the measurement locations with the highest power generation potential shows relatively good model fidelity, though, even with relatively high resolution, time-average power density disagreements of ± 50 % persist throughout the water column. Nonetheless, the simulation identifies several locations with power generation potential more than twice as high as locations in the measurement re-analysis. These sites were either not surveyed by NOS or excluded due to data quality issues. Altogether, these results demonstrate the complimentary of models and measurements for tidal energy site assessment and suggest that tidal currents could be an important source of electricity generation in the San Juan archipelago.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Insights on sources and formation mechanisms of liquid-bearing clouds over MOSAiC examined from a Lagrangian framework

Understanding Arctic stratiform liquid-bearing cloud life cycles and properly representing these life cycles in models is crucial for evaluations of cloud feedbacks as well as the faithfulness of climate projections for this rapidly warming region. Examination of cloud life cycles typically requires analyses of cloud evolution and origins on short time scales, on the order of hours to several days. Measurements from the recent Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition provide a unique view of the current state of the central Arctic over an annual cycle. Here, we use the MOSAiC radiosonde measurements to detect liquid-bearing cloud layers over full atmospheric columns and to examine the cloud-generating air masses’ properties. We perform 5-day (120 h) back-trajectory calculations for every detected cloud and cluster them using a unique set of variables extracted from these trajectories informed by ERA5 reanalysis data. This clustering method enables us to separate between the air mass source regions such as ice-covered Arctic and midlatitude open water. We find that moisture intrusions into the central Arctic typically result in multilayer liquid-bearing cloud structures and that more than half of multilayer profiles include overlying liquid-bearing clouds originating in different types of air masses. Finally, we conclude that Arctic cloud formation via prolonged radiative cooling of elevated stable subsaturated air masses circulating over the Arctic can occur frequently (up to 20% of detected clouds in the sounding data set) and may lead to a significant impact of ensuing clouds on the surface energy budget, including net surface warming in some cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comprehensive Review of DPF-based Flash Neutron Radiography Viability

MJOLNIR (MegaJOule Neutron Imaging Radiography) is a Dense Plasma Focus (DPF) being developed by LLNL as a prototype to assess the viability of a DPF as a flash neutron radiography source. To date, MJOLNIR has discharged up to 1.3 MJ of stored energy into a deuterium plasma load and achieved yields up to 1.2 × 10 12 neutrons per discharge. The MJOLNIR pulse length already meets the preliminary requirements for flash neutron radiography and we demonstrate in this report plausible paths forward to meet the yield requirement in a deuterium plasma. MJOLNIR’s neutron spot size has been characterized to be a factor of two greater than radiography requirements and presents the greatest challenge. A path forward to sufficiently shrinking spot size is outlined in this report, using a combination of gas dopants and a smaller anode implosion radius. A DPF is a plasma device with coaxial electrodes whose discharge ends with a stagnated hot and dense plasma column on-axis, at the tip of the central anode. Inside the MJOLNIR DPF plasma column, both thermal and beam target processes generate neutrons. To understand the underlying physics of neutron generation inside the DPF, we model the experiment using a combination of kinetic, fluid, and reduced-order models we have been developing for over ten years at LLNL. Our numerical tools are also pivotal to charting the path forward because they allow us to numerically test out which modifications would improve the DPF’s characteristic output.

42 ENGINEERING↗
DOE OSTI2024

Evaluation of Radiography for TRISO Buffer Layer Density Measurement

Tristructural isotropic (TRISO) fuel particles consist of a central uranium-bearing kernel and a series of coating layers designed to retain fission products and to ensure fuel performance. Several parameters such as thickness and density must be measured for these coating layers to show that they conform with fuel specifications. Current methods for measuring the density of pyrolytic carbon and silicon carbide layers (liquid gradient density column) and the buffer layer (mercury porosimetry) generate Resource Conservation and Recovery Act (RCRA) radiological-mixed waste. In addition, measurement of buffer and inner pyrolytic carbon layer densities require hot sampling or interrupted coating runs and the mercury porosimetry method used for buffer density measurement only measures the mean buffer density, not the interparticle distribution. A new approach has been evaluated to measure the density of coating layers in TRISO particles based on the dependence of x-ray attenuation in radiographs on material density. This method does not generate RCRA mixed waste, measures density on a particle-by-particle basis, and in principle is capable of measuring the density of all coating layers in a single process. Initial results using thinned TRISO particle sections to evaluate radiography measurement of density as a quality control characterization method are reported herein. In this work, the primary focus is on measurement of the density of the buffer layer; however, with appropriate calibration the method should be applicable to other coating layers. Improvements to the initial method and a full demonstration of the method on the remaining coating layers may be pursued as a future effort.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Low-synch Gram–Schmidt with delayed reorthogonalization for Krylov solvers

The parallel strong-scaling of iterative methods is often determined by the number of global reductions at each iteration. Low-synch Gram-Schmidt algorithms are applied here to the Arnoldi algorithm to reduce the number of global reductions and therefore to improve the parallel strong-scaling of iterative solvers for nonsymmetric matrices such as the GMRES and the Krylov-Schur iterative methods. In the Arnoldi context, the factorization is "left-looking" and processes one column at a time. Among the methods for generating an orthogonal basis for the Arnoldi algorithm, the classical Gram-Schmidt algorithm, with reorthogonalization (CGS2) requires three global reductions per iteration. A new variant of CGS2 that requires only one reduction per iteration is presented and applied to the Arnoldi algorithm. Delayed CGS2 (DCGS2) employs the minimum number of global reductions per iteration (one) for a one-column at-a-time algorithm. The main idea behind the new algorithm is to group global reductions by rearranging the order of operations. DCGS2 must be carefully integrated into an Arnoldi expansion or a GMRES solver. Numerical stability experiments assess robustness for Krylov-Schur eigenvalue computations. Performance experiments on the ORNL Summit supercomputer then establish the superiority of DCGS2 over CGS2.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories

In this work, we present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples. Program summary: Program Titles: LeptonInjector and LeptonWeighter CPC Library link to program files: https://doi.org/10.17632/662gkpjfd9.1 Developer's repository links: https://github.com/icecube/LeptonInjector and https://github.com/icecube/LeptonWeighter Licensing provisions: GNU Lesser General Public License, version 3. Programming Language: C++11 External Routines: • Boost • HDF5 • nuflux (https://github.com/icecube/nuflux) • nuSQuIDS (https://github.com/arguelles/nuSQuIDS) • Photospline (https://github.com/icecube/photospline) • SuiteSparse (https://github.com/DrTimothyAldenDavis/SuiteSparse) Nature of problem: LeptonInjector: Generate neutrino interaction events of all possible topologies and energies throughout and around a detector volume. LeptonWeighter: Reweight Monte Carlo events, generated by a set of LeptonInjector Generators, to any desired physical neutrino flux or cross section. Solution method: LeptonInjector: Projected ranges of generated leptons and the extent of the detector, in terms of column depth, are used to inject events in and around the detector volume. Event kinematics follow distributions provided in cross section files. LeptonWeighter: Event generation probabilities are calculated for each Generator, which are then combined into a generation weight and used to calculate an overall event weight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Demonstration of a 228 Ra/ 228 Ac isotope generator

An isotope generator to produce 228 Ac from 228 Ra was developed using a cation exchange resin column eluted with an acetate-diethylenetriaminpentaacetic acid buffer. Here, the elution behavior of 228 Ac on the column was studied with solutions of various pHs to select the optimal conditions. Two identical isotope generators were eluted for 47 days with over 100 mL of eluant with no detectable 228 Ra breakthrough and high 228 Ac yields (∼95%). The separation requires only biocompatible reagents and is performed at pH 4.6, conditions suitable for radiopharmaceutical studies of Ac.

and nuclear chemistry↗
DOE OSTI2020

Current and future neutrino oscillation constraints on leptonic unitarity

The unitarity of the lepton mixing matrix is a critical assumption underlying the standard neutrino-mixing paradigm. However, many models seeking to explain the as-yet-unknown origin of neutrino masses predict deviations from unitarity in the mixing of the active neutrino states. Motivated by the prospect that future experiments may provide a precise measurement of the lepton mixing matrix, we revisit current constraints on unitarity violation from oscillation measurements and project how next-generation experiments will improve our current knowledge. With the next-generation data, the normalizations of all rows and columns of the lepton mixing matrix will be constrained to ≲10% precision, with the e-row best measured at ≲1% and the τ-row worst measured at ~10% precision. The measurements of the mixing matrix elements themselves will be improved on average by a factor of 3. We highlight the complementarity of DUNE, T2HK, JUNO, and IceCube Upgrade for these improvements, as well as the importance of ντ appearance measurements and sterile neutrino searches for tests of leptonic unitarity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2023

A novel data gaps filling method for solar PV output forecasting

This study proposes a modified gaps filling method, expanding the column mean imputation method and evaluated using randomly generated missing values comprising 5%, 10%, 15%, and 20% of the original data on power output. The XGBoost algorithm was implemented as a forecasting model using the original and processed datasets and two sources of solar radiation data, namely, Shortwave Radiation (SWR) from Advanced Himawari Imager 8 (AHI-8) and Surface Solar Radiation Downward (SSRD) from ERA5 global reanalysis data. Further, the accuracy of the two sets of forecasted power output was evaluated using Root Mean Square Error (RMSE) and Mean Absolute Error (MAE). Results show that by applying the proposed gap filling method and using SWR in forecasting solar photovoltaic (PV) output, the improvement in the RMSE and MAE values range from 12.52% to 24.30% and from 21.10% to 31.31%, respectively. Meanwhile, using SSRD, the improvement in the RMSE values range from 14.01% to 28.54% and MAE values from 22.39% to 35.53%. To further evaluate the accuracy of the proposed gap-filling method, the proposed method could be validated using different datasets and other forecasting methods. Future studies could also consider applying the said method to datasets with data gaps higher than 20%.

Energy & Fuels↗
DOE OSTI2024

Plasma Bubble Column Reactor: A High Throughput Reactor Design for Water Treatment

Plasma-based water treatment (PWT) is a promising technology that can degrade various emerging contaminants. However, PWT application on an industrially viable scale is hindered by the lack of an efficient reactor design that combines enhanced plasma-liquid contact with high liquid throughput. This work investigates the applicability of a bubble column gas-liquid contactor to PWT. A pulsed plasma bubble column reactor (PBCR) with a concentric rod-cylinder electrode configuration was used to correlate contaminant removal performance with the gas-liquid contact parameters of the bubble column. A surfactant, rhodamine B dye, and a nonsurfactant, caffeine, were used as model contaminants at µM concentration levels. The bubble column characteristics, i.e., gas holdup, bubble size distribution, and gas-liquid area, were measured as a function of superficial gas velocity using image-based methods. Degradation rates of both contaminants increased with gas flowrate. For caffeine, the increase was attributed to intensified bulk liquid mixing, while dye degradation increased due to the increased gas-liquid area. Ultimately, we show that bubble column contactors significantly improve the utilization of plasma-generated reactive species toward contaminant degradation by distributing them over a large contact area. As a result, a better match between the plasma species interfacial flux and the interfacial contaminant concentration leads to improved treatment energy efficiency. Furthermore, typical degradation energy efficiencies were ~10 g/kWh for caffeine and ~60 g/kWh for rhodamine B.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗