Quasistatic propagation of steps along a phase boundary

We study quasistatic propagation of steps along a phase boundary in a two-dimensional lattice model of martensitic phase transitions. For analytical simplicity, the formulation is restricted to antiplane shear deformation of a cubic lattice with bi-stable interactions along one component of shear strain and harmonic interactions along the other. Energy landscapes connecting equilibrium configurations with periodic and non-periodic arrangements of steps are constructed, and the energy barriers separating metastable states are calculated. We show that a sequential one-by-one step propagation along a phase boundary requires smaller energy barriers than simultaneous motion of several steps.      

Effect of a latitudinal temperature gradient on convective motions arising in a rotating spherical layer

The hydrodynamics of planetary atmospheres and Interiors are frequently directly or indirectly connected with convective motions taking place in rotating liquid spherical layers in the field of a central force. Convective stability in a spherical layer at rest, in a central gravity field, was first discussed in [1, 2]. It was shown that the critical Rayleigh number Rao at which convective instability sets in and the wave number of the critical perturbations depend essentially on the thickness of the layer. As in the plane case, the problem of the convective stability of a spherical layer is found to be degenerate, and the form of the critical perturbations cannot be determined from the linear problem. In actuality, minimization of the Rayleigh number permits establishing only the wave numberl for the spherical harmonic Y _l ^m (θ, ?), realized at the limit of stability; the parameter m remains indeterminate and thus 2l+1 independent convective modes correspond to Rao. In [3] a study was made of the convective stability of a liquid in a slowly rotating thin spherical layer. It was shown that the presence of rotation eliminates the degeneracy; at the limit of stability there arise motions corresponding to the Y _l l (θ, ?) -harmonic with a degenerate maximum at the equator, and propagating in a wave manner toward the side opposite to the rotation. In the present work a study is made of the convective stability of a flow of liquid, arising in a rotating spherical layer due to a nonuniform distribution of the temperatures at one of the boundaries of the layer. In such a statement of the problem it is possible to model large-scale motions in the atmospheres of large planets having internal sources of heat and absorbing solar radiation near the cloud cover of the atmosphere. It is established that, depending on the relationships between the parameters imparting the rotation and the inhomogeneous distribution of the temperature, there is either stabilization or destabilization of the layer in comparison with a fixed layer of the same thickness and with the same, but uniformly distributed heat flux supplied to the layer. A study is made of the form of the corresponding critical perturbations.     

Pressure gradients due to gas expansion in the boundary layer combustion of a condensed fuel

Ignition and combustion of a condensed fuel in a gaseous high temperature oxidizing boundary layer flow is analyzed on the basis of higher order boundary layer theory. First order effects due to displacement thickness are taken into account and the pressure gradients generated in the outer potential flow are included in the numerical solution of the governing equations. The strong positive pressure gradients which are induced by expansion in front of the flame generate a low velocity region which facilitates a longitudinal diffusion of heat and mass.Es wird die Zündung und Verbrennung eines festen Brennstoffs in einer heißen gasförmigen oxidierenden Grenzschichtströmung auf der Grundlage der Grenzschichttheorie höherer Ordnung untersucht. Effekte erster Ordnung aufgrund der Verdrängungsdicke werden berücksichtigt und Druckgradienten, die in der äußeren Potentialströmung erzeugt werden, sind in der numerischen Lösung der Bestimmungsgleichungen eingeschlossen. Die starken positiven Druckgradienten, die durch die Expansion vor der Flamme induziert werden, erzeugen ein Gebiet niedriger Geschwindigkeit, was einen Wärme- und Stofftransport entgegen der Hauptströmungsrichtung ermöglicht.     

An analytical prediction of sliding motions along discontinuous boundary in non-smooth dynamical systems

