Search NASASearch

Engineering topics

Chan, K. L.

Publications and source records attributed to Chan, K. L..

At least 37 records · Page 2

Modeling the QBO and SAO Driven by Gravity Waves

Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW) is applied in a global scale numerical spectral model (NSM) to describe the semi-annual and quasi-biennial oscillations (SAO and QBO) as well as the long term interannual variations that are driven by wave mean flow interactions. This model has been successful in simulating the salient features observed near the equator at altitudes above 20 km, including the QBO extension into the upper mesosphere inferred from UARS measurements. The model has now been extended to describe also the mean zonal and meridional circulations of the upper troposphere and lower stratosphere that affect the equatorial QBO and its global scale extension. This is accomplished in part through tuning of the GW parameterization, and preliminary results lead to the following conclusions: (1) To reproduce the upwelling at equatorial latitudes associated with the Brewer/Dobson circulation that in part is modulated in the model by the vertical component of the Coriolis force, the eddy diffusivity in the lower stratosphere had to be enhanced and the related GW spectrum modified to bring it in closer agreement with the form recommended for the DSP. (2) To compensate for the required increase in the diffusivity, the observed QBO requires a larger GW source that is closer to the middle of the range recommended for the DSP. (3) Through global scale momentum redistribution, the above developments are conducive to extending the QBO and SAO oscillations to higher latitudes. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. (4) In a 3D version of the model, wave momentum is absorbed and dissipated by tides and planetary waves. Thus, a somewhat larger GW source is required to generate realistic amplitudes for the QBO and SAO.

Mayr, H. G.

The Role of Gravity Waves in Modulating Atmospheric Tides

We discuss results for the diurnal and semidiurnal tides obtained from our 3-D, time dependent numerical spectral model (NMS), extending from the ground up into the thermosphere, which incorporates Hines' Doppler spread parameterization of small scale gravity waves (GW). In the DSP, GW momentum (and energy) are conserved as the waves modulate the background flow and are filtered by the flow.As a consequence, the GW interaction tightly couples the dynamic components of the middle atmosphere with strong non-linear interactions between mean zonal circulation, tides and planetary waves to produce complicated patterns of variability much like those observed. The major conclusions are: (1) Since GW momentum is deposited in the altitude regime of increasing winds, the amplitude of the diurnal tide is amplified and its vertical wavelength is reduced at altitudes between 80 and 120 km. Wave filtering by the mean zonal circulation (with peak velocities during solstice) causes the GW flux to peak during equinox, and this produces a large semi-annual variation in the tide that has been observed on UARS. (2) Without the diurnal tide, the semidiurnal tide would also be modulated in this way. But the diurnal tide filters out the GW preferentially during equinox, so that the semidiurnal tide, at higher altitudes, tends to peak during solstice. (3) Under the influence of GW, the tides are modulated also significantly by planetary waves, with periods between 2 and 30 days, which are generated preferentially during solstice in part due to baroclinic instability.

Mayr, H. G.

What maintains the zonal circulation in planetary atmospheres?

The latest Voyager observations have shown large zonal velocities in Neptune's atmosphere, with some indication of alternating jets. Similar wind velocities have also been observed on Venus and are characteristic of planetary atmospheres in general, which is remarkable considering that the available solar or internal heating varies by more than a factor of 1000. A simplified model of the planetary circulation which provides some qualitative understanding is discussed. The basic assumption is that the source driving the circulation is also generating the dissipating eddies which are simulated by diffusion. Thus, the magnitude and structure of the zonal circulation are independent both of the source and the dissipation rate. The zonal velocities are related to the speed of sound and are of comparable magnitude in different atmospheres; although the available heating varies by a large factor, and the planetary parameters vary over a wide range. The alternating jets are described by a convective eigenmode which develops when energy transport out of the planetary interior is important, as is the case on Jupiter, Neptune, and Saturn.

Mayr, H. G.

Venus' superrotation, mixing length theory and eddy diffusion - A parametric study

