Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first-principles 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 37 records · Page 2

The 13th Technology of Deep Space One - Abstract

On October 24th, 1998, the Deep Space One (DS-1) spacecraft launched aboard a Delta II rocket as the first step towards the bold task of testing and validating 12 new technologies for future missions. This launch also represented yet another thrilling event; namely, the successful test and validation of a 13th heretofore undisclosed technology: model-based code-generation of the spacecraft's system-level fault-protection (FP) software from behavioral state diagrams and structural models.In this paper, we describe the process we used to leverage model-based code generation from state diagrams and structural specifications to better respond to the evolving requirements and scope of DS- I's system-level fault-protection design, development, test and operation. The evolution of the high-level design and the low-level changes in the flight software architecture and interfaces contributed to multiplying the number and frequency of fault-protection software releases thereby creating a multitude of software integration issues. To address the resulting software integration issues, we broadened the scope of code -eneration to other forms of model- based analysis techniques more traditionally associated with first-principle's reasoning about physical models. Additionally, we describe our in-flight launch and initial acquisition experience.

Rouquette, Nicolas↗

A Logical Account of Diagnosis with Multiple Theories

Model-based diagnosis is a powerful, first-principles approach to diagnosis. The primary drawback with model-based diagnosis is that it is based on a system model, and this model might be inappropriate. The inappropriateness of models usually stems from the fundamental tradeoff between completeness and efficiency. Recently, Struss has developed an elegant proposal for diagnosis with multiple models. Struss characterizes models as relations and develops a precise notion of abstraction. He defines relations between models and analyzes the effect of a model switch on the space of possible diagnoses. In this paper we extend Struss's proposal in three ways. First, our account of diagnosis with multiple models is based on representing models as more expressive first-order theories, rather than as relations. A key technical contribution is the use of a general notion of abstraction based on interpretations between theories. Second, Struss conflates component modes with models, requiring him to define models relations such as choices which result in non-relational models. We avoid this problem by differentiating component modes from models. Third, we present a more general account of simplifications that correctly handles situations where the simplification contradicts the base theory.

Pandurang, P.↗

Equivalent crystal theory of alloys

Equivalent Crystal Theory (ECT) is a new, semi-empirical approach to calculating the energetics of a solid with defects. The theory has successfully reproduced surface energies in metals and semiconductors. The theory of binary alloys to date, both with first-principles and semi-empirical models, has not been very successful in predicting the energetics of alloys. This procedure is used to predict the heats of formation, cohesive energy, and lattice parameter of binary alloys of Cu, Ni, Al, Ag, Au, Pd, and Pt as functions of composition. The procedure accurately reproduces the heats of formation versus composition curves for a variety of binary alloys. The results are then compared with other approaches such as the embedded atom and lattice parameters of alloys from pure metal properties more accurately than Vegard's law is presented.

Bozzolo, Guillermo↗

Adhesion at metal interfaces

A basic adhesion process is defined, the theory of the properties influencing metallic adhesion is outlined, and theoretical approaches to the interface problem are presented, with emphasis on first-principle calculations as well as jellium-model calculations. The computation of the energies of adhesion as a function of the interfacial separation is performed; fully three-dimensional calculations are presented, and universality in the shapes of the binding energy curves is considered. An embedded-atom method and equivalent-crystal theory are covered in the framework of issues involved in practical adhesion.

Banerjea, Amitava↗

Kuang's Semi-Classical Formalism for Electron Capture Cross-Sections in Ion-Ion Collisions at few MeV's/nucleon: Application to ENA Modeling

Accurate estimates of electron-capture cross sections at energies relevant to energetic neutral atom (ENA) modeling (approx few MeV per nucleon) and for multi-electron ions must rely on first-principles approaches and/or detailed quantum-mechanical simulation of the collision process. Kuang's semi-classical approach offers a middle-ground, elegant and efficient way to arrive at these estimates. We shall present a sample application and current progress in applying and extending Kuang's formalism to ENA modeling.

Barghouty, A. F.↗

CFD Extraction of Heat Transfer Coefficient in Cryogenic Propellant Tanks

Current reduced-order thermal model for cryogenic propellant tanks is based on correlations built for flat plates collected in the 1950's. The use of these correlations suffers from inaccurate geometry representation; inaccurate gravity orientation; ambiguous length scale; and lack of detailed validation. This study uses first-principles based CFD methodology to compute heat transfer from the tank wall to the cryogenic fluids and extracts and correlates the equivalent heat transfer coefficient to support reduced-order thermal model. The CFD tool was first validated against available experimental data and commonly used correlations for natural convection along a vertically heated wall. Good agreements between the present prediction and experimental data have been found for flows in laminar as well turbulent regimes. The convective heat transfer between the tank wall and cryogenic propellant, and that between the tank wall and ullage gas were then simulated. The results showed that the commonly used heat transfer correlations for either vertical or horizontal plate over-predict heat transfer rate for the cryogenic tank, in some cases by as much as one order of magnitude. A characteristic length scale has been defined that can correlate all heat transfer coefficients for different fill levels into a single curve. This curve can be used for the reduced-order heat transfer model analysis.

