Search NASA⌕ Search

SEARCH · Search NASA

Results for “outer approximation”

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

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Solving Coupled Cluster Equations by the Newton Krylov Method

We describe using the Newton Krylov method to solve the coupled cluster equation. The method uses a Krylov iterative method to compute the Newton correction to the approximate coupled cluster amplitude. The multiplication of the Jacobian with a vector, which is required in each step of a Krylov iterative method such as the Generalized Minimum Residual (GMRES) method, is carried out through a finite difference approximation, and requires an additional residual evaluation. The overall cost of the method is determined by the sum of the inner Krylov and outer Newton iterations. We discuss the termination criterion used for the inner iteration and show how to apply pre-conditioners to accelerate convergence. We will also examine the use of regularization technique to improve the stability of convergence and compare the method with the widely used direct inversion of iterative subspace (DIIS) methods through numerical examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Valence shell photoelectron angular distributions and vibrationally resolved spectra of imidazole: A combined experimental–theoretical study

Linearly polarized synchrotron radiation has been used to record polarization dependent valence shell photoelectron spectra of imidazole in the photon energy range 21–100 eV. These have allowed the photoelectron angular distributions, as characterized by the anisotropy parameter β, and the electronic state intensity branching ratios to be determined. Complementing these experimental data, theoretical photoionization partial cross sections and β-parameters have been calculated for the outer valence shell orbitals. The assignment of the structure appearing in the experimental photoelectron spectra has been guided by vertical ionization energies and spectral intensities calculated by various theoretical methods that incorporate electron correlation and orbital relaxation. Strong orbital relaxation effects have been found for the 15a', nitrogen lone-pair orbital. The calculations also predict that configuration mixing leads to the formation of several low-lying satellite states. The vibrational structure associated with ionization out of a particular orbital has been simulated within the Franck–Condon model using harmonic vibrational modes. The adiabatic approximation appears to be valid for the X 2A" state, with the β-parameter for this state being independent of the level of vibrational excitation. However, for all the other outer valence ionic states, a disparity occurs between the observed and the simulated vibrational structure, and the measured β-parameters are at variance with the behavior expected at the level of the Franck–Condon approximation. Furthermore, these inconsistencies suggest that the excited electronic states may be interacting vibronically such that the nuclear dynamics occur over coupled potential energy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Magnetic mesh generation and field line reconstruction for scrape-off layer and divertor modeling in stellarators

The design of divertor targets and baffles for optimal heat and particle exhaust from magnetically confined fusion plasmas requires a combination of fast, low-fidelity models (such as EMC3-Lite (Feng 2022 Plasma Phys. Control. Fusion 64 125012)) for scoping studies and high-fidelity ones (such as EMC3-EIRENE (Feng et al 2014 Contrib. Plasma Phys. 54 426–31)) for verification. Both of those approaches benefit from a magnetic flux tube mesh for fast interpolation and mapping of field line segments (Feng et al 2005 Phys. Plasmas 12 052505). A new automated mesh generator for unstructured quadrilateral flux tubes with adaptive refinement is presented and integrated into FLARE (Frerichs 2024 Nucl. Fusion 64 106034). For HSX with an extended first wall, it is found that several layers of flux tubes can span the entire half field period before splitting is required. This is an advantage over the traditional setup of the EMC3-EIRENE mesh where careful construction of several sub- domains is required already for the much tighter present first wall. In particular, there is no longer the need to manually construct a suitable outer boundary for the mesh. The divide and conquer paradigm with unstructured mesh layout offers a powerful alternative for fast head load approximation that is suitable for integration into optimization workflows. Further examples for W7-X and CTH demonstrate the versatile application range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling fast ion losses due to tearing and internal kink perturbations in MAST-U

