Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer model”

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,711 records · Page 95

Estimating short-period dynamics using an extended Kalman filter

An extended Kalman filter (EKF) is used to estimate the parameters of a low-order model from aircraft transient response data. The low-order model is a state space model derived from the short-period approximation of the longitudinal aircraft dynamics. The model corresponds to the pitch rate to stick force transfer function currently used in flying qualities analysis. Because of the model chosen, handling qualities information is also obtained. The parameters are estimated from flight data as well as from a six-degree-of-freedom, nonlinear simulation of the aircraft. These two estimates are then compared and the discrepancies noted. The low-order model is able to satisfactorily match both flight data and simulation data from a high-order computer simulation. The parameters obtained from the EKF analysis of flight data are compared to those obtained using frequency response analysis of the flight data. Time delays and damping ratios are compared and are in agreement. This technique demonstrates the potential to determine, in near real time, the extent of differences between computer models and the actual aircraft. Precise knowledge of these differences can help to determine the flying qualities of a test aircraft and lead to more efficient envelope expansion.

Bauer, Jeffrey E.↗

High-altitude balloon-launched aircraft - A piloted simulation study

A real-time piloted simulation at the NASA Dryden Flight Research Facility was used to study the feasibility of launching a research aircraft from a high-altitude balloon. In the study, the simulated aircraft was launched in a nosedown attitude at zero airspeed from 110,000 ft. After launch, the pilot flew the aircraft through a near-maximum-lift pullout and then through a zoom climb to a trimmed, l-g flight condition at the test altitude. The study included parametric variations to measure the effects of launch altitude, gross weight, Mach number limit, and parachute size on the test altitude attained. The aerodynamic model of the simulated aircraft was based on flight test results, low Reynolds number windtunnel tests, and computational models; the model included significant Mach number and Reynolds number effects at high altitude. A small parachute was included in the simulation to limit Mach number during the pullout to avoid adverse transonic effects and their resultant energy losses. A small rocket motor was included in the simulation and was investigated for boosting the aircraft to a higher test altitude. In the study, a test altitude of approximately 95,000 ft was attained without rocket boost, and a test altitude in excess of 100,000 ft was attained using small rocket boost.

Murray, James E.↗

Development and evaluation of a musculoskeletal model of the elbow joint complex

This paper describes the development and evaluation of a musculoskeletal model that represents human elbow flexion-extension and forearm pronation-supination. The length, velocity, and moment arm for each of the eight musculotendon actuators were based on skeletal anatomy and position. Musculotendon parameters were determined for each actuator and verified by comparing analytical torque-angle curves with experimental joint torque data. The parameters and skeletal geometry were also utilized in the musculoskeletal model for the analysis of ballistic elbow joint complex movements. The key objective was to develop a computational model, guided by parameterized optimal control, to investigate the relationship among patterns of muscle excitation, individual muscle forces, and movement kinematics. The model was verified using experimental kinematic, torque, and electromyographic data from volunteer subjects performing ballistic elbow joint complex movements.

Gonzalez, Roger V.↗

A musculoskeletal model of the elbow joint complex

This paper describes a musculoskeletal model that represents human elbow flexion-extension and forearm pronation-supination. Musculotendon parameters and the skeletal geometry were determined for the musculoskeletal model in the analysis of ballistic elbow joint complex movements. The key objective was to develop a computational model, guided by optimal control, to investigate the relationship among patterns of muscle excitation, individual muscle forces, and movement kinematics. The model was verified using experimental kinematic, torque, and electromyographic data from volunteer subjects performing both isometric and ballistic elbow joint complex movements. In general, the model predicted kinematic and muscle excitation patterns similar to what was experimentally measured.

Gonzalez, Roger V.↗

Operationally Efficient Propulsion System Study (OEPSS) Data Book: Integrated Booster Propulsion Module (BPM) Engine Start Dynamics - Volume 8

