Some properties of multi-fluid stokes flows

The solution of Stokes' equations for a rotating axisymmetric body which possesses reflection symmetry about a planar interface between two infinite immiscible quiescent viscous fluids is shown to be independent of the viscosities of the fluids and identical with the solution when the fluids have the same viscosity. The result is generalized to a rotating axisymmetric system of bodies which possesses reflection symmetry about each interface of a plane stratified system of fluids. An analogous result for two-fluid systems with a nonplanar static interface is also derived. The effect on torque reduction produced by the presence of a second fluid layer adjacent to a rotating axisymmetric body is considered and explicit calculations are given for the case of a sphere. A proof of uniqueness for unbounded multi-fluid Stokes' flow is given and the asymptotic far field structure of the velocity field is determined for axisymmetric flow caused by the rotation of axisymmetric bodies.     

Numerical modeling of slug flow initiation in a horizontal channels using a two-fluid model

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows that there are three zones associated with slug formation. In addition, long wavelength slugs are shown to initiate from short wavelength waves. These short waves are generated at the interface of the two phases by the Kelvin-Helmholtz hydrodynamic instability. The results obtained through numerical modeling show good agreement with experimental results.     

One-dimensional two-fluid equations for horizontal stratified two-phase flow

Advanced computer codes for water reactor loss-of-coolant analysis are based on the use of the two-fluid model of two-phase flow, in which conservation equations are solved for the gas and liquid phases separately. The standard two-fluid equations, however, sometimes predict the growth of instabilities in the flow, and occasionally become improperly posed. These difficulties have in the past led to the proposal of several different forms for the conservations equations.To help resolve these uncertainties a widely accepted form of the one-dimensional two-fluid equations is used to calculate wave propagation speeds, and stability limits, for the illustrative case of a frictionless horizontal stratified gas-liquid flow. Calculated propagation velocities are shown to agree with the appropriate limit of an exact solution, and the predicted stability limits are found consistent with available observations on the stability of the stratified flow regime.These comparisons help improve confidence in the ability of the two-fluid equations to analyse more complex problems in transient two-phase flow.     

Comparison of waveguide properties of curved versus straight planar elastic layers

Dispersion equations are solved for the in-plane and anti-plane wave propagation in planar elastic layer with constant curvature. The classical Lamé formulation of displacements via elastic potentials is applied and appropriate simplifications are employed. The dispersion diagrams in each case are compared with their counterparts for a straight layer, e.g., the classical Rayleigh–Lamb solution. The curvature-induced symmetry-breaking effects are investigated for layers with symmetric boundary conditions. The role of curvature is also investigated in the cases, when the boundary conditions are not symmetrical. The elementary Bernoulli–Euler theory is employed to analyze the wave guide properties of a curved planar elastic beam in its in-plane deformation. The validity range of the Bernoulli–Euler theory is assessed via comparison of dispersion diagrams.     

Elastica of a variable-arc-length circular curved beam subjected to an end follower force

This paper presents postbuckling behaviors of a variable-arc-length (VAL) circular curved beam subjected to an end follower force. One end of the VAL circular curved beam is hinged while the other end is supported by a frictionless slot, which is fixed horizontally and vertically but is allowed to rotate corresponding to loading direction. When the VAL circular curved beam is deformed, the total arc-length of the circular curved beam varies. Two approaches have been applied for the solution of this problem. The first approach is an elliptic integrals method based on elastica theory, which yields the exact closed-form solution in terms of the first and second kinds of elliptic integrals. For validation of the results, the shooting method is employed for a numerical solution by developing the set of nonlinear governing differential equations together with boundary conditions, and then integrating them by using the fourth-order Runge–Kutta algorithm. The results from both approaches are in very good agreement. From the results, it is found that the VAL circular curved beam subjected to an end follower force can be deformed in many mode shapes. For the first and third modes, the beam exhibits both stable and unstable configurations, whereas for the second mode only an unstable configuration exists. The influences of initial curvature on the critical load and the deformed configurations are highlighted.     

Prediction of liquid level distribution in horizontal gas-liquid stratified flows with interfacial level gradient

A one-dimensional momentum equation has been derived based on a two-fluid model and used to predict the axial distribution of liquid level or void fraction in steady, cocurrent, gas-liquid stratified flows in horizontal circular pipes and rectangular channels. The equation is carefully formulated to incorporate the effect of interfacial level gradient. Two different critical liquid levels are found from the momentum equation and are adopted as a boundary condition to calculate the liquid level or void fraction distribution upstream of the channel exit. The predicted void fraction distributions are compared with the experimental data obtained in a rectangular channel in this work and other data reported for large-diameter pipes. Good agreement is shown for air-kerosene, air-water and stream-water stratified flows with a smooth gas-liquid interface.     

