BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Bryan M. Johnson

Bryan M. Johnson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-phphysics.flu-dynphysics.plasm-phastro-ph.GAnlin.PS

Catalog footprint

What is connected

8works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2016arXiv

Three temperature plasma shock solutions with gray radiation diffusion

The effects of radiation on the structure of shocks in a fully-ionized plasma are investigated by solving the steady-state fluid equations for ions, electrons, and radiation. The electrons and ions are assumed to have the same bulk velocity but separate temperatures, and the radiation is modeled with the gray-diffusion approximation. Both electron and ion conduction are included, as well as ion viscosity. When the material is optically thin, three-temperature behavior occurs. When the diffusive flux of radiation is important but radiation pressure is not, two-temperature behavior occurs, with the electrons strongly coupled to the radiation. Since the radiation heats the electrons on length scales that are much longer than the electron-ion Coulomb coupling length scale, these solutions resemble radiative shock solutions rather than plasma shock solutions that neglect radiation. When radiation pressure is important, all three components are strongly coupled. Results with constant values for the transport and coupling coefficients are compared to a full numerical simulation with a good match between the two, demonstrating that steady shock solutions constitute a straightforward and comprehensive verification test methodology for multi-physics numerical algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Closed-form shock solutions

It is shown here that a subset of the implicit analytical shock solutions discovered by Becker and by Johnson can be inverted, yielding several exact closed-form solutions of the one-dimensional compressible Navier-Stokes equations for an ideal gas. For a constant dynamic viscosity and thermal conductivity, and at particular values of the shock Mach number, the velocity can be expressed in terms of a polynomial root. For a constant kinematic viscosity, independent of Mach number, the velocity can be expressed in terms of a hyperbolic tangent function. The remaining fluid variables are related to the velocity through simple algebraic expressions. The solutions derived here make excellent verification tests for numerical algorithms, since no source terms in the evolution equations are approximated, and the closed-form expressions are straightforward to implement. The solutions are also of some academic interest as they may provide insight into the non-linear character of the Navier-Stokes equations and may stimulate further analytical developments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Interaction between Turbulence and a Planar Rarefaction

The modeling of turbulence, whether it be numerical or analytical, is a difficult challenge. Turbulence is amenable to analysis with linear theory if it is subject to rapid distortions, i.e., motions occurring on a time scale that is short compared to the time scale for non-linear interactions. Such an approach could prove useful for understanding aspects of astrophysical turbulence, which is often subject to rapid distortions, such as supernova explosions or the free-fall associated with gravitational instability. As a proof of principle, a particularly simple problem is considered here: the evolution of vorticity due to a planar rarefaction in an ideal gas. Vorticity can either grow or decay in the wake of a rarefaction front, and there are two competing effects that determine which outcome occurs: entropy fluctuations couple to the mean pressure gradient to produce vorticity via baroclinic effects, whereas vorticity is damped due to the conservation of angular momentum as the fluid expands. In the limit of purely entropic fluctuations in the ambient fluid, a strong rarefaction generates vorticity with a turbulent Mach number on the order of the root-mean square of the ambient entropy fluctuations. The analytical results are shown to compare well with results from two- and three-dimensional numerical simulations. Analytical solutions are also derived in the linear regime of Reynolds-averaged turbulence models. This highlights an inconsistency in standard turbulence models that prevents them from accurately capturing the physics of rarefaction-turbulence interaction. Finally, dimensional analysis of the equations indicates that rapid distortion of turbulence can give rise to two distinct regimes in the turbulent spectrum: a distortion range at large scales where linear distortion effects dominate, and an inertial range at small scales where non-linear effects dominate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analytical shock solutions at large and small Prandtl number

Exact one-dimensional solutions to the equations of fluid dynamics are derived in the large-Pr and small-Pr limits (where Pr is the Prandtl number). The solutions are analogous to the Pr = 3/4 solution discovered by Becker and analytically capture the profile of shock fronts in ideal gases. The large-Pr solution is very similar to Becker's solution, differing only by a scale factor. The small-Pr solution is qualitatively different, with an embedded isothermal shock occurring above a critical Mach number. Solutions are derived for constant viscosity and conductivity as well as for the case in which conduction is provided by a radiation field. For a completely general density- and temperature-dependent viscosity and conductivity, the system of equations in all three limits can be reduced to quadrature. The maximum error in the analytical solutions when compared to a numerical integration of the finite-Pr equations is O(1/Pr) for large Pr and O(Pr) for small Pr.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Locality of MHD Turbulence in Isothermal Disks