This paper presents a method for the analytical prediction of sliding motions along discontinuous boundaries in non-smooth dynamical systems. The methodology is demonstrated through investigation of a periodically forced linear oscillator with dry friction. The switching conditions for sliding motions in non-smooth dynamical systems are given. The generic mappings for the friction-induced oscillator are introduced. From the generic mappings, the corresponding criteria for the sliding motions are presented through the force product conditions. The analytical prediction of the onset and vanishing of the sliding motions is illustrated. Finally, numerical simulations of sliding motions are carried out to verify the analytical prediction. This analytical prediction provides an accurate prediction of sliding motions in non-smooth dynamical systems. The switching conditions developed in this paper are expressed by the total force of the oscillator, and the nonlinearity and linearity of the spring and viscous damping forces in the oscillator cannot change such switching conditions. Therefore, the achieved force criteria can be applied to the other dynamical systems with nonlinear friction forces processing a C ⁰-discontinuity.     

On the convective instability of a thermal boundary layer

When the surface temperature of a liquid is a harmonic function of time with a frequency, a temperature wave propagates into the liquid. The amplitude of this wave decreases exponentially with distance from the surface. The temperature oscillation is essentially concentrated in a layer of the order of (2/)^1/2, where x is the thermal conductivity of the liquid (thermal boundary layer). Depending on the phase, at certain positions below the surface the temperature gradient is directed downwards and if its magnitude is sufficiently large (the magnitude is a function of the amplitude and frequency of the surface oscillations) the liquid can become unstable with respect to the onset of convection. In that case the convective motion may spread beyond the initial unstable layer. For low frequencies the stability condition can be derived from the usual static Rayleigh criterion, on the basis of the Rayleigh number and the average temperature gradient of the unstable layer. This quasi-static approach, used by Sal'nikov [1], is appropriate to those cases in which the period of the temperature oscillations is much larger than the characteristic time of the perturbations. But when these times are of the same order, the problem must be analyzed in dynamic terms. The stability problem must then be formulated as a problem of parametricresonance excitation of velocity oscillations due to the action of a variable parameter-the temperature gradient.In an earlier work [2] we considered the problem of the stability of a horizontal layer of liquid with a periodically varying temperature gradient. It was assumed that the thickness of the layer was much smaller than the penetration depth of the thermal wave, so that the temperature gradient could be assumed to be independent of position. In the present work we consider the opposite case, in which the liquid layer is assumed to be much larger than the penetration depth, i. e., a thermal boundary layer can be defined. The temperature gradient at equilibrium, which is a parameter in the equations determining the onset of perturbations, is here a periodic function of time and a relatively complicated function of the depth coordinate z. The periodic oscillations are solved by the Fourier method; the equations for the amplitudes are solved by the approximate method of KarmanPohlhausen.The authors are grateful to L. G. Loitsyanskii for helpful criticism.     

Velocities and structure of convective motions in a rotating fluid

The velocities in various rotating fluid regimes are investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 165–168, July–August, 1989.The author wishes to thank G. S. Golitsyn for his constant interest in the work.     

Numerical investigation of natural convective motions in a rotating cube

On the basis of a numerical integration of the Navier-Stokes equations in the Boussinesq approximation, the natural thermal convection in a rotating cube heated from below is investigated for three values of the Rayleigh number. The effect of the rotation velocity on the establishment of various forms of convection flows is studied. The existence of several types of convection flows is noted. Each of the motions is realized on a specific range of the rotation velocity of the cube. The domains of existence of the different types are determined. The results obtained are compared with experimental data.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 53–60, July–August, 1995.     

Double-diffusive convective motions for a saturated porous layer subject to modulated surface heating

Received December 22, 1999     

Effect of viscosity on the propagation of nonlinear Alfvén waves along a tangential magnetohydrodynamic discontinuity in an incompressible fluid

It is proposed to consider the propagation of surface waves along a tangential magnetohydrodynamic discontinuity in the particular case where the fluid velocities on both sides of the interface are equal to zero. In [1] it was shown that waves called surface Alfvén waves may be propagated along the surface separating a semi-infinite region without a field from a region with a uniform magnetic field. The linear theory of surface Alfvén waves in a compressible medium was considered in [2]. In [3] the damping of surface Alfvén waves as a result of viscosity and heat conduction was investigated. The propagation of low-amplitude nonlinear surface Alfvén waves in an incompressible fluid in the absence of dissipative processes is described by the integrodifferential equation obtained in [4]. By means of a numerical solution of this equation it was shown that a perturbation initially in the form of a sinusoidal wave will break. The breaking time was determined. In this paper the equation derived in [4] is extended to the case of a viscous fluid. It is shown that the equation obtained does not have steady-state solutions. The propagation of periodic disturbances is investigated numerically. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 94–104, November–December, 1986. The author wishes to thank L. S. Fedorov for assisting with the calculations.     

