Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bifurcation”

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 415 records · Page 23

Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations: Global Asymptotic Behavior of Time Discretizations - Part 2

The global asymptotic nonlinear behavior of 11 explicit and implicit time discretizations for four 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed. The objectives are to gain a basic understanding of the difference in the dynamics of numerics between the scalars and systems of nonlinear autonomous ODEs and to set a baseline global asymptotic solution behavior of these schemes for practical computations in computational fluid dynamics. We show how 'numerical' basins of attraction can complement the bifurcation diagrams in gaining more detailed global asymptotic behavior of time discretizations for nonlinear differential equations (DEs). We show how in the presence of spurious asymptotes the basins of the true stable steady states can be segmented by the basins of the spurious stable and unstable asymptotes. One major consequence of this phenomenon which is not commonly known is that this spurious behavior can result in a dramatic distortion and, in most cases, a dramatic shrinkage and segmentation of the basin of attraction of the true solution for finite time steps. Such distortion, shrinkage and segmentation of the numerical basins of attraction will occur regardless of the stability of the spurious asymptotes, and will occur for unconditionally stable implicit linear multistep methods. In other words, for the same (common) steady-state solution the associated basin of attraction of the DE might be very different from the discretized counterparts and the numerical basin of attraction can be very different from numerical method to numerical method. The results can be used as an explanation for possible causes of error, and slow convergence and nonconvergence of steady-state numerical solutions when using the time-dependent approach for nonlinear hyperbolic or parabolic PDEs.

Yee, H. C.↗

Comparison of the Aeroacoustics of Two Small-Scale Supersonic Inlets

An aerodynamic and acoustic investigation was performed on two small-scale supersonic inlets to determine which inlet would be more suitable for a High Speed Civil Transport (HSCT) aircraft during approach and takeoff flight conditions. The comparison was made between an axisymmetric supersonic P inlet and a bifurcated two-dimensional supersonic inlet. The 1/14 scale model supersonic inlets were used in conjunction with a 4.1 in (10.4 cm) turbofan engine simulator. A bellmouth was utilized on each inlet to eliminate lip separation commonly associated with airplane engine inlets that are tested under static conditions. Steady state measurements of the aerodynamic flowfield and acoustic farfield were made in order to evaluate the aeroacoustic performance of the inlets. The aerodynamic results show the total pressure recovery of the two inlets to be nearly identical, 99% at the approach condition and 98% at the takeoff condition. At the approach fan speed (60% design speed), there was no appreciable difference in the acoustic performance of either inlet over the entire 0 deg to 110 deg farfield measurement sector. The inlet flow field results at the takeoff fan speed (88% design speed), show the average inlet throat Mach number for the P inlet (Mach 0.52) to be approximately 2 times that of the 2D inlet (Mach 0.26). The difference in the throat Mach number is a result of the smaller throughflow area of the P inlet. This reduced area resulted in a 'soft choking' of the P inlet which lowered the tone and overall sound pressure levels of the simulator in the forward sector by an average of 9 dB and 3 dB, respectively, when compared to the 2D inlet.

Ng, Wing↗

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations that include transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Using several simplifying assumptions, linearized, stability equations are derived that describe the response of the plate just after bifurcation buckling occurs. These equations are then modified to allow the plate reference surface to be located a distance z(c), from the centroid surface which is convenient for modeling stiffened-plate assemblies. The implementation of the new theory into the VICONOPT buckling and vibration analysis and optimum design program code is described. Either classical plate theory (CPT) or first-order shear-deformation plate theory (SDPT) may be selected in VICONOPT. Comparisons of numerical results for several example problems with different loading states are made. Results from the new curved-plate analysis compare well with closed-form solution results and with results from known example problems in the literature. Finally, a design-optimization study of two different cylindrical shells subject to uniform axial compression is presented.

McGowan, David M.↗

