Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Motion models in attitude estimation

Attitude estimator use observations from different times to reduce the effects of noise. If the vehicle is rotating, the attitude at one time needs to be propagated to that at another time. If the vehicle measures its angular velocity, attitude propagating entails integrating a rotational kinematics equation only. If a measured angular velocity is not available, torques can be computed and an additional rotational dynamics equation integrated to give the angular velocity. Initial conditions for either of these integrations come from the estimation process. Sometimes additional quantities, such as gyro and torque parameters, are also solved for. Although the partial derivatives of attitude with respect to initial attitude and gyro parameters are well known, the corresponding partial derivatives with respect to initial angular velocity and torque parameters are less familiar. They can be derived and computed numerically in a way that is analogous to that used for the initial attitude and gyro parameters. Previous papers have demonstrated the feasibility of using dynamics models for attitude estimation but have not provided details of how each angular velocity and torque parameters can be estimated. This tutorial paper provides some of that detail, notably how to compute the state transition matrix when closed form expressions are not available. It also attempts to put dynamics estimation in perspective by showing the progression from constant to gyro-propagated to dynamics-propagated attitude motion models. Readers not already familiar with attitude estimation will find this paper an introduction to the subject, and attitude specialists may appreciate the collection of heretofore scattered results brought together in a single place.

Chu, D.↗
NASA NTRS1989

Motivating the sure bounds

Motivation is provided for a theorem that provides upper and lower bounds for the reliability of reconfigurable digital control systems. The reliability goals for these systems are too high to be established by natural life testing, which means the probability of system failure must be computed from mathematical models that capture the essential elements of fault occurence and system fault recovery. The upper and lower bound theorem shows that system recovery can be adequately described by its first two moments, provided component failure rate is low and system recovery is fast. This result greatly simplifies both the fault injection experiments that study system recovery and the numerical computations that estimate the probability of system failure from a mathematical model.

White, Allan L.↗
NASA NTRS1980

Comparison of experimental and computational shock structure in a transonic compressor rotor

Measurement of passage shock strength in a transonic compressor rotor using a gas fluorescent technique revealed an unexpected variation in shock strength in the radial direction. An axisymmetric idealization would normally predict that the passage shock strength would gradually weaken when moving radially inward until disappearing at the sonic radius. However, the measurements indicated a sharp peak in strength at the nominal sonic radius. Blade boundary layer separation originating at this point accounts for about one half of the total rotor losses. A numerical computation of the three-dimensional inviscid flow, using time-marching techniques, has accurately predicted in general the radial and tangential variations in passage shock strength and in particular the sharp pressure peak at the nominal sonic radius. The overall shock strength was somewhat overpredicted, but this overprediction may be the result of boundary layer separation in the experiment. This paper presents comparisons between the optical density measurements and computational results and in addition a short analytical discussion which demonstrates that the sharp shock strength rise may occur in many transonic compressor rotors.

Haymann-Haber, G.↗
NASA NTRS2019

Recommended Methods for Setting Mission Conjunction Analysis Hard Body Radii

For real-time conjunction assessment (CA) operations, computation of the Probability of Collision (P(sub c)) typically depends on the state vector, its covariance, and the combined hard body radius (HBR) of both the primary and secondary space-craft. However, most algorithmic approaches that compute the P(sub c) use generic conservatively valued HBRs that may tend to go beyond the physical limitations of both spacecraft, enough to drastically change the results of a conjunction assessment mitigation decision. On the other hand, if the attitude of the spacecraft is known and available, then a refined HBR can be obtained that could result in an improved and accurate numerically-computed P(sub c) value. The goal of this analysis is to demonstrate the various calculated P(sub c) values obtained based on a number of different HBR calculation techniques, oriented in the encounter or conjunction plane at the time of closest approach (TCA). Since in most conjunctions the secondary object is a debris object and thus orders of magnitude smaller than the primary, the greatest operational benefit is wrought by developing a better size estimate and representation for the primary object. We present an analysis that includes the attitude information of the primary object in the HBR calculation and assesses the resulting P(sub c) values for conjunction assessment decision making.

Mashiku, Alinda K.↗
NASA NTRS2019

Recommended Methods for Setting Mission Conjunction Analysis Hard Body Radii

For real-time conjunction assessment (CA) operations, computation of the Probability of Collision (P(sub c)) typically depends on the state vector, its covariance, and the combined hard body radius (HBR) of both the primary and secondary space-craft. However, most algorithmic approaches that compute the P(sub c) use generic conservatively valued HBRs that may tend to go beyond the physical limitations of both spacecraft, enough to drastically change the results of a conjunction assessment mitigation decision. On the other hand, if the attitude of the spacecraft is known and available, then a refined HBR can be obtained that could result in an improved and accurate numerically-computed P(sub c) value. The goal of this analysis is to demonstrate the various calculated P(sub c) values obtained based on a number of different HBR calculation techniques, oriented in the encounter or conjunction plane at the time of closest approach (TCA). Since in most conjunctions the secondary object is a debris object and thus orders of magnitude smaller than the primary, the greatest operational benefit is wrought by developing a better size estimate and representation for the primary object. We present an analysis that includes the attitude information of the primary object in the HBR calculation and assesses the resulting P(sub c) values for conjunction assessment decision making.

