Engineering topics

Toomre, Juri

Publications and source records attributed to Toomre, Juri.

At least 19 records

NASA NTRS2003

Global and Local Helioseismic Studies of Solar Convection Zone Dynamics Using SOI-MDI on SOHO

Our joint collaborative analyses of global mode data to characterize the solar differential rotation (e.g. Thompson et al. 1996, Schou et al. 1998), and most recently to detect and analyze temporal variations in angular velocity Omega profiles both within the convection zone and in the deeper radiative interior (e.g. Howe et al 2000a,b; Toomre et al. 2000), have led to a series of fascinating discoveries. These should be pursued further as the solar cycle continues. The physical deductions being made from these studies have been greatly strengthened by utilizing both SOI-MDI and GONG data in order to have two independent observational realizations of Doppler images spanning a five-year interval, using two separate procedures to determine global mode splittings, and then analyzing those splitting data sets using both RLS and SOLA inversion procedures. There are considerable subtleties in the effects of instrumental response functions and calibrations, sensitivity of peak finding algorithms and their mode leakage estimates, and stochastic variations in mode amplitudes that can all contribute to apparent changes in the Omega profiles being inferred from sequences of helioseismic data. We have come to understand the implications of many of these calibration and analysis steps, greatly aided by frequent multi-week collaborative working sessions in our Helioseismic Analysis Facility (HAF) at JILA involving many members of the SO1 dynamics and inversion team, including most of our Co-Is during the summer months when we hold intensive working sessions. Considerable further focused attention is required in a collaborative setting on such global mode issues as we continue studying the changing sun.

Toomre, Juri↗

NASA NTRS1999

Near-Surface Flow Fields Deduced Using Correlation Tracking and Time-Distance Analysis