Numerical investigation of convective motion in a closed cavity

The investigation of thermal convection in a closed cavity is of considerable interest in connection with the problem of heat transfer. The problem may be solved comparatively simply in the case of small characteristic temperature difference with heating from the side, when equilibrium is not possible and when slow movement is initiated for an arbitrarily small horizontal temperature gradient. In this case the motion may be studied using the small parameter method, based on expanding the velocity, temperature, and pressure in series in powers of the Grashof number—the dimensionless parameter which characterizes the intensity of the convection [1–4]. In the problems considered it has been possible to find only two or three terms of these series. The solutions obtained in this approximation describe only weak nonlinear effects and the region of their applicability is limited, naturally, to small values of the Grashof number (no larger than 10³).With increase of the temperature difference the nature of the motion gradually changes—at the boundaries of the cavity a convective boundary layer is formed, in which the primary temperature and velocity gradients are concentrated; the remaining portion of the liquid forms the flow core. On the basis of an analysis of the equations of motion for the plane case, Batchelor [4] suggested that the core is isothermal and rotates with constant and uniform vorticity. The value of the vorticity in the core must be determined as the eigenvalue of the problem of a closed boundary layer. A closed convective boundary layer in a horizontal cylinder and in a plane vertical stratum was considered in [5, 6] using the Batchelor scheme. The boundary layer parameters and the vorticity in the core were determined with the aid of an integral method. An attempt to solve the boundary layer equations analytically for a horizontal cylinder using the Oseen linearization method was made in [7].However, the results of experiments in which a study was made of the structure of the convective motion of various liquids and gases in closed cavities of different shapes [8–13] definitely contradict the Batchelor hypothesis. The measurements show that the core is not isothermal; on the contrary, there is a constant vertical temperature gradient directed upward in the core. Further, the core is practically motionless. In the core there are found retrograde motions with velocities much smaller than the velocities in the boundary layer.The use of numerical methods may be of assistance in clarifying the laws governing the convective motion in a closed cavity with large temperature differences. In [14] the two-dimensional problem of steady air convection in a square cavity was solved by expansion in orthogonal polynomials. The author was able to progress in the calculation only to a value of the Grashof numberG=10⁴. At these values of the Grashof numberG the formation of the boundary layer and the core has really only started, therefore the author's conclusion on the agreement of the numerical results with the Batchelor hypothesis is not justified. In addition, the bifurcation of the central isotherm (Fig. 3 of [14]), on the basis of which the conclusion was drawn concerning the formation of the isothermal core, is apparently the result of a misunderstanding, since an isotherm of this form obviously contradicts the symmetry of the solution.In [5] the method of finite differences is used to obtain the solution of the problem of strong convection of a gas in a horizontal cylinder whose lateral sides have different temperatures. According to the results of the calculation and in accordance with the experimental data [9], in the cavity there is a practically stationary core. However, since the authors started from the convection equations in the boundary layer approximation they did not obtain any detailed information on the core structure, in particular on the distribution of the temperature in the core.In the following we present the results of a finite difference solution of the complete nonlinear problem of plane convective motion in a square cavity. The vertical boundaries of the cavity are held at constant temperatures; the temperature varies linearly on the horizontal boundaries. The velocity and temperature distributions are obtained for values of the Grashof number in the range 0<G4·10⁵ and for a value of the Prandtl number P=1. The results of the calculation permit following the formation of the closed boundary layer and the very slowly moving core with a constant vertical temperature gradient. The heat flux through the cavity is found as a function of the Grashof number.     

