Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M CODES”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A framework for measuring weak-lensing magnification using the Fundamental Plane

ABSTRACT Galaxy–galaxy lensing is an essential tool for probing dark matter haloes and constraining cosmological parameters. While galaxy–galaxy lensing measurements usually rely on shear, weak-lensing magnification contains additional constraining information. Using the Fundamental Plane (FP) of elliptical galaxies to anchor the size distribution of a background population is one method that has been proposed for performing a magnification measurement. We present a formalism for using the FP residuals of elliptical galaxies to jointly estimate the foreground mass and background redshift errors for a stacked lens scenario. The FP residuals include information about weak-lensing magnification κ, and therefore foreground mass, since to first order, non-zero κ affects galaxy size but not other FP properties. We also present a modular, extensible code that implements the formalism using emulated galaxy catalogues of a photometric galaxy survey. We find that combining FP information with observed number counts of the source galaxies constrains mass and photo-z error parameters significantly better than an estimator that includes number counts only. In particular, the constraint on the mass is 17.0 per cent if FP residuals are included, as opposed to 27.7 per cent when only number counts are included. The effective size noise for a foreground lens of mass $M_\mathrm{ H}=10^{14}\, \mathrm{M}_\odot$, with a conservative selection function in size and surface brightness applied to the source population, is σκ, eff = 0.250. We discuss the improvements to our FP model necessary to make this formalism a practical companion to shear analyses in weak-lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Schneider Springs Fire Study 2023 for Ecosystem Respiration Rates: Surface Water Chemistry and Hydrologic Sensor Data across the Yakima River Basin, Washington, USA (v2)

This dataset supports a broader study examining the drivers of spatial variability in wildfire impacts across the Yakima River Basin. Data provided within this dataset were generated from sample collection across 17 total sites (8 sites affected by a recent wildfire, 9 sites unaffected by a recent wildfire) within multiple rivers throughout the Yakima River Basin in Washington, USA from May-July 2023. Fire affected sites are defined as those affected by the 2021 Schneider Springs Fire, based on the drainage area of the streams being within the 2021 Schneider Springs Fire burn perimeter or not (Figure 1, below). The contents include surface water geochemistry data (dissolved organic carbon; total dissolved nitrogen; total suspended solids); short-term sonde data (specific conductivity; turbidity; pH; chlorophyll A; temperature); stream depth data; stream velocity; manual chamber open channel respiration data; sensor time-series data (oxygen; water pressure; barometric pressure); field metadata (including qualitative information on in stream and river corridor characteristics); and environmental context photos taken in the field. The dataset also includes a summary file of the sensor data and plots of the sensor data. Sensors were only recovered at 15 out of the 17 sites, and not all sensors were recovered at all 15 sites (see Methods section for more details), therefore all data does not exist at all sites. Data from a 2022 study at the same sites, as well as additional sites, can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1969566. The data package was originally published in November 2023. It was updated in June 2025 (v2; modified files). See the change history section in the readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one folder with field photos and one main data folder with two subfolders. The main data folder consists of (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) field protocol; (5) readme; (6) international generic sample number (IGSN) mapping file; and (7) stream depth and averages. The sensor data subfolder consists of (1) sensor installation methods summary; (2) stream velocity; and (3) six subfolders. The BarotrollAtm (barometric pressure; temperature), DepthHOBO (water pressure; temperature), MantaRiver (specific conductivity; turbidity; pH; chlorophyll A; temperature), EXO (specific conductivity; pH; temperature), miniDOT (dissolved oxygen; temperature), and miniDOTManualChamber (dissolved oxygen; temperature) contain time-series data, plots, and summary files. The sample data subfolder consists of (1) total suspended solids (TSS) data; (2) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (3) total dissolved nitrogen (TN) data and averages; and (4) methods codes. All files are .csv, .pdf, .jpg, .jpeg, or .mov.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Dynamic Aperture Evaluation for the Hadron Storage Ring in the Electron-Ion Collider

The Electron-Ion Collider (EIC) is aiming at a design luminosity of 1e34 cm⁻²s⁻¹. To maintain such a high luminosity, both beams in the EIC need an acceptable beam lifetime in the presence of the beam-beam interaction. For this purpose, we carried out weak-strong element-by-element particle tracking to evaluate the long-term dynamic aperture for the hadron ring lattice design. We improved our simulation code SimTrack to treat some new lattice design features, such as radially offset on-momentum orbits, coordinate transformations in the interaction region, etc. In this article, we will present the preliminary dynamic aperture calculation results with β^{*}- function scan, radial orbit shift, crossing angle collision, and magnetic field errors.

Luo, Y.↗
DOE OSTI2024

FFA@CEBAF beam transport error and tolerance simulation studies