Radar Observations of Asteroids 7 Iris, 9 Metis, 12 Victoria, 216 Kleopatra, and 654 Zelinda

We report 13-cm wavelength radar observations of the main-belt asteroids 7 Iris, 9 Metis, 12 Victoria, 216 Kleopatra, and 654 Zelinda obtained at Arecibo between 1980 and 1989. The echoes are highly polarized yet broadly distributed in Doppler frequency, indicating that our targets are smooth on decimeter scales but very rough on some scale(s) larger than about I m. The echo spectra are generally consistent with existing size, shape, and spin information based on radiometric, lightcurve, and occultation data. All of our targets possess distinctive radar signatures that reveal large- scale topography. Reflectivity spikes within narrow ranges of rotation phase suggest large flat regions on Iris, Metis, and Zelinda, while bimodal spectra imply nonconvex, possibly bifurcated shapes for Kleopatra and Victoria. Kleopatra has the highest radar albedo yet measured for a main-belt asteroid, indicating a high metal concentration and making Kleopatra the best main-belt candidate for a core remnant of a differentiated and subsequently disrupted parent body. Upon completion of the Arecibo telescope upgrade, there will be several opportunities per year to resolve main-belt asteroids with hundreds of delay-Doppler cells, which can be inverted to provide estimates of both three-dimensional shape and radar scattering properties.

Mitchell, David L.↗

Radiative Ignition and Transition to Spread Investigation (RITSI)

A total of 25 experiments was conducted in the Glovebox Facility on the Third United States Microgravity payload (USMP3) mission to study the behavior of both flaming and smoldering ignition events, the transition from ignition to the flame/smoldering spread, and the flame/smoldering growth pattern in air. Twenty one of the 25 experiments were flaming experiments and the remainder were smoldering experiments. Ignition was initiated by a heated wire across a thermally thin filter paper in the middle of the sample for the two-dimensional configuration and by a focused beam from a halogen/tungsten lamp at the center of the sample for the three-dimensional configuration. The external air flow velocity was varied from 0 cm/s to 6.5 cm/s. The ignition and subsequent flame spread events were recorded by a video camera, a 35 mm camera, and 6 thermocouples (2 in the gas phase and 4 in the sample). The results indicate that non-piloted radiative ignition of the paper by external thermal radiation tends to occur more easily than in normal gravity. In the two-dimensional configuration, the transition from ignition to downstream flame spread never occurred; only the transition to upstream spread took place. In the three-dimensional configuration, a char growth pattern emerged in the shape of a fan, with the ignited area as the center of the fan and the fan spreading in the upstream direction. The internal angle of the fan increased with an increase in external flow velocity (higher flow velocity gave a more opened fan pattern). At an external flow velocity of 5 cm/s, the flame was horseshoe-shaped and the char pattern became an elongated semicircle toward upstream. It appears that the peak upstream flame spread rate was reached at a higher external velocity than for the two-dimensional flame front. Flame spread much faster along open edges of the sample paper than along the sample face, presumably due to there being larger supply of oxygen and energy feedback at the open edges than at the surface. These results indicate the possible importance of the effect of the flame front shape such as the curvature of the flame front on spread rate and flame strength. A complex, unexpected finger-shaped char growth pattern was observed during the smoldering experiment. Higher external flow velocity increased the number of localized smoldering fronts, of the char fingers they left behind, and frequency of bifurcations from the fingers. At present, it is not clear what caused such a complex char growth pattern.

Kashiwagi, Takashi↗

The Solar-Stellar Connection (NAG5-6124: SOHO Guest Investigator Program)

