Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation 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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

TOPOS: A new topometric patient positioning and tracking system for radiation therapy based on structured white light

Purpose: A patient positioning system for radiation therapy based on structured white light and using off-the-shelf hardware components for flexibility and cost-effectiveness has been developed in house. Increased accuracy, patient comfort, abandonment of any skin marks, accelerated workflow, objective reading/recording, better usability and robust sensor design, compared to other positioning approaches, were the main goals of this work. Another aim was the application of a 6 degrees of freedom tracking system working without dose deposition. Methods: Two optical sensors are the main parts of the TOPOS Registered-Sign system (Topometrical Positioning, cyberTECHNOLOGIES, Germany). The components: cameras, projectors, and computers are commercial off-the-shelf products, allowing for low production costs. The black/white cameras of the prototype are capable of taking up to 240 frames per second (resolution: 640 Multiplication-Sign 488 pixels). The projector has a resolution of 1024 Multiplication-Sign 768 and a refresh rate of 120 Hz. The patient's body surface is measured continuously and registered to a reference surface, providing a transformation to superimpose the patient's surface to the reference (planning CT) surface as best as possible. The execution of the calculated transformation provides the correct patient position before the treatment starts. Due to the high-speed acquisition of the surfaces, a surveillance of the patient's (respiration) motion during treatment is also accomplished. The accuracy of the system was determined using a male mannequin. Two treatment sites were evaluated: one simulating a head and neck treatment and the other simulating a thoracic wall treatment. The mannequin was moved to predefined positions, and shift vectors given by the surface registration were evaluated. Additionally manual positioning using a color-coding system was evaluated. Results: Two prototypes have been developed, each allowing a continuous high density scan of a 500 Multiplication-Sign 500 Multiplication-Sign 400 mm{sup 3} (L Multiplication-Sign W Multiplication-Sign D) large volume with a refresh rate of 10 Hz (extendible to 20 Hz for a single sensor system). Surface and position correction display, as well as respiratory motion, is shown in real-time (delay < 200 ms) using present graphical hardware acceleration. For an intuitive view of the patient's misalignment, a fast surface registration algorithm has been developed and tested and a real-time color-coding technique is proposed and verified that allows the user to easily verify the position of the patient. Using first the surface registration and then the color coding the best results were obtained: for the head and neck case, the mean difference between the actual zero position and the final match was 0.1 {+-} 0.4, -0.2 {+-} 0.7, and -0.1 {+-} 0.3 mm in vertical, longitudinal, and lateral direction. For the thoracic case, the mean differences were 0.3 {+-} 0.5, -0.6 {+-} 1.9, 0.0 {+-} 0.4 mm. Conclusions: The presented system copes with the increasing demand for more accurate patient positioning due to more precise irradiation technologies and minimizes the preparation times for correct patient alignment, therefore optimizing the treatment workflow. Moreover, TOPOS is a versatile and cost effective image guided radiation therapy device. It allows an objective rating of the patient's position before and during the irradiation and could also be used for respiratory gating or tracking.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-G-TeP1-06: Fast GPU Framework for Four-Dimensional Monte Carlo in Adaptive Intensity Modulated Proton Therapy (IMPT) for Mobile Tumors

Purpose: To demonstrate the feasibility of fast Monte Carlo (MC) treatment planning and verification using four-dimensional CT (4DCT) for adaptive IMPT for lung cancer patients. Methods: A validated GPU MC code, gPMC, has been linked to the patient database at our institution and employed to compute the dose-influence matrices (Dij) on the planning CT (pCT). The pCT is an average of the respiratory motion of the patient. The Dijs and patient structures were fed to the optimizer to calculate a treatment plan. To validate the plan against motion, a 4D dose distribution averaged over the possible starting phases is calculated using the 4DCT and a model of the time structure of the delivered spot map. The dose is accumulated using vector maps created by a GPU-accelerated deformable image registration program (DIR) from each phase of the 4DCT to the reference phase using the B-spline method. Calculation of the Dij matrices and the DIR are performed on a cluster, with each field and vector map calculated in parallel. Results: The Dij production takes ∼3.5s per beamlet for 10e6 protons, depending on the energy and the CT size. Generating a plan with 4D simulation of 1000 spots in 4 fields takes approximately 1h. To test the framework, IMPT plans for 10 lung cancer patients were generated for validation. Differences between the planned and the delivered dose of 19% in dose to some organs at risk and 1.4/21.1% in target mean dose/homogeneity with respect to the plan were observed, suggesting potential for improvement if adaptation is considered. Conclusion: A fast MC treatment planning framework has been developed that allows reliable plan design and verification for mobile targets and adaptation of treatment plans. This will significantly impact treatments for lung tumors, as 4D-MC dose calculations can now become part of planning strategies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Data-Driven Modeling of Spatially Explicit Air Pollutant Inventories for Agricultural and Forestry Feedstock Production

