Search NASA⌕ Search

Engineering topics

Frank, Adam

Publications and source records attributed to Frank, Adam.

The Radiative Magneto-Hydrodynamics of Colliding Flow: Instabilities, Reconnection and Exoplanet Atmosphere Connections

Hypersonic flows of magnetized plasmas naturally produce strong shocks when obstructions or self‐interactions (flow collisions) occur. Such complex flows occur in a wide variety of HEDP settings such as Z‐pinches and laser driven implosion experiments that include magnetic fields via explosive Helmholtz coils (MIFEDS) When thermal energy behind the shocks can be lost to optically thin radiation, the flow dynamics becomes even more complex, admitting a range of new instabilities such as the Nonlinear Thin Shell Instability (NTSI). Because of their inherent non‐ linearity, multi‐dimensionality and time‐dependence (i.e. turbulence), the ability to articulate the underlying physics of colliding hypersonic magnetized flows represents a significant frontier in HEDP studies. The addition of radiative cooling to these flows propels the study even further into the domain of extreme plasma physics. Exploration of these forms of HED plasmas is of interest to basic plasma physics as well as of interest for their application to Inertial Confinement Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Report DE-SC0020432 PI: Frank, Adam

This goal of this proposal was a focused study to understand multi‐dimensional Colliding Magnetized Radiative Plasma flows (or CMRPs) relevant to basic plasma physics and jets in the context of stellar evolution astrophysics.

Frank, Adam↗

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗

NGC 6302: The Tempestuous Life of a Butterfly

Abstract NGC 6302 (The Butterfly Nebula) is an extremely energetic and rapidly expanding bipolar planetary nebula (PN). If the central source is a single star, then its apparent location in an H-R diagram places it among the most massive, hottest, and presumably rapidly evolving of all central stars of PNe. Our proper motion study of NGC 6302, based on Hubble Space Telescope WFC3 images spanning 11 yr, has uncovered at least four different pairs of uniformly expanding internal lobes ejected at various times and orientations over the past two millennia at speeds ranging from 10–600 km s −1 . In addition, we find a pair of collimated off-axis flows in constant motion at ∼770 ± 100 km s −1 within which bright [Feii]feathersare conspicuous. Combining our results with those previously published, we find that the ensemble of flows has an ionized mass >0.1M ⊙ and its kinetic energy, between 10 46 and 10 48 erg, lies at the upper end of gravity-powered PNe ejection processes such as stellar mergers or mass accretion. We assemble our results into a plausible historical timeline of ejections from the nucleus and suggest that the ejections are powered by gravitational infall.

Astronomy & Astrophysics↗

Excitation and Excavation of the Claws of the Southern Crab

We show that the Southern Crab (aka Hen2–104) presents an auspicious opportunity to study the form and speed of the invisible winds that excavate and shock the lobes of various types of bipolar nebulae associated with close and highly evolved binary stars. A deep three-color image overlay of Hen2–104 reveals that the ionization state of its lobe edges, or “claws,” increases steadily from singly to doubly ionized values with increasing wall latitude. This “reverse” ionization pattern is unique among planetary nebulae (and similar objects) and incompatible with UV photoionization from a central source. We show that the most self-consistent explanation for the ionization pattern is shock ionization by a fast (~600 km s -1 ) “tapered” stellar wind in which the speed and momentum flux of the wind increase with equatorial latitude. We present a hydrodynamic simulation that places the latitude-dependent form, the knotty walls, and the reverse ionization of the outer lobes of Hen2–104 into a unified context.

79 ASTRONOMY AND ASTROPHYSICS↗

Jets from main sequence and white dwarf companions during common envelope evolution

ABSTRACT It has long been speculated that jet feedback from accretion on to the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model without a jet. Our binary consists of a 2-M⊙ red giant branch primary and a 1- or 0.5-M⊙ main sequence (MS) or white dwarf (WD) secondary companion modelled as a point particle. We run the simulations for 10 orbits (40 d). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of the order of the Eddington rate, with maximum velocity vj of the order of the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on the jet mass-loss rate, launch speed, companion mass, opening angle, and accretion rate. In line with our theoretical estimates, jets are choked around the time of first periastron passage and remain choked thereafter. Subsequent to choking, but not before, jets efficiently transfer energy to bound envelope material. This leads to increases in unbound mass of up to $\sim 10{{\ \rm per\ cent}}$, as compared to the simulations without jets. We also estimate the cumulative effects of jets over a full CE phase, finding that jets launched by MS and WD companions are unlikely to dominate envelope unbinding.

