Search NASA⌕ Search

SEARCH · Search NASA

Results for “model complexity”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Using CUBIT to Create Unstructured Mesh Models for MCNP Simulations

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because creating CSG models is a time-consuming and error-prone process as the complexities of geometries increase. An MCNP UM calculation requires UM geometry input files. The UM capability was originally designed to work with UM models created with the Abaqus/CAE software and ASCII input files that it generates. The Abaqus-formatted input files needed for MCNP UM calculations must have the correct Abaqus syntax and meet the additional MCNP requirements. Several software packages can generate a UM model formatted as an Abaqus input file. Cubit, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files created by Cubit cannot be used for MCNP simulations. A Python script has been developed to convert an Abaqus file created by Cubit to an Abaqus file that the MCNP code can process. This report describes the process of using Cubit to create UM models for MCNP calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Developing a robust strength model using physically-informed genetic programming

The strength of materials is influenced by a range of external conditions, such as temperature and deformation rate. Consequently, materials that demonstrate substantial variations in their mechanical behavior due to fluctuations in temperature and strain rate require complex strength models to accurately predict material performance in real-world applications. To predict such complex behavior, a robust and flexible strength model is necessary. In this work, we utilize genetic programming-based symbolic regression (GPSR) to develop data-driven strength models that accurately represent the measured stress–strain responses of tin across a wide range of strain, strain rate and temperature regimes. The GPSR models are constrained by physically-informed conditions, which leads to significant improvement in extrapolation. The best model is integrated into a multi-physics code to perform Taylor impact simulations, validating the model’s accuracy and robustness. In conclusion, the model predictions showed excellent agreement with experimental results, particularly when compared to predictions using traditional strength models.

Genetic programming↗
DOE OSTI2021

Complexity growth in integrable and chaotic models

We use the SYK family of models with N Majorana fermions to study the complexity of time evolution, formulated as the shortest geodesic length on the unitary group manifold between the identity and the time evolution operator, in free, integrable, and chaotic systems. Initially, the shortest geodesic follows the time evolution trajectory, and hence complexity grows linearly in time. We study how this linear growth is eventually truncated by the appearance and accumulation of conjugate points, which signal the presence of shorter geodesics intersecting the time evolution trajectory. By explicitly locating such “shortcuts” through analytical and numerical methods, we demonstrate that: (a) in the free theory, time evolution encounters conjugate points at a polynomial time; consequently complexity growth truncates at O($\sqrt{N}$), and we find an explicit operator which “fast-forwards” the free N-fermion time evolution with this complexity, (b) in a class of interacting integrable theories, the complexity is upper bounded by O(poly(N)), and (c) in chaotic theories, we argue that conjugate points do not occur until exponential times O(e N ), after which it becomes possible to find infinitesimally nearby geodesics which approximate the time evolution operator. Finally, we explore the notion of eigenstate complexity in free, integrable, and chaotic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy

A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.

58 GEOSCIENCES↗
DOE OSTI2021

Spectroscopic factors in dripline nuclei

Single-nucleon knockout reaction studies of the proton-dripline nuclei 9 C and 13 O suggest an appreciable suppression of spectroscopic factors. In this work, we calculate the one-neutron and one-proton spectroscopic factors for the mirror pair 9 C – 9 Li and 13 O using two variants of the continuum shell model: the complex-energy Gamow shell model and the real-energy shell model embedded in the continuum. Furthermore, our results indicate that the continuum effects strongly suppress the spectroscopic factors of well-bound orbits in the dripline systems but have less impact on the spectroscopic factors of weakly bound states.

6 ≤ A ≤ 19↗
DOE OSTI2023

Mobiliti v1.0

Mobiliti is a software platform designed to emulate the dynamics of a regional transportation road network. It is built on open-source software that provides parallel discrete-event simulation. The software is transformative in the area of transportation network simulation because of the geospatial scale and fidelity of the network model and the computational time it takes to model a full day of travel demand. For example, it runs a simulation of the entire San Francisco Bay Area, with a network model of ~1M links and a population that completes ~19M trips in ~5 minutes. This scale of simulation has not been attempted with existing simulation models due to the complexity of the model and the computational time it would take to complete. The intent of the software is to create a digital twin capability for cities to evaluate consequences of infrastructure or policy changes on road network dynamics.

