SEARCH · Engineering Papers

Results for “normalization method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Towards quantifying the impact of triaxiality on optical signatures of galaxy clusters: weak lensing and galaxy distributions

ABSTRACT We present observational evidence of the impact of triaxiality on radial profiles that extend to 40 Mpc from galaxy cluster centres in optical measurements. We perform a stacked profile analysis from a sample of thousands of nearly relaxed galaxy clusters from public data releases of the Dark Energy Survey and the Dark Energy Camera Legacy Survey. Using the central galaxy elliptical orientation angle as a proxy for galaxy cluster orientation, we measure cluster weak lensing and excess galaxy density axis-aligned profiles, extracted along the central galaxy’s major or minor axes on the plane of the sky. Our measurements show a ≳ 2σ–3σ difference per radial bin between the normalized axis-aligned profiles. The profile difference between each axis-aligned profile and the azimuthally averaged profile ($\sim \pm 10\,\rm per\ cent-20~{{\ \rm per\ cent}}$ along major/minor axis) appears inside the clusters (∼0.4 Mpc) and extends to the large-scale structure regime (∼10–20 Mpc). The magnitude of the difference appears to be relatively insensitive to cluster richness and redshift, and extends further out in the weak lensing surface mass density than in the galaxy overdensity. Looking forward, this measurement can easily be applied to other observational or simulation data sets and can inform the systematics in cluster mass modelling related to triaxiality. We expect imminent upcoming wide-area deep surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, to improve our quantification of optical signatures of cluster triaxiality.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

AAPM Truth‐based CT (TrueCT) reconstruction grand challenge

Background: This Special Report summarizes the 2022, AAPM grand challenge on Truth-based CT image reconstruction. Purpose: To provide an objective framework for evaluating CT reconstruction methods using virtual imaging resources consisting of a library of simulated CT projection images of a population of human models with various diseases. Methods: Two hundred unique anthropomorphic, computational models were created with varied diseases consisting of 67 emphysema, 67 lung lesions, and 66 liver lesions. The organs were modeled based on clinical CT images of real patients. The emphysematous regions were modeled using segmentations from patient CT cases in the COPDGene Phase I dataset. For the lung and liver lesion cases, 1–6 malignant lesions were created and inserted into the human models, with lesion diameters ranging from 5.6 to 21.9 mm for lung lesions and 3.9 to 14.9 mm for liver lesions. The contrast defined between the liver lesions and liver parenchyma was 82 ± 12 HU, ranging from 50 to 110 HU. Similarly, the contrast between the lung lesions and the lung parenchyma was defined as 781 ± 11 HU, ranging from 725 to 805 HU. For the emphysematous regions, the defined HU values were −950 ± 17 HU ranging from −918 to −979 HU. The developed human models were imaged with a validated CT simulator. The resulting CT sinograms were shared with the participants. The participants reconstructed CT images from the sinograms and sent back their reconstructed images. Further, the reconstructed images were then scored by comparing the results against the corresponding ground truth values. The scores included both task-generic (root mean square error [RMSE] and structural similarity matrix [SSIM]), and task-specific (detectability index [d’] and lesion volume accuracy) metrics. For the cases with multiple lesions, the measured metric was averaged across all the lesions. To combine the metrics with each other, each metric was normalized to a range of 0 to 1 per disease type, with “0” and “1” being the worst and best measured values across all cases of the disease type for all received reconstructions. Results: The True-CT challenge attracted 52 participants, out of which 5 successfully completed the challenge and submitted the requested 200 reconstructions. Across all participants and disease types, SSIM absolute values ranged from 0.22 to 0.90, RMSE from 77.6 to 490.5 HU, d’ from 0.1 to 64.6, and volume accuracy ranged from 1.2 to 753.1 mm3. The overall scores demonstrated that participant “A” had the best performance in all categories, except for the metrics of d’ for lung lesions and RMSE for liver lesions. Participant “A” had an average normalized score of 0.41 ± 0.22, 0.48 ± 0.32, and 0.42 ± 0.33 for the emphysema, lung lesion, and liver lesion cases, respectively. Conclusions: The True-CT challenge successfully enabled objective assessment of CT reconstructions with the unique advantage of access to a diverse population of diseased human models with known ground truth. This study highlights the significant potential of virtual imaging trials in objective assessment of medical imaging technologies.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

