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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Polyconvex neural network models of thermoelasticity

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. Here, in this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Deep neural network for the dielectric response of insulators

In this work, we introduce a deep neural network to model in a symmetry preserving way the environmental dependence of the centers of the electronic charge. The model learns from ab initio density functional theory, wherein the electronic centers are uniquely assigned by the maximally localized Wannier functions. When combined with the deep potential model of the atomic potential energy surface, the scheme predicts the dielectric response of insulators for trajectories inaccessible to direct ab initio simulation. The scheme is nonperturbative and can capture the response of a mutating chemical environment. We demonstrate the approach by calculating the infrared spectra of liquid water at standard conditions, and of ice under extreme pressure, when it transforms from a molecular to an ionic crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Computational Tools for Additive Manufacture of Tailored Microstructure and Properties

Additive manufacturing has the potential to revolutionize industrial hardware and unlock efficiency gains through the fabrication of geometries and architectures not possible by conventional processing. Currently most additive builds use a single set of process parameters (e.g. laser power and scan speed) which results in a part with a homogenous microstructure that provides a singular performance level. To move beyond this state, Raytheon Technologies Research Center worked to create a set of computational tools to track material evolution through each step of the additive process. Computational fluid dynamics and phase field models for microstructure evolution as a function of processing parameters, and crystal plasticity models fully coupling microstructure and mechanical properties for performance predictions were leveraged to establish a connection between additive parameters and the final microstructure. This framework was utilized to tailor spatially-varying mechanical properties in a part by appropriately controlling the microstructure evolution during the additive process. Specifically, a turbine blade was 3D printed from nickel superalloy IN718 using laser powder bed fusion with coarse grains in the airfoil section which experiences the highest temperatures and is creep limited while finer grains were printed in the root of the blade which experience higher stresses but at lower temperatures and is therefore fatigue limited. The benefit of being able to intentionally insert coarse grains in the high temperature region of the blade was showcased with a microstructure sensitive creep model that indicates longer creep life for coarser grains.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2021

mosartwmpy: A Python implementation of the MOSART-WM coupled hydrologic routing and water management model

mosartwmpy is a Python implementation of the Model for Scale Adaptive River Transport with Water Management (MOSART-WM). This new version retains the functionality of the legacy model (written in FORTRAN) while providing new features to enhance user experience and extensibility. MOSART is a large-scale river-routing model used to study riverine dynamics of water, energy, and biogeochemistry cycles across local, regional, and global scales (Li et al., 2013). The WM component introduced by Voisin et al. (2013) represents river regulation through reservoir storage and release operations, diversions from reservoir releases, and allocation to sectoral water demands. Each reservoir release is independently calibrated using long-term mean monthly inflow into the reservoir, long-term mean monthly demand associated with this reservoir, and reservoir goals (flood control, irrigation, recreation, etc.). Generic monthly pre-release rules and storage targets are set up for individual reservoirs; however, those releases are updated annually for inter-annual variability (dry or wet year) and daily for environmental constraints such as flow minimum release and minimum/maximum storage levels. The WM model allows an evaluation of the impact of water management over multiple river basins at once (global, continental scales) and with consistent representation of human operations over the full domain.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Accounting for herbaceous communities in process‐based models will advance our understanding of “grassy” ecosystems

Abstract Grassland and other herbaceous communities cover significant portions of Earth's terrestrial surface and provide many critical services, such as carbon sequestration, wildlife habitat, and food production. Forecasts of global change impacts on these services will require predictive tools, such as process‐based dynamic vegetation models. Yet, model representation of herbaceous communities and ecosystems lags substantially behind that of tree communities and forests. The limited representation of herbaceous communities within models arises from two important knowledge gaps: first, our empirical understanding of the principles governing herbaceous vegetation dynamics is either incomplete or does not provide mechanistic information necessary to drive herbaceous community processes with models; second, current model structure and parameterization of grass and other herbaceous plant functional types limits the ability of models to predict outcomes of competition and growth for herbaceous vegetation. In this review, we provide direction for addressing these gaps by: (1) presenting a brief history of how vegetation dynamics have been developed and incorporated into earth system models, (2) reporting on a model simulation activity to evaluate current model capability to represent herbaceous vegetation dynamics and ecosystem function, and (3) detailing several ecological properties and phenomena that should be a focus for both empiricists and modelers to improve representation of herbaceous vegetation in models. Together, empiricists and modelers can improve representation of herbaceous ecosystem processes within models. In so doing, we will greatly enhance our ability to forecast future states of the earth system, which is of high importance given the rapid rate of environmental change on our planet.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Modeling the Environment-Dependent Kinetics of Oxygen Reduction Reaction – a Continuum Model for Electric Double Layer

