SEARCH · Search NASA

Results for “dynamical phase transitions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Synthesis and Characterization of Ru(II) Tris(1,1O-phenanthroline)-Electron Acceptor Dyads Incorporating the 4-benzoyl-N-methylpyridinium Cation or N-Benzyl-N'-methyl-viologen. Improving the Dynamic Range, Sensitivity and Response Time of Sol-Gel Based Optical Oxygen Sensors

The title compounds (1 and 2, above) were synthesized by Sonogashira coupling reactions of appropriate Ru(1I) complexes with the electron a cceptors. Characterization was conducted in solution and in frozen ma trices. Finally, the title compounds were evaluated as dopants of sol-gel materials. It was found that the intramolecular quenching efficie ncy of 4-benzoyl-Nmethylpyridinium cation in solution depends on the solvent: photoluminescence is quenched completely in CH,CN, but not i n methanol or ethanol. On the other hand, intramolecular emission que nching by 4-benzyl-N-methyl viologen is complete in all solvents. The difference between the two quenchers is traced electrochemically to t he solvation of the 4-benzoyl-Nmethylpyridiniums by alcohol. In froze n matrices or adsorbed on the surfaces of silica aerogel, both Ru(I1) complex/electron acceptor dyads of this study are photoluminescent, and the absence of quenching has been traced to the environmental rigi dity. When doped aerogels are cooled at 77 K, the emission intensity increases by approximately 4x, and the spectra shift to the blue, analogous to what is observed with Ru(I1) complexes in solutions undergoi ng fluid-to-rigid transition. However, in contrast to frozen solution s, the luminescent moieties in the bulk of aerogels kept at low tempe ratures are still accessible to gas-phase quenchers diffusing through the mesopores, leading to more sensitive platforms for sensors than o ther room-temperature configurations. Thus the photoluminescence of o ur Ru(I1) complex dyads adsorbed on aerogel is quenchable by O2 both at room temperature and at 77 K. Furthermore, it was also found that O 2 modulates the photoluminescence of aerogels doped with 4-benzoyl -N -methylpyridinium-based dyads over a wider dynamic range compared wi th aerogels doped with either our vislogen-based dyads or with Ru(I1) tris(1,lO-phenanthroline) itself.

Leventis, Nicholas↗

NASA NTRS1996

Studies of Two-Phase Flow Dynamics and Heat Transfer at Reduced Gravity Conditions

The ability to predict gas-liquid flow patterns is crucial to the design and operation of two-phase flow systems in the microgravity environment. Flow pattern maps have been developed in this study which show the occurrence of flow patterns as a function of gas and liquid superficial velocities as well as tube diameter, liquid viscosity and surface tension. The results have demonstrated that the location of the bubble-slug transition is affected by the tube diameter for air-water systems and by surface tension, suggesting that turbulence-induced bubble fluctuations and coalescence mechanisms play a role in this transition. The location of the slug-annular transition on the flow pattern maps is largely unaffected by tube diameter, liquid viscosity or surface tension in the ranges tested. Void fraction-based transition criteria were developed which separate the flow patterns on the flow pattern maps with reasonable accuracy. Weber number transition criteria also show promise but further work is needed to improve these models. For annular gas-liquid flows of air-water and air- 50 percent glycerine under reduced gravity conditions, the pressure gradient agrees fairly well with a version of the Lockhart-Martinelli correlation but the measured film thickness deviates from published correlations at lower Reynolds numbers. Nusselt numbers, based on a film thickness obtained from standard normal-gravity correlations, follow the relation, Nu = A Re(sup n) Pr(exp l/3), but more experimental data in a reduced gravity environment are needed to increase the confidence in the estimated constants, A and n. In the slug flow regime, experimental pressure gradient does not correlate well with either the Lockhart-Martinelli or a homogeneous formulation, but does correlate nicely with a formulation based on a two-phase Reynolds number. Comparison with ground-based correlations implies that the heat transfer coefficients are lower at reduced gravity than at normal gravity under the same flow conditions. Nusselt numbers can be correlated in a fashion similar to Chu and Jones.

Witte, Larry C.↗

NASA NTRS2001

Experimental Evaluation of the Effect of Angle-of-attack on the External Aerodynamics and Mass Capture of a Symmetric Three-engine Air-breathing Launch Vehicle Configuration at Supersonic Speeds

