Pulsating waves in nonlinear magnetoconvection

Numerical experiments on compressible magnetoconvection reveal a new type of periodic oscillation, associated with alternating streaming motion. Analogous behaviour in a Boussinesq fluid is constrained by extra symmetry. A low-order model confirms that these pulsating waves appear via a pitchfork-Hopf-gluing bifurcation sequence from the steady state

Adding kinetics and hydrodynamics to the CHEETAH thermochemical code

In FY96 we released CHEETAH 1.40, which made extensive improvements on the stability and user friendliness of the code. CHEETAH now has over 175 users in government, academia, and industry. Efforts have also been focused on adding new advanced features to CHEETAH 2.0, which is scheduled for release in FY97. We have added a new chemical kinetics capability to CHEETAH. In the past, CHEETAH assumed complete thermodynamic equilibrium and independence of time. The addition of a chemical kinetic framework will allow for modeling of time-dependent phenomena, such as partial combustion and detonation in composite explosives with large reaction zones. We have implemented a Wood-Kirkwood detonation framework in CHEETAH, which allows for the treatment of nonideal detonations and explosive failure. A second major effort in the project this year has been linking CHEETAH to hydrodynamic codes to yield an improved HE product equation of state. We have linked CHEETAH to 1- and 2-D hydrodynamic codes, and have compared the code to experimental data. 15 refs., 13 figs., 1 tab

Solar Stereoscopy with STEREO/EUVI A and B spacecraft from small (6 deg) to large (170 deg) spacecraft separation angles

We performed for the first time stereoscopic triangulation of coronal loops in active regions over the entire range of spacecraft separation angles ($\alpha_{sep}\approx 6^\circ, 43^\circ, 89^\circ, 127^\circ$, and $170^\circ$). The accuracy of stereoscopic correlation depends mostly on the viewing angle with respect to the solar surface for each spacecraft, which affects the stereoscopic correspondence identification of loops in image pairs. From a simple theoretical model we predict an optimum range of $\alpha_{sep} \approx 22^\circ-125^\circ$, which is also experimentally confirmed. The best accuracy is generally obtained when an active region passes the central meridian (viewed from Earth), which yields a symmetric view for both STEREO spacecraft and causes minimum horizontal foreshortening. For the extended angular range of $\alpha_{sep}\approx 6^\circ-127^{\circ}$ we find a mean 3D misalignment angle of $\mu_{PF} \approx 21^\circ-39^\circ$ of stereoscopically triangulated loops with magnetic potential field models, and $\mu_{FFF} \approx 15^\circ-21^\circ$ for a force-free field model, which is partly caused by stereoscopic uncertainties $\mu_{SE} \approx 9^\circ$. We predict optimum conditions for solar stereoscopy during the time intervals of 2012--2014, 2016--2017, and 2021--2023.Comment: Solar Physics, (in press), 22 pages, 9 figure

Kang-Redner Anomaly in Cluster-Cluster Aggregation

The large time, small mass, asymptotic behavior of the average mass distribution \pb is studied in a $d$-dimensional system of diffusing aggregating particles for $1\leq d \leq 2$. By means of both a renormalization group computation as well as a direct re-summation of leading terms in the small reaction-rate expansion of the average mass distribution, it is shown that \pb \sim \frac{1}{t^d} (\frac{m^{1/d}}{\sqrt{t}})^{e_{KR}} for $m \ll t^{d/2}$, where $e_{KR}=\epsilon +O(\epsilon ^2)$ and $\epsilon =2-d$. In two dimensions, it is shown that \pb \sim \frac{\ln(m) \ln(t)}{t^2} for $m \ll t/ \ln(t)$. Numerical simulations in two dimensions supporting the analytical results are also presented.Comment: 11 pages, 6 figures, Revtex

Random Walks in Logarithmic and Power-Law Potentials, Nonuniversal Persistence, and Vortex Dynamics in the Two-Dimensional XY Model