Computation of normal modes of three-dimensional resonators by semiclassical and variational methods: seismological implications

The problem investigated in this paper is that of calculating the eigenfrequencies of the bouncing-ball modes of oscillation of scalar waves in a three-dimensiol, arbitrarily shaped enclosure or resonator. With a view to future applications in seismology, we consider wave resonators limited above by a flat reflector—the earth's surface—and below by a curved one—the sediment/rock interface. The surface of the curved reflector is smooth but otherwise of arbitrary shape. The EBK (Einstein-Brillouin-Keller) method of semiclassical quantization is used to compute the eigenfrequencies. Conventional ray tracing is used to determine the trajectories of the rays in configuration space as well as in the Poincaré surfaces of section needed for the application of the method. For a prescribed set of mode numbers the EBK conditions produce a set of algebraic equations whose solution gives the required eigenfrequency and the corresponding eigentrajectory. Additionally, a system of four nonlinearly coupled ODE's is derived that simulates accurately the trajectories of the bouncing points of rays and can be used as a fast alternative to ray tracing. A variational method based upon the Rayleigh-Ritz formulation is used to obtain bounds on the eigenfrequencies and to determine the resonator's eigenfunctions. For separable and smooth nonseparable domains both approaches give consistent results. If the curvature of the reflector is convex towards the approaching ray, repeated reflections will result in highly unstable ray trajectories which ultimately becomne chaotic. In this case the EBK method fails. A brief discussion of the inherent tendency to chaotic behavior of the multiply reflecting ray system inside arbitrarily shaped enclosures is presented.     

Lateral squeezing of horizontally stratified viscous layers

A fluid mechanics model for the slow deformations of layered geophysical materials is adopted. The problem consists of the slow viscous flow of two horizantally stratified fluid layers of different viscosities as being laterally squeezed between two frictionless planes. For solution a perturbation scheme is applied and the resulting approximation equations are solved numerically. The interface and free surface shapes, and the streamline patterns within the fluid layers are given.     

带集中质量的旋转柔性曲梁动力学特性分析

本文对带集中质量的平面内旋转柔性曲梁动力学特性进行了研究.基于绝对节点坐标法推导出曲梁单元,其中该曲梁单元采用Green-Lagrangian应变,并根据曲梁变形前后的曲率变化和曲率的精确表达式计算了曲梁单元弹性力所作的虚功.通过虚功原理,利用$\delta$函数和中心刚体与悬臂曲梁之间的固支边界条件,建立了带集中质量的旋转柔性曲梁非线性动力学模型.基于该模型,本文仿真计算了悬臂曲梁的纯弯曲问题和带有刚柔耦合效应的旋转柔性曲梁动力学响应问题,以此分别讨论了所提出曲梁单元的收敛性和动力学模型的正确性.进一步应用D'Alembert原理,将旋转曲梁等效为带离心力的无旋转曲梁,通过线性摄动处理得到系统的特征方程,以此分别研究了旋转角速度、初始曲率和集中质量对曲梁动力学特性的影响.最后重点分析了旋转曲梁的频率转向和振型切换问题,并阐述了两者之间的相互关系.研究结果表明:随着旋转角速度的增大,曲梁的频率特性与直梁的频率特性相近,以及曲梁拉伸变形占主导的模态振型会提前.      

Stability of a Strip Made of a Multilayer Composite with a Curved Structure

The stability of a long plate (strip) made of a layered composite with a curved structure is studied within the framework of the three-dimensional linearized theory of stability. The cases of simply supported and rigidly fixed plates compressed along the layers are considered. The curvature of the layers is allowed for within the framework of the continuum approach. Numerical results are obtained by the finite-element method for a strip consisting of alternating layers of two homogeneous isotropic materials     

A laboratory technique for generating a rotating, two-layer stratified channel flow

 Results are presented from a series of experiments in which the performance of a new type of rotating, stratified (two-layer) channel flow facility has been tested. The flow within the channel is driven by a source-sink system and is designed specifically to provide uniform, rectilinear, horizontal motion and relatively-quiescent conditions in the upper and lower layers respectively of the two-layer configuration. The channel is shown to perform well in this regard over a wide range of external forcing conditions, with negligible erosion of the diffusive interface separating the two constituent (miscible) homogeneous layers over periods that are long compared with the time scale of a typical laboratory experiment. Received: 22 January 1996/Accepted: 23 November 1998     

