BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. R. King

J. R. King contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRphysics.plasm-phphysics.comp-phmath-phmath.APmath.MPnlin.SIphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 published item(s)

preprint2021arXiv

Cross-Code verification and sensitivity analysis to effectively model the electrothermal instability

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. A sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Growing Neoclassical Tearing Modes Seeded via Transient-Induced-Multimode Interactions

Nonlinear extended MHD simulations demonstrating seeding of neoclassical tearing modes (NTMs) via MHD-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge [R.J. La Haye, et al., Proceedings IAEA FEC 2020]. The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase like that observed experimentally. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth which is well described by the nonlinear modified Rutherford equation. This work highlights the role of nonlinear mode coupling in seeding NTMs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multidimensional tests of a finite-volume solver for MHD with a real-gas equation of state

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation results where the spurious solution leads to an unphysical simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pattern formation in a pseudo-parabolic equation

We address the propagation into an unstable state of a localised disturbance in a forward-backward diffusion pseudo-parabolic equation. Three asymptotic regimes are distinguished as t tends to infinity, the first being a regime ahead of the propagating disturbance that is dominated by the linearised equation. The analysis of this leads to the determination of the speed of the leading edge of the propagating disturbance and implies that in the second, transition, regime the solution takes the form of a modulated travelling wave. In a third regime the solution approaches a nearly periodic steady state, where the period is obtained on matching with the modulated travelling wave. Detailed analysis of this pattern is also presented. The analysis is completed by contrasting the formal asymptotic description of the solution with numerical computations. It is assumed for the above analysis that the initial disturbance decays faster than an exponential rate; in this case a critical exponential decay rate at the leading edge of the front and propagation speed are found. We investigate the wave speed selection mechanism for exponentially decaying initial conditions. It is found that whenever the initial data behave as a real exponential (no matter how slow the rate of the decay) the speed selected is that selected by fast decaying initial conditions. However, for initial conditions with a complex exponential we find regimes of the decay rate and the wavelength for which the front propagates at a faster wave speed. This is investigated numerically and is worth emphasising since it gives a different scenario for wave speed behaviour than that exhibited by well-studied semilinear reaction-diffusion equations: there are initial conditions with exponential decay faster than the critical one for which the front propagates with a speed faster than the critical one.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lie symmetry properties of nonlinear reaction-diffusion equations with gradient-dependent diffusivity

Complete descriptions of the Lie symmetries of a class of nonlinear reaction-diffusion equations with gradient-dependent diffusivity in one and two space dimensions are obtained. A surprisingly rich set of Lie symmetry algebras depending on the form of diffusivity and source (sink) in the equations is derived. It is established that there exists a subclass in 1-D space admitting an infinite-dimensional Lie algebra of invariance so that it is linearisable. A special power-law diffusivity with a fixed exponent, which leads to wider Lie invariance of the equations in question in 2-D space, is also derived. However, it is shown that the diffusion equation without a source term (which often arises in applications and is sometimes called the Perona-Malik equation) possesses no rich variety of Lie symmetries depending on the form of gradient-dependent diffusivity. The results of the Lie symmetry classification for the reduction to lower dimensionality, and a search for exact solutions of the nonlinear 2-D equation with power-law diffusivity, also are included.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Detailed Abundances of the Solar Twins 16 Cygni A and B: Constraining Planet Formation Models

