Pulsating regime corresponding to an obstacle in a stationary inhomogeneous flow

The pulsating regime produced by the presence of a cylindrical cavity in a stationary inhomogeneous supersonic flow is simulated mathematically. The system of equations for an inviscid thermally nonconducting gas is solved by a numerical method based on a two-step difference scheme of second order of approximation. This method makes it possible to calculate in each time step the complete flow field at once, which makes it possible to follow the development of the nonstationary flow, which in the present case is a pulsating flow. The flow pattern in the pulsating regime is studied in detail. The pressure pulsations in the cavity are due to the alternating passage through it of shock waves and rarefaction waves, and the pulsations are nonlinear. The influence of the basic parameters on the characteristics of the pulsating flow is studied and some estimates are made.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 64–71, September–October, 1979.     

Separation of isochromatics and isopachics using a Faraday rotator in dynamic-holographic photoelasticity

A method is presented that allows the simultaneous separation of isochromatic- and isopachic-fringe patterns for transient-plane stress problems. Isopachic fringes are obtained by means of holography with a Faraday cell and a pulsed ruby laser flashing dual pulses. As usual isochromatic whole-order fringes are recorded in a circular-light polariscope. The shock generator (air-gun) and its synchronizing system with the ruby laser is described. The procedure is applied to the recording of the isochromatic- and isopachic-fringe patterns in a disk under radial dynamic loads.     

Experimental investigation of the development of motion of liquid in conduits

Only a few studies have been devoted to the experimental study of the initial stage of the motion of a liquid from the state of rest in a closed delivery conduit [1]. It can be concluded on the basis of the results of these studies that at the beginning of the process the mechanical energy losses are smaller than in quasistationary flow. These studies also suggest that the laminar nature of the flow persists in the nonstationary flow. However, investigations are of an integral nature and therefore in them the flow structure is not determined. In the present study the development of the motion of the liquid in a delivery conduit from the state of rest is investigated. The tangential frictional stresses at the wall of the conduit, measured by the thermal anemometric method, show that the transitional Reynolds number, at which the laminar flow regime changes into turbulent, depends on the acceleration of the flow and far exceeds the critical value for the case of the stationary flow. At maximum acceleration of the flow equal to 11.78 m/sec² the transition of the laminar regime to the turbulent at the wall of the conduit occurs at Re = 234, 500. The loss coefficients of mechanical energy have been computed from experimental results, which show that the use of the corresponding coefficient of quasistationary turbulent flow in the computation leads to appreciable errors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 79–85, November–December, 1977.     

Hydrodynamic instability in fluid layers with uniform volumetric energy sources

Linear and energy theory stability criteria are presented for fluid layers of infinite horizontal extent heated internally by a uniform volumetric energy source. The thermal coupling between the layer and its environment is modeled by a general mixed boundary condition in both the conduction state and the disturbance temperature. Rigid-rigid, free-free, free-rigid, and rigid-free boundaries are considered in the computations. For a fixed ratio of upper surface Biot number to that at the lower surface, decreasing the Biot number is strictly destabilizing for both linear and energy theory criteria. A region of possible subcritical instability is found; its size is strongly dependent on Biot number and becomes small for small values of lower surface Biot number and large Biot number ratio. For two rigid surfaces and an upper and lower surface Biot number of 47.5, mean energy transport measurementswithin the convecting layer indicate a critical Reyleigh number close to that predicted by linear theory. Subcritical instability is not observed when finite amplitude disturbances are introduced at a Rayleigh number between the critical values predicted by the linear theory and the energy theory.     

Comparison of Combustion Models and Assessment of Their Applicability to the Simulation of Premixed Turbulent Combustion in IC-Engines

