Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model assignment”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

A Systematic Methodology for Modeling and Attitude Control of Multi-body Space Telescopes

This paper derives a symbolic multi-body rigid nonlinear model for a space telescope using Stoneking’s implementation of Kane’s method. This symbolic nonlinear model is linearized using Matlab symbolic functions diff and inv because the analytic linearization is intractable for manual derivation. The linearized system model is then used to design the controllers using both linear quadratic regulator (LQR) and robust pole assignment methods. The closed-loop systems for the two designs are simulated using both the rigid model as well as a second model containing flexible modes. The performances of the two designs are compared based on the simulation testing results. Our conclusion is that the robust pole assignment design offers better performance than that of the LQR system in terms of actuator usage and pointing accuracy. However, the LQR approach remains an effective first design step that can inform the selection of real eigenvalues for robust pole assignment. The proposed method may be used for the modeling and controller designs for various multi-body systems.

Modeling↗
DOE OSTI2020

Structural modeling and validation of laminated stacks in magnetic gearing applications

Laminated structures are commonly employed in magnetic gears based on their ability to direct magnetic flux, while inhibiting eddy currents. Unfortunately, the stiffness of the laminated structure is much lower than a solid part, which can lead to unacceptable deflections, given the large magnetic forces and small air gaps in magnetic gears. Finite element analysis of laminations to accurately predict these deflections can be computationally expensive due to the small elements needed to model the thin lamination dimensions. These analyses also have some uncertainty due to the assignment of bonding stiffness. This paper presents a means of accurately modelling the deflection of a laminated stack with reduced complexity and computational time, which is compared to experimental results. Here, the resulting methods are then used to assess the magnetic and mechanical advantages of different designs for the modulating pieces of a magnetic gear.

42 ENGINEERING↗
DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗
DOE OSTI2022

Use of vibrational spectroscopy to identify the formation of neptunyl–neptunyl interactions: a paired density functional theory and Raman spectroscopy study

Actinyl–Actinyl interactions (AAIs) occur in pentavalent actinide systems, particularly for neptunium (Np), and lead to complex vibrational signals that are challenging to analyze and interpret. Previous studies have focused on neptunyl–neptunyl dimeric species, but trimers and tetramers have been identified as the primary motif for extended topologies observed in solid-state materials. Our hypothesis is that trimeric and tetrameric AAIs lead to the additional signals in the vibrational spectra, but this has yet to be explored systematically. Herein, we investigate three different neptunyl–neptunyl subunits (dimeric, trimeric, tetrameric) and determine the vibrational frequencies of the O=Np=O stretches using both computational and experimental approaches. Density Functional Theory (DFT) was used to identify distinct vibrational motions related to specific neptunyl oligomers and compared to previous literature precedent from Np(V) in HClO 4 and HCl systems. The vibrational behavior of Np(V) in HNO 3 was then evaluated via Raman spectroscopy. As the solution evaporated signals were linked to trimeric and tetrameric models. Solid phases produced in the evaporation include (NpO 2 ) 2 (NO 3 ) 2 (H 2 O) 5 and newly identified crystalline phase, Na(NpO 2 )(NO 3 ) 2 ·4H 2 O (NpNa). Furthermore, the combined computational studies and vibrational analysis provide evidence for unique observable vibrational bands for each polymerized subunit, allowing us to assign spectral features to trimeric and tetrameric models within three different simple anionic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗
NASA NTRS1987

An investigation of stellar coronae with AXAF

Questions concerning the phenomenon of stellar coronae are discussed together with the physical parameters essential to model stellar coronae and the different AXAF instruments assigned for measurements of these parameters. As an illustration of what AXAF will be able to accomplish, three specific examples of studies by AXAF are described. These include coronal modeling for a bright X-ray source such as HR 1099, using spectra obtained by AXAF; time-resolved spectroscopy of AR Lacertae and YY Geminormum during their eclipses, to use the spectra for the identification and modeling of individual geometrical physical structures on these stars; and time-resolved spectroscopy during stellar flares to obtain high-resolution spectra which will make it possible to infer the changes in the flaring plasma temperatures as a function of time.

Linsky, Jeffrey L.↗
NASA NTRS1972

