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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Design of model-reference adaptive control systems using Liapunov functions.

This paper considers the design of model-reference adaptive control systems using Liapunov functions. An adaption rule is developed analytically. This adaption rule results from the use of a Liapunov function which contains several positive semidefinite terms not included in previous references. The inclusion of these terms results in an adaption rule which is a function of the error, the plant states and the derivative and integral of these quantities. The adaption rule is applied to a simplified model of the Space Shuttle vehicle. Simulation results show that the maximum response error in the system is reduced by using this adaption rule.

Sutherlin, D. W.↗
DOE OSTI2022

A Framework for Inverse Prediction Using Functional Response Data

Inverse prediction models have commonly been developed to handle scalar data from physical experiments. However, it is not uncommon for data to be collected in functional form. When data are collected in functional form, it must be aggregated to fit the form of traditional methods, which often results in a loss of information. For expensive experiments, this loss of information can be costly. In this study, we introduce the functional inverse prediction (FIP) framework, a general approach which uses the full information in functional response data to provide inverse predictions with probabilistic prediction uncertainties obtained with the bootstrap. The FIP framework is a general methodology that can be modified by practitioners to accommodate many different applications and types of data. We demonstrate the framework, highlighting points of flexibility, with a simulation example and applications to weather data and to nuclear forensics. Results show how functional models can improve the accuracy and precision of predictions.

42 ENGINEERING↗
DOE OSTI2024

Potential reductions in fine particulate matter and premature mortality following implementation of air pollution controls on coal-fired power plants in India

Coal-fired power plants (CFPPs) account for > 70% of electricity generation in India, but < 5% of facilities have installed technologies for sulfur dioxide (SO 2 ) and nitrogen oxide (NO X ) removal. Emissions of these pollutants lead to the formation of fine particulate matter (PM 2.5 ) and an increased risk of premature mortality for exposed populations. Here, we use a nested version of the GEOS-Chem global chemical transport model (0.5° × 0.625° resolution) for India to estimate reductions in PM 2.5 concentrations that could have been achieved by implementing existing emission control technologies like flue-gas desulfurization (FGD) and/or selective catalytic reduction (SCR). We quantify the associated burden of disease using the integrated exposure response (IER) and global exposure mortality model (GEMM) functions and compare the costs of premature mortality to those for FGD installation. Model simulations for 2010 suggest installation of FGD would have reduced mean annual PM 2.5 concentrations across India by 8%, compared to 3% with SCR installation, and 11% with both FGD and SCR. A 7–28% reduction in PM 2.5 was simulated for local communities closest to CFPPs (same model grid cell), leading to up to 17% reduction in annual premature mortality. Overall, more than 0.21–0.48 million premature deaths would have been avoided over a 10-year period if FGD had been implemented on all CFPPs, compared to 0.09–0.21 million with SCR and 0.22–0.72 million with both FGD and SCR. Benefits associated with such actions are approximately $\$18.1$–$\$604$ billion USD per year, which is equivalent to ~ 0.44 to 10% of India’s GDP. These results suggest that monetary benefits from avoided premature mortality far outweigh the capital and operational costs of FGD and/or SCR installation of $\$19.5$ billion and/or $\$32.8$ billion per year, respectively. This information is essential because the high costs of installation and operation are often given as reasons for delaying installation and commissioning. Finally, we conclude that policy actions to control air pollution from CFPPs are economically justifiable.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Leveraging Artificial Intelligence to Predict Novel Eutectic Alloys

The goal of this project was to train an artificial neural network (ANN) to predict the fractional composition and melting point of eutectic alloys using fundamental atomic properties as inputs. The fundamental properties considered include atomic number, atomic weight, atomic radius, valence electron concentration, electronegativity, and electron affinity. The project involved several phases, starting with data preparation, where phase diagram data was harvested from the ASM International database. Approximately 1300 binary eutectics were collected and cleaned to ensure relevance and accuracy. A regression model was selected for training, utilizing a rectified linear unit as the activation function. Various model configurations were evaluated for predictive accuracy, with validation techniques employed to ensure robustness. The model demonstrated predictive capabilities above random guessing and was able to achieve up to 11% accuracy under certain conditions. An ablative test identified atomic radius and valence electron concentration as critical inputs for model performance. Incorporating the melting point of atomic constituents improved accuracy significantly, although ultimately the model’s predictive capability still fell short of the 80% target. This report details the methodology, results, and implications of the research, contributing to the understanding of employing artificial intelligence to predict the phase transition behavior of eutectic alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1993

