Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A microstructure-based modeling approach to predict the mechanical properties of Zr alloy with hydride precipitates

In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.

Kulkarni, Shank S.↗
DOE OSTI2022

Leaching Model of Radionuclides in Metal-Organic Framework Particles

Metal-organic frameworks (MOFs) have been used to sequester radionuclides and seal them inside of porous scaffolds using postsynthetic modification procedures. Experiments show that certain Zr-MOF with different capping linkers significantly affects the radionuclide release kinetics. In this work, we developed a leaching model of radionuclides in Zr-MOF particles. The model assumes that uranyl species occupy two energetically favored sites: the metal node and the MOF pores. For a given overall concentration of uranyl species, the partitions of uranyl species at the metal nodes and within the pores are determined by their chemical potentials. The model also considers the effect of particle surface and concentration on chemical potentials and diffusivity. The effects of spatial and structural dependent chemical potentials and diffusivity as well as particle sizes on leaching kinetics are investigated with the model. Predicted and measured uranyl leaching kinetics in Zr-MOF particles under batch experiments are compared. The results demonstrate the model’s capability for exploring the mechanisms of leaching and provide guidance for material design.

Metal-Organic Framework, Diffusion, uranyl, Leachi↗
DOE OSTI2015

Simulations of ecosystem hydrological processes using a unified multi-scale model

This paper presents a unified multi-scale model (UMSM) that we developed to simulate hydrological processes in an ecosystem containing both surface water and groundwater. The UMSM approach modifies the Navier–Stokes equation by adding a Darcy force term to formulate a single set of equations to describe fluid momentum and uses a generalized equation to describe fluid mass balance. The advantage of the approach is that the single set of the equations can describe hydrological processes in both surface water and groundwater where different models are traditionally required to simulate fluid flow. This feature of the UMSM significantly facilitates modelling of hydrological processes in ecosystems, especially at locations where soil/sediment may be frequently inundated and drained in response to precipitation, regional hydrological and climate changes. In this paper, the UMSM was benchmarked using WASH123D, a model commonly used for simulating coupled surface water and groundwater flow. Disney Wilderness Preserve (DWP) site at the Kissimmee, Florida, where active field monitoring and measurements are ongoing to understand hydrological and biogeochemical processes, was then used as an example to illustrate the UMSM modelling approach. The simulations results demonstrated that the DWP site is subject to the frequent changes in soil saturation, the geometry and volume of surface water bodies, and groundwater and surface water exchange. All the hydrological phenomena in surface water and groundwater components including inundation and draining, river bank flow, groundwater table change, soil saturation, hydrological interactions between groundwater and surface water, and the migration of surface water and groundwater interfaces can be simultaneously simulated using the UMSM. Overall, the UMSM offers a cross-scale approach that is particularly suitable to simulate coupled surface and ground water flow in ecosystems with strong surface water and groundwater interactions.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2017

Anisotropic single-particle dissipative particle dynamics model

We have developed a new single-particle dissipative particle dynamics (DPD) model for anisotropic particles with different shapes, e.g., prolate or oblate spheroids. In particular, the conservative and dissipative interactions between anisotropic single DPD particles are formulated using a linear mapping from the isotropic model of spherical particles. The proper mapping operator is constructed between each interacting pair of particles at every time step. Correspondingly, the random forces are properly formulated to satisfy the fluctuation-dissipation theorem (FDT). Notably, the model exactly conserves both linear and angular momentum. We demonstrate the proposed model's accuracy and efficiency by applying it for modeling colloidal ellipsoids. Specifically, we show it efficiently captures the static properties of suspensions of colloidal ellipsoids. The isotropic-nematic transition in an ellipsoidal suspension is reproduced by increasing its volume fraction or the aspect ratio of ellipsoid particles. Moreover, the hydrodynamics and diffusion of a single colloidal ellipsoid (prolate or oblate with moderate aspect ratios) are accurately captured. The calculated drag force on the ellipsoid and its diffusion coefficients (both translational and rotational) agree quantitatively with the theoretical predictions in the Stokes limit.

Computer Science↗
DOE OSTI2017

Eigendecomposition model of resistance temperature detector with applications to S-CO 2 cycle sensing

Super-critical carbon dioxide (S-CO2) is a promising thermodynamic cycle for advanced nuclear reactors and solar energy conversion applications. Dynamic control of the proposed recompression S-CO2 cycle is accomplished with input from resistance temperature detector (RTD) measurements of the process fluid. One of the challenges in practical implementation of S-CO2 cycle is high corrosion rate of component and sensor materials. In this paper, we develop a mathematical model of RTD sensing using eigendecomposition model of radial heat transfer in a layered long cylinder. We show that the value of RTD time constant primarily depends on the rate of heat transfer from the fluid to the outer wall of RTD. We also show that for typical material properties, RTD time constant can be calculated as the sum of reciprocal eigen-values of the heat transfer matrix. Using the computational model and a set of RTD and CO2 fluid thermophysical parameter values, we calculate the value of time constant of thermowell-mounted RTD sensor at the hot side of the precooler in the S-CO2 cycle. The eigendecomposition model of RTD will be used in future studies to model sensor degradation and its impact on control of S-CO2. (C) 2016 Elsevier B.V. All rights reserved.

