Search NASA⌕ Search

SEARCH · Search NASA

Results for “projected gradient 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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical Analysis of High-Order Modes in SRF Resonators for Particle Accelerators

Over the past decades, superconducting technology has rapidly evolved towards high accelerating gradients and low surface resistance, making it possible to operate particle accelerators with high average beam currents and large duty factors. However, RF losses due to coherent excitation of the HOM become the limiting factor for these regimes. Unlike the cavity operating mode, which is tuned separately, the HOM parameters can significantly vary from one cavity to another due to finite mechanical tolerances during the manufacturing process. Thus, it is of utmost importance to know the HOM parameter spread in advance in order to predict unexpected cryogenic losses, overheating of beam line components and maintain stable beam dynamics. In this paper, we present a method for generating cavity geometry with an arbitrary spread of mechanical imperfections and numerically evaluating HOM statistics. Knowing the spread of HOM parameters, we calculated the probability of resonant HOM losses in SRF accelerating cavities used in CW beam current machines such as the PIP-II and LCLS-II linacs, as well as for the SRF crab-cavity for the ILC project. Finally, we present experimental results of HOM spectra measurements in hundreds of 1.3 GHz cavities installed in LCLS-II cryomodules. Studying the effects of HOM excitation results in specifications of the SRF cavity and cryomodule and can significantly impact the efficiency and reliability of the machine operation.

Lunin, Andrei [Fermilab] (ORCID:0000000290096792)↗
DOE OSTI2024

Recent Advances in Assessing Environmental Effects of Marine Renewable Energy Around the World

Marine renewable energy (MRE) is increasingly of interest to coastal nations as a source of renewable energy that can support climate change mitigation goals as well as provide secure locally-produced energy for coastal and island communities. MRE extracts power from tidal streams, waves, ocean currents, run of rivers, and gradients in the ocean, with specialized devices developed and tested for each energy resource. Alongside development of MRE technologies and systems, first in Europe and then in North America, Australia, Asia, and other regions, it has been universally recognized that there is also a need to examine potential effects on marine animals, habitats, ecosystem processes, local communities and other sea users, to ensure that the MRE industry can be developed in a responsible and sustainable manner. This paper looks at the status of assessment and monitoring for potential environmental effects associated with MRE projects around the world. Over 80 projects were identified worldwide as having been tested, demonstrated, or commercially deployed with associated environmental monitoring. Five of the projects that represent tidal stream, wave, and run of river projects are examined in more detail to determine the types of data and information collected for those projects, the outputs of the monitoring campaigns, and the actions taken as a result of the data collection and analysis. Recommendations are provided for standardization of the monitoring approaches, instruments, and analysis methods at MRE project sites worldwide.

marine renewable energy, environmental effects, ri↗
DOE OSTI2016

TU-H-206-07: Assessment of Geometric Distortion in EPI with a SPAMM Tagged Acquisition

Purpose: Echo planar imaging (EPI) is known to exhibit gross geometric distortion caused by multiple factors, including B0 inhomgeneity and transient eddy currents. However, diffusion weighted (DW) EPI has become indispensable for diagnosis and therapy assessment. We propose a methodology for quantifying distortion in EPI sequences that does not require the use of dedicated spatial accuracy phantoms, enabling flexibility in phantom design for QA of distortion effects in EPI protocols. Methods: The proposed methodology utilizes a saturation technique known as Spatial Modulation of Magnetization (SPAMM) that tags the imaging subject with saturated grid lines. Originally intended for tracking cardiac motion, these grids are applied to assess differences between diffusion weighting directions and b-values, or against a more geometrically robust sequence such as fast spin echo (FSE). The saturation preparation sequence consists of binomially weighted (e.g. 1-3-3-1) pulses interleaved with gradient blips along the frequency encode direction, followed by the same sequence with gradient blips in the phase encode direction. Three phantoms were assessed with these sequences: a spherical head-sized phantom, a large shimming phantom, and a modified PET ACR phantom that included compartments of water, air, oil, and Teflon. Each phantom was acquired with three sequences using parameters from a clinically appropriate protocol (22 cm head or 46 cm abdomen): a conventional DW-EPI sequence (3 DW directions), and both the DW-EPI and FSE sequences with tagging. Differences in grid locations were visualized with minimum intensity projection between images, and measured using intersecting locations on the grids. Results: Grid lines were clearly visualized on tagged images and enabled quantification of distortions. Maximum eddy current induced errors of 10.8 to 14.8 mm were observed in areas away from isocenter with DW gradients applied in various directions. Conclusion: SPAMM tagging provides a promising mechanism for assessing spatial distribution of distortion in EPI sequences.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Radiopure ZnMoO{sub 4} scintillating bolometers for the LUMINEU double-beta experiment