The following is a final report from the SOHO Guest Investigator program to use the SUMER far-UV spectrometer to obtain imaging spectroscopy in support of the goals of the so-called "solar-stellar connection." In particular, a major emphasis was utilization of the long-slit time-resolved maps of the solar surface in bright far-UV emission lines to deduce how particular aspects of the temporally and spatially averaged line profiles trace back to individual structural features of the magnetically disturbed outer atmosphere; to help interpret the unresolved line profiles from high quality stellar observations (say, with the Space Telescope Imaging Spectrograph). The researchers served two tours of duty in the SOHO Operations Center as SUMER planners, during which time we conducted an extensive series of observing programs. These can be divided into three general categories: surface mapping, translimb spectroscopy, and active region diagnostics. We have analyzed some of the large volumes of data to the point where we have presented them in poster papers, and in invited papers at national and international meetings. Listed below are the titles of the preliminary publications we have written, including brief abstracts to indicate the main results. (1) Chromospheric structure and Dynamics-- Observations -- The chromosphere is a highly structured dynamic 'layer' of the solar outer atmosphere. Here, not only are the effects of mechanical heating first evident (moving upward in altitude from the deep photosphere), but also the amount of nonradiative energy deposited is far greater than in the albeit much hotter overlying transition region and corona. Further, the chromosphere is by far the thickest zone of the solar atmosphere with respect to the pressure scale height. A major goal of stellar astrophysics is to understand how the chromosphere is heated and why it adopts its peculiar structure. A cursory examination of solar filtergrams and high-resolution movies demonstrates that much of the chromospheric "action" must be occurring on fine spatial scales and short times; particularly in the cell interior transient brightenings, but also in the longer-lived network fragments. That regime of investigation is far removed from what one usually associates with "synoptic" measurements. Nevertheless, synoptic observations of chromospheric indices, filtergrams, and globally-averaged profile parameters not only can provide important insight concerning the crucial role of the cycle-variable part of the solar magnetic field; but they also can forge a key link with analogous measurements of the stars, where often the phenomena can be significantly exaggerated from the solar case, but high spatial resolution reconnaissance is not even a remote possibility. In addition to discussing the synoptic aspects of chromospheric structure and dynamics, I summarized new insights into the general problem of the solar chromosphere that have been obtained recently with the SUMER far-ultraviolet spectrometer on SOHO. (2) CO and the Temperature Structure of the Solar Atmosphere -- The surface layers of the Sun provide a crucial boundary condition for many of the processes that occur in the deep interior. The stratification of the outer solar atmosphere once was thought to be well understood. However, studies of thermally sensitive molecular absorptions in the infrared revealed puzzling anomalies. Strong lines of the CO fundamental vibration-rotation bands near 5 microns showed very cool temperatures at the extreme limb, and remarkable off-limb emissions extending well into the supposedly hot chromosphere. The conflicting pictures of the photosphere/chromosphere interface, from the widely separated wavelength regimes, has raised suspicions that those "layers" of the atmosphere are much more inhomogeneous than previously suspected. One proposal is that the low chromosphere is dominated by cool gas, the "COmosphere," which is threaded by a network of persistent small-scale hot magnetic filaments and occasionally disrupted by localized acoustic disturbances. The COmosphere is capped by the merged fields of the network elements in the chromospheric "canopy." (3) Translimb Spectroscopy with SOHO/SUMER -- We have used SUMER to obtain deep exposures of the 1300-1400 A spectrum, at the extreme limb and off-limb. Previous "translimb" studies in the thermal infrared had revealed remarkable extensions of cold material (T is approximately equal to 3000 K)-traced by carbon monoxide emission lines-into the heart of the hot chromosphere. A main objective of our program was to search for corresponding far-UV signatures of the "thermally-bifurcated" low chromosphere; for example, radiatively fluoresced emissions of the CO A-X 4th-positive system (collisional excitation would be negligible in cold gas). We conducted two separate observing programs with SUMER. Both made use of the 1 inch-diameter circular aperture, translated across the limb in the minimum motor step increments of 0."375, along the central meridian in the Northern polar coronal hole. The first program executed for nine hours beginning 19UT 25 Oct 1996. The full wavelength range was 1340-1400 A. It was recorded in two overlapping segments, placing key regions of the spectrum alternately on the KBr and bare parts of the detector, to help isolate 2nd-order features. Each segment was integrated for 500 s, and 32 pairs were obtained to span a 12" swath centered on the optical limb. The second program was conducted 00-09UT 01 Dec 1996. It consisted of a single wavelength setting (1300-1340 A) with exposure time 500 s, but twice the spatial coverage of the earlier series: 64 steps, for a total displacement of 24". The strong chromospheric resonance lines of atomic oxygen (1302-1305 A) and ionized carbon (1334-1335 A) were observed on the bare part of the MCP camera.

