SEARCH · Search NASA

Results for “Complete equation of state”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

NASA NTRS1976

Time scales for molecule formation by ion-molecule reactions

Analytical solutions are obtained for nonlinear differential equations governing the time-dependence of molecular abundances in interstellar clouds. Three gas-phase reaction schemes are considered separately for the regions where each dominates. The particular case of CO, and closely related members of the Oh and CH families of molecules, is studied for given values of temperature, density, and the radiation field. Nonlinear effects and couplings with particular ions are found to be important. The time scales for CO formation range from 100,000 to a few million years, depending on the chemistry and regime. The time required for essentially complete conversion of C(+) to CO in the region where the H3(+) chemistry dominates is several million years. Because this time is longer than or comparable to dynamical time scales for dense interstellar clouds, steady-state abundances may not be observed in such clouds.

Langer, W. D.↗

NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗

NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

NASA NTRS1976

The ionization of a low-density intercloud medium by a single O star

The ionization structure of hydrogen in a cloud-containing low-density gas surrounding an O star is calculated, emphasizing how UV photons cause the ionization of hydrogen in regions which are totally obscured from the star by neutral clouds. Particular consideration is given to the ionization of such obscured regions by diffuse Lyman continuum radiation produced when H(+) recombines directly to the ground state. The assumed physical structure of the cloud-containing neighborhood of an O star is discussed, a method is outlined for expressing the emission and absorption properties of the neutral clouds and their shadows in terms of effective volume emissivities and absorption coefficients on a line of sight, and equations of radiative transfer are derived for diffuse and stellar radiation in a cloud-containing H II region. Reradiation of Lyman continuum photons from the ionized boundary layer of a neutral cloud or shadow is analyzed. The expected emission measure across a large cloud-containing H II region is determined and found to compare well with that obtained from observations of H-alpha emission around runaway O stars. It is shown that the hydrogen in cloud shadows may be completely ionized by diffuse radiation within approximately half the Stroemgren radius, depending on the extent of that radius, the average frequency of the Lyman continuum radiation from the star, and the average cloud size.

Elmegreen, B. G.↗

NASA NTRS2008

Quantum-Inspired Maximizer

A report discusses an algorithm for a new kind of dynamics based on a quantum- classical hybrid-quantum-inspired maximizer. The model is represented by a modified Madelung equation in which the quantum potential is replaced by different, specially chosen 'computational' potential. As a result, the dynamics attains both quantum and classical properties: it preserves superposition and entanglement of random solutions, while allowing one to measure its state variables, using classical methods. Such optimal combination of characteristics is a perfect match for quantum-inspired computing. As an application, an algorithm for global maximum of an arbitrary integrable function is proposed. The idea of the proposed algorithm is very simple: based upon the Quantum-inspired Maximizer (QIM), introduce a positive function to be maximized as the probability density to which the solution is attracted. Then the larger value of this function will have the higher probability to appear. Special attention is paid to simulation of integer programming and NP-complete problems. It is demonstrated that the problem of global maximum of an integrable function can be found in polynomial time by using the proposed quantum- classical hybrid. The result is extended to a constrained maximum with applications to integer programming and TSP (Traveling Salesman Problem).

Zak, Michail↗

NASA NTRS2003

Theoretical Analysis of Thermodynamic Measurements near a Liquid-Gas Critical Point

Over the years, many ground-based studies have been performed near liquid-gas critical points to elucidate the expected divergences in thermodynamic quantities. The unambiguous interpretation of these studies very near the critical point is hindered by a gravity-induced density stratification. However, these ground-based measurements can give insight into the crossover behavior between the asymptotic critical region near the transition and the mean field region farther away. We have completed a detailed analysis of heat capacity, susceptibility and coexistence curve measurements near the He-3 liquid-gas critical point using the minimal-subtraction renormalization (MSR) scheme within the phi(exp 4) model. This MSR scheme, using only two adjustable parameters, provides a reasonable global fit to all of these experimental measurements in the gravity-free region out to a reduced temperature of |t| approx. 2x10(exp -2). Recently this approach has also been applied to the earlier microgravity measurements of Haupt and Straub in SF(sub 6) with surprising results. The conclusions drawn from the MSR analyses will be presented. Measurements in the gravity-affected region closer to the He-3 critical point have also been analyzed using the recent crossover parametric model (CPM) of the equation-of-state. The results of fitting heat capacity measurements to the CPM model along the He-3 critical isochore in the gravity-affected region will also be presented.

Barmatz, M.↗

NASA NTRS1960

Experimental Determination of the Effects of Frequency and Amplitude of Oscillation on the Roll-Stability Derivatives for a 60 Degree Delta-Wing Airplane Model

