SEARCH · Engineering Papers

Results for “flux analysis”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗

NASA NTRS1979

Neutron-gamma techniques for planetary exploration

Chemical analysis of planetary surfaces is necessary for the understanding of the origin and evolutionary processes of the solar system. Orbital lunar geochemical experiments performed during the Apollo 15 and 16 demonstrated that even with a low neutron flux (fast and thermal), reliable results could be extracted. The possibility of using a compact 14 MeV neutron generator for geochemical analysis of planetary surfaces, comets and asteroids is currently studied. This method allows the determination of bulk chemical composition, even in the presence of an atmosphere. This would be possible on the surface of Venus, for example, where alternate methods such as sample return are impractical. This method can be used for continuous monitoring of elemental abundances from a roving vehicle, for example, on the surface of Mars.

Johnson, R. G.↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

NASA NTRSDate not supplied

Radiative Flux Measurements from ARISE: A Comparison with CERES Top-of-Atmosphere Radiative Fluxes

Uncertainty in top-of-atmosphere (TOA) radiation fluxes observations are larger in the Arctic than in other regions. These uncertainties are due to the low sun angles and the highly reflective, anisotropic, and heterogeneous surface conditions. Quantifying, attributing, and reducing Arctic TOA radiative flux uncertainty enables a better understanding of the rapidly changing Arctic. To advance this goal, we compare the Cloud and Earth’s Radiant Energy System (CERES) TOA radiative fluxes with Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) campaign measurements collected in September 2014. We compare CERES TOA and aircraft radiative flux measurements using two complementary approaches: grid box average fluxes and instantaneously matched footprints. The grid box mean flux comparison indicates an agreement between CERES and aircraft measurements within 2 uncertainty (calibration and inversion) in the longwave for all five grid boxes and for four-of-five grid boxes in the shortwave; shortwave and longwave mean differences are -7.9 and +2.3 Wm‑2, respectively. The comparison of 36 instantaneously matched footprints with aircraft measurements reveals mean differences of -10.5 and 0.4 Wm‑2 in the shortwave and longwave, respectively. To further explore the persistent negative difference in the shortwave, we further quantify the effects of temporal and spatial sampling differences, angular distribution models, and scene identification to CERES-aircraft differences. Our analysis indicates that sampling differences (including scene evolution) account for an additional 1.8 and 1.7% uncertainty in the shortwave and longwave, respectively and indicates no bias. After accounting for this sampling uncertainty, all CERES-aircraft grid box mean fluxes agree within 2 uncertainty. Scene identification errors due to sea ice concentration data set differences exhibit no bias in the shortwave flux difference and indicate the possibility of substantial differences in the CERES fluxes in specific cases with large spatial heterogeneity. Considering the instantaneously matched footprints, we find that the angular distribution models account may account for up to ‑7.3 Wm-2 of the persistent CERES-aircraft shortwave flux difference due to systematic differences in the anisotropy for sea ice partly cloudy scenes. Additional analysis using a special scan model with the CERES FM2 instrument suggests a significant view zenith angle dependence of the CERES fluxes for sea ice partly cloudy scenes where shortwave fluxes systematically decrease with increasing view zenith angle; no dependence is found for other scene types. We conclude that (1) spatial heterogeneity and scene temporal evolution substantially limit our ability to use aircraft measurements to place strong constraints on CERES TOA fluxes and (2) that the representation of anisotropy in sea ice partly cloudy scenes is a significant factor contributing to the persistent negative CERES-aircraft shortwave flux difference in this comparison and require additional data to analysis fully quantify the potential bias.

Patrick C Taylor↗

NASA NTRS2004

MARIE Dose and Flux Measurements in Mars Orbit

