Search NASA⌕ Search

SEARCH · Search NASA

Results for “C Code”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Linear-Scaling Asymmetric Triples Correction through the Solution of the DLPNO–CCSD Lambda Equations: DLPNO–CCSD(T) Λ

In this research, we derive equations for solving for the stationary points of the DLPNO–CCSD Lagrangian, in the t 1 -transformed formalism introduced earlier and as currently implemented in the P SI 4 quantum chemistry software package. These lambda equations in the local pair natural orbital basis allow for the evaluation of CCSD(T) Λ energetics with linear-scaling computational effort, also known as the asymmetric triples correction. This DLPNO–CCSD(T) Λ method allows for accurate triples contributions to be computed for larger molecules, especially in cases that CCSD(T) is known to be insufficient, such as with multireference systems and bond-breaking systems. We showcase the accuracy of our code on reaction energies, barrier heights, and noncovalent interaction energies. Also showcased are the capabilities of our code by evaluating DLPNO–CCSD(T) Λ energetics on large noncovalent dimers up to 112 atoms, as well as a rhodium catalyst complex containing 66 atoms.

Cluster chemistry↗
DOE OSTI2020

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗
DOE OSTI2019

Simulation study of CO 2 laser-plasma interactions and self-modulated wakefield acceleration

3D numerical simulations of the interaction of a powerful CO 2 laser with hydrogen jets demonstrating the role of ionization on the characteristics of induced wakes are presented. Simulations using SPACE, a parallel relativistic Particle-in-Cell code, are performed in support of the plasma wakefield accelerator experiments being conducted at the Brookhaven National Laboratory (BNL) Accelerator Test Facility (ATF). A novelty of the SPACE code is its set of efficient atomic physics algorithms that compute ionization and recombination rates on the grid and transfer them to particles. The influence of ionization on the spectrum of pump laser has been studied for a range of gas densities.Simulations reproduce both Stokes and anti-Stokes shifts in the spectrum of pump laser, similar to those observed in experiments in the spectrum of probe. A good agreement has been achieved with the experiments on the effect of variation in gas density on Stokes/anti-Stokes intensity. In addition, self-injection and trapping of electrons into the self-modulated wakes has been observed and analyzed. Here, the experimentally validated code SPACE will be used for predictive simulations to guide future experiments at BNL ATF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Probing ion species separation and ion thermal decoupling in shock-driven implosions using multiple nuclear reaction histories

Simultaneously measured DD, DT, and D 3 He reaction histories are used to probe the impacts of multi-ion physics during the shock phase of inertial confinement fusion implosions. In these relatively hydrodynamiclike (burn-averaged Knudsen number ⟨N K ⟩ ~0.3) shock-driven implosions, average-ion hydrodynamic DUED simulations are able to reasonably match burnwidths, nuclear yields, and ion temperatures. However, kinetic-ion FPION simulations are able to better simulate the timing differences and time-resolved reaction rate ratios between DD, DT, and D 3 He reactions. FPION simulations suggest that the D 3 He/DT reaction rate ratio is most directly impacted by ion species separation between the 3 He and T ions, whereas the D 3 He/DD reaction rate ratio is affected by both ion species separation and ion temperature decoupling effects.

Hydrodynamic codes↗
DOE OSTI2021

Laser transport and backscatter in low-density SiO 2 and Ta2O 5 foams

Experiments using a single 527 nm wavelength beam interacting with sub- and supercritical density SiO 2 and Ta 2 O 5 foams examined laser propagation and backscatter from laser–plasma instabilities such as Stimulated Brillouin Scattering (SBS). Two densities of each material were examined, and multiple diagnostics were used to characterize the propagation and backscatter. For 5 mg/cc SiO 2 (n e /n c = 0.375), the laser propagation distance was well approximated by treating the foam as a gas. However, for the 2 mg/cc SiO 2 foam (n e /n c = 0.15), the same model over-predicts the propagation distance by ~40%. Furthermore, existing analytical theories on propagation through subcritical foams were able to account for this difference. The laser heat wave propagated ~1/2 as far in Ta 2 O 5 than SiO 2 foams with similar electron density. We showed that this difference is due to the increased radiation losses in the higher Z foam. The fraction of backscattered light scales linearly with incident laser intensity for the range of intensities examined. Ta 2 O 5 foams had significantly lower levels of backscatter (1–3%) than the SiO 2 (4–8%), which is consistent with estimates of large Landau damping due to the presence of the oxygen atoms. The measured fraction of SBS backscattered laser energy for a 2 mg/cc SiO 2 foam shot was ~4 times lower than predicted by simulations assuming a gas-like foam. We found that we needed to assume increased ion heating such that T i /T e ~1.2–1.5 in the plasma to agree with the measured SBS reflectivity. Analytical models of laser-heated foams predict preferential heating of the ions as has been observed in previous experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The Thermonuclear Runaway and the Classical Nova Outburst

