Numerical simulations of stellar convective dynamos. I. The model and method

Journal Article · Sat Sep 01 04:00:00 EDT 1984 · J. Comput. Phys.; (United States)

DOI:https://doi.org/10.1016/0021-9991(84)90033-0· OSTI ID:6302740

Glatzmaier, G A

A numerical model used to simulate global covection and magnetic field generation in stars is described. Nonlinear, three-dimensional, time-dependent solutions of the anelastic magnetohydrodynamic equations are presented for a strafified, rotating, spherical, fluid shell heated from below. The velocity, magnetic field, and thermodynamic perturbations are expanded in spherical harmonics to resolve their horizontal structure and in Chebyshev polynomials to resolve their radial structure. An explicit Adams--Bashforth time integration scheme is used with an implicit Crank--Nicolson treatment of the diffusion terms. Nonlinear terms are computed in physical space; and spatial derivatives are computed in spectral space. The resulting second-order differential equations are solved with a Chebyshev collocation method. Preliminary solutions for a Solar-like model are briefly discussed. Convective motions driven in the upper, superadiabatic part of the zone penetrate into the lower, subadiabatic part of the zone. Differential rotation, induced by the interaction of convection and rotation, is an equatorial acceleration at the surface as observed on the Sun; below the surface angular velocity decreases with depth. The meridional circulation and the equator-pole temperature excess are within Solar observational constraints; however, the giant-cell velocities at the surface are larger than observed. The large-scale magnetic field, induced by the differntial rotation and helical motions, peaks in the subadiabatic region below the convection zone; and, near the top of the convection zone, it is concentrated over the downdrafts of giant cells. A systematic drift in latitude of the magnetic field or a field reversal has not yet been seen.

• Research Organization: Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico
• OSTI ID: 6302740
• Journal Information: J. Comput. Phys.; (United States), Journal Name: J. Comput. Phys.; (United States) Vol. 55:3; ISSN JCTPA
• Country of Publication: United States
• Language: English

Related Subjects

640102* -- Astrophysics & Cosmology-- Stars & Quasi-Stellar
Radio & X-Ray Sources
71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
COMPUTERIZED SIMULATION
CONVECTION
ENERGY TRANSFER
ENTROPY
FLUID MECHANICS
HEAT TRANSFER
HYDRODYNAMICS
MAGNETIC FIELDS
MAGNETOHYDRODYNAMICS
MASS TRANSFER
MECHANICS
NONLINEAR PROBLEMS
PHYSICAL PROPERTIES
SIMULATION
SPHERICAL HARMONICS
STARS
STRATIFICATION
THERMODYNAMIC PROPERTIES
THREE-DIMENSIONAL CALCULATIONS
TIME DEPENDENCE