Yang, H. Q.↗

Computational Fluid Dynamics Based Extraction of Heat Transfer Coefficient in Cryogenic Propellant Tanks

Current reduced-order thermal model for cryogenic propellant tanks is based on correlations built for flat plates collected in the 1950's. The use of these correlations suffers from: inaccurate geometry representation; inaccurate gravity orientation; ambiguous length scale; and lack of detailed validation. The work presented under this task uses the first-principles based Computational Fluid Dynamics (CFD) technique to compute heat transfer from tank wall to the cryogenic fluids, and extracts and correlates the equivalent heat transfer coefficient to support reduced-order thermal model. The CFD tool was first validated against available experimental data and commonly used correlations for natural convection along a vertically heated wall. Good agreements between the present prediction and experimental data have been found for flows in laminar as well turbulent regimes. The convective heat transfer between tank wall and cryogenic propellant, and that between tank wall and ullage gas were then simulated. The results showed that commonly used heat transfer correlations for either vertical or horizontal plate over predict heat transfer rate for the cryogenic tank, in some cases by as much as one order of magnitude. A characteristic length scale has been defined that can correlate all heat transfer coefficients for different fill levels into a single curve. This curve can be used for the reduced-order heat transfer model analysis.

Yang, H. Q.↗

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

Atomistic Modeling of Co Growth on Cu(111]

The BFS method for alloys is applied to the study of Co growth on Cu(111). The parameterization of the Co-Cu system is obtained from first-principles calculations, and tested against known experimental features for low coverage Co deposition on Cu(100) and Cu(111). Atomistic simulations are performed to investigate the behavior of Co on Cu(111) as a function of coverage.

Khalil, Joseph↗

GOATS - Orbitology Component

The GOATS Orbitology Component software was developed to specifically address the concerns presented by orbit analysis tools that are often written as stand-alone applications. These applications do not easily interface with standard JPL first-principles analysis tools, and have a steep learning curve due to their complicated nature. This toolset is written as a series of MATLAB functions, allowing seamless integration into existing JPL optical systems engineering modeling and analysis modules. The functions are completely open, and allow for advanced users to delve into and modify the underlying physics being modeled. Additionally, this software module fills an analysis gap, allowing for quick, high-level mission analysis trades without the need for detailed and complicated orbit analysis using commercial stand-alone tools. This software consists of a series of MATLAB functions to provide for geometric orbit-related analysis. This includes propagation of orbits to varying levels of generalization. In the simplest case, geosynchronous orbits can be modeled by specifying a subset of three orbit elements. The next case is a circular orbit, which can be specified by a subset of four orbit elements. The most general case is an arbitrary elliptical orbit specified by all six orbit elements. These orbits are all solved geometrically, under the basic problem of an object in circular (or elliptical) orbit around a rotating spheroid. The orbit functions output time series ground tracks, which serve as the basis for more detailed orbit analysis. This software module also includes functions to track the positions of the Sun, Moon, and arbitrary celestial bodies specified by right ascension and declination. Also included are functions to calculate line-of-sight geometries to ground-based targets, angular rotations and decompositions, and other line-of-site calculations. The toolset allows for the rapid execution of orbit trade studies at the level of detail required for the early stage of mission concept development.

Haber, Benjamin M.↗

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗

Inverse isotope effect and the x dependence of the superconducting transition temperature in PdHx and PdDx

First-principles electron-band-structure calculations and measured phonon properties are employed in a calculation of the concentration dependence of the superconducting transition temperatures for PdHx and PdDx, including the inverse isotope effect. Excellent agreement with the experimental results is obtained. Our model also gives a qualitative understanding of the enhanced transition temperatures in the palladium-noble-metal hydrides.

Klein, B. M.↗

Rotor-vortex interaction noise

A theoretical and experimental study is carried out to develop a first-principles analysis for predicting noise generated by helicopter main-rotor shed vortices interacting with the tail rotor. The generalized prediction procedure calls for a knowledge of the incident vortex velocity field, the rotor geometry, and the rotor operating conditions. The analysis encompasses compressibility effects and chordwise and spanwise noncompactness and treats oblique intersections with the blade planform. In assessing the theory, a model rotor experiment which isolates the blade-vortex interaction noise from other rotor noise mechanisms is conducted. An isolated tip vortex, generated by an upstream semispan airfoil, is convected into the model rotor, and acoustic spectra, pressure signatures, and directivity are measured. Ingestion of the vortex by the rotor is experimentally observed to generate harmonic noise and impulsive waveforms.

