SEARCH · Engineering Papers

Results for “coding standards”

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

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Direct measurement of building transient and residual drift using an optical sensor system

Interstory drift (ID) is a key response parameter for buildings subjected to lateral loads and is used to define performance-based limit states, allowable deformations and damage states in a number of seismic design codes and standards. Here, an ability to rapidly and accurately measure both transient and residual ID during an earthquake would provide important observables for understanding the seismic demands and post-earthquake condition of a building. Accurate retrieval of ID from accelerometer-based instrumentation systems can be very challenging, if not impossible, as a result of instrumentation limitations and the post-processing associated with strong motion accelerometer data. This is particularly true for the case in which residual drifts occur during inelastic building response. In the study presented herein, a newly developed optical sensor system, designed specifically for directly measuring both transient and residual ID, was experimentally evaluated through shake table testing and computational simulations. The ability of the sensor to accurately measure ID is demonstrated and key operational considerations for sensor system deployment are examined through a model-based investigation.

42 ENGINEERING↗

DOE OSTI2023

A review of challenges, barriers, and opportunities for large-scale deployment of cool surfaces

Major urban centers are warming due to a combination of global and local phenomena. City governments are increasingly adopting strategies to mitigate the causes and impacts of extreme heat on their populations. Among these strategies are high solar-reflectance (cool) surfaces installed on building roofs and walls. Use of cool surfaces is a cost-effective and simple strategy that replaces conventional darker surfaces with surfaces that have a high reflectance to shortwave (solar) energy. This report reviews the recent history of cool-surface deployment efforts. This includes peer-reviewed literature, conference proceedings, and grey literature to identify challenges and barriers to wide-scale deployment of cool surfaces. We have also researched heat action plans and programs from cities and different codes and standards, as well as available incentive and rebate programs. The review identifies challenges, barriers, and opportunities associated with large-scale deployment of cool surfaces and categorize them broadly as being related to product development & performance or policies & mandates. It provides a foundation upon which we intend to build a roadmap for rapidly accelerating future deployments of cool surfaces. Finally, this roadmap will address identified challenges and incorporate lessons learned from historical efforts to generate a practical and actionable plan.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2017

Ignition and flame characteristics of cryogenic hydrogen releases

In this work, under-expanded cryogenic hydrogen jets were investigated experimentally for their ignition and flame characteristics. The test facility described herein, was designed and constructed to release hydrogen at a constant temperature and pressure, to study the dispersion and thermo-physical properties of cryogenic hydrogen releases and flames. In this study, a non-intrusive laser spark focused on the jet axis was used to measure the maximum ignition distance. The radiative power emitted by the corresponding jet flames was also measured for a range of release scenarios from 37 K to 295 K, 2–6 bar abs through nozzles with diameters from 0.75 to 1.25 mm. The maximum ignition distance scales linearly with the effective jet diameter (which scales as the square root of the stagnant fluid density). A 1-dimensional (stream-wise) cryogenic hydrogen release model developed previously at Sandia National Laboratories (although this model is not yet validated for cryogenic hydrogen) was exercised to predict that the mean mole fraction at the maximum ignition distance is approximately 0.14, and is not dependent on the release conditions. The flame length and width were extracted from visible and infra-red flame images for several test cases. The flame length and width both scale as the square root of jet exit Reynolds number, as reported in the literature for flames from atmospheric temperature hydrogen. As shown in previous studies for ignited atmospheric temperature hydrogen, the radiative power from the jet flames of cold hydrogen scales as a logarithmic function of the global flame residence time. The radiative heat flux from jet flames of cold hydrogen is higher than the jet flames of atmospheric temperature hydrogen, for a given mass flow rate, due to the lower choked flow velocity of low-temperature hydrogen. Lastly, this study provides critical information with regard to the development of models to inform the safety codes and standards of hydrogen infrastructure.

08 HYDROGEN↗

DOE OSTI2019

Modeling of underexpanded hydrogen jets through square and rectangular slot nozzles