Here, for proton-exchange-membrane fuel cells (PEMFCs) to achieve broad commercialization, improved energy-conversion efficiency with minimal Pt-based electrocatalyst is required. Because the sluggish rate of oxygen reduction reaction (ORR) limits the efficiency of PEMFCs, the efficiency improvement requires a better understanding of ORR kinetics and mechanism to design better catalyst. To understand the ORR mechanism, theoretical and experimental analyses have been conducted. While previous studies reasonably explained the catalyst-dependent activity on single crystal catalysts in 0.1 M perchloric acid solution, the explicit effect of electrolyte and related microenvironments is not thoroughly understood. The change in the electrolyte alters the electric-double-layer (EDL) structure and thus the local microenvironment at the electrode/electrolyte interface. Thus, the structure of the EDL should be carefully analyzed to uncover the electrolyte-dependent reaction kinetics. In this talk, we propose a multiscale continuum model to predict the EDL structure and examine the effect of perchloric acid concentration on ORR activity on Pt (111). The model includes Density Potential Functional Theory (DPFT) for electron density and Modified Poisson Boltzmann equation for species’ density and electric potential. Also, the interaction between adsorbents and electric field is taken into account by minimizing the grand potential. After model validation with experimentally measured double-layer capacity data as a function of applied potential and concentration, the effect of the perchloric acid concentration (0.02 M – 0.2 M) on ORR activity is analyzed and discussed. It is shown that the model reproduces the specific activity obtained in the experiments when assuming the oxygen adsorption is limiting the rate, which can be attributed to the large energetic barrier for solvent reorganization. Then, extension of the model to PEMFC ionomer electrolytes will be introduced. Overall, the model framework and findings provide insights into the ORR mechanism and guidance on how to tailor catalyst materials for increased PEMFC performance.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2025

Large-scale tearing-mode hazard function analysis with standard matched equilibrium reconstructions

The association between features from standard tokamak equilibrium reconstructions and the onset of n = 1 tearing modes (TMs) is analyzed at scale. The TM onset rate is directly modeled with a ‘hazard’ function which gives the expected number of onsets (per unit time spent) in a given equilibrium parameter region. In particular the different statistical modeling performance achieved for magnetics-only reconstructions and motional Stark effect (MSE) enhanced reconstructions is studied. It is observed that a better hazard model for the TM onset rate can be built with the MSE-enhanced equilibria compared to the matched magnetics-only situation. This advantage disappears if internal profile details are withheld from the matched analysis. Plausibility of the hazard function is further demonstrated with visualizations of global trends in the operational space, and time-traces from specific tokamak discharges. As a result, TMs typically degrade tokamak plasma performance and may lead to plasma termination, motivating this statistical study.

equilibrium↗

DOE OSTI2024

Large Divergence of Projected High Latitude Vegetation Composition and Productivity Due To Functional Trait Uncertainty

Abstract Vegetation distribution and composition are expected to change in northern high latitudes under rapid warming, which regulates ecosystem functions but remains challenging to predict. Vegetation change arises from the interplay of chronic climate trends such as warming and transient demographic processes of recruitment, growth, competition, and mortality. Most predictive models overlooked the role of demographic dynamics controlled by plant traits. Here, we simulate vegetation dynamics at the Kougarok Hillslope site in Alaska under historical and future climates using the E3SM Land Model coupled to the Functionally Assembled Terrestrial Simulator (ELM‐FATES). To evaluate the roles of plant traits, we parameterize the model with 5,265 trait configurations representing diverse physiological and demographic strategies. Results show current modeled biomass, composition, and productivity are most sensitive to traits controlling photosynthetic capacity, carbon allocation, allometry, and phenology. Among all trait configurations, ∼5% reproduce in situ biomass and plant functional type (PFT) composition measured in 2016, that are indistinguishable from these two observed ecosystem states. Notably, these same trait configurations produce diverging biomass, composition, and productivity under future climate, where the uncertainty attributable to traits is twice the change attributable to climate change. The variation of projected productivity arises from emerging PFT composition under novel climate regimes, primarily explained by traits controlling cold‐induced mortality, recruitment, and allometry. Our findings highlight the importance and uncertainty of demographic dynamics and its interaction with climate change in shaping Arctic vegetation change. Improved model predictions will likely benefit from explicit consideration of vegetation demography and better constraints of critical traits.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Younger facial looks are associate with a lower likelihood of several age-related morbidities in the middle-aged to elderly

