Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical analysis”

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 991 records · Page 55

Theoretical and experimental separation dynamics in capillary zone electrophoresis

The mathematical model of Bier et al. (1983) is used in a computer aided analysis of the conditions in capillary zone electrophoresis (ZE) under which sample zones migrate noninteractively with the carrier electrolyte. The monitoring of sample zones with a capillary analyzer that features both on-line conductivity and UV detection at the end of the separation trough is discussed. Data from a ZE analysis of a 5-component mixture are presented, and it is noted that all five components can be monitored via their conductivity change if enough sample is present. It is suggested from the results that the concentration ratio of background buffer to sample should be a minimum of 100:1 to effectively apply the plate concept to ZE.

Thormann, Wolfgang↗

Numerical, analytical, experimental study of fluid dynamic forces in seals

NASA/Lewis Research Center is sponsoring a program for providing computer codes for analyzing and designing turbomachinery seals for future aerospace and engine systems. The program is made up of three principal components: (1) the development of advanced three dimensional (3-D) computational fluid dynamics codes, (2) the production of simpler two dimensional (2-D) industrial codes, and (3) the development of a knowledge based system (KBS) that contains an expert system to assist in seal selection and design. The first task has been to concentrate on cylindrical geometries with straight, tapered, and stepped bores. Improvements have been made by adoption of a colocated grid formulation, incorporation of higher order, time accurate schemes for transient analysis and high order discretization schemes for spatial derivatives. This report describes the mathematical formulations and presents a variety of 2-D results, including labyrinth and brush seal flows. Extensions of 3-D are presently in progress.

Shapiro, William↗

Integrated design of structures, controls, and materials

In this talk we shall discuss algorithms and CAD tools for the design and analysis of structures for high performance applications using advanced composite materials. An extensive mathematical theory for optimal structural (e.g., shape) design was developed over the past thirty years. Aspects of this theory have been used in the design of components for hypersonic vehicles and thermal diffusion systems based on homogeneous materials. Enhancement of the design methods to include optimization of the microstructure of the component is a significant innovation which can lead to major enhancements in component performance. Our work is focused on the adaptation of existing theories of optimal structural design (e.g., optimal shape design) to treat the design of structures using advanced composite materials (e.g., fiber reinforced, resin matrix materials). In this talk we shall discuss models and algorithms for the design of simple structures from composite materials, focussing on a problem in thermal management. We shall also discuss methods for the integration of active structural controls into the design process.

Blankenship, G. L.↗

A High Order Finite Difference Scheme with Sharp Shock Resolution for the Euler Equations

We derive a high-order finite difference scheme for the Euler equations that satisfies a semi-discrete energy estimate, and present an efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semi-discrete energy estimate is based on a symmetrization of the Euler equations that preserves the homogeneity of the flux vector, a canonical splitting of the flux derivative vector, and the use of difference operators that satisfy a discrete analogue to the integration by parts procedure used in the continuous energy estimate. Around discontinuities or sharp gradients, refined grids are created on which the discrete equations are solved after adding a newly constructed artificial viscosity. The positioning of the sub-grids and computation of the viscosity are aided by a detection algorithm which is based on a multi-scale wavelet analysis of the pressure grid function. The wavelet theory provides easy to implement mathematical criteria to detect discontinuities, sharp gradients and spurious oscillations quickly and efficiently.

Gerritsen, Margot↗

Current Results and Proposed Activities in Microgravity Fluid Dynamics

The Institute for Problems in Mechanics' Laboratory work in mathematical and physical modelling of fluid mechanics develops models, methods, and software for analysis of fluid flow, instability analysis, direct numerical modelling and semi-empirical models of turbulence, as well as experimental research and verification of these models and their applications in technological fluid dynamics, microgravity fluid mechanics, geophysics, and a number of engineering problems. This paper presents an overview of the results in microgravity fluid dynamics research during the last two years. Nonlinear problems of weakly compressible and compressible fluid flows are discussed.

Polezhaev, V. I.↗

Control Design Strategies to Enhance Long-Term Aircraft Structural Integrity

