Layering and vertical transport in sheared double-diffusive convection in the diffusive regime

A sequence of two- and three-dimensional simulations are conducted for the double-diffusive convection (DDC) flows in the diffusive regime subjected to an imposed shear. For a wide range of control parameters, and for sufficiently strong perturbation of the conductive initial state, staircase-like structures spontaneously develop, with relatively well-mixed layers separated by sharp interfaces of enhanced scalar gradient. Such staircases appear to be robust even in the presence of strong shear over very long times, with early-time coarsening of the observed layers. For the same set of control parameters, different asymptotic layered states, with markedly different vertical scalar fluxes, can arise for different initial perturbation structures. The imposed shear significantly spatio-temporally modifies the vertical transport of the various scalars. The flux ratio $\gamma ^*$ (i.e. the ratio between the density fluxes due to the total salt flux and the total heat flux) is found, at steady state, to be essentially equal to the square root of the ratio of the salt diffusivity to the thermal diffusivity, consistent with the physical model proposed by Linden & Shirtcliffe (J. Fluid Mech., vol. 87, 1978, pp. 417–432) and the variational arguments presented by Stern (J. Fluid Mech., vol. 114, 1982, pp. 105–121) for unsheared double-diffusive convection.

Fractional order plasma modeling based on linear polarization of LASER light: an Atangana–Baleanu Caputo approach

PurposeThe purpose of this offered research is to articulate a multifaceted kind of highly unstable initial perturbation and further analyze the performance of the plasma particles for time-fractional order evaluation.Design/methodology/approachFor this purpose, the authors designed specific geometry and further interpreted it into the mathematical model using the concepts of the Vlasov Maxwell system. The suggested algorithm is based on the finite-difference and spectral estimation philosophy. The management of time and memory in generic code for computational purposes is also discussed.FindingsThe main purpose is to analyze the fractional behavior of plasma particles and also the capability of the suggested numerical algorithm. Due to initial perturbations, there are a lot of sudden variations that occurred in the formulated system. Graphical behavior shows that SR parameter produces devastation as compared to others. The variation of fractional parameter between the defend domain demonstrates the hidden pictures of plasma particles. The design scheme is efficient, convergent and has the capability to cover the better physics of the problem.Practical implicationsPlasma material is commonly used in different areas of science. Therefore, in this paper, the authors increase the capability of the mathematical plasma model with specific geometry, and further suitable numerical algorithm is suggested with detailed physical analysis of the outcomes. The authors gave a new direction to study the performance of plasma particles under the influence of LASER light.Originality/valueIn the recent era, science has produced a lot of advancements to study and analyze the physical natural process, which exist everywhere in the real word. On behalf of this current developments, it is now insufficient to study the first-order time evaluation of the plasma particles. One needs to be more precise and should move toward the bottomless state of it, that is, macroscopic and microscopic time-evaluation scales, and it is not wrong to say that there exits a huge gap, to study the time evaluation in this discussed manner. The presented study is entirely an advanced and efficient way to investigate the problem into the new directions. The capability of the proposed algorithm and model with fractional concepts can fascinate the reader to extend to the other dimensions.

Shock-induced interfacial instabilities of granular media

This paper investigates the shock-induced instability of the interfaces between gases and dense granular media with finite length via the coarse-grained compressible computational fluid dynamics–discrete parcel method. Despite generating a typical spike-bubble structure reminiscent of the Richtmyer–Meshkov instability (RMI), the shock-driven granular instability (SDGI) is governed by fundamentally different mechanisms. Unlike the RMI arising from baroclinic vorticity deposition on the interface, the SDGI is closely associated with the interfacial and bulk granular dynamics, which evolve with the transient coupling between particles and gases. Consequently, the SDGI follows a growth law distinctly different from that of the RMI, namely a semilinear slow regime followed by an exponentially expedited regime and a quadratic asymptotic regime. We further establish the instability criteria of the SDGI for granular media with infinite and finite lengths, which do not exist in the RMI. A scaling growth law of the SDGI for dense granular media with finite length is derived by normalizing the time with the rarefaction propagation time, which successfully collapses the data from cases with varying shock strength, particle column length and particle volume fraction and ought to hold for granular media with varying particle parameters. The effect of the initial perturbation magnitude can be properly considered in the scaling growth law by incorporating it into the length normalization.

Glacier geometry and flow speed determine how Arctic marine-terminating glaciers respond to lubricated beds

Basal conditions directly control the glacier sliding rate and the dynamic discharge of ice flow. Recent glacier destabilization events indicate that some marine-terminating glaciers quickly respond to lubricated beds with increased flow speed, but the underlying physics, especially how this vulnerability relates to glacier geometry and flow characteristics, remains unclear. This paper presents a 1-D physical framework for glacier dynamic vulnerability assuming sudden basal lubrication as an initial perturbation. In this new model, two quantities determine the scale and the areal extent of the subsequent thinning and acceleration after the bed is lubricated: Péclet number (Pe) and the product of glacier speed and thickness gradient (dubbed J0 in this study). To validate the model, this paper calculates Pe and J0 using multi-sourced data from 1996–1998 for outlet glaciers in Greenland and Austfonna Ice Cap, Svalbard, and compares the results with the glacier speed change during 1996/1998–2018. Glaciers with lower Pe and J0 are more likely to accelerate during this 20-year span than those with higher Pe and J0, which matches the model prediction. A combined factor of ice thickness, surface slope, and initial speed for ice flow physically determines how much and how fast glaciers respond to lubricated beds, as forms of speed, elevation, and terminus change.

