Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complexity theory”

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 883 records · Page 49

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

Geodynamics of Cenozoic deformation in central Asia

This paper presents a study of the tectonic stresses in central Asia based on an interpretation of satellite gravity data for mantle convection and supplemented with published fault plane solutions of earthquakes. Northwest-southeast to north-south compressional stresses exist in the Tien Shan region where reverse faulting dominates. The maximum compressive stress is oriented approximately northeast-southwest in the regions of Altai and southern Mongolia. Farther north, compressive stress gives way to tensional stress which causes normal faulting in the Baikal rift system. It is also shown that all of the tectonic stresses in the Tibetan plateau and Himalayan frontal thrust are related to the convection-generated stress patterns inferred from satellite gravity data. These results suggest that the complex crustal deformation in central Asia can be convincingly described by the deformation of the lithosphere on top of the up- and down-welling asthenospheric material beneath it. This observational fact may not only upset the simple view of the fluid crustal model of the Tibetan plateau, but also provide some useful constraints for the future development of deformation theory of continental crust.

Liu, H.-S.↗

Investigation of rotor wake interactions with a body in low speed forward flight

Wind tunnel tests are reported in which the unsteady aerodynamic interactions between rotor wake and a body are examined for the case of low-speed flight. A wide-field shadowgraph is employed to visualize the rotor-tip vortices during the interactions, and the pressure signatures are measured and correlated with wake trajectories. Consideration is given to unsteady potential-flow theory to estimate the pressure loads from measured wake trajectories. The pressure signatures are very different upstream of, downstream of, and at the wake-impingement point, at which point large transient loads are transferred to the body. Rotor wake is characterized, and the wake filaments impinge on the advancing side of the body first. The pressure signatures present a complex and varying range of patterns, but it is determined that vortex surface interactions can be modeled using potential-flow analysis.

Bi, Nai-Pei↗

Complexities of nadir-looking radiometric temperature measurements of plant canopies

Effective radiant temperatures (ERTs) of five wheat canopies in different stages of development were measured during morning and noon periods. The observed variability in nadir sensor response was quantitatively described as a function of canopy structure and the vertical temperature profile of canopy components. In many cases, the nadir sensor ERT was a poor measure of vegetation temperature due to effects of soil emissions. Strong vertical temperature profiles of vegetation components were also observed. The theory and measurements presented document that remote measurements of vegetation canopy temperatures cannot be made indiscriminately over large spatial regions without consideration of the underlying physical principles.

Kimes, D. S.↗

A first-order time-domain Green's function approach to supersonic unsteady flow

A time-domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. The Green's function method is employed in order to convert the potential-flow differential equation into an integral one. This integral equation is then discretized, in space through standard finite-element technique, and in time through finite-difference, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. Consistent with the spatial linear (first-order) finite-element approximations, the potential and co-normalwash are assumed to vary linearly in time. The long range goal of our research is to develop a comprehensive theory for unsteady supersonic potential aerodynamics which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗

Coronal heating by the resonant absorption of Alfven waves: Wavenumber scaling laws.

The importance of global modes in coronal loop heating is well established. In the present work the scaling of the global-mode resonant heating rate with the perturbation wavenumbers is studied with the numerical solution of the linearized time-dependent MHD equations for a full compressible, low-beta, resistive plasma using an implicit integration scheme. The numerical simulations demonstrate that the dissipation on inhomogeneties in the background Alfven speed occurs in narrow resonant layer with the highest heating rate at the global-mode frequency. The global-mode heating rate H (sub r) was found to scale as H (Sub r) approximately k (sub y) (exp 1.03) when k (sub z) = 0.1, and as H (sub r) approximately k (sub y) (exp -1.93) when k (sub z) = 0.75, where k (sub y) and k (sub z) are the wavenumbers in the perpendicular and parallel to the magnetic field directions, respectively, while the dependence of H (sub r) on k (sub z) is more complex. The quality factor Q of the MHD resonance cavity scales as Q approximately k (sub y) (exp -1.8) for k (sub z) = 0.75 and as Q approximately k (sub y) (exp -1.46) for k (sub z) = 0.1. The numerically determined heating rate scaling, the global-mode fequency, and the quality factor are in good agreement with the analytical linear theory. The magnitude of the perturbed velocities was found to decrease with k (sub y). Assuming typical coronal loop parameters (B (sub 0) = 100-200 G, upsilon (sub A) = 2000-4000 km/s), the Alfven waves can supply the required heating to a low-Q loops.