The development and revision of safety codes and standards for hydrogen infrastructure requires a solid scientific basis, including studies of unignited releases from high pressure systems for various scenarios. Most hydrogen releases are modeled as axisymmetric jets, but real leaks are more likely to be non-axisymmetric jets issuing from high aspect ratio cracks or slots. In the present study, underexpanded hydrogen jets from square and rectangular nozzles with aspect ratios of 1 to 16 were numerically modeled for stagnation pressures up to 20 MPa. The near and far flow fields were modeled separately using two sequential computational domains to accurately and efficiently capture the flow characteristics. The numerical models were first validated with experimental data from a previous experimental study and literature data. The mass fraction and velocity distributions show that the centerline decay rates increase as the nozzle aspect ratio increases, but this increase is dependent on the pressure. This means that the canonical decay law of round turbulent jets and plumes no longer applies to the slot nozzle jets for high pressures. The radial profiles collapse onto a Gaussian curve in the minor axis plane, but neither collapse, nor are they Gaussian in the major axis plane with peaks away from the jet centerline. Different shock patterns were identified along the major and minor axes and the axis switching phenomenon seen in the literature was also reproduced. The axis switching resulted in significantly wider flattened concentration distributions compared with the axisymmetric jet which may require consideration during safety analyses for non-circular nozzles. A scaling factor taking both the nozzle shape and pressure effects into account SAND2018-12056J This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government. was then devleoped to better scale the centerline decay rates for jets from both the square and rectangular nozzles. The present study demonstrates that the nozzle shape effects on the jet spreading should not be overlooked and proper scaling factors are required to collapse the data and calculate decay rates.

08 HYDROGEN↗

DOE OSTI2023

Liquid hydrogen storage system for heavy duty trucks: Configuration, performance, cost, and safety

In this work, we investigate the potential of liquid hydrogen storage (LH 2 ) on-board Class-8 heavy duty trucks to resolve many of the range, weight, volume, refueling time and cost issues associated with 350 or 700-bar compressed H 2 storage in Type-3 or Type-4 composite tanks. We present and discuss conceptual storage system configurations capable of supplying H 2 to fuel cells at 5-bar with or without on-board LH 2 pumps. Structural aspects of storing LH 2 in double walled, vacuum insulated, and low-pressure Type-1 tanks are investigated. Structural materials and insulation methods are discussed for service at cryogenic temperatures and mitigation of heat leak to prevent LH 2 boiloff. Failure modes of the liner and shell are identified and analyzed using the regulatory codes and detailed finite element (FE) methods. The conceptual systems are subjected to a Failure modes and effects analysis (FMEA) and a safety, codes, and standards (SCS) review to rank failures and identify safety gaps. The results indicate that the conceptual systems can reach 19.6% usable gravimetric capacity, 40.9 g-H 2 /L usable volumetric capacity and $174-183/kg-H 2 cost (2016 USD) when manufactured 100,000 systems annually.

08 HYDROGEN↗

DOE OSTI2024

Effects of hydrogen partial pressure on crack initiation and growth rate in vintage X52 steel

There is an increasing interest in ensuring compatibility of existing natural gas infrastructure with conveyance of hydrogen, especially for vintage (pre-1970s) pipeline steels. While hydrogen is known to affect the fatigue and fracture resistance of these materials, there is a lack of data on its effects on crack initiation at low hydrogen partial pressure. To address this gap, the current work presents a series of fatigue tests on vintage API Grade X52 steel in pure nitrogen, air and gaseous hydrogen at partial pressures of 1 and 207 bar. Circumferential notch tension (CNT) specimens with direct current potential difference (DCPD) measurements are used to get S–N curves for crack initiation and failure, and compact tension (CT) tests to assess crack growth rates. Results indicate that air accelerates crack initiation, low hydrogen partial pressure affects the crack growth rate, and higher pressures have the largest impact on both. Surprisingly, crack initiation is faster in air than at 1 bar hydrogen partial pressure. The CNT tests also reveal that the orientation of the banded microstructure has no effect on the crack initiation site, but it influences its growth behavior. Moreover, fractography shows that hydrogen also promotes instances of quasi-cleavage and intergranular fracture, linked to embrittlement of the material. Here, these results demonstrate that, while hydrogen may impact the entirety of the fatigue process, its effects on the early stages of damage accumulation become relevant only at larger pressures. More comprehensive studies of crack initiation in hydrogen atmospheres will be critical to future codes and standards for hydrogen infrastructure.