Fast ion (FI) loss properties in the presence of tearing mode and internal kink perturbations are numerically investigated for discharges in the MAST-U spherical tokamak, utilizing the MARS-F magnetohydrodynamic stability code and the REORBIT test particle guiding-center orbit-following module. Here, tracing about 100 000 particle markers sampled from the equilibrium distribution of the neutral-beam injection induced FIs, it is found that about 10% out of the total strike the limiting surface (including the divertor surface) in MAST-U discharge 46943, assuming a maximum perturbation of 100 G inside the plasma (corresponding to ~6 G at the Mirnov probe location at the outboard mid-plane). Detailed particle tracing, assuming a uniform initial distribution in the 2D phase space (at given radial locations), reveals that initially counter-current FIs launched near the plasma edge are subject to significant prompt losses, while almost all initially co-current ions remain well confined at the assumed perturbation level. Most lost FIs strike the lower-half of the limiting surface. Finite gyro-radius effects prevent lost ions from striking the top-outer corner of the super-X divertor chamber. A scan of the perturbation level (based on discharge 45163) reveals, not surprisingly, an approximately linear scaling of the particle loss fraction (for counter-current FIs) with respect to the perturbation amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Destructive Examination of a FeCrAl-UO 2 Irradiation Test

Several destructive postirradiation examinations were performed on an irradiation specimen that coupled an early iron-chromium-aluminum (FeCrAl) candidate alloy cladding with UO2 pellets. This irradiation test was designed to investigate the early life performance and compatibility of the FeCrAl-UO 2 cladding-fuel system under prototypic light water reactor neutronic conditions. Additionally, these tests were expected to provide neutron irradiated samples for severe accident testing of this fuel system. The irradiation studied in this work was part of the ATF-1 series of drop-in style irradiations performed in the Idaho National Laboratory Advanced Test Reactor. The rodlet studied in this work is one of three similar rodlets irradiated in ATF-1 that had approximately 7.6 cm of cladding machined from a wrought FeCrAl alloy and fueled with a 6.1 cm stack of UO 2 pellets. The outer diameter of the cladding was ~0.94 cm and the inner diameter was ~0.83 cm. After irradiation and non-destructive examination in Idaho this rodlet was shipped to the Oak Ridge National Laboratory hot-cells for further examination. The irradiated rodlet was sectioned into several samples for microstructural, micromechanical, and severe accident testing. During sectioning, it was noted that the fuel was not firmly bonded to the cladding and could be readily removed from small cladding slices. Microstructural characterization of fuel cross sections also revealed no significant interaction between the fuel and the cladding. Samples were also prepared for microhardness testing. To prepare for high temperature oxidation testing, the fuel was dissolved from segments of the cladding. High temperature oxidation testing of cladding segments was performed at 1200°C and 1300°C in a steam environment. Comparisons between the oxidation of this FeCrAl alloy in its neutron irradiated state, as-fabricated state, and the oxidation of Zircaloy-2 are made. The oxidation testing will be followed by ring compression testing to evaluate ductility. The microstructure of the samples after oxidation and ring compression testing will also be analyzed.

Harp, Jason↗
DOE OSTI2023

Setting up resistive wall modeling for the GA Compact Advanced Tokamak (GA FPP) (Final Report)

Resistive wall implementation for tokamaks has been achieved by both a 'thin wall' approximation in which time-varying fields are matched across a boundary for all modes; and, by meshing the region from the first wall to the outer boundary and setting certain fields to mimic a resistive wall. Both implementations were tested on existing tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Characterizing the seismic response of a molten salt nuclear reactor

