Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge instability”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Aero-acoustics of Drag Generating Swirling Exhaust Flows

Aircraft on approach in high-drag and high-lift configuration create unsteady flow structures which inherently generate noise. For devices such as flaps, spoilers and the undercarriage there is a strong correlation between overall noise and drag such that, in the quest for quieter aircraft, one challenge is to generate drag at low noise levels. This paper presents a rigorous aero-acoustic assessment of a novel drag concept. The idea is that a swirling exhaust flow can yield a steady, and thus relatively quiet, streamwise vortex which is supported by a radial pressure gradient responsible for pressure drag. Flows with swirl are naturally limited by instabilities such as vortex breakdown. The paper presents a first aero-acoustic assessment of ram pressure driven swirling exhaust flows and their associated instabilities. The technical approach combines an in-depth aerodynamic analysis, plausibility arguments to qualitatively describe the nature of acoustic sources, and detailed, quantitative acoustic measurements using a medium aperture directional microphone array in combination with a previously established Deconvolution Approach for Mapping of Acoustic Sources (DAMAS). A model scale engine nacelle with stationary swirl vanes was designed and tested in the NASA Langley Quiet Flow Facility at a full-scale approach Mach number of 0.17. The analysis shows that the acoustic signature is comprised of quadrupole-type turbulent mixing noise of the swirling core flow and scattering noise from vane boundary layers and turbulent eddies of the burst vortex structure near sharp edges. The exposed edges are the nacelle and pylon trailing edge and the centerbody supporting the vanes. For the highest stable swirl angle setting a nacelle area based drag coefficient of 0.8 was achieved with a full-scale Overall Sound Pressure Level (OASPL) of about 40dBA at the ICAO approach certification point.

Shah, P. N.↗
NASA NTRS2021

The Asymmetric Inner Disk of the Herbig Ae Star HD 163296 in the Eyes of VLTI/MATISSE: Evidence for a Vortex?

Context.A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observa-tions play a key role in our understanding of the disk structure and the dynamical processes at work.Aims.In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-massstar HD 163296 from early VLTI/MATISSE observations taken in theL- andN-bands. We put special emphasis on the detection ofpotential disk asymmetries.Methods.We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothedring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetriesas well. We also perform numerical hydrodynamical simulations of the inner edge of the disk.Results.Our modeling reveals a significant brightness asymmetry in theL-band disk emission. The brightness maximum of the asym-metry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio inthe azimuthal variation is3.5±0.2. Comparing our result on the location of the asymmetry with other interferometric measurements,we confirm that the morphology of ther<0.3au disk region is time-variable. We propose that this asymmetric structure, located in ornear the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortexis created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-lightradius of theL-band-emitting region is0.33±0.01au, the inclination is52◦+5◦−7◦, and the position angle is143◦±3◦. Our models predictthat a non-negligible fraction of theL-band disk emission originates inside the dust sublimation radius forμm-sized grains. Refractorygrains or large (&10μm-sized) grains could be the origin of this emission.N-band observations may also support a lack of smallsilicate grains in the innermost disk (r.0.6au), in agreement with our findings fromL-band data.

J Varga↗
NASA NTRS1986

Viscous vortex flows

Several computational studies are currently being pursued that focus on various aspects of representing the entire lifetime of the viscous trailing vortex wakes generated by an aircraft. The formulation and subsequent near-wing development of the leading-edge vortices formed by a delta wing are being calculated at modest Reynolds numbers using a three-dimensional, time-dependent Navier-Stokes code. Another computational code was developed to focus on the roll-up, trajectory, and mutual interaction of trailing vortices further downstream from the wing using a two-dimensional, time-dependent, Navier-Stokes algorithm. To investigate the effect of a cross-wind ground shear flow on the drift and decay of the far-field trailing vortices, a code was developed that employs Euler equations along with matched asymptotic solutions for the decaying vortex filaments. And finally, to simulate the conditions far down stream after the onset of the Crow instability in the vortex wake, a full three-dimensional, time-dependent Navier-Stokes code was developed to study the behavior of interacting vortex rings.

Weston, R. P.↗
NASA NTRS1975

On possible interactions between upper and lower atmosphere