Abstract Background Looking older for one’s chronological age is associated with a higher mortality rate. Yet it remains unclear how perceived facial age relates to morbidity and the degree to which facial ageing reflects systemic ageing of the human body. Objectives To investigate the association between ΔPA and age-related morbidities of different organ systems, where ΔPA represents the difference between perceived age (PA) and chronological age. Methods We performed a cross-sectional analysis on data from the Rotterdam Study, a population-based cohort study in the Netherlands. High-resolution facial photographs of 2679 men and women aged 51.5–87.8 years of European descent were used to assess PA. PA was estimated and scored in 5-year categories using these photographs by a panel of men and women who were blinded for chronological age and medical history. A linear mixed model was used to generate the mean PAs. The difference between the mean PA and chronological age was calculated (ΔPA), where a higher (positive) ΔPA means that the person looks younger for their age and a lower (negative) ΔPA that the person looks older. ΔPA was tested as a continuous variable for association with ageing-related morbidities including cardiovascular, pulmonary, ophthalmological, neurocognitive, renal, skeletal and auditory morbidities in separate regression analyses, adjusted for age and sex (model 1) and additionally for body mass index, smoking and sun exposure (model 2). Results We observed 5-year higher ΔPA (i.e. looking younger by 5 years for one’s age) to be associated with less osteoporosis [odds ratio (OR) 0.76, 95% confidence interval (CI) 0.62–0.93], less chronic obstructive pulmonary disease (OR 0.85, 95% CI 0.77–0.95), less age-related hearing loss (model 2; B = −0.76, 95% CI −1.35 to −0.17) and fewer cataracts (OR 0.84, 95% CI 0.73–0.97), but with better global cognitive functioning (g-factor; model 2; B = 0.07, 95% CI 0.04–0.10). Conclusions PA is associated with multiple morbidities and better cognitive function, suggesting that systemic ageing and cognitive ageing are, to an extent, externally visible in the human face.

Mekić, Selma↗

DOE OSTI2023

Modeling immunity in microphysiological systems

There is a need for better predictive models of the human immune system to evaluate safety and efficacy of immunomodulatory drugs and biologics for successful product development and regulatory approvals. Current in vitro models, which are often tested in two-dimensional (2D) tissue culture polystyrene, and preclinical animal models fail to fully recapitulate the function and physiology of the human immune system. Microphysiological systems (MPSs) that can model key microenvironment cues of the human immune system, as well as of specific organs and tissues, may be able to recapitulate specific features of the in vivo inflammatory response. This minireview provides an overview of MPS for modeling lymphatic tissues, immunity at tissue interfaces, inflammatory diseases, and the inflammatory tumor microenvironment in vitro and ex vivo. Broadly, these systems have utility in modeling how certain immunotherapies function in vivo, how dysfunctional immune responses can propagate diseases, and how our immune system can combat pathogens.

Research & Experimental Medicine↗

DOE OSTI2023

Empirical Study of the Effect of Thermal Loading on the Heating Efficiency of Variable-Speed Air Source Heat Pumps

Heating buildings with air source heat pumps (ASHPs) has the potential to save energy compared to utilizing conventional heat sources. Accurate understanding of the efficiency of ASHPs is important to maximize the energy savings. While it is well understood that, in general, ASHP efficiency decreases with decreasing outdoor temperature, it is not well understood how the ASHP efficiency changes with different levels of thermal loading, even though it is an important consideration for sizing and controlling ASHPs. The goal of this study was to create an empirical model of the ASHP efficiency as a function of two independent variables–outside temperature and level of thermal loading. Four ductless mini-split ASHPs were evaluated in a cold chamber where the temperature (representing the outdoor temperature) was varied over a wide range. For each temperature, the ASHP performance data were collected at several levels of thermal loading. The data for all four ASHPs were combined and approximated with an analytical function that can be used as a general model for the ASHP steady-state efficiency as a function of the outside temperature and level of thermal loading. To the knowledge of the authors, no such empirical model that is solely based on third-party test data has been published before. While limitations exist, the model can be used to help guide future selection and operation of ASHPs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

