Formation of large-scale structures in ablative Kelvin-Helmholtz instability

Wang, L. F.; Ye, W. H.; Don, Wai-Sun; Sheng, Z. M.; Li, Y. J.; He, X. T.
December 2010
Physics of Plasmas;Dec2010, Vol. 17 Issue 12, p122308
Academic Journal
In this research, we studied numerically nonlinear evolutions of the Kelvin-Helmholtz instability (KHI) with and without thermal conduction, aka, the ablative KHI (AKHI) and the classical KHI (CKHI). The second order thermal conduction term with a variable thermal conductivity coefficient is added to the energy equation in the Euler equations in the AKHI to investigate the effect of thermal conduction on the evolution of large and small scale structures within the shear layer which separate the fluids with different velocities. The inviscid hyperbolic flux of Euler equation is computed via the classical fifth order weighted essentially nonoscillatory finite difference scheme and the temperature is solved by an implicit fourth order finite difference scheme with variable coefficients in the second order parabolic term to avoid severe time step restriction imposed by the stability of the numerical scheme. As opposed to the CKHI, fine scale structures such as the vortical structures are suppressed from forming in the AKHI due to the dissipative nature of the second order thermal conduction term. With a single-mode sinusoidal interface perturbation, the results of simulations show that the growth of higher harmonics is effectively suppressed and the flow is stabilized by the thermal conduction. With a two-mode sinusoidal interface perturbation, the vortex pairing is strengthened by the thermal conduction which would allow the formation of large-scale structures and enhance the mixing of materials. In summary, our numerical studies show that thermal conduction can have strong influence on the nonlinear evolutions of the KHI. Thus, it should be included in applications where thermal conduction plays an important role, such as the formation of large-scale structures in the high energy density physics and astrophysics.


Related Articles

  • Preheating ablation effects on the Rayleigh-Taylor instability in the weakly nonlinear regime. Wang, L. F.; Ye, W. H.; Sheng, Z. M.; Don, Wai-Sun; Li, Y. J.; He, X. T. // Physics of Plasmas;Dec2010, Vol. 17 Issue 12, p122706 

    The two-dimensional Rayleigh-Taylor instability (RTI) with and without thermal conduction is investigated by numerical simulation in the weakly nonlinear regime. A preheat model κ(T)=κSH[1+f(T)] is introduced for the thermal conduction [W. H. Ye, W. Y. Zhang, and X. T. He, Phys. Rev. E 65,...

  • Efficient evaluation of analytic Fukui functions. Flores-Moreno, Roberto; Melin, Junia; Ortiz, J. V.; Merino, Gabriel // Journal of Chemical Physics;12/14/2008, Vol. 129 Issue 22, p224105 

    An efficient method for the analytic evaluation of Fukui functions is proposed. Working equations are derived and numerical results are used to validate the method on medium size set of molecules. In addition to the obvious advantages of analytic differentiation, the proposed method is efficient...

  • Spike deceleration and bubble acceleration in the ablative Rayleigh-Taylor instability. Ye, W. H.; Wang, L. F.; He, X. T. // Physics of Plasmas;Dec2010, Vol. 17 Issue 12, p122704 

    The nonlinear evolutions of the Rayleigh-Taylor instability (RTI) with preheat is investigated by numerical simulation (NS). A new phase of the spike deceleration evolution in the nonlinear ablative RTI (ARTI) is discovered. It is found that nonlinear evolution of the RTI can be divided into the...

  • Sensitivity analysis via simulation in the presence of discontinuities. Signahl, Mikael // Mathematical Methods of Operations Research;2004, Vol. 60 Issue 1, p29 

    In this paper we address the problem of estimating the mean derivative when the entity containing the parameter has jumps. The methods considered are finite differences, infinitesimal perturbation analysis and the likelihood ratio score function. We calculate the difference between the...

  • A New System of Multivalued Mixed Variational Inequality Problem. Xi Li; Xue-song Li // Abstract & Applied Analysis;2014, p1 

    We consider a new system of multivalued mixed variational inequality problem, which includes some known systems of variational inequalities as special cases. Under suitable conditions, the existence of solutions for the system of multivalued mixed variational inequality problem and the...

  • Numerical Investigation of the Nonlinear Schrödinger Equation with Saturation. Christou, M. A. // AIP Conference Proceedings;10/30/2008, Vol. 1067 Issue 1, p105 

    We construct an implicit finite difference scheme to investigate numerically the nonlinear Schrödinger equation with saturation (NLSS). We use Newton’s method to linearize the numerical scheme. We examine the propagation, interaction and overtake interaction of soliton solutions of the...

  • Propagation of nonlinear dust magnetoacoustic waves in cylindrical geometry. Hussain, S.; Mahmood, S. // Physics of Plasmas;Dec2011, Vol. 18 Issue 12, p123701 

    The cylindrical dust magnetoacoustic shocks in dissipative and solitons in non-dissipative dust-ion plasmas are studied. The dust particles are assumed to be fully negatively charged so that the density of electrons in comparison with dust and ions is negligible. The cylindrical Korteweg-de...

  • A UNIFORM NUMERICAL METHOD FOR A BOUNDARY-SHOCK PROBLEM. Vulanović, Relja // International Journal of Numerical Analysis & Modeling;2010, Vol. 7 Issue 3, p567 

    A singularly perturbed quasilinear boundary-value problem is considered in the case when its solution has a boundary shock. The problem is discretized by an upwind finite-difference scheme on a mesh of Shishkin type. It is proved that this numerical method has pointwise accuracy of almost first...

  • Parameter Iteration Method for Solving Nonlinear Problem. Lin, Jiang-guo // Applied Mathematics & Mechanics;Jul2001, Vol. 22 Issue 7, p836 

    The parameter iteration method was developed for solving the nonlinear problem in this paper. The results of several examples show that even the first iteration solution has very good accuracy.


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics