SEARCH · Engineering Papers

Results for “Internal Model Principle”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2022

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗

DOE OSTI2026

Solvent effects on triplet yields in BODIPY-based photosensitizers

Here, we employ molecular dynamics simulations and quantum rate theories to elucidate the complex condensed-phase dynamics underpinning triplet-state formation in organic photosensitizers. Using models informed by first-principles calculations complete with a molecular representation of solvents of different polarities, we elucidate the interplay of the internal and environmental interactions underlying triplet yield. We find that triplet yields depend sensitively on the dielectric stabilization of the charge transfer intermediate that facilitates a transition into the triplet manifold. Our results illustrate the importance of molecularly detailed models in understanding the excited-state internal charge-transfer dynamics of photochemically relevant organic molecules.

Coello Escalante, Leonardo [University of Californ↗

DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗

DOE OSTI2025

Optimal Design and Operation of Intensified Absorbers with 3D-Printed Packing for Solvent-Based CO 2 Capture

Many potential solvent-based carbon capture processes suffer from a high heat of absorption of CO 2 that adversely affects the thermodynamic driving force. While interstage coolers are often used for removing a portion of the generated heat by removing the solvent or a portion of the solvent from a stage and cooling and returning it back to the absorber, they can be placed only at discrete locations in the tower. This work investigates intensified absorbers with 3D-printed packing that includes an internal cooler and therefore can be potentially used for maximizing the operational efficiency of the absorbers for CO 2 capture. The intensified absorber is modeled by using a generic, first-principles, equation-oriented absorber column model. Since the placement of these intensified packings would cause a loss of area/volume used for mass transfer, optimization of the proposed intensified absorber is performed by optimally selecting the locations at which to place these devices and designing them such that the trade-off due to the addition of the heat removal area and the resulting loss in the mass transfer area is accounted for. Results show that optimally placed and designed intensified packings can lead to a significant increase in the capture efficiency of the process in comparison to a similar column with no internal cooling. It is also observed that by optimally placing and designing intensified packings, the lean solvent flow rate to the absorber can be decreased, and the CO 2 lean loading can be increased while still maintaining the same capture efficiency. These process changes can lead to a substantial reduction in the cost of capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Experimental Covariance Determination for Critical Integral Experiments

Integral benchmarks for criticality safety and nuclear data validation require expensive uncertainty quantification studies. In general, uncertainty quantification techniques ignore correlations between experiments and shared parts. Experiments, such as the TEX (Thermal/Epithermal eXperiments) campaigns, consist of many shared parts, such as the ‘Jemima’ HEU fuel plates, which create a strong correlation in their uncertainties. While these correlations are known to exist, they are often not estimated due to the complexity of such calculations. This paper describes an intuitive method of determining the covariance for each of the experimental components, providing a correlation for each family of parts across the multiple cases examined within a benchmark. A proof-of-principle study using the TEX-HEU experimental campaign was performed and verified that the correlations can be calculated with information commonly found in the ICSBEP (International Criticality Safety Benchmark Evaluation Project) benchmarks. This study showed that the introduction of model and experimental covariances reduces the χ 2 per degree of freedom from 2.203 to 1.179, indicating that the omission causes overly pessimistic bias quantifications. This technique can be seamlessly integrated to current benchmark evaluations as well as reevaluations of legacy benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2010

Modeling Social Influence via Combined Centralized and Distributed Planning Control

Real world events are driven by a mixture of both centralized and distributed control of individual agents based on their situational context and internal make up. For example, some people have partial allegiances to multiple, contradictory authorities, as well as to their own goals and principles. This can create a cognitive dissonance that can be exploited by an appropriately directed psychological influence operation (PSYOP). An Autonomous Dynamic Planning and Execution (ADP&E) approach is proposed for modeling both the unperturbed context as well as its reaction to various PSYOP interventions. As an illustrative example, the unrest surrounding the Iranian elections in the summer of 2009 is described in terms applicable to an ADP&E modeling approach. Aspects of the ADP&E modeling process are discussed to illustrate its application and advantages for this example.

Vaccaro, James↗

DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗

NASA NTRSDate not supplied

Effects of Rotation on Internal Structure and Dynamics of Main-Sequence Stars

Current state-of-the-art computer simulations allow us to build 3D dynamical and radiative models of stars from physical first principles with a high degree of realism. The radiative 3D dynamical stellar models obtained with the StellarBox code take into account the effects of turbulence, stellar abundances, a realistic equation of state, and radiative energy transport. In this talk, I will discuss the effects of rotation on the turbulent dynamics and surface structure for a 1.47Msun star for rotational periods of 1 and 14 days. The simulations are performed with the computational domain at various latitudes. The models reproduce stellar granulation, the subsurface shear layer, structural changes in convection, and the tachocline, which is the interface between the inner radiative zone and the outer convection zone and plays a crucial role in stellar variability. In particular, the model results reveal the formation of differential rotation and meridional circulation.

SMD↗

NASA NTRS1992

Interfacial shear modeling and flow predictions for internal film condesation flows

Internal flow of pure vapor experiencing film condesation on the walls of a straight duct is studied. The commonly occuring case of turbulent (or laminar) vapor flow in the core and laminar flow of the liquid condensate-with or without waves on the interface-is emphasized. We propose and implement a new first principle methodolgy which model interfacial shear with the help of reliable experimental data on heat transfer rates. Other details of the flow are predicted with the help of this model. These predictions are shown to be in agreement with relevant experimental data. Correlations for film thickness and heat transfer rates are also given.

Narain, A.↗

NASA NTRS1994

Interlaminar stress analysis of dropped-ply laminated plates and shells by a mixed method

A mixed method of approximation based on Reissner's variational principle is developed for the linear analysis of interlaminar stresses in laminated composites, with special interest in laminates that contain terminated internal plies (dropped-ply laminates). Two models are derived, one for problems of generalized plane deformation and the other for the axisymmetric response of shells of revolution. A layerwise approach is taken in which the stress field is assumed with an explicit dependence on the thickness coordinate in each layer. The dependence of the stress field on the thickness coordinate is determined such that the three-dimensional equilibrium equations are satisfied by the approximation. The solution domain is reduced to one dimension by integration through the thickness. Continuity of tractions and displacements between layers is imposed. The governing two-point boundary value problem is composed of a system of both differential and algebraic equations (DAE's) and their associated boundary conditions. Careful evaluation of the system of DAE's was required to arrive at a form that allowed application of a one-step finite difference approximation. A two-stage Gauss implicit Runge-Kutta finite difference scheme was used for the solution because of its relatively high degree of accuracy. Patch tests of the two models revealed problems with solution accuracy for the axisymmetric model of a cylindrical shell loaded by internal pressure. Parametric studies of dropped-ply laminate characteristics and their influence on the interlaminar stresses were performed using the generalized plane deformation model. Eccentricity of the middle surface of the laminate through the ply drop-off was found to have a minimal effect on the interlaminar stresses under longitudinal compression, transverse tension, and in-plane shear. A second study found the stiffness change across the ply termination to have a much greater influence on the interlaminar stresses.