A subscale aerodynamic model of the GTX air-breathing launch vehicle was tested at NASA Glenn Research Center's 10- by 10-Foot Supersonic Wind Tunnel from Mach 2.0 to 3.5 at various angles-of-attack. The objective of the test was to investigate the effect of angle-of-attack on inlet mass capture, inlet diverter effectiveness, and the flowfield at the cowl lip plane. The flow-through inlets were tested with and without boundary-layer diverters. Quantitative measurements such as inlet mass flow rates and pitot-pressure distributions in the cowl lip plane are presented. At a 3deg angle-of-attack, the flow rates for the top and side inlets were within 8 percent of the zero angle-of-attack value, and little distortion was evident at the cowl lip plane. Surface oil flow patterns showing the shock/boundary-layer interaction caused by the inlet spikes are shown. In addition to inlet data, vehicle forebody static pressure distributions, boundary-layer profiles, and temperature-sensitive paint images to evaluate the boundary-layer transition are presented. Three-dimensional parabolized Navier-Stokes computational fluid dynamics calculations of the forebody flowfield are presented and show good agreement with the experimental static pressure distributions and boundary-layer profiles. With the boundary-layer diverters installed, no adverse aerodynamic phenomena were found that would prevent the inlets from operating at the required angles-of-attack. We recommend that phase 2 of the test program be initiated, where inlet contraction ratio and diverter geometry variations will be tested.

Kim, Hyun D.↗

NASA NTRS2014

Physical and Dynamical Linkages Between Lightning Jumps and Storm Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. Currently, the lightning jump algorithm is being tested in two separate but important efforts. Schultz et al. (2014; this conference) is exploring the transition of the algorithm from its research based formulation to a fully objective algorithm that includes storm tracking, Geostationary Lightning Mapper (GLM) Proxy data and the lightning jump algorithm. Chronis et al. (2014) provides context for the transition to current operational forecasting using lightning mapping array based products. However, what remains is an end-to-end physical and dynamical basis for coupling total lightning flash rates to severe storm manifestation, so the forecaster has a reason beyond simple correlation to utilize the lightning jump algorithm within their severe storm conceptual models. Therefore, the physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relationship to updraft strength, updraft volume, precipitation-sized ice mass, etc.; however, their relationship specifically to lightning jumps is fragmented within the literature. Thus the goal of this study is to use multiple Doppler and polarimetric radar techniques to resolve the physical and dynamical storm characteristics specifically around the time of the lightning jump. This information will help forecasters anticipate lightning jump occurrence, or even be of use to determine future characteristics of a given storm (e.g., development of a mesocyclone, downdraft, or hail signature on radar), providing additional lead time/confidence in the severe storm warning paradigm.

Schultz, Christopher J.↗

NASA NTRSDate not supplied

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown↗

NASA NTRS2014

Physical and Dynamical Linkages between Lightning Jumps and Storm Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. Currently, the lightning jump algorithm is being tested in two separate but important efforts. Schultz et al. (2014; this conference) is exploring the transition of the algorithm from its research based formulation to a fully objective algorithm that includes storm tracking, Geostationary Lightning Mapper (GLM) Proxy data and the lightning jump algorithm. Chronis et al. (2014; this conference) provides context for the transition to current operational forecasting using lightning mapping array based products. However, what remains is an end-to-end physical and dynamical basis for coupling total lightning flash rates to severe storm manifestation, so the forecaster has a reason beyond simple correlation to utilize the lightning jump algorithm within their severe storm conceptual models. Therefore, the physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relationship to updraft strength, updraft volume, precipitation-sized ice mass, etc.; however, their relationship specifically to lightning jumps is fragmented within the literature. Thus the goal of this study is to use multiple Doppler and polarimetric radar techniques to resolve the physical and dynamical storm characteristics specifically around the time of the lightning jump. This information will help forecasters anticipate lightning jump occurrence, or even be of use to determine future characteristics of a given storm (e.g., development of a mesocyclone, downdraft, or hail signature on radar), providing additional lead time/confidence in the severe storm warning paradigm.

Schultz, Christopher J.↗

NASA NTRS2000

The Strategy for the Second Phase of Aerobraking Mars Global Surveyor