As biomass feedstock production expands, air pollutants generated from feedstock supply chain activities may impact air quality and human health. The National Renewable Energy Laboratory's Feedstock Production Emissions to Air Model (FPEAM) generates county-scale inventories of direct air pollutant emissions produced by feedstock supply chain activities. FPEAM can also be used to generate local, temporally resolved inventories for input into air quality models. This presentation provides an overview of FPEAM functionality and discusses several recent use cases prior to FPEAM's scheduled public release in September 2019. FPEAM is a modular, data-driven model. Each module contains calculations for an activity type, except for the feedstock transportation and on-farm equipment operation modules which connect to the U.S. Environmental Protection Agency's Motor Vehicle Emissions Simulator (MOVES) and Nonroad models, respectively. The data-driven aspect of FPEAM means that it is the input data provided for a scenario rather than the model itself that determines the pollutants and supply chain activities included in the scenario scope. Modules within FPEAM can also be turned on or off manually to adjust the scenario scope. Default input data packaged with the FPEAM code base can be used to calculate inventories of volatile organic compounds, particulate matter, ammonia, nitrogen oxides, sulfur oxides, and carbon monoxide generated by agricultural and forestry equipment operation, application of nitrogenous fertilizers, herbicide and insecticide, and feedstock transportation. Custom, user-developed data can be used in place of the default data to analyze scenarios with additional or alternate feedstocks, pollutants, and supply chain activities.

09 BIOMASS FUELS↗
DOE OSTI2020

Strategies for thorium fuel cycle transition in the SD-TMSR

Liquid-fueled Molten Salt Reactor (MSR) systems represent advances in safety, economics, and sustainability. The MSR has been designed to operate with a Th/ 233 U fuel cycle with 233 U used as startup fissile material. Since 233 U does not exist in nature, we must examine other available fissile materials to start up these reactor concepts. This work investigates the fuel cycle and neutronics performance of the Single-fluid Double-zone Thorium-based Molten Salt Reactor (SD-TMSR) with different fissile material loadings at startup: High Assay Low Enriched Uranium (HALEU) (19.79%), Pu mixed with HALEU (19.79%), reactor-grade Pu (a mixture of Pu isotopes chemically extracted from Pressurized Water Reactor (PWR) spent nuclear fuel (SNF) with 33 GW d/tHM burnup), transuranic elements (TRU) from Light Water Reactor (LWR) SNF, and 233 U. The MSR burnup routine provided by SERPENT-2 is used to simulate the online reprocessing and refueling in the SD-TMSR. The effective multiplication factor, fuel salt composition evolution, and net production of 233 U are studied in the present work. Additionally, the neutron spectrum shift during the reactor operation is calculated. The results show that the continuous flow of reactorgrade Pu helps transition to the thorium fuel cycle within a relatively short time (≈ 4.5 years) compared to 26 years for 233 U startup fuel. Finally, using TRU as the initial fuel materials offers the possibility of operating the SD-TMSR for an extended period of time (≈ 40 years) without any external feed of 233 U.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Scalar mixing in a Kelvin-Helmholtz shear layer and implications for Reynolds-averaged Navier-Stokes modeling of mixing layers

