Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling 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 379 records · Page 21

MFLOP to GFLOP: The Impact on High Fidelity Based Computational Aeroelasticity

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft which can experience complex flow/structure interactions. Even at moderate angles of attack supersonic aircraft can experience vortex induced aeroelastic oscillations. Near transonic speeds buffet associated structural oscillations are possible. Aircraft flying in transonic regime may experience a dip in the flutter speed. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overco me the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations. Modifications required in both fluids and structural solvers in order to run efficiently on parallel computers will be discussed. Implementation of moving grids and fluid/structural interface on parallel computers will be discussed.

Guruswamy, Guru P.↗

How the magnetosphere is driven into the substorm

This paper puts together several aspects of the substorm phenomena into a scenario subject to quantitative analysis, based on Ampere's law, magnetic stresses, and the effect of the parallel electric field. As the magnetospheric convection continuously shears the magnetic field, the parallel electric field suddenly switches on when the shear reaches a critical value. The parallel electric field not only accelerates the charged particles to cause the auroral display but also partially decouples the ionosphere from the magnetosphere so that the magnetospheric convection is enhanced significantly, causing the substorm expansion phase. Pi 2 pulsations are excited by the change of the magnetosphere-ionosphere coupling. We show a direct relationship between the parallel potential drop, the enhanced convection rate, and the amplitude of the Pi 2 pulsations. With no adjustable parameter, the relationship shows quantitative agreement with the diverse observations.

Zhu, Xiaoming↗

Extraction technique for the determination of oxygen-18 in water using preevacuated glass vials

The need for a rapid, inexpensive technique for routine 18O/16O extraction from water has arisen recently through applications in the medical sciences and in hydrology. The traditional experimental technique for determining the oxygen isotope composition of water, the CO2-H2O equilibration method, is tedious, time consuming, and involves the use of custom-made glass apparatus. Furthermore, because of potential memory effects from one sample to the next, the glassware needs to be thoroughly cleaned between runs. A few attempts have been made to improve upon the method. Attempts to analyze water directly in the source of the mass spectrometer produced large memory effects and questionable results. Commercially available apparatus for automated extraction of 18O/16O from water is generally prohibitively expensive and often is designed to interface only with the manufacturer's own mass spectrometer. The method described in this paper utilizes inexpensive, off-the-shelf, preevacuated, glass vials. Preevacuated vials have been used by others for the isotopic analysis of breath CO2 and are well tested. The vials can be purchased in bulk from scientific apparatus suppliers at a relatively low cost. These are coupled with a simplified extraction line consisting of a stainless steel syringe needle and a glass cold trap. Vials are filled with CO2 and H2O and shaken in a constant-temperature water bath for at least 90 min. Since the vials are discarded after use, no cleaning is necessary, essentially eliminating any memory effect. Reproducibility is generally better than +/- 0.05%. The only reagents required are gaseous CO2 for equilibration, a dry ice/alcohol mixture for trapping water, and liquid nitrogen for transferring the CO2.

NASA Discipline Number 52-10↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗

Progress in Computational Aeroelasticity Using High Fidelity Flow and Structural Equations on Parallel Computers

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENASERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru P.↗

NASTRAN solutions of problems described by simultaneous parabolic differential equations

NASTRAN solution techniques are shown for a numerical analysis of a class of coupled vector flow processes described by simultaneous parabolic differential equations. To define one physical problem type where equations of this form arise, the differential equations describing the coupled transfers of heat and mass in mechanical equilibrium with negligible mass average velocity are presented and discussed. Also shown are the equations describing seepage when both electrokinetic and hydrodynamic forces occur. Based on a variational statement of the general problem type, the concepts of scalar transfer elements and parallel element systems are introduced. It is shown that adoptation of these concepts allows the direct use of NASTRAN's existing Laplace type elements for uncoupled flow (the heat transfer elements) for treating multicomponent coupled transfer. Sample problems are included which demonstrate the application of these techniques for both steady-state and transient problems.

Mason, J. B.↗

An evaluation of the directed flow graph methodology

