Stochastic modeling of tire–snow interaction using a polynomial chaos approach

Developing accurate models to simulate the interaction between pneumatic tires and unprepared terrain is a demanding task. Such tire–terrain contact models are often used to analyze the mobility of a wheeled vehicle on a given type of soil, or to predict the vehicle performance under specified operational conditions (as related to the vehicle and tires, as well as to the running support). Due to the complex nature of the interaction between a tire and off-road environment, one usually needs to make simplifying assumptions when modeling such an interaction. It is often assumed that the tire–terrain interaction can be captured using a deterministic approach, which means that one assumes fixed values for several vehicle or tire parameters, and expects exact responses from the system. While this is rarely the case in real life, it is nevertheless a necessary step in the modeling process of a deterministic framework. In reality, the external excitations affecting the system, as well as the values of the vehicle and terrain parameters, do not have fixed values, but vary in time or space. Thus, although a deterministic model may capture the response of the system given one set of deterministic values for the system parameters, inputs, etc., this is in fact only one possible realization of the multitude of responses that could occur in reality. The goal of our study is to develop a mathematically sound methodology to improve the prediction of the tire–snow interaction by considering the variability of snow depth and snow density, which will lead to a significantly better understanding and a more realistic representation of tire–snow interaction. We constructed stochastic snow models using a polynomial chaos approach developed at Virginia Tech, to account for the variability of snow depth and of snow density. The stochastic tire–snow models developed are based on the extension of two representative deterministic tire–snow interaction models developed at the University of Alaska, including the pressure–stress deterministic model and the hybrid (on-road extended for off-road) deterministic model. Case studies of a select combination of uncertainties were conducted to quantify the uncertainties of the interfacial forces, sinkage, entry angle, and the friction ellipses as a function of wheel load, longitudinal slip, and slip angle. The simulation results of the stochastic pressure–stress model and the stochastic hybrid model are compared and analyzed to identify the most convenient tire design stage for which they are more suitable. The computational efficiency of the two models is also discussed.     

Parameters of jets out of rectangular nozzles when they are free or near a screen

The results are given of an experimental determination of the parameters of jets out of rectangular nozzles. The distributions of the mean velocity and an impurity concentration were measured. The rearrangement of the jet flow associated with the three-dimensional structure of the jet and interaction of the jet and a screen was investigated. A model that describes the occurrence of a pressure difference and curvature of the jet trajectory when it interacts with a screen is proposed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 39–48, July–August, 1979.     

Interaction of wall turbulence with a compliant surface

The effect of various parameters of the compliant surface on the interaction with wall turbulence is analyzed using the monoharmonic approximation. It is shown that the interaction is resonant in character and that for certain values of the parameters a considerable reduction in turbulent skin friction is possible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 67–72, July–August, 1990.     

Evaluation of link-track performances using DEM

A two-dimensional discrete-element model for the interaction between link-track and soil is presented. The model was developed using commercial PFC2D code. Two different particles, sphere and clump of two spheres, were used to represent the soil. The soil parameters of the model were determined using Hertzian contact theory. Based on the model and soil parameters, simulations of biaxial tests and calculations of the internal angle of friction and cohesion were preformed. The simulation results showed that the internal angle of friction should not exceed the value of 0.65 when using the spherical particles. Based on the clumped particles model, simulations of shear tests with two grouser plates (lengths 100 and 150 mm) were performed under different soil conditions, normal pressures, and cleat heights. A curve fitting of the simulation results was performed using three semi-empirical models from Bekker, Janosi, and Wong for representing the shear stress–displacement relationship. The best fitting was achieved using Wong’s approach. The simulation results of the cleat effects were compared with Bekker’s grouser approach and McKyes’s formulation for soil–blade interaction. In most of the cases, the results of Bekker’s model were the lower bound and McKyes’s model, the upper bound of the DEM simulation results. The properties of the soil model for the DEM were determined using simulation results of shear tests by grouser plate. In the range investigated, the size of the shearing grouser plate is not significant in determining the soil model properties.     

“Lateral” interaction of a supersonic underexpanded ideal-gas jet with surfaces of different shape