Ayres, Thomas R.↗

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

Constitutive and Stability Behavior of Granular Materials Under Very Low Effective Stresses

The study of the behavior of granular materials such as sand is of great importance to the understanding of the response of foundations resting on surface and near-surface soil deposits. A clear understanding of the behavior of such materials can also provide insights into other related issues where low effective stresses are encountered such as liquefaction. The main sources of the constitutive and stability properties of cohesionless granular materials is interparticle friction, which in turn under low confinement stress levels is strongly affected by gravitational body forces under terrestrial (1 gravity) conditions. Under moderate-to-high stress levels, the influence of gravity on the behavior of experiments may not be pronounced and therefore the test results in a terrestrial environment may be acceptable for engineering purposes. However, the conduct of experiments on granular materials under very low stress levels and under quasi-static conditions can only be performed in a microgravity environment. A series of displacement-controlled cyclic triaxial compression experiments were performed in a SPACEHAB module on the Space Shuttle during the STS-79 mission to Mir in September, 1996, and the STS-89 mission in January, 1998. The experiments were conducted on six right cylindrical specimens 75 mm in diameter and 150 mm long at effective confining pressures of 0.05, 0.52 and 1.30 kPa. The results show very high peak strength friction angles in the range of 47.6 to 70.0 degrees, which are mainly due to overconsolidation and grain interlocking effects. It was observed that the residual strength levels were in the same range as that observed at higher confining stress levels. The dilatancy angles were unusually high in the range of 30 to 31 degrees. All specimens display substantial initial stiffnesses and elastic moduli during unloading and reloading events, which are nearly an order of magnitude higher than conventional theories predict. A periodic instability phenomenon which appears to result from buckling of multiple internal arches and columnar systems, augmented by stick-slips was observed in the experiments. Computed Tomography (CT) measurements revealed valuable data about the internal fabric and the specimens deformation patterns. Uniform diffuse bifurcation with multiple radial shear bands was observed in the specimens tested in a microgravity environment. In the axial direction, two major conical surfaces were developed. Spatial nonsymmetrical deformations were observed in specimens tested in terrestrial laboratory.

Alshibli, Khalid A.↗

Ionospheric Instabilities and Shock Capturing Algorithms

We describe the basic physics of the gradient-drift instability and the collisional Rayleigh-Taylor instability, the two plasma instabilities which dominate the production of small-scale structure in the earth's ionosphere. We then focus our attention on the nonlinear evolution of these instabilities for two particular cases: 1) the recursive bifurcation and "freezing" of artificial ion clouds released in the ionosphere; and 2) equatorial spread F. Finally, we discuss the numerical challenges faced in attempting to numerically simulate these instabilities, and the role played by shock-capturing algorithms therein.

Zalesak, Steven T.↗

Statistical Approach for Determining the Onsets/Durations of ENSO Cycle Extremes

During the interval of 1950-mid 1998, some 16 El Nino and 10 La Nina have been identified on the basis of sea surface temperature in the Nino 3.4 region, these 26 events representing the extremes of the quasi-periodic ENSO cycle. Statistical aspects of these events are examined. Surprisingly, the durations of El Nino and La Nina appear to be strongly bifurcated into shorter and longer duration classes, as do the recurrence periods of El Nino. Moreover, the duration of an El Nino appears to provide a statistically meaningful indication as to when to expect the next onset of El Nino. Because the last El Nino had its onset in April 1997 and was of longer duration, onset of the next El Nino, probably, will not occur until after February 2000.

