Engineering topics

Cargill, P. J.

Publications and source records attributed to Cargill, P. J..

18 records

NASA NTRS2012

Enthalpy-Based Thermal Evolution of Loops: III. Comparison of Zero-Dimensional Models

Zero dimensional (0D) hydrodynamic models, provide a simple and quick way to study the thermal evolution of coronal loops subjected to time-dependent heating. This paper presents a comparison of a number of 0D models that have been published in the past and is intended to provide a guide for those interested in either using the old models or developing new ones. The principal difference between the models is the way the exchange of mass and energy between corona, transition region and chromosphere is treated, as plasma cycles into and out of a loop during a heating-cooling cycle. It is shown that models based on the principles of mass and energy conservation can give satisfactory results at some, or, in the case of the Enthalpy Based Thermal Evolution of Loops (EBTEL) model, all stages of the loop evolution. Empirical models can lead to low coronal densities, spurious delays between the peak density and temperature, and, for short heating pulses, overly short loop lifetimes.

Cargill, P. J.↗

NASA NTRS2011

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗

NASA NTRS1995

Geometry of interplanetary magnetic clouds

Two dimensional magnetohydrodynamic simulations are presented of the distortion of a magnetic flux rope that is being accelerated through ambient solar wind plasma. The flux rope magnetic field has an axial component parallel to the solar wind field and an azimuthal component, which lies in the simulation plane. As the flux rope moves through the solar wind plasma, vortices form on its trailing edge and couple strongly to its interior. If the flux rope azimuthal field is weak, it deforms into an elongated banana-like shape a few Alfven transit times. A strong azimuthal field component tends to inhibit this distortion. If the flux rope is taken to model a magnetic cloud, it is suggested that the shape of the cloud at 1 AU is determined by its distortion in the inner solar wind. Distortion timescales beyond 1 AU are estimated as many days. It is estimated that effective drag coefficients somewhat greater than unity are appropriate for modelling flux rope propagation.

Cargill, P. J.↗

NASA NTRS1994

Modulational instability of lower hybrid waves at the magnetopause

The role of lower bybrid waesat the magnetopause is reexamined. It is found that for the maximum observed wave power, the lower hybrid waves are unstable to a modulational instability on the magnetosheath side of the magnetopause. The modulational instabitlity leads to localized field structures oriented predominantly along the magnetic field. Such patchy lower hybrid turbulence has been observed by some spacecraft. As a result of the large T(sub i)/T(sub e) ratio, the waves saturate by ion heating; as a result, unlike other settings (e.g. comets, critical ionization phenomena) energetic electrons are not expected. The stochasitc electron transport in the presence of such turbulence is analyzed and results in strongly anistropic electron diffusion, with the dominant direction across the magnetic field. The diffusion rate exceeds significantly that expected from quasi-linear considerations and, for magnetospause parameters, also exceeds the rate discussed by Sonnerup (1980).

Shapiro, V. D.↗

NASA NTRS1991

The interaction of collisionless shocks in astrophysical plasmas

The interaction between collisionless shocks, such as may arise in a number of situations in space plasmas, is examined by means of hybrid numerical simulations. The production of energetic particles by shock collisions is investigated. When the shocks are quasi-perpendicular, ions with energies greater than 10E(0) are produced, where E(0) is the kinetic energy of the unshocked ions. The effect is optimized when the two colliding shocks have roughly equal strengths. As the shocks become closer to quasi-parallel, the number and energy of the accelerated ions increases. Energies in excess of 30E(0) are obtained in the quasi-parallel regime, with a small fall-off when the shocks are exactly parallel. For quasi-parallel collisions, the ion energization can be attributed to the interaction of hot ions with both strong electromagnetic waves and the motional electric field located at the shocks. These results are applied to shock collisions in the distant solar wind, at planetary bow shocks, as well as to the impulsive phase of solar flares.

Cargill, P. J.↗

NASA NTRS1991

Hybrid simulations of rotational discontinuities

1D hybrid simulations of rotational discontinuities (RDs) are presented. When the angle between the discontinuity normal and the magnetic field (theta-BN) is 30 deg, the RD broadens into a quasi-steady state of width 60-80 c/omega-i. The hodogram has a characteristic S-shape. When theta-BN = 60 deg, the RD is much narrower (10 c/omega-i). For right handed rotations, the results are similar to theta-BN = 30 deg. For left handed rotations, the RD does not evolve much from its initial conditions and the S-shape in the hodogram is much less visible. The results can be understood in terms of matching a fast mode wavelike structure upstream of the RD with an intermediate mode one downstream.

Goodrich, C. C.↗

NASA NTRS1990

The merging of quasiperpendicular collisionless shocks