The concept of the Hadley mechanism is adopted to describe the axisymmetric circulation of the Venus atmosphere. It is shown that, for the atmosphere of a slowly rotating planet such as Venus, a form of the nonliner 'closure' (self-consistent solution) of the fluid dynamics system which constrains the magnitude of the eddy diffusion coefficients can be postulated. A nonlinear one-layer spectral model of the zonally symmetric circulation was then used to establish the relationship between the heat source, the meridional circulation, and the eddy diffusion coefficients, yielding large zonal velocities. Computer experiments indicated that proportional changes in the heat source and eddy diffusion coefficients do not significantly change the zonal velocities. It was also found that, for large eddy diffusion coefficients, the meridional velocity is virtually constant; below a threshold in the diffusion rate, the meridional velocity decreases; and, for large eddy diffusion and small heating rates, the zonal velocities decrease with decreasing planetary rotation rates.

Mayr, H. G.

Mixing-length, shears, and differential rotation

Three aspects of convection in the solar envelope, the mixing length, the shears, and the differential rotation, are considered. Numerical experiments show that the mixing length is scaled by the pressure scale height rather than by the density scale height. Results indicate that the effective viscosity of convective turbulence on large scale shears may be approximated as 1/3 the product of the rms vertical velocity and the pressure scale height. It is suggested that the differential rotation of the solar convection zone may be viewed as an unstable axisymmetric mode which produces slower rotation at the surface.

Chan, K. L.

Turbulent compressible convection in a deep atmosphere. III - Tests on the validity and limitation of the numerical approach

The effects of various approximations on the results of three-dimensional computation of turbulent compressible convection are examined and related to the phenomenon of bubble breaking in stellar convection. The problem of viscosity and the numerical issue of efficiency versus accuracy in such computations are discussed, and a model is formulated and tested with regard to turbulent viscosity and accuracy. Results on the convergence to a statistically stationary state, the effects of time accuracy on the statistics, the effects of the adjustable constant parameter c(mu), and the effects of aspect ratio, horizontal grid size, and vertical grid distribution are discussed. A preliminary look at the three-dimensional flow structure is given.

Chan, K. L.

Conjecture on superrotation in planetary atmospheres - A diffusion model with mixing length theory

The component of rigid shell superrotation on Venus is discussed in the context of comparative planetary atmospheres. A simplified, heuristic analysis, utilizing mixed length theory to describe the small scale nonlinear advections of energy and angular momentum, thereby providing a closure of the dynamic system, is presented, on the basis of which a crude estimate of zonal velocity is made. The rigid shell (global average) component on Venus was calculated to be 105 m/sec.

Mayr, H. G.

Turbulent compressible convection in a deep atmosphere. II - Two-dimensional results for main-sequence A5 and F0 type envelopes

In the present two-dimensional numerical study of turbulent compressible convection in the A5 and F0 main-sequence envelope types, ionization effects are included in the equation of state of the gas, and radiative transfer is modeled in a diffusive process with tabulated gas opacities. It is noted that the thermal effects of ionization significantly affect the dynamics of the flows and that an inversion of the mean density can be created and sustained in a dynamical situation. The substantial differences in the flows of A5 and F0 indicate development trends in the transition from the radiative to the convective mode of energy transport. As convection becomes more effective, the flow becomes more turbulent and the scaling effects of local scale heights become more significant.

Sofia, S.

Differential rotation in a solar-driven quasi-axisymmetric circulation

The concept of a quasi-axisymmetric circulation is used to explore the global scale dynamics of planetary atmospheres. A numerical circulation model applicable to Jupiter is presented, and an analytical study is performed elucidating the conditions leading to differential rotation in an atmosphere which is convectively unstable. A linear system forced by solar differential heating is considered, with nonlinear effects arising from advection being represented in the form of eddy diffusion. An empirical, latitudinal spectrum of the observed zonal wind field on Jupiter is discussed. Numerical solutions are presented which reveal banded wind fields with alternating and equatorial zonal jets and a multicellular Ferrel-Thomson meridional circulation consistent with the observed cloud striations on Jupiter. The vertical derivatives are parameterized to construct a simplified one-layer model.