Mashiku, Alinda K.↗
NASA NTRS1976

Computation of turbulent flows

The essential ingredients and the effectiveness of several levels of turbulent-flow partial differential equation models are considered. Zero-equation models use only the partial differential equation (pde) for the mean velocity field and do not employ turbulence pde's. One-equation models make use of an additional pde relating to the turbulence velocity scale. Attention is also given to two-equation models, stress-equation models, and large eddy simulations. Large-eddy simulations are concerned with a three-dimensional time-dependent numerical computation of the large-scale turbulence.

Reynolds, W. C.↗
NASA NTRS1990

Resonant decay of global MHD modes at 'thick' interfaces

A global surface mode can exist on a 'surface' which is not a true discontinuity, but without sufficient dissipation it is not a normal mode and it decays in time via resonance absorption. The decay rate can be calculated analytically when the 'surface' is thin. The goal of this paper is to numerically estimate the decay rate when the surface is not thin. A cold plasma in a uniform magnetic field is considered in which density varies linearly across the 'surface'. In this linearized calculation, the global surface mode is driven in steady state by an antenna located in one of the uniform regions external to the surface. The frequency is a free parameter, and resonance curves are computed numerically without approximations. The widths of the resonance curves are used to estimate the free decay times of undriven surface modes, via the uncertainty principle. When applied to the active solar corona, the decay rates are large enough to account for the coronal heating, but it should be kept in mind that the role of nonlinearity in resonance absorption is still undetermined.

Hollweg, Joseph V.↗
NASA NTRS1990

Dynamic array processing for computationally intensive expert systems in CLIPS

This paper puts forth an architecture for implementing a loop for advanced data structure of arrays in CLIPS. An attempt is made to use multi-field variables in such an architecture to process a set of data during the decision making cycle. Also, current limitations on the expert system shells are discussed in brief in this paper. The resulting architecture is designed to circumvent the current limitations set by the expert system shell and also by the operating environment. Such advanced data structures are needed for tightly coupling symbolic and numeric computation modules.

Athavale, N. N.↗
NASA NTRS1989

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗
NASA NTRS2004

Theory, Modeling, Software and Hardware Development for Analytical and Computational Materials Science

The focus of this Cooperative Agreement between the Computational Materials Laboratory (CML) of the Processing Science and Technology Branch of the NASA Glenn Research Center (GRC) and the Department of Theoretical and Applied Mathematics at The University of Akron was in the areas of system development of the CML workstation environment, modeling of microgravity and earth-based material processing systems, and joint activities in laboratory projects. These efforts complement each other as the majority of the modeling work involves numerical computations to support laboratory investigations. Coordination and interaction between the modelers, system analysts, and laboratory personnel are essential toward providing the most effective simulations and communication of the simulation results. Toward these means, The University of Akron personnel involved in the agreement worked at the Applied Mathematics Research Laboratory (AMRL) in the Department of Theoretical and Applied Mathematics while maintaining a close relationship with the personnel of the Computational Materials Laboratory at GRC. Network communication between both sites has been established. A summary of the projects we undertook during the time period 9/1/03 - 6/30/04 is included.

Young, Gerald W.↗
NASA NTRS1980

Radiation boundary conditions for wave-like equations

In the numerical computation of hyperbolic equations it is not practical to use infinite domains; instead, the domain is truncated with an artificial boundary. In the present study, a sequence of radiating boundary conditions is constructed for wave-like equations. It is proved that as the artificial boundary is moved to infinity the solution approaches the solution of the infinite domain as O(r exp -m-1/2) for the m-th boundary condition. Numerical experiments with problems in jet acoustics verify the practical nature of the boundary conditions.

Bayliss, A.↗
NASA NTRS2004

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗
NASA NTRS1980

Comparison of undulation difference accuracies using gravity anomalies and gravity disturbances

Errors in the outer zone contribution to oceanic undulation differences computed from a finite set of potential coefficients based on satellite measurements of gravity anomalies and gravity disturbances are analyzed. Equations are derived for the truncation errors resulting from the lack of high-degree coefficients and the commission errors arising from errors in the available lower-degree coefficients, and it is assumed that the inner zone (spherical cap) is sufficiently covered by surface gravity measurements in conjunction with altimetry or by gravity anomaly data. Numerical computations of error for various observational conditions reveal undulation difference errors ranging from 13 to 15 cm and from 6 to 36 cm in the cases of gravity anomaly and gravity disturbance data, respectively for a cap radius of 10 deg and mean anomalies accurate to 10 mgal, with a reduction of errors in both cases to less than 10 cm as mean anomaly accuracy is increased to 1 mgal. In the absence of a spherical cap, both cases yield error estimates of 68 cm for an accuracy of 1 mgal and between 93 and 160 cm for the lesser accuracy, which can be reduced to about 110 cm by the introduction of a perfect 30-deg reference field.