Nova explosions occur on the white dwarf component of a cataclysmic variable binary stellar system that is accreting matter lost by its companion. When sufficient material has been accreted by the white dwarf, a thermonuclear runaway occurs and ejects material in what is observed as a classical nova explosion. Here, we describe both the recent advances in our understanding of the progress of the outburst and outline some of the puzzles that are still outstanding. We report on the effects of improving both the nuclear reaction rate library and including a modern nuclear reaction network in our one-dimensional, fully implicit, hydrodynamic computer code. In addition, there has been progress in observational studies of supernovae Ia with implications about the progenitors, and we discuss that in this review.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Simulations of tokamak boundary plasma turbulence transport in setting the divertor heat flux width

The BOUT++ code has been used to simulate edge plasma electromagnetic (EM) turbulence and transport, and to study the role of EM turbulence in setting the scrape-off layer (SOL) heat flux width λ q . More than a dozen tokamak discharges from C-Mod, DIII-D, EAST, ITER and CFETR have been simulated with encouraging success. The parallel electron heat fluxes onto the target from the BOUT++ simulations of C-Mod, DIII-D, and EAST follow the experimental heat flux width scaling of the inverse dependence on the poloidal magnetic field. Further turbulence statistics analysis shows that the blobs are generated near the pedestal pressure peak gradient region inside the separatrix and contribute to the transport of the particle and heat in the SOL region. Transport simulations indicate two distinct regimes: drift dominant regime and turbulence dominant regime. Goldston's heuristic drift-based (HD) model yields a consistent divertor heat flux width in the drift dominant regime. For C-Mod enhanced D α H-mode discharges, drifts and turbulence are competing in setting the divertor heat flux width, possibly due to its compact machine size and good pedestal confinement. The simulations for ITER and CFETR indicate that divertor heat flux width λ q of the future machines may no longer follows the 1/B pol,OMP HD-based empirical (Eich) scalings and the HD model gives a pessimistic limit of divertor heat flux width. The simulation results show a transition from a drift dominant regime to a turbulence dominant regime from current machines to future machines such as ITER and CFETR for two reasons. (1) The magnetic drift-based radial transport decreases due to large CFETR and ITER machine sizes. (2) The SOL turbulence thermal diffusivity increases due to larger turbulent fluxes ejected from the pedestal into the SOL when operating in a different pedestal structure, from an ELM-free H-mode pedestal regime to a small and grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Study of shock waves and related phenomena motivated by astrophysics

This study discusses the recent research in High-Energy-Density Physics at our Center. Our work in complex hydrodynamics is now focused on mode coupling in the Richtmyer-Meshkov process and on the supersonic Kelvin-Helmholtz instability. These processes are believed to occur in a wide range of astrophysical circumstances. In radiation hydrodynamics, we are studying radiative reverse shocks relevant to cataclysmic variable stars. Our work on magnetized flows seeks to produce magnetized jets and study their interactions. We build the targets for all these experiments, and simulate them using our CRASH code. We also conduct diagnostic research, focused primarily on imaging x-ray spectroscopy and its applications to scattering and fluorescence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

NuGrid stellar data set – II. Stellar yields from H to Bi for stellar models with MZAMS = 1–25 M⊙ and Z = 0.0001–0.02

We provide here a significant extension of the NuGrid Set 1 models in mass coverage and towards lower metallicity, adopting the same physics assumptions. The combined data set now includes the initial masses M ZAMS /M ⊙ = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with α-enhanced composition for the lowest three metallicities. These models are computed with the MESA stellar evolution code and are evolved up to the AGB, the white dwarf stage, or until core collapse. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid MPPNP code. Explosive nucleosynthesis is based on semi-analytic 1D shock models. Metallicity-dependent mass-loss, convective boundary mixing in low- and intermediate-mass models and H and He core burning massive star models are included. Convective O-C shell mergers in some stellar models lead to the strong production of odd-Z elements P, Cl, K, and Sc. In AGB models with hot dredge-up, the convective boundary mixing efficiency is reduced to accommodate for its energetic feedback. In both low-mass and massive star models at the lowest metallicity, H-ingestion events are observed and lead to i-process nucleosynthesis and substantial 15 N production. Finally, complete yield data tables, derived data products and online analytic data access are provided.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Observation of a magneto-Rayleigh-Taylor instability in magnetically collimated plasma jets

