SEARCH · Search NASA

Results for “TENSOR 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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Aerodynamic Sensitivity of a Novel Data-Driven Airfoil Shape Representation Framework

We explore the aerodynamic implications of a novel data-driven separable shape tensor framework used to represent discrete airfoil shapes. In this study, we construct a data-driven parameter space defined by separable shape tensors and informed by tens of thousands of distinct airfoils. We use this design space to generate new airfoil designs to study parametric sensitivities with respect to various aerodynamic responses. We use a HAM2D RANS solver to approximate the lift, drag, and moment coefficients for the generated airfoils at two different angles-of-attack. We analyze the robustness and sensitivities of using the separable shape tensor design space by examining the coverage of the aerodynamic response space, uncovering low-dimensional polynomial ridge approximations, and computing various sensitivity metrics. The results show that the data-driven design space produce significant variation in target aerodynamic quantities and facilitate highly accurate approximations (R^2 > 0.96) of one- and two-dimensional structures in each aerodynamic response. This further reduces the effective dimension to enable simplified design and optimization tasks.

aerodynamics↗

DOE OSTI2021

Few-nucleon matrix elements in pionless effective field theory in a finite volume

Here, pionless effective field theory in a finite volume (FVEFT $\notπ$ ) is investigated as a framework for the analysis of multinucleon spectra and matrix elements calculated in lattice QCD (LQCD). By combining FVEFT π with the stochastic variational method, the spectra of nuclei with atomic number A ∈ {2, 3} are matched to existing finite-volume LQCD calculations at heavier-than-physical quark masses corresponding to a pion mass m $\notπ$ = 806 MeV , thereby enabling infinite-volume binding energies to be determined using infinite-volume variational calculations. Based on the variational wave functions that are constructed in this approach, the finite-volume matrix elements of various local operators are computed in FVEFT $\notπ$ and matched to LQCD calculations of the corresponding QCD operators in the same volume, thereby determining the relevant one- and two-body effective field theory counterterms and enabling an extrapolation of the LQCD matrix elements to infinite volume. As examples, the scalar, tensor, and axial matrix elements are considered, in addition to the magnetic moments and the isovector longitudinal momentum fraction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr, down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2023

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B -mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr , down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

U. Fuskeland↗

DOE OSTI2022

Primary donor triplet states of Photosystem I and II studied by Q-band pulse ENDOR spectroscopy

The photoexcited triplet state of the “primary donors” in the two photosystems of oxygenic photosynthesis has been investigated by means of electron-nuclear double resonance (ENDOR) at Q-band (34 GHz). The data obtained represent the first set of 1 H hyperfine coupling tensors of the 3 P700 triplet state in PSI and expand the existing data set for 3 P680. We achieved an extensive assignment of the observed electron-nuclear hyperfine coupling constants (hfcs) corresponding to the methine α-protons and the methyl group β-protons of the chlorophyll (Chl) macrocycle. The data clearly confirm that in both photosystems the primary donor triplet is located on one specific monomeric Chl at cryogenic temperature. In comparison to previous transient ENDOR and pulse ENDOR experiments at standard X-band (9–10 GHz), the pulse Q-band ENDOR spectra demonstrate both improved signal-to-noise ratio and increased resolution. The observed ENDOR spectra for 3 P700 and 3 P680 differ in terms of the intensity loss of lines from specific methyl group protons, which is explained by hindered methyl group rotation produced by binding site effects. Contact analysis of the methyl groups in the PSI crystal structure in combination with the ENDOR analysis of 3 P700 suggests that the triplet is located on the Chl a' (P A ) in PSI. The results also provide additional evidence for the localization of 3 P680 on the accessory Chl D1 in PSII.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

DOE OSTI2020

Optical Characterization of the Keck Array and BICEP3 CMB Polarimeters from 2016 to 2019

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale cosmic microwave background (CMB) polarization from the South Pole in search of a primordial B-mode signature. This B-mode signal arises from primordial gravitational waves interacting with the CMB and has amplitude parametrized by the tensor-to-scalar ratio r. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier transform spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on r induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

Effects of Adiabatic Heating on the High Strain Rate Deformation of Polymer Matrix Composites