A fluid-dynamic, digital-transient computer model of an integrated, parallel propulsion system was developed for the CDC mainframe and the SUN workstation computers. Since all STME component designs were used for the integrated system, computer subroutines were written characterizing the performance and geometry of all the components used in the system, including the manifolds. Three transient analysis reports were completed. The first report evaluated the feasibility of integrated engine systems in regards to the start and cutoff transient behavior. The second report evaluated turbopump out and combined thrust chamber/turbopump out conditions. The third report presented sensitivity study results in staggered gas generator spin start and in pump performance characteristics.

Kemp, Victoria R.↗

Magnetic Shield Design Modeling and Validation on SWOT Spacecraft Applied to Mars Flux Pinning Orbiting Sample Design

A modeling methodology was developed for use on the Surface Water and Ocean Topography (SWOT) mission for the computational modeling and design of a magnetic shield for a 63 A-sq m source. Shield options were modeled and tested across various design parameters and validated with measurement. Measurement results fell within 10% of simulation in most cases of concern with sources of error well understood. These methods results informed a subsequent modeling activity for a future Mars Sample Return (MSR) mission concept with stringent magnetic cleanliness requirements to retain geologic integrity of the samples. One option for capturing the sample from orbit requires rare earth magnets in close proximity to the Martian samples. Magnetic modeling with a finite-element method solver estimated the magnetic environment and established the need for shielding. Further modeling then determined the shielding necessary to meet magnetic requirements, followed by the design of a mass-optimized solution.

Gonzales, Edward↗

Ig-VAE: Generative modeling of protein structure by direct 3D coordinate generation

While deep learning models have seen increasing applications in protein science, few have been implemented for protein backbone generation—an important task in structure-based problems such as active site and interface design. We present a new approach to building class-specific backbones, using a variational auto-encoder to directly generate the 3D coordinates of immunoglobulins. Our model is torsion- and distance-aware, learns a high-resolution embedding of the dataset, and generates novel, high-quality structures compatible with existing design tools. We show that the Ig-VAE can be used with Rosetta to create a computational model of a SARS-CoV2-RBD binder via latent space sampling. We further demonstrate that the model’s generative prior is a powerful tool for guiding computational protein design, motivating a new paradigm under which backbone design is solved as constrained optimization problem in the latent space of a generative model.

59 BASIC BIOLOGICAL SCIENCES↗

Cathode Modeling of Solid-State Batteries

The search for safe, reliable, and compact high-capacity energy storage devices has led to increased interest in all-solid-state battery research. The use of solid electrolytes provides enhanced safety and durability due to their reduced flammability and increased mechanical strength compared to organic liquid electrolytes. Still, the use of solid electrolytes remains challenging. Computational modeling plays a substantial role in addressing these challenges. A particle dynamics electromechanical model for simulating electrochemical processes in a solid-state battery cathode will be presented. The model presents cathode microstructure at the particle level as a mixture of ionically conductive solid electrolyte particles, electrically conductive carbon additives, and cathodic reactant particles. After densification, the particle connectivity is analyzed to reconstruct the complex electric network connecting reactant particles with an anodic and cathodic current collectors through the electrolyte and carbon particles. The Kirchhoff’s matrix equation describing this electric network, is solved to obtain values of various critical parameters, such as the overall conductivity of the cathode for lithium ions and electrons, cathodic reactant material utilization, and the distribution of the electric current and voltages within the cathode. In addition, by representing the reactant particles as electrolyte or galvanic microcells governed by the Butler-Volmer electrochemical equation, the overall performance of battery cells during charge or discharge processes, respectively, can be predicted for a given cathodic powder composition. The presented model, executed on a high-performance computing architecture, essentially provides a valuable guidance in designing and developing future solid-state batteries.

solid-state battery↗

Variable-Complexity Multidisciplinary Optimization on Parallel Computers