Over the operational lifetime of both military and civil aircraft, structural components are exposed to hundreds of thousands of low-stress repetitive load cycles and less frequent but higher-stress transient loads originating from maneuvering flight and atmospheric gusts. Micro-material imperfections in the structure, such as cracks and debonded laminates, expand and grow in this environment, reducing the structural integrity and shortening the life of the airframe. Extreme costs associated with refurbishment of critical load-bearing structural components in a large fleet, or altogether reinventoring the fleet with newer models, indicate alternative solutions for life extension of the airframe structure are highly desirable. Increased levels of operational safety and reliability are also important factors influencing the desirability of such solutions. One area having significant potential for impacting crack growth/fatigue damage reduction and structural life extension is flight control. To modify the airframe response dynamics arising from command inputs and gust disturbances, feedback loops are routinely applied to vehicles. A dexterous flight control system architecture senses key vehicle motions and generates critical forces/moments at multiple points distributed throughout the airframe to elicit the desired motion characteristics. In principle, these same control loops can be utilized to influence the level of exposure to harmful loads during flight on structural components. Project objectives are to investigate and/or assess the leverage control has on reducing fatigue damage and enhancing long-term structural integrity, without degrading attitude control and trajectory guidance performance levels. In particular, efforts have focused on the effects inner loop control parameters and architectures have on fatigue damage rate. To complete this research, an actively controlled flexible aircraft model and a new state space modeling procedure for crack growth have been utilized. Analysis of the analytical state space model for crack growth revealed the critical mathematical factors, and hence the physical mechanism they represent, that influenced high rates of airframe crack growth. The crack model was then exercised with simple load inputs to uncover and expose key crack growth behavior. To characterize crack growth behavior, both "short-term" laboratory specimen test type inputs and "long-term" operational flight type inputs were considered. Harmonic loading with a single overload revealed typical exponential crack growth behavior until the overload application, after which time the crack growth was retarded for a period of time depending on the overload strength. An optimum overload strength was identified which leads to maximum retardation of crack growth. Harmonic loading with a repeated overload of varying strength and frequency again revealed an optimum overload trait for maximizing growth retardation. The optimum overload strength ratio lies near the range of 2 to 3 with dependency on frequency. Experimental data was found to correlate well with the analytical predictions.

Newman, Brett A.↗

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗

Assessment of a One-Dimensional Finite Element Charring Ablation Material Response Model for Phenolic-Impregnated Carbon Ablator

In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.

Thermal Protection System↗

Decomposition-aggregation stability analysis

This report presents the development and description of the decomposition aggregation approach to stability investigations of high dimension mathematical models of dynamic systems. The high dimension vector differential equation describing a large dynamic system is decomposed into a number of lower dimension vector differential equations which represent interconnected subsystems. Then a method is described by which the stability properties of each subsystem are aggregated into a single vector Liapunov function, representing the aggregate system model, consisting of subsystem Liapunov functions as components. A linear vector differential inequality is then formed in terms of the vector Liapunov function. The matrix of the model, which reflects the stability properties of the subsystems and the nature of their interconnections, is analyzed to conclude over-all system stability characteristics. The technique is applied in detail to investigate the stability characteristics of a dynamic model of a hypothetical spinning Skylab.

Siljak, D. D.↗

Computer graphics in aerodynamic analysis

The use of computer graphics and its application to aerodynamic analyses on a routine basis is outlined. The mathematical modelling of the aircraft geometries and the shading technique implemented are discussed. Examples of computer graphics used to display aerodynamic flow field data and aircraft geometries are shown. A future need in computer graphics for aerodynamic analyses is addressed.

Cozzolongo, J. V.↗

Subspace Iteration Method for Complex Eigenvalue Problems with Nonsymmetric Matrices in Aeroelastic System

Modern airplane design is a multidisciplinary task which combines several disciplines such as structures, aerodynamics, flight controls, and sometimes heat transfer. Historically, analytical and experimental investigations concerning the interaction of the elastic airframe with aerodynamic and in retia loads have been conducted during the design phase to determine the existence of aeroelastic instabilities, so called flutter .With the advent and increased usage of flight control systems, there is also a likelihood of instabilities caused by the interaction of the flight control system and the aeroelastic response of the airplane, known as aeroservoelastic instabilities. An in -house code MPASES (Ref. 1), modified from PASES (Ref. 2), is a general purpose digital computer program for the analysis of the closed-loop stability problem. This program used subroutines given in the International Mathematical and Statistical Library (IMSL) (Ref. 3) to compute all of the real and/or complex conjugate pairs of eigenvalues of the Hessenberg matrix. For high fidelity configuration, these aeroelastic system matrices are large and compute all eigenvalues will be time consuming. A subspace iteration method (Ref. 4) for complex eigenvalues problems with nonsymmetric matrices has been formulated and incorporated into the modified program for aeroservoelastic stability (MPASES code). Subspace iteration method only solve for the lowest p eigenvalues and corresponding eigenvectors for aeroelastic and aeroservoelastic analysis. In general, the selection of p is ranging from 10 for wing flutter analysis to 50 for an entire aircraft flutter analysis. The application of this newly incorporated code is an experiment known as the Aerostructures Test Wing (ATW) which was designed by the National Aeronautic and Space Administration (NASA) Dryden Flight Research Center, Edwards, California to research aeroelastic instabilities. Specifically, this experiment was used to study an instability known as flutter. ATW was a small-scale airplane wing comprised of an airfoil and wing tip boom. This wing was formulated based on a NACA-65A004 airfoil shape with a 3.28 aspect ratio. The wing had a span of 18 inch with root chord length of 13.2 inch and tip chord length of 8.7 inch. The total area of this wing was 197 square inch. The wing tip boom was a 1 inch diameter hollow tube of length 21.5 inch. The total weight of the wing was 2.66 lbs.