A Reynolds stress model for near-wall turbulence

The paper formulates a tensorially consistent near-wall second-order closure model. Redistributive terms in the Reynolds stress equations are modeled by an elliptic relaxation equation in order to represent strongly nonhomogeneous effects produced by the presence of walls; this replaces the quasi-homogeneous algebraic models that are usually employed, and avoids the need for ad hoc damping functions. The model is solved for channel flow and boundary layers with zero and adverse pressure gradients. Good predictions of Reynolds stress components, mean flow, skin friction, and displacement thickness are obtained in various comparisons to experimental and direct numerical simulation data. The model is also applied to a boundary layer flowing along a wall with a 90-deg, constant-radius, convex bend.

Durbin, P. A.↗
NASA NTRS1999

Computational Models of the Cardiovascular System and Its Response to Microgravity

Computational models of the cardiovascular system are powerful adjuncts to ground-based and in-flight experiments. We will provide NSBRI with a model capable of simulating the short-term effects of gravity on cardiovascular function. The model from this project will: (1) provide a rational framework which quantitatively defines interactions among complex cardiovascular parameters and which supports the critical interpretation of experimental results and testing of hypotheses. (2) permit predictions of the impact of specific countermeasures in the context of various hypothetical cardiovascular abnormalities induced by microgravity. Major progress has been made during the first 18 months of the program: (1) We have developed an operational first-order computer model capable of simulating the cardiovascular response to orthostatic stress. The model consists of a lumped parameter hemodynamic model and a complete reflex control system. The latter includes cardiopulmonary and carotid sinus reflex limbs and interactions between the two. (2) We have modeled the physiologic stress of tilt table experiments and lower body negative pressure procedures (LBNP). We have verified our model's predictions by comparing them with experimental findings from the literature. (3) We have established collaborative efforts with leading investigators interested in experimental studies of orthostatic intolerance, cardiovascular control, and physiologic responses to space flight. (4) We have established a standardized method of transferring data to our laboratory from the ongoing NSBRI bedrest studies. We use this data to estimate input parameters to our model and compare our model predictions to actual data to further verify our model. (5) We are in the process of systematically simulating current hypotheses concerning the mechanism underlying orthostatic intolerance by matching our simulations to stand test data from astronauts pre- and post-flight. (6) We are in the process of developing a JAVA version of the simulator which will be distributed amongst the cardiovascular team members. Future work on this project involves modifications of the model to represent a rodent (rat) model, further evaluation of the bedrest astronaut and animal data, and systematic investigation of specific countermeasures.

Kamm, Roger D.↗
NASA NTRS1978

The clouds of Venus - An anisotropic scattering model

The phase function for anisotropic scattering is applied to a homogeneous atmospheric model to ascertain the effects of anisotropy on the near-infrared spectrum of Venus. Young's (1971) equivalent widths for the 7820-A CO2 band are analyzed to derive allowed combinations of CO2 specific abundance, continuum albedo, pressure, and degree of anisotropy. From these combinations, values are derived for the photon mean free path and total optical thickness of the clouds. The Venera 10 measurement of the transmitted flux at 7200 A is then used in conjunction with spherical-albedo requirements to reduce the range of possible solutions to the equivalent-width analysis. Through the application of approximate similarity relations, the 'true' (i.e., anisotropic) atmospheric parameters are derived from the isotropic values. For an assumed Mie scatterer with an average anisotropy factor of 0.7, the results indicate a total optical thickness of about 55 with a single-scattering albedo of 0.9992 to 0.9995.

Lestrade, J. P.↗
NASA NTRS2016

A Method for Improving Hotspot Directional Signatures in BRDF Models Used for MODIS

