Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Stochastic simulation of vertically nonhomogeneous gusts

The small-scale horizontal gust structure of detailed wind profiles along the vertical in the first 20 km of the atmosphere is a vertically nonhomogeneous process. A linear stochastic model was developed based on the process covariance function. This model was formulated through the use of a scaling hypothesis which transforms the nonhomogeneous gust process into a nondimensional gust process which is homogeneous in a nondimensional gust height coordinate. The velocity scaling parameter for the gust process is the gust standard deviation, and the length scale used to nondimensionalize the altitude is the vertical space lag associated with the first zero of the gust covariance function. State space theory was used to derive a digital filter from the model, which can be readily used to simulate gusts for space vehicle design applications.

Fichtl, G. H.↗
DOE OSTI2022

Atoms and interatomic bonding synergism inherent in molecular electronic wave functions

The chemical model of matter consists of atoms held together by bonds. The success of this model implies that the physical interactions of the electrons and nuclei in molecules combine into compound interactions that create the bonding. In the quantum mechanical description, the modified atoms in molecules and the bonding synergism are contained in the molecular electronic wave function. So far, only part of this information has been recovered from the wave function. Notably, the atoms have remained unidentified in the wave function. One reason is that conventional energy decomposition analyses formulate separate model wave functions, independent of the actual wave function, to represent “prepared atoms” and preconceived interactions and, then, intuitively catenate the parts. In the present work, the embedded modified atoms and the inherent physical synergisms between them are determined by a unified derivation entirely from the actual molecular valence space wave function. By means of a series of intrinsic orbital and configurational transformations of the wave function, the energy of formation of a molecule is additively resolved in terms of intra-atomic energy changes, interference energies, and quasi-classical, non-classical, and charge-transfer Coulombic interactions. Here, the analysis furnishes an algorithm for the quantitative resolution of the energy of formation, which enables analyses elucidating reaction energies.

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

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d -Block

First-principles electronic structure simulations are an invaluable tool for understanding chemical bonding and reactions. While machine-learning models such as interatomic potentials significantly accelerate the exploration of potential energy surfaces, electronic structure information is generally lost. Particularly in the field of heterogeneous catalysis, simulated electron band structures provide fundamental insights into catalytic reactivity. This ab initio knowledge is preserved in semiempirical methods such as density functional tight binding (DFTB), which extend the accessible computational length and time scales beyond first-principles approaches. In this paper here we present Shell-Optimized Atomic Confinement (SOAC) DFTB electronic-part-only parametrizations for bulk and surface band structures of all d-block transition metals that enable efficient predictions of electronic descriptors for large structures or high-throughput studies on complex systems outside the computational reach of density functional theory.

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

Informing Robust Functional Relationship Benchmarks: An Evaluation of the Temperature Sensitivity of Ecosystem Respiration Across the Arctic-Boreal Region

During land model development, simulated carbon dynamics are often benchmarked against observational data sets to evaluate model performance. Functional relationship benchmarks are the relationship between a driving variable (e.g., temperature) and a response variable (e.g., ecosystem respiration) and are a promising tool for assessing model performance by evaluating modeled sensitivities to changing environmental conditions. However, observed functional relationships can be influenced by choices made during data collection and throughout the benchmarking process, impacting the inferred skill of land models. To avoid misrepresenting a model's true performance, it is necessary to systematically evaluate best practices when constructing functional relationship benchmarks. We developed a set of guidelines for constructing functional relationship benchmarks, considering the choice of data set, number of daily observations, temporal extent, and temporal resolution across Alaska and Canada over a 20-year period from 2001 to 2020. The temperature sensitivity of ecosystem respiration from observations, evaluated through an apparent Q 10 , is highly variable both spatially and as a result of the data processing approach applied in the benchmark formation. When benchmarking 13 models from the Warming Permafrost Model Intercomparison Project (WrPMIP), the range in inferred model skill is substantially impacted by the choices applied in constructing functional relationship benchmarks. The inferred performance of a given model is most sensitive to the number of daily observations and temporal extent, followed by choice of benchmark data set and temporal averaging. Results from this analysis can guide the development of consistent and robust functional relationships for future model evaluation studies.

