Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical mathematics”

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 199 records · Page 11

Ballistic transport in planetary ring systems due to particle erosion mechanisms. I - Theory, numerical methods, and illustrative examples

Ballistic transport, defined as the net radial transport of mass and angular momentum due to exchanges of meteoroid hypersonic-impact ejecta by neighboring planetary ring regions on time-scales orders-of-magnitude shorter than the age of the solar system, is presently considered as a problem in mathematical physics. The preliminary results of a numerical scheme for following the combined effects of ballistic transport and viscous diffusion demonstrate that ballistic transport generates structure near sharp edges already present in the ring-mass distribution; the entire ring system ultimately develops an undulatory structure whose length scale is typically of the order of the radial excursion of the impact ejecta.

Durisen, Richard H.↗

Numerical Simulations of Plasma Based Flow Control Applications

A mathematical model was developed to simulate flow control applications using plasma actuators. The effects of the plasma actuators on the external flow are incorporated into Navier Stokes computations as a body force vector. In order to compute this body force vector, the model solves two additional equations: one for the electric field due to the applied AC voltage at the electrodes and the other for the charge density representing the ionized air. The model is calibrated against an experiment having plasma-driven flow in a quiescent environment and is then applied to simulate a low pressure turbine flow with large flow separation. The effects of the plasma actuator on control of flow separation are demonstrated numerically.

Suzen, Y. B.↗

Recovery of thrust magnitude from minitrack data for the SERT-2 spacecraft

An estimation of the thrust magnitude of the Space Electric Rocket Test (SERT) engine is obtained. This estimation is based on changes in the SERT-2 spacecraft orbit as determined from tracking data by using a revised definitive orbit determination system which includes the effect of thrust upon the motion of a spacecraft in general and the SERT-2 satellite in particular. The results of these computations are found to be in good agreement with values of acceleration as determined by an on-board accelerometer. A discussion of the mathematical model used to obtain the numerical results is also presented.

Jones, E. M.↗

Effect of gravity on methane-air combustion

Analytical and numerical techniques dealing with the theoretical description of the influence of zero and reduced gravitational acceleration on diffusion flames, with a view to improving understanding of fires in space vehicles, were developed in support of experimental work performed in this area. This was done in order to confirm qualitative understanding of the process, to determine the quantitative accuracy of numerical predictions, and to establish a mathematical model of the process for subsequent use as a predictive and exploratory tool. The following results were accomplished: (1) derivation of differential equations and boundary conditions describing the system, (2) details of the computations, using a FORTRAN computer program, for calculating the flow and heat and mass transfer in two dimensions (both steady and unsteady). It was shown that the experimental behavior can be reproduced with fair accuracy, provided that the time step is sufficiently short.

Elghobashi, S.↗

Transonic flow calculations

The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.

Jameson, A.↗

Recent advances in the solution of three-dimensional flow over wings with leading edge vortex separation

Recent advances in a panel method for the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation are presented. These advances were achieved as part of an ultimately successful assault on two shortcomings of the method, namely convergence failures in seemingly random cases, and overprediction of lift coefficient for high aspect-ratio wings. Advances include the implementation of improved panel numerics for the purpose of eliminating the highly non-linear effects of ring vortices around doublet panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A variety of cases generated by the computer program implementing the method are presented. These cases are of two types. The first type consists of numerical studies, which verify the underlying mathematical assumptions of the method and moreover show that the results are strongly invariant with respect to such user dependent input as wing panel layout, initial sheet shape, sheet rollup, etc. The second type consists of cases run for the purpose of comparing computed results with experimental data, and these comparisons verify the underlying physical assumptions made by the method.

Johnson, F. T.↗

Hypersonic merged-layer flow on a sphere

With the objective to achieve the desired reduction in velocity by aerodynamic forces, an aeroassisted orbital transfer vehicle (AOTV) is expected to fly in the higher region of the atmosphere for a sustained period of time. This aeroassist maneuver will occur in the transitional regime, if the vehicle is designed for a low ballistic coefficient. In this case, a merged-layer (ML) is formed around the space vehicle. The flow characteristics of the vehicle can be studied by making use of the full Navier-Stokes (NS) equations, taking into account slip and temperature jump conditions. The present investigation has the objective to obtain a two-term series solution of the full Navier-Stokes equations with surface slip and temperature jump conditions for the ML flow on a sphere. Attention is given to the mathematical formulation of the problem, the numerical method of integration, and the results.

Jain, A. C.↗

Analysis of boundary conditions for SSME subsonic internal viscous flow analysis

A study was completed of mathematically proper boundary conditions for unique numerical solution of internal, viscous, subsonic flows in the space shuttle main engine. The study has concentrated on well posed considerations, with emphasis on computational efficiency and numerically stable boundary condition statements. The method of implementing the established boundary conditions is applicable to a wide variety of finite difference and finite element codes, as demonstrated.

Baker, A. J.↗

Artificial Dissipation In Computations Of Hypersonic Flows

Four artificial-dissipation mathematical models introduced to suppress spurious numerical oscillations in finite-difference computations of hypersonic external flows. Models tested for their effects in capturing details of shocks in hypersonic flows about bodies of various shapes.

Yoon, S.↗

Mathematical Models Of Turbulence In Hypersonic Flow

Report discusses mathematical models of turbulence used in numerical simulations of complicated viscous, hypersonic flows. Includes survey of essential features of models and their statuses in applications.