We present the direct experimental observation of the formation of a diamagnetic cavity and magneto-Rayleigh-Taylor (MRT) instability in a β ≈ 1 high energy density plasma. Proton radiography is used to measure the two dimensional path-integrated magnetic field in a laser-produced plasma propagating parallel to a preimposed magnetic field. Flutelike structures, associated with the MRT instability, are observed to grow at the surface of the cavity, with a measured wavelength of 1.2 mm and growth time of 4 ns. These measurements are in good agreement with predictions of three dimensional magnetohydrodynamic simulations using the GORGON code. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Evaluation of Thin Plate Hydrodynamic Stability through a Combined Numerical Modeling and Experimental Effort

An experimental and computational effort was undertaken in order to evaluate the capability of the fluid-structure interaction (FSI) simulation tools to describe the deflection of a Missouri University Research Reactor (MURR) fuel element plate redesigned for conversion to lowenriched uranium (LEU) fuel due to hydrodynamic forces. Experiments involving both flat plates and curved plates were conducted in a water flow test loop located at the University of Missouri (MU), at conditions and geometries that can be related to the MURR LEU fuel element. A wider channel gap on one side of the test plate, and a narrower on the other represent the differences that could be encountered in a MURR element due to allowed fabrication variability. The difference in the channel gaps leads to a pressure differential across the plate, leading to plate deflection. The induced plate deflection the pressure difference induces in the plate was measured at specified locations using a laser measurement technique. High fidelity 3-D simulations of the experiments were performed at MU using the computational fluid dynamics code STAR-CCM+ coupled with the structural mechanics code ABAQUS. Independent simulations of the experiments were performed at Argonne National Laboratory (ANL) using the STAR-CCM+ code and its built-in structural mechanics solver. The simulation results obtained at MU and ANL were compared with the corresponding measured plate deflections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Review of V&V Documentation

Consortium for Advance Simulation of Light Water Reactors (CASL) is a Department of Energy Innovation Hub whose mission is the following, "CASL is a collaboration of the nation's leading scientists, institutions, and supercomputers, with an aggressive 10-year mission to confidently predict the performance of existing and next-generation commercial nuclear reactors through comprehensive, science-based modeling and simulation." The CASL program to date has focused on developing the necessary predictive capability. Rightly so, it is characterized by many as a research project. As a matter of intent, the first 6 years of CASL focused on developing and demonstrating the prediction capability of a suite of independent physics codes: MPACT (neutronics), CTH (thermal hydraulics in the core), BISON (fuel performance), and MAMBA (CRUD and boron uptake on fuel rod surfaces) The last 4 years focused on initial attempts to couple the codes and to demonstrate those capabilities through a series of challenge problems aligned to 3 key issues of interest to the nuclear power industry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Sample problems for Section III, Division 5 design by inelastic analysis of Grade 91 components

This report works through two sample design analyses of representative Grade 91 components using the ASME Boiler & Pressure Vessel Code Section III, Division 5, Appendix HBB-T rules for the strain and deformation and creep-fatigue design criteria using a design by inelastic analysis approach. The component geometries and loadings were selected to span a wide variety of potential applications. The sample components are an endplate geometry, which is fairly representative of typical pressure vessel construction, and a core block, representing a critical component in a heat pipe microreactor, which has a complicated geometry and complex thermal stress history. These sample problems serve several functions. First, they are complete worked examples of how to apply the current design by inelastic analysis rules. Second, they serve as test cases to evaluate the current rules and provide feedback on how they could be optimized to deal with complicated component geometries or to be better automated in modern finite element analysis software. Finally, the analysis here uses an inelastic constitutive model for Grade 91 developed by Argonne National Laboratory and proposed for inclusion in the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules as a reference constitutive model. Using this constitutive model in the sample problem analyses tests the model in the context of the ASME Code rules. The final conclusions of this report are that the material model and the current Code rules are satisfactory, but that the Code rules could be optimized to simplify automation and reduce over conservatism. Future work should address the specific optimization topics discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Annual Report on NCSP Technical Support task in BNL during FY24