This report covers work conducted under grant NAG1-1562 for the NASA High Performance Computing and Communications Program (HPCCP) from December 7, 1993, to December 31, 1997. The objective of the research was to develop new multidisciplinary design optimization (MDO) techniques which exploit parallel computing to reduce the computational burden of aircraft MDO. The design of the High-Speed Civil Transport (HSCT) air-craft was selected as a test case to demonstrate the utility of our MDO methods. The three major tasks of this research grant included: development of parallel multipoint approximation methods for the aerodynamic design of the HSCT, use of parallel multipoint approximation methods for structural optimization of the HSCT, mathematical and algorithmic development including support in the integration of parallel computation for items (1) and (2). These tasks have been accomplished with the development of a response surface methodology that incorporates multi-fidelity models. For the aerodynamic design we were able to optimize with up to 20 design variables using hundreds of expensive Euler analyses together with thousands of inexpensive linear theory simulations. We have thereby demonstrated the application of CFD to a large aerodynamic design problem. For the predicting structural weight we were able to combine hundreds of structural optimizations of refined finite element models with thousands of optimizations based on coarse models. Computations have been carried out on the Intel Paragon with up to 128 nodes. The parallel computation allowed us to perform combined aerodynamic-structural optimization using state of the art models of a complex aircraft configurations.

Grossman, Bernard↗

Improved dynamic analysis method using load-dependent Ritz vectors

The dynamic analysis of large space structures is important in order to predict their behavior under operating conditions. Computer models of large space structures are characterized by having a large number of degrees of freedom, and the computational effort required to carry out the analysis is very large. Conventional methods of solution utilize a subset of the eigenvectors of the system, but for systems with many degrees of freedom, the solution of the eigenproblem is in many cases the most costly phase of the analysis. For this reason, alternate solution methods need to be considered. It is important that the method chosen for the analysis be efficient and that accurate results be obtainable. It is important that the method chosen for the analysis be efficient and that accurate results be obtainable. The load dependent Ritz vector method is presented as an alternative to the classical normal mode methods for obtaining dynamic responses of large space structures. A simplified model of a space station is used to compare results. Results show that the load dependent Ritz vector method predicts the dynamic response better than the classical normal mode method. Even though this alternate method is very promising, further studies are necessary to fully understand its attributes and limitations.

Escobedo-Torres, J.↗

An Instrumented Capsule Design to Measure Thermal Conductivity in Miniature UO2 Specimens

Numerous separate effects irradiations of miniature nuclear fuel specimens have been conducted in the High Flux Isotope Reactor (HFIR) under the experimental platform designated as MiniFuel. MiniFuel is a static irradiation capability in which microstructural evolution and fuel performance phenomena are observed during postirradiation examination thereby offering a snapshot of the terminal fuel characteristics. This approach inherently requires fielding an irradiation where the experimental conditions are determined using predictive models and the pertinent outcomes are measured at the end of the test. Static irradiations can provide useful insights to the relationships between fuel performance and the pivotal irradiation conditions, namely temperature and burnup, but the ability to monitor fuel performance in situ would further support fuel development and qualification. To this end, an instrumented experiment design is being developed at Oak Ridge National Laboratory to capture thermal conductivity degradation and fission gas release during HFIR irradiation. These phenomena will be monitored using unique capsule designs that each target a different phenomenon. This paper details the thermal conductivity capsule (TCC) design and its expected performance envelope as determined using computer models. Each TCC will contain a miniature UO2 disc specimen (~0.5 mm thick × 5 mm diameter) sandwiched between metallic slugs with embedded thermocouples. The coupling of in situ temperature measurements, known thermal conductivity of the metallic components, and heat generation rates computed using high-fidelity neutronics models make the thermal conductivity measurement possible. This paper describes the reactor physics and heat transfer models used to predict the capsule’s performance and the methodology for calculating the fuel specimen’s thermal conductivity from the thermocouple measurements.

Gorton, Jacob [ORNL] (ORCID:0000000269806083)↗