The Continuous Electron Beam Accelerator Facility (CEBAF) is a 12 GeV recirculating electron accelerator at the Thomas Jefferson National Accelerator Facility (JLAB). Major upgrades to the accelerator are being investigated which include a new 650 MeV injection beamline and state-of-the-art fixed-field alternating (FFA) gradient recirculation arcs. The upgrade will extend the energy of the electron beam to over 20 GeV. In this paper, we provide an error and tolerance simulation study of the amended beam optics transport of the existing accelerator tuned for 22 GeV operation. The study is conducted with the particle tracking codes elegant and Bmad in two parts. In the first part, we treat each section of the accelerator (electromagnetic arcs and linacs) modularly with ideal conditions at the beginning. The second part is a pseudo start-to-end (S2E) simulation with accumulated errors propagating from one beamline to the next.

Accelerator Physics↗
DOE OSTI2020

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.3 is a program that obtains a steady-state flow and temperature solution for a nuclear reactor core, or for a single fuel assembly. It is based on an evolutionary sequence of codes originally used for plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over several decades. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional (Appendix XV). It begins with a solution across all coolant channels and fuel plates or tubes within a given fuel assembly, at the entrance to the assembly. The temperature solution is repeated for each axial node along the length of the fuel assembly. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset-of-nucleate boiling ratio(ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability is available (Appendix XII) to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A User’s Guide to the PLTEMP/ANL Code (V.4.3)

PLTEMP/ANL V4.3 is a program that obtains a steady-state flow and temperature solution for a nuclear reactor core, or for a single fuel assembly. It is based on an evolutionary sequence of codes originally used for plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over several decades. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional (Appendix XV). It should be noted that for a 3-dimensional solution the axial power shape is modeled without accounting for axial heat conduction. It begins with a solution across all coolant channels and fuel plates or tubes within a given fuel assembly, at the entrance to the assembly. The temperature solution is repeated for each axial node along the length of the fuel assembly. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset-of-nucleate boiling ratio(ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability is available (Appendix XII) to greatly reduce the reactor analyst’s time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Effects of the U.S. inflation reduction act on SMR economics

The U.S. Inflation Reduction Act (IRA) of 2022 provides a wide array of tax credits and other incentives for low-carbon energy. The technology-neutral clean generation production tax credit (PTC) (Section 45Y of the U.S. Internal Revenue Code) and the technology-neutral investment tax credit (ITC) (Section 48E) lower the net cost of new electricity generation projects with zero or negative greenhouse gas emission rates. We evaluate the impact of the IRA legislation—specifically the PTC and ITC—on the cost-competitiveness of small modular reactors (SMRs). We use the Argonne Low-carbon Energy Analysis Framework (A-LEAF) model to calculate the capacity factor of an SMR with a range of hypothetical variable operating and maintenance (O&M) costs in the Electric Reliability Council of Texas (ERCOT) electricity market. We selected ERCOT for market modeling because of its competitive structure, available data, and extensive use in prior literature. We use a discounted cash flow model to calculate the SMR’s net present value based on the market prices and capacity factors from A-LEAF, hypothetical ranges of capital and variable O&M costs, and other input parameters, with or without the IRA tax credits. We determine the SMR owner’s optimal choice of PTC or ITC for the hypothetical ranges of capital and variable O&M costs. We also evaluate potential shifts in the SMR owner’s optimal choice of PTC or ITC based on historical patterns of nuclear capital cost overruns in the United States. We also assess the sensitivity of our results to longer PTC period and electricity prices from the New England market, which tend to be higher than electricity prices in ERCOT. We find that even with the IRA tax credits, only SMRs with low capital and variable O&M costs would be economically feasible in the low-price ERCOT market scenario modeled. A longer PTC period and higher-price market such as New England, however, would significantly expand the economic feasibility of SMRs in the United States.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Overview on Shielding Analyses for the VENUS Instrument at SNS

VENUS, a world‐class versatile neutron imaging instrument, is under construction and is expected to be completed and ready to start commissioning in 2023. The range of cold to epithermal neutrons at SNS will give users of VENUS access to novel imaging methods, as well as to significantly improved existing methods. The instrument is being built on beam line 10 at Spallation Neutron Source (SNS) First Target Station (FTS) facing a decoupled poisoned hydrogen moderator. Instrument design, which includes optics, Front-End components and instrument enclosure components and shape started over 10 years ago and during this time had significant changes. All changes were supported by neutronics analyses to provide adequate shielding for both Front-End and instrument enclosure.Final optics design will provide the field of view (FOV) at the detector position (where the image is formed) to be as high as 0.20 by 0.20 m. Both the beam, which contains a large fraction of high-energy neutrons, and the desire for a large footprint on the detector are challenges for shielding design and budget, because the driving cost for the instrument is the beam line and enclosure shielding. For cost reason VENUS baseline design is optimized with respect to both the instrument cave footprint and its wall thickness, and Front-End shielding is tailored along the beam. Analyses are performed with the Monte Carlo particle transport code MCNPX version 2.7.0 to make choices on materials, thicknesses and configurations of the shielding. Numerous calculations are performed to meet space and coast constrain and to satisfy instrument physics needs and to comply with radiation protection requirements.

Gallmeier, Franz X.↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2020

SDSS-IV MaNGA: The link between bars and the early cessation of star formation in spiral galaxies

ABSTRACT Bars are common in low-redshift disc galaxies, and hence quantifying their influence on their host is of importance to the field of galaxy evolution. We determine the stellar populations and star formation histories of 245 barred galaxies from the Mapping Nearby Galaxies at APO (MaNGA) galaxy survey, and compare them to a mass- and morphology-matched comparison sample of unbarred galaxies. At fixed stellar mass and morphology, barred galaxies are optically redder than their unbarred counterparts. From stellar population analysis using the full spectral fitting code starlight, we attribute this difference to both older and more metal-rich stellar populations. Dust attenuation however, is lower in the barred sample. The star formation histories of barred galaxies peak earlier than their non-barred counterparts, and the galaxies build up their mass at earlier times. We can detect no significant differences in the local environment of barred and unbarred galaxies in this sample, but find that the H i gas mass fraction is significantly lower in high-mass ($\rm {M}_{\star } \gt 10^{10}~\rm {M}_{\odot }$) barred galaxies than their non-barred counterparts. We speculate on the mechanisms that have allowed barred galaxies to be older, more metal-rich and more gas-poor today, including the efficient redistribution of galactic fountain byproducts, and a runaway bar formation scenario in gas-poor discs. While it is not possible to fully determine the effect of the bar on galaxy quenching, we conclude that the presence of a bar and the early cessation of star formation within a galaxy are intimately linked.

Fraser-McKelvie, Amelia↗
DOE OSTI2026

RHOD Site - NOAA PSL Wind Retrievals WINDoe / Derived Data

This dataset contains daily NetCDF files with horizontal wind profiles retrieved with the WINDoe retrieval (Gebauer and Bell 2024) at Rhode Island (RHOD). WINDoe retrievals datasets are also available at Nantucket Island (NANT, nant.windoe.z01.c1) and Block Island (BLOC, bloc.windoe.z01.c1). WINDoe is an optimal estimation algorithm to retrieve wind profiles combining multiple instruments. The code is available in this github repository (https://github.com/OAR-atmospheric-observations/WINDoe/tree/main) and the retrieval is described by Gebauer and Bell (2024). WINDoe allows combining the individual datasets and outputs into one profile taking into account the information and uncertainties of each dataset. The use of WINDoe minimizes data gaps and maximizes data availability, compared to using wind profiles from only one of the instruments. The regular height grid eases comparisons to numerical weather prediction models. Code modifications have been made that include reading in WFIP3 specific instruments, averaging Doppler lidar radial velocities at various azimuth angles to avoid overfitting, and allowing the user to define a height grid by the user in the vipfile. The instruments used as input to the retrieval are a radar wind profiler (low- and high resolution mode) providing data in and above the boundary layer, a scanning Doppler lidar usually providing data throughout the boundary layer, a profiling lidar providing data from 50 to 200 m at BLOC and NANT, and from 10 to 280 m at Rhode Island, and a surface tower (4 m at NANT and RHOD and 10 m at BLOC). From the scanning lidars, we used radial velocity measurements at 60 deg elevation angle at six different azimuth angles with a resolution of approximately 30 m along the line of sight and the lowest range gate at approximately 70 m. The wind profiles are retrieved with WINDoe up to 3.74 km with 10 m vertical resolution. The profiles are retrieved every 15 min at BLOC and NANT and every 60 min at RHOD.

17 WIND ENERGY↗
DOE OSTI2026

BLOC Site - NOAA PSL Wind Retrievals WINDoe / Derived Data

This dataset contains daily netcdf files with horizontal wind profiles retrieved with the WINDoe retrieval (Gebauer and Bell 2024) at Block Island (BLOC). WINDoe retrievals datasets are also available at Nantucket Island (NANT, nant.windoe.z01.c1) and Rhode Island (RHOD, rhod.windoe.z01.c1). WINDoe is an optimal estimation algorithm to retrieve wind profiles combining multiple instruments. The code is available in this github repository (https://github.com/OAR-atmospheric-observations/WINDoe/tree/main), and the retrieval is described by Gebauer and Bell (2024). WINDoe allows combining the individual datasets and outputs into one profile taking into account the information and uncertainties of each dataset. The use of WINDoe minimizes data gaps and maximizes data availability, compared to using wind profiles from only one of the instruments. The regular height grid eases comparisons to numerical weather prediction models. Code modifications have been made that include reading in WFIP3 specific instruments, averaging Doppler lidar radial velocities at various azimuth angles to avoid overfitting, and allowing the user to define a height grid by the user in the vipfile. The instruments used as input to the retrieval are a radar wind profiler (low- and high resolution mode) providing data in and above the boundary layer, a scanning Doppler lidar usually providing data throughout the boundary layer, a profiling lidar providing data from 50 to 200 m at BLOC and NANT, and from 10 to 280 m at Rhode Island, and a surface tower (4 m at NANT and RHOD and 10 m at BLOC). From the scanning lidars, we used radial velocity measurements at 60 deg elevation angle at six different azimuth angles with a resolution of approximately 30 m along the line of sight and the lowest range gate at approximately 70 m. The wind profiles are retrieved with WINDoe up to 3.74 km with 10 m vertical resolution. The profiles are retrieved every 15 min at BLOC and NANT and every 60 min at RHOD.

17 WIND ENERGY↗
DOE OSTI2026

NANT Site - NOAA PSL Wind Retrievals WINDoe / Derived Data

This dataset contains daily NetCDF files with horizontal wind profiles retrieved with the WINDoe retrieval (Gebauer and Bell 2024) at Nantucket Island (NANT). WINDoe retrievals datasets are also available at Block Island (BLOC, bloc.windoe.z01.c1) and Rhode Island (RHOD, rhod.windoe.z01.c1). WINDoe is an optimal estimation algorithm to retrieve wind profiles combining multiple instruments. The code is available in this github repository (https://github.com/OAR-atmospheric-observations/WINDoe/tree/main), and the retrieval is described by Gebauer and Bell (2024). WINDoe allows combining the individual datasets and outputs into one profile taking into account the information and uncertainties of each dataset. The use of WINDoe minimizes data gaps and maximizes data availability, compared to using wind profiles from only one of the instruments. The regular height grid eases comparisons to numerical weather prediction models. Code modifications have been made that include reading in WFIP3 specific instruments, averaging Doppler lidar radial velocities at various azimuth angles to avoid overfitting, and allowing the user to define a height grid by the user in the vipfile. The instruments used as input to the retrieval are a radar wind profiler (low- and high resolution mode) providing data in and above the boundary layer, a scanning Doppler lidar usually providing data throughout the boundary layer, a profiling lidar providing data from 50 to 200 m at BLOC and NANT, and from 10 to 280 m at Rhode Island, and a surface tower (4 m at NANT and RHOD and 10 m at BLOC). From the scanning lidars, we used radial velocity measurements at 60 deg elevation angle at six different azimuth angles with a resolution of approximately 30 m along the line of sight and the lowest range gate at approximately 70 m. The wind profiles are retrieved with WINDoe up to 3.74 km with 10 m vertical resolution. The profiles are retrieved every 15 min at BLOC and NANT and every 60 min at RHOD.

17 WIND ENERGY↗
DOE OSTI2020

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in Its Infancy by the Zwicky Transient Facility

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. Here, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of $R={744.0}_{-128.0}^{+183.0}\,{R}_{\odot }$ and an ejected mass of ${M}_{\mathrm{ej}}={9.3}_{-5.8}^{+0.4}\,{M}_{\odot }$. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Simulation studies of the interaction of laser radiation with additively manufactured foams

The interaction of laser radiation with foams of various porosities and low densities has been the subject of several numerical and experimental studies. In all cases, the modeling of low-Z under-dense foams as uniform gases of equivalent average density using standard radiation-hydrodynamics codes has resulted in heat-front velocities that are considerably faster than those observed experimentally. It has been theoretically conjectured that this difference may be attributed to the breakdown of the foam’s morphology, leading to a dynamics of filament expansion where the ion and electron energy partitions are significantly different from those calculated using the uniform gas model. We found that 3D computer simulations employing a disconnected representation of the foam’s microstructure which allowed for the dynamics of foam element heating, expansion, and stagnation largely supported the theoretical picture. Simulations using this model for laser experiments on under-dense 2 mg cc -1 SiO 2 aerogel foams reproduced the experimental data fairly well. We used the validated model in simulations of low-density structured foam-like materials (produced via additive manufacturing) with a variety of morphologies. We found that the log-pile configurations were consistent with the analytical propagation model of Gus’kov et al. Further validation of the model was obtained by simulating experiments performed at the Jupiter Laser Facility using the log-pile and octet-truss foam morphologies. Simulations of the foam-laser interaction using a wave propagation code showed that the microstructure was able to enhance stimulated Brillouin scattering (SBS) by concentrating the light energy into density holes. In turn, this promotes laser filamentation, reducing SBS and bringing the predicted values closer to the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