This work reports on the investigations done during FY 24 in collaboration with two 2024 Summer undergraduate interns through the DOE SULI and BNL SURP programs. These efforts were partially funded by the Nuclear Criticality Safety Program through the Technical Support Succession Plan task. The first project, developed by Ian Snider, focused on testing the impact of uncertainties in thermal cross sections for multiple materials in nuclear reactor simulations. The second project aimed implement concurrent into BNL’s machine-learning code to correct spin mis-assignments in neutron resonances, the Bayesian Resonance Reclassifier (BRR).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Expected Gamma-Ray Burst Detection Rates and Redshift Distributions for the BlackCAT CubeSat Mission

We report the results of an extensive set of simulations exploring the sensitivity of the BlackCAT CubeSat to long-duration gamma-ray bursts (GRBs). BlackCAT is a NASA APRA-funded CubeSat mission for the detection and real-time subarcminute localization of high-redshift (z ≳ 3.5) GRBs. Thanks to their luminous and long-lived afterglow emissions, GRBs are uniquely valuable probes of high-redshift star-forming galaxies and the intergalactic medium. In addition, each detected GRB with a known redshift serves to localize a region of high-redshift star formation in three dimensions, enabling deep follow-on searches for host galaxies and associated local and large-scale structures. We explore two distinct models for the GRB redshift distribution and luminosity function, both consistent with Swift observations. We find that, for either model, BlackCAT is expected to detect a mean of 42 bursts per year on orbit, with 6.7% to 10% of these at z > 3.5. BlackCAT bursts will be localized to an r 90 ≲ 55'' precision and reported to the community within seconds. Due to the mission orbit and pointing scheme, bursts will be located in the night sky and well placed for deep multiwavelength follow-up observations. BlackCAT is on schedule to achieve launch readiness in 2025.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Residential Energy Efficiency and the Three Rs: Resistance, Resilience, and Recovery

The Building America Solution Center focuses primarily on providing best practices information on energy efficiency, indoor air quality, and water efficiency. While gathering information to address building practices that can improve disaster resistance, resilience, and recovery, it became apparent that there is a significant amount of overlap in areas of home construction related to energy efficiency and in areas related to resiliency. This paper will explore three areas where mutual goals and/or experience may be beneficial to both camps in meeting objectives related to efficient, resistant, resilient, and recoverable homes. The first is the common basis in building science that each industry must draw on to design and build adequate houses. This portion of the paper will include a discussion of photovoltaics. The second is the shared experiences each industry faces in market transformation and achieving changes in building codes. We will explore lessons that these industries may share around economic criteria, incentives, and prescriptive measures in code. And finally, this paper will address the role that insurance companies play in achieving changes in building practices related to energy efficiency and disaster planning.

residential, energy efficiency, disaster resistanc↗
DOE OSTI2018

Numerical Modelling of Radionuclide Migration in the Context of a Near-Surface LILW Disposal Facility - 18618

The reference solution for the management of low and intermediate level short-lived waste (LILW) in France consists in their emplacement in concrete vaults built at the near-surface. The corresponding disposal facility under operation by the French waste management organization (Andra), is located in North-eastern France and planned to be closed in 2060. The safety case developed by Andra has to show the evidence that, amongst other aspects, the radionuclide transfer through the confinement barriers will occur in such a way that the human and environmental exposure will be as low as reasonably achievable (ALARA) and in accordance with the dose constraints. For this purpose, modelling is used to integrate uncertainties about the characteristics of concrete components, waste and water table and their evolution in time and space. At last, modelling is finally used to investigate the radiological impact of specific scenarios in which the failure of safety barriers occurs. In the frame of its technical review of the safety assessment, the French Institute for Radiological Protection and Nuclear Safety (IRSN) also performs such kind of numerical modelling of radionuclide fate from the waste to the hydrogeological outlets, in order to check whether the assessment performed by Andra takes into account sufficiently conservative assumptions. The present work aims to evaluate the role of components of the vaults on radionuclide transport during the post-closure period over a few thousand years. For this purpose, IRSN carries out numerical calculations with MELODIE software, which simulates water flow and transport of species, in order to simulate radionuclide migration from waste packages emplaced in a concrete vault to the surrounding groundwater. The developed transport model focuses on 3-H radionuclide and integrates Andra's main assumptions (release model from the waste packages, chemical and hydraulic properties of the vault components, shrinkage and cracks in concrete components). An adaptive mesh method, based on an a posteriori error estimator, was used to model the cracks and the consideration of a strong heterogeneity in the hydraulic and transport parameters values in the different components of the model. This method shows an improvement of results of about 40 times, whereas the mesh size was increased by only 2%. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