Macfarlane, Jane↗
DOE OSTI2025

Pathways for Nucleation and Growth in Confined Spaces and at Interfaces

Mineral crystallization is central to myriad natural processes from the formation of snowflakes to stalagmites, but the molecular-scale mechanisms are often far more complex than models reflect. Feedbacks between the hydro-, bio-, and geo-spheres drive complex crystallization processes that challenge our ability to observe and quantify them, motivating an expansion of crystallization theories. Here, in this article, we discuss how the driving forces and timescales of nucleation are influenced by factors ranging from simple geometric confinement to distinct interfacial solution structures involving solvent organization, electrical double layers, and surface charging effects. Taken together, these ubiquitous natural phenomena can preserve metastable intermediates, drive precipitation of undersaturated phases, and modulate crystallization in time and space.

biosphere↗
DOE OSTI2020

Generating MCNP Input Files for Unstructured Mesh Geometries

Los Alamos National Laboratory's (LANL) Monte Carlo N-Particle (MCNP) transport code version 6 has the capability for tracking particles on unstructured mesh (UM) geometry models. The MCNP UM feature has been developed for performing calculations of complex geometry models. This capability tracks particles on hybrid geometries where finite element meshes are embedded into constructive solid geometry (CSG) cells. The MCNP UM feature was originally designed to read UM models created by Abaqus/CAE software suite. MCNP versions 6.2.0 and later can process UM models read from Abaqus input les or MCNPUM les converted from Abaqus input les. Sandia National Laboratory's Cubit Toolkit and other finite element analysis software packages may generate UM models and then convert these models into Abaqus input formats.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy

Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.

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

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

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

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Large Eddy Simulations of Turbulent and Buoyant Flows in Urban and Complex Terrain Areas Using the Aeolus Model

Fast and accurate predictions of the flow and transport of materials in urban and complex terrain areas are challenging because of the heterogeneity of buildings and land features of different shapes and sizes connected by canyons and channels, which results in complex patterns of turbulence that can enhance material concentrations in certain regions. To address this challenge, we have developed an efficient three-dimensional computational fluid dynamics (CFD) code called Aeolus that is based on first principles for predicting transport and dispersion of materials in complex terrain and urban areas. The model can be run in a very efficient Reynolds average Navier–Stokes (RANS) mode or a detailed large eddy simulation (LES) mode. The RANS version of Aeolus was previously validated against field data for tracer gas and radiological dispersal releases. As a part of this work, we have validated the Aeolus model in LES mode against two different sets of data: (1) turbulence quantities measured in complex terrain at Askervein Hill; and (2) wind and tracer data from the Joint Urban 2003 field campaign for urban topography. As a third set-up, we have applied Aeolus to simulate cloud rise dynamics for buoyant plumes from high-temperature explosions. For all three cases, Aeolus LES predictions compare well to observations and other models. These results indicate that Aeolus LES can be used to accurately simulate turbulent flow and transport for a wide range of applications and scales.

42 ENGINEERING↗
DOE OSTI2023

Comparison of Zirconium Redistribution in BISON EBR-II Models Using FIPD and IMIS Databases with Experimental Post Irradiation Examination

Metallic fuels have seen increased interest for future sodium fast reactors due to their material properties: high thermal conductivities and advantageous neutronic properties allow for greater fission densities. One drawback to typical metallic fuels is zirconium redistribution, which impacts this advantageous material and its neutronic properties. Unfortunately, the processes behind zirconium migration behavior are understood using first principles, so before these fuels are implemented in future fast reactors, characterization and fuel qualification regimes must be completed. These activities can be supported through the use of robust modeling using the most accurate empirical models currently available to fuel researchers around the world. The tool that allows researchers to model this complex coupled thermo-mechanical behavior and nuclear properties is BISON. Additionally, BISON model parameters need to be compared against PIE measurements. The current work utilizes two fuel pins from EBR-II experiment X441 to optimize various model parameters, including porosity correction factor, thermal conductivity, phase transition temperature, and diffusion coefficient multipliers, before implementing the final model for seven fuel pins with differing characteristics. To properly evaluate the BISON simulations, the results are compared to PIE metallography data for each fuel pin, to ensure the zirconium redistribution is properly reflected in the simulation results. Six out of seven analyzed fuel pins demonstrate good agreement between the metallography images and BISON results, showing alignment of the Zr-rich, Zr-depleted, and moderately Zr-enriched zones at various axial heights along the fuel pins. Further work is needed to refine the model parameters for general pin use.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Particle Receiver Models for Systems Analysis