Circumferential notch tension↗

DOE OSTI2024

Recent advances in combustion science related to hydrogen safety

Hydrogen is a key pillar in the global Net Zero strategy. Rapid scaling up of hydrogen production, transport, distribution and utilization is expected. This entails that hydrogen, which is traditionally an industrial gas, will come into proximity of populated urban areas and in some situations handled by the untrained public. To realize all their benefits, hydrogen and its technologies must be safely developed and deployed. The specific properties of hydrogen involving wide flammability range, low ignition energy and fast flame speed implies that any accidental release of hydrogen can be easily ignited. Comparing with conventional fuels, combustion systems fueled by hydrogen are also more prone to flame instability and abnormal combustion. This paper aims to provide a comprehensive review about combustion research related to hydrogen safety. It starts with a brief introduction which includes some overview about risk analysis, codes and standards. The core content covers ignition, fire, explosions and deflagration to detonation transition (DDT). Considering that DDT leads to detonation, and that detonation may also be induced directly under special circumstances, the subject of detonation is also included for completeness. The review covers laboratory, medium and large-scale experiments, as well as theoretical analysis and numerical simulation results. While highlights are provided at the end of each section, the paper closes with some concluding remarks highlighting the achievements and key knowledge gaps.

08 HYDROGEN↗

DOE OSTI2024

Ambient energy for buildings: Beyond energy efficiency

The following Key Messages comprise the salient findings of this study: 1. Ambient energy (from sun, air, ground, and sky) can heat and cool buildings; provide hot water, ventilation and daylighting; dry clothes; and cook food. These services account for about three-quarters of building energy consumption and a third of total US demand. Biophilic design (direct and indirect connections with nature) is an intrinsic adjunct to ambient energy systems, and improves wellness and human performance. 2. The current strategy of electrification and energy efficiency for buildings will not meet our climate goals, because the transition to an all-renewable electric grid is too slow. Widespread adoption of ambient energy is needed. Solar-heated buildings also flatten the seasonal demand for electricity compared to all-electric buildings, reducing required production capacity and long-term energy storage. In addition, ambient-conditioned buildings improve resilience by remaining livable during power outages. 3. National policies, incentives, and marketing should be enacted to promote ambient energy use. Federal administrative priorities should reflect the importance of ambient energy for buildings. Use of ambient energy should be encouraged through existing and new building codes and standards. 4. Ambient energy system design tools are needed for architects, engineers, builders, building scientists, realtors, appraisers, and consumers. PVWatts is used over 100 million times per year for photovoltaic system design. A similar, simple, and accessible tool for ambient design is crucial. 5. Training on ambient energy is needed throughout secondary, post-secondary, and continuing education for workforce development. Currently, only about 10% of colleges teach courses on passive heating and cooling systems. 6. Ambient-conditioned buildings should be demonstrated in all US climate zones. Performance should be monitored and reported, with quantitative case studies made widely available. 7. While current technology is sufficient to build high-performance ambient buildings now, research is needed to develop new technologies to harness ambient energy more effectively and more economically. Such advancements will facilitate adoption of ambient energy technologies in a wider range of buildings, including retrofits. Examples include windows with much lower thermal losses, use of the building shell as thermal storage, alternative light-weight thermal storage systems, sky radiation cooling systems, automated controls for solar gains and passive cooling, and ground coupling.

14 SOLAR ENERGY↗

DOE OSTI2021