Ofman, L.↗

Aberrations for Grazing Incidence Optics

Large number of grazing incidence telescope configurations have been designed and studied. Wolte1 telescopes are commonly used in astronomical applications. Wolter telescopes consist of a paraboloidal primary mirror and a hyperboloidal or an ellipsoidal secondary mirror. There are 8 possible combinations of Wolter telescopes. Out of these possible designs only type 1 and type 2 telescopes are widely used. Type 1 telescope is typically used for x-ray applications and type 2 telescopes are used for EUV applications. Wolter-Schwarzshild (WS) telescopes offer improved image quality over a small field of view. The WS designs are stigmatic and free of third order coma and, therefore, the PSF is significantly better over a small field of view. Typically the image is more symmetric about its centroid. As for the Wolter telescopes there are 8 possible combinations of WS telescopes. These designs have not been widely used because the surface equations are complex parametric equations complicating the analysis and typically the resolution requirements are too low to take full advantage of the WS designs. There are several other design options. Most notable are wide field x-ray telescope designs. Polynomial designs were originally suggested by Burrows4 and hyperboloid-hyperboloid designs for solar physics applications were designed by Harvey5. No general aberration theory exists for grazing incidence telescopes that would cover all the design options. Several authors have studied the aberrations of grazing incidence telescopes. A comprehensive theory of Wolter type 1 and 2 telescopes has been developed. Later this theory was expanded to include all possible combinations of grazing incidence and also normal incidence paraboloid-hyperboloid and paraboloid-ellipsoid telescopes. In this article the aberration theory of Wolter type telescopes is briefly reviewed.

Saha, Timo T.↗

A new method to real-normalize measured complex modes

A time domain subspace iteration technique is presented to compute a set of normal modes from the measured complex modes. By using the proposed method, a large number of physical coordinates are reduced to a smaller number of model or principal coordinates. Subspace free decay time responses are computed using properly scaled complex modal vectors. Companion matrix for the general case of nonproportional damping is then derived in the selected vector subspace. Subspace normal modes are obtained through eigenvalue solution of the (M sub N) sup -1 (K sub N) matrix and transformed back to the physical coordinates to get a set of normal modes. A numerical example is presented to demonstrate the outlined theory.

Wei, Max L.↗

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗

Theoretical Study of Infrared and Raman Spectra of Hydrated Magnesium Sulfate Salts

Harmonic and anharmonic vibrational frequencies, as well as infrared and Raman intensities, are calculated for MgSO4.nH20 (n=1-3). Electronic structure theory at the second order Moller-Plesset perturbation theory (MP2) level with a triple-zeta + polarization (TZP) basis set is used to determine the geometry, properties, and vibrational spectra of pure and hydrated MgSO4 salts. The direct vibrational self-consistent field (VSCF) method and its correlation corrected (CC-VSCF) extension are used to determine anharmonic corrections to vibrational frequencies and intensities for the pure MgSO4 and its complex with one water molecule. Very significant differences are found between vibrational of water molecules in complexes with MgSO4 and pure water. Some of the O-H stretching frequencies are shifted to the red very significantly (by up to 1500-2000/cm) upon complexation with magnesium sulfate. They should be observed between 1700 and 3000/cm in a region very different from the corresponding O-H stretch frequency region of pure water (3700-3800/cm). In addition, the SO2 stretching vibrations are found at lower frequency regions than the water vibrations. They can serve as unique identifiers for the presence of sulfate salts. The predicted infrared and Raman spectra should be of valuable help in the design of future missions and analysis of observed data from the ice surface of Jupiter's moon Europa that possibly contains hydrated MgSO4 salts.