The semi-empirical, kernel-driven, linear RossThick-LiSparseReciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is used to generate the routine MODIS BRDFAlbedo product due to its global applicability and the underlying physics. A challenge of this model in regard to surface reflectance anisotropy effects comes from its underestimation of the directional reflectance signatures near the Sun illumination direction; also known as the hotspot effect. In this study, a method has been developed for improving the ability of the RTLSR model to simulate the magnitude and width of the hotspot effect. The method corrects the volumetric scattering component of the RTLSR model using an exponential approximation of a physical hotspot kernel, which recreates the hotspot magnitude and width using two free parameters (C(sub 1) and C(sub 2), respectively). The approach allows one to reconstruct, with reasonable accuracy, the hotspot effect by adjusting or using the prior values of these two hotspot variables. Our results demonstrate that: (1) significant improvements in capturing hotspot effect can be made to this method by using the inverted hotspot parameters; (2) the reciprocal nature allow this method to be more adaptive for simulating the hotspot height and width with high accuracy, especially in cases where hotspot signatures are available; and (3) while the new approach is consistent with the heritage RTLSR model inversion used to estimate intrinsic narrowband and broadband albedos, it presents some differences for vegetation clumping index (CI) retrievals. With the hotspot-related model parameters determined a priori, this method offers improved performance for various ecological remote sensing applications; including the estimation of canopy structure parameters.

airborne measurements↗
NASA NTRS1988

Shot noise cross-correlation functions and cross spectra - Implications for models of QPO X-ray sources

The cross-correlation functions (CCFs) and cross spectra expected for quasi-periodic oscillation (QPO) shot noise models are calculated under various assumptions, and the results are compared to observations. Effects due to possible coherence of the QPO oscillations are included. General formulas for the cross spectrum, the cross-phase spectrum, and the time-delay spectrum for QPO shot models are calculated and discussed. It is shown that the CCFs, cross spectra, and power spectra observed for Cyg X-e2 imply that the spectrum of the shots evolves with time, with important implications for the interpretation of these functions as well as of observed average energy spectra. The possible origins for the observed hard lags are discussed, and some physical difficulties for the Comptonization model are described. Classes of physical models for QPO sources are briefly addressed, and it is concluded that models involving shot formation at the surface of neutron stars are favored by observation.

Shibazaki, N.↗
DOE OSTI2025

In vivo mapping of mutagenesis sensitivity of human enhancers

Distant-acting enhancers are central to human development1. However, our limited understanding of their functional sequence features prevents the interpretation of enhancer mutations in disease2. Here we determined the functional sensitivity to mutagenesis of human developmental enhancers in vivo. Focusing on seven enhancers that are active in the developing brain, heart, limb and face, we created over 1,700 transgenic mice for over 260 mutagenized enhancer alleles. Systematic mutation of 12-base-pair blocks collectively altered each sequence feature in each enhancer at least once. We show that 69% of all blocks are required for normal in vivo activity, with mutations more commonly resulting in loss (60%) than in gain (9%) of function. Using predictive modelling, we annotated critical nucleotides at the base-pair resolution. The vast majority of motifs predicted by these machine learning models (88%) coincided with changes in in vivo function, and the models showed considerable sensitivity, identifying 59% of all functional blocks. Taken together, our results reveal that human enhancers contain a high density of sequence features that are required for their normal in vivo function and provide a rich resource for further exploration of human enhancer logic.

Kosicki, Michael↗
DOE OSTI2024

Functionally Assembled Terrestrial Ecosystem Simulator (FATES) for Hurricane Disturbance and Recovery

Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post-hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES) to study the processes and the key factors regulating post-hurricane forest recovery. We implemented hurricane-induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post-hurricane forest recovery, and quantified the sensitivity of key model parameters to the post-hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre-hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post-hurricane forest recovery in simulations. This research improves understanding of the ELM-FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

EvoProtGrad (Directed Evolution for Proteins with Gradients) [SWR-23-48]

A Python package for directed evolution on a protein sequence with gradient-based discrete Markov chain monte carlo (MCMC). Users are able to compose custom models that map sequence to function with pretrained models, including protein language models (PLMs), to guide and constrain search. Our package natively integrates with the HuggingFace platform and supports PLMs from transformers. Our MCMC sampler identifies promising amino acids to mutate via model gradients taken with respect to the input (i.e., sensitivity analysis). We allow users to compose their own custom target function for MCMC by leveraging the Product of Experts MCMC paradigm. Each model is an "expert" that contributes its own knowledge about the protein's fitness landscape to the overall target function. The sampler is designed to be more efficient and effective than brute force and random search while maintaining most of the generality and flexibility. Additional information can be found in the related publication: https://iopscience.iop.org/article/10.1088/2632-2153/accacd