We present results from the Martian Radiation Environment Experiment (MARIE), aboard the 2001 Mars Odyssey spacecraft in orbit around Mars. MARIE operated successfully from March 2002 through October 2003. At the time of this writing, the instrument is off due to a loss of communications during an extremely intense Solar Particle Event. Efforts to revive MARIE are planned for Spring 2004, when Odyssey's role as a communications relay for the MER rovers is completed. During the period of successful operation, MARIE returned the first detailed energetic charged particle data from Mars. Due to limitations of the instrument, normalizing MARIE data to flux or dose is not straightforward - several large corrections are needed. Thus normalized results (like dose or flux) have large uncertainties and/or significant model-dependence. The problems in normalization are mainly due to inefficiency in detecting high-energy protons (signal-to-noise problems force the trigger threshold to be higher than optimal), to the excessively high gains employed in the signal processing electronics (many ions deposit energy sufficient to saturate the electronics, and dE/dx information is lost), and to artifacts associated with the two trigger detectors (incomplete registration of dE/dx). Despite these problems, MARIE is efficient for detecting helium ions with kinetic energies above about 30 MeV/nucleon, and for detecting high-energy ions (energies above about 400 MeV/nucleon) with charges from 5 to 10. Fluxes of these heavier ions can be compared to fluxes obtained from the ACE/CRIS instrument, providing at least one area of direct comparison between data obtained at Earth and at Mars; this analysis will be presented as a work in progress. We will also present dose-rate data, with a detailed explanation of the many sources of uncertainty in normalization. The results for both flux and dose will be compared to predictions of the HZETRN model of the GCR.

Zeitlin, C.↗

DOE OSTI2022

Two-step neutronics calculations with Shift and Griffin for advanced reactor systems

This research develops the initial coupling of the Shift Monte Carlo (MC) code and the Griffin reactor physics code for reactor analysis of non–light-water reactor systems. The novelty of this work is twofold. It is the first application of Shift to produce the multigroup cross sections needed for Griffin as applied to a non–light-water reactor system; and, the first investigation and analysis of characteristics of the Empire microreactor benchmark that should be considered for steady state and transient reactor physics calculations. This application uses the previously developed two-step neutronics analysis workflow to demonstrate this initial coupling. Here, we outline the two-step neutronics analysis workflow in which the Shift MC code is used to generate the multigroup cross sections and fluxes needed by the Griffin deterministic solver. Details on how these multigroup cross sections are generated using MC tallies are given, as well as the practicalities and limitations of the two-step neutronics workflow. The Empire microreactor benchmark was used to investigate and validate this coupling. Results using this benchmark show good agreement between Griffin calculations using Serpent-generated cross sections and Shift-generated cross sections. Analysis of the characteristics of this Empire benchmark show larger eigenvalue differences between heterogeneous and pin–homogenized solutions compared to those of traditional light-water reactor (LWR) designs, thus requiring super homogenization factor corrections for accurate eigenvalue and power distribution predictions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1990

Sounding rocket measurement of the absolute solar EUV flux utilizing a silicon photodiode

A newly developed stable and high quantum efficiency silicon photodiode was used to obtain an accurate measurement of the integrated absolute magnitude of the solar extreme UV photon flux in the spectral region between 50 and 800 A. The adjusted daily 10.7-cm solar radio flux and sunspot number were 168.4 and 121, respectively. The unattenuated absolute value of the solar EUV flux at 1 AU in the specified wavelength region was 6.81 x 10 to the 10th photons/sq cm per s. Based on a nominal probable error of 7 percent for National Institute of Standards and Technology detector efficiency measurements in the 50- to 500-A region (5 percent on longer wavelength measurements between 500 and 1216 A), and based on experimental errors associated with the present rocket instrumentation and analysis, a conservative total error estimate of about 14 percent is assigned to the absolute integral solar flux obtained.

Ogawa, H. S.↗

DOE OSTI2024

Evaluation of Additively Manufactured Monolithic SiC and SiC-Matrix Ceramic Matrix Composites for Concentrating Solar Receiver Applications - Fabrication and Testing of Receiver Design Feature Specimens in a Simulating Lab Test via Laser Heating

Increasing operating temperatures of solar receivers is paramount to the efficiency of concentrated solar thermal (CST) and solar power (CSP) systems. Owing to its high temperature stability combined with excellent thermal and optical properties, SiC has been the material of choice for application in high-temperature solar receivers. We report the results of our study of the effective heat transfer characteristics of several candidate SiC structure motifs, or feature geometries, which are fabricated via additive manufacturing. The SiC structure motifs studied include different permutations of three-dimensional periodic lattices and defined shapes. A solar-thermal simulating laboratory test setup is constructed using a 4kW CO2 laser system with beam shaping optics to apply concurrent radiative heating power on one face of 2”-diameter cylindrical feature specimens, representing the structure motifs of interest for receiver element design, while flowing through the sample as heat transfer fluid. Using the test setup, simulative test conditions representative of a concentrated solar flux of up to ~2000 suns could be achieved in the lab tests under varying air flow through the test structure. A simple 1D numerical analysis scheme is developed to extract an effective or compound heat transfer coefficient representative of the test structure under steady-state heat flow conditions. The test results and their use to guide the selection and optimization of SiC material and structure motifs for the receiver element design fabrication are discussed.

