
In this paper the problem is posed of the prescription of the socalled
BoltzmannGrad (BG) limit ($\mathcal{L}_{BG}$) for the $N$body system of
smooth hardspheres 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 hardsphere systems recently developed (Tessarotto \textit{et al.},
20132017). The issue addressed here concerns the prescription of the BGlimit
operator and specifically the noncommutative property of $\mathcal{L}_{BG}$
with the freestreaming operator which enters the same kinetic equation.\ It is
shown\ that the form of the resulting limit equation remains in principle
nonunique, 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 weaklydetuned 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\'enCoriolis 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, verticallyisothermal 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
inplane 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
pathconnected, bounded, and Lipschitz domains, we establish the notion that a
multilevel nonlinear explicit scheme is stable iff (if and only if) the
corresponding scheme in variations is stable. A new modification of the central
LaxFriedrichs (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
operatorsplitting techniques.

The present study concerns the numerical modeling of seawave instability
under the effect of logarithmicwind 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 weaklysheared
(logarithmic) wind to the wave exponentially varying with time. It is shown for
hurricane conditions that the Milestype 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 seawave instability
under the effect of logarithmicwind 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 weaklysheared
(logarithmic) wind to the wave exponentially varying with time. It is shown for
hurricane conditions that the Milestype 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, nonresonant excitation of vertically confined
sound waves by stable planar inertiacoriolis 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 stratorotational 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 nonlinear 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 nonlinear scheme. It is established the notion that a nonlinear
scheme is stable if and only if the corresponding scheme in variations is
stable. A new modification of the central LaxFriedrichs (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
operatorsplitting 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 nonlinear
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 nonlinear scheme. The main
methodological innovation is the theorems establishing the notion that a
nonlinear 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
nonsymmetric matrices. A new modification of the central LaxFriedrichs (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
operatorsplitting 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 timeadvance exactly the same
fluid pressure \textit{without solving the Poisson equation}% . In fact, the
incompressible NavierStokes 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
(20042007). 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 socalled NavierStokes (NS) dynamical system. As a
consequence it is found that the fluid pressure obeys a welldefined 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 KantLaplace theory of
formation of the planetary system remains valid also for contemporary
cosmogony. Traditionally, axisymmetric magnetohydrodynamic (MHD), and recently
HallMHD 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 diffusiveshockacceleration on
the stability of the accelerating shock is investigated. A twofluid 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 steadystate shocks that contain a gas subshock. The downstream
conditions of the latter are parametrized by the ratio of the upstream
accleratedparticles 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 entropyvortex waves.