Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model based development”

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,477 records · Page 82

Detection of DoS Attacks Using ARFIMA Modeling of GOOSE Communication in IEC 61850 Substations

Integration of Information and Communication Technology (ICT) in modern smart grids (SGs) offers many advantages including the use of renewables and an effective way to protect, control and monitor the energy transmission and distribution. To reach an optimal operation of future energy systems, availability, integrity and confidentiality of data should be guaranteed. Research on the cyber-physical security of electrical substations based on IEC 61850 is still at an early stage. In the present work, we first model the network traffic data in electrical substations, then, we present a statistical Anomaly Detection (AD) method to detect Denial of Service (DoS) attacks against the Generic Object Oriented Substation Event (GOOSE) network communication. According to interpretations on the self-similarity and the Long-Range Dependency (LRD) of the data, an Auto-Regressive Fractionally Integrated Moving Average (ARFIMA) model was shown to describe well the GOOSE communication in the substation process network. Based on this ARFIMA-model and in view of cyber-physical security, an effective model-based AD method is developed and analyzed. Two variants of the statistical AD considering statistical hypothesis testing based on the Generalized Likelihood Ratio Test (GLRT) and the cumulative sum (CUSUM) are presented to detect flooding attacks that might affect the availability of the data. Our work presents a novel AD method, with two different variants, tailored to the specific features of the GOOSE traffic in IEC 61850 substations. The statistical AD is capable of detecting anomalies at unknown change times under the realistic assumption of unknown model parameters. The performance of both variants of the AD method is validated and assessed using data collected from a simulation case study. We perform several Monte-Carlo simulations under different noise variances. The detection delay is provided for each detector and it represents the number of discrete time samples after which an anomaly is detected. In fact, our statistical AD method with both variants (CUSUM and GLRT) has around half the false positive rate and a smaller detection delay when compared with two of the closest works found in the literature. Our AD approach based on the GLRT detector has the smallest false positive rate among all considered approaches. Whereas, our AD approach based on the CUSUM test has the lowest false negative rate thus the best detection rate. Depending on the requirements as well as the costs of false alarms or missed anomalies, both variants of our statistical detection method can be used and are further analyzed using composite detection metrics.

IEC 61850 electrical substations↗

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗

Modeling human decision making behavior in supervisory control

An optimal decision control model was developed, which is based primarily on a dynamic programming algorithm which looks at all the available task possibilities, charts an optimal trajectory, and commits itself to do the first step (i.e., follow the optimal trajectory during the next time period), and then iterates the calculation. A Bayesian estimator was included which estimates the tasks which might occur in the immediate future and provides this information to the dynamic programming routine. Preliminary trials comparing the human subject's performance to that of the optimal model show a great similarity, but indicate that the human skips certain movements which require quick change in strategy.

Tulga, M. K.↗

A Coupled Theoretical and Experimental Approach to Elucidating the Mechanisms of Methyl Esters

The goal of this project was to probe fundamental kinetics questions regarding the gas phase reactive behavior of oxygenates. In order to succeed, our program necessitated both the development of new experimental tools and the development of theory based kinetic models to explain the chemistry of oxygenate fuels. By using advanced hybrid additive manufacturing/ traditional manufacturing techniques to create specialized ceramic micro-reactors, the team was able to simulate extreme high-temperature environments with much greater precision and durability than was previously possible. These tiny, high-tech tubes allow researchers to capture and identify "fleeting" chemical species—molecules that exist for only a fraction of a millisecond—using sophisticated light sources and mass spectrometry. Additionally, the team supported the development of a tabletop VUV laser system for isomer detection without a synchrotron. Through this work, we have successfully mapped out the specific chemical pathways of various oxygen-rich fuels, solving long-standing mysteries about how these substances break down kinetically. Ultimately, this research provides the fundamental knowledge needed to design next-generation engines and fuels that are better for the environment.

09 BIOMASS FUELS↗

An integrated experimental-computational approach for predicting virulence in New Zealand white rabbits and humans following inhalation exposure to Bacillus anthracis spores

Inhalation of Bacillus anthracis spores can lead to an anthrax infection that can be fatal. Previously published mathematical models have extrapolated kinetic rates associated with bacterial growth in New Zealand White (NZW) rabbits to humans, but to date, actual measurements of the underlying processes associated with anthrax virulence between species have not been conducted. To address this knowledge gap, we have quantified species-specific rate constants associated with germination, proliferation, and immune cell inactivation of B. anthracis Sterne using an in vitro test platform that includes primary lung epithelial and immune cells. The generated data was then used to develop a physiologically based biokinetic model (PBBK) which quantitatively compares bacterial growth and mean time to death under lethal conditions in rabbits and humans. Simulations based upon our in vitro data and previously published in vivo data from rabbits indicate that disease progression is likely to be faster in humans than in NZW rabbits under comparable total deposited dose conditions. With the computational framework established, PBBK parameters can now be refined using experimental data for lethal B. anthracis strains (e.g. Ames) under identical conditions in future studies. The PBBK model can also be linked to existing aerosol dosimetry models that account for species-specific differences in aerosol deposition patterns to further improve the human health risk assessment of inhalation anthrax.