Abstract One Generation IV nuclear reactor, which uses a fluoride salt as a coolant, graphite reflector blocks as a moderator, and circulating buoyant TRISO pebbles as fuel is at an advanced stage of development. To characterize the seismic behavior of components of this reactor, validate numerical models for analysis, and develop recommendations for design, a set of earthquake‐simulator experiments on a scaled model of the reactor vessel and its internals was executed on a six‐degree‐of‐freedom earthquake simulator. The model was seismically isolated at its base using two types of spherical sliding bearings. The scaled model involved representations of the prototype reactor vessel, core barrel, reflector blocks, coolant, and spherical fuel pebbles. The material and geometric properties of different test components were selected based on a dynamic similitude scaling analysis and an approximate length scale of 0.4. Four sets of three‐component earthquake motions were used as inputs for testing. Instrumentation on the test specimen recorded the dynamic responses of the outer vessel, core barrel, and reflector‐block assembly, the hydrodynamic responses (sloshing and hydrodynamic pressure) of the liquid coolant, pebble consolidation under earthquake shaking, and the behavior of the isolation systems. This paper describes the design of the experiments and presents key results from the tests. The dynamic responses of the outer vessel, core barrel, and the reflector blocks revealed that the components responded as a unit for the intense shaking used in the experiments. The sloshing response of the fluid in a thin annulus near the perimeter of the vessel was heavily damped. The change in the packing fraction of the pebble bed under repeated, intense 3D earthquake shaking was less than 3%. Seismically isolating the vessel substantially reduced demands on its internal components.

Engineering↗
DOE OSTI2023

Direct stellarator coil optimization for nested magnetic surfaces with precise quasi-symmetry

We present a robust optimization algorithm for the design of electromagnetic coils that generate vacuum magnetic fields with nested flux surfaces and precise quasi-symmetry. The method is based on a bilevel optimization problem, where the outer coil optimization is constrained by a set of inner least squares optimization problems whose solutions describe magnetic surfaces. The outer optimization objective targets coils that generate a field with nested magnetic surfaces and good quasi-symmetry. The inner optimization problems identify magnetic surfaces when they exist, and approximate surfaces in the presence of magnetic islands or chaos. We show that this formulation can be used to heal islands and chaos, thus producing coils that result in magnetic fields with precise quasi-symmetry. We show that the method can be initialized with coils from the traditional two-stage coil design process, as well as coils from a near-axis expansion optimization. We present a numerical example where island chains are healed and quasi-symmetry is optimized up to surfaces with aspect ratio 6. Another numerical example illustrates that the aspect ratio of nested flux surfaces with optimized quasi-symmetry can be decreased from 6 to approximately 4. In conclusion, the last example shows that our approach is robust and a cold-start using coils from a near-axis expansion optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nucleosynthesis of an 11.8 M ⊙ Supernova with 3D Simulation of the Inner Ejecta: Overall Yields and Implications for Short-lived Radionuclides in the Early Solar System

Based on a 3D supernova simulation of an 11.8 M ⊙ progenitor model with an initial solar composition, we study the nucleosynthesis using tracers covering the innermost 0.1 M ⊙ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of 45 Sc and 64 Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emission from the proto-neutron star. The yields of these isotopes are rather uncertain due to the approximate neutrino transport used in the simulation. In order to obtain the total yields for the whole supernova, we combine the results from the tracers with those for the outer layers from a suitable 1D model. Using the yields of short-lived radionuclides (SLRs), we explore the possibility that an 11.8 M ⊙ supernova might have triggered the formation of the solar system and provided some of the SLRs measured in meteorites. In particular, we discuss two new scenarios that can account for at least the data on 41 Ca, 53 Mn, and 60 Fe without exceeding those on the other SLRs.

Astronomy & Astrophysics↗
DOE OSTI2020

Transportation analysis and related design optimization of the Fermilab high-beta 650 MHz cryomodule