Geomagnetic data was compared with data on tropospheric and stratospheric circulation characteristics; a statistically significant shrinking was found in areal extent of the stratospheric vortex from the third to the eighth day following a 'geomagnetic storm' The meridionality of the 30 640-m contour line at 10 millibars increases markedly from 5 to 8 days after the storm. During the contraction of the polar vortex edge, the mean height of the vortex central contour decreases only slightly. This indicates that a stratospheric warming event is associated with a steepening of the contour gradient rather than a warming over the entire area of the stratospheric polar vortex. The troposphere reacts to these weak, but significant, stratospheric warming events by a shrinkage of the area of the 500-millibar cold air pool. This shrinkage commences about 3 days after the stratospheric warming. The investigation indicates that the energy input into the stratosphere that is received in conjunction with the geomagnetic disturbance has to come at a propitious time, when the stratospheric-tropospheric circulation system is not already undergoing a major readjustment because of an inherent dynamic instability.

Macdonald, B. C.↗
NASA NTRS1942

Control-Surface Instability on High-Speed Airplanes

Tests of several modern airplanes indicate that control surfaces with a high degree of aerodynamic balance are likely to possess characteristics which make them unsatisfactory or dangerous in high-speed flight. Dive tests made in the spring of 1940 at the NACA on a naval fighter-type airplane illustrate one form of instability that may be encountered. During a dive at an indicated airspeed of 365 miles per hour, the ailerons suddenly overbalanced. The efforts of the pilot to bring the ailerons back to neutral resulted in a violent oscillation of the control stick from side to side. Fortunately, the force required to return the ailerons to neutral was within the pilot's capabilities. A time history of the maneuver is given in figure1 and typical frames from motion pictures of the cockpit and of the wing, taken during the maneuver, are given in figure 2. In the illustrated case, the occurrence of aerodynamic overbalance was attributed to a slight bulge, approximately 1/16 inch thick, on the lower surface of the leading edges of the ailerons, caused by the installation of additional mass balance ahead of the hinge line. A drawing showing the shape of the bulge is given in figure 3. After this slight protuberance had been eliminated, dives were successfully made at higher speeds.

Phillips, William H.↗
NASA NTRS1992

Discrete-vortex simulation of pulsating flow on a turbulent leading-edge separation bubble

Studies are made of the turbulent separation bubble in a two-dimensional semi-infinite blunt plate aligned to a uniform free stream with a pulsating component. The discrete-vortex method is applied to simulate this flow situation because this approach is effective for representing the unsteady motions of the turbulent shear layer and the effect of viscosity near the solid surface. The numerical simulation provides reasonable predictions when compared with the experimental results. A particular frequency with a minimum reattachment is related to the drag reduction. The most effective frequency is dependent on the amplified shedding frequency. The turbulent flow structure is scrutinized. This includes the time-mean and fluctuations of the velocity and the surface pressure, together with correlations between the fluctuating components. A comparison between the pulsating flow and the non-pulsating flow at the particular frequency of the minimum reattachment length of the separation bubble suggests that the large-scale vortical structure is associated with the shedding frequency and the flow instabilities.

Sung, Hyung Jin↗
NASA NTRS1998

Inner Magnetospheric Superthermal Electron Transport: Photoelectron and Plasma Sheet Electron Sources

Two time-dependent kinetic models of superthermal electron transport are combined to conduct global calculations of the nonthermal electron distribution function throughout the inner magnetosphere. It is shown that the energy range of validity for this combined model extends down to the superthermal-thermal intersection at a few eV, allowing for the calculation of the en- tire distribution function and thus an accurate heating rate to the thermal plasma. Because of the linearity of the formulas, the source terms are separated to calculate the distributions from the various populations, namely photoelectrons (PEs) and plasma sheet electrons (PSEs). These distributions are discussed in detail, examining the processes responsible for their formation in the various regions of the inner magnetosphere. It is shown that convection, corotation, and Coulomb collisions are the dominant processes in the formation of the PE distribution function and that PSEs are dominated by the interplay between the drift terms. Of note is that the PEs propagate around the nightside in a narrow channel at the edge of the plasmasphere as Coulomb collisions reduce the fluxes inside of this and convection compresses the flux tubes inward. These distributions are then recombined to show the development of the total superthermal electron distribution function in the inner magnetosphere and their influence on the thermal plasma. PEs usually dominate the dayside heating, with integral energy fluxes to the ionosphere reaching 10(exp 10) eV/sq cm/s in the plasmasphere, while heating from the PSEs typically does not exceed 10(exp 8) eV/sq cm/s. On the nightside, the inner plasmasphere is usually unheated by superthermal electrons. A feature of these combined spectra is that the distribution often has upward slopes with energy, particularly at the crossover from PE to PSE dominance, indicating that instabilities are possible.