Urban microclimate and its impact on building performance: A case study of San Francisco

Urban microclimate exerts an increasing influence on urban buildings, energy, and sustainability. This study uses 10-year measured hourly weather data at 27 sites in San Francisco, California, to (1) analyze and visualize the urban microclimate patterns and urban heat island effect; (2) simulate annual energy use and peak electricity demand of typical large office buildings and large hotels to investigate the influence of urban microclimate on building performance; (3) simulate indoor air temperature of a single-family house without air-conditioning during the record three-day heatwave of 2017, to quantify the divergence of climate resilience due to urban microclimate effect. Results show significant microclimate effects in San Francisco with up to 11 °C outdoor air temperature difference between the coastal and downtown areas on September 1, 2017, during the record three-day heatwave. The simulated energy results of the prototype large office and large hotel buildings using the 2017 weather data show over 100% difference in annual heating energy use and 65% difference in annual cooling energy use across different stations; as well as up to 30% difference in peak cooling electricity demand. The impacts on annual site or source energy use are minimal (less than 5%) as cooling and heating in a mild climate are a relatively small portion of overall building energy use in San Francisco. Results also show the microclimate effects influence indoor air temperature of unconditioned homes by up to 5 °C. Newer buildings and homes are much less affected by microclimate effects due to more stringent performance requirements of the building envelope and energy systems. Furthermore, these findings inform that San Francisco microclimate variations should be considered in urban energy planning, building energy codes and standards, as well as heat resilience policymaking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Thoroughly testing and integrating hundreds of Pull Requests per month: ROOT’s new Cost-efficient and Feature Rich GitHub-based CI

ROOT is an open source framework, freely available on GitHub, at the heart of data acquisition, processing and analysis of HE(N)P experiments, and beyond. It is developed collaboratively: contributions are not authored only by ROOT team members, but also by the user community at large: developers and scientists from universities, labs as well as the private sector. More than 1500 GitHub Pull Requests are merged on average per year. It is in this context that code integration acquires a primary role. The review of code contributions isn’t enough: not only they need to be thoroughly reviewed, they also need to be thoroughly tested through a powerful CI infrastructure on several different platforms to comply with the high code quality standards of the project. Since the end of 2023, ROOT moved its continuous integration system from Jenkins to GitHub Actions. In this contribution, we characterise the transition to the GitHub CI, focussing on our strategy, its implementation and the lessons learned, as well as the advantages the new system offers with respect to the previous one. Particular emphasis will be given to the evaluation of the cost-benefit ratio for Jenkins and GitHub Actions for the ROOT project. We also describe how we manage to run in less than one hour thousands of unit, integration, functional and end-to-end tests on different flavours of Windows, four versions of macOS, as well as about ten of the most used Linux distributions, taking advantage of the CERN computing infrastructure.

Piparo, Danilo [CERN]↗

DOE OSTI2016

Ion-kinetic simulations of D- 3 He gas-filled inertial confinement fusion target implosions with moderate to large Knudsen number

Experiments designed to investigate the transition to non-collisional behavior in D 3 He-gas inertial confinement fusion target implosions display increasingly large discrepancies with respect to simulations by standard hydrodynamics codes as the expected ion mean-free-paths λ c increase with respect to the target radius R (i.e., when the Knudsen number N K = λ c /R grows). To take properly into account large N K 's, multi-ion-species Vlasov-Fokker-Planck computations of the inner gas in the capsules have been performed, for two different values of N K , one moderate and one large. The results, including nuclear yield, reactivity-weighted ion temperatures, nuclear emissivities, and surface brightness, have been compared with the experimental data and with the results of hydrodynamical simulations, some of which include an ad hocmodeling of kinetic effects. The experimental results are quite accurately rendered by the kinetic calculations in the smaller-N K case, much better than by the hydrodynamical calculations. The kinetic effects at play in this case are thus correctly understood. However, in the higher-N K case, the agreement is much worse. Furthermore, the remaining discrepancies are shown to arise from kinetic phenomena (e.g., inter-species diffusion) occurring at the gas-pusher interface, which should be investigated in the future work.

