Engineering topics

Hawley, John F.

Publications and source records attributed to Hawley, John F..

9 records

NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗

NASA NTRS1992

A powerful local shear instability in weakly magnetized disks. III - Long-term evolution in a shearing sheet. IV - Nonaxisymmetric perturbations

The nonlinear evolution of the recently identified accretion disk magnetic shear instability is investigated through a series of numerical simulations. Finite-difference computations of the equations of compressible MHD are carried out on an axisymmetric shearing sheet system with periodic boundary conditions designed to approximate a local region within an accretion disk. Initial field configurations that include some net vertical component evolve into a nonlinear, exponentially growing solution with large poloidal velocities and magnetic fields with energies comparable to the thermal energy density. The stability of a purely azimuthal field configuration is examined, and it is found that nonaxisymmetric instability is present, but with a growth time measured in tens of orbital periods. In general, the most rapid growth occurs for very small radial and azimuthal wavenumbers, leading to coherent magnetic field structure in planes parallel to the disk. It is suggested that this instability is a key ingredient for the generation of magnetic fields in disks.

Hawley, John F.↗

NASA NTRS1992

Is the Oort A-value a universal growth rate limit for accretion disk shear instabilities?

A weak-field local MHD instability that is of importance to accretion disks is examined. The maximum growth rate of the instability is found to be not only independent of the magnetic field strength but independent of field geometry as well. In particular, all Keplerian disks are unstable in the presence of any weak poloidal field, with the ratio of the maximum growth rate to disk angular velocity given by 3/4. The maximum growth rate of any weak field configuration that is not purely toroidal is given by the local Oort A-value of the disk. The behavior is studied by using a form of the dynamical Hill equations. It is conjectured that the Oort A-value is an upper bound to the growth rate of any instability feeding upon the free energy of differential rotation.

Balbus, Steven A.↗

NASA NTRS1992

A test suite for magnetohydrodynamical simulations

A collection is presented of MHD problems that will provide researchers who have interests in modeling MHD phenomena with a battery of tests (a 'test suite') for calibrating their numerical algorithms. Use of these tests will provide a common reference for comparison of different numerical MHD algorithms. The test suite includes both 1D and 2D problems. Taken together, these problems test the abilities of a numerical method to propagate accurately all of the MHD wave families in moving and stationary media, to capture MHD shocks and contact discontinuities, and to model accurately Lorentz force terms in multiple dimensions. An example solution of each test problem is presented. Diagnostic convergence-testing procedures that provide a quantitative evaluation of MHD algorithms are demonstrated.

Stone, James M.↗

NASA NTRS1991

Three-dimensional simulations of black hole tori

The global hydrodynamic nonaxisymmetric instabilities in thick constant angular momentum (l) accretion gas tori in orbit around a Schwarzschild black hole are investigated via 2D and 3D numerical simulations. A radially-wide torus is found to develop and m = 1 nonaxisymmetric density perturbation near the pressure maximum and a trailing spiral wave that extends through the outer part of the torus. This global wave transports angular momentum through the torus, causing the average angular momentum distribution to evolve slowly away from l = constant. The unstable mode drives an accretion flow from the torus into the black hole. The role of accretion is investigated by modeling a wide torus with an inner boundary at the cusp of the Schwarzschild effective potential. In such a torus, accretion is present throughout the evolution; only modest unstable mode growth is observed, and saturation occurs at low amplitude. Comparison simulations performed in the torus equatorial plane show qualitative similarities between the 2D and 3D system, although the 3D modes have a function dependence on height.

Hawley, John F.↗

NASA NTRS1991

A powerful local shear instability in weakly magnetized disks. I - Linear analysis. II - Nonlinear evolution

A broad class of astronomical accretion disks is presently shown to be dynamically unstable to axisymmetric disturbances in the presence of a weak magnetic field, an insight with consequently broad applicability to gaseous, differentially-rotating systems. In the first part of this work, a linear analysis is presented of the instability, which is local and extremely powerful; the maximum growth rate, which is of the order of the angular rotation velocity, is independent of the strength of the magnetic field. Fluid motions associated with the instability directly generate both poloidal and toroidal field components. In the second part of this investigation, the scaling relation between the instability's wavenumber and the Alfven velocity is demonstrated, and the independence of the maximum growth rate from magnetic field strength is confirmed.

Balbus, Steven A.↗

NASA NTRS1990

Nonaxisymmetric instabilities in a slender torus - Two- and three-dimensional simulations

Nonaxisymmetric instabilities in accretion disks are investigated in terms of the slender, non-self-gravitating, ideal fluid torus model. Nonlinear simulations are carried out in two and tree dimensions, using Cartesian-grid finite differencing. The fastest growing instability is the principal mode. This mode saturates with the formation of ellipsoidal density distributions, designated as 'planets'. The simulations provide evidence for nonlinear mode-mode coupling. When several wavenumbers m are unstable, multiple 'planets' form that subsequently merge. Sources of numerical error are examined, and the effects are contrasted with those of the physical instability. The same qualitative evolution is seen in both the two- and the three-dimensional simulations, even though strict vertical hydrostatic equilibrium is no longer rigorously maintained after mode saturation.

Hawley, John F.↗

NASA NTRS1990

Hollow conical jet models for SS 433 - A paradigm lost?

A precessing jet such as that in SS 433 may be approximated as an axisymmetric flow, if the precession time is short by comparison to the propagation time. A series of simulations has been conducted for precessing jets using an axisymmetric finite-difference hydrodynamics code. Examinations are made of hollow cylindrical jets, which lack the complication of a growing interior volume, and conical jets, which model the behavior of a precessing jet propagating on the surface of its precession cone.

Kochanek, Christopher S.↗

NASA NTRS1989

Vortex paradigm for shock-accelerated density-stratified interfaces

Numerical simulations are reported that capture the observations of Haas and Sturtevant (1985) when a planar shock accelerates an inclined planar interface between two media. A vortex paradigm is presented for interpreting the evolution of shock-accelerated density-stratified interfaces beyond early times. The simulations illustrate the paradigm and are in good agreement with the results of Haas and Sturtevant. Color images of space-time diagrams of the one-space integrated vorticity function are given to show the primary and secondary features of the emerging vortex structures.

Hawley, John F.↗