14 SOLAR ENERGY↗

NASA NTRS1991

Evidence of chaotic pattern in solar flux through a reproducible sequence of period-doubling-type bifurcations

A preliminary study of the limits to solar flux intensity prediction, and of whether the general lack of predictability in the solar flux arises from the nonlinear chaotic nature of the Sun's physical activity is presented. Statistical analysis of a chaotic signal can extract only its most gross features, and detailed physical models fail, since even the simplest equations of motion for a nonlinear system can exhibit chaotic behavior. A recent theory by Feigenbaum suggests that nonlinear systems that can be led into chaotic behavior through a sequence of period-doubling bifurcations will exhibit a universal behavior. As the control parameter is increased, the bifurcation points occur in such a way that a proper ratio of these will approach the universal Feigenbaum number. Experimental evidence supporting the applicability of the Feigenbaum scenario to solar flux data is sparse. However, given the hypothesis that the Sun's convection zones are similar to a Rayleigh-Bernard mechanism, we can learn a great deal from the remarkable agreement observed between the prediction by theory (period doubling - a universal route to chaos) and the amplitude decrease of the signal's regular subharmonics. It is shown that period-doubling-type bifurcation is a possible route to a chaotic pattern of solar flux that is distinguishable from the logarithm of its power spectral density. This conclusion is the first positive step toward a reformulation of solar flux by a nonlinear chaotic approach. The ultimate goal of this research is to be able to predict an estimate of the upper and lower bounds for solar flux within its predictable zones. Naturally, it is an important task to identify the time horizons beyond which predictability becomes incompatible with computability.

Ashrafi, S.↗

NASA NTRSDate not supplied

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

NASA NTRSDate not supplied

A comparison between CERES TOA Radiative Fluxes and Airborne Radiative Flux Measurements from ARISE

Uncertainty in top-of-atmosphere (TOA) radiation fluxes observations are larger in the Arctic than in other regions. These uncertainties are due to the low sun angles and the highly reflective, anisotropic, and heterogeneous surface conditions. Quantifying, attributing, and reducing Arctic TOA radiative flux uncertainty enables a better understanding of the rapidly changing Arctic. To advance this goal, we compare the Cloud and Earth’s Radiant Energy System (CERES) TOA radiative fluxes with Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) campaign measurements collected in September 2014. We compare CERES TOA and aircraft radiative flux measurements using two complementary approaches: grid box average fluxes and instantaneously matched footprints. The grid box mean flux comparison indicates an agreement between CERES and aircraft measurements within 2σ uncertainty (calibration and inversion) in the longwave for all five grid boxes and for four-of-five grid boxes in the shortwave; shortwave and longwave mean differences are -7.9 and +2.3 Wm 2, respectively. The comparison of 36 40 instantaneously matched footprints with aircraft measurements reveals mean differences of -120.25 and -1.00.4 Wm 2 in the shortwave and longwave, respectively. To further explore the persistent negative difference in the shortwave, Wwe further quantify the effects of temporal and spatial sampling differences, angular distribution models, and scene identification to CERES-aircraft differences. Our analysis indicates that sampling differences (including scene evolution) account for an additional 1.8 and 1.7% uncertainty in the shortwave and longwave, respectively and , but indicates no bias. After accounting for this sampling uncertainty, all CERES-aircraft grid box mean fluxes agree within 2σ uncertainty. Scene identification errors due to sea ice concentration data set differences exhibit no bias in the shortwave flux difference and indicate the possibility of substantial differences in the CERES fluxes in specific cases with large spatial heterogeneity. Considering the instantaneously matched footprints, we find that the angular distribution models account may account for up to 7.3 Wm-2 of the persistent CERES-aircraft shortwave flux difference due to systematic differences in the anisotropy for sea ice partly cloudy scenes. Additional analysis using a special programmable scan model with the CERES FM2 instrument suggests a significant view zenith angle dependence of the CERES fluxes for sea ice partly cloudy scenes where shortwave fluxes systematically decrease with increasing view zenith angle; no dependence is found for other scene types. We conclude that (1) spatial heterogeneity and scene temporal evolution substantially limit our ability to use aircraft measurements to place strong constraints on CERES TOA fluxes and (2) that the representation of anisotropy in sea ice partly cloudy scenes is likely a significant factor contributing to the persistent negative CERES-aircraft shortwave flux difference in this comparison and require additional data to analysis fully quantify the potential bias.