In order to simulate the turbulent combustion process occurring in spark-ignition (IC) engines, it is necessary to provide suitable and numerically economical models, the latter being particularly important in the application to industrial problems. Moreover, these models must deliver sufficiently accurate results for the unsteady operation of spark combustion engines, concerning variable geometries, temperatures, pressures and charge development in different configurations. In this work different turbulent combustion models for premixed hydrocarbon combustion are compared with respect to their ability to accurately predict the propagation of turbulent perfectly premixed flames. As a first configuration a cylinder of constant volume was studied. Transient calculations were used to simulate the propagation of the turbulent flame and to evaluate the resulting turbulent burning velocity. These calculations were performed for a perfect mixture of air and hydrocarbons at stoichiometric mixture and different initial conditions concerning pressure, temperature and turbulence intensity. As a second configuration a stationary turbulent bunsen-type flame with methane fuel was used to validate the turbulent combustion model of [Lindstedt and Vaos, Combust. Flame 116 (1999) 461] at different pressures. Calculated results were then compared to experimental data of [Kobayashi, Tamura, Maruta and Niioka. In: Proceedings of the 26th Symposium on Combustion, 1996, p. 389] and show excellent agreement for the turbulent burning velocity at several pressure levels using only a single set of model parameters.     

Experimental measurements and computational modeling of the flow field in an idealized human oropharynx

The velocity field in the central sagittal plane of an idealized representation of the human oropharynx (HOP) during steady inspiration, simulating oral inhalation through an inhaler mouthpiece, was measured experimentally using endoscopic particle image velocimetry (PIV). Measurements were made at three flow rates: 15, 30, and 90 L/min, which correspond to a wide range of physiological conditions. Extensive tests were performed to verify the veracity of the PIV data. The flow was also modeled computationally using Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) methods. The PIV data clearly indicate the complex nature of HOP flow, with three-dimensionality and several regions of separation and recirculation evident. Comparison of the experimental and computational results shows that, although the RANS CFD reproduces the basic features of the flow, it does not adequately capture the increased viscous effects at lower Reynolds numbers. The results demonstrate the need for more development and validation of CFD modeling, in particular RANS methods, in these flows.     

Deformation of a Peridynamic Bar

The deformation of an infinite bar subjected to a self-equilibrated load distribution is investigated using the peridynamic formulation of elasticity theory. The peridynamic theory differs from the classical theory and other nonlocal theories in that it does not involve spatial derivatives of the displacement field. The bar problem is formulated as a linear Fredholm integral equation and solved using Fourier transform methods. The solution is shown to exhibit, in general, features that are not found in the classical result. Among these are decaying oscillations in the displacement field and progressively weakening discontinuities that propagate outside of the loading region. These features, when present, are guaranteed to decay provided that the wave speeds are real. This leads to a one-dimensional version of St. Venant's principle for peridynamic materials that ensures the increasing smoothness of the displacement field remotely from the loading region. The peridynamic result converges to the classical result in the limit of short-range forces. An example gives the solution to the concentrated load problem, and hence provides the Green's function for general loading problems.     

Some features of flexural-gravity wave propagation in the presence of a crack in sheet ice

The effect of a crack in an ice sheet on the propagation of surface flexural-gravity waves in a basin of constant depth is analyzed. The ice sheet is simulated by two floating semi-infinite fragments of a thin elastic isotropic plate. As the boundary-contact conditions on the line of contact between the parts of the plate the conditions of continuity of displacements for arbitrary slopes simulating one ice-floe overlying on another and free-edge conditions (crack) are considered. The dependence of the amplitude characteristics of the perturbations on the thickness of the ice, its degree of compression, the incident wave frequency, the depth of the basin, and the form of the boundary-contact conditions is investigated. Problems of wave diffraction on inhomogeneities of an elastic plate were solved in [1, 2], and on a crack in the ice sheet in [3, 4].Sevastopol. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 144–150, March–April, 1996.     

Dynamic compression-shear response of brittle materials with specimen recovery

A new configuration for compression-shear soft-recovery experiments is presented. This technique is used to investigate various failure mechanisms during dynamic multiaxial loading of an Al₂O₃/SiC nanocomposite and TiB₂. Velocity profiles of the target surface are measured with a variable sensitivity displacement interferometer, yielding normal and transverse velocity-time histories. A dynamic shear stress of approximately 280 MPa is obtained, in the Al₂O₃/SiC nanocomposite, for an imposed axial stress of about 3.45 GPa on a 540 m thick sample. This dynamic shear stress is well below the value predicted by elastic wave propagation theory. This could be the result of stress-induced damage and inelasticity in the bulk of the sample or inelasticity on the sample surface due to frictional sliding. To gain further insight into the possible failure mechanisms, an investigation of compression-shear recovery techniques, with simultaneous trapping of longitudinal and lateral release waves, is conducted.