A 60' delta-wing airplane model was oscillated in roll for several frequencies and amplitudes of oscillation to determine the effects of the oscillatory motion on the roll-stability derivatives for the model. The derivatives were measured at a Reynolds number of 1,600,000 for the wing alone, the wing-fuselage combination, and the complete model which included a triangular-plan-form vertical tail. Both rolling and yawing moments due to rolling velocity exhibited large frequency effects for angles of attack higher than 16 degrees. Variations in these derivatives were measured for the lowest frequencies of oscillation; as the frequency increased, the derivatives because more nearly linear with angle of attack. Both velocity derivatives were considerably different at high angles of attack from the corresponding derivatives measured by the steady-state rolling-flow technique. Rolling and yawing moments due to rolling acceleration were measured and similarly found to be highly dependent on frequency at high angles of attack. Some period and time-to-damp computations, which were made to reveal the significance of the acceleration derivatives, indicated that inclusion of the measured derivatives in the equations of motion lengthened the period of the lateral oscillation by 10 percent for a typical delta-wing airplane and increased the time to damp to one-half amplitude by 50 percent.

Fisher, Lewis R.↗

NASA NTRS2016

Trajectory Design Employing Convex Optimization for Landing on Irregularly Shaped Asteroids

Mission proposals that land on asteroids are becoming popular. However, in order to have a successful mission the spacecraft must reliably and softly land at the intended landing site. The problem under investigation is how to design a fuel-optimal powered descent trajectory that can be quickly computed on- board the spacecraft, without interaction from ground control. An optimal trajectory designed immediately prior to the descent burn has many advantages. These advantages include the ability to use the actual vehicle starting state as the initial condition in the trajectory design and the ease of updating the landing target site if the original landing site is no longer viable. For long trajectories, the trajectory can be updated periodically by a redesign of the optimal trajectory based on current vehicle conditions to improve the guidance performance. One of the key drivers for being completely autonomous is the infrequent and delayed communication between ground control and the vehicle. Challenges that arise from designing an asteroid powered descent trajectory include complicated nonlinear gravity fields, small rotating bodies and low thrust vehicles. There are two previous studies that form the background to the current investigation. The first set looked in-depth at applying convex optimization to a powered descent trajectory on Mars with promising results.1, 2 This showed that the powered descent equations of motion can be relaxed and formed into a convex optimization problem and that the optimal solution of the relaxed problem is indeed a feasible solution to the original problem. This analysis used a constant gravity field. The second area applied a successive solution process to formulate a second order cone program that designs rendezvous and proximity operations trajectories.3, 4 These trajectories included a Newtonian gravity model. The equivalence of the solutions between the relaxed and the original problem is theoretically established. The proposed solution for designing the asteroid powered descent trajectory is to use convex optimization, a gravity model with higher fidelity than Newtonian, and an iterative solution process to design the fuel optimal trajectory. The solution to the convex optimization problem is the thrust profile, magnitude and direction, that will yield the minimum fuel trajectory for a soft landing at the target site, subject to various mission and operational constraints. The equations of motion are formulated in a rotating coordinate system and includes a high fidelity gravity model. The vehicle's thrust magnitude can vary between maximum and minimum bounds during the burn. Also, constraints are included to ensure that the vehicle does not run out of propellant, or go below the asteroid's surface, and any vehicle pointing requirements. The equations of motion are discretized and propagated with the trapezoidal rule in order to produce equality constraints for the optimization problem. These equality constraints allow the optimization algorithm to solve the entire problem, without including a propagator inside the optimization algorithm.

Pinson, Robin M.↗

NASA NTRS1991

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗

NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

NASA NTRS2009

Simulating Bone Loss in Microgravity Using Mathematical Formulations of Bone Remodeling

Most mathematical models of bone remodeling are used to simulate a specific bone disease, by disrupting the steady state or balance in the normal remodeling process, and to simulate a therapeutic strategy. In this work, the ability of a mathematical model of bone remodeling to simulate bone loss as a function of time under the conditions of microgravity is investigated. The model is formed by combining a previously developed set of biochemical, cellular dynamics, and mechanical stimulus equations in the literature with two newly proposed equations; one governing the rate of change of the area of cortical bone tissue in a cross section of a cylindrical section of bone and one governing the rate of change of calcium in the bone fluid. The mechanical stimulus comes from a simple model of stress due to a compressive force on a cylindrical section of bone which can be reduced to zero to mimic the effects of skeletal unloading in microgravity. The complete set of equations formed is a system of first order ordinary differential equations. The results of selected simulations are displayed and discussed. Limitations and deficiencies of the model are also discussed as well as suggestions for further research.

Pennline, James A.↗

NASA NTRS1968

Predicted vibration responses of Apollo structure and effects of pressure correlation lengths on response

