(Brief Overview)

In [J1], [J4], [J5] [J15], [J16], [J17], [P4] the steady and oscillatory flows and heat transfer in toroidal tubes were investigated numerically. The algorithm developed was applied also to the oscillatory flow in spherical annuli [J20], [J27].

In [J14], [J33], [J44], [P6] we proposed to use scaled coordinates in order to keep the boundaries of the computational domain fixed. Based on this notion, a difference scheme was devised in [J33] for solving the Navier-Stokes equations with a free boundary. As a featuring example the stationary rising of bubbles in ideal and viscous liquid was considered. Results were obtained for wide range of governing parameters: Reynolds and Weber numbers.

The scheme of Helmholtz-Kirchhoff was implemented numerically for ideal flows with separation treating the separation line as an unknown free boundary at which both the kinematic condition and Bernoulli integral hold. Flows around blunt 2D profiles (cylinder, airfoils) were treated by the said technique and in each case two distinct solutions for the shape of free boundary were found: open and cusp-ended zone ([J44], [J52], [P16], [P18], [P25]). A further development of this approach involved quasi-parabolic coordinates [P32].

In [P45] discriminant and cluster analysis were used to develop a procedure for pattern recognition of macrosynoptic processes. The adaptive grids from [P8], [P9], were applied to the 2D limited-area barotropic models for numerical weather prediction in [P44]. In [P43] we developed an algorithm for approximation of 4D geopotential field with multivariate B-splines in least-square sense.

A new coupled system of two strongly nonlinear elliptic equations is derived for 2D orthogonal meshes with a manageable Jacobian [P8], [P9]. It is shown how the boundary conditions are to be posed in order to have coordinate lines that are orthogonal to the domain boundaries. A difference scheme of splitting type is devised in a manner to preserve the coupling of the system at each fractional time-step.

A block-three-diagonal algorithm of Gaussian elimination was developed in [P5]. For multidiagonal systems, the algorithm for Gaussian elimination with pivoting was implemented in FORTRAN in [P33]. Some specialized minimization procedures are developed in [J8], [P2].

The application of longitudinal waves for evaluation of surface dilational viscosity was examined critically in [J22] and the limitations due to the unsteady nature of the flow were shown. A new experimental setup is proposed [J23] employing a steady flow and the respective viscometric formula was derived in which the surface elasticity does not take part.

In [J6] is discussed the application of the principle of minimal dissipation to turbulent flows with nearly plane-parallel mean velocity profile. The critical revision of the experimental data [J24], showed that it governs the lower Reynolds number of transition to turbulence in Poiseuille flows. It was instrumental in investigating flows with random structure [P11], [P15], [P17], [P22], [B1].

In many different situations the random physical fields can be quantitatively very well approximated with random point functions. The technique of Volterra-Wiener functional series is corroborated in [J21], [J28], [J31], [J32], [P10], [P12], [B1].

The most straightforward application of the random point approximation was in particulate continua with random constitution where the random point nature is obvious [J31], [J32], [J35], [J36], [J59], [P10], [P12]. The famous Maxwell formula for overall conductivity and Einstein formula for the effective viscosity were rigorously derived in the linear order of the random point approximation.

On a heuristic basis, the random point approximation was used also as an alternative approach to the closure problem in turbulence. Application to Burgers turbulence was demonstrated in [J7], [J13], [J21] (see also [J40], [J57], [P24], [B1]). The first kernel of the functional expansion (the "linear" term) represents the shape of localized structures (called by the experimentalists "coherent structures") that are randomly dispersed throughout the space or in time . When the system is near the threshold of the instability, the first term accurately represents the random flow [P38], [B1]. It turned out that a first-order truncation (neglecting the pair-, triple-, etc., interactions among the structures) can be used on a heuristic basis even for systems that are far from the threshold of instability. The quantitative agreement of the predicted statistical characteristics with turbulence measurements was encouraging.

The method was also applied to Poiseuille flows [J26], [J29], [P11], [P15], [P22]; Lorenz attractor [J42], [P23], [P24]; Homsy-Kuramoto-Sivashinsky equation [J47], [J48], [J57] and self-similar random solutions of Navier-Stokes equations [J49], [J54], [P13], [P17], [P21], [P38], e.g., the plane mixing layer (see also the monograph [B1]).

In [J53]
the generalization of the Boussinesq approximation for large-scale
atmospheric motions is put on a consistent asymptotic basis. In
[J19]
an integro-differential equation for the deflection of the cupula in
human vestibular apparatus was derived. TE CO_{2} lasers are
modeled in
[J12],
[J18].
The works
[J3],
[J9],
[J12]
were devoted to some problems in plasticity. In
[J10],
[J11],
[J38],
[J39],
inhomogeneous (e.g., reacting) convective flows were investigated
numerically. In
[P1],
[P3]
the stability of elastic bars subjected to non-conservative loading
was studied.

In [J37], [J51] a singular asymptotic solution for the high frequency oscillatory flow in eccentric spherical annuli is developed. In [J30] we elaborated a Fourier-series expansion for the gradient linear temperature field around two non-equal and non-overlapping spheres. In the 2D case (gradient flow around two cylinders) the solution is found numerically [P7]. Two-dimensional model of convective cloud was proposed in [J82].