The results of R&D of radiopure zinc molybdate (ZnMoO{sub 4}) based scintillating bolometers for the LUMINEU (Luminescent Underground Molybdenum Investigation for NEUtrino mass and nature) double-beta decay experiment are presented. A dedicated two-stage molybdenum purification technique (sublimation in vacuum and recrystallization from aqueous solutions) and an advanced directional solidification method (the low-thermal-gradient Czochralski technique) were utilized to produce high optical quality large mass (∼1 kg) ZnMoO{sub 4} crystal boules and first {sup 100}Mo (99.5%) enriched Zn{sup 100}MoO{sub 4} crystal scintillator (mass of ∼0.2 kg). Scintillating bolometers based on ZnMoO{sub 4} (≈ 0.33 kg) and Zn{sup 100}MoO{sub 4} (≈ 0.06 kg) scintillation elements and high purity Ge wafers were tested in the EDELWEISS set-up at the Modane Underground Laboratory (France). Long term low temperature tests demonstrate excellent detectors’ performance and effectiveness of the purification and solidification procedures for the achievement of high radiopurity of the material, in particular with a bulk activity of {sup 228}Th and {sup 226}Ra below 4 µBq/kg. The adopted protocol was used to produce for the first time a large volume Zn{sup 100}MoO{sub 4} crystal scintillator (mass of ∼1.4 kg, {sup 100}Mo enrichment is 99.5%) to search for neutrinoless double-beta decay of {sup 100}Mo in the framework of the LUMINEU project.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗
NASA NTRS1999

Project Description and Publications List for UAH CMMR

This research combines a state of the art X-ray Transmission Microscope, XTM, with a specially designed x-ray transparent horizontal Bridgman furnace to image (with resolutions up to 3 micrometers) the solidification of metal alloys in real-time. The objective is to obtain real-time dynamic data to provide direct measure of the solute profile in the liquid, phase coalescence and growth in the liquid, and the detailed interface morphology (e,g., dendrites and cells) during solidification. We are also enhancing the XTM data with precise solid-liquid interfacial temperature and the thermal gradient measurement techniques, and working on the application of this technology to the study of the fundamentals of solidification in microgravity. Over the last several years we have successfully imaged in real-time: interfacial-morphologies, phase growth, coalescence, incorporation of phases into the growing interface, and the solute boundary layer in the liquid at the solid-liquid interface. We have also measured true local growth rates and can evaluate segregation structures in the solid. Interfacial undercoolings are being measured either with a special Seebeck furnace or with micro-thermocouple arrays we are developing. These later techniques are presently being incorporated with the XTM furnace. This last year emphasized the investigation of the solute layer in the melt during solidification. Methods were developed to quantify the solute concentrations using x-ray absorption and to compare to predictions from simulations. In addition, work is being completed on a brass-board portable XTM that incorporates a vertical Bridgman furnace.

Kaukler, William F.↗
DOE OSTI2017

Quantum Mechanical Simulations of Complex Nanostructures for Photovoltaic Applications