Schlinker, R. H.↗

Modelling of Surfaces. Part 1: Monatomic Metallic Surfaces Using Equivalent Crystal Theory

We present a detailed description of equivalent crystal theory focusing on its application to the study of surface structure. While the emphasis is in the structure of the algorithm and its computational aspects, we also present a comprehensive discussion on the calculation of surface energies of metallic systems with equivalent crystal theory and other approaches. Our results are compared to experiment and other semiempirical as well as first-principles calculations for a variety of fcc and bcc metals.

Bozzolo, Guillermo↗

Prediction and measurement of low-frequency harmonic noise of a hovering model helicopter rotor

Far-field acoustic data for a model helicopter rotor have been gathered in a large open-jet, acoustically treated wind tunnel with the rotor operating in hover and out of ground-effect. The four-bladed Boeing 360 model rotor with advanced airfoils, planform, and tip shape was run over a range of conditions typical of today's modern helicopter main rotor. Near in-plane acoustic measurements were compared with two independent implementations of classical linear theory. Measured steady thrust and torque were used together with a free-wake analysis (to predict the thrust and drag distributions along the rotor radius) as input to this first-principles theoretical approach. Good agreement between theory and experiment was shown for both amplitude and phase for measurements made in those positions that minimized distortion of the radiated acoustic signature at low-frequencies.

Aggarawal, H. R.↗

Resolution Assessment in Dynamic Image Formation

Remote sensing and astronomical image formation is often complicated by deficiencies in measurement quality, density, or diversity. Penalized likelihood methods can incorporate additional first-principles physical prior knowledge and improve the image reconstructions, but a systematic bias is unavoidable as a consequence. This work derives theory to understand the bias and develops a computational tool to probe its effect on the reconstructed image and bound resolution limits. Though the focus is on image formation, the contributions of this paper apply to any inference problem that can be expressed under the linear state-space signal model.

multidimensional signal processing↗

Noise Prediction for Maneuvering Rotorcraft

This paper presents the initial work toward first-principles noise prediction for maneuvering rotors. Both the aeromechanical and acoustics aspects of the maneuver noise problem are discussed. The comprehensive analysis code, CAMRAD 2. was utilized to predict the time-dependent aircraft position and attitude, along - with the rotor blade airloads and motion. The major focus of this effort was the enhancement of the acoustic code WOPWOP necessary to compute the noise from a maneuvering rotorcraft. Full aircraft motion, including arbitrary transient motion, is modeled together with arbitrary rotor blade motions. Noise from a rotorcraft in turning and descending flight is compared to level flight. A substantial increase in the rotor noise is found both for turning flight and during a transient maneuver. Additional enhancements to take advantage of parallel computers and clusters of workstations, in addition to a new compact-chordwise loading formulation, are also described.

Brentner, Kenneth S.↗

The Importance of the Cathode Plume and Its Interactions with the Ion Beam in Numerical Simulations of Hall Thrusters

Hall2De is a first-principles, 2-D axisymmetric code that solves the equations of motion for ions, electrons, and neutrals on a magnetic-field-aligned grid. The computational domain downstream of the acceleration channel exit plane is large enough to include self-consistently the cathode boundary. In this paper, we present results from numerical simulations of the H6 laboratory thruster with an internally mounted cathode, with the aim of highlighting the importance of properly accounting for the interactions between the ion beam and cathode plume. The anomalous transport of electrons across magnetic field lines in Hall2De is modelled using an anomalous collision frequency, νanom, yielding νanom approximately equal to omega ce (i.e., the electron cyclotron frequency) in the plume. We first show that restricting the anomalous collision frequency to only regions where the current density of ions is large does not alter the plasma discharge in the Hall thruster as long as the interaction between the ion beam and the cathode plume is captured properly in the computational domain. This implies that the boundary conditions must be placed sufficiently far as to not interfere with the electron transport in this region. These simulation results suggest that electron transport across magnetic field lines occurs largely inside the beam and may be driven by the interactions between beam ions and electrons. A second finding that puts in relevance the importance of including the cathode plume in numerical simulations is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in the Hall2De simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume obtained during the operation of the H6 thruster at nominal conditions (300V, 20A) and chamber background pressure of approximately 1.5 x 10(exp -5) Torr. In addition, the numerical results obtained with the latter approach exhibit less sensitivity to background pressure than previous attempts at explaining the features of the plasma properties in the near plume.

Lopez Ortega, Alejandro↗