The Proton Improvement Plan-II (PIP-II) at Fermi National Accelerator Laboratory (FNAL) will create a new and vastly improved accelerator, which will be the source of high-energy particles for the experiments taking place at FNAL. The new linear accelerator (LINAC) contains several types of cryomodules, which are individual particle accelerators. The last cryomodule in the LINAC will be the High-Beta 650 (HB650), which will operate at 650 MHz. Each module is approximately 15 meters in length and 1.5 meters wide, weighs 13 tonnes, and shares many design features with three of the other cryomodules. The HB650 consists of two primary sections, the outer vacuum vessel and the components that reside within it, which will be cooled to cryogenic temperatures – the cold mass. Once assembly is complete at FNAL, it will be transported to one of three places: another location at FNAL, another national laboratory within the continental United States, or to a scientific partner in Europe. Any excitation the module experiences, such as going over rough roads when on a semitrailer, can create high stresses in components and cause failure if severe enough. Additionally, any delicate components with low resonant frequencies that match the excitation spectrum could achieve resonance, potentially causing a fatigue failure by repeated flexing. To ensure the successful transportation of the HB650, a transport analysis utilizing ANSYS simulation software has been performed on major subassemblies, as well as analysis of the subassemblies as a combined system. The optimization of components, design of support structures, and overall increasing of the lowest resonant frequencies have resulted in a cryomodule that will be better equipped to handle transportation and any large impact loads that come with it.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

The HOSTS Survey: Evidence for an Extended Dust Disk and Constraints on the Presence of Giant Planets in the Habitable Zone of β Leo

The young (50–400 Myr) A3V star β Leo is a primary target to study the formation history and evolution of extrasolar planetary systems as one of the few stars with known hot (∼1600 K), warm (∼600 K), and cold (∼120 K) dust belt components. In this paper, we present deep mid-infrared measurements of the warm dust brightness obtained with the Large Binocular Telescope Interferometer (LBTI) as part of its exozodiacal dust survey (HOSTS). The measured excess is 0.47% ± 0.050% within the central 1.5 au, rising to 0.81% ± 0.026% within 4.5 au, outside the habitable zone of β Leo. This dust level is 50 ± 10 times greater than in the solar system’s zodiacal cloud. Poynting–Robertson drag on the cold dust detected by Spitzer, and Herschel underpredicts the dust present in the habitable zone of β Leo, suggesting an additional delivery mechanism (e.g., comets) or an additional belt at ∼5.5 au. A model of these dust components is provided that implies the absence of planets more than a few Saturn masses between ∼5 au and the outer belt at ∼40 au. We also observationally constrain giant planets with the LBTI imaging channel at 3.8 μm wavelength. Assuming an age of 50 Myr, any planet in the system between approximately 5–50 au must be less than a few Jupiter masses, consistent with our dust model. Taken together, these observations showcase the deep contrasts and detection capabilities attainable by the LBTI for both warm exozodiacal dust and giant exoplanets in or near the habitable zone of nearby stars.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Surplus Plutonium Disposition Sphincter Seal Development & Testing

The Surplus Plutonium Disposition Project (SPD, Project Y744) requested the Savannah River National Laboratory (SRNL) to perform a series of leak and durability tests on sphincter seals that were designed to process blend cans and shield cans into an SPD glovebox. Each sphincter seal assembly consisted of a stainless-steel outer sleeve weldment that was bolted onto the test enclosure and an aluminum inner sleeve that housed seven to nine Neoprene seals. All tests were executed to determine long term effects from processing approximately a year’s worth of blend and shield cans through the sphincter seals.

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

Calibration of the CMS hadron calorimeters using proton-proton collision data at √ s = 13 TeV

Methods are presented for calibrating the hadron calorimeter system of the CMS detector at the LHC. The hadron calorimeters of the CMS experiment are sampling calorimeters of brass and scintillator, and are in the form of one central detector and two endcaps. These calorimeters cover pseudorapidities |η| < 3 and are positioned inside the solenoidal magnet. An outer calorimeter, outside the magnet coil, covers |η| < 1.26, and a steel and quartz-fiber Cherenkov forward calorimeter extends the coverage to |η| < 5.19. The initial calibration of the calorimeters was based on results from test beams, augmented with the use of radioactive sources and lasers. The calibration was improved substantially using proton-proton collision data collected at √ s = 7 , 8, and 13 TeV, as well as cosmic ray muon data collected during the periods when the LHC beams were not present. The present calibration is performed using the 13 TeV data collected during 2016 corresponding to an integrated luminosity of 35.9 fb -1 . The intercalibration of channels exploits the approximate uniformity of energy collection over the azimuthal angle. The absolute energy scale of the central and endcap calorimeters is set using isolated charged hadrons. The energy scale for the electromagnetic portion of the forward calorimeters is set using Z→ ee data. The energy scale of the outer calorimeters has been determined with test beam data and is confirmed through data with high transverse momentum jets. In this paper, we present the details of the calibration methods and accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Comparison of reduced model predictions for divertor detachment onset and reattachment timescales in ASDEX Upgrade and JET experiments