Thermomechanical peeling in multilayer beams and plates—a solution from first principles

The purpose of this paper is to present a solution for the peeling moment arising from the peeling stress distribution in any interface in a multilayer beam or plate. The solution is implemented for three- and four- layer beams; it is shown that it can readily be implemented for any desired number of layers. The solution is derived from first principles, and is evolved from the well known [Timoshenko, A. 1925. Analysis of bi-metal thermostats. Journal of the Optical Society of America, 11, 233–255] bimetal thermostat analysis. A physical interpretation of the factors that make up the peeling moment is given, enabling quick identification of how any layer property may be changed in order to resist delamination at any interface of interest. The concept of moments being transferred across the interface to cause equal radii of curvature is helpful in understanding the factors that influence the magnitude and direction of the peeling moment at any interface. This analytical method applies equally to multilayer stack-ups cured simultaneously at a common temperature, and to products such as integrated circuit wafers where the layers are formed sequentially at different temperatures.     

Three-dimensional studies on bicomponent extrusion

The present work is concerned with the mathematical modelling and numerical simulation of three-dimensional (3-D) bicomponent extrusion. The objective is to provide an understanding of the flow phenomena involved and to investigate their impact on the free surface shape and interface configuration of the extruded article. A finite element algorithm for the 3-D numerical simulation of bicomponent stratified free surface flows is described. The presence of multiple free surfaces (layer interface and external free surfaces) requires special free surface update schemes. The pressure and viscous stress discontinuity due to viscosity mismatch at the interface between the two stratified components is handled with both a double node (u–v–w–P ₁ –P ₂ –h ₁ –h ₂) formulation and a penalty function (u–v–w–P–h ₁ –h ₂) formulation.The experimentally observed tendency of the less viscous layer to encapsulate the more viscous layer in stratified bicomponent flows of side-by-side configuration is established with the aid of a fully 3-D analysis in agreement with experimental evidence. The direction and degree of encapsulation depend directly on the viscosity ratio of the two melts. For shear thinning melts exhibiting a viscosity crossover point, it is demonstrated that interface curvature reversal may occur if the shearing level is such that the crossover point is exceeded. Extrudate bending and distortion of the bicomponent system because of the viscosity mismatch is shown. For flows in a sheath-core configuration it is shown that the viscosity ratio may have a severe effect on the swelling ratio of the bicomponent system.Modelling of the die section showed that the boundary condition imposed at the fluid/fluid/wall contact point is critical to the accuracy of the overall solution.     

A consideration of curved interfaces in stress-assisted martensite formation

A purely mechanical, sharp interface model is developed to consider curved interfaces that have been observed between martensite phase variants. The approach is based on a theory of small strains as distinct from small displacement gradients. It admits a realistic characterization of each phase with standard elasticity tensors and allows for inhomogeneous states of strain within each phase including inhomogeneous, finite rotations. The model indicates that any signficant interface curvature must be due to material rotation because interfaces cannot be finitely curved with respect to the material lattice. It is also found that the interface driving traction is not influenced by local lattice rotations unless inertia affects the reaction.     

Dynamics of a fluid sheet on a curved rigid wall with allowance for phase transformation on the free surface

The flow of a three-dimensional sheet on a curved wall is considered. Gravity and surface tension forces act on the sheet while a droplet stream falls on its free surface. The systems of equations of viscous incompressible fluid dynamics on a curved rigid surface and the boundary conditions with allowance for the falling droplet stream are formulated. The problems of steady axisymmetric motion of the sheet on cylindrical and conical surfaces are considered. The effect of the curvature of the rigid wall on the solution is examined. Kharkov. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 42–50, July–August, 1994.     

Flow pattern transition,pressure gradient,hold-up predictions in gas/non-Newtonian power-law fluid stratified flow

This work focuses on gas/non-Newtonian power-law fluid stratified pipe flow. Two different theoretical approaches to obtain pressure gradient and hold-up predictions are presented: the steady fully developed two-fluid model and the pre-integrated model. The theoretical predictions are compared with experimental data available for horizontal and for slightly downward inclined air/shear thinning fluid stratified flow taken from literature. The predictions of the pre-integrated model are validated showing a good agreement when compared with experimental data. The criteria for the transition from the stratified flow pattern are applied to gas/non-Newtonian stratified flow. The neutral stability analysis (smooth/wavy stratified flow) and the well-posedness (existence region of stratified flow) of governing equations are carry out. The predicted transition boundaries are obtained using the steady fully developed two-fluid model and the pre-integrated model, where the shape factors and their derivatives are accounted for. A comparison between the predicted boundaries and experimental flow pattern maps is presented and shows a good agreement. A comment on the shear stress modeling by the pre-integrated model is provided.     

