Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface 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,261 records · Page 70

Stellarator Design Exploration Using Symbolic-Regression Neutronics Surrogates

Systems codes require fast, simplified models to rapidly evaluate fusion power plant concepts, but neutronics analyses are often a computational bottleneck. Here, to address this, surrogate models for key neutronics responses have been developed using 3-D neutronics-ready models built with the open-source code ParaStell from a database of stellarator equilibria. Neutronics responses such as tritium breeding ratio (TBR), nuclear heating, and neutron-induced radiation damage displacements per atom (dpa) were simulated using OpenMC. Through sensitivity analysis and symbolic regression (SR), simple power-law formulas were derived connecting these neutronics responses to global stellarator parameters, including fusion power, plasma surface area, and plasma elongation. Validation shows these formulas can predict the simulation results with low error, enabling quick and accurate assessment of neutronics requirements in stellarator design exploration activities with systems codes.

Modeling↗

NEQAIR96,Nonequilibrium and Equilibrium Radiative Transport and Spectra Program: User's Manual

This document is the User's Manual for a new version of the NEQAIR computer program, NEQAIR96. The program is a line-by-line and a line-of-sight code. It calculates the emission and absorption spectra for atomic and diatomic molecules and the transport of radiation through a nonuniform gas mixture to a surface. The program has been rewritten to make it easy to use, run faster, and include many run-time options that tailor a calculation to the user's requirements. The accuracy and capability have also been improved by including the rotational Hamiltonian matrix formalism for calculating rotational energy levels and Hoenl-London factors for dipole and spin-allowed singlet, doublet, triplet, and quartet transitions. Three sample cases are also included to help the user become familiar with the steps taken to produce a spectrum. A new user interface is included that uses check location, to select run-time options and to enter selected run data, making NEQAIR96 easier to use than the older versions of the code. The ease of its use and the speed of its algorithms make NEQAIR96 a valuable educational code as well as a practical spectroscopic prediction and diagnostic code.

Whiting, Ellis E.↗

Electric-field enhancement caused by subwavelength-sized particles located on the surface of multilayer dielectric mirrors

A laser pulse impinging on the surface of an optical component can interact with particles, such as contamination debris, to produce a scattered electric field, which, either by itself or combined with the incident laser field, coherently can significantly increase the local field intensity. This effect can be of critical importance as it can reduce the laser-induced–damage threshold of the affected component. In this work, we use a field-propagation code to improve understanding regarding the factors that determine the magnitude and location of the electric-field enhancement for the case of subwavelength-sized particles located on the surface of multilayer dielectric mirrors.

47 OTHER INSTRUMENTATION↗

Radiometric Calibration of A Non-Imaging Airborne Spectrometer To Measure the Greenland Ice Sheet Surface

Methods to radiometrically calibrate a non-imaging airborne visible-to-shortwave infrared (VSWIR) spectrometer to measure the Greenland ice sheet surface are presented. Airborne VSWIR measurement performance for bright Greenland ice and dark bare rock/soil targets is compared against the MODerate resolution atmospheric TRANsmission (MODTRAN ® ) radiative transfer code (version 6.0), and a coincident Landsat 8 Operational Land Imager (OLI) acquisition on 29 July 2015 during an in-flight radiometric calibration experiment. Airborne remote sensing flights were carried out in northwestern Greenland in preparation for the Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimeter mission. A total of nine science flights were conducted over the Greenland ice sheet, sea ice, and open-ocean water. The campaign's primary purpose was to correlate green laser pulse penetration into snow and ice with spectroscopic-derived surface properties. An experimental airborne instrument configuration that included a nadir-viewing (looking downward at the surface) non-imaging Analytical Spectral Devices (ASD) Inc. spectrometer that measured upwelling VSWIR (0.35 to 2.5 µm) spectral radiance (Wm −2 sr −1 µm −1 ) in the two-color Slope Imaging Multi-polarization Photon-Counting Lidar's (SIMPL) ground instantaneous field of view, and a zenith-viewing (looking upward at the sky) ASD spectrometer that measured VSWIR spectral irradiance (W m −2 nm −1 ) was flown. National Institute of Standards and Technology (NIST) traceable radiometric calibration procedures for laboratory, in-flight, and field environments are described in detail to achieve a targeted VSWIR measurement requirement of within 5 % to support calibration/validation efforts and remote sensing algorithm development. Our MODTRAN predictions for the 29 July flight line over dark and bright targets indicate that the airborne nadir-viewing spectrometer spectral radiance measurement uncertainty was between 0.6 % and 4.7 % for VSWIR wavelengths (0.4 to 2.0 µm) with atmospheric transmittance greater than 80 %. MODTRAN predictions for Landsat 8 OLI relative spectral response functions suggest that OLI is measuring 6 % to 16 % more top-of-atmosphere (TOA) spectral radiance from the Greenland ice sheet surface than was predicted using apparent reflectance spectra from the nadir-viewing spectrometer. While more investigation is required to convert airborne VSWIR spectral radiance into atmospherically corrected airborne surface reflectance, it is expected that airborne science flight data products will contribute to spectroscopic determination of Greenland ice sheet surface optical properties to improve understanding of their potential influence on ICESat-2 measurements.