The overtaking of one collisionless shock by another is studied by means of hybrid numerical simulations. The two shocks merge into a stronger shock and trailing nonshock discontinuities. The strong shock continues to propagate in the same direction as the two weaker shocks. The merging is shown to occur by a self-consistent process involving the interaction of ions reflected at the overtaking shock with the plasma upstream of the leading shock. The characteristic time scale for the merging is typically 1/Omega(i), where Omega(i) is the ion gyrofrequency. For exactly perpendicular shocks, the trailing discontinuity is a tangential discontinuity. It has a width of 2-3 ion Larmor radii. For oblique shocks, a contact discontinuity is present in the downstream plasma state. These results are of relevance to shock interactions in the very distant solar wind as well as in other energetic astrophysical situations such as solar flares.

Cargill, P. J.↗

NASA NTRS1990

Hybrid simulations of tangential discontinuities

One-dimensional hybrid simulations of tangential discontinuities (TDs) are presented. Ion kinetic effects are shown to play a crucial role in determining the TD structure. When the density ratio across a TD is large, magnetic profiles with either overshoots or monotonic gradients are obtained. When there is a bulk plasma flow perpendicular to the magnetic field, so that the motional electric field is directed into the moving plasma, the field profiles exhibit similar behavior when the density ratio differs from unity. In each of these cases, the final width is between 2 and 5 ion Larmor radii, roughly what is seen at the earth's magnetopause.

Cargill, P. J.↗

NASA NTRS1988

Lower hybrid waves upstream of comets and their implications for the Comet Halley 'bow wave'

Observed and theoretical features concerning the nature of so-called cometary 'bow shocks' or 'bow waves' are discussed. Collective plasma effects associated with the presence of pickup ring ions (protons and water ions) in the vicinity of the supermagnetosonic to submagnetosonic transition region in the quasi-perpendicular limit are considered; the linear and nonlinear evolution of instabilities around the lower hybrid frequency is emphasized. It is shown that lower hybrid waves can lead to heating and produce distributions with magnitudes in reasonable agreement with Giotto data. The implications to the existence and structure of cometary bow shocks are discussed.

Hizanidis, K.↗

NASA NTRS1988

Resonance absorption of Alfven waves at comet-solar wind interaction regions

The interaction of the low-frequency (0.01 Hz) MHD waves, observed upstream of comets, with the structured plasma near the cometary bow wave is examined. It is suggested that the waves undergo resonant absorption due to either ambient density gradients or localized shear in the background magnetic field. The absorption process can give rise to rapid heating of the solar wind protons, in agreement with observations from Comet Halley. Since the free energy for the generation of MHD waves came from deceleration (without accompanying heating) of the solar wind protons during the pick-up of cometary ions, the subsequent reabsorption of the energy is equivalent to a nonlocal transformation of ordered to random energy and can be described as nonlocal viscosity.

Sharma, A. S.↗

NASA NTRS1988

A mechanism for strong shock electron heating in supernova remnants

It is shown that collisionless shock waves propagating away from a supernova may be directly responsible for the 10 keV X-ray emission seen in supernova remnants. A sequence of plasma instabilities (Buneman and ion acoustic) between the reflected and/or transmitted ions and the background electrons at the foot of the shock front can give rise to rapid anomalous heating of electrons. Hybrid simulations of a perpendicular collisionless shock are presented to demonstrate that this heating can arise within a self-consistently computed shock structure.

Cargill, P. J.↗

NASA NTRS1988

Collisionless shock formation and the prompt acceleration of solar flare ions

The formation mechanisms of collisionless shocks in solar flare plasmas are investigated. The priamry flare energy release is assumed to arise in the coronal portion of a flare loop as many small regions or 'hot spots' where the plasma beta locally exceeds unity. One dimensional hybrid numerical simulations show that the expansion of these 'hot spots' in a direction either perpendicular or oblique to the ambient magnetic field gives rise to collisionless shocks in a few Omega(i), where Omega(i) is the local ion cyclotron frequency. For solar parameters, this is less than 1 second. The local shocks are then subsequently able to accelerate particles to 10 MeV in less than 1 second by a combined drift-diffusive process. The formation mechanism may also give rise to energetic ions of 100 keV in the shock vicinity. The presence of these energetic ions is due either to ion heating or ion beam instabilities and they may act as a seed population for further acceleration. The prompt acceleration of ions inferred from the Gamma Ray Spectrometer on the Solar Maximum Mission can thus be explained by this mechanism.

Cargill, P. J.↗

NASA NTRS1987

The interaction of quasiperpendicular shock waves in a collisionless plasma

The interaction between collisionless shocks is examined from a kinetic viewpoint using computer simulations of shock collisions spanning the range of ambient plasma and shock parameters typical of the interplanetary medium. It is found that the interaction between two subcritical, perpendicular shocks is largely fluidlike with no evidence for accelerated particles. The electromagnetic fields associated with each shock, but not the plasmas, pass through each other during the collision. This interaction is due to self-consistently generated electric fields which adjust to the value needed to bring to rest the two incoming streams. The interaction between two supercritical perpendicular collisionless shocks is similar to the subcritical case, in that the electromagnetic fields associated with each shock are transmitted whereas particles are not. However, the interaction is also highly dynamic with substantial acceleration of ions. Increasing the Mach number gives rise to an increasing number of energetic particles.