Large-eddy simulation of a temporally evolving Kelvin-Helmholtz (KH) mixing layer is performed with the tenth-order compact difference code miranda to examine the steady-state behavior of a passive scalar in a shear-driven mixing layer. It is shown that the integral behavior of scalar variance in a KH mixing layer behaves similarly to the integral behavior of scalar variance in a Rayleigh-Taylor (RT) mixing layer, and mixedness of the simulated KH shear layer tends towards a value of about 0.8. It is further shown that if the k-L-a-V Reynolds-averaged Navier-Stokes (RANS) model [B. E. Morgan et al., Phys. Rev. E 98, 033111 (2018)], calibrated to reproduce steady-state mixing in an RT layer, is applied to simulate a KH mixing layer, the RANS model will significantly overpredict the magnitude of scalar variance in the KH layer. A straightforward addition to the k-L-a-V model is then suggested, and self-similarity analysis is applied to determine constraints on model coefficients. Furthermore, it is shown that with the addition of a buoyancy production term in the model equation for scalar variance, it becomes possible to eliminate the model deficiency and match steady-state mixedness in simulations of both RT and KH mixing layers with a single model calibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Requirements Description of GAMSOR and GAMSRC

This report presents the modeling and simulation capabilities of Argonne National Laboratory’s GAMSOR and GAMSRC code that is used in present reactor analysis activities. These capabilities will then be used to establish the set of verification tasks necessary to verify GAMSOR for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The GAMSOR program is a modified version of DIF3D used to obtain the coupled neutron gamma flux distributions for a given reactor system. GAMSOR was first developed in the 1980s and its work has progressed in parallel with DIF3D throughout its various development phases. While GAMSOR is compatible with DIF3D-FD and DIF3D-Nodal, DIF3D-VARIANT is most commonly used today and will be the focus of the present verification effort. As covered in [1], GAMSOR was built to provide more accurate heating distributions by explicitly solving for both the gamma and neutron heating distributions. Without this consideration, gamma heat is normally assumed to be deposited locally at the site of gamma production. This is acceptable for modeling some phenomena because the total energy production is preserved, but it can lead to inaccurate approximations of fluence-dependent behaviors such as fuel performance, swelling, cladding embrittlement, radiation induced creep, and power peaking. The outputs from GAMSOR are contained in multiple intermediate and final output files: FIXSRC, VARSRC, GTFLUX, GHFLUX, PWDINT, NPDINT, and GPDINT. DIF3D normally produces a PWDINT file containing power density results, with the assumption that gamma rays are absorbed where they are emitted (i.e., energy deposition is local with no modeled transport). With GAMSOR, neutron power density (NPDINT) and gamma power density (GPDINT) are both produced, with their summed power densities replacing the results in PWDINT for total heating (i.e., recoverable) power density.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Countercurrent Actinide Lanthanide Separation Process (ALSEP) Demonstration Test with a Simulated PUREX Raffinate in Centrifugal Contactors on the Laboratory Scale

An Actinide Lanthanide Separation Process (ALSEP) for the separation of trivalent actinides (An(III)) from simulated raffinate solution was successfully demonstrated using a 32-stage 1 cm annular centrifugal contactor setup. The ALSEP solvent was composed of a mixture of 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) and N,N,N',N'-tetra-(2-ethylhexyl)-diglycolamide (T2EHDGA) in n-dodecane. Flowsheet calculations and evaluation of the results were done using the Argonne’s Model for Universal Solvent Extraction (AMUSE) code using single-stage distribution data. The co-extraction of Zr(IV) and Pd(II) was prevented using CDTA (trans-1,2-diaminocyclohexane-N,N,N',N'-tetraacetic acid) as a masking agent in the feed. For the scrubbing of co-extracted Mo; citrate-buffered acetohydroxamic acid was used. The separation of An(III) from the trivalent lanthanides (Ln(III)) was achieved using citrate-buffered diethylene-triamine-N,N,N',N",N"-pentaacetic acid (DTPA), and Ln(III) were efficiently back extracted using N,N,N',N'-tetraethyl-diglycolamide (TEDGA). A clean An(III) product was obtained with a recovery of 95% americium and curium. The Ln(III) were efficiently stripped; but the Ln(III) product contained 5% of the co-stripped An(III). The carryover of Am and Cm into the Ln(III) product is attributed to too few actinide stripping stages, which was constrained by the number of centrifugal contactors available. Improved separation would be achieved by increasing the number of An strip stages. The heavier lanthanides (Pr, Nd, Sm, Eu, and Gd) and yttrium were mainly routed to the Ln product, whereas the lighter lanthanides (La and Ce) were mostly routed to the raffinate.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Neutron and Gamma Correlations using CGM in MCNP 6.2.0