Pak, Chan-gi↗

Discrete Deterministic and Stochastic Petri Nets

Petri nets augmented with timing specifications gained a wide acceptance in the area of performance and reliability evaluation of complex systems exhibiting concurrency, synchronization, and conflicts. The state space of time-extended Petri nets is mapped onto its basic underlying stochastic process, which can be shown to be Markovian under the assumption of exponentially distributed firing times. The integration of exponentially and non-exponentially distributed timing is still one of the major problems for the analysis and was first attacked for continuous time Petri nets at the cost of structural or analytical restrictions. We propose a discrete deterministic and stochastic Petri net (DDSPN) formalism with no imposed structural or analytical restrictions where transitions can fire either in zero time or according to arbitrary firing times that can be represented as the time to absorption in a finite absorbing discrete time Markov chain (DTMC). Exponentially distributed firing times are then approximated arbitrarily well by geometric distributions. Deterministic firing times are a special case of the geometric distribution. The underlying stochastic process of a DDSPN is then also a DTMC, from which the transient and stationary solution can be obtained by standard techniques. A comprehensive algorithm and some state space reduction techniques for the analysis of DDSPNs are presented comprising the automatic detection of conflicts and confusions, which removes a major obstacle for the analysis of discrete time models.

MATHEMATICAL MODELS↗

Investigation of crew motion disturbances on Skylab-Experiment T-013

Astronaut crew motions can produce some of the largest disturbances acting on a manned spacecraft which can affect vehicle attitude and pointing. Skylab Experiment T-013 was developed to investigate the magnitude and effects of some of these disturbances on the Skylab spacecraft. The methods and techniques used to carry out this experiment are discussed, and preliminary results of data analysis presented. Initial findings indicate that forces on the order of 300 N were exerted during vigorous soaring activities, and that certain experiment activities produced spacecraft angular rate excursions 0.03 to 0.07 deg/sec. Results of Experiment T-013 will be incorporated into mathematical models of crew-motion disturbances, and are expected to be of significant aid in the sizing, design, and analysis of stabilization and control systems for future manned spacecraft.

Conway, B. A.↗

Flight flutter testing technology at Grumman

Analysis techniques used in the automated telemetry station (ATS) for on line data reduction are encompassed in a broad range of software programs. Concepts that form the basis for the algorithms used are mathematically described. The control the user has in interfacing with various on line programs is discussed. The various programs are applied to an analysis of flight data which includes unimodal and bimodal response signals excited via a swept frequency shaker and/or random aerodynamic forces. A nonlinear response error modeling analysis approach is described. Preliminary results in the analysis of a hard spring nonlinear resonant system are also included.

Perangelo, H. J.↗

Predicting Tail Buffet Loads of a Fighter Airplane

Buffet loads on aft aerodynamic surfaces pose a recurring problem on most twin-tailed fighter airplanes: During maneuvers at high angles of attack, vortices emanating from various surfaces on the forward parts of such an airplane (engine inlets, wings, or other fuselage appendages) often burst, immersing the tails in their wakes. Although these vortices increase lift, the frequency contents of the burst vortices become so low as to cause the aft surfaces to vibrate destructively. Now, there exists a new analysis capability for predicting buffet loads during the earliest design phase of a fighter-aircraft program. In effect, buffet pressures are applied to mathematical models in the framework of a finite-element code, complete with aeroelastic properties and working knowledge of the spatiality of the buffet pressures for all flight conditions. The results of analysis performed by use of this capability illustrate those vibratory modes of a tail fin that are most likely to be affected by buffet loads. Hence, the results help in identifying the flight conditions during which to expect problems. Using this capability, an aircraft designer can make adjustments to the airframe and possibly the aerodynamics, leading to a more robust design.

Moses, Robert W.↗