Hydraulic analysis of free-surface flows with solidification

Summary The paper is concerned with a one-dimensional analysis of plane open-channel flow with continuous solidification. The process is of relevance for recent developments in the casting of steel and other metals. The bottom of the channel consists of a rotating casting roll and a horizontal cooling table, where the solidified material is withdrawn with given velocity. The study is restricted to the region downstream of the top of the casting roll. Surface tension is neglected. In the main part of the analysis inviscid fluid flow is considered since the Reynolds number is very large in the applications. It is found that the steady-state solutions are nonunique in a certain parameter range. In addition to a continuous solution, there are two solutions including hydraulic jumps, with one hydraulic jump being located on the casting roll, the other one on the cooling table. Regarding the stability of the non unique solutions, the evolution of disturbances is investigated numerically as an initial-value problem. It is concluded that the hydraulic jump on the cooling table is unstable, while the other discontinuous solution as well as the continuous solution are stable for sufficiently small disturbances. Which stable solution is attained in the steady state, depends on the history of the process. Friction at the liquid/solid interface is taken into account in the last part of the analysis. A constant friction coefficient is assumed. It is found that the history of the process determines the steady-state solution if, and only if, the friction coefficient is sufficiently small. For larger values of the friction coefficient, the steady-state solution is unique and independent of the history of the transient process. Furthermore, for sufficiently large friction coefficients, stable hydraulic jumps are found, in contrast to the inviscid case, also on the cooling table. Received 19 March 1999; accepted for publication 3 May 1999     

Indentation of an Elastic Half-space by a Rigid Flat Punch as a Model Problem for Analysing Contact Problems with Coulomb Friction

One important problem which still remains to be solved today is the uniqueness of the solution of contact problems in linearized elastostatics with small Coulomb friction. This difficult question is addressed here in the case of the indentation of a two-dimensional elastic half-space by a rigid flat punch of finite width, which has been previously studied by Spence in Proc. Camb. Philos. Soc. 73, 249–268 (1973). It is proved that all the solutions have the same simple structure, involving active contact everywhere below the punch and a sticking interval surrounded by two inward slipping intervals. All these solutions show the same local asymptotics for surface traction and displacement at a border between a sticking and a slipping zone. These asymptotics describe (soft) singularities, which are universal (they hold with any geometry) and are explicitly given. It is also proved that in cases where the friction coefficient is small enough, the sticking intervals present in two distinct solutions, if two distinct solutions exist, cannot overlap.     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

Thermo-mechanical contact behavior of a finite graded layer under a sliding punch with heat generation

The problem of a rigid punch contacting with a finite graded layer on a rigid substrate is investigated within the framework of steady-state plane strain thermoelasticity, in which heat generated by contact friction is considered with a constant friction coefficient and inertia effects are neglected. The material properties of the graded layer vary according to an exponential function in the thickness direction. Fourier integral transform method and transform matrix approach are employed to reduce the current thermocontact problem to the second kind of Cauchy-type singular integral equation. Distributions of the contact pressure and the in-plane stress under the prescribed thermoelastic environment with different parameter combinations, including ratio of shear moduli, relative sliding speed, friction coefficient and thermal parameters are obtained and analyzed, as well as the stress singularity and the stress intensity factors near the contact edges. The results should be helpful for the design of surfaces with strong wear resistance and novel graded materials for real applications.     

渤、黄海的底摩擦系数

利用伴随法，根据32个验潮站的水位观测资料来优化渤、黄海的底摩擦系数．所作的一系列数值模拟实验表明；利用伴随方法优化渤、黄海的底摩擦系数能够有效地减少数值模拟结果与观测结果的差异．     

Dam-break flows at a jump in the width of a rectangular channel

The problem of dam-break flow at a jump in the width of a rectangular channel is studied in the first shallow-water approximation. Two cases where the upstream channel width is greater or smaller than the downstream channel width are considered. It is shown that in the first case, the problem is uniquely solvable under the assumption that the total energy of the flow is conserved at the jump in the channel width, and in the second case, a solution of the problem for some initial data exists only provided that the total energy of the flow is lost at this jump.     

An Apparent Interface Location as a Tool to Solve the Porous Interface Flow Problem