Emami, Patrick↗
DOE OSTI2021

Verification Testing For Solid-Element Material Models 11-19 in DYNA3D/ParaDyn

This technical report documents the creation and implementation of verification tests for solid-element material models 11 through 19 available in DYNA3D/ParaDyn. The verification tests covered all aspects of each material model, except for the Weibull distribution functionality in material models 15 and 19. General test cases were created to verify the elastic and plastic behavior of the material models. Other additional tests were developed to examine the intricacies of each material model. Each test involved the use a kinematic load case and specification of material parameters necessary to activate corresponding features of the material model. The load cases prescribed the full time history of the kinematic motion for the solid elements, and these loads are independent of the material model or element formulation. When possible, closed form solutions were then derived for each verification test in a continuum setting. The DYNA3D simulations for each test were carried out over a unit time interval, t ϵ [0, 1], and the as implemented DYNA3D response was compared to the closed form solutions evaluated at discrete points in time. A relative error measure was determined for each test to justify the proper implementation of the material model. The relative errors comparing the DYNA3D solution to the analytical solution were, in general, on the order of machine precision except where noted. This signifies the proper implementation of solid-element material models 11-19. In the development of these verification tests, six bugs were found and fixed in the source code. Additionally, this work generated eighteen DYNA3D input decks and answer extraction scripts in the DYNA3D/ParaDyn Software Quality Assurance test suite, which are comprised of a total of 285 solid-element tests. Testing for each material model utilizes two input decks and answer extraction scripts, where one focuses on the linear elastic response and the other examines the inelastic and remaining functionalities of the material model. In total, this work added 285 individual verification test problems in the DYNA3D/ParaDyn test suite.

42 ENGINEERING↗
DOE OSTI2023

Prediction of nonlinear specific heat during single crystal HMX phase transition

Here, we develop a thermodynamically consistent chemo-thermo-mechanical model for the β → δ phase transition of energetic HMX crystals. In contrast to previous models, which either considered specific heat to be a constant or utilized a calibrated function, this model provides novel expressions for the specific heats at constant volume and constant elastic strains derived directly from continuum mechanics. In addition, the model provides a novel prediction for the critical temperature at which the chemical heating rate achieves its extremum for Arrhenius kinetics. The numerical solution predicts highly nonlinear specific heat behavior including order of magnitude changes.

42 ENGINEERING↗
DOE OSTI2022

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗
DOE OSTI2023

RouteE-Powertrain [SWR-19-19]

RouteE-Powertrain is a tool for predicting energy usage over a set of road links. RouteE-Powertrain is a Python package that allows users to work with a set of pre-trained mesoscopic vehicle energy prediction models for a varity of vehicle types. Additionally, users can train their own models if "ground truth" energy consumption and driving data are available. RouteE-Powertrain models predict vehicle energy consumption over links in a road network, so the features considered for prediction often include traffic speeds, road grade, turns, etc. The typical user will utilize RouteE's catalog of pre-trained models. Currently, the catalog consists of light-duty vehicle models, including conventional gasoline, diesel, hybrid electric (HEV), and battery electric (BEV). These models can be applied to link-level driving data (in the form of pandas dataframes) to output energy consumption predictions. Users that wish to train new RouteE models can do so. The model training function of RouteE enables users to use their own drive-cycle data, powertrain modeling system, and road network data to train custom models. https://pypi.org/project/nrel.routee.powertrain/ pip install nrel.routee.powertrain

Holden, Jacob↗
NASA NTRS1991

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed in this work. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional CFD flow solver with the Monte Carlo simulation of the pdf evolution equation has been developed. The algorithm has been validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames has been carried out using this algorithm. Numerical results compared favorably with experimental data. The computuations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗
NASA NTRS1991

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional computational fluid dynamic (CFD) flow solver with the Monte Carlo simulation of the pdf evolution equation was developed. The algorithm was validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames was carried out using this algorithm. Numerical results compared favorably with experimental data. The computations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