Progressive wove and reverberant acoustic fields ore often used in the laboratory for qualification testing of aerospace structures that are subjected to unsteady aerodynamic environments which occur during flight through the atmosphere. Acoustic simulation of these in-flight environments i s achieved i f the structural vibration response to acoustic excitation is equivalent, over the frequency range of interest, to the response caused by aerodynamic turbulence. To determine the degree of response simulation that could be expected for such qualification testing of the Apollo Spacecraft, responses of a segment of the Apollo structure were estimated theoretically for both acoustic and flight environments. The vibration analysis presented in this report was performed for the Spacecraft Lunar Module Adaptor; and the three environments treated are turbulence at Mach 1 and Mach 2, a reverberant acoustic field, and a modified acoustic progressive wave field. The progressive wave field is assumed to be generated within a specially designed shroud of sixteen axially oriented ducts which is constructed around and contoured to the external skin of the vehicle. Each duct is open on the side adjacent to the skin and is driven acoustically at one end by an independent noise source, thus allowing for different pressure correlation patterns around the circumference of the vehicle. The purpose of the analysis was to determine the effects of surface pressure correlation lengths on the response characteristics of representative Apollo structure, and to select on this basis, the optimum acoustic environment for structural vibration qualification tests of the Apollo in the Spacecraft Acoustic Laboratory at MSC. The method of analysis employed consists of determining, independently, the space average, mean-square acceleration spectral density for each of the significant modes of vibration of an equivalent cylindrical shell, and then expressing the total structural response as a linear summation of these modal responses. Approximately 570 classical modes of a pinned-end cylinder were used in order to find the spectral response characteristics of the structure throughout the frequency range of 10 - 1,000 Hz. The effective forcing functions or joint acceptances of the various modes of vibration were computed and are presented graphically for each of the pressure excitation environments. A brief development of the response equations used is presented in the report. Space average acceleration response spectra were computed for twelve different cases including aerodynamic turbulence at Mach 1 and Mach 2, reverberant acoustic field, five duct correlation patterns, two cases of an axially damped progressive wave field, and two cases in which the structural damping constant was varied from the expected value. The response spectra were computed and plotted by a high speed digital computer far a flat excitation pressure spectrum; and the response spectra were numerically integrated to give one-third octave bond average responses. A complete set of digital computer programs has been developed for both flat plate and cylindrical shell structures for the above environments, and each requires run times of three minutes or less. Thus, the work presented herein, along with the results, constitutes a practical advancement to the existing state-of-the-art of vibration prediction.

R W White↗

NASA NTRS1994

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.

Murthy, V. R.↗

NASA NTRS1986

Combustion research in the Internal Fluid Mechanics Division

The goal of this research is to bring computational fluid dynamics to a state of practical application for the aircraft engine industry. The approach is to have a strongly integrated computational and experimental program for all the disciplines associated with the gas turbine and other aeropropulsion systems by advancing the understanding of flow physics, heat transfer, and combustion processes. The computational and experimental research is integrated in the following way: the experiments that are performed provide an empirical data set so that physical models can be formulated to describe the processes that are occurring - for example, turbulence or chemical reaction. These experiments also form a data base for those who are doing code development by providing experimental data against which the codes can be verified and assesed. Models are generated as closure to some of the numerical codes, and they also provide physical insight for experiments. At the same time, codes which solve the complete Navier-Stokes equations can be used as a kind of numerical experiment from which far more extensive data can be obtained than ever could be obtained experimentally. This could provide physical insight into the complex processes that are taking place. These codes are also exercised against experimental data to assess the accuracy and applicability of models.

Mularz, Edward J.↗

NASA NTRS2011

Introduction to: Atlantic Meridional Overturning Circulation(AMOC)

