Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ring model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3

Robust inference of the Galactic Centre gamma-ray excess spatial properties

ABSTRACT The gamma-ray Fermi-LAT Galactic Centre excess (GCE) has puzzled scientists for over 15 yr. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology – potentially arising from the annihilation of dark matter (DM) particles – with a boxy morphology – expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VISTA Variables in the Via Lactea survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a galprop-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology for the GCE, regardless of the mask applied to the Galactic plane.

Astronomy & Astrophysics↗

Rapid Tuning of Synchrotron Surrogate Model at the Recycler Ring

The 8 GeV proton-storage Recycler Ring (RR) is essential for reaching megawatt beam intensity goals for the DUNE neutrino beam at Fermilab. Custom shims on each RR permanent magnet were designed to cancel manufacturing defects and bring magnetic fields to the design values. Remaining imperfections cause the observed tune variation vs energy to deviate from what is calculated using the design fields. Using the POUNDERS (“Practical Optimization Using No Derivatives for sums of Squares”) optimization method with Synergia in the loop, we demonstrate rapid convergence to a set of additive, higher-order multipole moments of these magnetic shims which reproduce that observed variation, and show that the convergence advantage grows with the parameter-space dimensionality.

43 PARTICLE ACCELERATORS↗

Multifrequency Models of Black Hole Photon Rings from Low-luminosity Accretion Disks

Images of black holes encode both astrophysical and gravitational properties. Detecting highly lensed features in images can differentiate between these two effects. We present an accretion disk emission model coupled to the Adaptive Analytical Ray Tracing (AART) code that allows a fast parameter space exploration of black hole photon ring images produced from synchrotron emission from 10 to 670 GHz. As an application, we systematically study several disk models and compute their total flux density, average radii, and optical depth. The model parameters are chosen around fiducial values calibrated to general relativistic magnetohydrodynamic (GRMHD) simulations and observations of M87*. For the parameter space studied, we characterize the transition between optically thin and thick regimes and the frequency at which the first photon ring is observable. Our results highlight the need for careful definitions of photon ring radius in the image domain, as in certain models the highly lensed photon ring is dimmer than the direct emission at certain angles. We find that at low frequencies the ring radii are set by the electron temperature, while at higher frequencies the magnetic field strength plays a more significant role, demonstrating how multifrequency analysis can also be used to infer plasma parameters. Lastly, we show how our implementation can qualitatively reproduce multifrequency black hole images from GRMHD simulations when adding time variability to our disk model through Gaussian random fields. This approach provides a new method for simulating observations from the Event Horizon Telescope and the proposed Black Hole Explorer space mission.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerator Physics at NSLS-II: Research Accomplishments in 2022

NSLS-II accelerator physicists are working on scientific support of the NSLS-II operations and start-up of the accelerator systems including Storage Ring, Injector, beam diagnostics, and top-off injection; improvement of beam stability and orbit feedbacks; commissioning of a superconducting wiggler for the HEX beamline; development of a new operation mode for time-resolved experiments. The routine lattice characterization is continued to keep tracking the machine, the storage ring lattice model including strong coupling is developed, and seasonal variations of the ring circumference are studied. The collaboration with the EIC project is ongoing. In the frameworks of three LDRD projects, the physicists are working on the lattice development, studies of collective effects, and specifications of alignment and stability for the NSLS-II accelerator upgrade; novel methods of lattice design and optimization; studies of free-electron laser options; ultrafast electron diffraction and microscopy. Research accomplishments achieved in 2022, are summarized in this report.

43 PARTICLE ACCELERATORS↗

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Data for "Depth of nutrient uptake by deep-rooted plants is regulated by water availability"