Khazanov, G. V.↗
NASA NTRS1998

Inner Magnetospheric Superthermal Electron Transport: Photoelectron and Plasma Sheet Electron Sources

Two time-dependent kinetic models of superthermal electron transport are combined to conduct global calculations of the nonthermal electron distribution function throughout the inner magnetosphere. It is shown that the energy range of validity for this combined model extends down to the superthermal-thermal intersection at a few eV, allowing for the calculation of the entire distribution function and thus an accurate heating rate to the thermal plasma. Because of the linearity of the formulas, the source terms are separated to calculate the distributions from the various populations, namely photoelectrons (PEs) and plasma sheet electrons (PSEs). These distributions are discussed in detail, examining the processes responsible for their formation in the various regions of the inner magnetosphere. It is shown that convection, corotation, and Coulomb collisions are the dominant processes in the formation of the PE distribution function, and that PSEs are dominated by the interplay between the drift terms. Of note is that the PEs propagate around the nightside in a narrow channel at the edge of the plasmasphere as Coulomb collisions reduce the fluxes inside of this and convection compresses the flux tubes inward. These distributions are then recombined to show the development of the total superthermal electron distribution function in the inner magnetosphere and their influence on the thermal plasma. PEs usually dominate the dayside heating, with integral energy fluxes to the ionosphere reaching 10(exp 10) eV/sq cm/s in the plasmasphere, while heating from the PSEs typically does not exceed 10(exp 8)eV/sq cm/s. On the nightside, the inner plasmasphere is usually unheated by superthermal electrons. A feature of these combined spectra is that the distribution often has upward slopes with energy, particularly at the crossover from PE to PSE dominance, indicating that instabilities are possible.

Khazanov, G. V.↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRS2005

Parametric Inlet Tested in Glenn's 10- by 10-Foot Supersonic Wind Tunnel

The Parametric Inlet is an innovative concept for the inlet of a gas-turbine propulsion system for supersonic aircraft. The concept approaches the performance of past inlet concepts, but with less mechanical complexity, lower weight, and greater aerodynamic stability and safety. Potential applications include supersonic cruise aircraft and missiles. The Parametric Inlet uses tailored surfaces to turn the incoming supersonic flow inward toward an axis of symmetry. The terminal shock spans the opening of the subsonic diffuser leading to the engine. The external cowl area is smaller, which reduces cowl drag. The use of only external supersonic compression avoids inlet unstart--an unsafe shock instability present in previous inlet designs that use internal supersonic compression. This eliminates the need for complex mechanical systems to control unstart, which reduces weight. The conceptual design was conceived by TechLand Research, Inc. (North Olmsted, OH), which received funding through NASA s Small-Business Innovation Research program. The Boeing Company (Seattle, WA) also participated in the conceptual design. The NASA Glenn Research Center became involved starting with the preliminary design of a model for testing in Glenn s 10- by 10-Foot Supersonic Wind Tunnel (10 10 SWT). The inlet was sized for a speed of Mach 2.35 while matching requirements of an existing cold pipe used in previous inlet tests. The parametric aspects of the model included interchangeable components for different cowl lip, throat slot, and sidewall leading-edge shapes and different vortex generator configurations. Glenn researchers used computational fluid dynamics (CFD) tools for three-dimensional, turbulent flow analysis to further refine the aerodynamic design.

Slater, John W.↗
NASA NTRS1992

Midlatitude weather systems on Mars: Is there a hemispheric asymmetry?