The applicability of the Directed Graph Methodology (DGM) to the design and analysis of special purpose image and signal processing hardware was evaluated. A special purpose image processing system was designed and described using DGM. The design, suitable for very large scale integration (VLSI) implements a region labeling technique. Two computer chips were designed, both using metal-nitride-oxide-silicon (MNOS) technology, as well as a functional system utilizing those chips to perform real time region labeling. The system is described in terms of DGM primitives. As it is currently implemented, DGM is inappropriate for describing synchronous, tightly coupled, special purpose systems. The nature of the DGM formalism lends itself more readily to modeling networks of general purpose processors.

Snyder, W. E.↗

Surface temperature/heat transfer measurement using a quantitative phosphor thermography system

A relative-intensity phosphor thermography technique developed for surface heating studies in hypersonic wind tunnels is described. A direct relationship between relative emission intensity and phosphor temperature is used for quantitative surface temperature measurements in time. The technique provides global surface temperature-time histories using a 3-CCD (Charge Coupled Device) video camera and digital recording system. A current history of technique development at Langley is discussed. Latest developments include a phosphor mixture for a greater range of temperature sensitivity and use of castable ceramics for inexpensive test models. A method of calculating surface heat-transfer from thermal image data in blowdown wind tunnels is included in an appendix, with an analysis of material thermal heat-transfer properties. Results from tests in the Langley 31-Inch Mach 10 Tunnel are presented for a ceramic orbiter configuration and a four-inch diameter hemisphere model. Data include windward heating for bow-shock/wing-shock interactions on the orbiter wing surface, and a comparison with prediction for hemisphere heating distribution.

Buck, G. M.↗

A reciprocity formulation for the EM scattering by an obstacle within a large open cavity

A formulation based on a generalized reciprocity theorem is developed for analyzing the external high frequency EM scattering by a complex obstacle inside a relatively arbitrary open-ended waveguide cavity when it is illuminated by an external source. This formulation is also extended to include EM fields whose time dependence may be nonperiodic. A significant advantage of this formulation is that it allows one to break up the analysis into two independent parts; one deals with the waveguide cavity shape alone and the other with the obstacle alone. The external scattered field produced by the obstacle (in the presence of the waveguide cavity structure) is given in terms of a generalized reciprocity integral over a surface S(T) corresponding to the interior waveguide cavity cross section located conveniently but sufficiently close to the obstacle. Furthermore, the fields coupled into the cavity from the source in the exterior region generally need to propagate only one-way via the open front end (which is directly illuminated) to the interior surface S(T) in this approach, and not back, in order to find the external field scattered by the obstacle.

Pathak, Prabhakar H.↗

Fluid/Structure Interaction Studies of Aircraft Using High Fidelity Equations on Parallel Computers

Abstract Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional super computers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru↗

Prediction and measurement of the proportionality constant in statistical energy analysis of structures.

The fundamental equation of statistical energy analysis (SEA) states that the average power flow between two coupled vibrating systems is proportional to the difference in their average modal energies. Under certain circumstances it is possible to estimate the proportionality constant by modifying system boundary conditions on the separated systems and calculating or measuring changes in the systems. Newland's estimate, based upon blocking part of the system, is reexamined, and limitations are discussed. Three alternative methods which circumvent blocking are presented. These were applied to predict power flow in experiments on coupled beams and on coupled plates wherein power flow through the coupling was measured directly as a product of force times velocity. The measurements support the fundamental SEA relation, including the null power point where the average modal energies are equal.

Lotz, R.↗

Calculation of unsteady airfoil loads with and without flap deflection at -90 degrees incidence

A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This unique method provides for the direct solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body-fitted computational mesh incorporating a staggered grid method. The vorticity is determined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and for the conservation of mass at the mesh-cell centers. The method provides for the direct solution of the flow field and satisfies the conservation of mass to machine zero at each time-step. The results of the present analysis and experimental results obtained for a XV-15 airfoil are compared. The comparisons indicate that the calculated drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results. Comparisons of the numerical results of the present method for several airfoils demonstrate the significant influence of airfoil curvature and flap deflection on the predicted download.

Stremel, Paul M.↗