The data set consists of strontium (Sr) isotope ratios (87Sr/86Sr), water isotopes, soil cation concentrations, soil water potential sensor data, and results of 87Sr/86Sr mixing model. The plant canopy size files include the dataset of canopy dimension of sagebrush, lupine, and sunflower. The soil and plant ICPMS (Inductively Coupled Plasma Mass Spectrometry) data file includes both of 87Sr/86Sr, and cation concentration dataset from soil exchangeable pool, apatite pool, silicate extract, atmospheric rain deposition, and plant leaf and stem tissues. The plant dendrochronology file includes the dendrochronogical ring width of several sagebrush, and dendrochemical sample data includes the 87Sr/86Sr for each separated growth ring. The modeling result gives the proportion of nutrient sources of each plants (based on their 87Sr/86Sr in leaf tissues and growth rings) from atmospheric deposition and mineral weathering. Soil water potential data includes continuous collection of soil water potential dataset at 2 depths (30 cm and 60 cm, from Nov 24 - Jun 25) of the sampling site. All the samples were collected from 2 sampling campaign June and July 2023, and rain water is a separate sampling from Aug - Sept 2023, at north-facing hillslope near pumphouse site. The data showed that the depth of cation nutrient acquisition is thus tightly coupled with, and likely determined by, water availability in soil, saprolite and bedrock. The enhanced uptake of cations and water from regions of mineral weathering could confer plant and ecosystem resilience during low water years and may impact the rate of bedrock weathering and watershed chemistry during drought. This dataset includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type; a location metadata file (locations.csv); and a samples metadata file (samples.csv). All files are provided as comma-separated values (CSV) files (.csv). This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Compare low-gain 1D and 3D models toward a storage ring-based X-ray free-electron-laser oscillator

Fourth-generation diffraction-limited synchrotron light sources (DL-SLS) present a compelling platform for advanced light source applications due to their ultra-low emittance (around tens of picometers), sub-micrometer beam stability, and high repetition rates (greater than 1 MHz). Despite these advantages, fundamental limitations, such as long bunch durations (tens of picoseconds), large relative energy spread (~10⁻³), and low peak current (<300 A), have hindered their effectiveness in free-electron laser (FEL) applications. Using a low-gain three-dimensional (3D) model based on NSLS-II upgrade lattice parameters, we demonstrate that a 6% single-pass gain is achievable with an optimized transverse gradient undulator (TGU). Concurrently, recent advances in X-ray optics—specifically, a four-mirror bow-tie resonator utilizing diamond mirrors and beryllium compound refractive lenses (CRLs)—enable total round-trip losses below 4% and an output coupling efficiency near 1%, rendering a Storage Ring based X-ray FEL oscillator (XFELO) a realistic possibility. In the 3D framework, the TGU effectively mitigates the detrimental effects of large energy spread by spatially dispersing the beam, leading to enhanced gain relative to conventional one-dimensional (1D) models. We establish both analytical and numerical connections between the 1D and 3D descriptions, leveraging GENESIS simulations to map equivalent gain conditions. While 1D resonance is defined in the energy domain, the 3D TGU system shifts the resonance dependence to transverse gradient and dispersion parameters. Once optimized, this framework allows predictive modeling using 1D formulas, which can then be validated through full 3D simulations.

43 PARTICLE ACCELERATORS↗

The Effects of Field Line Curvature (FLC) Scattering on Ring Current Dynamics and Isotropic Boundary

Abstract In the ring current dynamics, various loss mechanisms contribute to the ring current decay, including losses to the upper atmosphere through particle precipitation. This study implements the field line curvature (FLC) scattering mechanism in a kinetic ring current model and investigates its role in precipitating ions into the ionosphere during the 17 March 2013 storm. Simulation results indicate that (1) the FLC scattering process exerts on energetic ions on the nightside where the magnetospheric configuration is more stretching. It is more effective on heavy ions (e.g., O + ). These ion losses thereafter lead to a faster recovery of the ring current. (2) The FLC‐associated ion precipitation mainly occurs in the outer region (L > 5 for protons and L > 4.5 for oxygen ions) on the nightside. The O + precipitation takes places in a wider region than protons although its intensity is much lower. Comparisons with POES observations suggest that more proton precipitation is needed in the inner region. This is probably caused by the less stretched configuration in the simulation that prevents more precipitation. It may also imply that other loss process is required in the model such as wave‐particle interactions. (3) The storm time precipitating proton flux of tens of keV due to the FLC scattering sometimes becomes comparable to that of electrons on the nightside, although electrons usually dominate the ionospheric energy deposition from the midnight eastward toward the dayside. (4) The FLC scattering process seems to be capable of explaining the formation of isotropic boundary in the ionosphere during the investigated event.