DOE OSTI2022

LAI, EVI, NDVI, and kNDVI in 23 pantropical forests affected by 21 cyclones

Statement of purpose: Cyclones alter the function and composition of tropical forests, making effects of intensifying cyclones on carbon-rich forests a critical topic of study. Here, we quantified cyclone-induced damage and recovery of 21 cyclone disturbances affecting 23 pantropical forest sites between 1988-2017 utilizing leaf area index (LAI), enhanced vegetation index (EVI), normalized difference vegetation index (NDVI), and transformed NDVI (kNDVI) values from Google Earth Engine. Field observations collected in a meta-analysis (Bomfim et al., 2022, in review) were used to ground-truth and test effects of soil resource availability and disturbance factors on damage and recovery. This meta-analysis also served as the basis to begin vegetation index extraction, utilizing unique site and date combinations, from tropical forests effect by cyclone disturbances. We began collecting NDVI (5km resolution) from the NOAA Climate Data Record (CDR) of AVHRR Normalized Difference Vegetation Index (NDVI), Version 5 data product (Vermote, 2019) for all case studies included, 42. Next, we began extracting Landsat data from Landsat 4, 5, and 8, courtesy of the U.S. Geological Survey, in search of higher resolution data. We selected a 3 by 3 Landsat pixel area, leading to a 90m resolution data extraction. The specific imagery used includes Landsat 4 USGS Landsat 4 TM Collection 1 Tier 1 TOA (top of atmosphere) Reflectance, Landsat 5 USGS Landsat 5 TM (thematic mapper) Collection 1 Tier 1 TOA Reflectance, and Landsat 8 USGS Landsat 8 Collection 1 Tier 1 TOA Reflectance. Within Google Earth Engine, we selected the date and location (latitude and longitude), calculated NDVI, kNDVI, and EVI utilizing Landsat bands (see metadata_NGEE-tropics_cyclones), and extracted post- and pre-cyclone values for each case study to calculate cyclone-induced change in the vegetative index. Due to limited spatial resolution of Landsat remote sensing data, MODIS products were investigated next. First, the MOD13Q1.006 Terra Vegetation Indices 16-Day Global 250m product was used to extract 250m EVI and NDVI (Didan, 2015) and then the MCD15A3H.006 MODIS Leaf Area Index/FPAR 4-Day Global 500m product product was used to extract LAI 500m (Myneni et al., 2015). Pre- and post-cyclone values, change in the vegetative index, and standard deviation for all values are included in the main csv (see case_study_data.csv) for all vegetative indices collected, including LAI 500m, EVI 250m, NDVI 250m, NDVI 90m, kNDVI 90m, EVI 90m, and NDVI 5km. Lastly, recovery values were calculated utilizing a standardization method (see metadata_NGEE-tropics_cyclones) and recovery values for MODIS (see MODIS_recovery.csv) and Landsat (Landsat_recovery.csv) data are included.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Flammability and dispersion of tritium in confined release scenarios

Ignition of a flammable tritium-air mixture is the most probable means to produce the water form (T 2 O or HTO), which is more easily absorbed by living tissue and is hence ~10,000 times more hazardous to human health when uptake occurs compared to the gaseous form (T 2 or HT; per Mishima and Steele, 2002). Tritium-air mixtures with T 2 concentrations below 4 mol% are considered sub-flammable and will not readily convert to the more hazardous water form. It is therefore desirable from a safety perspective to understand the dispersion behavior of tritium under different release conditions, especially since tritium is often stored in quantities and pressures much lower than is typical for normal hydrogen. The formation of a flammable layer at the ceiling is a scenario of particular concern because the rate of dispersion to nonflammable conditions is slowest in this configuration, which maximizes the time window over which the flammable tritium may encounter an ignition source. This report describes the processes of buoyant rise and dispersion of tritium. Accumulation of flammable concentrations of tritium next to the ceiling is a common safety concern for hydrogen, but this situation can only occur if dispersion rates are slow with respect to rates of release and rise. Theory and simulations demonstrate that buoyancy does not cause regions with flammable concentrations to form within buildings from sources that have previously been mixed to sub-flammable concentrations. A simulated series of tritium release events with their associated dispersion behavior are reported herein; these simulations apply computational fluid dynamics to rooms with three different ceiling heights and a variety of tritium release rates. Safety related quantities from these simulations are reported, including the mass and volume of tritium occurring in a flammable mixture, the presence or absence of a flammable layer at the ceiling, and the time required for dispersion to nonflammable conditions after the end of the tritium release event. These safety metrics are influenced by the magnitude and rate of the tritium release with respect to the air volume in the room and also the momentum of the plume or jet with respect to the ceiling height. Several screening criteria are recommended to assess whether a specific tritium release scenario is likely to form a flammable layer at the ceiling. The methods and results in this modeling study have applicability to explosion safety analysis for other buoyant flammable gases, including the lighter isotopes of hydrogen.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2023