A quantitative understanding of the electronic excitations in nanostructures, especially complex nanostructures, is crucial for making new-generation photovoltaic (PV) cells based on nanotechnology, which have high efficiency and low cost. Yet current quantum mechanical simulation methods are either computationally too expensive or not accurate and reliable enough, hindering the rational design of the nanoscale PV cells. The PI seeks to develop new methodologies to overcome the challenges in this very difficult and long-lasting problem, pushing the field forward so that electronic excitations can be accurately predicted for systems involving thousands of atoms. The primary objective of this project is to develop new approaches for electronic excitation calculations that are more accurate than traditional density functional theory (DFT) and are applicable to systems larger than what current beyond-DFT methods can treat. In this proposal, the PI will first address the excited-state problem within the DFT framework to obtain quasiparticle energies from both Kohn-Sham (KS) eigenvalues and orbitals; and the electron-hole binding energy will be computed based on screened Coulomb interaction of corresponding DFT orbitals. The accuracy of these approaches will be examined against many-body methods of GW/BSE and quantum Monte Carlo (QMC). The PI will also work on improving the accuracy and efficiency of the GW/BSE and QMC methods in electronic excitation computations by using better KS orbitals obtained from orbital-dependent DFT as inputs. Then an extended QMC database of ground- and excited-state properties will be developed, and this will be spot checked and supplemented with data from GW/BSE calculations. The investigation will subsequently focus on the development of an improved exchange-correlation (XC) density functional beyond the current generalized gradient approximation (GGA) level of parameterization, with parameters fitted to the QMC database. This will allow the ground-state properties of focus systems to be more precisely predicted using DFT. These new developments will then be applied to investigate a chosen set of complex nanostructures that have great potential for opening new routes in designing materials with improved transport, electronic, and optical properties for PV and other optoelectronic usages: (1) Hybrid interfaces between materials with distinct electronic and optical properties, such as organic molecules (conjugated polymers, e.g. P3HT) and inorganic semiconducting materials (Si and ZnO). Complicated interface structures, including interface bonding configurations, compositional and geometrical blending patterns, interfacial defects, and various sizes and shapes of inorganic nanomaterials, will be considered for the purpose of understanding the working mechanisms of present organic/nano PV systems and designing optimum interface structures for fast charge separation and injection. (2) Complex-structured semiconducting nanomaterials that could induce charge separation without pn- or hetero-junctions. The new methodology will allow the PI to investigate the performance of realistic semiconducting nanomaterials of internal (impurities, defects, etc.) and external (uneven surface, mechanical twisting and bending, surface chemistry, etc.) complexities on optical absorption and charge transport against charge trapping and recombination. Of particular interest is whether such structural complexity in a single material could even be beneficial for PV usage, for example, charge separation through morphology control. Successful completion of the proposed DFT methodology would have a far-reaching impact on our ability to study and exploit the nature of electronic excitations in complex materials, advancing the design of next-generation electronic and optoelectronic devices in all facets of renewable energy conversion and storage, including photovoltaics, thermoelectricity, photochemistry, etc.

14 SOLAR ENERGY↗
DOE OSTI2019

IRNet: A General Purpose Deep Risidual Regression Framework for Materials Discovery

Materials discovery is crucial for making scientific advances in many domains. Collections of data from experiments and first-principle computations have spurred interest in applying machine learning methods to create predictive models capable of mapping from composition and crystal structures to materials properties. Generally, these are regression problems with the input being a 1D vector composed of numerical attributes representing the material composition and/or crystal structure. While neural networks consisting of fully connected layers have been applied to such problems, their performance often suffers from the vanishing gradient problem when network depth is increased. Hence, predictive modeling for such tasks has been mainly limited to traditional machine learning techniques such as Random Forest. In this paper, we study and propose design principles for building deep regression networks composed of fully connected layers with numerical vectors as input. We introduce a novel deep regression network with individual residual learning, IRNet, that places shortcut connections after each layer so that each layer learns the residual mapping between its output and input. We use the problem of learning properties of inorganic materials from numerical attributes derived from material composition and/or crystal structure to compare IRNet’s performance against that of other machine learning techniques. Using multiple datasets from the Open Quantum Materials Database (OQMD) and Materials Project for training and evaluation, we show that IRNet provides significantly better prediction performance than the state-of-the-art machine learning approaches currently used by domain scientists. We also show that IRNet’s use of individual residual learning leads to better convergence during the training phase than when shortcut connections are between multi-layer stacks while maintaining the same number of parameters.

deep learning↗
DOE OSTI2022

Comparisons Between Array Derived Dynamic Strain Rate (ADDS) and Fiber‐Optic Distributed Acoustic Sensing (DAS) Strain Rate