Wilson, Robert M.↗

Mechanics of Granular Materials (MGM) Microgravity Experiment

The second series of MGM experiment was conducted during the STS-89 mission in January 1998. The experiment was previously flow on Atlantis's STS-79 mission in September 1996. Six displacement-controlled, drained triaxial compression experiments were performed at very low effective confining stresses. The confining stresses were in the ranges 0.05, 0.52 and 1.30 kPa. Three experiments were subjected to monotonic loading and unloading cycles while the other three experiments were subjected to cyclic loading. The results show very high peak strength friction angles in the range of 47.6 to 70.0 degrees, which are mainly due to overconsolidation and grain interlocking effects. It was observed that the residual strength levels in the monotonic loading experiments were in the same range as that observed at higher confining stress levels. The dilatancy angles were unusually high in the range of 30 to 31 degrees. All specimens display substantial initial stiffnesses and elastic moduli during unloading and reloading events, which are nearly an order of magnitude higher than conventional theories predict. A periodic instability phenomenon which appears to result from buckling of multiple internal arches and columnar systems, augmented by stick-slips was observed in the experiments. Computed Tomography (CT) measurements revealed valuable data about the internal fabric and the specimens deformation patterns. Uniform diffuse bifurcation with multiple radial shear bands was observed in the specimens tested in a microgravity environment. In the axial direction, two major conical surfaces were developed. Spatial nonsymmetrical deformations were observed in specimens tested in terrestrial laboratory.

Alshibli, Khalid A.↗

Translimb Spectroscopy with SOHO/SUMER

We have used SUMER to obtain deep exposures of the 1300-1400 A spectrum, at the extreme limb and off-limb. Previous "translimb" studies in the thermal infrared had revealed remarkable extensions of cold material (T approximately 3000 K)--traced by carbon monoxide emission lines--into the heart of the hot chromosphere. A main objective of our program was to search for corresponding far-UV signatures of the "thermally-bifurcated" low chromosphere; for example, radiatively fluoresced emissions of the CO A-X 4th-positive system (collisional excitation would be negligible in cold gas). We conducted two separate observing programs with SUMER. Both made use of the 1"-diameter circular aperture, translated across the limb in the minimum motor step increments of 0.375", along the central meridian in the Northern polar coronal hole. The first program executed for nine hours beginning 19UT 25 Oct 1996. The full wavelength range was 1340-1400 A. It was recorded in two overlapping segments, placing key regions of the spectrum alternately on the KBr and bare parts of the detector, to help isolate 2nd-order features. Each segment was integrated for 500 s, and 32 pairs were obtained to span a 12" swath centered on the optical limb. The second program was conducted 00-09UT 01 Dec 1996. It consisted of a single wavelength setting (1300-1340 A) with exposure time 500 s, but twice the spatial coverage of the earlier series: 64 steps, for a total displacement of 24". The strong chromospheric resonance lines of atomic oxygen (1302-1305 A) and ionized carbon (1334-1335 A) were observed on the bare part of the MCP camera. We report our progress in cataloging the rich, diverse translimb emission spectrum; and our efforts to deduce fundamental properties of the thermally heterogeneous chromosphere.

Ayer, Thomas R.↗

Radar Images of Asteroid 4179 Toutatis

Delay-Doppler images of the Earth-crossing asteroid 4179 Toutatis achieve resolutions as fine as 125 nanoseconds (19 meters in range) and 8.3 millihertz (0.15 millimeter per second in radial velocity) and place hundreds to thousands of pixels on the asteroid, which appears to be several kilometers long, topographically bifurcated, and heavily cratered. The image sequence reveals Toutatis to be in an extremely slow, non-principal axis rotation state.

Ostro, Steven J.↗

Isolated Nozzle Analysis Using OVERFLOW