X-ray absorption spectroscopy study of Mn reference compounds for Mn speciation in terrestrial surface environments

Abstract X-ray absorption spectroscopy (XAS) offers great potential to identify and quantify Mn species in surface environments by means of linear combination fit (LCF), fingerprint, and shell-fit analyses of bulk Mn XAS spectra. However, these approaches are complicated by the lack of a comprehensive and accessible spectrum library. Additionally, molecular-level information on Mn coordination in some potentially important Mn species occurring in soils and sediments is missing. Therefore, we investigated a suite of 32 natural and synthetic Mn reference compounds, including Mn oxide, oxyhydroxide, carbonate, phosphate, and silicate minerals, as well as organic and adsorbed Mn species, by Mn K-edge X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopy. The ability of XAS to infer the average oxidation state (AOS) of Mn was assessed by comparing XANES-derived AOS with the AOS obtained from redox titrations. All reference compounds were studied for their local (<5 Å) Mn coordination environment using EXAFS shell-fit analysis. Statistical analyses were employed to clarify how well and to what extent individual Mn species (groups) can be distinguished by XAS based on spectral uniqueness. Our results show that LCF analysis of normalized XANES spectra can reliably quantify the Mn AOS within ~0.1 v.u. in the range +2 to +4. These spectra are diagnostic for most Mn species investigated, but unsuitable to identify and quantify members of the manganate and Mn(III)-oxyhydroxide groups. First-derivative XANES fingerprinting allows the unique identification of pyrolusite, ramsdellite, and potentially lithiophorite within the manganate group. However, XANES spectra of individual Mn compounds can vary significantly depending on chemical composition and/or crystallinity, which limits the accuracy of XANES-based speciation analyses. In contrast, EXAFS spectra provide a much better discriminatory power to identify and quantify Mn species. Principal component and cluster analyses of k2-weighted EXAFS spectra of Mn reference compounds implied that EXAFS LCF analysis of environmental samples can identify and quantify at least the following primary Mn species groups: (1) Phyllo- and tectomanganates with large tunnel sizes (2 × 2 and larger; hollandite sensu stricto, romanèchite, todorokite); (2) tectomanganates with small tunnel sizes (2 × 2 and smaller; cryptomelane, pyrolusite, ramsdellite); (3) Mn(III)-dominated species (nesosilicates, oxyhydroxides, organic compounds, spinels); (4) Mn(II) species (carbonate, phosphate, and phyllosilicate minerals, adsorbed and organic species); and (5) manganosite. All Mn compounds, except for members of the manganate group (excluding pyrolusite) and adsorbed Mn(II) species, exhibit unique EXAFS spectra that would allow their identification and quantification in mixtures. Therefore, our results highlight the potential of Mn K-edge EXAFS spectroscopy to assess bulk Mn speciation in soils and sediments. A complete XAS-based speciation analysis of bulk Mn in environmental samples should preferably include the determination of Mn valences following the “Combo” method of Manceau et al. (2012), EXAFS LCF analyses based on principal component and target transformation results, as well as EXAFS shell-fit analyses for the validation of LCF results. For this purpose, all 32 XAS reference spectra are provided in the Online Materials1 for further use by the scientific community.

Geochemistry & Geophysics↗

DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Pairing glue in cuprate superconductors from the self-energy revealed via machine learning

