No Arabic abstract
Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently rapidly the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one and three dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.
We use the Fokker-Planck equation and model the dispersive dynamics of solid particles in annular protoplanetary disks whose gas component is more massive than the particle phase. We model particle--gas interactions as hard sphere collisions, determine the functional form of diffusion coefficients, and show the existence of two global unstable modes in the particle phase. These modes have spiral patterns with the azimuthal wavenumber $m=1$ and rotate slowly. We show that in ring-shaped disks, the phase space density of solid particles increases linearly in time towards an accumulation point near the location of pressure maximum, while instabilities grow exponentially. Therefore, planetesimals and planetary cores can be efficiently produced near the peaks of unstable density waves. In this mechanism, particles migrating towards the accumulation point will not participate in the formation of planets, and should eventually form a debris ring like the main asteroid belt or classical Kuiper belt objects. We present the implications of global instabilities to the formation of ice giants and terrestrial planets in the solar system.
Aims. We show how the build-up of magnetic gradients in the Suns corona may be inferred directly from photospheric velocity data. This enables computation of magnetic connectivity measures such as the squashing factor without recourse to magnetic field extrapolation. Methods.Assuming an ideal evolution in the corona, and an initially uniform magnetic field, the subsequent field line mapping is computed by integrating trajectories of the (time-dependent) horizontal photospheric velocity field. The method is applied to a 12 hour high-resolution sequence of photospheric flows derived from Hinode/SOT magnetograms. Results. We find the generation of a network of quasi-separatrix layers in the magnetic field, which correspond to Lagrangian coherent structures in the photospheric velocity. The visual pattern of these structures arises primarily from the diverging part of the photospheric flow, hiding the effect of the rotational flow component: this is demonstrated by a simple analytical model of photospheric convection. We separate the diverging and rotational components from the observed flow and show qualitative agreement with purely diverging and rotational models respectively. Increasing the flow speeds in the model suggests that our observational results are likely to give a lower bound for the rate at which magnetic gradients are built up by real photospheric flows. Finally, we construct a hypothetical magnetic field with the inferred topology, that can be used for future investigations of reconnection and energy release.
Complex mixing and magnetic field generation occurs within stellar interiors particularly where there is a strong shear flow. To obtain a comprehensive understanding of these processes, it is necessary to study the complex dynamics of shear regions. Due to current observational limitations, it is necessary to investigate the inevitable small-scale dynamics via numerical calculations. Here, we examine direct numerical calculations of a local model of unstable shear flows in a compressible polytropic fluid primarily in a two-dimensional domain, where we focus on determining how key parameters affect the global properties and characteristics of the resulting saturated turbulent phase. We consider the effect of varying both the viscosity and the thermal diffusivity on the non-linear evolution. Moreover, our main focus is to understand the global properties of the saturated phase, in particular estimating for the first time the spread of the shear region from an initially hyperbolic tangent velocity profile. We find that the vertical extent of the mixing region in the saturated regime is generally determined by the initial Richardson number of the system. Further, the characteristic quantities of the turbulence, i.e. typical length-scale and the root-mean-square velocity are found to depend on both the Richardson number, and the thermal diffusivity. Finally, we present our findings of our investigation into saturated flows of a `secular shear instability in the low Peclet number regime with large Richardson numbers.
A reactive fluid dissolving the surrounding rock matrix can trigger an instability in the dissolution front, leading to spontaneous formation of pronounced channels or wormholes. Theoretical investigations of this instability have typically focused on a steadily propagating dissolution front that separates regions of high and low porosity. In this paper we show that this is not the only possible dissolutional instability in porous rocks; there is another instability that operates instantaneously on any initial porosity field, including an entirely uniform one. The relative importance of the two mechanisms depends on the ratio of the porosity increase to the initial porosity. We show that the inlet instability is likely to be important in limestone formations where the initial porosity is small and there is the possibility of a large increase in permeability. In quartz-rich sandstones, where the proportion of easily soluble material (e.g. carbonate cements) is small, the instability in the steady-state equations is dominant.
A dynamic mitigation mechanism for instability growth was proposed and discussed in the paper [Phys. Plasmas 19, 024503 (2012)]. In the present paper the robustness of the dynamic instability mitigation mechanism is discussed further. The results presented here show that the mechanism of the dynamic instability mitigation is rather robust against changes in the phase, the amplitude and the wavelength of the wobbling perturbation applied. Generally instability would emerge from the perturbation of the physical quantity. Normally the perturbation phase is unknown so that the instability growth rate is discussed. However, if the perturbation phase is known, the instability growth can be controlled by a superposition of perturbations imposed actively: if the perturbation is induced by, for example, a driving beam axis oscillation or wobbling, the perturbation phase could be controlled and the instability growth is mitigated by the superposition of the growing perturbations.