Results of a detailed abundance analysis of the solar twins 16 Cyg A and 16 Cyg B based on high-resolution, high signal-to-noise ratio echelle spectroscopy are presented. 16 Cyg B is known to host a giant planet while no planets have yet been detected around 16 Cyg A. Stellar parameters are derived directly from our high-quality spectra, and the stars are found to be physically similar, with Delta_T_eff = +43 K, Delta_logg = -0.02 dex, and Delta_xi = +0.10 km/s (in the sense of A - B), consistent with previous findings. Abundances of 15 elements are derived and are found to be indistinguishable between the two stars. The abundances of each element differ by <= 0.026 dex, and the mean difference is +0.003 +/- 0.015 (sigma) dex. Aside from Li, which has been previously shown to be depleted by a factor of at least 4.5 in 16 Cyg B relative to 16 Cyg A, the two stars appear to be chemically identical. The abundances of each star demonstrate a positive correlation with the condensation temperature of the elements (T_c); the slopes of the trends are also indistinguishable. In accordance with recent suggestions, the positive slopes of the [m/H]-T_c relations may imply that terrestrial planets have not formed around either 16 Cyg A or 16 Cyg B. The physical characteristics of the 16 Cyg system are discussed in terms of planet formation models, and plausible mechanisms that can account for the lack of detected planets around 16 Cyg A, the disparate Li abundances of 16 Cyg A and B, and the eccentricity of the planet 16 Cyg B b are suggested.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fe I and Fe II Abundances of Solar-Type Dwarfs in the Pleiades Open Cluster

We have derived Fe abundances of 16 solar-type Pleiades dwarfs by means of an equivalent width analysis of Fe I and Fe II lines in high-resolution spectra obtained with the Hobby - Eberly Telescope and High Resolution Spectrograph. Abundances derived from Fe II lines are larger than those derived from Fe I lines (herein referred to as over-ionization) for stars with Teff < 5400 K, and the discrepancy (deltaFe = [Fe II/H] - [Fe I/H]) increases dramatically with decreasing Teff, reaching over 0.8 dex for the coolest stars of our sample. The Pleiades joins the open clusters M 34, the Hyades, IC 2602, and IC 2391, and the Ursa Major moving group, demonstrating ostensible over-ionization trends. The Pleiades deltaFe abundances are correlated with Ca II infrared triplet and Halpha chromospheric emission indicators and relative differences therein. Oxygen abundances of our Pleiades sample derived from the high-excitation O I triplet have been previously shown to increase with decreasing Teff, and a comparison with the deltaFe abundances suggests that the over-excitation (larger abundances derived from high excitation lines relative to low excitation lines) and over-ionization effects that have been observed in cool open cluster and disk field main sequence (MS) dwarfs share a common origin. Star-to-star Fe I abundances have low internal scatter, but the abundances of stars with Teff < 5400 K are systematically higher compared to the warmer stars. The cool star [Fe I/H] abundances cannot be connected directly to over-excitation effects, but similarities with the deltaFe and O I triplet trends suggest the abundances are dubious. Using the [Fe I/H] abundances of five stars with Teff > 5400 K, we derive a mean Pleiades cluster metallicity of [Fe/H] = +0.01 +/- 0.02.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Li I and K I Scatter in Cool Pleiades Dwarfs

We utilize high-resolution (R~60,000), high S/N (~100) spectroscopy of 17 cool Pleiades dwarfs to examine the confounding star-to-star scatter in the 6707 Li I line strengths in this young cluster. Our Pleiads, selected for their small projected rotational velocity and modest chromospheric emission, evince substantial scatter in the linestrengths of 6707 Li I feature that is absent in the 7699 K I resonance line. The Li I scatter is not correlated with that in the high-excitation 7774 O I feature, and the magnitude of the former is greater than the latter despite the larger temperature sensitivity of the O I feature. These results suggest that systematic errors in linestrength measurements due to blending, color (or color-based T_eff) errors, or line formation effects related to an overlying chromosphere are not the principal source of Li I scatter in our stars. There do exist analytic spot models that can produce the observed Li scatter without introducing scatter in the K I line strengths or the color-magnitude diagram. However, these models predict factor of >3 differences in abundances derived from the subordinate 6104 and resonance 6707 Li I features; we find no difference in the abundances determined from these two features. These analytic spot models also predict CN line strengths significantly larger than we observe in our spectra. The simplest explanation of the Li, K, CN, and photometric data is that there must be a real abundance component to the Pleiades Li dispersion. We suggest that this real abundance component is the manifestation of relic differences in erstwhile pre-main-sequence Li burning caused by effects of surface activity on stellar structure. We discuss observational predictions of these effects.

0 citations0 reviewsNo code linkedOpen access