Solving Modal Equations of Motion with Initial Conditions Using MSC/NASTRAN DMAP: Coupled Versus Uncoupled Integration - Part 2

By utilizing MSC/NASTRAN DMAP (Direct Matrix Abstraction Program) in an existing NASA Lewis Research Center coupled loads methodology, solving modal equations of motion with initial conditions is possible using either coupled (Newmark-Beta) or uncoupled (exact mode superposition) integration available within module TRD1. Both the coupled and newly developed exact mode superposition methods have been used to perform transient analyses of various space systems. However, experience has shown that in most cases, significant time savings are realized when the equations of motion are integrated using the uncoupled solver instead of the coupled solver. Through the results of a real-world engineering analysis, advantages of using the exact mode superposition methodology are illustrated.

Barnett, Alan R.↗

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗

Integrated Modeling for the Next Generation Space Telescope "Yardstick" Concept

The so-called NASA "Yardstick" design concept for the Next Generation Space Telescope presents unique challenges for systems-level analysis. Simulations that integrate controls, optics, thermal, and structural models are required to evaluate baseline performance, study design sensitivities, and perform design optimization. An integrated modeling approach was chosen using a combination of commercial off-the-shelf and "in-house" developed codes. The resulting capability provides a foundation for linear and non-linear analysis, using both the time and frequency-domain methods. It readily allows various combinations of design parameters and environmental loads to be evaluated directly in terms of key science-related metrics, in this case the scalar RMS (root mean square) line-of-sight and RMS wavefront errors. This presentation first addresses the development of the component, or discipline, models for the Yardstick design. It will then proceed to present the integration of the component models, using linear-systems approaches, in order to support two of the most critical baseline performance analyses: jitter and thermal-elastic stability of the optical telescope assembly (OTA). The results of the jitter analysis indicate that disturbances from the reaction wheels coupled with the lightly-damped and highly-flexible structure present significant challenges to the baseline line-of-sight control architecture. Vibration isolation will be required to meet jitter error requirements. The results of the thermal-elastic analysis indicate that the mirror segment displacements due to ground-to-orbit cool-down of the telescope are within the expected capture range of the segment rigid-body control actuators. This means we will be able to align and phase the primary mirror. However, the results for the analysis of the thermal transient response following an attitude maneuver (slew) show that this telescope design is not sufficiently stable, passively, to meet the wavefront error requirements. Structural re-design is one possibility; alternatively, active thermal control of the OTA may be considered. The Yardstick integrated models were successfully used to demonstrate the feasibility of two thermal control strategies.

Mosier, Gary E.↗

Crater ejecta scaling laws - Fundamental forms based on dimensional analysis

Self-consistent scaling laws are developed for meteoroid impact crater ejecta. Attention is given to the ejection velocity of material as a function of the impact point, the volume of ejecta with a threshold velocity, and the thickness of ejecta deposit in terms of the distance from the impact. Use is made of recently developed equations for energy and momentum coupling in cratering events. Consideration is given to scaling of laboratory trials up to real-world events and formulations are developed for calculating the ejection velocities and ejecta blanket profiles in the gravity and strength regimes of crater formation. It is concluded that, in the gravity regime, the thickness of an ejecta blanket is the same in all directions if the thickness and range are expressed in terms of the crater radius. In the strength regime, however, the ejecta velocities are independent of crater size, thereby allowing for asymmetric ejecta blankets. Controlled experiments are recommended for the gravity/strength transition.

Housen, K. R.↗

Active control of sound transmission/radiation from elastic plates by vibration inputs. I - Analysis

Active control of sound radiation from vibrating plates by oscillating forces applied directly to the structure is analytically studied. The model consists of a plane acoustic wave incident on a clamped elastic circular thin plate. Control is achieved by point forces, and quadratic optimization is used to calculate the optimal control gains necessary to minimize a cost function proportional to the radiated acoustic power (the transmitted field). The results show that global attenuation of broadband radiated sound levels for low to mid-range frequencies can be achieved with one or two control forces, irrespective of whether the system is on or off resonance. The efficiency of the control strategy is demonstrated to be related to the nature of the coupling between the plate modes of response and the radiated field.

Fuller, C. R.↗