Recently, machine learning was applied to extract both the normal and the anomalous components of the self-energy from photoemission data at the antinodal points in Bi– based cuprate high-temperature superconductors. It was argued that both components do show prominent peaks near 50 meV, which hold information about the pairing glue, but the peaks are hidden in the actual data, which measure only the total self-energy. We analyze the self-energy within an effective fermion-boson theory. We show that soft thermal fluctuations give rise to peaks in both components of the self-energy at a frequency comparable to the superconducting gap, while they cancel in the total self-energy, all irrespective of the nature of the pairing boson. Furthermore, in the quantum limit T → 0 prominent peaks survive only for a very restricted subclass of pairing interactions. We argue that the way to potentially nail down the pairing boson is to determine the thermal evolution of the peaks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Response of thermal and fast-ion transport to beam ion population, rotation and T e / T i in the DIII-D steady state hybrid scenario

The thermal and fast-ion transport properties of DIII-D steady-state hybrid discharges with normalized beta β N ≳ 3 are studied at low injected torque and an increased electron to ion temperature ratio T e / T i . Linear stability analysis performed with the TGLF turbulent code indicates that a high-k mode is usually dominant at smaller radii, whereas a low-k mode is usually dominant at larger radii in these plasmas. A reduction in the net injected torque from 8.6 to 4.3 N-m leads to reduced E × B shear and hence, an enhanced turbulence that was observed on the Doppler backscattering diagnostic and was also computed with TGLF. As T e / T i in the core was increased from 0.57 to 0.66 by adding electron cyclotron current drive (ECCD) to these plasmas, higher levels of transport are observed with increased high-k modes indicated by TGLF. The fast-ion transport level varied over an order of magnitude in these discharges depending on whether Alfvén eigenmodes, fishbones, or no instabilities were observed. Hybrid plasmas with fishbones have decreased fast-ion transport, compared to plasmas with Alfvén eigenmodes, since they are resonant with a smaller portion of phase space and their resonance is farther from the wall. This reduction in fast-ion transport with ECCD mitigates the increase in turbulent transport, resulting in higher performance than expected during strong electron heating. Similarly, the lowest fast-ion transport was observed in the low torque plasma, which also led to better than expected performance at this torque value. Finally, the thermal and fast-ion transport changes observed as the torque/rotation and T e / T i are varied indicate possible methods for transferring this scenario to a reactor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Chiral Topological Surface States on a Finite Square Photonic Crystal Bounded by Air

Chiral, topologically protected, photonic surface states can be found at the boundary between gyrotropic photonic crystals where a changing magnetic field induces different topology across the interface. Typically, photonic crystals with either a suitable band structure on both sides of the interface to provide a band gap and evanescent decay of the surface states away from the interface, or an outer layer with engineered material properties is required. In this paper, we show the emergence of topological, unidirectional surface states at the termination of finite gyrotropic photonic crystals with a simple square lattice and C 4 rotational symmetry bounded by a vacuum, eliminating the need for an outside layer to enable chiral surface modes. Here, we start from an infinite, time-reversal-symmetry-breaking photonic crystal with a band gap associated with bands with nonzero Chern numbers, different from all-zero Chern numbers in air. We then modify the photonic crystal to move this band gap below the light line, while maintaining the Chern-number discontinuities. Band-structure calculations for a supercell approximating a photonic crystal finite in the direction normal to the surface demonstrate the existence, dispersion, and chirality of the surface mode. Extensive direct scattering calculations for a point source and spatial Fourier analysis further reveal a unidirectional free-space topological surface state, which propagates counterclockwise around the surface of a finite photonic crystal, providing a nearly foolproof way to cross-check the surface-mode band structure unaffected by backscattering from local defects. Additionally, scattering simulations allow an independent characterization of the state dispersion and unveil the robustness of the topological plasmonic mode propagation around the 90° bends of the structure, being due to only radiation leakage. In contrast to buried topological surface states, the observed surface modes at the photonic crystal–air interface have the advantage of being accessible to the outside world, allowing one to take advantage of the defect-tolerant backscattering-free surface modes to engineer emission from photonic crystal surfaces into arbitrary free-space beam shapes and directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Density of States and Spectral Function of a Superconductor out of a Quantum-Critical Metal