Abstract Distributed acoustic sensing (DAS) strain rate and particle velocity can be compared through approximate scaling with medium velocity. We instead performed a direct comparison between array derived dynamic strain (ADDS) rate and DAS strain rate for six frequency bands. The PoroTomo project at Brady's Hot Springs, Nevada, deployed a 240‐geophone 3C array co‐located with fiber‐optic DAS system and 8.7 km of buried cable. We selected subsets of the geophone array to create four smaller arrays and computed ADDS. The horizontal components of the ADDS were rotated into the direction of the fiber‐optic cable and then compared with the observed DAS strain rates. From three example regional earthquakes of local magnitudes 2.9, 4.1, and 4.3, the ADDS are found to be coherent with DAS for frequencies ≤1 Hz. For frequencies >1‐Hz, this correlation decays quickly. Small differences between linear and areal dynamic strains at 1‐Hz suggest poor signal‐to‐noise or localized strain that is perturbed by shallow heterogeneities compare to the average strain propagating across the geophone array. The implication is that around 1‐Hz, straight fiber DAS is measuring axial strain along the fiber and can provide good approximations to translational particle motions. However, above 1‐Hz, DAS becomes more sensitive to shallow velocity gradients that can be beneficial for geophysical imaging yet becomes a limitation for traditional seismic analysis methods depending on absolute amplitude and phase from translational particle motions.

58 GEOSCIENCES↗
DOE OSTI2015

Back to the future: using historical climate variation to project near-term shifts in habitat suitable for coast redwood

Studies that mimic the effect of climate change on terrestrial ecosystems often use climate projections from downscaled global climate models (GCMs). These simulations are generally too coarse to capture patterns of fine-scale climate variation, such as the sharp coastal energy and moisture gradients associated with wind-driven upwelling of cold water. Coastal upwelling may limit future increases in coastal temperatures, compromising GCMs' ability to provide realistic scenarios of future climate in these coastal ecosystems. Taking advantage of naturally occurring variability in the high-resolution historic climatic record, we developed multiple fine-scale scenarios of California climate that maintain coherent relationships between regional climate and coastal upwelling. We compared these scenarios against coarse resolution GCM projections at a regional scale to evaluate their temporal equivalency. We utilized these historically based scenarios to estimate potential suitable habitat for coast redwood (Sequoia sempervirens D. Don) under 'normal' combinations of temperature and precipitation, and under anomalous combinations representative of potential future climates. We found that a scenario of warmer temperature with historically normal precipitation is equivalent to climate projected by GCMs for California by 2020-2030 and that under these conditions, climatically suitable habitat for coast redwood significantly contracts at the southern end of its current range. Our results indicate that historical climate data provide a high-resolution alternative to downscaled GCM outputs for near-term ecological forecasts. This method may be particularly useful in other regions where local climate is strongly influenced by ocean-atmosphere dynamics that are not represented by coarse-scale GCMs.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2001

International Space Station Radiation Shielding Model Development

The projected radiation levels within the International Space Station (ISS) have been criticized by the Aerospace Safety Advisory Panel in their report to the NASA Administrator. Methods for optimal reconfiguration and augmentation of the ISS shielding are now being developed. The initial steps are to develop reconfigurable and realistic radiation shield models of the ISS modules, develop computational procedures for the highly anisotropic radiation environment, and implement parametric and organizational optimization procedures. The targets of the redesign process are the crew quarters where the astronauts sleep and determining the effects of ISS shadow shielding of an astronaut in a spacesuit. The ISS model as developed will be reconfigurable to follow the ISS. Swapping internal equipment rack assemblies via location mapping tables will be one option for shield optimization. Lightweight shield augmentation materials will be optimally fit to crew quarter areas using parametric optimization procedures to minimize the augmentation shield mass. The optimization process is being integrated into the Intelligence Synthesis Environment s (ISE s) immersive simulation facility at the Langley Research Center and will rely on High Performance Computing and Communication (HPCC) for rapid evaluation of shield parameter gradients.

Qualls, G. D.↗
NASA NTRS1971

Angular momentum desaturation for Skylab using gravity gradient torques