LaRC conducted a code validation study for the OVERFLOW code to ascertain its accuracy for boattail drag prediction. The OVERFLOW results compared favorably with the LaRC 16-ft. Transonic Wind Tunnel (TWT) data, and prior CFD solutions from PAB3D and CFL3D. The ultimate goal is to investigate the installation drag of the nacelle boattails with powered nozzles at transonic mach numbers. The OVERFLOW solver was chosen because of its ability to accept volume overlapping structured grid for very complex airframe configurations. Structured grid components for representing the transonic nozzle boattail can be added to the BCAG grid for a TCA airframe with 2D bifurcated inlet and flow through nacelle without alteration. The focus of this research was to determine the suitability of the OVERFLOW solver for accomplishing this ultimate goal. This presentation will first introduce the transonic nozzle boattail wind-tunnel model geometry, followed by an examination of aerodynamic features based on the current OVERFLOW solutions and the solutions obtained previously using PAB3D, comparisons of Cp on the flap surface between the OVERFLOW solutions, wind tunnel data, and solutions from other CFD codes, an assessment of boattail drag count prediction, and a work plan for FY99.

Wong, Tin-Chee↗

Electromagnetics and Fluid Dynamics

Previous efforts focused on developing tools for design of low observables were sustained. The final product was the maturation of a high-order accurate finite-volume based code to solve Maxwell's equations. One of the primary achievements was the development and implementation of efficient filtering techniques which enhance the robustness of high-order and optimized schemes without significant adverse impact on accuracy. This has eliminated the stability barrier which restrains the common use of high-order schemes for conservative wave propagation phenomena on curvilinear meshes. A study was performed of crossing shock interactions under conditions of increasing interaction strength and asymmetry. In the first category, the observed computed topological bifurcations were correlated with the formation of various lines of coalescence and divergence evident in experimental and computed surf-ace oil maps. ne flow structure arising from asymmetric interactions was investigated with particular emphasis on: 1) vorticity dynamics, 2) shock-structure and 3) sidewall vortex loading. Several efforts of the prior year were successfully published in archival journals. The high-order algorithms developed for CEM have been implemented into the FDL3DI CFD code are presently undergoing extensive testing. Preliminary results are highly encouraging.

Gaitonde, Datta↗

The Effect of a Rotating Magnetic Field on Rayleigh-Benard Convection

A destabilizing vertical temperature gradient and a rotating magnetic field have been applied to a cylindrical column of liquid gallium. The convective flows which arise as a function of these driving parameters are identified. For small magnetic field strengths, a region of stationary flow is observed. This region is bounded by critical values of the Rayleigh and Hartmann numbers. As the rotating magnetic field is increased, the critical Rayleigh number can increase by more than a factor of 10. The rotating magnetic field itself induces an instability at a critical value of the Hartmann number independent of the Rayleigh number. The nature of the bifurcations (whether subcritical or supercritical) and the convective flows occurring at the critical Rayleigh numbers are dependent upon the Hartmann number. In the stationary flow region, the rotating magnetic field can induce fluid motion with velocities on the order of cm/s. Thus, a rotating magnetic field can be used to achieve the benefits of forced convection during crystal growth processes without triggering deleterious instabilities.

Volz, Martin P.↗

Preliminary Chaotic Model of Snapover on High Voltage Solar Cells

High voltage power systems in space will interact with the space plasma in a variety of ways. One of these, Snapover, is characterized by a sudden enlargement of the electron current collection area across normally insulating surfaces. A power drain on solar array power systems will results from this enhanced current collection. Optical observations of the snapover phenomena in the laboratory indicates a functional relation between bia potential and surface glow area. This paper shall explore the potential benefits of modeling the relation between current and bia potential as an aspect of bifurcation analysis in chaos theory. Successful characterizations of snapover as a chaotic phenomena may provide a means of snapover prevention and control through chaotic synchronization.

Mackey, Willie R.↗