A new approach for solving the laminar flow problem above a porous medium is presented here, using an apparent interface for which both superficial velocity and intrinsic shear stress are continuous. The derivation of this approach is based on a detailed investigation of the Ochoa-Tapia and Whitaker (Int. J. Heat Mass Transfer 38:2635–2646, 1995a) jump condition and its sensitivity to the value of β (the jump condition coefficient) and to an error in the interface location. While the value of the jump condition coefficient is highly sensitive to the interface location, the new apparent interface approach does not require an a priori information about the location of the interface. This approach can be easily used knowing only one measurable parameter—the maximum velocity or the flow rate. Validation of the apparent interface approach against measurements from other works shows that it can be successfully used to predict the velocity profile for different geometrical models.     

Subsonic semi-infinite crack with a finite friction zone in a bimaterial

Propagation of a semi-infinite crack along the interface between an elastic half-plane and a rigid half-plane is analyzed. The crack advances at constant subsonic speed. It is assumed that, ahead of the crack, there is a finite segment where the conditions of Coulomb friction law are satisfied. The contact zone of unknown a priori length propagates with the same speed as the crack. The problem reduces to a vector Riemann–Hilbert problem with a piece-wise constant matrix coefficient discontinuous at three points, 0, 1, and ∞. The problem is solved exactly in terms of Kummer's solutions of the associated hypergeometric differential equation. Numerical results are reported for the length of the contact friction zone, the stress singularity factor, the normal displacement u₂, and the dynamic energy release rate G. It is found that in the case of frictionless contact for both the sub-Rayleigh and super-Rayleigh regimes, G is positive and the stress intensity factor K_II does not vanish. In the sub-Rayleigh case, the normal displacement is positive everywhere in the opening zone. In the super-Rayleigh regime, there is a small neighborhood of the ending point of the open zone where the normal displacement is negative.     

Strain rate intensity factors for a plastic mass flow between two conical surfaces

The concept of strain rate intensity factor was introduced in [1], where the asymptotic expansion of the velocity field in a perfectly rigid-plastic material was obtained near the maximum friction surface, which is determined by the condition that the specific friction forces on this surface are equal to the simple shear yield strength. In particular, it was shown in this paper that near the maximum friction surface the equivalent strain rate (the second invariant of the strain rate tensor) tends to infinity inversely proportional to the square root of the distance to this surface. We note that the same result was obtained in the case of plane flow in [2]. The strain rate intensity factor is defined to be the coefficient of the leading singular number in the series expansion of the equivalent strain rate near the maximum friction surface. It was shown in [3] that there is a sufficiently complete formal analogy between the strain rate intensity factor and the stress intensity factor in mechanics of cracks [4]. In [5], it was suggested to use the concept of strain rate intensity factor to estimate the thickness of the layer near the friction surface where one should take into account viscosity effects. (Thus, this is an intensive strain layer formed as a result of a very large equivalent strain rate.) Therefore, the problem of calculating the strain rate intensity factor in specific processes is topical in the development of the general concept based on the use of the strain rate intensity factor and its applications in the theory of metal forming processes. These factors have already been calculated for several processes such as plane upsetting and drawing [3]. In the present paper, we calculate the distribution of the strain rate intensity factor in a plastic mass flow through an infinite converging channel formed by two conical surfaces on which the law of maximum friction acts (Fig. 1). A specific characteristic of this problem is the existence of two maximum friction surfaces and, accordingly, two distributions of the strain rate intensity factor. Since, according to the theory [5], the strain rate intensity factor is related to the thickness of the intensive strain layer near the friction surface, the solution of this problem may serve as a starting point for experimental confirmations of the theory. Note that the intensive strain layer thickness can be determined experimentally without any difficulties [6, 7] and the flow in an infinite channel of the shape under study can successfully model the tube drawing process [8].     

带控制片方柱绕流的非定常数值模拟

对高雷诺数下带控制片的方体(在方柱之前放置一小尺度的柱形薄片)非定常绕流进行了数值模拟。数值模拟方法采用RNG重整化群紊流模型，SIMPLE算法，在非定常计算中引入双重时间步方法，将方体和控制片区域作为求解域的一部分作整体求解。结果表明，控制片与方柱之间距离的不同，使得控制片所形成的涡被抑止或产生卡门涡街，从而对柱体侧面边界层的分离发生影响，产生三种典型的流态，柱体时均阻力系数相对于无控制片的情况迅速减小，当控制片偏离柱体轴心时，柱体升力系数迅速增大，偏心至柱体外壁面附近达到极大值。阻力系数达到最小值后将产生跳跃式的增大，而升力系数在达到最大值后也将产生跳跃式的减小。     