On the possibility of explaining features of the latitude variations of Sq(H) and Sq(Z) in the equatorial region by electrical modelling

Curves of H (the horizontal component) and Z (the vertical component) of the magnetic field from a strip with an assigned current-density distribution compiled from the results of electrical modelling are calculated. The calculated curves are compared qualitatively with latitude curves obtained from ground magnetic observatory data. It is shown that the forms of the theoretical H and Z and the experimental latitude variations are qualitatively similar.

Pogrebnoy, V. N.↗
DOE OSTI2024

CovTransformer: A transformer model for SARS-CoV-2 lineage frequency forecasting

With hundreds of SARS-CoV-2 lineages circulating in the global population, there is an ongoing need for predicting and forecasting lineage frequencies and thus identifying rapidly expanding lineages. Accurate prediction would allow for more focused experimental efforts to understand pathogenicity of future dominating lineages and characterize the extent of their immune escape. Here, we first show that the inherent noise and biases in lineage frequency data make a commonly-used regression-based approach unreliable. To address this weakness, we constructed a machine learning model for SARS-CoV-2 lineage frequency forecasting, called CovTransformer, based on the transformer architecture. We designed our model to navigate challenges such as a limited amount of data with high levels of noise and bias. We first trained and tested the model using data from the UK and the USA, and then tested the generalization ability of the model to many other countries and US states. Remarkably, the trained model makes accurate predictions two months into the future with high levels of accuracy both globally (in 31 countries with high levels of sequencing effort) and at the US-state level. Our model performed substantially better than a widely used forecasting tool, the multinomial regression model implemented in Nextstrain, demonstrating its utility in SARS-CoV-2 monitoring. Assuming a newly emerged lineage is identified and assigned, our test using retrospective data shows that our model is able to identify the dominating lineages 7 weeks in advance on average before they became dominant. Overall, our work demonstrates that transformer models represent a promising approach for SARS-CoV-2 forecasting and pandemic monitoring.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

General Protocol for the Accurate Prediction of Molecular 13 C/ 1 H NMR Chemical Shifts via Machine Learning Augmented DFT

An accurate prediction of NMR chemical shifts at affordable computational cost is very important for different types of structural assignments in experimental studies. Density functional theory (DFT) and gauge-including atomic orbital (GIAO) are two of the most popular computational methods for NMR calculation, yet they often fail to resolve ambiguities in structural assignments. In this work, we present a new method that uses machine learning (ML) techniques (DFT + ML) that significantly increases the accuracy of 13 C/ 1 H NMR chemical shift prediction for a variety of organic molecules. The input of the generalizable DFT + ML model contains two critical parts: one is a vector providing insights into chemical environments, which can be evaluated without knowing the exact geometry of the molecule; the other one is the DFT-calculated isotropic shielding constant. The DFT + ML model was trained with a data set containing 476 13 C and 270 1 H experimental chemical shifts. For the DFT methods used here, the root mean square deviations (RMSDs) for the errors between predicted and experimental 13 C/ 1 H chemical shifts can be as small as 2.10/0.18 ppm, which is much lower than those from simple DFT (5.54/0.25 ppm), or DFT + linear regression (LR) (4.77/0.23 ppm) approaches. It also has a smaller maximum absolute error than two previously proposed NMR-predicting ML models. The robustness of the DFT + ML model is tested on two classes of organic molecules (TIC10 and hyacinthacines), where the correct isomers were unambiguously assigned to the experimental ones. Overall, the DFT + ML model shows promise for structural assignments in a variety of systems, including stereoisomers, that are often challenging to determine experimentally.

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

Machine Learning Based Resilience Testing of an Address Randomization Cyber Defense

