Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 271 records · Page 15
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 OSTI2012

Protein subcellular location pattern classification in cellular images using latent discriminative models

Motivation: Knowledge of the subcellular location of a protein is crucial for understanding its functions. The subcellular pattern of a protein is typically represented as the set of cellular components in which it is located, and an important task is to determine this set from microscope images. In this article, we address this classification problem using confocal immunofluorescence images from the Human Protein Atlas (HPA) project. The HPA contains images of cells stained for many proteins; each is also stained for three reference components, but there are many other components that are invisible. Given one such cell, the task is to classify the pattern type of the stained protein. We first randomly select local image regions within the cells, and then extract various carefully designed features from these regions. This region-based approach enables us to explicitly study the relationship between proteins and different cell components, as well as the interactions between these components. To achieve these two goals, we propose two discriminative models that extend logistic regression with structured latent variables. The first model allows the same protein pattern class to be expressed differently according to the underlying components in different regions. The second model further captures the spatial dependencies between the components within the same cell so that we can better infer these components. To learn these models, we propose a fast approximate algorithm for inference, and then use gradient-based methods to maximize the data likelihood. Results: In the experiments, we show that the proposed models help improve the classification accuracies on synthetic data and real cellular images. The best overall accuracy we report in this article for classifying 942 proteins into 13 classes of patterns is about 84.6%, which to our knowledge is the best so far. In addition, the dependencies learned are consistent with prior knowledge of cell organization.

59 BASIC BIOLOGICAL SCIENCES↗
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↗
DOE OSTI2012

Two self-referencing methods for the measurement of beam spot position

Purpose: Two quantitative methods of measuring electron beam spot position with respect to the collimator axis of rotation (CAOR) are described. Methods: Method 1 uses a cylindrical ion chamber (IC) mounted on a jig corotational with the collimator making the relationship among the chamber, jaws, and CAOR fixed and independent of collimator angle. A jaw parallel to the IC axis is set to zero and the IC position adjusted so that the IC signal is approximately 50% of the open field value, providing a large dose gradient in the region of the IC. The cGy/MU value is measured as a function of collimator rotation, e.g., every 30 Degree-Sign . If the beam spot does not lie on the CAOR, the signal from the ion chamber will vary with collimator rotation. Based on a measured spatial sensitivity, the distance of the beam spot from the CAOR can be calculated from the IC signal variation with rotation. The 2nd method is image based. Two stainless steel rods, 3 mm in diameter, are mounted to a jig attached to the Linac collimator. The rods, offset from the CAOR, lay in different planes normal to the CAOR, one at 158 cm SSD and the other at 70 cm SSD. As the collimator rotates the rods move tangent along an envelope circle, the centers of which are on the CAOR in their respective planes. Three images, each at a different collimator rotation, containing the shadows of both rods, are acquired on the Linac EPID. At each angle the shadow of the rods on the EPID defines lines tangent to the projection of the envelope circles. From these the authors determine the projected centers of the two circles at different heights. From the distance of these two points using the two heights and the source to EPID distance, the authors calculate the distance of the beam spot from the CAOR. Measurements with all two techniques were performed on an Elekta Linac. Measurements were performed with the beam spot in nominal clinical position and in a deliberately offset position. Measurements were also performed using the Flexmap image registration/ball-bearing test. Results: Within their uncertainties, both methods report the same beam spot displacement. In clinical use, a total of 203 monthly beam spot measurements on 14 different beams showed an average displacement of 0.11 mm ({sigma}= 0.07 mm) in-plane and 0.10 mm ({sigma}= 0.07 mm) cross-plane with maximum displacement of 0.37 mm in-plane and 0.34 mm cross-plane. Conclusions: The methods described provide a quantitative measure of beam spot position, are easy to use, and provide another tool for Linac setup and quality assurance. Fundamental to the techniques is that they are self-referencing-i.e., they do not require the user to independently define the CAOR.

43 PARTICLE ACCELERATORS↗
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↗
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 OSTI2011

Transverse Emittance Reduction with Tapered Foil

The idea of reducing transverse emittance with tapered energy-loss foil is proposed by J.M. Peterson in 1980s and recently by B. Carlsten. In this paper, we present the physical model of tapered energy-loss foil and analyze the emittance reduction using the concept of eigen emittance. The study shows that, to reduce transverse emittance, one should collimate at least 4% of particles which has either much low energy or large transverse divergence. The multiple coulomb scattering is not trivial, leading to a limited emittance reduction ratio. Small transverse emittances are of essential importance for the accelerator facilities generating free electron lasers, especially in hard X-ray region. The idea of reducing transverse emittance with tapered energy-loss foil is recently proposed by B. Carlsten [1], and can be traced back to J.M. Peterson's work in 1980s [2]. Peterson illustrated that a transverse energy gradient can be produced with a tapered energy-loss foil which in turn leads to transverse emittance reduction, and also analyzed the emittance growth from the associated multiple coulomb scattering. However, what Peterson proposed was rather a conceptual than a practical design. In this paper, we build a more complete physical model of the tapered foil based on Ref. [2], including the analysis of the transverse emittance reduction using the concept of eigen emittance and confirming the results by various numerical simulations. The eigen emittance equals to the projected emittance when there is no cross correlation in beam's second order moments matrix [3]. To calculate the eigen emittances, it requires only to know the beam distribution at the foil exit. Thus, the analysis of emittance reduction and the optics design of the subsequent beam line section can be separated. In addition, we can combine the effects of multiple coulomb scattering and transverse energy gradient together in the beam matrix and analyze their net effect. We find that,when applied to an electron linac or electron beam line, the energy spread increase and angular growth due to multiple scattering are not trivial; as a result, the transverse emittance can only be reduced with a limited ratio, e.g. down to about 65% the original value. The contents of this paper are arranged as follows. In Sec. II, we build the physical model of the tapered foil, derive the transverse eigen emittance and discuss the emittance reduction criteria. In Sec. III, we implement numerical simulations to verify the physical model; and in Sec. IV, we present numerical experiments and subsequent beam line to remove the transverse energy gradient to demonstrate the applicability of such method. Conclusions are given in the last section.

43 PARTICLE ACCELERATORS↗
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↗
DOE OSTI2017

Numerical and Experimental Investigation of Turbulent Transport Control via Shaping of Radial Plasma Flow Profiles

Turbulence, and turbulence-driven transport are ubiquitous in magnetically confined plasmas, where there is an intimate relationship between turbulence, transport, instability driving mechanisms (such as gradients), plasma flows, and flow shear. Though many of the detailed physics of the interrelationship between turbulence, transport, drive mechanisms, and flow remain unclear, there have been many demonstrations that transport and/or turbulence can be suppressed or reduced via manipulations of plasma flow profiles. This is well known in magnetic fusion plasmas [e.g., high confinement mode (H-mode) and internal transport barriers (ITB’s)], and has also been demonstrated in laboratory plasmas. However, it may be that the levels of particle transport obtained in such cases [e.g. H-mode, ITB’s] are actually lower than is desirable for a practical fusion device. Ideally, one would be able to actively feedback control the turbulent transport, via manipulation of the flow profiles. The purpose of this research was to investigate the feasibility of using both advanced model-based control algorithms, as well as non-model-based algorithms, to control cross-field turbulence-driven particle transport through appropriate manipulation of radial plasma flow profiles. The University of New Mexico was responsible for the experimental portion of the project, while our collaborators at the University of Montana provided plasma transport modeling, and collaborators at Lehigh University developed and explored control methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