Christopher J. Crawford↗

Understanding detachment of the W7-X island divertor

The fundamental behavior of the W7-X island divertor under detached conditions, which has been theoretically predicted with the EMC3-Eirene code, is re-examined here under the experimental conditions achieved so far and compared with the first experimental results. Both simulations and experiments cover a range of divertor configurations and plasma parameters, and show the following common trends: (1) with rising impurity radiation, the target heat load decreases 'uniformly' over the entire target surface in the sense that both the peak and average heat loads can drop by an order of magnitude. Impurity radiation (mainly from intrinsic carbon) occurs primarily at the plasma edge and the resulting negative impact on the stored energy is less than 10%. (2) When the total radiation exceeds a critical level, the target particle flux (the recycling flux Γrecy) begins to fall and can drop by a factor of 3–5 at high radiation levels without an obvious indication of significant volume recombination. (3) While Γ recy decreases, the divertor neutral pressure continues to build up and reaches a maximum, at which point Γ recy has declined significantly. (4) During detachment, the electron temperature at the last closed flux surface falls in a way that is not quantitatively understandable from parallel classical heat conduction processes. This paper presents a physical explanation of the numerical/experimental results described above. Furthermore, using the EMC3-Eirene code as a diagnostic tool, we are able, apparently for the first time, to provide a full quantitative analysis of each transport channel in the island divertor, aiming to clarify how the island divertor plasma self-regulates to maintain particle, energy, and momentum balance under detached conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Upgrade of the Materials Analysis Particle Probe (MAPP-U) to decipher the impact of lithium-based surfaces on NSTX-U plasma behavior

Understanding the plasma wall interaction (PWI) remains a critical issue for the feasibility of thermonuclear magnetic fusion energy solutions. Key issues with PWI mechanisms in fusion tokamak reactors include: evolution of surface chemistry and its role on hydrogen retention. In particular how low-Z coatings such as lithium can impact the behavior of plasma at the edge and in the core. PMI (plasma-material interactions) are particularly important for strategies that involve low-recycling regimes and the use of lithium PFS (plasma-facing surfaces) to attain them, as in the case of NSTX-U. Recent reports have indicated the importance of access to the evolving plasma-facing surface during and in-between plasma discharges. Changes in surface chemistry and morphology due to ion bombardment and the difficulty of diagnosing plasma-facing surfaces, especially reactive surfaces, complicate the development of a predictive understanding of the wall and its interaction with the plasma. Consequently, this impairs the ability to design advanced PFC materials for future plasma-burning fusion reactors and appropriate PMI code validation. The Materials Analysis Particle Probe (MAPP) is an established and on-going PMI probe diagnostic system compatible with the highly chemically reactive system of lithium and boron coatings adopted by the NSTX-U research program. MAPP is the first PMI diagnostic to capture the surface physics and chemistry in-vacuo in a fusion tokamak system and correlate this data to controlled plasma shots. Currently MAPP captures this information at a fixed radial location at the NSTX-U outboard divertor region. The MAPP diagnostic has enabled understanding of the near-surface and surface chemistry of complex evolving lithiated and borated carbon-based PFC surfaces retention and transport of hydrogen. Coupled to atomistic simulations in collaboration with P. Krstic of Stony Brook U. MAPP has been very successful in achieving high-impact scientific research in its current grant period evidenced by two invited review articles and over 20 peer-reviewed manuscripts and over 40 contributed and invited presentations at both national and international conferences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗

Comparative Study of Active Flow Control Strategies for Lift Enhancement of a Simplified High-Lift Configuration

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of the conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-AFC configuration to recover the aerodynamic performance (lift) of the CRM-HL configuration. In the current paper, a Lattice Boltzmann method-based computational fluid dynamics (CFD) code, known as PowerFLOWQ® is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with several different types of AFC devices. The transonic version of the PowerFLOWQ® code that has been validated for high speed flows is used to accurately simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient AFC systems and actuator arrangements emerged based on the parametric studies performed prior to a Fall 2018 wind tunnel test. Preliminary comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations.

Vatsa, Veer N.↗

Impact of T - and ρ -dependent decay rates and new (n, γ ) cross-sections on the s process in low-mass asymptotic giant branch stars

Aims. We study the impact of nuclear input related to weak-decay rates and neutron-capture reactions on predictions for the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars. We provide the first database of surface abundances and stellar yields of the isotopes heavier than iron from the Monash models. Methods. We ran nucleosynthesis calculations with the Monash post-processing code for seven stellar structure evolution models of low-mass AGB stars with three different sets of nuclear inputs. The reference set has constant decay rates and represents the set used in the previous Monash publications. The second set contains the temperature and density dependence of β decays and electron captures based on the default rates of nuclear NETwork GENerator (NETGEN). In the third set, we further update 92 neutron-capture rates based on re-evaluated experimental cross sections from the ASTrophysical Rate and rAw data Library. We compare and discuss the predictions of the sets relative to each other in terms of isotopic surface abundances and total stellar yields. We also compare the results to isotopic ratios measured in presolar stardust silicon carbide (SiC) grains from AGB stars. Results. The new sets of models result in a ∼66% solar s-process contribution to the p-nucleus 152 Gd, confirming that this isotope is predominantly made by the s process. The nuclear input updates result in predictions for the 80 Kr/ 82 Kr ratio in the He intershell and surface 64 Ni/ 58 Ni, 94 Mo/ 96 Mo, and 137 Ba/ 136 Ba ratios that are more consistent with the corresponding ratios measured in stardust; however, the new predicted 138 Ba/ 136 Ba ratios are higher than the typical values of the SiC grains. The W isotopic anomalies are in agreement with data from the analyses of other meteoritic inclusions. We confirm that the production of 176 Lu and 205 Pb is affected by too large uncertainties in their decay rates from NETGEN.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Yield degradation mechanisms for two-shock capsules evaluated through simulations

An investigation of twenty two-shock campaign indirectly driven capsules on the National Ignition Facility was conducted using the xRAGE computer code. The two-shock platform was developed to look at the sensitivity of fuel–ablator mix with shock timing, asymmetry, surface roughness, and convergence on roughly ignition size scale capsules. This platform used CH/CD (plastic/deuterated plastic) shell capsules that were about 685-μm outer radius and filled with D2 or hydrogen-tritium (HT) gas. The experimental radius and velocity vs time, neutron yield, burn averaged ion temperature (Tion), burn width, and self-emission image size were compared to one-dimensional (1D) and two-dimensional (2D) simulations. Our 2D simulations suggest that the mixing of glass from the fill tube was the dominant source of impurity in the gas region of the capsule during burn, along with fuel–ablator mix. The mass of glass mixed in is about 5–10 ng. Our 2D simulations capture most of the yield trends from different degradation mechanisms, and they match the observed burn width and Tion measurements. Our 2D models match all the available data to within 2.5 times the normalized experimental error for 19 of 20 capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current advances on Talbot–Lau x-ray imaging diagnostics for high energy density experiments (invited)