A numerical investigation is made of the interaction of an underexpanded jet of an inviscid and nonheat-conducting gas issuing from an axisymmetric conical nozzle with plane, cylindrical, and spherical surfaces. It is assumed that the flow turning angle for flow about a barrier is smaller than the critical angle, and subsonic regions are absent in the flow field studied. The effect of the characteristic parameters (Mach number at the nozzle exit, jet underexpansion) on the flow pattern and jet forces is analyzed. The results of numerical calculations are compared to the results of approximate theories and experimental data. A theoretical solution of the problem of the effect of a supersonic jet on a surface of given shape, even in the approximation of an inviscid, nonheat-conducting gas, is quite difficult. A reason for this is that the flow region contains shock waves interacting with each other, contact discontinuities, and zones of mixed sub-and supersonic flow. As far as is known to the authors, the results obtained for three-dimensional problems for the interaction of supersonic jets with each other or with barriers are primarily experimental (for example, [1–6]). A numerical analysis of the interaction of axisymmetric ideal-gas jets was carried out in [7–10]. In [7] a three-dimensional form of the method of characteristics was used to calculate the initial interaction region for two supersonic cylindrical jets (with Mach number M=10) intersecting at an angle of 60. The interaction of several jets has been considered in [8, 9], where the solution was obtained according to the Lax—Wendroff method without elimination of the discontinuity lines of flow parameters. In [10] the lateral interaction of axisymmetric supersonic jets with each other and with a plate is investigated by means of a straight-through calculationTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–8, November–December, 1974.The authors thank A. N. Kraiko for useful discussions of the results, and A. L. Isakov and É. N. Gasparyan for kindly providing the experimental data.     

Parameters of flow near the axis in the shock layer in the interaction of a supersonic jet with a barrier

The axisymmetric interaction between a supersonic jet with a finite expansion ratio and a barrier is accompanied by the formation of complex sub- and supersonic flow in a shock layer whose thickness depends on the parameters of the jet and the position of the barrier. The main relationships of the interaction process have been established experimentally ([1–3] and others) and individual results of numerical calculations of such flows are known [4]. An analytical investigation of the parameters in the shock layer formed ahead of a plane barrier when an underexpanded jet impinges on it is presented below. The results of [5], where the region near the axis of a shock layer of arbitrary thickness is analyzed within the framework of a model of flow with a constant density, is placed at the basis of the analysis.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 63–70, September–October, 1978.The author thanks Yu. M. Tsirkunov for useful discussions.     

Shock wave formation as a result of interaction with a weak discontinuity on the boundary of a local subsonic zone

A self-similar solution, which explains the formation of a strong-family shock wave (Mach number behind the wave less than unity) on the sonic line, is obtained for the Tricomi equation of plane potential flow in hodograph variables. A characteristic with a discontinuity of the derivatives of the gas dynamic parameters arrives at the formation (interaction) point, while the characteristic of the other family leaving this point does not contain a singularity. The intensity of the shock wave varies along its generator in accordance with a power law with an exponent close to unity. At the interaction point the discontinuity of the derivatives along the streamline is equal to infinity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 152–158, July–August, 1990.The results were presented at the G. G. Chernyi seminar. The author is grateful to the seminar director and the participants for useful discussions.     

Stationary supersonic conducting gas flow in a channel with nonconducting walls under weak magnetohydrodynamic interaction

A solution of the problem of flow in a channel with nonconducting walls for a small magnetohydrodynamic interaction parameter N is obtained by numerical methods. In the 0–10 range of variation of the Hall and magnetic Reynolds number parameters the distributions of the electrical parameters and the average (over the cross section) and local gasdynamic flow parameters are computed for two different geometries of the applied magnetic field. It is shown that an increase in the Hall and magnetic Reynolds number parameters is accompanied by a diminution in the Joule dissipation and the perturbation of the average (over the cross section) gasdynamic flow characteristics. It is disclosed that the distribution of the gasdynamic parameters over the channel cross section is extremely nonmonotonic in the end current zones.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 20–29, July–August, 1970.In conclusion, the author is grateful to A. B. Vatazhin for useful comments and constant attention to the research and to I. U. Tolmach for valuable comments.     

Calculation of flow of highly attenuated gas around a body,taking into account interactions with the surface

The results of Monte Carlo calculations are given for the aerodynamic characteristics of different bodies in free-molecular and near-free-molecular conditions. The calculations are made for different models of the interaction between the attenuated-gas flow and the surface of the body. Analysis and comparison of the results obtained with available experimental and natural data serves as the basis for a discussion of the applicability of individual parameters of these models in practical calculations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 117–123, September–October, 1978.It remains to thank M. N. Kogan, V. A. Perepukhov, A. I. Erofeev, A. I. Bunimovich, L. L. Zvorykin, A. A. Pyarnpuu, and V. F. Kameko for discussions of the results and for valuable comments.     

Formation of secondary drops during impact interaction of a drop with a liquid surface

This paper reports results from experimental studies of the formation of secondary drops during impact interaction of a drop with a liquid surface. The experimental data are compared with analytical estimates of the parameters of the cavern formed and the Rayleigh column.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 55–62, January–February, 2005.     

Nonregular interaction of shock waves in magnetohydrodynamics