Mayr, H. G.

What causes the buoyancy reversal in compressible convection?

The problem posed by the existence of a negative buoyancy work region at the top of cellular type convection in a deeply stratified superadiabatic layer (Massaguer and Zahn, 1980) is addressed. It is approached by studying two-dimensional cellular compressible convection with different physical parameters. The results suggest that a large viscosity, together with density stratification, is responsible for the buoyancy reversal. The numerical results obtained are analyzed. It is pointed out, however, that in an astrophysical situation a fluid involved in convection will generally have very small viscosity. It is therefore thought unlikely that buoyancy reversal occurs in this way.

Chan, K. L.

Turbulent compressible convection in a deep atmosphere. I - Preliminary two-dimensional results

The turbulent convection of a compressible fluid in a deep atmosphere is simulated by two-dimensional numerical computations, displaying coexisting 'cells' whose sizes range from the total depth of the convection zone to the smallest scale height at the top. While the largest cells traverse the entire zone, smaller ones are clustered near the top. The vertical correlation length of the vertical velocity is proportional to the local pressure or density scale height, in a way reminiscent of the concept of mixing length.

Chan, K. L.

ADI on staggered mesh - A method for the calculation of compressible convection

An alternating direction implicit (ADI) method has been applied to a staggered grid for the computation of convection in a highly stratified fluid. Since artificial viscosity is not needed, subtle effects like the onset of convection can be studied. These computations compare well with the 2-D results by Graham (1975) and also agree with standard Boussinesq results when taken to that limit. Good efficiency has been achieved with a time step hundreds of times larger than the stability limit imposed by the explicit treatment of diffusion and the Courant number is not restricted to be below 1. The Navier-Stokes equation contains cross spatial derivatives which are treated explicitly in most ADI schemes. The destabilizing effect of such a practice on a 2-D model system with second-order spatial derivative terms only was analyzed and found to be not excessive.

Chan, K. L.

Quick matching technique to study the relationship between solar radius and luminosity variations

A simple matching technique is developed which makes it possible to compute the response of the solar envelope to perturbations which occur within the solar convective region in timescales of importance to climate. The technique is applied to perturbation of the convective efficiency (alpha-mechanism) and of the nongas component of the pressure in different regions of the convection zone (beta-mechanism). The results indicate that whereas either perturbation affects the solar luminosity, the alpha-mechanism has almost no effect on the solar outer radius, regardless of the affected region. The beta-mechanism, however, produces radius changes which may be large if the location of perturbation is deep enough.

Sofia, S.

Estimating short-term solar variations by a simple envelope matching technique

A simple matching technique is explained which allows the computation of the response of the solar surface to perturbations which occur at any depth within the convective envelope of the Sun. This technique was applied to a perturbation of the convective efficiency (alpha-mechanism), and of the non-gas component of the pressure (beta-mechanism) in different regions of the convection zone. The results indicate that either perturbation affects the solar luminosity. However, the alpha-mechanism has little effect in the solar radius, regardless of the location of the perturbed region, whereas the beta-mechanism produces radius changes that become quite large if the location of the perturbed region is deep within the solar convection zone.

Chan, K. L.

Synchrotron brightness distribution of turbulent radio jets

Radio jets are considered as turbulent mixing regions and it is proposed that the essential small scale viscous dissipation in these jets is by emission of MHD waves and by their subsequent strong damping due, at least partly, to gyro-resonant acceleration of supra-thermal particles. A formula relating the synchrotron surface brightness of a radio jet to the turbulent power input is deduced from physical postulates, and is tested against the data for NGC315 and 3C31 (NGC383). The predicted brightness depends essentially on the collimation behavior of the jet, and, to a lesser extent, on the CH picture of a 'high' nozzle with accelerating flow. The conditions for forming a large scale jet at a high nozzle from a much smaller scale jet are discussed. The effect of entrainment on the prediction is discussed with the use of similarity solutions. Although entrainment is inevitably associated with the turbulent jet, it may or may not be a dominant factor depending on the ambient density profile.

Henriksen, R. N.