Three-dimensional modelling of a self-sustained atmospheric pressure glow discharge

The atmospheric pressure glow discharge (APGD) is a relatively simple and versatile plasma source used in diverse applications. Stable APGD operation at high currents, generally a challenge due to instabilities leading to glow-to-arc transition, has been demonstrated using actively-controlled cathodic cooling. This article presents the computational modelling and simulation of a self-sustained direct-current APGD in helium within a 10 mm pin-to-plate inter-electrode gap for currents ranging from 4 to 40 mA. The APGD model is comprised of the conservation equations for total mass, chemical species, momentum, thermal energy of heavy-species and of free electrons, and electric charge. The model equations are discretized using a nonlinear variational multi-scale finite element method that has demonstrated superior accuracy in other plasma flow problems, on a temporal and three-dimensional computational domain suitable to unveil the potential occurrence of instabilities. Modelling results show good agreement with experimental measurements of voltage drop and the same trend but higher values of temperature. The higher temperatures obtained by the simulations appear to be due to the absence of a near-cathode heat dissipation model. Here, the results also reveal that the distribution of electron density and of the ratio of atomic helium ions to total ions transitions from monotonically increasing away from the cathode to presenting a minimum near the centre of the gap with increasing current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Emergence of Complexity in Protein Functions and Metabolic Networks

In modern organisms proteins perform a majority of cellular functions, such as chemical catalysis, energy transduction and transport of material across cell walls. Although great strides have been made towards understanding protein evolution, a meaningful extrapolation from contemporary proteins to their earliest ancestors is virtually impossible. In an alternative approach, the origin of water-soluble proteins was probed through the synthesis of very large libraries of random amino acid sequences and subsequently subjecting them to in vitro evolution. In combination with computer modeling and simulations, these experiments allow us to address a number of fundamental questions about the origins of proteins. Can functionality emerge from random sequences of proteins? How did the initial repertoire of functional proteins diversify to facilitate new functions? Did this diversification proceed primarily through drawing novel functionalities from random sequences or through evolution of already existing proto-enzymes? Did protein evolution start from a pool of proteins defined by a frozen accident and other collections of proteins could start a different evolutionary pathway? Although we do not have definitive answers to these questions, important clues have been uncovered. Considerable progress has been also achieved in understanding the origins of membrane proteins. We will address this issue in the example of ion channels - proteins that mediate transport of ions across cell walls. Remarkably, despite overall complexity of these proteins in contemporary cells, their structural motifs are quite simple, with -helices being most common. By combining results of experimental and computer simulation studies on synthetic models and simple, natural channels, I will show that, even though architectures of membrane proteins are not nearly as diverse as those of water-soluble proteins, they are sufficiently flexible to adapt readily to the functional demands arising during evolution.

Pohorille, Andzej↗

A performance analysis method for distributed real-time robotic systems: A case study of remote teleoperation

Robot coordination and control systems for remote teleoperation applications are by necessity implemented on distributed computers. Modeling and performance analysis of these distributed robotic systems is difficult, but important for economic system design. Performance analysis methods originally developed for conventional distributed computer systems are often unsatisfactory for evaluating real-time systems. The paper introduces a formal model of distributed robotic control systems; and a performance analysis method, based on scheduling theory, which can handle concurrent hard-real-time response specifications. Use of the method is illustrated by a case of remote teleoperation which assesses the effect of communication delays and the allocation of robot control functions on control system hardware requirements.

Lefebvre, D. R.↗

Structure of 3-mercaptopropionic acid dioxygenase with a substrate analog reveals bidentate substrate binding at the iron center