The nonregular quasisteady interaction of plane-polarized shock waves in a magnetic field is considered. Within the framework of the magnetohydrodynamic approximation, a Mach interaction pattern is proposed and investigated. This pattern contains the Mach shock wave moving between the initial shocks and connecting their branch points. Unlike the gasdynamic case, in the vicinity of both branch points the flow contains, in addition to the Mach shock, up to five waves whose qualitative structure cannot be predicted in advance. Numerical methods of solving the problem are elaborated. The numerical simulations of the phenomenon under consideration were carried out by means of a specially designed software package with the entire range of all the upstream parameters being covered. It was found that the velocity and inclination of the Mach shock depend significantly on the vector of the magnetic induction in the initial state. The wave pattern of the developing flow is a function of all the key parameters, being particularly sensitive to the inclination of the magnetic field. The values of the parameters for which catastrophic restructurings of the wave pattern occur are determined. Sudden changes in the flow and the magnetic field occur in a zone behind the colliding shock waves. The interaction of the shock waves at different angles and for different magnetic fields are analyzed in detail.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 184–198, May–June, 1993.The authors wish to thank G. G. Chernyi, A. B. Vatazhin, A. G. Kulikovskii, and G. A. Lyubimov for helpful discussions of their results.     

Oblique interaction of alfvén and contact discontinuities in magnetohydrodynamics

A general method of solving problems of the interaction of stationary discontinuities is proposed. The problem of the oblique incidence of an Alfvén plane-polarized discontinuity on a contact discontinuity is examined in the general formulation. A solution is constructed numerically over the entire range of variation of the governing parameters. A number of effects associated with the magnetohydrodynamic nature of the interaction are explored. For example, the formation in space of sectors in which the density falls by several orders (almost to a vacuum) is detected. The solutions obtained are of interest, for example, for investigating the interaction between Alfvén discontinuities in the solar wind and the magnetopause, plasmopause and other inhomogeneities whose boundary can be approximated by a contact discontinuity [13–15].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 131–142, January–February, 1990.     

Unsteady reflection of a shock wave from a body with a cylindrical recess

In a number of cases of supersonic flow past bodies with recesses pulsations in the flow arise [1–3]. Experiments [4, 5] indicate that stabilization of the steady supersonic flow past the body with a recess on which a shock wave is incident takes place after a series of oscillations of the bow wave. Numerical calculation of the interaction of a supersonic jet with a cylindrical cavity [6] reveals that damped pressure pulsations arise inside the cavity if the jet is homogeneous, and undamped pulsations it is inhomogeneous. The authors explain the damping of the pulsations by the influence of artificial viscosity. This paper investigates experimentally and theoretically (by numerical methods) the oscillations of the bow shock wave and the parameters of the flow behind it in the case of unsteady reflection of a shock wave from a body with a cylindrical recess turned towards the flow. The problem is posed as follows. A plane shock wave with constant parameters impinges on a cylinder with a cavity. The unsteady flow originating from this interaction is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 199–202, September–October, 1984.     

Some results of a numerical analysis of transient waves in gas suspensions

Special aspects of the transmission of transient waves through gas mixtures carrying suspended solid particles of chemically inert substances are examined. The influence of the parameters of the gas suspension on the conditions governing the occurrence of transient processes is discussed. The interaction of shock waves with a dust-laden half-space is considered. The results of calculations relating to the decay of an arbitrary discontinuity during the reflection of a shock wave from a wall are presented.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 64–69, September–October, 1976.     

Interference of steady shock waves traveling in the same direction

A classification of the possible types of shock-wave structures formed as a result of the interference between overtaking shocks in a homogeneous flow is developed on the basis of a previous study [1]. A series of analytic and numerical interaction type criteria is obtained, which makes it possible to justify and supplement the analysis, carried out in [2], of the regions of the governing flow parameters in which steady-state solutions for shock-wave structures of different types exist. The calculations are found to be consistent with the known experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 143–152, July–August, 1987.     

Nonlinear steady-state oscillations of an elastic plate floating on the surface of a liquid of finite depth

In addition to obtaining solutions by the perturbation method it is shown that in the case of nonlinear wave interaction given a certain relationship between the parameters of the interacting waves steady-state compound waves may exist.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 146–154, May–June, 1989.     

Interaction of unsteady pressure waves and perforated baffles

The results of experiments to investigate the interaction of plane shock waves with a triangular pressure profile and perforated baffles are presented. The parameters of the reflected waves are determined and the amplitude and impulse attenuation coefficients are measured as functions of the permeability of the barrier. The quenching coefficients are determined for the case of two baffles arranged in series.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 174–176, November–December, 1987.