Polymer matrix composites (PMCs) are increasingly being used in aerospace structures that are expected to experience complex dynamic loading conditions throughout their lifetime. As such, a detailed understanding of the high strain rate behavior of the constituents, particularly the strain rate, temperature, and pressure dependent polymer matrix, is paramount. In this paper, preliminary efforts in modeling experimentally observed temperature rises due to plastic deformation in PMCs subjected to dynamic loading are presented. To this end, an existing isothermal viscoplastic polymer constitutive formulation is extended to model adiabatic conditions by incorporating temperature dependent elastic properties and modifying the components of the inelastic strain rate tensor to explicitly depend on temperature. It is demonstrated that the modified polymer constitutive model is capable of capturing strain rate and temperature dependent yield as well as thermal softening associated with the conversion of plastic work to heat at high rates of strain. The modified constitutive model is then embedded within a strength of materials based micromechanics framework to investigate the manifestation of matrix thermal softening, due to the conversion of plastic work to heat, on the high strain rate response of a T700Epon 862 (T700E862) unidirectional composite. Adiabatic model predictions for high strain rate composite longitudinal tensile, transverse tensile, and in-plane shear loading are presented. Results show a substantial deviation from isothermal conditions; significant thermal softening is observed for matrix dominated deformation modes (transverse tension and in-plane shear), highlighting the importance of accounting for the conversion of plastic work to heat in the polymer matrix in the high strain rate analysis of PMC structures.

polymer matrix composites↗

DOE OSTI2025

Proton internal pressure from deeply virtual Compton scattering on collider kinematics

Abstract The unique experimental connection to the QCD energy–momentum tensor offered by generalised parton distributions has been strongly highlighted in the past few years with attempts to extract the pressure and shear forces distributions within the nucleon. If, in principle, this can be performed in a model independent way from experimental data, in practice, the current limited precision and kinematic coverage make such an extraction very challenging. Moreover, the limitation to a leading-order description in the strong coupling of the data has provided only an indirect and weakly sensitive access to gluon degrees of freedom, solely through their mixing to quarks via evolution. In this paper we address this issue by providing a next-to-leading order formalism allowing a reanalysis of global fits with genuine gluonic degrees of freedom. In addition, we provide an estimate of the reduction in uncertainty that could stem from the extended kinematic range relevant for the future Electron Ion Collider. Finally, we stress the connection between the analysis of the dispersion relation in terms of generalised parton distributions and the deconvolution problem.

Dutrieux, H. (ORCID:0000000183344885)↗

DOE OSTI2024

Standard model Higgs inflation supplemented by minimal dark matter

Renormalization group analysis with the present measurements of the top quark mass m t = 172.69 ± 0.30 GeV [1] indicates that the Standard Model (SM) Higgs potential becomes unstable at energy scales ∼ 10 10 GeV . This may be interpreted as hinting at new particles at high energy. The minimal extension of the SM that can avoid this instability while leaving the SM Higgs as the sole scalar particle of the theory is obtained by adding suitable fermions to the SM. These fermions are good dark matter candidates and the model is known as the minimal dark matter model. We revisit the inflationary scenario based on the minimal dark matter model, taking into account updated parameter constraints and recent understanding of reheating dynamics. We explore the model with different values of the right-handed neutrino mass and find that the cosmological prediction is insensitive to such details. We obtained a spectral index of the cosmic microwave background n s = 0.9672 and a tensor-to-scalar ratio r = 0.0031 as a robust prediction of this scenario. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

NASA NTRS1991

Monte Carlo simulation of errors in the anisotropy of magnetic susceptibility - A second-rank symmetric tensor

Monte Carlo perturbations of synthetic tensors to evaluate the Hext/Jelinek elliptical confidence regions for anisotropy of magnetic susceptibility (AMS) eigenvectors are used. When the perturbations are 33 percent of the minimum anisotropy, both the shapes and probability densities of the resulting eigenvector distributions agree with the elliptical distributions predicted by the Hext/Jelinek equations. When the perturbation size is increased to 100 percent of the minimum eigenvalue difference, the major axis of the 95 percent confidence ellipse underestimates the observed eigenvector dispersion by about 10 deg. The observed distributions of the principal susceptibilities (eigenvalues) are close to being normal, with standard errors that agree well with the calculated Hext/Jelinek errors. The Hext/Jelinek ellipses are also able to describe the AMS dispersions due to instrumental noise and provide reasonable limits for the AMS dispersions observed in two Hawaiian basaltic dikes. It is concluded that the Hext/Jelinek method provides a satisfactory description of the errors in AMS data and should be a standard part of any AMS data analysis.

Lienert, Barry R.↗