The Langevin equation for a particle (`random walker') moving in d-dimensional space under an attractive central force, and driven by a Gaussian white noise, is considered for the case of a power-law force, F(r) = - Ar^{-sigma}. The `persistence probability', P_0(t), that the particle has not visited the origin up to time t, is calculated. For sigma > 1, the force is asymptotically irrelevant (with respect to the noise), and the asymptotics of P_0(t) are those of a free random walker. For sigma < 1, the noise is (dangerously) irrelevant and the asymptotics of P_0(t) can be extracted from a weak noise limit within a path-integral formalism. For the case sigma=1, corresponding to a logarithmic potential, the noise is exactly marginal. In this case, P_0(t) decays as a power-law, P_0(t) \sim t^{-theta}, with an exponent theta that depends continuously on the ratio of the strength of the potential to the strength of the noise. This case, with d=2, is relevant to the annihilation dynamics of a vortex-antivortex pair in the two-dimensional XY model. Although the noise is multiplicative in the latter case, the relevant Langevin equation can be transformed to the standard form discussed in the first part of the paper. The mean annihilation time for a pair initially separated by r is given by t(r) \sim r^2 ln(r/a) where a is a microscopic cut-off (the vortex core size). Implications for the nonequilibrium critical dynamics of the system are discussed and compared to numerical simulation results.Comment: 10 pages, 1 figur

Signatures of the slow solar wind streams from active regions in the inner corona

Some of local sources of the slow solar wind can be associated with spectroscopically detected plasma outflows at edges of active regions accompanied with specific signatures in the inner corona. The EUV telescopes (e.g. SPIRIT/CORONAS-F, TESIS/CORONAS-Photon and SWAP/PROBA2) sometimes observed extended ray-like structures seen at the limb above active regions in 1MK iron emission lines and described as "coronal rays". To verify the relationship between coronal rays and plasma outflows, we analyze an isolated active region (AR) adjacent to small coronal hole (CH) observed by different EUV instruments in the end of July - beginning of August 2009. On August 1 EIS revealed in the AR two compact outflows with the Doppler velocities V =10-30 km/s accompanied with fan loops diverging from their regions. At the limb the ARCH interface region produced coronal rays observed by EUVI/STEREO-A on July 31 as well as by TESIS on August 7. The rays were co-aligned with open magnetic field lines expanded to the streamer stalks. Using the DEM analysis, it was found that the fan loops diverged from the outflow regions had the dominant temperature of ~1 MK, which is similar to that of the outgoing plasma streams. Parameters of the solar wind measured by STEREO-B, ACE, WIND, STEREO-A were conformed with identification of the ARCH as a source region at the Wang-Sheeley-Arge map of derived coronal holes for CR 2086. The results of the study support the suggestion that coronal rays can represent signatures of outflows from ARs propagating in the inner corona along open field lines into the heliosphere.Comment: Accepted for publication in Solar Physics; 31 Pages; 13 Figure

On the structure and evolution of a polar crown prominence/filament system

Polar crown prominences are made of chromospheric plasma partially circling the Suns poles between 60 and 70 degree latitude. We aim to diagnose the 3D dynamics of a polar crown prominence using high cadence EUV images from the Solar Dynamics Observatory (SDO)/AIA at 304 and 171A and the Ahead spacecraft of the Solar Terrestrial Relations Observatory (STEREO-A)/EUVI at 195A. Using time series across specific structures we compare flows across the disk in 195A with the prominence dynamics seen on the limb. The densest prominence material forms vertical columns which are separated by many tens of Mm and connected by dynamic bridges of plasma that are clearly visible in 304/171A two-color images. We also observe intermittent but repetitious flows with velocity 15 km/s in the prominence that appear to be associated with EUV bright points on the solar disk. The boundary between the prominence and the overlying cavity appears as a sharp edge. We discuss the structure of the coronal cavity seen both above and around the prominence. SDO/HMI and GONG magnetograms are used to infer the underlying magnetic topology. The evolution and structure of the prominence with respect to the magnetic field seems to agree with the filament linkage model.Comment: 24 pages, 14 figures, Accepted for publication in Solar Physics Journal, Movies can be found at http://www2.mps.mpg.de/data/outgoing/panesar

Non-monotonic variation with salt concentration of the second virial coefficient in protein solutions

The osmotic virial coefficient $B_2$ of globular protein solutions is calculated as a function of added salt concentration at fixed pH by computer simulations of the ``primitive model''. The salt and counter-ions as well as a discrete charge pattern on the protein surface are explicitly incorporated. For parameters roughly corresponding to lysozyme, we find that $B_2$ first decreases with added salt concentration up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength. Our studies demonstrate that the existence of a discrete charge pattern on the protein surface profoundly influences the effective interactions and that non-linear Poisson Boltzmann and Derjaguin-Landau-Verwey-Overbeek (DLVO) theory fail for large ionic strength. The observed non-monotonicity of $B_2$ is compared to experiments. Implications for protein crystallization are discussed.Comment: 43 pages, including 17 figure