The Thomas theorem and the Efimov States within a generalized Bohr model

Using a generalized Bohr model and the hyper-spherical formalism for a three-body system, we derive the universal energy functions. We apply our model to 4He atom where only the Coulomb potential is dominant and find that its binding energy is well reproduced (less than 5.5% difference). Later, we focus on the equal mass three-body systems and derive the Thomas theorem assuming a simple interaction depending on the range of the potential. We discuss the conditions for which an unbound two-body system produces a bound three-body system and apply our model to 4He and triton atoms as well as to the triton nucleus. Using their scattering lengths and effective ranges, we are able to reproduce the two-body or the three-body binding energies (less than 5% difference) with only one parameter fitted. Prediction for excited (Efimov) levels are also given and in particular we demonstrate that for some hyper-angles two equal minima appear which indicate a phase (shape) transition similar to the Landau's theory of phase transition. We suggest that the observed excited levels in two different experiments for the triton nucleus are indeed Efimov levels and there may be more surprises.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Fault Contribution of Grid-Following and Grid-Forming Inverters Considering Generic Inverter Controls and Ride-through Requirements

Gradual transformation of the synchronous generator (SG)-dominated power systems to inverter-based resource (IBR)-dominated power systems is bringing new challenges regarding power system protection. Increasing IBR’s penetration will reduce the reliability of protection coordination, originally developed for SG-dominated power systems, potentially leading to less accurate calculation or even mis-operation. Thus, new protection schemes compatible with IBR-dominated power systems should be established. Prior to establishing the protection schemes, research on the fault contribution of generic IBR models should be conducted considering interconnection requirements. This paper conducts fault current contribution analysis of generic grid-following (GFL) and grid-forming (GFM) inverters considering interconnection requirements using an electromagnetic transient (EMT) simulation tool. To do this, EMT-based generic GFL models have been developed with the ride-through requirements from IEEE-1547 standards. In addition, the major functions of the generic models are explained. Then, the fault current analysis was conducted using GFL and GFM inverters in a single-inverter-infinite-bus system and a modified IEEE 34-node system considering different grid, fault, inverter control parameters.

Kim, Jinho [BATTELLE (PACIFIC NW LAB)]↗

DOE OSTI2024

Laboratory Upgrade Point Absorber WEC-Sim Model with MoorDyn Moorings

This dataset includes a WEC-Sim and MoorDyn model of the Laboratory Upgrade Point Absorber (LUPA). LUPA is an open-source wave energy converter designed and tested by Oregon State University. The files provided here constitute a stable LUPA configuration with three mooring lines. This model is 1/20 scale, optimized for the O.H. Hinsdale Wave Lab at Oregon State University. This model of LUPA adds MoorDyn functionality for more accurate mooring predictions and uses a more stable, updated version of LUPA's current physical configuration. This model is for WEC-Sim Version 6.0. A recent update of WEC-Sim has changed some functionality of MoorDyn such that this model will not work with WEC-Sim Version 6.1.

16 TIDAL AND WAVE POWER↗

DOE OSTI2020

HPC-enabled computation of demand models at scale

The purpose of this project is to examine the energy impact of urban-scale traffic for the Los Angeles Basin by developing and implementing a scalable traffic assignment model. An energy optimization function will be posed and when integrated into the optimization code for travel assignment it can be mathematically proven to converge. The energy optimization function can then be compared to the typical travel time optimization that is traditionally used in traffic assignment models. The analysis will begin with static traffic assignment models with the routing for all origin and destinations computed in parallel on high performance computing facilities. Convergence of the numerical methods rely on the solution of convex programs (or extensions of these). This step will mostly consist of demonstrating the ability to parallelize the Frank Wolfe algorithm on various platforms.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2020

HPC4Mobilty w/ UCB

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