A striking conclusion of the Intergovernmental Panel on Climate Change 2007 report is the crucial role that the Atlantic Meridional Overturning Circulation (AMOC) may play in anthropogenic climate change. However, these IPCC coupled climate simulations show a broad range of uncertainty in the magnitude and timing of AMOC transport change ranging from none to nearly complete collapse within the 21st century. The potential consequences of large changes in the characteristics of AMOC have motivated the creation in the United States of an interagency program and implementation plan to develop monitoring and prediction capabilities for the AMOC This program parallels the development of substantial monitoring efforts by European, South American and African countries -- notably the UK Rapid and Rapid-Watch programs. The papers contained in this volume are derived from presentations at the First U.S. Atlantic Meridional Overturning Circulation (AMOC) Meeting held 4 - 6 May, 2009 to review the US implementation plan and its coordination with other monitoring activities. The Atlantic Meridional Overturning Circulation consists of multiple components illustrated in an attached figure. Water enters the South Atlantic at upper and intermediate depths through both western and eastern routes (where eddy transport is especially important) and is transported northward across the equator, where it recirculates within the northern subtropical and subpolar gyres. The northern end is defined by the sinking regions of the Nordic Seas and the Labrador Sea where the waters that eventually form the upper and lower branches of North Atlantic Deep Water are conditioned. High surface salinities, the result of high net evaporation in the tropics and subtropics (including the Mediterranean Sea), and presence of regions of the Arctic Ocean that remain ice-free even in winter allow for the rapid cooling and thus densification of surface water. This dense surface water becomes the source of deep water formation in the sinking regions. In addition to transporting mass, the AMOC transports roughly half of the total amount of heat carried northward through the northern subtropics (down the temperature-gradient) by the ocean. In contrast in the Southern Hemisphere AMOC transports heat up-gradient from the cool Circumpolar Current to the warm tropics. Paleoevidence suggests that AMOC heat transport in the two hemispheres has varied over time in ways intimately tied to millennial changes in the Earth's climate. In one example, the abrupt Younger Dryas spell of cold weather over the North Atlantic, which began 13,000 years ago, has generally been linked to a millennial shutdown of the AMOC as a result of massive freshwater discharge from the North American continent. The current AMOC monitoring array consists of a series of instrumented transects located across key passages (see Cunningham et al., 2010 for a recent review). In the Arctic and sub-Arctic, transects cross Fram Strait, Denmark Strait and the Faroe Channel (connecting Greenland, Iceland, and the United Kingdom), as well as the entrance to the Labrador Sea. Further south and extending outwards from the east coast of North America there are a series of monitoring arrays including arrays of the Canadian Atlantic Zone Monitoring Program, deployments of the Rapid Western Atlantic Variability Experiment (WAVE), Line W at 39 N, as well as the Rapid-MOC moored array. The latter spans the entire Atlantic basin along 26.5 N. At tropical latitudes we have the Meridional Overturning Variability Experiment (MOVE) array at 16 N, while in the Southern Hemisphere a corresponding basin-spanning transect is being established at the latitude of Cape of Good Hope, complemented by arrays at Drake Passage.

Hakkinen, Sirpa↗

NASA NTRS1979

Flight test evaluation of predicted light aircraft drag, performance, and stability

A technique was developed which permits simultaneous extraction of complete lift, drag, and thrust power curves from time histories of a single aircraft maneuver such as a pullup (from V sub max to V sub stall) and pushover (to sub V max for level flight.) The technique is an extension to non-linear equations of motion of the parameter identification methods of lliff and Taylor and includes provisions for internal data compatibility improvement as well. The technique was show to be capable of correcting random errors in the most sensitive data channel and yielding highly accurate results. This technique was applied to flight data taken on the ATLIT aircraft. The drag and power values obtained from the initial least squares estimate are about 15% less than the 'true' values. If one takes into account the rather dirty wing and fuselage existing at the time of the tests, however, the predictions are reasonably accurate. The steady state lift measurements agree well with the extracted values only for small values of alpha. The predicted value of the lift at alpha = 0 is about 33% below that found in steady state tests while the predicted lift slope is 13% below the steady state value.

Smetana, F. O.↗

NASA NTRS2019

Pulsed Plasma Acceleration Modeling in Detonation and Deflagration Modes

Gas-fed electromagnetic pulsed plasma accelerators discharge electrical energy into a gas to ionize and electromagnetically accelerate propellant in the domain. Efforts to model pulsed accelerators have either assumed that the discharge is short and completely transient, accelerating the gas by entraining it in a moving current sheet, or that the discharge is relatively long, establishing a stable quasi-steady current distribution pattern that accelerates a plasma flowing through it. There have been pulsed plasma accelerator tests that appear to fit somewhere between these two bounds, exhibiting some properties that are associated with the purely transient devices while also showing others that are associated with quasi-steady-state plasma acceleration. A model is presented based upon the premise that all pulsed plasma accelerators first form an accelerating current sheet (detonation mode accelerator). Depending upon the pulse length and the gas conditions in the dis- charge channel, the plasma sheet may reach the end of the accelerator before the discharge has completed a full half-cycle, wherein the proposed model transitions to a quasi-steady description of the acceleration process (deflagration mode accelerator). A review of the entire model is presented, highlighting improvements and upgrades implemented to aid in the stability of the numerical scheme used to model the gas flow in the channel and an improved treatment of current sheet mass shedding, which adds gas to the wake of the sheet. The assumptions employed to model the transition from detonation to deflagration mode are presented and used to generate solutions to the governing equations. The modeling of the deflagration-mode under the present assumptions results in very low deflagration impulse bits relative to those obtained at the end of the detonation mode, implying that the present assumptions and modeling of the deflagration mode may not represent a fruitful approach.

Polzin, Kurt A.↗