Particle receivers are gaining importance in the field of Concentrating Solar Power (CSP) due to the high temperature that particles can achieve without degradation. Several researchers are studying the potential of this technology by means of system analyses, which need simple and light models of the system. This study presents two simple models for the particle receiver. The simplest model is a correlation obtained by fitting the results calculated with a more complex receiver model simulated in CFD. The other model is a 1D model, which is benchmarked against the same CFD results. Although both models achieve high coefficient of determination, R2, when compared to CFD results, the 1D model seems to provide more accurate results (especially during sunsets and sunrises). Both models are integrated into a tecno-economic model developed in previous work. The LCOE obtained with the 1D model is between 7% and 10% greater than the one obtained with the correlation.

González-Portillo, Luis F. (ORCID:0000000348643825↗
DOE OSTI2021

A Model for the Signal Initiation Complex Between Arrestin-3 and the Src Family Kinase Fgr

Arrestins regulate a wide range of signaling events, most notably when bound to active G protein-coupled receptors (GPCRs). Among the known effectors recruited by GPCR-bound arrestins are Src family kinases, which regulate cellular growth and proliferation. Here, we focus on arrestin-3 interactions with Fgr kinase, a member of the Src family. Previous reports demonstrated that Fgr exhibits high constitutive activity, but can be further activated by both arrestin-dependent and arrestin-independent pathways. We report that arrestin-3 modulates Fgr activity with a hallmark bell-shaped concentration-dependence, consistent with a role as a signaling scaffold. Here, we further demonstrate using NMR spectroscopy that a polyproline motif within arrestin-3 interacts directly with the SH3 domain of Fgr. To provide a framework for this interaction, we determined the crystal structure of the Fgr SH3 domain at 1.9 Å resolution and developed a model for the GPCR-arrestin-3-Fgr complex that is supported by mutagenesis. This model suggests that Fgr interacts with arrestin-3 at multiple sites and is consistent with the locations of disease-associated Fgr mutations. Collectively, these studies provide a structural framework for arrestin-dependent activation of Fgr.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Standoff Infrared Measurements of Chemical Plume Dynamics in Complex Terrain Using a Combination of Active Swept-ECQCL Laser Spectroscopy with Passive Hyperspectral Imaging

Chemical plume detection and modeling in complex terrain present numerous challenges. We present experimental results from outdoor releases of two chemical tracers (sulfur hexafluoride and Freon-152a) from different locations in mountainous terrain. Chemical plumes were detected using two standoff instruments collocated at a distance of 1.5 km from the plume releases. A passive long-wave infrared hyperspectral imaging system was used to show time- and space-resolved plume transport in regions near the source. An active infrared swept-wavelength external cavity quantum cascade laser system was used in a standoff configuration to measure quantitative chemical column densities with high time resolution and high sensitivity along a single measurement path. Both instruments provided chemical-specific detection of the plumes and provided complementary information over different temporal and spatial scales. The results show highly variable plume propagation dynamics near the release points, strongly dependent on the local topography and winds. Effects of plume stagnation, plume splitting, and plume mixing were all observed and are explained based on local topographic and wind conditions. Measured plume column densities at distances ~100 m from the release point show temporal fluctuations over ~1 s time scales and spatial variations over ~1 m length scales. The results highlight the need for high-speed and spatially resolved measurement techniques to provide validation data at the relevant spatial and temporal scales required for high-fidelity terrain-aware microscale plume propagation models.

47 OTHER INSTRUMENTATION↗