Moving target defenses (MTDs) are widely used as an active defense strategy for thwarting cyberattacks on cyber-physical systems by increasing diversity of software and network paths. Recently, machine Learning (ML) and deep Learning (DL) models have been demonstrated to defeat some of the cyber defenses by learning attack detection patterns and defense strategies. It raises concerns about the susceptibility of MTD to ML and DL methods. Here, in this article, we analyze the effectiveness of ML and DL models when it comes to deciphering MTD methods and ultimately evade MTD-based protections in real-time systems. Specifically, we consider a MTD algorithm that periodically randomizes address assignments within the MIL-STD-1553 protocol—a military standard serial data bus. Two ML and DL-based tasks are performed on MIL-STD-1553 protocol to measure the effectiveness of the learning models in deciphering the MTD algorithm: 1) determining whether there is an address assignments change i.e., whether the given system employs a MTD protocol and if it does 2) predicting the future address assignments. The supervised learning models (random forest and k-nearest neighbors) effectively detected the address assignment changes and classified whether the given system is equipped with a specified MTD protocol. On the other hand, the unsupervised learning model (K-means) was significantly less effective. The DL model (long short-term memory) was able to predict the future addresses with varied effectiveness based on MTD algorithm's settings.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1993

Features of a SINDA/FLUINT model of a liquid oxygen supply line

The modeling features used in a steady-state heat transfer problem using SINDA/FLUINT are described. The problem modeled is a 125 feet long, 3 inch diameter pipe, filled with liquid oxygen flow driven by a given pressure gradient. The pipe is fully insulated in five sections. Three sections of 1 inch thick spray-on foam and two sections of vacuum jacket. The model evaluates friction, turns losses and convection heat transfer between the fluid and the pipe wall. There is conduction through the foam insulation with temperature dependent thermal conductivity. The vacuum space is modeled with radiation and gas molecular conduction, if present, in the annular gap. Heat is transferred between the outer surface and surrounding ambient by natural convection and radiation; and, by axial conduction along the pipe and through the vacuum jacket spacers and welded seal flanges. The model makes extensive use of SINDA/FLUINT basic capabilities such as the GEN option for nodes and conductors (to generate groups of nodes or conductors), the SIV option (to generate single, temperature varying conductors), the SIM option (for multiple, temperature varying conductors) and the M HX macros for fluids (to generate strings of lumps, paths, and ties representing a diabatic duct). It calls subroutine CONTRN (returns the relative location in the G-array of a network conductor, given an actual conductor number) enabling an extensive manipulation of conductor (calculates an assignment of their values) with DO loops. Models like this illustrate to the new and even to the old SINDA/FLUINT user, features of the program that are not so obvious or known, and that are extremely handy when trying to take advantage of both, the automation of the DATA headers and make surgical modifications to specific parameters of the thermal or fluid elements in the OPERATIONS portion of the model.

Simmonds, Boris G.↗
DOE OSTI2022

Hidden vibronic and excitonic structure and vibronic coherence transfer in the bacterial reaction center

We report two-dimensional electronic spectroscopy (2DES) experiments on the bacterial reaction center (BRC) from purple bacteria, revealing hidden vibronic and excitonic structure. Through analysis of the coherent dynamics of the BRC, we identify multiple quasi-resonances between pigment vibrations and excitonic energy gaps, and vibronic coherence transfer processes that are typically neglected in standard models of photosynthetic energy transfer and charge separation. We support our assignment with control experiments on bacteriochlorophyll and simulations of the coherent dynamics using a reduced excitonic model of the BRC. We find that specific vibronic coherence processes can readily reveal weak exciton transitions. While the functional relevance of such processes is unclear, they provide a spectroscopic tool that uses vibrations as a window for observing excited state structure and dynamics elsewhere in the BRC via vibronic coupling. Vibronic coherence transfer reveals the upper exciton of the “special pair” that was weakly visible in previous 2DES experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Computational investigation of hydrogen-induced phonon changes in carbon fiber

Optical vibrational spectroscopy has shown promise as a noninvasive means of monitoring the mechanical properties of carbon fiber (CF), which is increasingly used for industrial and consumer purposes. However, interpretation of optical vibrational spectra for solid materials is inferential, particularly when defects are present. Because inelastic neutron scattering (INS) spectroscopy is not subject to selection rules, the full vibrational spectra can be measured. And, identifying correlations between INS features and tensile properties can assist in the interpretation of spectra from more commonly used optical vibrational spectroscopic techniques, such as Raman and infrared (IR) spectroscopy. Recent INS experiments on high-performance commercial carbon fibers showed features near 900 and 1100 cm –1 in addition to a broad feature near 3000 cm –1 that increased in intensity with decreasing tensile strength. These features were assigned to hydrogen defects. Here, we use density functional theory to simulate the INS spectra of several hydrogen defect geometries in graphite as a model for carbon fiber structure units, confirming the experimental assignment of these peaks to hydrogen modes and providing insights into the structure and lattice dynamics of the defects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibration of cloud and aerosol related parameters for solar irradiance forecasts in WRF-solar