Using Conditional Nonlinear Optimal Perturbation to Generate Initial Perturbations in ENSO Ensemble Forecasts

Using the latest operational version of the ENSO forecast system from the National Marine Environmental Forecasting Center (NMEFC) of China, ensemble forecasting experiments are performed for El Niño-Southern Oscillation (ENSO) events that occurred from 1997 to 2017 by generating initial perturbations of the conditional nonlinear optimal perturbation (CNOP) and Climatically relevant Singular Vector (CSV) structures. It is shown that when the initial perturbation of the leading CSV structure in the ensemble forecast of the CSVs-scheme is replaced by those of the CNOP structure, the resulted ensemble ENSO forecasts of the CNOP+CSVs-scheme tend to possess a larger spread than the forecasts obtained with the CSVs-scheme alone, leading to a better match between the root mean square error and the ensemble spread, a more reasonable Talagrand diagram and an improved Brier skill score (BSS). All these results indicate that the ensemble forecasts generated by the CNOP+CSVs-scheme can improve both the accuracy of ENSO forecasting and the reliability of the ensemble forecasting system. Therefore, ENSO ensemble forecasting should consider the effect of nonlinearity on the ensemble initial perturbations to achieve a much higher skill. It is expected that fully nonlinear ensemble initial perturbations can be sufficiently yielded to produce ensemble forecasts for ENSO, finally improving the ENSO forecast skill to the greatest possible extent. The CNOP will be a useful method to yield fully nonlinear optimal initial perturbations for ensemble forecasting.

Statistical characteristics of turbulent mixing in spherical and cylindrical converging Richtmyer–Meshkov instabilities

In this paper, the Richtmyer–Meshkov instabilities in spherical and cylindrical converging geometries with a Mach number of approximately 1.5 are investigated by using the high resolution implicit large eddy simulation method, and the influence of the geometric effect on the turbulent mixing is investigated. The heavy fluid is sulphur hexafluoride (SF6), and the light fluid is nitrogen (N2). The shock wave converges from the heavy fluid into the light fluid. The Atwood number is 0.678. The total structured and uniform Cartesian grid node number in the main computational domain is 20483. In addition, to avoid the influence of boundary reflection, a sufficiently long sponge layer with 50 non-uniform coarse grids is added for each non-periodic boundary. Present numerical simulations have high and nonlinear initial perturbation levels, which rapidly lead to turbulent mixing in the mixing layers. Firstly, some physical-variable mean profiles, including mass fraction, Taylor Reynolds number, turbulent kinetic energy, enstrophy and helicity, are provided. Second, the mixing characteristics in the spherical and cylindrical turbulent mixing layers are investigated, such as molecular mixing fraction, efficiency Atwood number, turbulent mass-flux velocity and density self-correlation. Then, Reynolds stress and anisotropy are also investigated. Finally, the radial velocity, velocity divergence and enstrophy in the spherical and cylindrical turbulent mixing layers are studied using the method of conditional statistical analysis. Present numerical results show that the geometric effect has a great influence on the converging Richtmyer–Meshkov instability mixing layers.

Universal self-similar attractor in the bending-driven levelling of thin viscous films

We study theoretically and numerically the bending-driven levelling of thin viscous films within the lubrication approximation. We derive Green’s function of the linearized thin-film equation and further show that it represents a universal self-similar attractor at long times. As such, the rescaled perturbation of the film profile converges in time towards the rescaled Green’s function, for any summable initial perturbation profile. In addition, for stepped axisymmetric initial conditions, we demonstrate the existence of another, short-term and one-dimensional-like self-similar regime. We also characterize the convergence time towards the long-term universal attractor in terms of the relevant physical and geometrical parameters, and provide the local hydrodynamic fields and global elastic energy in the universal regime as functions of time. Finally, we extend our analysis to the nonlinear thin-film equation through numerical simulations.

New Types of Dissipative Streaming Instabilities

Two new, previously unknown types of dissipative streaming instabilities (DSI) are substantiated. They follow from new approach, which allows solving in general form the classical problem of an initial perturbation development for streaming instabilities (SI). SI is caused by relative motion of the streams of plasma components. With an increase in level of dissipation SI transforms into a DSI. The transformation occurs because dissipation serves as a channel for energy removal for the growth of the negative energy wave of the stream. Until recently, only one type of DSI was known. Its maximal growth rate depends on the beam density nb and the collision frequency ν in the plasma as ∼nb/ν. All types of conventional beam-plasma instabilities (Cherenkov, cyclotron, etc.) transform into it. The solution of the problem of the initial perturbation development in systems with weak beam-plasma coupling leads to a new type of DSI. With an increase in the level of dissipation, the instability in these systems transforms to the new DSI. Its maximal growth rate is ∼nb/ν. The second new DSI develops in beam-plasma waveguide with over-limiting current of e-beam. Its growth rate ∼nb/ν. In addition, the solutions of abovementioned problem provide much information about SI and DSI, significant part of which is unavailable by other methods.