The general-purpose Monte Carlo radiation-transport code MCNP6{sup TM} came from the integration of MCNP5{sup TM} and MCNPX{sup TM}. The most recent release of MCNP, version 6.2.0, is scheduled for release in early 2017, which includes a variety of bug fixes and new features to the code. MCNP traditionally models nuclear reactions using Monte Carlo sampling techniques on measured and modeled cross-sectional data contained within ACE (A Compact Evaluated Nuclear Data File) libraries. In addition to reaction cross-sections, ACE libraries contain information relating to the production of secondary particles, such as secondary neutrons and photons. For example, an incident neutron on a nucleus can produce a nuclear reaction in the form (n, M{sub n}n'), where n is the incident neutron and M{sub n} is the number of secondary neutrons n' produced. Using this notation, different reactions can be described based on their neutron multiplicity, such as neutron capture (M{sub n} = 0), elastic and inelastic scattering (M{sub n} = 1), and reactions with multiple secondary neutrons (M{sub n} > 1). Accounting for secondary γ-rays emitted from the residual nuclei, this reaction can further be generalized as n, Mnn'M to describe a reaction producing M{sub n} neutrons and M -rays. These values M{sub n} and M are referred to as the neutron and multiplicities, respectively. ACE libraries contain the cross-sectional data for each reaction, and their corresponding statistically averaged multiplicities for neutrons and photons, M-bar{sub n} and M-bar{sub γ}. However, there are several limitations to reaction sampling using ACE libraries. First, the data contained within ACE libraries is limited by the accuracy of the model and/or experiment, which out of necessity can greatly simplify the true physics behind a reaction. For example, the value M{sub γ} in a (n, M{sub n}n'M{sub γ}γ) reaction is determined by the statistically averaged multiplicity M-bar{sub γ} instead of a distribution. Consequently, MCNP, unable to produce a continuous distribution, will utilizes a binary sampling around M-bar{sub γ}. While this sampling method provides a statistically average M{sub γ} over a large simulation, it is an inaccurate representation of an individual reaction and its corresponding secondary particles. Secondary neutrons are affected in a similar manner, such that different kinds of secondary neutrons can be sampled from the same interaction (eg., an inelastic neutron with an [n, 2n] neutron). Lastly, the random sampling technique for secondary particles remains completely independent of the sampled reaction, eliminating the ability to correlate secondary particles for a given reaction. To supplement these data libraries, additional codes and models describing certain nuclear reactions are employed. These codes can fill in the information missing from ACE, provide a better model to the true physics distribution, and can provide true correlated secondary particles, which has many useful applications in research and industry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Time dependent supervisory control update with FARM using rolling window

This report describes improvements to the Feasible Actuator Range Modifier (FARM) component of the RAVEN-based HYBRID framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON plug-in that solves the power dispatch problem. The solution to the dispatch problem involves economically optimal dispatches that satisfy limits on production variables and their rates of variation (explicit constraints) as well as process variables tied to the service life of equipment (implicit constraints). FARM serves to validate or confirm that a HERON solution for explicit constraints also satisfies the implicit constraints. FARM-alpha was released by Argonne National Laboratory in January 2021 followed by FARM-Beta in January 2022 with the latter providing increased flexibility for the user. In this report, FARM-Gamma, the latest version of the code, is described. The major improvement is the implementation of a system identification algorithm based on the Dynamic Mode Decomposition with Control (DMDc) coupled with a “Rolling Window” scheme that allows obtaining linear time-varying state-space models. This feature equips FARM with the most accurate approximation of system dynamics, and it relieves the user from the burden of performing an exhaustive off-line characterization of the dynamics. FARM-Gamma capabilities are assessed by solving the power dispatch problem for a representative IES unit. The simulation times corresponding to the different releases are estimated and compared. These values capture the increasing computational burden of the successively higher fidelity state-space models adopted by FARM-Alpha, FARM-Beta and FARM-Gamma. The code implementation provides significant flexibility, i.e., the user can always select the most suitable version of FARM according to the problem to be solved and the available computational resources. It is anticipated that FARM will play a role in addressing several future IES applications. We outline how it can support the coordinated management and safe operation of a nuclear plant coupled to industrial processes to produce hydrogen and synfuels.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