Building on prior analysis of ASDEX Upgrade (AUG) experiments (Henderson et al 2023 Nucl. Fusion63 086024), this study compares simple analytical formula predictions for divertor detachment onset and reattachment timescales in JET experiments. Detachment onset primarily scales with divertor neutral pressure, impurity concentration, power directed to the targets, machine size, and integral perpendicular power decay length. JET experiments, focusing on seeding mixtures of Ne and Ar, align with the detachment onset predictions. Radiation efficiencies among the impurities show good agreement with the model predictions, contrasting with AUG observations which suggested higher efficiency for Ar and lower efficiency for Ne. The time taken to re-ionise the neutral volume in front of the outer target in fully detached divertor conditions was measured following both abrupt increases in injected neutral beam power and, separately, cutting of the impurity gas flow. Re-ionisation of the neutrals occurs within approximately 1 s on JET, which aligns with the simple model prediction derived from AUG data. While the AUG results are not new, their comparison with the JET results enhances understanding, reinforcing confidence in using simple models to predict future reactor scenarios.

JET↗
DOE OSTI2022

Machine-Learning Accelerated Studies of Materials with High Performance and Edge Computing

In the studies of materials, experimental measurements often serve as the reference to verify physics theory and modeling; while theory and modeling provide a fundamental understanding of the physics and principles behind. However, the interactions and cross validation between them have long been a challenge even to-date. Not only that inferring a physics model from experimental data is itself a difficult inverse problem, another major challenge is the orders-of-magnitude longer wall-clock time required to carry out high-fidelity computer modeling to match the timescale of experiments. We envisage that by combining high performance computing, data science, and edge computing technology, the current predicament can be alleviated, and a new paradigm of data-driven physics research will open up. For example, we can accelerate computer simulations by first performing the large-scale modeling on high performance computers and train a machine-learned surrogate model. This computationally inexpensive surrogate model can then be transferred to the computing units residing closely to the experimental facilities to perform high-fidelity simulations at a much higher throughout. The model will also be more amenable to analyzing and validating experimental observations in comparable time scales at a much lower computational cost. Further integration of these accelerated computer simulations with an outer machine learning loop can also inform and direct future experiments, while making the inverse problem of physics model inference more tractable. We will demonstrate a proof-of-concept by using a quantum Monte Carlo application, Dynamical Cluster Approximation (DCA++), to machine-learn a surrogate model and accelerate the study of quantum correlated materials.

Li, Ying Wai↗
DOE OSTI2023

Validation of a detonation product equation of state for an insensitive high explosive via slab geometry expansion tests