79 ASTRONOMY AND ASTROPHYSICS↗

Advancements in NEAMS Tool Capabilities for Multiphysics Simulation of Fast Reactor Core Bowing and Identification of Validation Test Data

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves Multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Last year, thermal fluids and mechanics were coupled on a multi-assembly benchmark problem based on ABR-1000 design. This year, the reactor physics code Griffin was assessed for readiness of core bowing calculations. Preliminary integration of Griffin’s ring-heterogeneous model with thermal fluids and thermal mechanics solvers was performed. Specifically, thermal-mechanics and reactor physics were coupled for single- and multi-assembly problems, and reactor physics and subchannel methods were coupled for a single assembly model. Finally, the workflow of all three physics was preliminarily demonstrated on a single assembly model. Caveats and future development needed have been identified. To supplement the multiphysics demonstration, verification and assessment efforts of thermos-mechanical capabilities for modeling thermo-mechanical core bowing behavior were continued by analyzing IAEA Verification Problem 5 which includes radiation swelling and creep. Additionally, a small core reactor physics benchmark defined by Japan Atomic Energy Agency (JAEA) was performed to assess neutronics models for estimating reactivity feedback. Finally, Fast Flux Test Facility (FFTF) validation test data for core bowing phenomena has been identified and summarized, with a recommended path forward for validation once this capability is mature.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗

On the Ion Precipitation due to Field Line Curvature (FLC) and EMIC Wave Scattering and Their Subsequent Impact on Ionospheric Electrodynamics

Both field line curvature (FLC) and electromagnetic ion cyclotron (EMIC) wave scattering are believed to be associated with energetic ion precipitation, but their relative contributions to ionospheric ion precipitation and subsequent effects are still unclear. Here, in this study, by using a kinetic ring current model, we investigate their impact on the ionosphere from two aspects: the global distribution of ion precipitation and resulting ionospheric conductance. Our results show that the intensity and coverage of ion precipitation due to EMIC waves are larger than that due to FLC scattering, while the latter mostly contributes to ion precipitation in outer regions (L > 4–5). We then estimate the conductance with empirical models using the simulated ion precipitation energy flux. When the EMIC wave associated proton precipitation is included, the conductance is significantly enhanced in the dusk-to-midnight sector and has a wide magnetic latitude (MLAT) range from around 52° to 62°, considerably altering the electric potential in the dusk sector, which further influences particle dynamics in the magnetosphere. On the contrary, the proton precipitation caused by FLC scattering only occurs at higher latitudes above MLAT = 60° and the corresponding conductance is slightly enhanced at midnight, with negligible influence on the convective electric potential. Although electron precipitation-associated conductance is predominant globally, the results show that proton precipitation can also play an important role in ionospheric electrodynamics, especially when EMIC wave-associated precipitating proton energy flux exceeds that of electrons in the dusk sector, a region where many subauroral coupling processes occur.

79 ASTRONOMY AND ASTROPHYSICS↗

Multimaterial Aerosol Jet Printed Magnetic Nanocomposites for Microwave Circuits

Here this article demonstrates the process and characterization of additively manufactured polymer-matrix magnetic nanocomposites. These nanocomposites are composed of nickel-zinc ferrite (Ni 0.5 Zn 0.5 Fe 2 O 4 ) nanoparticles in a polyimide matrix. These materials are printed using multimaterial aerosol jet printing where aerosols of nickel-zinc ferrite and polyimide inks are mixed in place to fabricate structures with nonunity relative permeability. To characterize these materials, we use energy-dispersive X-ray spectroscopy to estimate the content of ferrite in the composite, fitted simulation models, and ring resonator measurements to estimate the permeability. In fully aerosol jet printed 1-mm 2 inductors, we demonstrate an increase in inductance density of 40% in a 4.5-turn microstrip spiral inductor using this printed composite. These fully printed inductors demonstrate an inductance density of 4.2-nH/mm 2 . EDS analysis shows films with 18 wt.% Ni 0.5 Zn 0.5 Fe 2 O 4 . Finally, we show microstrip transmission lines on printed composite compared to the same on polyimide subjected to a dc magnetic bias, displaying a magnetic resonance behavior.

42 ENGINEERING↗