finite element↗
DOE OSTI2013

Modeling mixtures of thyroid gland function disruptors in a vertebrate alternative model, the zebrafish eleutheroembryo

Maternal thyroxine (T4) plays an essential role in fetal brain development, and even mild and transitory deficits in free-T4 in pregnant women can produce irreversible neurological effects in their offspring. Women of childbearing age are daily exposed to mixtures of chemicals disrupting the thyroid gland function (TGFDs) through the diet, drinking water, air and pharmaceuticals, which has raised the highest concern for the potential additive or synergic effects on the development of mild hypothyroxinemia during early pregnancy. Recently we demonstrated that zebrafish eleutheroembryos provide a suitable alternative model for screening chemicals impairing the thyroid hormone synthesis. The present study used the intrafollicular T4-content (IT4C) of zebrafish eleutheroembryos as integrative endpoint for testing the hypotheses that the effect of mixtures of TGFDs with a similar mode of action [inhibition of thyroid peroxidase (TPO)] was well predicted by a concentration addition concept (CA) model, whereas the response addition concept (RA) model predicted better the effect of dissimilarly acting binary mixtures of TGFDs [TPO-inhibitors and sodium-iodide symporter (NIS)-inhibitors]. However, CA model provided better prediction of joint effects than RA in five out of the six tested mixtures. The exception being the mixture MMI (TPO-inhibitor)-KClO{sub 4} (NIS-inhibitor) dosed at a fixed ratio of EC{sub 10} that provided similar CA and RA predictions and hence it was difficult to get any conclusive result. There results support the phenomenological similarity criterion stating that the concept of concentration addition could be extended to mixture constituents having common apical endpoints or common adverse outcomes. - Highlights: • Potential synergic or additive effect of mixtures of chemicals on thyroid function. • Zebrafish as alternative model for testing the effect of mixtures of goitrogens. • Concentration addition seems to predict better the effect of mixtures of goitrogens.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to stop unphysical rearrangements during geometry relaxation. A simple fixed-atom or “coordinate-lock” approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

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

Coarse-grained modeling of polyethylene melts: Effect on dynamics

The distinctive viscoelastic behavior of polymers results from a coupled interplay of motion on multiple length and time scales. Capturing the broad time and length scales of polymer motion remains a challenge. Using polyethylene (PE) as a model macromolecule, we construct coarse-grained (CG) models of PE with three to six methyl groups per CG bead and probe two critical aspects of the technique: pressure corrections required after iterative Boltzmann inversion (IBI) to generate CG potentials that match the pressure of reference fully atomistic melt simulations and the transferability of CG potentials across temperatures. While IBI produces nonbonded pair potentials that give excellent agreement between the atomistic and CG pair correlation functions, the resulting pressure for the CG models is large compared with the pressure of the atomistic system. We find that correcting the potential to match the reference pressure leads to nonbonded interactions with much deeper minima and slightly smaller effective bead diameter. However, simulations with potentials generated by IBI and pressure-corrected IBI result in similar mean-square displacements (MSDs) and stress autocorrelation functions G( t ) for PE melts. While the time rescaling factor required to match CG and atomistic models is the same for pressure- and non-pressure-corrected CG models, it strongly depends on temperature. Furthermore, transferability was investigated by comparing the MSDs and stress autocorrelation functions for potentials developed at different temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Improved Models for Metallic Nanoparticle Cores from Atomic Pair Distribution Function (PDF) Analysis

X-ray atomic pair distribution functions (PDFs) were collected from a range of canonical metallic nanomaterials, both elemental and alloyed, prepared using different synthesis methods and exhibiting drastically different morphological properties. Widely applied shape-tuned attenuated crystal (AC) fcc models proved inadequate, yielding structured, coherent, and correlated fit residuals. However, equally simple discrete cluster models could account for the largest amplitude features in these difference signals. A hypothesis testing based approach to nanoparticle structure modeling systematically ruled out effects from crystallite size, composition, shape, and surface faceting as primary factors contributing to the AC misfit. On the other hand, decahedrally twinned cluster cores were found to be the origin of the AC structure misfits for a majority of the nanomaterials reported here. It is further motivated that the PDF can readily differentiate between the arrangement of domains in these multiply twinned motifs. Most of the nanomaterials surveyed also fall within the sub-5 nm size regime where traditional electron microscopy cannot easily detect and quantify domain structures, with sampling representative of the average nanocrystal synthesized. Here, the results demonstrate that PDF analysis is a powerful method for understanding internal atomic interfaces in small noble metallic nanomaterials. Such core cluster models, easily built algorithmically, should serve as starting structures for more advanced models able to capture atomic positional disorder, ligand induced or otherwise, near nanocrystal surfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