Poe, Jeralyn [Northern Arizona University, Flagsta↗
DOE OSTI2024

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

The E3SM version 1 single-column model

Abstract. The single-column model (SCM) functionality of the Energy Exascale Earth System Model version 1 (E3SMv1) is described in this paper. The E3SM SCM was adopted from the SCM used in the Community Atmosphere Model (CAM) but has evolved significantly since then. We describe changes made to the aerosol specification in the SCM, idealizations, and developments made so that the SCM uses the same dynamical core as the full general circulation model (GCM) component. Based on these changes, we describe and demonstrate the seamless capability to “replay” a GCM column using the SCM. We give an overview of the E3SM case library and briefly describe which cases may serve as useful proxies for replicating and investigate some long-standing biases in the full GCM runs while demonstrating that the E3SM SCM is an efficient tool for both model development and evaluation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Meteor velocity distribution and an optimum monitoring mathematical model

At present, there are a great number of radio meteor, ionosphere and rocket observation data for the altitude range of 80 to 100 km which indicate the existence of large scale circulation systems in the mesopause to low thermosphere range which change regularly with season and latitude. But the existing observation network and observation programs are not optimal for revealing the main factors forming the circulation mode at these altitudes. A generalized optimum monitoring mathematical model is offered for consideration. The model input data are distribution density, response function, individual measurement root mean square uncertainty and detection effectiveness function. The model makes it possible to obtain the observation distribution density, the minimal possible dispersion and optimized system effectiveness.

Volkov, N. G.↗
NASA NTRS1988

Extremely luminous far-infrared colliding galaxies - Computation of colliding-galaxy luminosity function

A computational model used to determine the luminosity function for gaseous galaxies in collision is described. It is based on the assumption of high-velocity encounters of two thin, gas-rich disks having identical radii and thickness, able to approach each other from arbitrary directions, with arbitrary orientations and with arbitrary degree of overlap. The distribution of expected luminosities generated in the dissipation of translational energy is then found to closely fit the observed distribution of luminosities of extremely luminous IRAS galaxies. Bursts of star formation resulting from such collisions would not have the correct luminosity function.

Harwit, Martin↗
NASA NTRS1988

A test of synthetic accretion disk spectra using ultraviolet flux distributions of novalike variables

Ultraviolet (UV) fluxes and other data for a sample of nine novalike cataclysmic variables are assembled from the literature. The UV fluxes and colors are compared to spectra of steady state model accretion disks, constructed using either Planck functions or model stellar atmosphere spectra as elementary radiators. Deficiencies are found with both sets of models. The restriction to steady state radial temperature profiles is relaxed for the 'stellar atmosphere' disks, but these generalized models still fail to account simultaneously for the absolute UV flux and the UV color of the observed spectra. The conclusion is that stellar atmosphere models do not reflect the physics of disks. More appropriate models probably can and definitely should be made.

Wade, Richard A.↗
NASA NTRS1988

Electron-impact excitation of atomic oxygen

A close-coupling approximation is employed to theoretically investigate the excitation of the 3s3S0 and 3p3P states of atomic oxygen for electron impact energies ranging from 13.87-100 eV. The target states included in the expansion of the total wave function are modeled by accurate configuration-interaction wave functions. The results are consistent with recent experimental findings for the dipole-allowed 3P-3s3S0 transition.

Tayal, S. S.↗
NASA NTRS2023

Distorted-Toroidal Flux Rope Model

The 3D characterization of magnetic flux ropes observed in the heliosphere has been a challenging task for decades. This is mainly due to the limitations on inferring the 3D global topology and physical properties from the 1D time series from any spacecraft. To advance our understanding of magnetic flux ropes whose configuration departs from the typical stiff geometries, here we present an analytical solution for a 3D flux rope model with an arbitrary cross section and a toroidal global shape. This constitutes the next level of complexity following the elliptic-cylindrical (EC) geometry. The mathematical framework was established by Nieves-Chinchilla et al. with the EC flux rope model, which describes a magnetic topology with an elliptical cross section as a first approach to changes in the cross section. In the distorted-toroidal flux rope model, the cross section is described by a general function. The model is completely described by a nonorthogonal geometry and the Maxwell equations can be consistently solved to obtain the magnetic field and relevant physical quantities. As a proof of concept, this model is generalized in terms of the radial dependence of current density components. The last part of this paper is dedicated to a specific function, F (φ) = δ (1 - λ cosφ), to illustrate possibilities of the model. This model paves the way toward the investigation of complex distortions of magnetic structures in the solar wind. Future investigations will explore these distortions in depth by analyzing specific events; studying implications for physical quantities, such as magnetic fluxes, helicity, or energy; and evaluating the force balance with the ambient solar wind that allows such distortions.

Solar wind↗
NASA NTRS1986

Reflection and emission models for deserts derived from Nimbus-7 ERB scanner measurements

Broadband shortwave and longwave radiance measurements obtained from the Nimbus-7 Earth Radiation Budget scanner were used to develop reflectance and emittance models for the Sahara-Arabian, Gibson, and Saudi Deserts. The models were established by fitting the satellite measurements to analytic functions. For the shortwave, the model function is based on an approximate solution to the radiative transfer equation. The bidirectional-reflectance function was obtained from a single-scattering approximation with a Rayleigh-like phase function. The directional-reflectance model followed from integration of the bidirectional model and is a function of the sum and product of cosine solar and viewing zenith angles, thus satisfying reciprocity between these angles. The emittance model was based on a simple power-law of cosine viewing zenith angle.

Staylor, W. F.↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
NASA NTRS1996

Stress and Damage in Polymer Matrix Composite Materials Due to Material Degradation at High Temperatures

This report describes analytical methods for calculating stresses and damage caused by degradation of the matrix constituent in polymer matrix composite materials. Laminate geometry, material properties, and matrix degradation states are specified as functions of position and time. Matrix shrinkage and property changes are modeled as functions of the degradation states. The model is incorporated into an existing composite mechanics computer code. Stresses, strains, and deformations at the laminate, ply, and micro levels are calculated, and from these calculations it is determined if there is failure of any kind. The rationale for the model (based on published experimental work) is presented, its integration into the laminate analysis code is outlined, and example results are given, with comparisons to existing material and structural data. The mechanisms behind the changes in properties and in surface cracking during long-term aging of polyimide matrix composites are clarified. High-temperature-material test methods are also evaluated.

McManus, Hugh L.↗
DOE OSTI2023

Multi-objective Bayesian alloy design using multi-task Gaussian processes

In design applications, correlations among material properties (such as the tendency for stronger materials to be less ductile) are often neglected. This approach is echoed in multi-objective optimization techniques which treat each performance characteristic as an independent objective, aiming to optimize scalar functions and find optimal Pareto fronts. However, this overlooks the statistical relationships between performance characteristics inherent in a material system. To address this, we propose the use of Bayesian optimization, a highly efficient black-box optimization algorithm known for constructing Gaussian processes (GPs) – uncorrelated surrogates - to model objective functions. Rather than evaluating multiple GPs for each objective function separately, we argue for a shift towards jointly modeling these objective functions, considering their statistical correlations. This integrated approach utilizes naturally occurring relationships among material properties, providing additional information to enhance the performance of the design framework. This requires the replacement of multiple independent GPs with a single multi-task GP, employing a correlation matrix to construct a multi-task kernel function, wherein each task corresponds to a single objective function. Here, we anticipate this refined methodology will better leverage material correlations, improving design optimization results.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Kinetic modeling of microbial growth, enzyme activity, and gene deletions: An integrated model of β-glucosidase function in Cellvibrio japonicus

Understanding the complex growth and metabolic dynamics in microorganisms requires advanced kinetic models containing both metabolic reactions and enzymatic regulation to predict phenotypic behaviors under different conditions and perturbations. Most current kinetic models lack gene expression dynamics and are separately calibrated to distinct media, which consequently makes them unable to account for genetic perturbations or multiple substrates. This challenge limits our ability to gain a comprehensive understanding of microbial processes towards advanced metabolic optimizations that are desired for many biotechnology applications. Here, we present an integrated computational and experimental approach for the development and optimization of mechanistic kinetic models for microbial growth and metabolic and enzymatic dynamics. Our approach integrates growth dynamics, gene expression, protein secretion, and gene-deletion phenotypes. We applied this methodology to build a dynamic model of the growth kinetics in batch culture of the bacterium Cellvibrio japonicus grown using either cellobiose or glucose media. The model parameters were inferred from an experimental data set using an evolutionary computation method. The resulting model was able to explain the growth dynamics of C. japonicus using either cellobiose or glucose media and was also able to accurately predict the metabolite concentrations in the wild-type strain as well as in β-glucosidase gene deletion mutant strains. We validated the model by correctly predicting the non-diauxic growth and metabolite consumptions of the wild-type strain in a mixed medium containing both cellobiose and glucose, made further predictions of mutant strains growth phenotypes when using cellobiose and glucose media, and demonstrated the utility of the model for designing industrially-useful strains. Importantly, the model is able to explain the role of the different β-glucosidases and their behavior under genetic perturbations. Here, this integrated approach can be extended to other metabolic pathways to produce mechanistic models for the comprehensive understanding of enzymatic functions in multiple substrates.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1993

Starburst modeling of M82 - Test case for a biased initial mass function

We compute starburst models for M82, making use of recent theoretical tracks of stellar evolution. Detailed comparisons of our models and those of others demonstrate this technique to be quite reliable, with relatively little change in output parameters as a function of the selection of theoretical tracks or of estimates of the observational characteristics of the stars along these tracks. The models are matched to the observational constraints for M82 summarized by McLeod et al. (1993). The rate of star formation and time of observation were thoroughly optimized to produce the most favorable fit to the observations, but we still found that the recently proposed forms for the solar-neighborhood IMF cannot produce starbursts adequate to fit the observations of this galaxy. We then explored adjustments to the shape of the IMF to improve the fit to M82. We find (1) the shape of the IMF for high-mass stars need not be different from that observed locally; and (2) the most likely modification to the IMF in M82 is that stars with masses below a few solar masses form much less commonly than in the solar neighborhood.

Rieke, G. H.↗