Analytical estimates of beam intensity limitations in the EIC Beam Accumulator Ring

This report summarizes analytical estimates of beam intensity limitations in the Beam Accumulator Ring (BAR), a part of the Electron-Ion Collider injector complex under development at Brookhaven National Laboratory. The BAR is designed to accumulate a polarized electron beam with a single-bunch charge up to 28 nC for subsequent acceleration in the Rapid Cycling Synchrotron and then injection into the Electron Storage Ring. Analytical models of longitudinal and transverse collective effects are developed using a broadband impedance model calibrated with the beam-based measurement results from NSLS VUV. The study evaluates the impact of resistive-wall and geometric impedances, assessing microwave and transverse mode coupling instabilities, bunch lengthening, and beam-induced heating.

43 PARTICLE ACCELERATORS↗

The Variability of the Black Hole Image in M87 at the Dynamical Timescale

The black hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5–61 days) is comparable to the 6 day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure-phase measurements on all six linearly independent nontrivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of ∼3°–5°. The only triangles that exhibit substantially higher variability (∼90°–180°) are the ones with baselines that cross the visibility amplitude minima on the u–v plane, as expected from theoretical modeling. We used two sets of general relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

Surrogate Modelling of 3rd Integer Resonant Extraction at Fermilab Delivery Ring

We present an ongoing work in which a surrogate model is being developed to reproduce the response dynamics of the third-integer resonant extraction process in the Delivery Ring (DR) at Fermilab. This effort is in pursuit of smoothly extracting circulating beam to the Mu2e Experiment s production target, wherein the goal is to extract a uniform slice of the circulating $1e12$ protons in the DR over 25,000 turns (43~ms). The DR contains 3 harmonic sextupoles which excite a third-integer resonance as well as three fast, tune-ramping quadrupole magnets which drive the horizontal tune towards the $29/3$ resonance. In our initial work the surrogate model trains on a semi-analytical simulation provided in the same format as live data. Using Reinforcement Learning (and other potential ML methods), the trained surrogate acts as the environment in which a simple ML control agent could learn to dynamically adjust the quadrupole ramp at 430 break points within the 43 microsecond spill window. The control agent will be hosted on a dedicated Arria 10 FPGA, introducing its own requirements on control agent architecture. In this work we report the accuracy and fidelity of surrogate models in comparison to the response dynamics of the physics simulator.

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

Constraining long‐term model predictions for woody growth using tropical tree rings

The strength and persistence of the tropical carbon sink hinges on the long-term responses of woody growth to climatic variations and increasing CO 2 . However, the sensitivity of tropical woody growth to these environmental changes is poorly understood, leading to large uncertainties in growth predictions. Here, we used tree ring records from a Southeast Asian tropical forest to constrain ED2.2-hydro, a terrestrial biosphere model with explicit vegetation demography. Specifically, we assessed individual-level woody growth responses to historical climate variability and increases in atmospheric CO 2 (C a ). When forced with historical C a , ED2.2-hydro reproduced the magnitude of increases in intercellular CO 2 concentration (a major determinant of photosynthesis) estimated from tree ring carbon isotope records. In contrast, simulated growth trends were considerably larger than those obtained from tree rings, suggesting that woody biomass production efficiency (WBPE = woody biomass production:gross primary productivity) was overestimated by the model. The estimated WBPE decline under increasing C a based on model-data discrepancy was comparable to or stronger than (depending on tree species and size) the observed WBPE changes from a multi-year mature-forest CO 2 fertilization experiment. In addition, we found that ED2.2-hydro generally overestimated climatic sensitivity of woody growth, especially for late-successional plant functional types. The model-data discrepancy in growth sensitivity to climate was likely caused by underestimating WBPE in hot and dry years due to commonly used model assumptions on carbon use efficiency and allocation. To our knowledge, this is the first study to constrain model predictions of individual tree-level growth sensitivity to C a and climate against tropical tree-ring data. Our results suggest that improving model processes related to WBPE is crucial to obtain better predictions of tropical forest responses to droughts and increasing C a . Finally, more accurate parameterization of WBPE will likely reduce the stimulation of woody growth by C a rise predicted by biosphere models.

54 ENVIRONMENTAL SCIENCES↗