Wave induced sliding at the interface between a layered elastic medium and a half-space

The sliding interface between an unrestrained elastic half-space and a grounded layered half-space excited by an incident harmonic wave is investigated. The problem is formulated considering various possible boundary conditions and boundary inequalities at the sliding interface. The Coulomb friction model without distinction between the static and kinetic coefficients of friction is considered to govern the sliding condition. Three possible bands at the interface, namely slip, stick, and separation, are considered. The interface is assumed to be preloaded under normal and shear stresses. The solution is developed by modifying the problem of welded interface, which then is reduced to a set of algebraic equations. The effects of the incident angle, layer thickness, friction coefficient and externally applied stresses on the drift velocity of the unrestrained half-space are studied numerically for a pair of materials. It is shown that the sliding interface, and hence the drift velocity of unrestrained half-space is noticeably influenced by the layered medium. These results are expected to be useful for the development of a new kind of ultrasonic drive in future.     

Experiments on the flow of a thin liquid film over a horizontal stationary and rotating disk surface

Experiments on characterization of thin liquid films flowing over stationary and rotating disk surfaces are described. The thin liquid film was created by introducing deionized water from a flow collar at the center of an aluminum disk with a known initial film thickness and uniform radial velocity. Radial film thickness distribution was measured using a non-intrusive laser light interface reflection technique that enabled the measurement of the instantaneous film thickness over a finite segment of the disk. Experiments were performed for a range of flow rates between 3.0 lpm and 15.0 lpm, corresponding to Reynolds numbers based on the liquid inlet gap height and velocity between 238 and 1,188. The angular speed of the disk was varied from 0 rpm to 300 rpm. When the disk was stationary, a circular hydraulic jump was present in the liquid film. The liquid-film thickness in the subcritical region (downstream of the hydraulic jump) was an order of magnitude greater than that in the supercritical region (upstream of the hydraulic jump) which was of the order of 0.3 mm. As the Reynolds number increased, the hydraulic jump migrated toward the edge of the disk. In the case of rotation, the liquid-film thickness exhibited a maximum on the disk surface. The liquid-film inertia and friction influenced the inner region where the film thickness progressively increased. The outer region where the film thickness decreased was primarily affected by the centrifugal forces. A flow visualization study of the thin film was also performed to determine the characteristics of the waves on the free surface. At high rotational speeds, spiral waves were observed on the liquid film. It was also determined that the angle of the waves which form on the liquid surface was a function of the ratio of local radial to tangential velocity.     

Friction measurement at low sliding speed using a servohydraulic tension-torsion machine

A method is presented for the determination of the coefficient of friction at low sliding speeds down to 3 nm/s. It is also possible to vary the contact pressure. The maximum average contact pressure is 40 MPa in the present test configuration, where a servohydraulic tension-torsion machine is used. The method, involving test equipment, measurement and the principles of evaluation, is exemplified by a study where the variation of the coefficient of friction with sliding velocity contact pressure and surface roughness is investigated for the material combination steel and concrete. The measurements have been performed as a 2³ factorial design. An error analysis was performed and the relative measurement uncertainty when determining the coefficient of friction was estinated to be 1.2 percent.     

Turbulent thermal boundary layer on a plate. Reynolds analogy and heat transfer law over the entire range of prandtl numbers

Arational asymptotic theory is proposed,which describes the turbulent dynamic and thermal boundary layer on a flat plate under zero pressure gradient. The fact that the flow depends on a finite number of governing parameters makes it possible to formulate algebraic closure conditions relating the turbulent shear stress and heat flux with the gradients of the averaged velocity and temperature. As a result of constructing an exact asymptotic solution of the boundary layer equations, the known laws of the wall for velocity and temperature, the velocity and temperature defect laws, and the expressions for the skin friction coefficient, Stanton number, and Reynolds analogy factor are obtained. The latter makes it possible to give two new formulations of the temperature defect law, one of which is identical to the velocity defect law and contains neither the Stanton number nor the turbulent Prandtl number, and the second formulation does not contain the skin friction coefficient. The heat transfer law is first obtained in the form of a universal functional relationship between three parameters: the Stanton number, the Reynolds number, and the molecular Prandtl number. The conclusions of the theory agree well with the known experimental data.     