We analyze the validity of a quasiparticle description of a superconducting state above a metallic quantum-critical point (QCP). A normal state at a QCP is a non-Fermi liquid with no coherent quasiparticles. A superconducting order gaps out low-energy excitations, except for a sliver of states for non-s-wave gap symmetry, and at a first glance, restores coherent quasiparticle behavior. We argue that this does not necessarily hold as the fermionic self-energy may remain singular above the gap edge. This singularity gives rise to markedly non-BCS behavior of the density of states and to the appearance of a nondispersing mode at the gap edge in the spectral function. We analyze the set of quantum-critical models with an effective dynamical four-fermion interaction V(Ω) ∝ 1/Ω γ , where Ω is a frequency of a boson, which mediates the interaction. We show that coherent quasiparticle behavior in a superconducting state holds for γ < 1/2, but breaks down for larger γ. Here, we discuss signatures of quasiparticle breakdown and compare our results with the photoemission data for Bi2201 and Bi2212.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Advancement of Commercial Intumescent Expanding Foams for Deactivation and Decommissioning in the Nuclear Sector - 20198

Florida International University (FIU), in collaboration with The Department of Energy's Office of Environmental Management (DOE-EM), Savannah River National Laboratory (SRNL), and sites across the Savannah River complex, have identified an operational requirement for a fixative technology that is intended to immobilize and/or isolate residual contamination within a 3-dimensional space. Fixation of radiological contamination can reduce worker risk and mitigate potentially hazardous conditions, however nearly every marketed contamination fixative has been found to be flammable; a significant concern in radiological facilities. Coupled with this, industry fixatives are normally used as a thin coating which can present problems when attempting to stabilize irregular geometry or areas that are difficult to access whilst ensuring full coverage. The technical evaluation and advancement of commercial-off-the-shelf (COTS) polyurethane foams has yielded a down-selected candidate that shows potential in meeting the requirements to support deactivation and decommissioning activities. Several performance criteria have been established and tested to progress the technology readiness level towards an active field demonstration (TRL-7). Such criteria include: mechanical failure limits, adhesive and cohesive properties, thermal/fire resilience, determining thermal behavior, ability to immobilize contamination, and a means of non-destructive evaluation of applications. The test scenario examined was targeted towards an application for decommissioning nuclear pipework, in which the down-selected polyurethane foam would act as a barrier to segregate pipework and mitigate the potential for release during cutting, packaging, and storage operations. Testing carried out at SRNL included: mechanical evaluation of tensile, compressive, and adhesion strength by dynamic mechanical analysis (DMA), as well as thermogravimetric analysis (TGA). FIU examined the foam's fixative properties by utilizing phosphorescent europium-dysprosium doped strontium aluminate powder to investigating the extent to which contamination can be immobilized. FIU has also exploited previous successes in the field of intumescent technologies to assess the down-selected foam's tolerance to an extreme fire scenario, while maintaining the ability to effectively mitigate a contamination release. Parallel to this, extensive thermal investigations were carried out to determine the upper boundary of anticipated heat generation during the curing process as heat generation has the potential to compromise rubber parts of contaminated enclosures. These investigations subsequently yielded a promising method for a non-destructive application evaluation by means of infrared thermography. Utilizing the high sensitivity of modern IR cameras, coupled with the heat generated during the curing process of the polyurethane foam, FIU has been exploring the concept of monitoring the external pipe surface for indications of an irregular or abnormal application, thus informing operational decision making. The testing carried out utilized several current 'best fit' ASTM standards, which serve as helpful guidelines for testing, however, a precise definition of the operational parameters and requirements is still necessary. With continued collaboration with SRNL, FIU aims refine said definitions and develop new standards by which this, and other decommissioning technologies, can be accredited by relevant standards based testing. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Chromatographic Modeling of 1-Foot Detritiation Test Column Results for the Optimization of Regeneration Time

This report documents the methods and results of a modeling study on a series of data generated by 1-foot test columns for the Water Detritiation project. In order to compliment modeling work already in progress on this project, a purely empirical curve-fitting approach is used in this study to determine chromatographic performance indicators such as retention time, curve width, and the number of chromatographic stages. The data used in this study comes from four sets of chromatographic test separations performed in March through May of 2017 using 1-foot columns. All are using the same flow rates and temperature (room temperature, nominally 25 C°). The first set of data (referred to in this report as the Column 1 data) is a series of twelve breakthrough curve measurements. Prior to each of these twelve runs, the column was regenerated by flowing hydrogen gas through the column, for increasing amounts of time, ranging from 2 hours to 25 hours. The second set of data (referred to as the Column 2 data) is comprised of five measurements of breakthrough curves measured in similar fashion as the Column 1 data, across a range of regeneration times from 2 to 21 hours. The third set of data (referred to as the Purge data) is of three breakthrough curves, measured after a regeneration time of 8 hours, which takes place after a nitrogen purge of 15, 30, or 45 minutes. The fourth data set is a series of five breakthrough curve measurements using a 2 nd batch of stationary phase, taken after regeneration times ranging from 2 to 20 hours. The overall modeling approach is as follows. First, the data is pre-processed for formatting, artifact removal, and is normalized. The data is then fit by integrated versions of four common chromatographic models. The best model is chosen based on the predictive performance of the model, and chromatographic performance values are calculated based on the best model’s fit parameters. Each of these processes are detailed below.