Marvin, J. G.↗

Identification of linear systems by an asymptotically stable observer

A formulation is presented for the identification of a linear multivariable system from single or multiple sets of input-output data. The system input-output relationship is expressed in terms of an observer, which is made asymptotically stable by an embedded eigenvalue assignment procedure. The prescribed eigenvalues for the observer may be real, complex, mixed real and complex, or zero. In this formulation, the Markov parameters of the observer are identified from input-output data. The Markov parameters of the actual system are then recovered from those of the observer and used to obtain a state space model of the system by standard realization techniques. The basic mathematical formulation is derived, and extensive numerical examples using simulated noise-free data are presented to illustrate the proposed method.

Phan, Minh Q.↗

A mathematical model of physiological processes and its application to the study of aging

The behavior of a physiological system which, after displacement, returns by homeostatic mechanisms to its original condition can be described by a simple differential equation in which the "recovery time" is a parameter. Two such systems, which influence one another, can be linked mathematically by the use of "coupling" or "feedback" coefficients. These concepts are the basis for many mathematical models of physiological behavior, and we describe the general nature of such models. Next, we introduce the concept of a "fatal limit" for the displacement of a physiological system, and show how measures of such limits can be included in mathematical models. We show how the numerical values of such limits depend on the values of other system parameters, i.e., recovery times and coupling coefficients, and suggest ways of measuring all these parameters experimentally, for example by monitoring changes induced by X-irradiation. Next, we discuss age-related changes in these parameters, and show how the parameters of mortality statistics, such as the famous Gompertz parameters, can be derived from experimentally measurable changes. Concepts of onset-of-aging, critical or fatal limits, equilibrium value (homeostasis), recovery times and coupling constants are involved. Illustrations are given using published data from mouse and rat populations. We believe that this method of deriving survival patterns from model that is experimentally testable is unique.

Aging/physiology↗

Forward flight trim calculation and frequency response validation of a high-order helicopter simulation model

This paper describes a new trim procedure, that includes the calculation of the steady-state response of the rotor blades, and that is applicable to straight flight and steady coordinated turns. This paper also describes the results of a validation study for a high order linearized model of helicopter flight dynamics, that includes rotor, inflow, and actuator dynamics. The model is obtained by numerical perturbations of a nonlinear, blade element type mathematical model. Predicted responses are compared with flight test data for two values of flight speed. The comparison is carried out in the frequency domain. Numerical simulations show that the trim algorithm is very accurate, and preserves the periodicity of the aircraft states. The results also indicate that the predictions of the linearized model are in good agreement with flight test data, especially at medium and high frequencies.

Kim, Frederick D.↗

Modeling and control of flexible space platforms with articulated payloads

The first steps in developing a methodology for spacecraft control-structure interaction (CSI) optimization are identification and classification of anticipated missions, and the development of tractable mathematical models in each mission class. A mathematical model of a generic large flexible space platform (LFSP) with multiple independently pointed rigid payloads is considered. The objective is not to develop a general purpose numerical simulation, but rather to develop an analytically tractable mathematical model of such composite systems. The equations of motion for a single payload case are derived, and are linearized about zero steady-state. The resulting model is then extended to include multiple rigid payloads, yielding the desired analytical form. The mathematical models developed clearly show the internal inertial/elastic couplings, and are therefore suitable for analytical and numerical studies. A simple decentralized control law is proposed for fine pointing the payloads and LFSP attitude control, and simulation results are presented for an example problem. The decentralized controller is shown to be adequate for the example problem chosen, but does not, in general, guarantee stability. A centralized dissipative controller is then proposed, requiring a symmetric form of the composite system equations. Such a controller guarantees robust closed loop stability despite unmodeled elastic dynamics and parameter uncertainties.

Graves, Philip C.↗

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We describe a variation of the Optimal Estimation (OE) method for greenhouse gas remote sensing retrievals using a singular value decomposition (SVD) and an uninformative prior. The SVD method is capable of discerning vertical information in column integrated absorption measurements. While traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We discuss the method, illustrate its use on an integrated path differential absorption CO2 lidar measurement model, and compare it to traditional optimal estimation using numerical simulations. In the absence of forward model error, the mathematics behind the SVD method guarantee it to be bias-free, which is confirmed by the numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. While the SVD approach can be used for most trace gas retrievals, it is particularly useful for situations where prior knowledge of the trace gas profile is poor. The SVD analysis is also useful in illustrating how vertical information is treated by the traditional OE approach.

Ramanathan, Anand K.↗

A comparison of two recent approaches to machine parameterization for mathematical software

Two proposals for expressing the dependence of numerical software on the environment are examined. The approaches are briefly summarized and the differences and similarities are discussed. One proposal characterizes the dependence of the software on the machine architecture essentially in terms of the representation of floating point entities in the machine. The second characterizes the dependence of the software on the environment in terms of models both of the floating point numbers and of the behavior of the arithmetic unit of the machine.

Smith, B. T.↗

Wing shape optimization for maximum cross-country speed, with mathematical programming

A computer program was developed to calculate numerically the speed and circling polars of an aircraft when the lift and drag characteristics of the wing airfoils are known. The planform of the wing is described by variables which are optimized so that the cross-country speed of the glider is maximum for the particular type of thermal model. Two thermal models are compared and shown that with a greater wing area than now normally used the performance can be increased.

Helwig, G.↗