What Processes are Defining the Ionospheric Conductivity and its Variability During Geomagnetic Disturbances?

Modeling of electrodynamic coupling between the magnetosphere, ionosphere, and upper atmosphere (MIA) depends on accurate specification of ionospheric conductance produced by auroral precipitation of high-energy electrons and ions. The precipitation of energetic electrons into the ionosphere is the result of a three-step process that relies on the proper selections of the simulation tools for the ionospheric conductivity studies, while observation can only measure the results of the three steps. In the region of diffuse aurora, the first step is the initiation of electron precipitation into both magnetically conjugate foot points from the Earth’s magnetosphere via wave-particle interactions. The second step is the multiple atmospheric backscatters (or reflections) of electrons at the two magnetic conjugate points, which produces secondary superthermal electron fluxes. The third step is namely the self-consistent electric and magnetic fields that influence magnetospheric particle transport and re-distribute precipitating electrons and ions through the ionospheric electrodynamics. These steps are especially important for revealing electron precipitation dynamics that carry most of the energy in the aurora, resulting also from ion precipitation production and the formation of ionospheric conductance during geomagnetic disturbances. We demonstrate all above results based on SuperThermal Electron Transport (STET), Superthermal Proton, Electron and Atomic Hydrogen tRansport in the Ionosphere and Thermosphere (SPEAH-RIT), and Comprehensive Inner Magnetosphere and Ionospere (CIMI) codes developed at NASA Goddard Space Flight Center.

George V. Khazanov↗
DOE OSTI2021

Cognitive simulation models for inertial confinement fusion: Combining simulation and experimental data

The design space for inertial confinement fusion (ICF) experiments is vast, and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design and investigation, simulations require input from past experimental data to better predict future performance. In this work, we describe a cognitive simulation method for combining simulation and experimental data into a common, predictive model. This method leverages a machine learning technique called “transfer learning,” the process of taking a model trained to solve one task, and partially retraining it on a sparse dataset to solve a different, but related task. In the context of ICF design, neural network models are trained on large simulation databases and partially retrained on experimental data, producing models that are far more accurate than simulations alone. Here, we demonstrate improved model performance for a range of ICF experiments at the National Ignition Facility and predict the outcome of recent experiments with less than 10% error for several key observables. We discuss how the methods might be used to carry out a data-driven experimental campaign to optimize performance, illustrating the key product—models that become increasingly accurate as data are acquired.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Fault-tolerant resource comparison of qudit and qubit encodings for diagonal quadratic operators

Finite local Hilbert-space truncations arise naturally in quantum simulations of lattice field theories and motivate qudit encodings, but their fault-tolerant advantage over qubit encodings remains unclear. We compare the non-Clifford cost of implementing quadratic diagonal evolutions, exemplified by 𝑈 = 𝑒$^{−𝑖⁢𝑡⁢𝜙^2_𝑥}$ in a uniform field-amplitude discretization of a real scalar field, using either one logical 𝑑-level qudit or 𝑛 𝑏 = ⌈log 2⁡ 𝑑⌉ logical qubits. We analyze two standard settings: product-formula simulation and linear combination of unitaries (LCU) per block encoding, taking the resource metric to be the number of non-Clifford gates after synthesis into a discrete logical gate set. Because tight synthesis bounds for general single-qudit rotations are not known, we express the qudit constructions in terms of embedded two-level SU⁡(2) rotations and derive explicit finite-𝑑 break-even conditions for their synthesis cost; these serve as compiler targets for when qudit encodings can outperform the qubit baseline. Within the constructive models studied here, product-formula implementations would require an exponentially stronger per-primitive synthesis advantage for qudits to win asymptotically, while in the LCU setting the qubit encoding is asymptotically cheaper in 𝑑. Nevertheless, the finite-𝑑 threshold analysis identifies low-dimensional regions in which qudits can yield meaningful constant-factor savings, particularly for LCU-based implementations. As a secondary analysis of the LCU construction, we use an idealized negligible-overhead qubit-qudit code-switching model to give an absolute 𝑇-count comparison and reinterpret the savings as an allowable per-switch overhead budget.