Harrison, Peter N.↗

DOE OSTI2022

Vibrational Entropy Contribution to Mixing Free Energy of Ni-Rich LiNi 1- y Co y O 2

First principles density-functional-theory calculations are performed of the mixing free energy of the layered oxide LiNi 1-y Co y O 2 , a model Li-ion-battery cathode material, at 1000 K. Analysis of internal energy, LN exchange, configurational entropy and vibrational entropy are combined to obtain the net mixing free energy. The vibrational entropy requires a more indirect treatment than the other terms; analysis based on heat capacity is presented here, to show that the vibrational entropy is of comparable magnitude to the configurational entropy, and it tends to promote solid-solution behavior in LiNi 1-y Co y O 2 .

25 ENERGY STORAGE↗

DOE OSTI2023

3D CAD Modeling and Mechanical Integration at the Mu2e Experiment at Fermilab

This abstract encapsulates a transformative summer internship at Fermilab, focusing on the Mu2e Experiment within the Mechanical Integration and 3D CAD Modeling group. The internship was characterized by a two-phase approach, with the initial month dedicated to the design of an aluminum platform tailored for DS Trench Maintenance. This phase encompassed a comprehensive exploration of structural integrity, material selection, and precision engineering principles, demonstrating an adeptness in CAD modeling and a keen eye for detail. The subsequent month saw a seamless transition into a dynamic phase of the project, wherein the intern played a pivotal role in effecting numerous adjustments across diverse sections of the 3D experiment structure. This phase demanded a nuanced understanding of the intricacies of the Mu2e Experiment and a high level of proficiency in CAD modeling.

3D CAD Modeling↗

DOE OSTI2023

Scientific machine learning for modeling and simulating complex fluids

The formulation of rheological constitutive equations—models that relate internal stresses and deformations in complex fluids—is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these “digital fluid twins” to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation—a task that is not achievable using any previously developed data-driven rheological equation of state.

Science & Technology - Other Topics↗

DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗

DOE OSTI2025

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A computational approach to complex junctions and interfaces

Molecular- and nano- junctions, surfaces, interfaces, and interfacial processes have been the subjects of this DOE supported program for about 20 years. We have investigated electronic and magnetic structure, spin-dependent charge transport, and emerging phenomena at surfaces and interfaces in a number of physical systems ranging from nano-clusters, to nanowires, to two-dimensional (2D) crystals and their interactions, to bulk matter. Eighty-four refereed papers were published [1-84], with a majority effort from the group led by the PI, and over a hundred presentations were made at international conferences. We use first-principles methods in the framework of density functional theory in conjunction with nano-equilibrium Green function techniques, semi-empirical Boltzmann transport theory, tight-binding models, and the Hubbard model and many-body wavefunction approaches. In pursuing our scientific goals, we have also developed a number of computational algorithms and computational packages.

2D materials↗

NASA NTRS1989

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗

NASA NTRS1989

Mixed finite element models for free vibrations of thin-walled beams

Noor, Ahmed K.↗