Cargill, P. J.↗

NASA NTRS1986

Flare energetics

In this investigation of flare energetics, researchers sought to establish a comprehensive and self-consistent picture of the sources and transport of energy within a flare. To achieve this goal, they chose five flares in 1980 that were well observed with instruments on the Solar Maximum Mission, and with other space-borne and ground-based instruments. The events were chosen to represent various types of flares. Details of the observations available for them and the corresponding physical parameters derived from these data are presented. The flares were studied from two perspectives, the impulsive and gradual phases, and then the results were compared to obtain the overall picture of the energics of these flares. The role that modeling can play in estimating the total energy of a flare when the observationally determined parameters are used as the input to a numerical model is discussed. Finally, a critique of the current understanding of flare energetics and the methods used to determine various energetics terms is outlined, and possible future directions of research in this area are suggested.

Wu, S. T.↗

NASA NTRS1986

The energetics of the gradual phase

Reseachers compare results with those in the chapter by Moore et al. (1980), who reached five main conclusions about the gradual phase: (1) the typical density of the soft X-ray emitting plasma is between 10 to the 11th power and 10 to the 12th power cm-3 for compact flares and between 10 to the 10th power and 10 to the 11th power cm-3 for a large-area flare; (2) cooling is by conduction and radiation in roughly equal proportions; (3) continual heating is needed in the decay phase of two-ribbon flares; (4) continual heating is probably not needed in compact events; (5) most of the soft-X-ray-emitting plasma results from chromospheric evaporation. The goal was to reexamine these problems with the data from the Solar Maximum Mission (SMM) and other supporting instruments as well as to take advantage of recent theoretical advances. SMM is capable of measuring coronal temperatures more accurately and with a better cadence than has been possible before. The SMM data set is also unique in that the complete transit of an active region was observed, with soft X-ray and UV images being taken every few minutes. Researcher's were therefore able to establish the pre-flare conditions of the region and see whether anything has changed as a result of the flare. The assumptions made in attempting to determine the required plasma parameters are described. The derived parameters for the five prime flares are presented, and the role of numerical simulations is discussed.

Strong, K. T.↗

NASA NTRS1986

The energy balance of plasmoids in the solar atmosphere

The properties of an isolated magnetized plasmoid in a nonuniform magnetic field such as arises in stellar atmospheres are studied. The work of Pneuman and Cargill (1985) on the so-called melon-seed effect is extended to include an equation describing the energy balance, so giving a unified picture of the shape, motion, and energetics of the plasmoid. Three treatments of plasmoid energy balance are considered: (1) a polytropic law, (P = about N to the gamma); (2) one in which the plasmoid cools radiatively; and (3) one in which a heating function proportional to the local density balances the radiation. For a gamma = 4/3 polytrope the evolution is self-similar, so that the plasmoid maintains its shape as it moves out from the stellar surface. If gamma is less than 4/3, the final shape is a long thin cigar-shaped body, whereas if gamma is greater than or equal to 4/3, it ultimately becomes self-similar. In cases with radiation and also with heating, the ultimate shape of the plasmoid is determined by whether its gas or magnetic pressure dominate. The former is equivalent to the gamma-less-than-4/3 polytrope, and the latter to the gamma-greater-than-4/3 one. If radiation alone is present, the plasmoid cools rapidly and subsequently evolves self-similarly. If heating balances radiation initially, then the plasmoid heats up as it moves out, but, if the ratio of the transit of time of Alfven waves across it is much less than the radiative cooling time, it ultimately evolves as a gamma = 5/3 polytrope. In each case the plasmoid can be ejected to large distances (several radii) in a stellar atmosphere, for a reasonable choice of surface parameters.

Cargill, P. J.↗

NASA NTRS1986

Interaction of two collisionless shocks

Kinetic simulations of the interaction between two collisionless shocks are presented. During the collision of two perpendicular shocks, the shock electromagnetic field structures pass through each other, while the previously shocked ions are kept separate by the electric field arising in the collision. When two supercritical shocks collide, a fraction of ions are accelerated up to an order of magnitude in energy by first being reflected at one shock, then interacting with the electric fields of the other shock.

Cargill, P. J.↗

NASA NTRS1984

The development and cooling of a solar limb-flare

Observations of a flare that began in soft X-rays at 20:37 UT on April 12, 1980 at the west limb of the sun are discussed. The observations of the flare and postflare loops are first described, and the Solar Maximum Mission data are interpreted in terms of the temperature, density, and geometry of the emitting regions. The observed postflare cooling time is compared with that expected from radiation, conduction, and enthalpy flux. The loop model is discussed, and the observed events are compared with other proposed models.

Veck, N. J.↗