We numerically evolve turbulence driven by the magnetorotational instability (MRI) in a 3D, unstratified shearing box and study its structure using two-point correlation functions. We confirm Fromang and Papaloizou's result that shearing box models with zero net magnetic flux are not converged; the dimensionless shear stress $α$ is proportional to the grid scale. We find that the two-point correlation of the magnetic field shows that it is composed of narrow filaments that are swept back by differential rotation into a trailing spiral. The correlation lengths along each of the correlation function principal axes decrease monotonically with the grid scale. For mean azimuthal field models, which we argue are more relevant to astrophysical disks than the zero net field models, we find that: $α$ increases weakly with increasing resolution at fixed box size; $α$ increases slightly as the box size is increased; $α$ increases linearly with net field strength, confirming earlier results; the two-point correlation function of the magnetic field is resolved and converged, and is composed of narrow filaments swept back by the shear; the major axis of the two-point increases slightly as the box size is increased; these results are code independent, based on a comparison of ATHENA and ZEUS runs. The velocity, density, and magnetic fields decorrelate over scales larger than $\sim H$, as do the dynamical terms in the magnetic energy evolution equations. We conclude that MHD turbulence in disks is localized, subject to the limitations imposed by the absence of vertical stratification, the use of an isothermal equation of state, finite box size, finite run time, and finite resolution

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Simple Waves in Ideal Radiation Hydrodynamics

In the dynamic diffusion limit of radiation hydrodynamics, advection dominates diffusion; the latter primarily affects small scales and has negligible impact on the large scale flow. The radiation can thus be accurately regarded as an ideal fluid, i.e., radiative diffusion can be neglected along with other forms of dissipation. This viewpoint is applied here to an analysis of simple waves in an ideal radiating fluid. It is shown that much of the hydrodynamic analysis carries over by simply replacing the material sound speed, pressure and index with the values appropriate for a radiating fluid. A complete analysis is performed for a centered rarefaction wave, and expressions are provided for the Riemann invariants and characteristic curves of the one-dimensional system of equations. The analytical solution is checked for consistency against a finite difference numerical integration, and the validity of neglecting the diffusion operator is demonstrated. An interesting physical result is that for a material component with a large number of internal degrees of freedom and an internal energy greater than that of the radiation, the sound speed increases as the fluid is rarefied. These solutions are an excellent test for radiation hydrodynamic codes operating in the dynamic diffusion regime. The general approach may be useful in the development of Godunov numerical schemes for radiation hydrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Magnetohydrodynamic Shearing Waves

I consider the nonaxisymmetric linear theory of a rotating, isothermal magnetohydrodynamic (MHD) shear flow. The analysis is performed in the shearing box, a local model of a thin disk, using a decomposition in terms of shearing waves, i.e., plane waves in a frame comoving with the shear. These waves do not have a definite frequency as in a normal mode decomposition, and numerical integration of a coupled set of amplitude equations is required to characterize their time dependence. Their generic time dependence, however, is oscillatory with slowly-varying frequency and amplitude, and one can construct accurate analytical solutions by applying the Wentzel-Kramers-Brillouin method to the full set of amplitude equations. The solutions have the following properties: 1) Their accuracy increases with wavenumber, so that most perturbations that fit within the disk are well-approximated as modes with time-dependent frequencies and amplitudes. 2) They can be broadly classed as incompressive and compressive perturbations, the former including the nonaxisymmetric extension of magnetorotationally unstable modes, and the latter being the extension of fast and slow modes to a differentially-rotating medium. 3) Wave action is conserved, implying that their energy varies with frequency. 4) Their shear stress is proportional to the slope of their frequency, so that they transport angular momentum outward (inward) when their frequency increases (decreases). The complete set of solutions constitutes a comprehensive linear test suite for numerical MHD algorithms that incorporate a background shear flow. I conclude with a brief discussion of possible astrophysical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Effects of Thermal Conduction on Radiatively-Inefficient Accretion Flows

We quantify the effects of electron thermal conduction on the properties of hot accretion flows, under the assumption of spherical symmetry. Electron heat conduction is important for low accretion rate systems where the electron cooling time is longer than the conduction time of the plasma, such as Sgr A* in the Galactic Center. For accretion flows with density profiles similar to the Bondi solution (n[r] ~ r^[-3/2]), we show that heat conduction leads to super-virial temperatures, implying that conduction significantly modifies the dynamics of the accretion flow. We then self-consistently solve for the dynamics of spherical accretion in the presence of saturated conduction and electron heating. We find that the accretion rate onto the central object can be reduced by ~1-3 orders of magnitude relative to the canonical Bondi rate. Electron conduction may thus be an important ingredient in explaining the low radiative efficiencies and low accretion rates inferred from observations of low-luminosity galactic nuclei. The solutions presented in this paper may also describe the nonlinear saturation of the magnetothermal instability in hot accretion flows.

0 citations0 reviewsNo code linkedOpen access