Chaban, Galina M.↗

Modeling Distributed Situation Awareness in Resilience-Based Design of Complex Engineered Systems

Human operators play a major role in the resilience of complex systems–while human error is one of the biggest contributors to hazardous events, operators additionally play a critical role in mitigating hazardous events. A key factor underlying this operator resilience is situation awareness–the ability of operators to understand their environment and each other to achieve desired system functions. In contrast to situation awareness-related accident models in the literature, which are largely conceptual in nature, this work proposes the use of a dynamic simulation framework to concretely model both the effects of situation awareness-related human errors and situation awareness-related hazard-mitigating properties using the distributed situation awareness theory. This work then presents specialized model constructs to enable agents’ individual perceptions of the system state and transactions with other agents (and thus distributed situation awareness) to be represented in simulation. To demonstrate this framework, it is then adapted to an aircraft taxiway case study, where it is used to model aircraft conflicts due to lack of vision and poor communications from the air traffic controller. This demonstration shows the potential of using simulation models to rigorously understand situation awareness-related human errors and thus inform the design of resilience.

Resilience Modeling↗

Modeling Distributed Situation Awareness in Resilience-based Design of Complex Systems

Human operators play a major role in the resilience of complex systems–while human error is one of the biggest contributors to hazardous events, operators additionally play a critical role in mitigating hazardous events. A key factor underlying this operator resilience is situation awareness–the ability of operators to understand their environment and each other to achieve desired system functions. In contrast to situation awareness-related accident models in the literature, which are largely conceptual in nature, this work proposes the use of a dynamic simulation framework to concretely model both the effects of situation awareness-related human errors and situation awareness-related hazard-mitigating properties using the distributed situation awareness theory. This work then presents specialized model constructs to enable agents’ individual perceptions of the system state and transactions with other agents (and thus distributed situation awareness) to be represented in simulation. To demonstrate this framework, it is then adapted to an aircraft taxiway case study, where it is used to model aircraft conflicts due to lack of vision and poor communications from the air traffic controller. This demonstration shows the potential of using simulation models to rigorously understand situation awareness-related human errors and thus inform the design of resilience.

Resilience Modeling↗

High-Tc Superconducting Bolometer Noise Measurement Using Low Noise Transformers - Theory and Optimization

Care must always be taken when performing noise measurements on high-Tc superconducting materials to ensure that the results are not from the measurement system itself. One situation likely to occur is with low noise transformers. One of the least understood devices, it provides voltage gain for low impedance inputs (< 100 ), e.g., YBaCuO and GdBaCuO thin films, with comparatively lower noise levels than other devices for instance field effect and bipolar junction transistors. An essential point made in this paper is that because of the complex relationships between the transformer ports, input impedance variance alters the transformer s transfer function in particular, the low frequency cutoff shift. The transfer of external and intrinsic transformer noise to the output along with optimization and precautions are treated; all the while, we will cohesively connect the transfer function shift, the load impedance, and the actual noise at the transformer output.

Aslam, Shahid↗

Probabilistic structural analysis to evaluate the structural durability of SSME critical components

NASA Lewis Research Center is currently developing probabilistic structural analysis methods for select Space Shuttle Main Engine (SSME) structural components. Briefly, the deterministic, three-dimensional, inelastic analysis methodology developed under the Hot Section Technology (HOST) and R and T Base Programs is being augmented to accommodate the complex probabilistic loading spectra, the thermoviscoplastic material behavior, and the material degradation associated with the environment of space propulsion system structural components representative of the SSME such as turbine blades, transfer ducts, and liquid-oxygen posts. The development of probabilistic structural analysis methodology consists of the following program elements: (1) composite load spectra; (2) probabilistic structural analysis methods; (3) probabilistic finite element theory - new variational principles; and (4) probabilistic structural analysis application. In addition, the program includes deterministic analysis elements: (1) development of structural tailoring computer codes (SSME/STAEBL); (2) development of dynamic creep buckling/ratcheting theory; (3) evaluation of the dynamic characteristics of single-crystal SSME blades; (4) development of SSME blade damper technology; and (5) development of integrated boundary elements for hotfluid structure interaction.