On February 19, 1999, the Mars Global Surveyor (MGS) spacecraft was able to propulsively establish its mapping orbit. This event followed the completion of the second phase of aerobraking for the MGS spacecraft on February 4, 1999. For the first time, a spacecraft at Mars had successfully employed aerobraking methods in order to reach its desired pre-launch mapping orbit. This was accomplished despite a damaged spacecraft solar array. The MGS spacecraft was launched on November 7, 1996, and after a ten month interplanetary transit was inserted into a highly elliptical capture orbit at Mars on September 12, 1997. Unlike other interplanetary missions, the MGS spacecraft was launched with a planned mission delta-V ((Delta)V) deficit of nearly 1250 m/s. To overcome this AV deficit, aerobraking techniques were employed. However, damage discovered to one of the spacecraft's two solar arrays after launch forced major revisions to the original aerobraking planning of the MGS mission. In order to avoid a complete structural failure of the array, peak dynamic pressure levels for the spacecraft were established at a major spacecraft health review in November 1997. These peak dynamic pressure levels were roughly one-third of the original mission design values. Incorporating the new dynamic pressure limitations into mission replanning efforts resulted in an 'extended' orbit insertion phase for the mission. This 'extended' orbit insertion phase was characterized by two distinct periods of aerobraking separated by an aerobraking hiatus that would last for several months in an intermediate orbit called the "Science Phasing Orbit" (SPO). This paper describes and focuses on the strategy for the second phase of aerobraking for the MGS mission called "Aerobraking Phase 2." This description will include the baseline aerobraking flight profile, the trajectory control methodology, as well as the key trajectory metrics that were monitored in order to successfully "guide' the spacecraft to its desired mapping orbit. Additionally, the actual aerobraking progress is contrasted to the planned aerobraking flight profile. (A separate paper will describe the navigation aspects of MGS aerobraking in detail.) Key to the success of the MGS mission is the delivery of the spacecraft to its final mapping orbit and the synergy the instrument complement provides to its scientific investigators when science data is returned from that orbit. The MGS mapping orbit is characterized as a low altitude, near-circular, near-polar orbit that is Sun-synchronous with the descending equatorial crossing at 2:00 AM local mean solar time (LMST).

Johnston, M. D.↗

NASA NTRS1989

The influence of velocity-changing collisions on resonant degenerate four-wave mixing

The phase-conjugate signal observed in resonant degenerate four-wave mixing on the 6 3P2 to 7 3S1 transition of atomic Hg in an Hg-Ar discharge is investigated. At a fixed Ar pressure the variation of the signal with pump powers is explained by a model that includes the effects of velocity-changing collisions (VCCs). As the Ar pressure was varied from 0 to 1 torr, an increase in the phase-conjugate signal was observed and is ascribed to a change in the discharge dynamics with Ar pressure and to the influence of VCCs. To further clarify the role of collisions and optical pumping, degenerate four-wave mixing spectra are examined as a function of pump power. Line shapes are briefly discussed.

Richardson, W. H.↗

NASA NTRS2021

Combined Cycle Engine Large-Scale Inlet for Mode Transition Studies: System Identification Rack Software Design

This report describes the development of the system identification (SysID) rack custom application code. This code is used to automate the system identification experimental processes for the National Aeronautics and Space Administration (NASA) Combined Cycle Engine Large-Scale Inlet for Mode Transition Experiments testbed. This series of experiments took place in the NASA Glenn Research Center (GRC) 10 foot by 10 foot (10x10) Supersonic Wind Tunnel (SWT) test facility. The SysID code was developed to apply command signals to testbed actuators, receive actuator feedback signals, and acquire sensor data which is used to identify the dynamics of the experimental processes. This code development and the process identification experiments were conducted during Phase-2 testing at various operating points. These points were predefined from data collected during the Phase-1 testbed characterization experiments. The developed SysID code was executed on a NASA designed and built real-time data acquisition (DAQ) and controls (SysID) rack. This report focuses on the development of the SysID rack real-time DAQ and control hardware and the custom application code used to automate the system identification process. The implementation of this reliable and efficient data acquisition and controls development tool is demonstrated by presenting an example of experimental data collected to be used for system identification.

R Thomas↗

NASA NTRS2000