42 ENGINEERING↗

DOE OSTI2018

Speeding up simulations by slowing down particles: Speed-limited particle-in-cell simulation

Based on the particle-in-cell (PIC) plasma simulation method, the speed-limited PIC (SLPIC) method delivers faster kinetic plasma simulation in cases where the particle distributions evolve slowly compared with the maximum stable PIC timestep. SLPIC thus offers more feasible, fully kinetic simulation in regimes that historically have required fluid approaches, such as magnetohydrodynamic (MHD), two-fluid, or Boltzmann electron treatments. In particular, SLPIC allows an explicit time advance with steps much larger than the inverse plasma frequency, avoiding the instability explicit PIC faces with large timesteps. SLPIC avoids this instability by slowing down fast particles (e.g., electrons) in a way that is rigorously underpinned by an approximate Vlasov equation; unlike MHD, two-fluid, and Boltzmann electron approaches, SLPIC does not fundamentally neglect any first-principles plasma physics, although the choices of the grid cell size, timestep, and number of macroparticles per cell naturally limit the physical phenomena that can be accurately represented. SLPIC can be implemented with minor modifications of a standard PIC code and does not require an implicit time advance. It enables large timesteps in first-principles kinetic plasma simulation of appropriately slow phenomena, and it can handle many of the same complications as PIC, such as boundary conditions and collisions. Finally, in an argon plasma sheath test problem, a SLPIC simulation achieved a speed-up of a factor of 160 over the corresponding PIC simulation, without loss of accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Design and Analysis of the US ITER Central Solenoid Module Lifting Fixture

The central solenoid (CS) consists of six large high field superconducting magnets (also known as modules) approximately 4 m in diameter and 2 m tall that weigh approximately 120 tonnes each. These large and complex modules create challenges during assembly of the CS that require the development of custom assembly tooling such as the CS lifting fixture. The CS module lifting fixture is designed to lift and stack the six CS modules in the assembly building on the ITER site. Because of its unique design, fabrication, and assembly features, no lifting attachments could be incorporated within or under the CS modules. This limitation motivated the development of a friction-based lift fixture. The design and evaluation of the CS module lifting fixture considered both worker safety and investment protection, and the assessments were performed to international codes and standards. The CS module lifting fixture consists of two principal subassemblies: spider assembly and ring weldment. These subassemblies allow the frictional force to be augmented by the mechanical advantage of shallow-angle wedges. Large radial preloads created by both screw jack assemblies and the weight of a CS module develop frictional forces capable of performing a lift with a safety factor of at least 2. The design effort resulted in the use of low friction linear bearings on angle surfaces to ensure constant pressure, integrated jacks for pretensioning the fixture prior to lifting, and load pin strain gauges for monitoring the normal force. Testing of various materials and surface treatments led to the selection of laminated aluminum and rubber pads as the friction interface on the CS lifting fixture side and a grit-blasted Nitronic 50 stainless steel band on the CS module side. A redundant lifting method using the module slings between the spider and module is also utilized after the initial friction lift. The CS lifting fixture provides a safe and reliable solution for lifting and stacking the CS module during assembly.

Morris, Eric↗

DOE OSTI2021

Structural Analysis of the US ITER ECH Transmission Line System

