ALBERT

Description: The 10.7 cm flux patrols in Canada recorded 4 Great Bursts (peaks greater than 500 sfu) during the disk passage of AR 5395 in March 1989. The Great Bursts of 16 and 17 March were simple events of great amplitude and with half-life durations of only several minutes. Earlier Great Bursts, originating on 6 March towards the NE limb and on 10 March closer to the central meridian, belong to an entirely different category of event. Each started with a very strong impulsive event lasting just minutes. After an initial recovery, however, the emission climbed back to level as greater or greater than the initial impulsive burst. The events of 6 and 10 March stayed above the Great Burst threshold for at least 100 minutes. The second component of long duration in these cases is associated with Type 4 continuum emission and thus very likely with CMEs. Major geomagnetic disturbances did not occur as a result of the massive complex event of 6 March or the two simple but strong events of 16 and 17 March. But some 55 hours after the peak in the long-enduring burst of 10 March, a storm began which qualifies as the fourth strongest geomagnetic storm in Canada since 1932. The vertical component of the earth's field measured during the storm by a fluxgate magnetometer at a station in Manitoba is presented. Within a minute of the sudden commencement of this storm, a series of breakdowns began in the transmission system of Hydro-Quebec which resulted in a total loss of power, on a bitterly cold winter's day, for at least 10 hours. The loss of power provoked an enormous outcry from the public resulting in the power utilities being more receptive to the need to monitor solar as well as geomagnetic activity.

Description: A detailed study was made of face and annular seals under conditions where boiling, i.e., phase change of the leaking fluid, occurs within the seal. Many seals operate in this mode because of flashing due to pressure drop and/or heat input from frictional heating. High pressure, water pumps, industrial chemical pumps, and cryogenic pumps are mentioned as a few of many applications. The initial motivation was the LOX-GOX seals for the space shuttle main engine, but the study was expanded to include any face or annular seal where boiling occurs. Some of the distinctive behavior characteristics of two-phase seals were discussed, particularly their axial stability. While two-phase seals probably exhibit instability to disturbances of other degrees of freedom such as wobble, etc., under certain conditions, such analyses are too complex to be treated at present. Since an all liquid seal (with parallel faces) has a neutral axial stiffness curve, and is stabilized axially by convergent coning, other degrees of freedom stability analyses are necessary. However, the axial stability behavior of the two-phase seal is always a consideration no matter how well the seal is aligned and regardless of the speed. Hence, axial stability is thought of as the primary design consideration for two-phase seals and indeed the stability behavior under sub-cooling variations probably overshadows other concerns. The main thrust was the dynamic analysis of axial motion of two-phase face seals, principally the determination of axial stiffness, and the steady behavior of two-phase annular seals. The main conclusions are that seals with two-phase flow may be unstable if improperly balanced. Detailed theoretical analyses of low (laminar) and high (turbulent) leakage seals are presented along with computer codes, parametric studies, and in particular a simplified PC based code that allows for rapid performance prediction. A simplified combined computer code for the performance prediction over the laminar and turbulent ranges of a two-phase seal is described and documented. The analyses, results, and computer codes are summarized.

Description: The present investigation is concerned with a fully discrete accuracy and stability analysis of the one-dimensional heat equation, taking into account the evaluation of two-pass explicit schemes which simultaneously employ lumped and coupled capacity matrices. Schemes of the considered characteristics are not amenable to uncoupled semidiscrete and ordinary differential equation analyses. The obtained results illustrate that superior behavior may be achieved by schemes of the employed type when compared with the performance of the standard one-pass explicit schemes. The key idea in the considered approach is related to the utilization of a reduced-quadrature capacity matrix in the evaluation of the right-hand-side residual.

Description: In order that magnetic flux be confined within the solar interior for times comparable to the solar cycle period it has been suggested that the bulk of the solar toroidal field is stored in the convectively stable overshoot region situated beneath the convection zone proper. Such a magnetic field, though, is still buoyant and is therefore subject to Rayleigh-Taylor type instabilities. The model problem of an isolated region of magnetic field embedded in a convectively stable atmosphere is considered. The fully nonlinear evolution of the two dimensional interchange of modes is studied, thereby shedding some light on one of the processes responsible for the escape of flux from the solar interior.

Description: Motivated by considerations of the solar toroidal magnetic field, the behavior of a layer of uniform magnetic field embedded in a convectively stable atmosphere is studied. Since the field can support extra mass, such a configuration is top-heavy and thus instabilities of the Rayleigh-Taylor type can occur. For both static and rotating basic states, the evolution of the interchange modes (no bending of the field lines) is followed by integrating numerically the nonlinear compressible MHD equations. The initial Rayleigh-Taylor instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The subsequent motions are determined primarily by the strong interactions between vortices which are responsible for the rapid disruption of the magnetic layer.

Description: The current status of streamline-upwind/Petrov-Galerkin (SUPG) methods for the analysis of flow problems is surveyed in an analytical review. Problem areas addressed include classical Galerkin, upwind, artificial-diffusion, SUPG, discontinuous Galerkin, space-time FEM, and discontinuity-capturing approaches to the scalar advection-diffusion equation; incompressible flows; advective-diffusive systems; and the compressible Euler and Navier-Stokes equations. Graphs and diagrams are provided, and the good stability properties of state-of-the-art SUPG methods are pointed out.

Description: Galerkin/least-squares finite-element methods are presented for advective-diffusive equations. Galerkin/least-squares represents a conceptual simplification of streamline-upwind Petrov-Galerkin methods, and is in fact applicable to a wide variety of other problem types. A convergence analysis and error estimates are presented. Some numerical results for compressible Navier-Stokes flows are presented.

Description: There are good reasons for believing that the sun has a strong toroidal magnetic field in the stably stratified region of convective overshoot sandwiched between the radiative zone and convective zone proper. The magnetic field in this region is modeled by studying the behavior of a layer of uniform field embedded in a subadiabatic atmosphere. Since the field can support extra mass, such a configuration is top-heavy, and instabilities of the Rayleigh-Taylor type can occur. Numerical integration of the two-dimensional compressible MHD equations makes it possible to follow the evolution of this instability into the nonlinear regime. The initial buoyancy-driven instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The somewhat surprising subsequent motions are determined primarily by the strong interactions between vortices.

Description: Adding to the classical Hellinger-Reissner formulation, a residual form of the equilibrium equation, a new Galerkin/least-squares finite element method is derived. It fits within the framework of a mixed finite element method and is stable for rather general combinations of stress and velocity interpolations, including equal-order discontinuous stress and continuous velocity interpolations which are unstable within the Galerkin approach. Error estimates are presented based on a generalization of the Babuska-Brezzi theory. Numerical results (not presented herein) have confirmed these estimates as well as the good accuracy and stability of the method.

Description: Symmetric finite-element formulations are presented for the primitive-variables form of the Stokes equations. These formulations are convergent for any combination of pressure and velocity interpolations. Various boundary conditions, such as pressure, are accommodated.