Slab expansion (SLABEX) tests are conducted to validate a process for calibrating the detonation products equation of state (EOS) of a high explosive (HE). The SLABEX tests use rectangular slabs of PBX 9502, a polymer-bonded HE formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene bound with Kel F-800, a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Three PBX 9502 slab thicknesses are examined, each confined symmetrically by two rectangular copper (Cu) plates approximately one-tenth the thickness of the HE slab. For the duration of each experiment, the detonation flow along the central axis of the PBX 9502 slab remains two-dimensional. The lateral flow velocity component of the outer surfaces of the expanding Cu plates is measured, along with the steady axial detonation speed along the central axis of the SLABEX. Hydrodynamic simulations of the Cu plate expansion in the SLABEX geometry, driven by the energy stored in the detonation products by the detonation combustion event, are conducted using a Jones–Wilkins–Lee EOS for the detonation products. This EOS form was recently parameterized for PBX 9502 in the cylinder expansion test geometry using a newly developed calibration technique [Voelkel et al., Combust. Flame 246, 112373 (2022)]. Good agreement between the experiment and prediction is found in each SLABEX test, demonstrating that the detonation product EOS calibration technique produces EOSs that are predictive when applied to other geometries.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Deep Neural Network Assisted Distributed Strain and Temperature Fiber Sensor System for Natural Gas Pipeline Monitoring

Natural gas pipeline integrity monitoring is crucial to detect potential leaks, find structural issues, and prevent environmental damage. This article presents a system of natural gas pipeline monitoring that uses a specialized double Brillouin peak sensing fiber along with the Brillouin optical time domain analysis (BOTDAs) technique. The calibrated sensing fiber coefficients for strain and temperature are 41.8 kHz/ με and 0.9 MHz/°C for peak 1; and 47.2 kHz/ με , and 1.11 MHz/°C for peak 2, respectively. Initially, lab tests were performed by installing a short section of double Brillouin peak fiber (DBPF) on a 1-in steel pipe under pressure up to 1000 per square inch (psi) at elevated temperatures. Simultaneous distributed measurements of temperature and pressure-induced hoop strain were successfully measured. Considering the long processing speed to extract Brillouin frequency shift (BFS), we employ a novel probabilistic deep neural network (PDNN) framework for rapid BFS prediction. Additionally, using the Finite Element Method, the effects of the pipeline pressure on hoop strain were modeled and compared to the experimental hoop strain under the same set of pipeline conditions. Finally, an actual 4-in outer diameter steel natural gas pipeline was used for pilot-scale tests, where hoop strain was measured at various pressure levels. Leaks were simulated to demonstrate accurate pipeline integrity monitoring. At an internal pipe pressure of 1000 psi, hoop strain of approximately 300 με was observed, and the sensitivity was calculated as 0.28 με /psi. The results of this pilot-scale study demonstrated that the system is capable of performing distributed monitoring sufficient to detect pipeline pressure and the presence of leaks to ensure the safe operation of gas pipelines in the field.

03 NATURAL GAS↗
DOE OSTI2020

The breakBRD Breakdown: Using IllustrisTNG to Track the Quenching of an Observationally Motivated Sample of Centrally Star-forming Galaxies

The observed breakBRD ("break bulges in red disks") galaxies are a nearby sample of face-on disk galaxies with particularly centrally concentrated star formation: they have red disks but recent star formation in their centers as measured by the D n 4000 spectral index. In this paper, we search for breakBRD analogs in the IllustrisTNG simulation and describe their history and future. We find that a small fraction (~4% at z = 0; ~1% at z = 0.5) of galaxies fulfill the breakBRD criteria, in agreement with observations. In comparison with the mass-weighted parent IllustrisTNG sample, these galaxies tend to consist of a higher fraction of satellite and splashback galaxies. However, the central, non-splashback breakBRD galaxies show similar environments, black hole masses, and merger rates, indicating that there is not a single formation trigger for inner star formation and outer quenching. We determine that breakBRD analog galaxies as a whole are in the process of quenching. The breakBRD state, with its highly centrally concentrated star formation, is uncommon in the history of either currently quiescent or star-forming galaxies; however, approximately 10% of 10 < M*/M⊙ < 10 11 quiescent galaxies at z = 0 have experienced SFR concentrations comparable to those of the breakBRDs in their past. Additionally, the breakBRD state is short lived, lasting a few hundred Myr up to ~2 Gyr. The observed breakBRD galaxies may therefore be a unique sample of outside-in quenching galaxies.

Astronomy & Astrophysics↗