Talbot–Lau x-ray interferometry is a refraction-based diagnostic that can map electron density gradients through phase-contrast methods. The Talbot–Lau x-ray deflectometry (TXD) diagnostics have been deployed in several high energy density experiments. To improve diagnostic performance, a monochromatic TXD was implemented on the Multi-Tera Watt (MTW) laser using 8 keV multilayer mirrors (Δθ/θ = 4.5%-5.6%). Copper foil and wire targets were irradiated at 10 14 –10 15 W/cm 2 . Laser pulse length (~10 to 80 ps) and backlighter target configurations were explored in the context of Moiré fringe contrast and spatial resolution. Foil and wire targets delivered increased contrast <30%. The best spatial resolution (<6 μm) was measured for foils irradiated 80° from the surface. Further TXD diagnostic capability enhancement was achieved through the development of advanced data postprocessing tools. The Talbot Interferometry Analysis (TIA) code enabled x-ray refraction measurements from the MTW monochromatic TXD. Additionally, phase, attenuation, and dark-field maps of an ablating x-pinch load were retrieved through TXD. The images show a dense wire core of ~60 μm diameter surrounded by low-density material of ~40 μm thickness with an outer diameter ratio of ~2.3. Attenuation at 8 keV was measured at ~20% for the dense core and ~10% for the low-density material. Instrumental and experimental limitations for monochromatic TXD diagnostics are presented. Enhanced postprocessing capabilities enabled by TIA are demonstrated in the context of high-intensity laser and pulsed power experimental data analysis. Significant advances in TXD diagnostic capabilities are presented. Furthermore, these results inform future diagnostic technique upgrades that will improve the accuracy of plasma characterization through TXD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗

Advanced beaded and tubular structural panels

A program to develop lightweight beaded and tubular structural panels is described. Applications include external surfaces, where aerodynamically acceptable, and primary structure protected by heat shields. The design configurations were optimized and selected with a computer code which iterates geometric parameters to satisfy strength, stability, and weight constraints. Methods of fabricating these configurations are discussed. Nondestructive testing produced extensive combined compression, shear, and bending test data on local buckling specimens and large panels. The optimized design concepts offer 25 to 30% weight savings compared to conventional stiffened sheet construction.

Musgrove, M. D.↗

Reeds computer code

The REEDS (rocket exhaust effluent diffusion single layer) computer code is used for the estimation of certain rocket exhaust effluent concentrations and dosages and their distributions near the Earth's surface following a rocket launch event. Output from REEDS is used in producing near real time air quality and environmental assessments of the effects of certain potentially harmful effluents, namely HCl, Al2O3, CO, and NO.

Bjork, C.↗

Flow field description for the reaction control system of the Space Shuttle Orbiter

The flow field for the Reaction Control System (RCS) jets of the Space Shuttle Orbiter has been calculated from the nozzle throat to the far-field flow region. The calculations involved the use of recently improved rocket engine nozzle/plume codes and the source flow analytic model. The flow solution is discussed and used to produce impingement forces and moments for surfaces immersed in the flow. The impingement results are compared to calculations involving previous flow solutions in order to evaluate the impingement effects due to the RCS jets during orbital operations.

Alred, J. W.↗

An experimental and computational study of turbulent boundary layers in oscillating flows

An experimental-computational study of the behavior of turbulent boundary layers for oscillating air flows over a plane surface with a small favorable mean pressure gradient is described. Experimental data were obtained for boundary layers generated on the test section wall of a new type of facility that produces a mean flow velocity with a superposed sinusoidal component over a wide frequency range. Experimental results for unsteady velocity profiles and turbulence intensity obtained from hot wire measurements at a fixed stream-wise position on the surface are presented for frequencies ranging from 5 to 21 Hz. Numerical predictions were obtained using an unsteady boundary layer computational code and the Cebeci-Smith and Glushko turbulence models. Predicted quantities related to unsteady velocity profiles exhibit fair agreement with experiment when the Cebeci-Smith model is used.

Cook, W. J.↗