• In this paper the problem is posed of the prescription of the so-called Boltzmann-Grad (BG) limit ($\mathcal{L}_{BG}$) for the $N-$body system of smooth hard-spheres which undergo unary, binary as well as multiple elastic instantaneous collisions. The statistical description is couched in terms of the Master kinetic equation, i.e., the kinetic equation which realizes the axiomatic "\textit{ab initio}" approach to the classical statistical mechanics of finite hard-sphere systems recently developed (Tessarotto \textit{et al.}, 2013-2017). The issue addressed here concerns the prescription of the BG-limit operator and specifically the non-commutative property of $\mathcal{L}_{BG}$ with the free-streaming operator which enters the same kinetic equation.\ It is shown\ that the form of the resulting limit equation remains in principle non-unique, its precise realization depending critically on the way the action of the same operator is prescribed. Implications for the global prescription of the Boltzmann equation are pointed out.
  • Differentially rotating disks under the effect of axial magnetic field are prone to a nonlinear explosive magnetorotational instability (EMRI). The dynamic equations that govern the temporal evolution of the amplitudes of three weakly-detuned resonantly interacting modes are derived. As distinct from exponential growth in the strict resonance triads EMRI occurs due to the resonant interactions of a MRI mode with stable Alfv\'en-Coriolis and magnetosonic modes. Numerical solutions of the dynamic equations for amplitudes of a triad indicate that two types of perturbations behavior can be excited for resonance conditions: (i) EMRI which leads to infinite values of the three amplitudes within a finite time, and (ii) bounded irregular oscillations of all three amplitudes. Asymptotic explicit solutions of the dynamic equations are obtained for EMRI regimes and are shown to match the numerical solutions near the explosion time.
  • The linear instability of thin, vertically-isothermal Keplerian discs, under the influence of axial magnetic field is investigated. Solutions of the stability problem are found explicitly by asymptotic expansions in the small aspect ratio of the disc. It is shown that the perturbations are decoupled into in-plane and vertical modes. Exact expressions for the growth rates as well as the number of unstable modes are derived. Those are the discrete counterpart of the continuous infinite homogeneous cylinder magnetorotational (MRI) spectrum. In addition, a weakly nonlinear analysis of the MRI is performed. It is shown that near the instability threshold the latter is saturated by the stable magnetoacoustic modes.
  • The stability of nonlinear explicit difference schemes with not, in general, open domains of the scheme operators are studied. For the case of path-connected, bounded, and Lipschitz domains, we establish the notion that a multi-level nonlinear explicit scheme is stable iff (if and only if) the corresponding scheme in variations is stable. A new modification of the central Lax-Friedrichs (LxF) scheme is developed to be of the second order accuracy. The modified scheme is based on nonstaggered grids. A monotone piecewise cubic interpolation is used in the central scheme to give an accurate approximation for the model in question. The stability of the modified scheme is investigated. Some versions of the modified scheme are tested on several conservation laws, and the scheme is found to be accurate and robust. As applied to hyperbolic conservation laws with, in general, stiff source terms, it is constructed a second order nonstaggered central scheme based on operator-splitting techniques.
  • The present study concerns the numerical modeling of sea-wave instability under the effect of logarithmic-wind profile in hurricane conditions. The central point of the study is the calculation of the wave growth rate, which is proportional to the fractional input energy from the weakly-sheared (logarithmic) wind to the wave exponentially varying with time. It is shown for hurricane conditions that the Miles-type stability model based on the Charnock's formula with the standard constant coefficient underestimates the growth rate ~5 to 50 times as compared with the model employing the roughness adopted from experimental data for hurricane winds. The drag reduction with wind speed at hurricane conditions coupled with the similar behavior of the dimensionless gravity acceleration, leads to the minimum in the maximal growth rate and the maximum in the most unstable wavelength.
  • The present study concerns the numerical modeling of sea-wave instability under the effect of logarithmic-wind profile in hurricane conditions. The central point of the study is the calculation of the wave growth rate, which is proportional to the fractional input energy from the weakly-sheared (logarithmic) wind to the wave exponentially varying with time. It is shown for hurricane conditions that the Miles-type stability model based on the Charnock's formula with the standard constant coefficient underestimates the growth rate ~5 to 50 times as compared with the model employing the roughness adopted from experimental data for hurricane winds. The drag reduction with wind speed at hurricane conditions coupled with the similar behavior of the dimensionless gravity acceleration, leads to the minimum in the maximal growth rate and the maximum in the most unstable wavelength.
  • The short time evolution of three dimensional small perturbations is studied. Exhibiting spectral asymptotic stability, thin discs are nonetheless shown to host intensive hydrodynamical activity in the shape of non modal growth of initial small perturbations. Two mechanisms that lead to such behavior are identified and studied, namely, non-resonant excitation of vertically confined sound waves by stable planar inertia-coriolis modes that results in linear growth with time, as well as resonant coupling of those two modes that leads to a quadratic growth of the initial perturbations. It is further speculated that the non modal growth can give rise to secondary strato-rotational instabilities and thus lead to a new route to turbulence generation in thin discs.
  • The stability of difference schemes for, in general, hyperbolic systems of conservation laws with source terms are studied. The basic approach is to investigate the stability of a non-linear scheme in terms of its cor- responding scheme in variations. Such an approach leads to application of the stability theory for linear equation systems to establish stability of the corresponding non-linear scheme. It is established the notion that a non-linear scheme is stable if and only if the corresponding scheme in variations is stable. A new modification of the central Lax-Friedrichs (LxF) scheme is developed to be of the second order accuracy. A monotone piecewise cubic interpolation is used in the central schemes to give an accurate approximation for the model in question. The stability of the modified scheme are investigated. Some versions of the modified scheme are tested on several conservation laws, and the scheme is found to be accurate and robust. As applied to hyperbolic conservation laws with, in general, stiff source terms, it is constructed a second order nonstaggered central scheme based on operator-splitting techniques.
  • The monotonicity and stability of difference schemes for, in general, hyperbolic systems of conservation laws with source terms are studied. The basic approach is to investigate the stability and monotonicity of a non-linear scheme in terms of its corresponding scheme in variations. Such an approach leads to application of the stability theory for linear equation systems to establish stability of the corresponding non-linear scheme. The main methodological innovation is the theorems establishing the notion that a non-linear scheme is stable (and monotone) if the corresponding scheme in variations is stable (and, respectively, monotone). Criteria are developed for monotonicity and stability of difference schemes associated with the numerical analysis of systems of partial differential equations. The theorem of Friedrichs (1954) is generalized to be applicable to variational schemes with non-symmetric matrices. A new modification of the central Lax-Friedrichs (LxF) scheme is developed to be of the second order accuracy. A monotone piecewise cubic interpolation is used in the central schemes to give an accurate approximation for the model in question. The stability and monotonicity of the modified scheme are investigated. Some versions of the modified scheme are tested on several conservation laws, and the scheme is found to be accurate and robust. As applied to hyperbolic conservation laws with, in general, stiff source terms, it is constructed a second order scheme based on operator-splitting techniques.
  • An important aspect of computational fluid dynamics is related to the determination of the fluid pressure in isothermal incompressible fluids. In particular this concerns the construction of an exact evolution equation for the fluid pressure which replaces the Poisson equation and yields an algorithm which is a Poisson solver, i.e., it permits to time-advance exactly the same fluid pressure \textit{without solving the Poisson equation}% . In fact, the incompressible Navier-Stokes equations represent a mixture of hyperbolic and elliptic pde's, which are extremely hard to study both analytically and numerically. In this paper we intend to show that an exact solution to this problem can be achieved adopting the approach based on inverse kinetic theory (IKT) recently developed for incompressible fluids by Ellero and Tessarotto (2004-2007). In particular we intend to prove that the evolution of the fluid fields can be achieved by means of a suitable dynamical system, to be identified with the so-called Navier-Stokes (N-S) dynamical system. As a consequence it is found that the fluid pressure obeys a well-defined evolution equation. The result appears relevant for the construction of Lagrangian approaches to fluid dynamics.
  • Cool weakly ionized gaseous rotating disk, are considered by many models as the origin of the evolution of protoplanetary clouds. Instabilities against perturbations in such disks play an important role in the theory of the formation of stars and planets. Thus, a hierarchy of successive fragmentations into smaller and smaller pieces as a part of the Kant-Laplace theory of formation of the planetary system remains valid also for contemporary cosmogony. Traditionally, axisymmetric magnetohydrodynamic (MHD), and recently Hall-MHD instabilities have been thoroughly studied as providers of an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, the Hall instability against nonaxisymmetric perturbations in compressible rotating fluids in external magnetic field is proposed as a viable mechanism for the azimuthal fragmentation of the protoplanetary disk and thus perhaps initiating the road to planet formation. The Hall instability is excited due to the combined effect of the radial stratification of the disk and the Hall electric field, and its growth rate is of the order of the rotation period.
  • The effect of particles that undergo strong diffusive-shock-acceleration on the stability of the accelerating shock is investigated. A two-fluid model is employed in which the accelerated particles are treated as a fluid whose effect is incorporated as an additional pressure in the momentum equation. The Dyakov and Kontorovich stability criteria are used in order to study the stability of those steady-state shocks that contain a gas sub-shock. The downstream conditions of the latter are parametrized by the ratio of the upstream acclerated-particles pressure to the total pressure. For some range of values of that parameter, three possible downstream states are possible for each upstream state. It is shown that in that range of parameters the shocks are either corrugationally unstable or lose their energy by spontaneous emission of acoustic as well as entropy-vortex waves.