Model parameters are a major source of uncertainty in numerical weather prediction. Recently, the Weather Research and Forecasting model with Solar extensions (WRF-Solar) has been upgraded by enhancing the treatment of sub-grid scale cloud and aerosols with augmentations of a sub-grid scale cloud scheme (CLD3) and an upgraded aerosol-aware Thompson-Eidhammer scheme (TE14). However, the value of model parameters associated with these parameterizations are assigned based on limited measurements or theoretical calculations. Calibrating the most sensitive parameters has the potential to improve solar irradiance predictions. Here, we adopted a multiobjective surrogate-based optimization (SBO) framework to calibrate nine parameters used in CLD3 and TE14 that lead to the largest sensitivity in simulated irradiance. The normalized mean-absolute-error (NMAE) of global horizontal irradiance (GHI) and direct normal irradiance (DNI) are minimized by calibrating WRF-Solar over two regions including the Southern Great Plains (SGP) and Central California, in order to focus on parameter calibration under cloudy conditions with different aerosol loading. The results show that generalized linear model (GLM)-based surrogate models approximate physical models well, particularly when the third order and three-way interaction terms are considered. The SBO framework efficiently searches the parameter space for optimal solutions with less computational costs than directly calibrating the physical model. We first calibrate CLD3 parameters over the less-polluted SGP region. Optimized CLD3 parameters alone result in NMAE reduction by 14% for the site-mean and up to 33% for individual cases over the SGP region. With further calibration of TE14 parameters over the Central California during active fire periods, the optimized parameters lead to over 20% reductions of NMAE. Our investigation reveals, however, that optimizing TE14 has a limited impact on irradiance simulations under less-polluted conditions in the SGP.

14 SOLAR ENERGY↗
NASA NTRS1982

A Markov chain model for reliability growth and decay

A mathematical model is developed to describe a complex system undergoing a sequence of trials in which there is interaction between the internal states of the system and the outcomes of the trials. For example, the model might describe a system undergoing testing that is redesigned after each failure. The basic assumptions for the model are that the state of the system after a trial depends probabilistically only on the state before the trial and on the outcome of the trial and that the outcome of a trial depends probabilistically only on the state of the system before the trial. It is shown that under these basic assumptions, the successive states form a Markov chain and the successive states and outcomes jointly form a Markov chain. General results are obtained for the transition probabilities, steady-state distributions, etc. A special case studied in detail describes a system that has two possible state ('repaired' and 'unrepaired') undergoing trials that have three possible outcomes ('inherent failure', 'assignable-cause' 'failure' and 'success'). For this model, the reliability function is computed explicitly and an optimal repair policy is obtained.

Siegrist, K.↗
NASA NTRS2010

Optimization of Regression Models of Experimental Data Using Confirmation Points

A new search metric is discussed that may be used to better assess the predictive capability of different math term combinations during the optimization of a regression model of experimental data. The new search metric can be determined for each tested math term combination if the given experimental data set is split into two subsets. The first subset consists of data points that are only used to determine the coefficients of the regression model. The second subset consists of confirmation points that are exclusively used to test the regression model. The new search metric value is assigned after comparing two values that describe the quality of the fit of each subset. The first value is the standard deviation of the PRESS residuals of the data points. The second value is the standard deviation of the response residuals of the confirmation points. The greater of the two values is used as the new search metric value. This choice guarantees that both standard deviations are always less or equal to the value that is used during the optimization. Experimental data from the calibration of a wind tunnel strain-gage balance is used to illustrate the application of the new search metric. The new search metric ultimately generates an optimized regression model that was already tested at regression model independent confirmation points before it is ever used to predict an unknown response from a set of regressors.

Ulbrich, N.↗