An angular momentum desaturation method for momentum exchange devices of orbiting spacecraft is described. The specific application of the method is to the Skylab which contains three double-gimbaled control moment gyros for precise attitude control and maneuvering. It is assumed that the attitude reference is inertially fixed and that two of the vehicle principal moments of inertia are much larger than the third. Gravity gradient torques and resultant angular momentum accumulation are developed for small deviations from the reference. The assumed moment-of-inertia distribution allows desaturation about all axes with only two attitude angles each for the two axes with large moments of inertia. The necessary desaturation maneuvers can be decoupled for a special set of orbital coordinates. All maneuvers are made during the night portion of the orbit, and the percentage utilized for desaturation is selectable. Expressions for the attitude angle commands are developed assuming infinite vehicle rates. The effect of finite rates introduces an efficiency into the desaturation. Expressions for this efficiency are developed and means for compensation are treated. Arbitrary misalignments between geometric vehicle axes and principal moment-of-inertia axes are permissible. An angle bias about the sun line minimizes the angular momentum accumulation about the sun line projection into the orbital plane. Adaptive desaturation maneuver limiting consistent with the available maneuver momentum is included.

Kennel, H. F.↗
DOE OSTI2018

Closures for Two-fluid Five-moment Equations for Magnetized Plasmas (Final Scientific Report)

In the proposed work, we will develop a complete set of closures for two-fluid five-moment equations for magnetized plasmas. Specifically, we propose to obtain explicit formulas for electrons and ions to express heat flow, viscosity, frictional force, and collisional heating in terms of density, and temperature, flow velocity. For this purpose, we solve a set of general moment equations for the relevant closure moments. An infinite hierarchy of general moment equations is completely equivalent to the original Landau-Fokker-Planck kinetic equation. In the small gyro-period expansion, the parallel components of the general moment equations agree with the moments of the drift-kinetic equation (DKE). The advantages of solving the moment equations include (i) an exact treatment of the collision operator, (ii) closure and transport relations with no flux-surface average, and (iii) a unified form of closures for general collisionality. Practically we solve a truncated set of moment equations and hence must verify the convergence of the solutions by increasing number of moments. The parallel integral (nonlocal) closures obtained by solving a set of linearized parallel moment equations are useful only when the nonlinear coupling of moments to the temperature and magnetic field gradient terms is negligible. When the nonlinear coupling terms are substantial, we propose a Fourier transform method to solve the parallel moment equations. By Fourier-expanding the magnetic field, we can convert the differential equations of parallel moments to linear algebraic equations for Fourier components of the moments. The convergence of the solutions can be confirmed by increasing the number of Fourier modes. Our preliminary study in axisymmetric circular magnetic geometry with a large aspect ratio shows that (i) the solutions for Fourier modes n=0 and n=1 converge when modes up to n=2 are included in the expansion, and (ii) lower collisionality requires more moments for convergent closures, as in the case of the integral closures. The proposed work will start from this simple magnetic geometry to find general parallel closures and then be generalized to explore a non-circular magnetic field with finite aspect-ratio. With parallel moments known, the perpendicular component of moment closures can be obtained by inverting the generalized cross product operator in the general moment equations. These efforts will be performed in conjunction with several multi-institutional projects, NIMROD, the Center for Extended MHD Modeling (CEMM), and the Plasma Science and Innovation (PSI) Center.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
DOE OSTI2017

Describing rainfall in northern Australia using multiple climate indices

Savanna landscapes are globally extensive and highly sensitive to climate change, yet the physical processes and climate phenomena which affect them remain poorly understood and therefore poorly represented in climate models. Both human populations and natural ecosystems are highly susceptible to precipitation variation in these regions due to the effects on water and food availability and atmosphere–biosphere energy fluxes. Here we quantify the relationship between climate phenomena and historical rainfall variability in Australian savannas and, in particular, how these relationships changed across a strong rainfall gradient, namely the North Australian Tropical Transect (NATT). Climate phenomena were described by 16 relevant climate indices and correlated against precipitation from 1900 to 2010 to determine the relative importance of each climate index on seasonal, annual and decadal timescales. Precipitation trends, climate index trends and wet season characteristics have also been investigated using linear statistical methods. In general, climate index–rainfall correlations were stronger in the north of the NATT where annual rainfall variability was lower and a high proportion of rainfall fell during the wet season. This is consistent with a decreased influence of the Indian–Australian monsoon from the north to the south. Seasonal variation was most strongly correlated with the Australian Monsoon Index, whereas yearly variability was related to a greater number of climate indices, predominately the Tasman Sea and Indonesian sea surface temperature indices (both of which experienced a linear increase over the duration of the study) and the El Niño–Southern Oscillation indices. These findings highlight the importance of understanding the climatic processes driving variability and, subsequently, the importance of understanding the relationships between rainfall and climatic phenomena in the Northern Territory in order to project future rainfall patterns in the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Measurement of Atmospheric Circular Polarization with POLARBEAR