Modal Analysis and Model Correlation of the Mir Space Station

This paper will discuss on-orbit dynamic tests, modal analysis, and model refinement studies performed as part of the Mir Structural Dynamics Experiment (MiSDE). Mir is the Russian permanently manned Space Station whose construction first started in 1986. The MiSDE was sponsored by the NASA International Space Station (ISS) Phase 1 Office and was part of the Shuttle-Mir Risk Mitigation Experiment (RME). One of the main objectives for MiSDE is to demonstrate the feasibility of performing on-orbit modal testing on large space structures to extract modal parameters that will be used to correlate mathematical models. The experiment was performed over a one-year span on the Mir-alone and Mir with a Shuttle docked. A total of 45 test sessions were performed including: Shuttle and Mir thruster firings, Shuttle-Mir and Progress-Mir dockings, crew exercise and pushoffs, and ambient noise during night-to-day and day-to-night orbital transitions. Test data were recorded with a variety of existing and new instrumentation systems that included: the MiSDE Mir Auxiliary Sensor Unit (MASU), the Space Acceleration Measurement System (SAMS), the Russian Mir Structural Dynamic Measurement System (SDMS), the Mir and Shuttle Inertial Measurement Units (IMUs), and the Shuttle payload bay video cameras. Modal analysis was performed on the collected test data to extract modal parameters, i.e. frequencies, damping factors, and mode shapes. A special time-domain modal identification procedure was used on free-decay structural responses. The results from this study show that modal testing and analysis of large space structures is feasible within operational constraints. Model refinements were performed on both the Mir alone and the Shuttle-Mir mated configurations. The design sensitivity approach was used for refinement, which adjusts structural properties in order to match analytical and test modal parameters. To verify the refinement results, the analytical responses calculated using original and refined math models were compared with the measured responses. The results from this study show that the refined models predict more accurately the dynamics of the actual structure. The MiSDE test and analysis results provided the information and experience on test design, flight testing, and data analysis for the ISS and other future spacecraft, which are critical to the verification of analytical models and structural loads.

Kim, Hyoung M.↗

NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗

NASA NTRS1986

Dynamic dielectric analysis for nondestructive cure monitoring and process control

Dynamic dielectric analysis (DDA) is an effective in situ NDE method that can monitor the reaction status in thermosets and the phase changes in thermoplastics, including slow reactions occuring late in the cure cycle and recrystallization during annealing. The effects of moisture and resin history on reaction rate can also be determined, as can ionic and dipolar contributions. The ionic mobility parameter is noted to be an excellent monitor of viscosity above the glass transition temperature. The ability of DDA to monitor cure rate variations in a thick section during autoclaving has been demonstrated.

Kranbuehl, D. E.↗

NASA NTRS1993

Dynamical instability as the cause of the massive outbursts in Eta Carinae and other luminous blue variables

A new type of stellar envelope structure has been computationally discovered at very high stellar masses. The outer part of the envelope resembles a nearly detached, diffusely filled shell overlying an ultrahot surface of small radius. This structural anomaly is caused by a large iron bump occurring in the new opacities of Iglesias et al. (1992). The new stellar models with normal metallicity encounter a strong ionization-induced dynamical instability in the outer envelope as they rapidly transit the H-R diagram after the end of central hydrogen burning. Preliminary evolutionary and hydrodynamical calculations successfully mimic the most basic observed properties of Eta Carinae and other very luminous blue variables. The Humphreys-Davidson sloped line in the H-R diagram, however, seems to be unrelated to these variables, and is instead the observed terminus of the main-sequence phase of evolution if convective core overshooting is insignificant.

Stothers, Richard B.↗

NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

NASA NTRS2016

Experimental Study of Tip Vortex Flow from a Periodically Pitched Airfoil Section