Patrick C. Taylor↗

NASA NTRS1993

A study of the Merritt Island, Florida sea breeze flow regimes and their effect on surface heat and moisture fluxes

Data collected during the Convective and Precipitation/Electrification Experiment were analyzed as part of an investigation of the sea breeze in the vicinity of Merritt Island, Florida. Analysis of near-surface divergence fields shows that the classical 24-hour oscillation in divergence over the island due to the direct sea breeze circulation is frequently disrupted and exhibits two distinct modes: the classical sea breeze pattern and deviations from that pattern. A comparison of clear day surface energy fluxes with fluxes on other days indicates that changes in magnitudes were dominated by the presence or absence of clouds. Non-classical sea breeze days tended to lose more available energy in the morning than classical sea breeze days due to earlier development of small cumulus over the island. A composite storm of surface winds, surface energy fluxes, rainfall, and satellite visible data was constructed. A spectral transmittance over the visible wavelengths for the cloud cover resulting from the composite storm was calculated. It is shown that pre-storm transmittances of 0.8 fall to values near 0.1 as the downdraft moves directly over the site. It is also found that under post-composite storm conditions of continuous clear sky days, 3.5 days are required to evaporate back into the atmosphere the latent heat energy lost to the surface by rainfall.

Rubes, M. T.↗

DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Analysis of WIRE-21 SPND and Optical Fiber Sensor Measurements

The Wireless Instrumented RB Experiment 2021 (WIRE-21) was a highly instrumented experiment designed to test prototype wireless temperature and pressure sensors developed by Westinghouse Electric Corporation (WEC) in the High Flux Isotope Reactor (HFIR). In addition to these sensors, a suite of other sensors was also integrated into the experiment capsule, including various types of distributed optical fiber sensors and self-powered neutron detectors (SPNDs). This report presents an analysis of the data generated by the SPNDs and optical fiber sensors that were tested under the highest reported neutron flux. The distributed optical fiber sensor results demonstrate that F-doped optical fibers, particularly those inscribed with fs fiber Bragg gratings (FBGs), are capable of surviving fast neutron fluences on the order of 1021 n fast /cm 2 at temperatures relevant to light-water reactors (between 200 and 400°C). However, a significant blue-shift in the optical spectra of these sensors was observed over the course of irradiation which cannot be explained based on the current understanding of radiation-induced compaction in fused silica glass. A mechanistic understanding of this drift has not yet been developed and is proposed as future scope under the Advanced Sensors and Instrumentation program. Of the four SPNDs, only one appeared to operate normally during three cycles of irradiation. A radiation transport model of the experiment in HFIR was used to determine time-dependent neutron flux in this SPND and to calculate the neutron sensitivity of the device. Results showed a cycle averaged sensitivity of 1.5 ×10 -22 and 1.4 ×10 -22 A/nν for the first and third cycles of irradiation, respectively. Additional details pertaining to the optical fibers and SPNDs are included herein.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Experimental evidence of enhanced radial transport in small ELM regimes at DIII-D

Small/type-II edge-localized-modes (ELMs), carrying 1% of the plasma stored energy, are found to deposit only 45 ± 5% of the ELM power near the strike point, and the remaining 55 ± 5% to the far scrape-off-layer (SOL). Small ELMs spread their power over a larger area compared to type-I ELMs, where such a ratio is about 60% and 40% for near- and far-SOL regions, respectively. The larger spread is reflected in the heat flux width (λq) in the SOL for the intra-small ELMs profile of 6.0 mm, almost a factor 2 larger than that of type-I ELMs of 3.15 mm, for similar plasma conditions and magnetic configuration. At the ELM peak, the small ELMs λq is found to be up to 4 times larger than for the type-I ELMs, going from 2 to 7.9 mm, indicating enhanced radial transport in the neon-seeded small ELM scenario. Inter-ELM λqs have been also calculated at the secondary outer divertor in quasi-double-null (QDN) discharges. It is found that, on average, λq is 2.2 times larger in the high-separatrix-density small ELM regime, compared to a reference type-I ELM one. These findings are supported by small ELMs radial velocity profiles, measured at the outer midplane with a fast reciprocating probe, showing a decay length (λvr) in the SOL of 12.8 cm, which is 3.3 times larger than that for the type-I ELMs of 3.9 cm. This analysis shows that small ELMs, although attractive for future machines due to low peak heat flux and large λq, might be of concern for the larger flux to the outer wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Error analysis of numerical methods for thick diffusive neutron transport problems on Shishkin mesh