Steady Sliding Frictional Contact Problems in Linear Elasticity

The existence and uniqueness of an equilibrium solution to frictional contact problems involving a class of moving rigid obstacles is studied. At low friction coefficient values, the steady sliding frictional contact problem is uniquely solvable, thanks to the Lions-Stampacchia theorem on variational inequalities associated with a nonsymmetric coercive bilinear form. It is proved that the coerciveness of the bilinear form can be lost at some positive critical value of the friction coefficient, depending only on the geometry and the elastic properties of the body. An example presented here, shows that infinitely many solutions can be obtained when the friction coefficient is larger than the critical value. This result is paving the road towards a theory of jamming in terms of bifurcation in variational inequality. The particular case where the elastic body is an isotropic half-space is studied. The corresponding value of the critical friction coefficient is proved to be infinite in this case. In the particular case of the frictionless situation, our analysis incidentally unifies the approaches developed by Lions-Stampacchia (variational inequalities) and Hertz (harmonic analysis on the half-space) to contact problems in linear elasticity.     

The bouncing motion appearing in a robotic system with unilateral constraint

The hopping or bouncing motion can be observed when robotic manipulators are sliding on a rough surface. Making clear the reason of generating such phenomenon is important for the control and dynamical analysis for mechanical systems. In particular, such phenomenon may be related to the problem of Painlevé paradox. By using LCP theory, a general criterion for identifying the bouncing motion appearing in a planar multibody system subject to single unilateral constraint is established, and found its application to a two-link robotic manipulator that comes in contact with a rough constantly moving belt. The admissible set in state space that can assure the manipulator keeping contact with the rough surface is investigated, and found which is influenced by the value of the friction coefficient and the configuration of the system. Painlevé paradox can cause either multiple solutions or non-existence of solutions in calculating contact force. Developing some methods to fill in the flaw is also important for perfecting the theory of rigid-body dynamics. The properties of the tangential impact relating to the inconsistent case of Painlevé paradox have been discovered in this paper, and a jump rule for determining the post-states after the tangential impact finishes is developed. Finally, the comprehensively numerical simulation for the two-link robotic manipulator is carried out, and its dynamical behaviors such as stick-slip, the bouncing motion due to the tangential impact at contact point or the external forces, are exhibited.     

Identification of In Situ Frictional Properties of Bolted Assemblies with Digital Image Correlation

The present paper is devoted to the identification of frictional properties in bolted assemblies. It is shown that kinematic data provided by digital image correlation can be used to analyze the change of the friction coefficient with the number of cycles. Two approaches are followed. The first one is based on the displacement jump between the assembled plates, and the second one relies on displacement fields measured on the same surface.     

Discontinuity of quasi-static solution in the two-node Coulomb frictional system

A single-node system introduced by Klarbring has provided insight into the non-uniqueness of solution in the quasi-static contact problem at the high coefficient of Coulomb friction. Here, we explore this issue for the two-node system under the slip displacement space in which the instantaneous condition is efficiently represented. In the paper, we identify a qualitatively different failure of the quasi-static evolution algorithm in which a more complex dynamic transition may occur. When the system evolves from the point where both-node discontinuity occurs, the transient evolution behavior involving a damping matrix is explored in order to investigate a final state of the two-node system. It is demonstrated that the final state is uniquely determined which is independent of the damping matrix.     

Estimating the dynamic coefficient of friction from experimental data: the motion of a rocket inside a smoothbore launcher

The motion of a spinning rocket inside a smoothbore launcher has drawn some attention in recent years. The inaccuracy of spin-stabilized rockets can be attributed, in part, to the initial motion of the rocket inside the smoothbore launcher. Recently, tests have indicated that a closer examination of the problem is required. The tests were performed to determine the time history of the spin and forward motion of the rocket by Fastax cameras. The evaluation of the data with respect to engineering-design parameters is difficult since the experimental results must be compared with an analytical model of the system. In this paper, the results of one test are reported and a nonlinear analytical model is used which includes both the rotary and forward motion of the rocket. The two motions are coupled through a dynamic friction coefficient. The intent of this paper is to show how an analytical model can be made to fit the experimental data, that is, the initial conditions and the coefficient of friction are found which define the solution giving the best least-square fit to the data.