An experimental investigation of a tip vortex from a NACA0012 airfoil is conducted in a low-speed wind tunnel at a chord Reynolds number of 4x10(exp 4). Initially, data for a stationary airfoil held at various angles-of-attack (alpha) are gathered. Detailed surveys are done for two cases: alpha=10 deg with attached flow and alpha=25 deg with massive flow separation on the upper surface. Distributions of various properties are obtained using hot-wire anemometry. Data include mean velocity, streamwise vorticity and turbulent stresses at various streamwise locations. For all cases, the vortex core is seen to involve a mean velocity deficit. The deficit apparently traces to the airfoil wake, part of which gets wrapped by the tip vortex. At small alpha, the vortex is laminar within the measurement domain. The strength of the vortex increases with increasing alpha but undergoes a sudden drop around alpha (is) greater than 16 deg. The drop in peak vorticity level is accompanied by transition and a sharp rise in turbulence within the core. Data are also acquired with the airfoil pitched sinusoidally. All oscillation cases pertain to a mean alpha=15 deg while the amplitude and frequency are varied. An example of phase-averaged data for an amplitude of +/-10 deg and a reduced frequency of k=0.2 is discussed. All results are compared with available data from the literature shedding further light on the complex dynamics of the tip vortex.

vortex flows↗

NASA NTRS2022

Planet Hunters TESS IV: A Massive, Compact Hierarchical Triple Star System TIC 470710327

We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA’s Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse-timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14–17 M ꙩ ) is more massive than the combined mass of the inner binary (10.9–13.2 M ꙩ ). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne–Żytkow object. Further observational characterization of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products.

NASA's Transiting Exoplanet Survey Satellite (TESS↗

NASA NTRS1989

Lessons learned in the transition to Ada from FORTRAN at NASA/Goddard

Two dynamics satellite simulators are developed from the same requirements, one in Ada and the other in FORTRAN. The purpose of the research was to find out how well the prescriptive Ada development model worked to develop the Ada simulator. The FORTRAN simulator development, as well as past FORTRAN developments, provided a baseline for comparison. Since this was the first simulator developed, the prescriptive Ada development model had many similarities to the usual FORTRAN development model. However, it was modified to include longer design and shorter testing phases, which is generally expected with Ada developments. One result was that the percentage of time the Ada project spent in the various development activities was very similar to the percentage of time spent in these activities when doing a FORTRAN project. Another finding was the difficulty the Ada team had with unit testing as well as with integration. It was realized that adding additional steps to the design phase, such as an abstract data type analysis, and certain guidelines to the implementation phase, such as to use primarily library units and nest sparingly, would have made development easier. These are among the recommendations made to be incorporated in a new Ada development model next time.

Brophy, Carolyn Elizabeth↗

NASA NTRS1993

Basic studies of baroclinic flows

A fully nonlinear 3-dimensional numerical model (GEOSIM), previously developed and validated for several cases of geophysical fluid flow, has been used to investigate the dynamical behavior of laboratory experiments of fluid flows similar to those of the Earth's atmosphere. The phenomena investigated are amplitude vacillation, and the response of the fluid system to uneven heating and cooling. The previous year's work included hysteresis in the transition between axisymmetric and wave flow. Investigation is also continuing of the flows in the Geophysical Fluid Flow Cell (GFFC), a low-gravity Spacelab experiment. Much of the effort in the past year has been spent in validation of the model under a wide range of external parameters including nonlinear flow regimes. With the implementation of a 3-dimensional upwind differencing scheme, higher spectral resolution, and a shorter time step, the model has been found capable of predicting the majority of flow regimes observed in one complete series of baroclinic annulus experiments of Pfeffer and co-workers. Detailed analysis of amplitude vacillation has revealed that the phase splitting described in the laboratory experiments occurs in some but not all cases. Through the use of animation of the models output, a vivid 3-dimensional view of the phase splitting was shown to the audience of the Southeastern Geophysical Fluid Dynamics Conference in March of this year. A study on interannual variability was made using GEOSIM with periodic variations in the thermal forcing. Thus far, the model has not predicted a chaotic behavior as observed in the experiments, although there is a sensitivity in the wavenumber selection to the initial conditions. Work on this subject, and on annulus experiments with non-axisymmetric thermal heating, will continue. The comparison of GEOSIM's predictions will result from the Spacelab 3 GFFC experiments continued over the past year, on a 'back-burner' basis. At this point, the study (in the form of a draft of a journal article) is nearly completed. The results from GEOSIM compared very well with the experiments, and the use of the model allows the demonstration of flow mechanics that were not possible with the experimental data. For example, animation of the model output shows that the forking of the spiral bands is a transient phenomenon, due to the differential east-west propagation of convection bands from different latitudes.

Miller, Timothy L.↗