A thin layer will develop at the boundary if the incoming angular flux is anisotropic in thick diffusive neutron transport problems. Solving such singularly perturbed problems, which have non-smooth solutions with singularity near the boundary, is computationally challenging. Standard finite difference schemes on a uniform mesh cannot yield ε-uniform convergence, where ε is a small parameter, while it can be achieved on a suitable piecewise-uniform Shishkin mesh. We present a formal error analysis of the diamond difference (DD) method and step difference (SD) method for solving the S{sub N} neutron transport equation. The analysis can be extended to other finite difference methods. Numerical results are presented to confirm the error estimates and the advantages of the Shishkin mesh. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Assess and Benchmark Magneto-Inertial Fusion (MIF) Scaling. Final Technical Report for the SNL/Rochester ALPHA Follow-on Project

This project was a follow-on to the Sandia National Laboratories (SNL) and the Laboratory for Laser Energetics (LLE) ARPA-E ALPHA project entitled “Demonstrating Fuel Magnetization and Laser Heating Tools for Low-Cost Fusion Energy”. The primary purpose of this follow-on project was to obtain additional data at the OMEGA facility to help better understand how MagLIF, a platform that has already demonstrated the scientific viability of magneto-inertial fusion, scales across a factor of 1000 in driver energy. A secondary aspect of this project was to extend simulations and analysis at SNL to cover a wider magneto-inertial fusion (MIF) parameter space and test scaling of those models across this wide range of input energies and conditions of the target. This work was successful in improving understanding of how key physics elements of MIF scales and improves confidence in setting requirements for fusion gain with larger drivers. The OMEGA experiments at the smaller scale verified the hypothesis that preheating the fuel plays a significant role in introducing wall contaminants that mix into the fuel and significantly degrade fusion performance. This contamination not only impacts target performance but the optimal input conditions for the target. However, analysis at the Z-scale showed that target performance at high preheat levels is limited by the Nernst effect, which advects magnetic flux from the hot spot, reducing magnetic insulation and consequently reduces the temperature of the fuel. The combination of MagLIF experiments at the disparate scales of OMEGA and Z along with a multiscale 3D simulation analysis has led to new insight into the physical mechanisms responsible for limiting target performance and provides important benchmarks to assess target scaling more generally for MIF schemes. Finally, in addition to the MagLIF related work, a semi-analytic model of liner driven Field Reversed Configuration (FRC) was developed that predicts the fusion gain for such systems. This model was also validated with 2D radiation magneto-hydrodynamic simulations and predicts that fusion gains of near unity could be driven by the Z machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1983

Experiments using atmospheric forcing from a FGGE analysis to drive an upper ocean model

Results are presented of a series of numerical experiments in which an upper ocean model is driven by surface heat fluxes and stress fields derived from the FGGE SOP-1 GLAS analysis/forecast system (Halem, et al., 1982). The model results show that most changes in the mixed layer height and horizontal velocity occurs in the first days. On the other hand, changes in the temperature field take a longer time to develop. In the most realistic case (real initial conditions, instantaneous forcing fields from the atmospheric analysis), the resulting changes in temperature were larger than observed and the correlation between observed and predicted changes was poor. The deficiency in the forecast of SST changes may be due to several factors: lack of sufficient ocean resolution, improper initialization, lack of feedback between the ocean and the atmosphere and the absence of transports by the strong boundary currents, and perhaps unrealistic surface fluxes of heat and momentum. Unless these problems are alleviated it is not reasonable to perform coupled atmospheric ocean forecasts.

Camerlengo, A.↗