3-D tissue culture, B. anthracis Sterne, biothreat↗

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed.

Impact↗

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed

Impact↗

A critique of some recent second-order turbulence closure models for compressible boundary layers

Computations based on two recently developed second-order turbulence closure models are compared with a series of boundary-layer experiments and with predictions of these experiments using an algebraic mixing length model. One of the models employs an eddy viscosity, whereas the other evaluates components of the Reynolds stress tensor. For flat plates, the computations are compared with the van Driest skin-friction transformation to assess the handling of compressibility. For boundary layers in pressure gradients, four experiments at Mach 4 and one at Mach 6.7 are used as the bases for comparison. In general, both models represent mean velocities and skin friction reasonably well, but represent the turbulence shear stress less accurately.

Rubesin, M. W.↗

NASA's Quiet Aircraft Technology Project

NASA's Quiet Aircraft Technology Project is developing physics-based understanding, models and concepts to discover and realize technology that will, when implemented, achieve the goals of a reduction of one-half in perceived community noise (relative to 1997) by 2007 and a further one-half in the far term. Noise sources generated by both the engine and the airframe are considered, and the effects of engine/airframe integration are accounted for through the propulsion airframe aeroacoustics element. Assessments of the contribution of individual source noise reductions to the reduction in community noise are developed to guide the work and the development of new tools for evaluation of unconventional aircraft is underway. Life in the real world is taken into account with the development of more accurate airport noise models and flight guidance methodology, and in addition, technology is being developed that will further reduce interior noise at current weight levels or enable the use of lighter-weight structures at current noise levels.

Whitfield, Charlotte E.↗

ROM-Based Surrogate Systems Modeling of EBR-II

We report the System Analysis Module (SAM), developed and maintained by Argonne National Laboratory, is designed to provide whole-plant transient safety analysis capabilities for a number of advanced non-light water reactors, including sodium-cooled fast reactor (SFR), lead-cooled fast reactor (LFR), and molten salt reactor (MSR)/fluoride-salt-cooled high-temperature reactor (FHR) designs. SAM is primarily constructed as a systems-level analysis tool, with the potential to incorporate reduced order models from three-dimensional computational fluid dynamics (CFD) simulations to improve characterization of complex, multidimensional physics. It is recognized that the computational expense associated with CFD can be intractable for various engineering analyses, such as uncertainty quantification, inference, and design optimization. This paper explores the reducibility of a SAM model using recent advances in randomized linear algebra techniques, which attempt to find recurring patterns in the various realizations generated by a model after randomly perturbing all its input parameters. The reduction is described in terms of fewer degrees of freedom (DOFs), referred to as the active DOFs, for the model variables such as input model parameters and model responses. The results indicate that there is significant room for additional reduction that may be leveraged for additional computational gains when employing SAM for engineering-intensive analyses that require repeated model executions. Different from physics-based reduction approaches, the proposed approach allows one to estimate upper bounds on the reduction errors, which are rigorously developed in this work. Finally, different methods for surrogate model construction, such as regression and neural network-based training, are employed to correlate the input and output active DOFs, which are related back to the original variables using matrix-based linear transformations.

42 ENGINEERING↗

Models to study gravitational biology of Mammalian reproduction

Mammalian reproduction evolved within Earth's 1-g gravitational field. As we move closer to the reality of space habitation, there is growing scientific interest in how different gravitational states influence reproduction in mammals. Habitation of space and extended spaceflight missions require prolonged exposure to decreased gravity (hypogravity, i.e., weightlessness). Lift-off and re-entry of the spacecraft are associated with exposure to increased gravity (hypergravity). Existing data suggest that spaceflight is associated with a constellation of changes in reproductive physiology and function. However, limited spaceflight opportunities and confounding effects of various nongravitational factors associated with spaceflight (i.e., radiation, stress) have led to the development of ground-based models for studying the effects of altered gravity on biological systems. Human bed rest and rodent hindlimb unloading paradigms are used to study exposure to hypogravity. Centrifugation is used to study hypergravity. Here, we review the results of spaceflight and ground-based models of altered gravity on reproductive physiology. Studies utilizing ground-based models that simulate hyper- and hypogravity have produced reproductive results similar to those obtained from spaceflight and are contributing new information on biological responses across the gravity continuum, thereby confirming the appropriateness of these models for studying reproductive responses to altered gravity and the underlying mechanisms of these responses. Together, these unique tools are yielding new insights into the gravitational biology of reproduction in mammals.

Review, Tutorial↗

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗

Finite state aeroelastic model for use in rotor design optimization

In this article, a rotor aeroelastic model based on a newly developed finite state dynamic wake, coupled with blade finite element analysis, is described. The analysis is intended for application in rotor blade design optimization. A coupled simultaneous system of differential equations combining blade structural dynamics and aerodynamics is established in a formulation well-suited for design sensitivity computation. Each blade is assumed to be an elastic beam undergoing flap bending, lead-lag bending, elastic twist, and axial deflections. Aerodynamic loads are computed from unsteady blade element theory where the rotor three-dimensional unsteady wake is described by a generalized dynamic wake model. Correlation of results obtained from the analysis with flight test data is provided to assess model accuracy.