A coupled model for simulation of the gas–liquid two-phase flow with complex flow patterns

A new model coupling two basic models, the model based on interface tracking method and the two-fluid model, for simulating gas–liquid two-phase flow is presented. The new model can be used to simulate complex multiphase flow in which both large-length-scale interface and small-length-scale gas–liquid interface coexist. By the physical state and the length scale of interface, three phases are divided, including the liquid phase, the large-length-scale-interface phase (LSI phase) and the small-length-scale-interface phase (SSI phase). A unified solution framework shared by the two basic models is built, which makes it convenient to perform the solution process. Based on the unified solution framework, the modified MCBA–SIMPLE algorithm is employed to solve the Navier–Stokes equations for the proposed model. A special treatment called “volume fraction redistribution” is adopted for the special grids containing all three phases. Another treatment is proposed for the advection of large-length-scale interface when some portion of SSI phase coalesces into LSI phase. The movement of the large-length-scale interface is evaluated using VOF/PLIC method. The proposed model is equivalent to the two-fluid model in the zone where only the liquid phase and the SSI phase are present and to the model based on interface tracking method in the zone where only the liquid phase and the LSI phase are present. The characteristics of the proposed model are shown by four problems.     

A finite volume–multigrid method for flow simulation on stratified porous media on curvilinear co‐ordinate systems

This paper presents a numerical study of infiltration processes on stratified porous media. The study is carried out to examine the performance of a finite volume method on problems with discontinuous solutions due to the transmission conditions in the interfaces. To discretize the problem, a curvilinear co‐ordinate system is used. This permits matching the interface with the boundary of the control volumes that interchange fluxes between layers. The use of the multigrid algorithm for the resulting systems of equations allows problems involving a large number of nodes with low computational cost to be solved. Finally, some numerical experiments, which show the capillary barrier behaviour depending on the material used for the different layers and the geometric design of the interface, are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Electrokinetic flow and energy conversion in a curved microtube

Curved channels are ubiquitous in microfluidic systems. The pressuredriven electrokinetic flow and energy conversion in a curved microtube are investigated analytically by using a perturbation analysis method under the assumptions of the small curvature ratio and the Reynolds number. The results indicate that the curvature of the microtube leads to a skewed pattern in the distribution of the electrical double layer (EDL) potential. The EDL potential at the outer side of the bend is larger than that at the inner side of the bend. The curvature shows an inhibitory effect on the magnitude of the streaming potential field induced by the pressure-driven flow. Since the spanwise pressure gradient is dominant over the inertial force, the resulting axial velocity profile is skewed into the inner region of the curved channel. Furthermore, the flow rate in a curved microtube could be larger than that in a straight one with the same pressure gradient and shape of cross section. The asymptotic solutions of the axial velocity and flow rate in the absence of the electrokinetic effect are in agreement with the classical results for low Reynolds number flows. Remarkably, the curved geometry could be beneficial to improving the electrokinetic energy conversion (EKEC) efficiency.     

Simulations of free surface flows with implementation of surface tension and interface sharpening in the two-fluid model

Two-fluid model used for free surface flows with large characteristic scales is improved; the smeared interface is sharpened with conservative level set method and the surface tension force with wetting angle is implemented. Surface tension force is split between two phases with several models. Detailed analysis showed the splitting of surface tension force with volume averaging as the most appropriate. The improved two-fluid model with interface sharpening and implemented surface tension is validated on several test cases. The pressure jump over a droplet interface test case showed that the pressure jump in simulation converges with grid refinement to the analytical one. The parasitic currents in simulation are one order of magnitude larger than in simulation with volume of fluid model. In the oscillating droplet test case the time period of oscillating droplet with initially ellipsoid or square shape is similar to the analytical time period. In the rising bubble test case, the rising bubble position, terminal velocity, and circularity are similar to the one observed in simulations with level set model. The wetting angle is implemented in the two-fluid model with interface sharpening and surface tension force. Model is tested in the simulation of droplet in contact with wall with different wetting angles.