Near-photospheric flow fields on the Sun are deduced using two independent methods applied to the same time series of velocity images observed by SOI-MDI on SOHO. Differences in travel times between f modes entering and leaving each pixel measured using time-distance helioseismology are used to determine sites of supergranular outflows. Alternatively, correlation tracking analysis of mesogranular scales of motion applied to the same time series is used to deduce the near-surface flow field. These two approaches provide the means to assess the patterns and evolution of horizontal flows on supergranular scales even near disk center, which is not feasible with direct line-of-sight Doppler measurements. We find that the locations of the supergranular outflows seen in flow fields generated from correlation tracking coincide well with the locations of the outflows determined from the time-distance analysis, with a mean correlation coefficient after smoothing of bar-r(sub s) = 0.840. Near-surface velocity field measurements can used to study the evolution of the supergranular network, as merging and splitting events are observed to occur in these images. The data consist of one 2048-minute time series of high-resolution (0.6" pixels) line-of-sight velocity images taken by MDI on 1997 January 16-18 at a cadence of one minute.

DeRosa, Marc↗

NASA NTRS1998

Phase Inversion: Inferring Solar Subphotospheric Flow and Other Asphericity from the Distortion of Acoustic Waves

A method is proposed for analyzing an almost monochromatic train of waves propagating in a single direction in an inhomogeneous medium that is not otherwise changing in time. An effective phase is defined in terms of the Hilbert transform of the wave function, which is related, via the JWKB approximation, to the spatial variation of the background state against which the wave is propagating. The contaminating effect of interference between the truly monochromatic components of the train is eliminated using its propagation properties. Measurement errors, provided they are uncorrelated, are manifest as rapidly varying noise; although that noise can dominate the raw phase-processed signal, it can largely be removed by low-pass filtering. The intended purpose of the analysis is to determine the distortion of solar oscillations induced by horizontal structural variation and material flow. It should be possible to apply the method directly to sectoral modes. The horizontal phase distortion provides a measure of longitudinally averaged properties of the Sun in the vicinity of the equator, averaged also in radius down to the depth to which the modes penetrate. By combining such averages from different modes, the two-dimensional variation can be inferred by standard inversion techniques. After taking due account of horizontal refraction, it should be possible to apply the technique also to locally sectoral modes that propagate obliquely to the equator and thereby build a network of lateral averages at each radius, from which the full three-dimensional structure of the Sun can, in principle, be determined as an inverse Radon transform.

Gough, Douglas↗

NASA NTRS1998

Turbulent Compressible Convection with Rotation: Mean Flows and Differential Rotation - 2

The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Brummell, Nicholas H.↗

NASA NTRS1998

Confined States in Large-Aspect-Ratio Thermosolutal Convection

Two-dimensional thermosolutal convection with no-slip boundary conditions is studied using numerical simulations in a periodic domain. The domain is large enough to follow the evolution of phase instabilities of fully nonlinear traveling waves. In the parameter regime studied these instabilities evolve, without loss of phase or hysteresis, into a series of confined states or pulses characterized by locally enhanced heat and solute transport. The wavelength and phase velocity of the traveling rolls within a pulse differ substantially from those in the background. The pulses drift in the same direction as the convection rolls on which they ride but more slowly, and are characterized by an exponential leading front and an oscillatory trailing end. Multiple, apparently stable, states are found for identical parameter values. The qualitative properties of the pulses are in good agreement with the predictions of a third-order phase equation which accounts for the relation between wave number and phase velocity, the oscillatory tails and the multiplicity of states. These properties of the pulses are shown to be a consequence of Shil'nikov dynamics in the spatial domain.

Spina, Alejandro↗

NASA NTRS1996

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗

NASA NTRS1995

The sensitivity of various mode sets for probing differential rotation shear zones

The potential of future data sets from global oscillations network group (GONG) and solar oscillations investigation (SOI) for resolving long-lived azimuthal jets and shearing flows, is investigated. Various artificial data sets are constructed, containing noise resembling that of a one-year observation run. These are inverted using a two dimensional regularized least squares inversion. The ability of this method to form well localized averages of the rotation rate, as measured by the averaging kernels, is investigated using an extensive mode set and subsets. It is shown that it is possible to keep the noise in the solution down to a few nHz in much of the solar interior, while obtaining a reasonable resolution for a GONG-like data set. At low latitudes in the middle of the convection zone, an angular resolution of less than 10 deg and a radial resolution of about 0.04 solar radii, are obtained. The averaging kernels depend on the mode set, and a reduction in the number of modes tends to introduce small-scale near surface structures into the averaging kernels which would adversely affect the inferred rotation rate.

Christensen-Dalsgaard, Joergen↗

NASA NTRS1995

Local area analysis of high-degree solar oscillations: New ring fitting procedures

The local area analysis of five-minute solar oscillations using ring diagrams to determine subphotospheric velocity flows is a tool for convection zone dynamics. In relation to the problem of the large computational task of fitting the rings, a faster method is presented that carries out the ring fitting using data obtained with a high l helioseismometer. Noise sources are eliminated, and a perturbation approach is used to fit the azimuthally averaged spectrum. The parameters determined in this way are held constant while the ring diagram is fitted. The results obtained are presented and discussed.

Haber, Deborah A.↗

NASA NTRS1995

Plane-wave analysis of solar acoustic-gravity waves: A (slightly) new approach

The plane-wave decomposition of the acoustic-gravity wave effects observed in the photosphere provides a computationally efficient technique that probes the structure of the upper convective zone and boundary. In this region, the flat sun approximation is considered as being reasonably accurate. A technique to be used for the systematic plane-wave analysis of Michelson Doppler imager data, as part of the solar oscillations investigation, is described. Estimates of sensitivity are presented, and the effects of using different planar mappings are discussed. The technique is compared with previous approaches to the three dimensional plane-wave problem.

Bogart, Richard S.↗

NASA NTRS1995

Analysis of helioseismic wave fields to examine horizontal structures

A preliminary inversion procedure that carries out a local area analysis on simulated oscillation data is presented and discussed. The procedure is carried out in order to deduce two dimensional subsurface structures in the horizontal plane, representative of thermal variations, as a function of depth. The aim is the evaluation of an inversion procedure that utilizes information gained from the phase distortion occurring in artificially generated acoustic waves in order to determine the subsurface thermal structure. These distortions would naturally occur as a direct consequence of convective motions in the solar interior.

Julien, K. A.↗

NASA NTRS1995

Local-area simulations of rotating compressible convection and associated mean flows

The dynamics of compressible convection within a curved local segment of a rotating spherical shell are considered in relation to the turbulent redistribution of angular momentum within the solar convection zone. Current supercomputers permit fully turbulent flows to be considered within the restricted geometry of local area models. By considering motions in a curvilinear geometry in which the Coriolos parameters vary with latitude, Rossby waves which couple with the turbulent convection are thought of as being possible. Simulations of rotating convection are presented in such a curved local segment of a spherical shell using a newly developed, sixth-order accurate code based on compact finite differences.

Hurlburt, Neal E.↗

NASA NTRS1995

Non-local transport in turbulent MHD convection

The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Miesch, Mark↗

NASA NTRS1995

Solar convection zone dynamics and rotation

Helioseismology reveals that there appears to be a discrepancy between the differential rotation profile with the radius and latitude deduced from the inversion of the observed frequency splitting of p-modes and the rotation profile anticipated from various global simulations of convection in rotating shells. The problem may be caused by deductions in numerical simulations of convection which was considered near laminar flows. A high performance computing offers paths for studying the properties of such astrophysical turbulence. A range of approaches to study the basic dynamics of convection is reviewed, along with its ability to redistribute angular momentum in rotating systems. The results of recent three-dimensional simulations of turbulent compressible convection constrained by the effects of rotation are described for the cases in spherical shells and local area f-planes. The turbulence possesses intense vortex tubes with intricate interactions and instabilities. The theoretical studies reveal that the transition to turbulent states has contributed to significant changes in flow symmetries and transports, and that such a turbulent convection can drive substantial mean flows, which are distinctly different from those in which the convection is dominantly laminar.

Toomre, Juri↗

NASA NTRS1994

Penetration below a convective zone

Two-dimensional numerical simulations are used to investigate how fully compressible nonlinear convection penetrates into a stably stratified zone beneath a stellar convection zone. Estimates are obtained of the extent of penetration as the relative stability S of the stable to the unstable zone is varied over a broad range. The model deals with a perfect gas possessing a constant dynamic viscosity. The dynamics is dominated by downward-directed plumes which can extend far into the stable material and which can lead to the excitation of a broad spectrum of internal gravity waves in the lower stable zone. The convection is highly time dependent, with the close coupling between the lateral swaying of the plumes and the internal gravity waves they generate serving to modulate the strength of the convection. The depth of penetration delta, determined by the position where the time-averaged kinetic flux has its first zero in the stable layer, is controlled by a balance between the kinetic energy carried into the stable layer by the plumes and the buoyancy braking they experience there. A passive scalar is introduced into the unstable layer to evaluate the transport of chemical species downward. Such a tracer is effectively mixed within a few convective overturning times down to a depth of delta within the stable layer. Analytical estimates based on simple scaling laws are used to interpret the variation of delta with S, showing that it first involves an interval of adiabatic penetration if the local Peclet number of the convection exceeds unity, followed by a further thermal adjustment layer, the depths of each interval scaling in turn as S(exp -1) and S(exp -1/4). These estimates are in accord with the penetration results from the simulations.

Hurlburt, Neal E.↗

NASA NTRS1993

Ionization effects in three-dimensional solar granulation simulations

These numerical studies show that ionization influences both the transport and dynamical properties of compressible convection near the surface of the Sun. About two-thirds of the enthalpy transported by convective motions in the region of partial hydrogen ionization is carried as latent heat. The role of fast downflow plumes in total convective transport is substantially elevated by this contribution. Instability of the thermal boundary layer is strongly enhanced by temperature sensitive variations in the radiative properties of the fluid, and this provides a mechanism for plume initiation and cell fragmentation in the surface layers. As the plumes descend, temperature fluctuations and associated buoyancy forces are maintained because of the increased specific heat of the partially ionized material. This can result is supersonic vertical flows. At greater depths, ionization effects diminish, and the plumes are decelerated by significant entrainment of surrounding fluid.

Rast, Mark P.↗

NASA NTRS1993

Large-eddy simulations of compressible convection on massively parallel computers

We report preliminary implementation of the large-eddy simulation (LES) technique in 2D simulations of compressible convection carried out on the CM-2 massively parallel computer. The convective flow fields in our simulations possess structures similar to those found in a number of direct simulations, with roll-like flows coherent across the entire depth of the layer that spans several density scale heights. Our detailed assessment of the effects of various subgrid scale (SGS) terms reveals that they may affect the gross character of convection. Yet, somewhat surprisingly, we find that our LES solutions, and another in which the SGS terms are turned off, only show modest differences. The resulting 2D flows realized here are rather laminar in character, and achieving substantial turbulence may require stronger forcing and less dissipation.

Xie, Xin↗

NASA NTRS1991

Laboratory and theoretical models of planetary-scale instabilities and waves

Meteorologists and planetary astronomers interested in large-scale planetary and solar circulations recognize the importance of rotation and stratification in determining the character of these flows. The two outstanding problems of interest are: (1) the origins and nature of chaos in baroclinically unstable flows; and (2) the physical mechanisms responsible for high speed zonal winds and banding on the giant planets. The methods used to study these problems, and the insights gained, are useful in more general atmospheric and climate dynamic settings. Because the planetary curvature or beta-effect is crucial in the large scale nonlinear dynamics, the motions of rotating convecting liquids in spherical shells were studied using electrohydrodynamic polarization forces to generate radial gravity and centrally directed buoyancy forces in the laboratory. The Geophysical Fluid Flow Cell (GFFC) experiments performed on Spacelab 3 in 1985 were analyzed. The interpretation and extension of these results have led to the construction of efficient numerical models of rotating convection with an aim to understand the possible generation of zonal banding on Jupiter and the fate of banana cells in rapidly rotating convection as the heating is made strongly supercritical. Efforts to pose baroclinic wave experiments for future space missions using a modified version of the 1985 instrument have led us to develop theoretical and numerical models of baroclinic instability. Some surprising properties of both these models were discovered.

Hart, John E.↗

NASA NTRS1991

Nonlinear behavior of solar gravity modes driven by He-3 in the core. II - Numerical simulations

The nonlinear behavior of gravity-mode oscillations driven by He-3-destroying reactions in the solar core has been examined by numerically integrating equations describing a very simplified model. The results of a previous bifurcation analysis, which suggest that such oscillations are unlikely to attain amplitudes sufficient to trigger core convection, are verified. These results are extended to models whose nuclear reaction rates and thermal stratification represent the core somewhat more accurately. Nonlinear processes give rise to a preference for the oscillations to develop as standing waves rather than traveling waves, thus breaking the degeneracy between these two types of motion which exists in linearized theory. Study of the large-amplitude behavior of the oscillations is hindered by a tendency for the model to become thermally unstable.

Merryfield, William J.↗