79 ASTRONOMY AND ASTROPHYSICS↗

New Clues to the Evolution of Dwarf Carbon Stars From Their Variability and X-Ray Emission

As main-sequence stars with C > O, dwarf carbon (dC) stars are never born alone but inherit carbon-enriched material from a former asymptotic giant branch (AGB) companion. In contrast to M dwarfs in post-mass-transfer binaries, C 2 and/or CN molecular bands allow dCs to be identified with modest-resolution optical spectroscopy, even after the AGB remnant has cooled beyond detectability. Accretion of substantial material from the AGB stars should spin up the dCs, potentially causing a rejuvenation of activity detectable in X-rays. Indeed, a few dozen dCs have recently been found to have photometric variability with periods under a day. However, most of those are likely post-common-envelope binaries, spin–orbit locked by tidal forces, rather than solely spun-up by accretion. Here, we study the X-ray properties of a sample of the five nearest-known dCs with Chandra. Two are detected in X-rays, the only two for which we also detected short-period photometric variability. We suggest that the coronal activity detected so far in dCs is attributable to rapid rotation due to tidal locking in short binary orbits after a common-envelope phase, late in the thermally pulsing (TP) phase of the former C-AGB primary (TP-AGB).

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of charge exchange on the evaporative wind of HD 209458b

ABSTRACT The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet’s host star and neutral hydrogen from the planet’s wind has been proposed as a mechanism to create ‘energetic neutral atoms’ (ENAs), which could explain the high absorption line velocities observed in systems where mass-loss is occurring. In this paper, we present results from three-dimensional hydrodynamic simulations of the mass-loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind–wind interaction analytically, and then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities.

79 ASTRONOMY AND ASTROPHYSICS↗

X-ray Emission from the Pre-planetary Nebula Henize 3-1475

We report the first detection of X-ray emission in a pre-planetary nebula, He 3-1475. Pre-planetary nebulae are rare objects in the short transition stage between the asymptotic giant branch (AGB) and planetary nebula evolutionary phases, and He 3-1475, characterized by a remarkable S-shaped chain of optical knots, is one of the most noteworthy members of this class. Observations with the Advanced CCD Imaging Spectrometer on board the Chandra X-Ray Observatory show the presence of compact emission coincident with the brightest optical knot in this bipolar object, which is displaced from the central star by 2'.7 along the polar axis. Model fits to the X-ray spectrum indicate an X-ray temperature and luminosity, respectively, of (4.3-5.7) x 10(exp 6) K and (4 +/- 1.4) x 10(exp 31) (D/5 kpc)(exp 2) ergs s(exp -1) respectively. Our 3 sigma upper limit on the luminosity of compact X-ray emission from the central star in He 3-1475 is approximately equal to 5 x 10(exp 31) (D/5 kpc)(exp 2) ergs s(exp -1). The detection of X-rays in He 3-1475 is consistent with models in which fast collimated post-AGB outflows are crucial to the shaping of nebulae; we discuss such models in the context of our observations.

circumstellar matter↗

The homunculus of Eta Carinae: An interacting stellar winds paradigm

We simulate the origin and evolution of the bipolar nebula surrounding Eta Car using numerical two-dimensional gasdynamic models. The generalized interacting stellar winds scenario, wherein a stellar wind interacts with an aspherical circumstellar environment, is adopted. The eruption wind of 1840-1860, which is taken to be spherically symmetric, interacts with a preeruption toroidal density environment. Using reasonable assumptions of initial conditions and eruption parameters based on archival data, we have performed over 30 simulations in an effort to bracket the initial parameters which produce models that best match observations. We find that models with high pole-to-equator density contrasts (greater than 100) and toroidal density configurations nicely account for the observed morphology and kinematics of the homunculus.

Frank, Adam↗

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