He, Chengjian↗

A non-local, ordinary-state-based viscoelasticity model for peridynamics.

A non-local, ordinary-state-based, peridynamics viscoelasticity model is developed. In this model, viscous effects are added to deviatoric deformations and the bulk response remains elastic. The model uses internal state variables and is conceptually similar to linearized isotropic viscolelasticity in the local theory. The modulus state, which is used to form the Jacobian matrix in Newton-Raphson algorithms, is presented. The model is shown to satisfy the 2nd law of thermodynamics and is applicable to problems in solid continuum mechanics where fracture and rate effects are important; it inherits all the advantages for modeling fracture associated with peridynamics. By combining this work with the previously published ordinary-state-based plasticity model, the model may be amenable to viscoplasticity problems where plasticity and rate effects are simultaneously important. Also, the model may be extended to include viscous effects for spherical deformations as well. The later two extensions are not presented and may be the subject of further work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fifteen-minute Extravehicular Activity Prebreathe Protocol Using NASA's Exploration Atmosphere (8.2 psia/ 34% 02)

A TBDM DCS probability model based on an existing biophysical model of inert gas bubble growth provides significant prediction and goodness-of-fit with 84 cases of DCS in 668 human altitude exposures. 2. Model predictions suggest that 15-minute O2 prebreathe protocols used in conjunction with suit ports and an 8.2 psi, 34% O2, 66% N2 atmosphere may enable rapid EVA capability for future exploration missions with the risk of DCS ≤ 12%.  EVA could begin immediately at 6.0 psi, with crewmembers decreasing suit pressure to 4.3 psi after completing the 15-minute in-suit prebreathe. 3. Model predictions suggest that intermittent recompression during exploration EVA may reduce decompression stress by 1.8% to 2.3% for 6 hours of total EVA time. Savings in gas consumables and crew time may be accumulated by abbreviating the EVA suit N2 purge to 2 minutes (20% N2) compared with 8 minutes (5% N2) at the expense of an increase in estimated decompression risk of up to 2.4% for an 8-hour EVA.  Increased DCS risk could be offset by IR or by spending additional time at 6 psi at the beginning of the EVA.  Savings of 0.48 lb of gas and 6 minutes per person per EVA corresponds to more than 31 hours of crew time and 1800 lb of gas and tankage under the Constellation lunar architecture. 6. Further research is needed to characterize and optimize breathing mixtures and intermittent recompression across the range of environments and operational conditions in which astronauts will live and work during future exploration missions. 7. Development of exploration prebreathe protocols will begin with definition of acceptable risk, followed by development of protocols based on models such as ours, and, ultimately, validation of protocols through ground trials before operational implementation.

Abercromby, Andrew F. J.↗

A probabilistic graphical model approach to stochastic multiscale partial differential equations

We develop a probabilistic graphical model based methodology to efficiently perform uncertainty quantification in the presence of both stochastic input and multiple scales. Both the stochastic input and model responses are treated as random variables in this framework. Their relationships are modeled by graphical models which give explicit factorization of a high-dimensional joint probability distribution. The hyperparameters in the probabilistic model are learned using sequential Monte Carlo (SMC) method, which is superior to standard Markov chain Monte Carlo (MCMC) methods for multi-modal distributions. Finally, we make predictions from the probabilistic graphical model using the belief propagation algorithm. Numerical examples are presented to show the accuracy and efficiency of the predictive capability of the developed graphical model.

97 MATHEMATICS AND COMPUTING↗

A Physics-Based Modeling Framework for Prognostic Studies

Prognostics and Health Management (PHM) methodologies have emerged as one of the key enablers for achieving efficient system level maintenance as part of a busy operations schedule, and lowering overall life cycle costs. PHM is also emerging as a high-priority issue in critical applications, where the focus is on conducting fundamental research in the field of integrated systems health management. The term diagnostics relates to the ability to detect and isolate faults or failures in a system. Prognostics on the other hand is the process of predicting health condition and remaining useful life based on current state, previous conditions and future operating conditions. PHM methods combine sensing, data collection, interpretation of environmental, operational, and performance related parameters to indicate systems health under its actual application conditions. The development of prognostics methodologies for the electronics field has become more important as more electrical systems are being used to replace traditional systems in several applications in the aeronautics, maritime, and automotive fields. The development of prognostics methods for electronics presents several challenges due to the great variety of components used in a system, a continuous development of new electronics technologies, and a general lack of understanding of how electronics fail. Similarly with electric unmanned aerial vehicles, electrichybrid cars, and commercial passenger aircraft, we are witnessing a drastic increase in the usage of batteries to power vehicles. However, for battery-powered vehicles to operate at maximum efficiency and reliability, it becomes crucial to both monitor battery health and performance and to predict end of discharge (EOD) and end of useful life (EOL) events. We develop an electrochemistry-based model of Li-ion batteries that capture the significant electrochemical processes, are computationally efficient, capture the effects of aging, and are of suitable accuracy for reliable EOD prediction in a variety of usage profiles.

Li-ion Batteries↗