Chamis, Christos C.↗

Nonlinear flight test trajectory controllers for aircraft

Flight test trajectory control systems are designed to enable the pilot to follow complex trajectories for evaluating an aircraft within its known flight envelope and to explore the boundaries of its capabilities. Previous design approaches were baed on linearized aircraft models necessitating a large amount of data storage along with gain schedules. In this paper, the synthesis of nonlinear flight test trajectory controllers for a fixed-wing aircraft is described. This approach uses singular perturbation theory and the recently developed theory of prelinearizing transformations. These controllers do not require gain scheduling for satisfactory operation, can be used in arbitrarily nonlinear maneuvers, and are mechanized with a direct, noniterative analytic solution.

Menon, P. K. A.↗

Membership generation using multilayer neural network

There has been intensive research in neural network applications to pattern recognition problems. Particularly, the back-propagation network has attracted many researchers because of its outstanding performance in pattern recognition applications. In this section, we describe a new method to generate membership functions from training data using a multilayer neural network. The basic idea behind the approach is as follows. The output values of a sigmoid activation function of a neuron bear remarkable resemblance to membership values. Therefore, we can regard the sigmoid activation values as the membership values in fuzzy set theory. Thus, in order to generate class membership values, we first train a suitable multilayer network using a training algorithm such as the back-propagation algorithm. After the training procedure converges, the resulting network can be treated as a membership generation network, where the inputs are feature values and the outputs are membership values in the different classes. This method allows fairly complex membership functions to be generated because the network is highly nonlinear in general. Also, it is to be noted that the membership functions are generated from a classification point of view. For pattern recognition applications, this is highly desirable, although the membership values may not be indicative of the degree of typicality of a feature value in a particular class.

Kim, Jaeseok↗

Improving Rotor-Stator Interaction Noise Code Through Analysis of Input Parameters

There are two major sources of aircraft noise. The first is from the airframe and the second is from the engines. The focus of the acoustics branch at NASA Glenn is on the engine noise sources. There are two major sources of engine noise; fan noise and jet noise. Fan noise, produced by rotating machinery of the engine, consists of both tonal noise, which occurs at discrete frequencies, and broadband noise, which occurs across a wide range of frequencies. The focus of my assignment is on the broadband noise generated by the interaction of fan flow turbulence and the stator blades. such as the sweep and stagger angles and blade count, as well as the flow parameters such as intensity of turbulence in the flow. The tool I employed in this work is a computer program that predicts broadband noise from fans. The program assumes that the complex shape of the curved blade can be represented as a single flat plate, allowing it to use fairly simple equations that can be solved in a reasonable amount of time. While the results from such representation provided reasonable estimates of the broadband noise levels, they did not usually represent the entire spectrum accurately. My investigation found that the discrepancy between data and theory can be improved if the leading edge and the trailing edge of the blade are treated separately. Using this approach, I reduced the maximum error in noise level from a high of 30% to less than 5% for the cases investigated. Detailed results of this investigation will be discussed at my presentation. The objective of this study is to investigate the influence of geometric parameters

Unton, Timothy J.↗

Far-field Lorenz-Mie Scattering in an Absorbing Host Medium. II: Improved Stability of the Numerical Algorithm

A recently developed FORTRAN program computing far-field optical observables for spherical particles in an absorbing medium has exhibited numerical instability arising when the product of the particle vacuum size parameter and the imaginary part of the refractive index of the host becomes sufficiently large. We offer a simple analytical explanation of this instability and propose a compact numerical algorithm for the stable computation of Lorenz–Mie coefficients based on an upward recursion formula for spherical Hankel functions of a complex argument. Extensive tests confirm an excellent accuracy of this algorithm approaching machine precision. The improved public-domain FORTRAN program is available at https://www.giss.nasa.gov/staff/mmishchenko/Lorenz-Mie.html.

Lorenz-Mie theory↗