A laser interferometric measurement of the temperature distribution in convective thermal transfer boundary layers

By means of a double mirror interferometry a two-dimensional temperature distribution measurement in convective thermal boundary layers is presented. When the cold air flows along a hot plate model, the interferometric fringe inside the boundary layer will bend. According to the displacement of the fringe and the relation between temperature and index of refraction, a two-dimensional temperature profile is obtained. All is accomplished by optical device with the help of micro-computer without any contact with the flow field. The project supported by the National Natural Science Foundation of China     

Influence of a strong discontinuity on the propagation of second order discontinuities

We study the propagation of second order weak discontinuities in quasi-linear hyperbolic systems of equations with discontinuous coefficients. The general theory is applied to shallow water waves.     

Secondary vibrational convective motions in a flat vertical layer of liquid

Investigations of the stability of steady-state plane-parallel convective motion between vertical planes heated to different temperatures [1–5] have shown that this motion, depending on the value of the Prandtl number P, exhibits instability of two types. With small and moderate Prandtl numbers, the instability is of a hydrodynamic nature. It is brought about by monotonic perturbations which, in the supercritical region, develop into a periodic, with respect to the vertical, system of steady-state vortices at the interface between the opposing convective flows. Articles [6, 7] are devoted to the numerical investigation of nonlinear secondary steady-state flows. If P>11.4, there appears a new mode of instability, i.e., running thermal waves increasing in the flow; with P>12, this mode becomes more dangerous [4]. This instability is connected with the development of vibrational perturbations, and it can be considered that in the supercritical region the perturbations lead to the establishment of steady-state vibrations. Linear theory has made it possible to determine the boundaries of the regions of stability. In the present article a numerical investigation is made of nonlinear supercritical conditions developing as a result of a loss of stability of the steady-state flow with respect to vibrational perturbations.     

Effects of temperature dependent fluid properties and variable Prandtl number on the transient convective flow due to a porous rotating disk

In this paper we have studied the effects of temperature dependent fluid properties such as density, viscosity and thermal conductivity and variable Prandtl number on unsteady convective heat transfer flow over a porous rotating disk. Using similarity transformations we reduce the governing nonlinear partial differential equations for flow and heat transfer into a system of ordinary differential equations which are then solved numerically by applying Nachtsheim–Swigert shooting iteration technique along with sixth-order Runge–Kutta integration scheme. Comparison with previously published work for steady case of the problem were performed and found to be in very good agreement. The obtained numerical results show that the rate of heat transfer in a fluid of constant properties is higher than in a fluid of variable properties. The results further show that consideration of Prandtl number as constant within the boundary layer for variable fluid properties lead unrealistic results. Therefore, modeling thermal boundary layers with temperature dependent fluid properties Prandtl number must treated as variable inside the boundary layer.     

A finite element solution to turbulent diffusion in a convective boundary layer

A second-order closure turbulence model is used to simulate the plume behaviour of a passive contaminant dispersed in a convective boundary layer. A time-splitting finite element scheme is used to solve the set of partial differential equations. It is shown that the second-order closure model compares favourably with recent findings from laboratories, wind-tunnel experiments and large-eddy simulations. We also compare the second-order closure model with the commonly used K-diffusion model for the same meteorological conditions. Case studies also show the effects of model parameters and turbulence variables on the plume behaviour.     

On a magneto-elastic system with discontinuous coefficients and the propagation of a weak discontinuity

Summary In this paper we deduce the magneto-elastic system of equations in the three-dimensional case. As an application we study the propagation of a weak discontinuity when a strong discontinuity also occurs.

---

Sommario In questa Nota si deduce un sistema di equazioni per la Magneto-elasticità con deformazioni finite nel caso tridimensionale. Come Applicazione si studia la propagazione di una discontinuità debole nella ipotesi che i coefficienti delle equazioni risultano essi stessi discontinui. Tale ipotesi risulta verificata, ad es., quando si considerano due differenti mezzi magneto-elastici uno in contatto con l'altro.

---

---

Work supported by the CNR through the Gruppo Nazionale per la Fisica-Matematica.