At millimeter wavelengths, the atmospheric emission is circularly polarized owing to the Zeeman splitting of molecular oxygen by the Earth's magnetic field. We report a measurement of the signal in the 150 GHz band using 3 yr of observational data with the POLARBEAR project. Nonidealities of a continuously rotating half-wave plate (HWP) partially convert circularly polarized light to linearly polarized light. While POLARBEAR detectors are sensitive to linear polarization, this effect makes them sensitive to circular polarization. Although this was not the intended use, we utilized this conversion to measure circular polarization. We reconstruct the azimuthal gradient of the circular polarization signal and measure its dependency from the scanning direction and the detector bandpass. We compare the signal with a simulation based on atmospheric emission theory, the detector bandpass, and the HWP leakage spectrum model. We find the ratio of the observed azimuthal slope to the simulated slope is 0.92 ± 0.01(stat) ± 0.07(sys). This ratio corresponds to a brightness temperature of 3.8 mK at the effective band center of 121.8 GHz and bandwidth of 3.5 GHz estimated from representative detector bandpass and the spectrum of Zeeman emission. This result validates our understanding of the instrument and reinforces the feasibility of measuring the circular polarization using the imperfection of continuously rotating HWP. Continuously rotating HWP is popular in ongoing and future cosmic microwave background experiments to modulate the polarized signal. This work shows a method for signal extraction and leakage subtraction that can help measure circular polarization in such experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Stellar populations with MEGARA: The inner regions of NGC 7025

This paper aims to determine the capabilities of the MEGARA spectrograph at the Gran Telescopio Canarias (GTC), which is an optical integral-field unit, for studying stellar populations. We also aim to exploit its combination of high spectral (R ~ 6000, 12 000 and 20 000) and spatial (0.62") resolution within its 12".5 × 11"​​.3 field of view. We do this by analysing the commissioning data of the nearby S0a galaxy NGC 7025. We establish a systematic method through which we can determine the properties of the stellar populations in the observations made with MEGARA, more specifically, within the MEGADES legacy project. For this paper in particular, we determine the properties of the stellar populations of NGC 7025. We used MEGARA observations of galaxy NGC 7025 that were taken during the commissioning phase of the instrument. We applied different approaches to estimate the properties of the stellar populations with the highest possible certainty. In addition to the specific study of NGC 7025 and in the context of the MEGADES survey, we have carried out a number of tests to determine the expected errors (including potential biases) in these star formation history (SFH) derivations as a function of these parameters, namely spectral setup, signal-to-noise ratio, σ, and the SFH itself. All the studies we conduct (both full spectral fitting and absorption line indices) of the stellar populations of NGC 7025 indicate that the stars that form its bulge have supersolar metallicity and considerably old ages (~10 Gyr) in general. Using three different combinations of MEGARA spectral setups, we determined that the bulge of NGC 7025 has smild negative mass-weighted age gradient. For the more detailed SFH, our results indicate that in addition to a rather constant star formation at early epochs, a peak in the formation history of the stars in the bulge is also found 3.5–4.5 Gyr ago. This partly explains the mass-weighted age gradients we measured. The scenario presented in NGC 7025 is that of an isolated galaxy under secular evolution that about 3.5–4.5 Gyr ago likely experienced a minor merger (mass ratio 1/10) that induced an increase in star formation and also perturbed the morphology of its outer disc. In addition to these specific results for NGC 7025, we report different lessons learned for the ongoing exploitation of the MEGADES survey with the GTC, such as the need to obtain combined observations in the LR-B + LR-V setups and a signal-to-noise ratio of at least 20 per Å.

79 ASTRONOMY AND ASTROPHYSICS↗