Search NASA⌕ Search

SEARCH · Search NASA

Results for “reference model”

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 847 records · Page 47

A multifidelity method for a nonlocal diffusion model

Nonlocal models feature a finite length scale, referred to as the horizon, such that points separated by a distance smaller than the horizon interact with each other. Such models have proven to be useful in a variety of settings. However, due to the reduced sparsity of discretizations, they are also generally computationally more expensive compared to their local differential equation counterparts. In this work, we introduce a multifidelity Monte Carlo method that combines the high-fidelity nonlocal model of interest with surrogate models that use coarser grids and/or smaller horizons and thus have lower fidelities and lower costs. Using the multifidelity method, the overall computational cost of uncertainty quantification is reduced without compromising accuracy. It is shown for a one-dimensional nonlocal diffusion example that speedups of up to two orders of magnitude can be achieved using the multifidelity method to estimate the expectation of an output of interest.

97 MATHEMATICS AND COMPUTING↗

Investigation of assumptions and approximations in the virtual fields method for a viscoplastic material model

The Virtual Fields Method (VFM) is an inverse technique used for parameter estimation and calibration of constitutive models. Many assumptions and approximations—such as plane stress, incompressible plasticity, and spatial and temporal derivative calculations—are required to use VFM with full-field deformation data, for example, from Digital Image Correlation (DIC). Here, this work presents a comprehensive discussion of the effects of these assumptions and approximations on parameters identified by VFM for a viscoplastic material model for 304L stainless steel. We generated synthetic data from a Finite-Element Analysis (FEA) in order to have a reference solution with a known material model and known model parameters, and we investigated four cases in which successively more assumptions and approximations were included in the data. We found that VFM is tolerant to small deviations from the plane stress condition in a small region of the sample, and that the incompressible plasticity assumption can be used to estimate thickness changes with little error. A local polynomial fit to the displacement data was successfully employed to compute the spatial displacement gradients. The choice of temporal derivative approximation (i.e., backwards difference versus central difference) was found to have a significant influence on the computed rate of deformation and on the VFM results for the rate-dependent model used in this work. Finally, the noise introduced into the displacement data from a stereo-DIC simulator was found to have negligible influence on the VFM results. Evaluating the effects of assumptions and approximations using synthetic data is a critical first step for verifying and validating VFM for specific applications. Finally, the results of this work provide the foundation for confidently using VFM for experimental data.

304L stainless steel↗

The Effect of Protoplanetary Disk Cooling Times on the Formation of Gas Giant Planets by Gravitational Instability

Observational evidence exists for the formation of gas giant planets on wide orbits around young stars by disk gravitational instability, but the roles of disk instability and core accretion for forming gas giants on shorter period orbits are less clear. The controversy extends to population synthesis models of exoplanet demographics and to hydrodynamical models of the fragmentation process. The latter refers largely to the handling of radiative transfer in three-dimensional (3D) hydrodynamical models, which controls heating and cooling processes in gravitationally unstable disks, and hence dense clump formation. A suite of models using the β cooling approximation is presented here. The initial disks have masses of 0.091 M {sub ⊙} and extend from 4 to 20 au around a 1 M {sub ⊙} protostar. The initial minimum Toomre Qi values range from 1.3 to 2.7, while β ranges from 1 to 100. We show that the choice of Q {sub i} is equal in importance to the β value assumed: high Q{sub i} disks can be stable for small β , when the initial disk temperature is taken as a lower bound, while low Q{sub i} disks can fragment for high β . These results imply that the evolution of disks toward low Q{sub i} must be taken into account in assessing disk fragmentation possibilities, at least in the inner disk, i.e., inside about 20 au. The models suggest that if low Q{sub i} disks can form, there should be an as yet largely undetected population of gas giants orbiting G dwarfs between about 6 au and 16 au.

79 ASTRONOMY AND ASTROPHYSICS↗

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗

On multichannel film dosimetry with channel-independent perturbations