The Viking Lander meteorology observations confirmed that midlatitude weather systems are present in the Northern Hemisphere of Mars during the autumn, winter, and early spring. These systems have properties consistent with a development through baroclinic instability of the wintertime zonal-mean circulation. It is known that the weather systems must be of importance for the zonal-mean circulation by virtue of their heat and momentum transports. Observations show that they are associated with dust raising, and they must be involved in the transport of dust; the weather systems almost certainly must act to produce significant water transports as well. Recent simulations with the ARC Mars General Circulation Model (GCM) show that the heat transports by the weather system (and by the zonal-mean circulation) during northern winter can be very substantial: large enough to significantly reduce the CO2 condensation at polar latitudes. This could be largely responsible for the observed early spring halt in the north polar cap recession and for the presence of westerly winds near the cap edge during this season.

Barnes, J. R.↗
NASA NTRS1988

Plasma waves in the range of the lower hybrid frequency - ISEE 1 and 2 observations at the earth's bow shock

This report presents a characterization of plasma wave noise in the range of the lower hybrid frequency associated with 65 crossings of earth's bow shock observed by the ISEE 1 and 2 satellites. Wave growth generally becomes detectable at the upstream edge of the shock foot, increases at the upstream edge of the shock ramp, peaks within the ramp, and then quickly decays to steady downstream values. The upstream extent of the noise is on the same order as that of specularly reflected gyrating ions. Similar profiles were observed in subcritical and supercritical shocks, and no special behavior was associated with the first critical Mach number. Spectra in the foot and ramp were similar in shape, although the noise was 1 to 2 orders of magnitude more intense in the shock ramps than in the feet. Electric field intensities are positively correlated with solar wind speed and inversely related to electron beta and Mach number. Magnetic components are positively correlated with Mach number and beta. The results are generally consistent with suggestions that the noise consists of lower hybrid waves driven by reflected gyrating ions in the foot, and by additional instabilities, such as the cross-field current, in the shock ramp.

Mellott, M. M.↗
NASA NTRS2009

Impulsive Reconnection in the Sun's Atmosphere

Recent high-resolution observations from the Hinode mission show dramatically that the Sun's atmosphere is filled with explosive activity ranging from chromospheric explosions that reach heights of Mm, to coronal jets that can extend to solar radii, to giant coronal mass ejections (CME) that reach the edge of the heliosphere. The driver for all this activity is believed to be 3D magnetic reconnection. From the large variation observed in the temporal behavior of solar activity, it is clear that reconnection in the corona must take on a variety of distinct forms. The explosive nature of jets and CMEs requires that the reconnection be impulsive in that it stays off until a substantial store of free energy has been accumulated, but then turns on abruptly and stays on until much of this free energy is released. The key question, therefore, is what determines whether the reconnection is impulsive or not. We present some of the latest observations and numerical models of explosive and non-explosive solar activity. We argue that, in order for the reconnection to be impulsive, it must be driven by a quasi-ideal instability. We discuss the generality of our results for understanding 31) reconnection in other contexts.

Antiochos, Spiro K.↗
NASA NTRS1997

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗
NASA NTRS1992

Preliminary study of the CRRES magnetospheric barium releases

Preliminary theoretical and computational analyses of the Combined Release and Radiation Effects Satellite (CRRES) magnetospheric barium releases are presented. The focus of the studies is on the evolution of the diamagnetic cavity which is formed by the barium ions as they expand outward, and on the structuring of the density and magnetic field during the expansion phase of the releases. Two sets of simulation studies are discussed. The first set is based upon a 2D ideal MHD code and provides estimates of the time and length scales associated with the formation and collapse of the diamagnetic cavity. The second set uses a nonideal MHD code; specifically, the Hall term is included. This additional term is critical to the dynamics of sub-Alfvenic plasma expansions, such as the CRRES barium releases, because it leads to instability of the expanding plasma. Detailed simulations of the G4 and G10 releases were performed. In both cases the expanding plasma rapidly structured: the G4 release structured at time t less than about 3 s and developed scale sizes of about 1-2 km, while the G10 release structured at time t less than about 22 s and developed scale sizes of about 10-15 km. It is also found that the diamagnetic cavity size is reduced from those obtained from the ideal MHD results because of the structure. On the other hand, the structuring allows the formation of plasma blobs which appear to free stream across the magnetic field; thus, the barium plasma can propagate to larger distances traverse to the magnetic field than the case where no structuring occurs. Finally, a new normal mode of the system was discovered which may be excited at the leading edge of the expanding barium plasma.