Jekeli, C.↗
NASA NTRS1988

Symbolic derivation of material property matrices in finite element analysis

The principles and operation of MMAX, a symbolic-computation program which automates the process of generating property matrices for structural materials, are briefly described and illustrated with sample analyses of a rubberlike material and an elastoplastic material. MMAX is written in LISP under the symbolic finite-element generator FINGER and the general symbolic manipulator MACSYMA; it first derives the formulas required by mathematical manipulation, and then translates the formulas into FORTRAN code, adapted to the particular type of machine to be used for the numerical calculations. This approach is shown to combine efficiently the advantages of symbolic and numerical computation for engineering applications.

Tan, H. Q.↗
NASA NTRS2019

Efficient On-Orbit Singularity-Free Geopotential Estimation

For vehicles in orbital cruise, either in orbit around a body or in a transfer orbit, gravity is the primary external force; therefore, the complexity of the geopotential model heavily influences the accuracy of the navigation state. As this type of spacecraft requires a highly accurate geopotential model and performs the computation numerous times a second, this calculation needs to be efficient. This calculation needs to include the first and second derivatives of the geopotential because navigation systems can require both the gravitational acceleration and gradient for inertial integration and state filtering. The most efficient method for calculating the geopotential, the forward column method, is discussed in detail, as well as a method to avoid the singularities that exist when using this method. This is shown to decrease the computation time of the geopotential compared to other popular methods. In addition, methods for first and second order propagation are discussed, which decrease the rate at which the full geopotential model needs to be calculated while maintaining the accuracy required of navigation systems. These estimation methods are shown to decrease the necessary computation of the gravity model by multiple orders of magnitude. This method decreases the computation time of the Exploration Upper Stage navigation system geopotential model by potentially an order of magnitude compared with the previous model without affecting the navigation error. A properly implemented geopotential model can have the accuracy of an 8 by 8 model while approaching the computational requirement as using only the J234 coefficients.

Amert, Joel↗
NASA NTRS1974

Application of a numerical differencing analyzer computer program to a Modular Integrated Utility System

This paper describes the application of a numerical differencing analyzer computer program to the thermal analyzation of a MIUS model. The MIUS model which was evaluated is one which would be required to support a 648-unit Garden Apartment Complex. This computer program was capable of predicting the thermal performance of this MIUS from the impressed electrical, heating, and cooling loads.

Brandli, A. E.↗
NASA NTRS1999

Magnetic Photon Splitting: The S-Matrix Formulation in the Landau Representation

Calculations of reaction rates for the third-order QED process of photon splitting gamma yields gamma.gamma in strong magnetic fields traditionally have employed either the effective Lagrangian method or variants of Schwinger's proper-time technique. Recently, Mentzel, Berg and Wunner [1] presented an alternative derivation via an S-matrix formulation in the Landau representation. Advantages of such a formulation include the ability to compute rates near pair resonances above pair threshold. This paper presents new developments of the Landau representation formalism as applied to photon splitting, providing significant, advances beyond the work of [1] by summing over the spin quantum numbers of the electron propagators, and analytically integrating over the component of momentum of the intermediate states that is parallel to field. The ensuing tractable expressions for the scattering amplitudes are satisfyingly compact, and of an appearance familiar to S-matrix theory applications. Such developments can facilitate numerical computations of splitting considerably both below and above pair threshold. Specializations to two regimes of interest are obtained, namely the limit of highly supercritical fields and the domain where photon energies are far inferior to that for the threshold of single-photon pair creation. In particular, for the first time the low-frequency amplitudes are simply expressed in terms of the Gamma function, its integral and its derivatives. In addition, the equivalence of the asymptotic forms in these two domains to extant results from effective Lagrangian/proper- time formulations is demonstrated.

Baring, Matthew G.↗
NASA NTRS1999

Parallel Tetrahedral Mesh Adaptation with Dynamic Load Balancing

The ability to dynamically adapt an unstructured grid is a powerful tool for efficiently solving computational problems with evolving physical features. In this paper, we report on our experience parallelizing an edge-based adaptation scheme, called 3D_TAG. using message passing. Results show excellent speedup when a realistic helicopter rotor mesh is randomly refined. However. performance deteriorates when the mesh is refined using a solution-based error indicator since mesh adaptation for practical problems occurs in a localized region., creating a severe load imbalance. To address this problem, we have developed PLUM, a global dynamic load balancing framework for adaptive numerical computations. Even though PLUM primarily balances processor workloads for the solution phase, it reduces the load imbalance problem within mesh adaptation by repartitioning the mesh after targeting edges for refinement but before the actual subdivision. This dramatically improves the performance of parallel 3D_TAG since refinement occurs in a more load balanced fashion. We also present optimal and heuristic algorithms that, when applied to the default mapping of a parallel repartitioner, significantly reduce the data redistribution overhead. Finally, portability is examined by comparing performance on three state-of-the-art parallel machines.

Oliker, Leonid↗