Purpose: Different multichannel methods for film dosimetry have been proposed in the literature. Two of them are the weighted mean method and the method put forth byMicke et al. [“Multichannel film dosimetry with nonuniformity correction,” Med. Phys. 38, 2523–2534 (2011)] and Mayer et al. [“Enhanced dosimetry procedures and assessment for EBT2 radiochromic film,” Med. Phys. 39, 2147–2155 (2012)]. The purpose of this work was to compare their results and to develop a generalized channel-independent perturbations framework in which both methods enter as special cases. Methods: Four models of channel-independent perturbations were compared: weighted mean, Micke–Mayer method, uniform distribution, and truncated normal distribution. A closed-form formula to calculate film doses and the associated type B uncertainty for all four models was deduced. To evaluate the models, film dose distributions were compared with planned and measured dose distributions. At the same time, several elements of the dosimetry process were compared: film type EBT2 versus EBT3, different waiting-time windows, reflection mode versus transmission mode scanning, and planned versus measured dose distribution for film calibration and for γ-index analysis. The methods and the models described in this study are publicly accessible through IRISEU. Alpha 1.1 ( http://www.iriseu.com ). IRISEU. is a cloud computing web application for calibration and dosimetry of radiochromic films. Results: The truncated normal distribution model provided the best agreement between film and reference doses, both for calibration and γ-index verification, and proved itself superior to both the weighted mean model, which neglects correlations between the channels, and the Micke–Mayer model, whose accuracy depends on the properties of the sensitometric curves. With respect to the selection of dosimetry protocol, no significant differences were found between transmission and reflection mode scanning, between 75 ± 5 min and 20 ± 1 h waiting-time windows or between employing EBT2 or EBT3 films. Significantly better results were obtained when a measured dose distribution was used instead of a planned one as reference for the calibration, and when a planned dose distribution was used instead of a measured one as evaluation for the γ-analysis. Conclusions: The truncated normal distribution model of channel-independent perturbations was found superior to the other three models under comparison and the authors propose its use for multichannel dosimetry.

61 RADIATION PROTECTION AND DOSIMETRY↗

Two-Step Procedure for Liquid-Salt-Cooled-Reactor Analysis

Here, this paper presents the development of a lattice physics–to–core simulator two-step procedure for the rapid analysis of the Advanced High Temperature Reactor (AHTR). Lattice physics, reflector, and control blade models were developed from which cross-section libraries could be generated for a nodal core simulator. Few-group structures for the core simulator were also generated to account for the neutronic characteristics of AHTR. After developing the AHTR two-step procedure, cross-section libraries were generated using the SERPENT continuous-energy Monte Carlo code. These libraries were then used in the core simulator NESTLE to perform full-core calculations, which were in turn benchmarked against reference SERPENT full-core models. Benchmarking results showed reasonable accuracy of the developed two-step procedure but revealed an inherent inadequacy in the one-dimensional radial reflector model and showed a likely need for a greater number of energy groups than were used in this study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A bilayer Double Semion model with symmetry-enriched topological order

We construct a new model of two-dimensional quantum spin systems that combines intrinsic topological orders and a global symmetry called flavour symmetry. It is referred as the bilayer Doubled Semion model (bDS) and is an instance of symmetry-enriched topological order. A honeycomb bilayer lattice is introduced to combine a Double Semion Topological Order with a global spin–flavour symmetry to get the fractionalization of its quasiparticles. The bDS model exhibits non-trivial braiding self-statistics of excitations and its dual model constitutes a Symmetry-Protected Topological Order with novel edge states. This dual model gives rise to a bilayer Non-Trivial Paramagnet that is invariant under the flavour symmetry and the well-known spin flip symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Report on FY 2023 Research and Development on Specially Designed Creep-fatigue Experiments on Alloy 617 in Support of Improving Creep-fatigue Evaluation Approaches

Experimental and numerical studies in support of developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY 2023. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue life curves at long hold times and low strain ranges at elevated temperatures. In this work, the available uniaxial creep-fatigue failure data on Alloy 617 at temperatures of 950°C and 850°C were analyzed to provide a guidance on the development of the extrapolation method and the creep-fatigue failure criteria. A viscoplastic constitutive model for Alloy 617 was adopted to extrapolate the mechanical responses to low strain ranges and long hold times. A set of design curves of Alloy 617 at temperatures of 950°C, 850°C, and 800°C with tensile hold times of 1 hr, 100 hr, and 1,000 hr are developed. Furthermore, creep-fatigue testing on two notch specimen geometries, shallow-notch and sharp V-notch, on Alloy 617 was performed 950°C to understand the multiaxial stress relaxation behavior. The experimental and numerical results on the notch specimens were compared with those on the standard uniaxial smooth bar specimens. The effect of multi-axial stress state combined with elastic follow-up on the stress relaxation behavior was investigated in this report.

36 MATERIALS SCIENCE↗

SU-F-303-17: Real Time Dose Calculation of MRI Guided Co-60 Radiotherapy Treatments On Free Breathing Patients, Using a Motion Model and Fast Monte Carlo Dose Calculation

Purpose: To demonstrate real-time dose calculation of free-breathing MRI guided Co−60 treatments, using a motion model and Monte-Carlo dose calculation to accurately account for the interplay between irregular breathing motion and an IMRT delivery. Methods: ViewRay Co-60 dose distributions were optimized on ITVs contoured from free-breathing CT images of lung cancer patients. Each treatment plan was separated into 0.25s segments, accounting for the MLC positions and beam angles at each time point. A voxel-specific motion model derived from multiple fast-helical free-breathing CTs and deformable registration was calculated for each patient. 3D images for every 0.25s of a simulated treatment were generated in real time, here using a bellows signal as a surrogate to accurately account for breathing irregularities. Monte-Carlo dose calculation was performed every 0.25s of the treatment, with the number of histories in each calculation scaled to give an overall 1% statistical uncertainty. Each dose calculation was deformed back to the reference image using the motion model and accumulated. The static and real-time dose calculations were compared. Results: Image generation was performed in real time at 4 frames per second (GPU). Monte-Carlo dose calculation was performed at approximately 1frame per second (CPU), giving a total calculation time of approximately 30 minutes per treatment. Results show both cold- and hot-spots in and around the ITV, and increased dose to contralateral lung as the tumor moves in and out of the beam during treatment. Conclusion: An accurate motion model combined with a fast Monte-Carlo dose calculation allows almost real-time dose calculation of a free-breathing treatment. When combined with sagittal 2D-cine-mode MRI during treatment to update the motion model in real time, this will allow the true delivered dose of a treatment to be calculated, providing a useful tool for adaptive planning and assessing the effectiveness of gated treatments.

60 APPLIED LIFE SCIENCES↗

Considerations for improving the biosphere assessment methodology and the associated conceptual and mathematical models for application to LLW/ILW radioactive wastes in Spain - 15194

We give consideration to future needs in assessment studies relating to the disposal of very low level radioactive waste (VLLW), low level radioactive waste (LLW) and intermediate level radioactive waste (ILW) in Spain. As these wastes are primarily disposed at El Cabril, particular attention is given to the characteristics of that site. The focus is on post-closure performance assessment studies and specifically on the biosphere aspect of those studies. In particular, it gives consideration to the assessment methodology that should be used, taking into account the recommendations made in the international BIOMASS and BIOCLIM programs, as well as more recent developments. It then goes on to consider the conceptual models of the biosphere that should be developed and comments on how these should be transformed into mathematical models. The issues arising in conceptual and mathematical model development are illustrated by reference to El Cabril, but several of these issues would also be applicable to other engineered facilities excavated from the surface in a variety of geographical and geological contexts. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nuclear Archaeology for Gaseous Diffusion Enrichment Plants

Gaseous diffusion was historically the most widely used technology for military production of highly enriched uranium. Since June 2013, all gaseous diffusion enrichment plants worldwide are permanently shut down. The experience with decommissioning some of these plants has shown that they contain large amounts of uranium particles deposited in the cascade equipment. This article evaluates the potential of using uranium particle deposition to understand and reconstruct the operating histories of gaseous diffusion enrichment plants. First, a squared-off cascade enrichment model is derived to estimate the enrichment capacity of a reference plant. Then, using a cross-flow filtration model, the mass of solid uranium particles deposited over time in the tubular separation membranes of the stage diffusers is calculated. Lastly, potential techniques to characterize these uranium deposits and help reconstruct the operating history of the plant are assessed.

07 ISOTOPE AND RADIATION SOURCES↗

Coda Calibration Tool

CCT is a Java based application for calibrating 10 shear wave coda measurement models to observed data using a much smaller set of reference moment magnitudes (MWs) calculated from other means (waveform modeling, etc.). These calibrated measurement models can then be used in other tools to generate coda moment magnitude measurements, source spectra, estimated stress drop, and other useful measurements for any additional events and any new data collected in the calibrated region.

Addair, Travis↗

The future evolution of energy-water-agriculture interconnectivity across the US

Abstract Energy, water, and agricultural resources across the globe are highly interconnected. This interconnectivity poses science challenges, such as understanding and modeling interconnections, as well as practical challenges, such as efficiently managing interdependent resource systems. Using the US as an example, this study seeks to define and explore how interconnectivity evolves over space and time under a range of influences. Concepts from graph theory and input–output analysis are used to visualize and quantify key intersectoral linkages using two new indices: the ‘Interconnectivity Magnitude Index’ and the ‘Interconnectivity Spread Index’. Using the Global Change Analysis Model (GCAM-USA), we explore the future evolution of these indices under four scenarios that explore a range of forces, including socioeconomic and technological change. Analysis is conducted at both national and state level spatial scales from 2015 to 2100. Results from a Reference scenario show that resource interconnectivity in the US is primarily driven by water use amongst different sectors, while changes in interconnectivity are driven by a decoupling of the water and electricity systems, as power plants become more water-efficient over time. High population and GDP growth results in relatively more decoupling of sectors, as a larger share of water and energy is used outside of interconnected sector feedback loops. Lower socioeconomic growth results in the opposite trend. Transitioning to a low-carbon economy increases interconnectivity because of the expansion of purpose-grown biomass, which strengthens the connections between water and energy. The results highlight that while some regions may experience similar sectoral stress projections, the composition of the intersectoral connectivity leading to that sectoral stress may call for distinctly different multi-sector co-management strategies. The methodology we introduce here can be applied in diverse geographical and sectoral contexts to enable better understanding of where, when, and how coupling or decoupling between sectors could evolve and be better managed.

Khan, Zarrar (ORCID:0000000281478553)↗

Definition of the IEA Wind 22-Megawatt Offshore Reference Wind Turbine

This technical report describes the design of a new 22 Megawatt reference wind turbine (RWT). The turbine model was designed collaboratively by two teams at the Denmark Technical University and at the National Renewable Energy Laboratory within the International Energy Agency (IEA) Wind Technology Commercialization Programme (TCP) Task 55 on Reference Wind Turbines and Farms. RWTs serve an important purpose in the wind energy community, since they provide openly available data for models representative of current wind turbine technology, which can be used by practitioners for a variety of modeling purposes, ranging from aerodynamic, structural, and aeroelastic turbine modeling to wind farm flow modeling, across a range of fidelities. The IEA 22 RWT aims to model machines with projected installation in the 2025-2030 time frame. The turbine has a rotor diameter of 284 meters and a hub height of 170 meters. It is a class 1-B machine with a rotor specific power nearing 350 W m -2 and it is mounted on either a fixed-bottom offshore foundation or a semi-submersible floating platform.

17 WIND ENERGY↗

Disorder effect on the anisotropic resistivity of phosphorene determined by a tight-binding model

In this work we develop a compact multiorbital tight-binding model for phosphorene that accurately describes states near the main band gap. The model parameters are adjusted using as reference the band structure obtained by a density functional theory calculation with the hybrid HSE06 functional. We use the optimized tight-binding model to study the effects of disorder on the anisotropic transport properties of phosphorene. In particular, we evaluate how the longitudinal resistivity depends on the lattice orientation for two typical disorder models: dilute scatterers with high potential fluctuation amplitudes, mimicking screened charges in the substrate, and dense scatterers with lower amplitudes, simulating weakly bounded adsorbates. We show that the intrinsic anisotropy associated to the band structure of this material, although sensitive to the type and intensity of the disorder, is robust.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Retrospective Analysis of Backwater Habitat Availability Using Remote Sensing

Low-velocity channel-margin habitats serve as important nursery habitats for the endangered Colorado pikeminnow (Ptychocheilus lucius) in the middle Green River between Jensen and Ouray, Utah. These habitats, known as backwaters, are associated with emergent sandbars, and are shaped and reformed annually by peak flows. Our recent knowledge about backwater characteristics and dynamics that is summarized in the synthesis report (Grippo et al. 2015) was based on detailed annual survey data from a relatively small sample of backwaters that were collected from 2003 to 2014 and reach-wide evaluations of backwater surface area based on manual interpretation of aerial and satellite imagery. Methods that bridge the gap between the detailed surveys from a small number of backwaters and the reach-wide assessment of their surface area would enable an assessment of the availability of backwater habitats that meet the minimum depth requirements for suitable habitat for Colorado pikeminnow. In 2015 Argonne National Laboratory (Argonne) tested three regression models—linear, multiple, and partial least square (PLS) regression models—for estimating backwater depth using National Agriculture Imagery Program (NAIP) imagery collected in July 2006 that covered the Jensen-Ouray reach of the Green River (Hamada and LaGory 2016). The results suggested that a PLS regression model showed high correlation with the reference depth (R 2 = 0.69) and had the most unbiased and consistent estimate of backwater depth. The results also indicated that the PLS model would provide reasonable estimates of the amount of habitat providing a minimum suitable depth of 30 cm for young-of-the-year Colorado pikeminnow, even though absolute depth estimates may be uncertain for backwater areas deeper than approximately 40 cm. The study also provided insights regarding the amount and the selection of calibration and validation (cal-val) data needed for improving the accuracy of depth prediction.

54 ENVIRONMENTAL SCIENCES↗

Comparing the Performance of Three Land Models in Global C Cycle Simulations: A Detailed Structural Analysis: Structural Analysis of Land Models

Land models are valuable tools to understand the dynamics of global carbon (C) cycle. Various models have been developed and used for predictions of future C dynamics but uncertainties still exist. Diagnosing the models’ behaviors in terms of structures can help to narrow down the uncertainties in prediction of C dynamics. In this study three widely used land surface models, namely CSIRO’s Atmosphere Biosphere Land Exchange (CABLE) with 9 C pools, Community Land Model (version 3.5) combined with Carnegie-Ames-Stanford Approach (CLM-CASA) with 12 C pools and Community Land Model (version 4) (CLM4) with 26 C pools were driven by the observed meteorological forcing. The simulated C storage and residence time were used for analysis. The C storage and residence time were computed globally for all individual soil and plant pools, as well as net primary productivity (NPP) and its allocation to different plant components’ based on these models. Remotely sensed NPP and statistically derived HWSD, and GLC2000 datasets were used as a reference to evaluate the performance of these models. Results showed that CABLE exhibited better agreement with referenced C storage and residence time for plant and soil pools, as compared with CLM-CASA and CLM4. CABLE had longer bulk residence time for soil C pools and stored more C in roots, whereas, CLM-CASA and CLM4 stored more C in woody pools due to differential NPP allocation. Overall, these results indicate that the differences in C storage and residence times in three models are largely due to the differences in their fundamental structures (number of C pools), NPP allocation and C transfer rates. Our results have implications in model development and provide a general framework to explain the bias/uncertainties in simulation of C storage and residence times from the perspectives of model structures.

54 ENVIRONMENTAL SCIENCES↗

Deployment of neural-network-based neutron microscopic cross sections in the Griffin reactor physics application

The capability to utilize neural networks to predict macroscopic and microscopic cross section parametric spaces has been developed for the Griffin reactor physics application. The LibTorch interface enables Griffin's MOOSE-based materials to interact with LibTorch-trained models, allowing for the evaluation of complex macroscopic or microscopic cross section spaces, which are then used to evaluate the neutronic properties of the Griffin finite element model. This study benchmarks traditional ISOXML-formatted tabulation libraries against neural network-based models for 279 nuclides on 20,160 grid points for zero-dimensional and two-dimensional reactor models. Benchmark metrics include the fundamental mode eigenvalue, fission and absorption rates, and various temperature coefficients of reactivity (isothermal, fuel, and moderator). From the perspective of storage space, the complete set of LibTorch models uses 11 MB on disk, compared to the 10 GB for the ISOXML multigroup library that covers the same grid space. For the two-dimensional performance case considered in Griffin, the Torch model uses 97% less RAM than the reference ISOXML dataset while runtime increases by a factor of 3 when using the LibTorch model compared to the ISOXML dataset with multi-linear interpolation. The LibTorch model consistently yields errors within 0.01% for most analyzed quantities except for the temperature coefficients of reactivity where the maximum discrepancies are up to 0.3 $\frac{pcm}{K}$. Due to the neural network attempting to best predict quantities with no regard for a positive or negative bias for any given quantity, predictions may experience random fluctuations, resulting in both positive and negative errors. Future work will entail both depletion and coupled transient analysis to determine the predictive capabilities of Griffin with neural network-based cross sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