The electron cyclotron heating (ECH) and current drive system is one of the main plasma heating systems for ITER. It uses high-power microwave beams with the power deposition location steerable across the plasma cross section. Microwave power is conveyed via transmission lines (TLs) that run from the gyrotrons in the radio frequency building through the assembly hall and tokamak building to the ECH launchers within the tokamak vacuum vessel. The ECH system includes a vast array of interconnected TL waveguides, in-line components, and support structures.Finite element (FE) modeling provides an essential means of simulating the system, applying loads and determining deflections, rotations, forces, moments, and stresses in order to evaluate various structural and microwave transmission performance metrics. A representative FE model of the overall ECH TL system is developed in ANSYS®. This top-level model defines the centerline of the waveguide system. Waveguide segments are represented by line elements (beams and pipes) with equivalent section properties, and support structures are represented by boundary conditions. A systematic approach is used to model each ECH component with lumped masses and structurally equivalent stiffness matrices or ANSYS superelements.The top-level TL FE model is used to evaluate the various loads (thermal, vacuum, seismic, etc.) and operating scenarios. The top-level model directly calculates stresses in the straight aluminum waveguide segments. The model provides the forces and moments acting on the in-line components for detailed submodel assessments. Displacement results from the top-level analysis feed into a separate microwave performance model to help determine operational efficiency. All TL performance and thermal-structural requirements are met, as specified by the applicable codes and standards, and successfully documented in numerous technical reports and demonstrated at the final design review.

Shanmugasundaram, Aravind↗

DOE OSTI2015

ON VERTICALLY GLOBAL, HORIZONTALLY LOCAL MODELS FOR ASTROPHYSICAL DISKS

Disks with a barotropic equilibrium structure, for which the pressure is only a function of the density, rotate on cylinders in the presence of a gravitational potential, so that the angular frequency of such a disk is independent of height. Such disks with barotropic equilibria can be approximately modeled using the shearing box framework, representing a small disk volume with height-independent angular frequency. If the disk is in baroclinic equilibrium, the angular frequency does generally depend on height, and it is thus necessary to go beyond the standard shearing box approach. In this paper, we show that given a global disk model, it is possible to develop approximate models that are local in horizontal planes without an expansion in height with shearing-periodic boundary conditions. We refer to the resulting framework as the vertically global shearing box (VGSB). These models can be non-axisymmetric for globally barotropic equilibria but should be axisymmetric for globally baroclinic equilibria. We provide explicit equations for this VGSB which can be implemented in standard magnetohydrodynamic codes by generalizing the shearing-periodic boundary conditions to allow for a height-dependent angular frequency and shear rate. We also discuss the limitations that result from the radial approximations that are needed in order to impose height-dependent shearing periodic boundary conditions. We illustrate the potential of this framework by studying a vertical shear instability and examining the modes associated with the magnetorotational instability.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Simulation of site and topographic effects on ground motion in Los Alamos, NM mesas

Mesas, which consist of an elevated area of land with a flat top and steep cliffs at the sides, are one of the common geological formations present in the Los Alamos region. Previous research has shown that geological formations such as sedimentary canyons can amplify soil response during earthquakes. There have also been parametric studies to understand the response of an idealized and isolated mountain/canyon under inclined plane waves. In this study, a 2-D linear anelastic soil domain, with topography modelled after the Los Alamos region with many mesas and canyons, is considered to understand site-specific topographic effects in the presence of non-isolated topographical features. Various earthquake sources configurations ranging from point sources to finite faults with varying rupture length, dip angles and spatial slip distributions are considered. To isolate the effect of topography, three other soil domains—a homogeneous half-space, homogeneous half-space with mesas and canyons, and a computational domain with just the soil properties from the Los Alamos region on a flat soil domain is also considered. 2-D site-response analyses of these soil domains under earthquake excitation show that the free-field response of the soil can be amplified or de-amplified depending on the topography of the region and the location of the station relative to the fault. These studies also show that even relatively small mesas with height less than 100 m can significantly amplify the response (by a factor of 2 or more), which differ from the much smaller amplification factors (≤1.2) specified by standard building codes such as Eurocode-8 (2000) for topographical features with a similar slope. De-amplifications are also fairly common, especially close to the canyons. The results from this study agree qualitatively with those from the ambient vibration study conducted by Stolte et al. on a mesa from the Los Alamos region. Furthermore, such site-specific studies provide important insights into the variability of the topographic amplification factors within a region of interest. Such knowledge is important in the design of safety-related critical infrastructure located within that region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