Thiol dioxygenases are a subset of nonheme iron oxygenases that catalyze the formation of sulfinic acids from sulfhydryl-containing substrates and dioxygen. Among this class, cysteine dioxygenases (CDOs) and 3-mercaptopropionic acid dioxygenases (3MDOs) are the best characterized, and the mode of substrate binding for CDOs is well understood. However, the manner in which 3-mercaptopropionic acid (3MPA) coordinates to the nonheme iron site in 3MDO remains a matter of debate. A model for bidentate 3MPA coordination at the 3MDO Fe-site has been proposed on the basis of computational docking, whereas steady-state kinetics and EPR spectroscopic measurements suggest a thiolate-only coordination of the substrate. To address this gap in knowledge, we determined the structure of Azobacter vinelandii 3MDO (Av3MDO) in complex with the substrate analog and competitive inhibitor, 3-hydroxypropionic acid (3HPA). The structure together with DFT computational modeling demonstrates that 3HPA and 3MPA associate with iron as chelate complexes with the substrate-carboxylate group forming an additional interaction with Arg168 and the thiol bound at the same position as in CDO. A chloride ligand was bound to iron in the coordination site assigned as the O 2 -binding site. Supporting HYSCORE spectroscopic experiments were performed on the (3MPA/NO)-bound Av3MDO iron nitrosyl (S = 3/2) site. In combination with spectroscopic simulations and optimized DFT models, this work provides an experimentally verified model of the Av3MDO enzyme–substrate complex, effectively resolving a debate in the literature regarding the preferred substrate-binding denticity. These results elegantly explain the observed 3MDO substrate specificity, but leave unanswered questions regarding the mechanism of substrate-gated reactivity with dioxygen.

36 MATERIALS SCIENCE↗

NEAMS Problem Draft (FY2020 Assessment)

This report documents the physical configuration required to model the assessment problem posed by the Reactor Physics focus area within the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. This model is based on the Empire model developed by staff at Los Alamos National Laboratory under the MEITNER (Modeling-Enhanced Innovations Trailblazing Nuclear Energy Reinvigoration) program. The reference report does not provide dimensions for development of computational models; those data are provided elsewhere in published. Thus, this report has been prepared to provide a single reference that can be used to create consistent models by assessment participants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Partially Ionized Plasma Physics and Technological Applications

Partially ionized plasma physics has attracted increased attention recently due to numerous technological applications made possible by the increased sophistication of computer modelling, the depth of the theoretical analysis, and the technological applications to a vast field of manufacturing for computer components. Partially ionized plasma is characterized by a significant presence of neutral particles in contrast to the fully ionized plasma. The theoretical analysis is based upon solutions of the kinetic Boltzmann equation, yielding the non-Maxwellian electron energy distribution function (EEDF), thereby emphasizing the difference with a fully ionized plasma. The impact of the effect on discharges in inert and molecular gases is described in detail, yielding the complex nonlinear phenomena resulting in plasma selforganization. A few examples of such phenomena are given, including the non-monotonic EEDFs in the discharge afterglow in a mixture of argon with the molecular gas NF3; the explosive generation of cold electron populations in capacitive discharges, hysteresis of EEDF in inductively coupled plasmas. Recently, highly advanced computer codes were developed in order to address the outstanding challenges in plasma technology. These developments are briefly described in general terms.

non-Maxwellian electron energy distribution functi↗

Magnetorotational iron core collapse

During its final evolutionary stages, a massive star, as considered in current astrophysical theory, undergoes rapid collapse, thereby triggering a sequence of a catastrophic event which results in a Type II supernova explosion. A remnant neutron star or a black hole is left after the explosion. Stellar collapse occurs, when thermonuclear fusion has consumed the lighter elements present. At this stage, the core consists of iron. Difficulties arise regarding an appropriate model with respect to the core collapse. The present investigation is concerned with the evolution of a Type II supernova core including the effects of rotation and magnetic fields. A simple neutrino model is developed which reproduced the spherically symmetric results of Bowers and Wilson (1982). Several two-dimensional computational models of stellar collapse are studied, taking into account a case in which a 15 solar masses iron core was artificially given rotational and magnetic energy.

Symbalisty, E. M. D.↗