Godwood, Samuel [Univ. of Liverpool (United Kingdo↗
DOE OSTI2012

Experimental Definition and Validation of Protein Coding Transcripts in Chlamydomonas reinhardtii

Algal fuel sources promise unsurpassed yields in a carbon neutral manner that minimizes resource competition between agriculture and fuel crops. Many challenges must be addressed before algal biofuels can be accepted as a component of the fossil fuel replacement strategy. One significant challenge is that the cost of algal fuel production must become competitive with existing fuel alternatives. Algal biofuel production presents the opportunity to fine-tune microbial metabolic machinery for an optimal blend of biomass constituents and desired fuel molecules. Genome-scale model-driven algal metabolic design promises to facilitate both goals by directing the utilization of metabolites in the complex, interconnected metabolic networks to optimize production of the compounds of interest. Using Chlamydomonas reinhardtii as a model, we developed a systems-level methodology bridging metabolic network reconstruction with annotation and experimental verification of enzyme encoding open reading frames. We reconstructed a genome-scale metabolic network for this alga and devised a novel light-modeling approach that enables quantitative growth prediction for a given light source, resolving wavelength and photon flux. We experimentally verified transcripts accounted for in the network and physiologically validated model function through simulation and generation of new experimental growth data, providing high confidence in network contents and predictive applications. The network offers insight into algal metabolism and potential for genetic engineering and efficient light source design, a pioneering resource for studying light-driven metabolism and quantitative systems biology. Our approach to generate a predictive metabolic model integrated with cloned open reading frames, provides a cost-effective platform to generate metabolic engineering resources. While the generated resources are specific to algal systems, the approach that we have developed is not specific to algae and can be readily expanded to other microbial systems as well as higher plants and animals.

09 BIOMASS FUELS↗
DOE OSTI2023

Advanced Photon Acceleration Schemes for Tunable XUV/Soft X-Ray Sources

We report on the advances made under the purview of DOE award number DE-SC0019135 that was active during the period: 09/01/2018–08/31/2021. The grant investigated the application of a “flying focus” to the problem of photon acceleration—in which a dynamic refractive index gradient is used to continuously upshift a probe beam’s frequency. Codes were written to describe the creation of ionization waves of arbitrary velocity (IWAVs) for use as a photon-accelerating medium, as well as the behavior of a witness pulse residing in said medium. Experiments first verified the spatiotemporal control over laser intensity provided by a chromatic flying focus, then used that ability to produce small-diameter IWAVs in the far-field with the expected dynamics, and finally demonstrated even further flexibility by producing large-diameter IWAVs in the laser quasi-far-field that maintained the beneficial dynamics. Multiple innovative diagnostics—spectrally resolved Schlieren and spectrally resolved interferometry—were pioneered in order to diagnose the IWAVs. For IWAV production in the laboratory, however, beam quality was identified as a key limitation in the quasi-far-field. Since the original chromatic flying focus was found to result in relatively long (ps duration) intensity peaks, which could limit some applications including photon acceleration, additional techniques were invented to provide similar spatiotemporal control while also retaining ultrashort intensity peaks. While simulations have identified several interesting regimes for photon acceleration—first predicting the upshift of a counterpropagating witness pulse from the optical to the extreme ultraviolet in less than 1 cm, and later obtaining similar shifts in less than 100 μm in a simpler self-seeded configuration—experimental demonstration is left for future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Utah FORGE Project 3-2417: Simulations for Distributed Acoustic Sensing Strain Signatures as an Indicator of Fracture Connectivity

This dataset encompasses simulations of strain signatures from both hydraulically connected and "near-miss" fractures in enhanced geothermal systems (EGS). The files and results are presented from the perspective of digital acoustic sensing's (DAS) potential to differentiate the two fracture types. This dataset was acquired by the FOGMORE R&D project (Fiber Optic Geophysical MOnitoring of Reservoir Evolution), under Utah FORGE R&D Project 3-2417. Included are simulation and results via MatLab and COMSOL files, as well as a thesis and paper summarizing the results. Some stimulated fractures may be incomplete, approaching but not intersecting the production well. These "near-miss" fractures can be addressed in future stimulation stages or re-stimulated to complete the connection. We propose the use of fiber optic distributed acoustic sensing (DAS) as a method by which near-miss stimulated fractures may be identified and distinguished from hydraulically connected fractures. The low-frequency sub-nanostrain signatures of both complete and near-miss fractures in DAS data are simulated in this study using a hydrogeomechanical discrete fracture network model. The spatial distribution of strain was found to be an accurate indicator. However, this indicator must be evaluated in the context of DAS gauge length and spatial sampling. These simulations are a precursor to tests conducted at FORGE in 2023.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Developments in Performance and Portability for MadGraph5_aMC@NLO

Event generators simulate particle interactions using Monte Carlo techniques, providing the primary connection between experiment and theory in experimental high energy physics. These software packages, which are the first step in the simulation worflow of collider experiments, represent approximately 5 to 20% of the annual WLCG usage for the ATLAS and CMS experiments. With computing architectures becoming more heterogeneous, it is important to ensure that these key software frameworks can be run on future systems, large and small. In this contribution, recent progress on porting and speeding up the Madgraph5_aMC@NLO event generator on hybrid architectures, i.e. CPU with GPU accelerators, is discussed. The main focus of this work has been in the calculation of scattering amplitudes and "matrix elements", which is the computational bottleneck of an event generation application. For physics processes limited to QCD leading order, the code generation toolkit has been expanded to produce matrix element calculations using C++ vector instructions on CPUs and using CUDA for NVidia GPUs, as well as using Alpaka, Kokkos and SYCL for multiple CPU and GPU architectures. Performance is reported in terms of matrix element calculations per time on NVidia, Intel, and AMD devices. The status and outlook for the integration of this work into a production release usable by the LHC experiments, with the same functionalities and very similar user interfaces as the current Fortran version, is also described.

Valassi, Andrea↗
DOE OSTI2015

2015 IECC: Energy Savings Analysis

Section 304(a) of the Energy Conservation and Production Act, as amended, requires the Secretary of Energy to make a determination each time a revised edition of the 1992 Model Energy Code (MEC), or any successor thereof, is published with respect to whether the revised code would improve energy efficiency in residential buildings. The International Energy Conservation Code (IECC), as administered by the International Code Council (ICC), establishes the national model code for energy efficiency requirements for residential buildings. The latest edition of the IECC, the 2015 IECC, was published on June 3, 2014 and forms the basis of this analysis. To meet these statutory requirements, as well as to assist states and adopting entities in understanding associated savings, the DOE Building Energy Codes Program and Pacific Northwest National Laboratory (PNNL) conduct analyses to evaluate the differences between the latest edition of the IECC and its immediate predecessor. A qualitative analysis is conducted, identifying all changes made to the previous edition of the IECC, and characterizing these changes in terms of their anticipated impact on residential building energy consumption. A quantitative analysis is then modeled through building energy simulation to estimate the resulting energy impacts. This report documents the technical analysis used to evaluate whether residential buildings constructed to meet the requirements of the 2015 IECC would result in energy efficiency improvements over residential buildings constructed to meet the requirements of the previous edition, the 2012 IECC. PNNL considered all code change proposals approved for inclusion in the 2015 IECC during the ICC code development cycle, and evaluated their combined impact on a suite of prototypical residential building energy models across all U.S. climate zones.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