DOE OSTI2023

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Planck 2018 results

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = -0.9 ± 5.1; f NL equil = -26 ± 47; and f NL ortho = -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 σ . Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is g NL local = (-5.8 ± 6.5) × 10 4 (68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the ΛCDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Planck constraints on the tensor-to-scalar ratio

In this work, we present constraints on the tensor-to-scalar ratio r using Planck data. We use the latest release of Planck maps (PR4), processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarization for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck's noise level, the TT spectrum gives constraints on r that are cosmic-variance limited (with $\sigma$ r = 0.093), but we show that the marginalized posterior peaks towards negative values of r at about the 1.2$\sigma$ level. We derived Planck constraints using the BB power spectrum at both large angular scales (the 'reionization bump') and intermediate angular scales (the 'recombination bump') from $\ell$=2 to 150, and find a stronger constraint than that from TT, with $\sigma$ r = 0.069. The Planck BB spectrum shows no systematic bias, and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the constraint r < 0.158 at 95% confidence using more than 50% of the sky. This upper limit tightens to r < 0.069 when Planck EE, BB, and EB power spectra are combined consistently, and it tightens further to r < 0.056 when the Planck TT power spectrum is included in the combination. Finally, combining Planck with BICEP2/Keck 2015 data yields an upper limit of r < 0.044.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

The Physics of Turbulence in the Boundary Layer

The geometry of the velocity field in a numerically simulated incompressible turbulent boundary layer over a flat plate at Re theta=670 has been studied using the invariants of the velocity gradient tensor. These invariants are computed at every grid point in the flow and used to form the discriminant. Of primary interest are those regions in the flow where the discriminant is positive; regions where, according to the characteristic equation, the eigenvalues of the velocity gradient tensor are complex. An observer moving with a frame of reference which is attached to a fluid particle lying within such a region would see a local flow pattern of the type stable-focus-stretching or unstable-focus-compressing. When the flow is visualized this way, continuous, connected, large-scale structures are revealed that extend from the point just below the buffer layer out to the beginning of the wake region. These structures are aligned with the mean shear close to the wall and arch in the cross-stream direction away from the wall. In some cases the structures observed are very similar to to the hairpin eddy vision of boundary layer structure proposed by Theodorsen. That the structure of the flow is revealed more effectively by the discriminant rather than by the vorticity is important and adds support to recent observations of the discriminant in a channel flow simulation. Of particular importance is the fact that the procedure does not require the use of an arbitrary threshold in the discriminant. Further analysis using computer flow visualization shows a high degree of spatial correlation between regions of positive discriminant, extreme negative pressure fluctuations and large instantaneous values of Reynolds shear stress.

Kline, Stephen↗

NASA NTRSDate not supplied

Assessing the Relationships Between Sensorimotor Biomarkers and Post-Landing Functional Task Performance

Spaceflight drives adaptive changes in healthy individuals appropriate for sensorimotor function in a microgravity environment. These changes are maladaptive for return to earth's gravity. The inter-individual variability of sensorimotor decrements is striking, although poorly understood. The goal of this study was to identify a set of behavioral, neuroimaging and genetic measures that can be used to predict early post-flight performance on a set of sensorimotor tasks. Astronauts were recruited who previously participated in sensorimotor field tests and/or posturography soon after long-duration spaceflight. Behavioral tests included assessments of sensory dependency and adaptability. Visual dependency involved treadmill walking while viewing a moving virtual visual scene. Vestibular perceptual thresholds were measured while seated during lateral translations. Proprioception dependency was measured during one-legged stance on a horizontal air-bearing surface. Ground assessment of adaptability was performed (1) during treadmill walking with a virtual linear hallway and a moving walking surface, and (2) during multiple trials of navigating an obstacle course while wearing reversing prisms (adaptive Functional Mobility Test, aFMT). The neuroimaging tests characterized individual differences in regional brain volumes (using Structural MRI) and white matter microstructure (using Diffusion Tensor Imaging) to serve as potential predictors of adaptive capacity. The genetic tests utilized saliva samples to examine variations in four genes chosen because of their ability to differentiate sensorimotor adaptation ability in a normative population, including Catechol-O-methyltransferase (COMT), Dopamine Receptor D2 (DRD2), Brain-derived neurotrophic factor (BDNF) and the α2-adrenergic receptor. Thirty ISS crewmembers were recruited, including 12 first-time fliers, 6F, and mission durations lasting 185.5 ± 45.5 days, mean ± std. We utilized a combination of three post-flight functional task outcomes: tandem walk, recovery from fall and dynamic posturography. There was considerable variability among all post-flight performance outcomes. Based on an ordinal scale survey, 72% indicated their ability to perform functional tasks were more impacted postflight relative to inflight with 50% indicating they needed to restrict movements for a longer period postflight relative to inflight. While there is a strong association within tests obtained at different R+0 timepoints, by R+24 hr performance on one post-flight test does not necessarily correlate with performance on other post-flight tests. Preliminary statistical analysis indicates combining biomarkers will increase predictive power and this will be explored with future analyses. Our preliminary findings underscore the importance of a comprehensive post-flight test battery including different types of tasks with varying sensory feedback. Understanding the relationships between these sensorimotor biomarkers and post-flight functional task performance improve both our understanding of the individual variability and our strategy to optimize sensorimotor countermeasures.

biomarkers↗

DOE OSTI2023

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