42 ENGINEERING↗

DOE OSTI2020

Target normal sheath acceleration with a large laser focal diameter

In this work, the dependence of the laser-driven ion acceleration from thin titanium foils in the Target Normal Sheath Acceleration (TNSA) regime on target and laser parameters is explored using two dimensional particle-in-cell simulations. The oblique incidence (θ L = 45°) and large focal spot size (w 0 = 40μm) are chosen to take an advantage of quasi one-dimensional geometry of sheath fields and effective electron heating. This interaction setup also reveals low and achromatic angular divergence of a proton beam. It is shown that the hot electron temperature deviates from the ponderomotive scaling for short laser pulses and small pre-plasmas. This deviation is mainly due to the laser sweeping, as the short duration laser pulse each moment in time effectively heats only a fraction of a focal spot on the foil. This instantaneous partial heating results in an electron temperature deviation from the ponderomotive scaling and, thus, lower maximum proton energies than it could have been expected from the TNSA theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗

DOE OSTI2024

Tuning superconductivity and spin-vortex instabilities in CaKFe 4 As 4 through in-plane antisymmetric strains

Lattice strains of appropriate symmetry have served as an excellent tool to explore the interaction of superconductivity in the iron-based superconductors with orthorhombic-nematic and stripe spin-density-wave (SSDW) order. In this Letter, we contribute to a broader understanding of the coupling of strain to superconductivity and competing normal-state orders by studying CaKFe 4 As 4 under large, in-plane strains of B 1 g and B 2 g symmetry. In contrast to the majority of iron-based superconductors, pure CaKFe 4 As 4 exhibits superconductivity with a relatively high transition temperature of T c ∼ 35 K in proximity of a noncollinear, tetragonal, hedgehog spin-vortex crystal (SVC) order. Through experiments and calculations, we demonstrate an anisotropic in-plane strain response of T c and the favored SVC configuration in CaKFe 4 As 4 . This supports a scenario, in which the change in spin fluctuations dominates the strain response of superconducting T c . Overall, by suggesting moderate B 2 g strains as an effective parameter to change the stability of SVC and SSDW, we outline a pathway to a unified phase diagram of iron-based superconductivity. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

On Anomaly Detection for Transactive Energy Systems with Competitive Market

Two anomaly-detection criteria are proposed for transactive energy systems with competitive markets. Participants of transactive energy systems seek an optimal power allocation through hybrid economic-control methods to facilitate the integration of various types of distributed energy resources to power distribution systems. In transactive energy systems, every participant is assumed to be a rational entity, and consumers have diminishing marginal utility and suppliers have increasing marginal cost. With the first proposed anomaly-detection criterion, the monotonicity of marginal cost and marginal utility are examined. The impact of line flow constraints is also taken into consideration. Then, the second anomaly-detection criterion is proposed for TESs with marginal cost and marginal utility which change faster than a certain rate. The second criterion is more accurate than the first one for TES with marginal cost and marginal utility which change faster than a certain rate, but it requires the knowledge of that rate. Neither criteria requires more data than those necessary to find the optimal power allocation and the market-clear price in a transactive energy system. Therefore, the proposed criteria do not disclose any more data than necessary. As the monotonicity of marginal cost and marginal utility in a TES with competitive markets results in convex objective functions in an optimization problem and strongly convex ones when marginal cost and marginal utility changes faster than a certain rate, the two detection criteria are also applicable to anomaly detection of general convex optimization problems. Simulations are carried out to show the efficacy of the proposed criteria to detect anomalies caused by cyberattacks.

Wang, Peng↗

DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