Huba, J. D.↗
NASA NTRS2012

Variability of Jupiter's Five-Micron Hot Spot Inventory

Global upheavals on Jupiter involve changes in the albedo of entire axisymmetric regions, lasting several years, with the last two occurring in 1989 and 2006. Against this backdrop of planetary-scale changes, discrete features such as the Great Red Spot (GRS), and other vortices exhibit changes on shorter spatial- and time-scales. We track the variability of the discrete equatorial 5-micron hot spots, semi-evenly spaced in longitude and confined to a narrow latitude band centered at 6.5degN (southern edge of the North Equatorial Belt, NEB), abundant in Voyager images. Tantalizingly similar patterns were observed in the visible (bright plumes and blue-gray regions), where reflectivity in the red is anti-correlated with 5-microns thermal radiance. Ortiz et al. (1998, GRL, 103) characterized the latitude and drift rates of the hot spots, including the descent of the Galileo probe at the southern edge of a 5-micron hot spot, as the superposition of equatorial Rossby waves, with phase speeds between 99 - 103m/s, relative to System III. We note that the high 5-micron radiances correlate well but not perfectly with high 8.57-micron radiances. Because the latter are modulated primarily by changes in the upper ammonia (NH3) ice cloud opacity, this correlation implies that changes in the ammonia ice cloud field may be responsible for the variability seen in the 5-~m maps. During the NEB fade (2011 - early 2012), however, these otherwise ubiquitous features were absent, an atmospheric state not seen in decades. The ongoing NEB revival indicates nascent 5-~m hot spots as early as April 2012, with corresponding visible dark spots. Their continuing growth through July 2012 indicates the possit.le re-establishment of Rossby waves. The South Equatorial Belt (SEB) and NEB revivals began similarly with an instability that developed into a major outbreak, and many similarities in the observed propagation of clear regions.

Yanamandra-Fisher, Padma A.↗
NASA NTRS2016

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗
NASA NTRS2000

Test Cases for Flutter of the Benchmark Models Rectangular Wings on the Pitch and Plunge Apparatus

The supercritical airfoil was chosen as a relatively modem airfoil for comparison. The BOO12 model was tested first. Three different types of flutter instability boundaries were encountered, a classical flutter boundary, a transonic stall flutter boundary at angle of attack, and a plunge instability near M = 0.9 and for zero angle of attack. This test was made in air and was Transonic Dynamics Tunnel (TDT) Test 468. The BSCW model (for Benchmark SuperCritical Wing) was tested next as TDT Test 470. It was tested using both with air and a heavy gas, R-12, as a test medium. The effect of a transition strip on flutter was evaluated in air. The B64AOlO model was subsequently tested as TDT Test 493. Some further analysis of the experimental data for the BOO12 wing is presented. Transonic calculations using the parameters for the BOO12 wing in a two-dimensional typical section flutter analysis are given. These data are supplemented with data from the Benchmark Active Controls Technology model (BACT) given and in the next chapter of this document. The BACT model was of the same planform and airfoil as the BOO12 model, but with spoilers and a trailing edge control. It was tested in the heavy gas R-12, and was instrumented mostly at the 60 per cent span. The flutter data obtained on PAPA and the static aerodynamic test cases from BACT serve as additional data for the BOO12 model. All three types of flutter are included in the BACT Test Cases. In this report several test cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flutter. Cases are selected for classical and stall flutter for the BSCW model, for classical and plunge for the B64AOlO model, and for classical flutter for the BOO12 model. Test Cases are also presented for BSCW for static angles of attack. Only the mean pressures and the real and imaginary parts of the first harmonic of the pressures are included in the data for the test cases, but digitized time histories have been archived. The data for the test cases are available as separate electronic files. An overview of the model and tests is given, the standard formulary for these data is listed, and some sample results are presented.

Robert M Bennett↗