A two-parameter model of wave regimes for viscous liquid film flows

Wave regimes of viscous liquid film flows are considered when the viscosity coefficients vary in a wide range. An approximate model system of differential equations with two external governing parameters for the film layer thickness and the local flow rate is derived. The viscous dissipation of a film layer is taken into account in this system more accurately than in the well-known one-parameter Shkadov model. New properties of linear and nonlinear waves caused by the hydrodynamic instability of high-viscous liquid flows under gravity and surface tension are found.     

Effective and efficient characterization of lubrication flow over soft coatings

This work models the fluid-structure interactions associated with separating a solid from a soft elastic film in a liquid environment. One side of the soft film is permanently attached to a rigid substrate. The ensuing liquid flow and elastic deformations are derived by considering a system of partial differential equations, that govern the mechanics of the separation process. A finite element based computational scheme is developed to solve the governing equations and predict the resultant forces acting on the solid. It is shown that the resulting forces are influenced by the elasticity of the film for an initial duration and by the viscosity of the liquid at larger times. The proposed model is utilized to shed insights into the mechanics of the separation process in constrained-surface stereolithography process.

     

Development of a PLIC-VOF method for the dynamic simulation of entry region flow in a laminar falling film

The present work describes a numerical procedure to simulate the development of hydrodynamic entry region in a gravity-driven laminar liquid film flow over an inclined plane. It provides a better insight into the physics of developing film in entry region. A novel numerical approach is proposed which has the potential to provide solutions for the complex physics of liquid film spreading on solid walls. The method employs an incompressible flow algorithm to solve the governing equations, a PLIC-VOF method to capture the free surface evolution and a continuum surface force (CSF) model to include the effect of surface tension. To account for the moving contact line on the solid substrate, a precursor film model based wall treatment is implemented. Liquid film flow has been simulated for the Reynolds number range of 5 ≤ Re ≤ 37.5, and the predicted results are found to agree well with the available analytical and experimental data.     

Numerical modeling of elastic waves in micropolar plates and shells taking into account inertial characteristics

Mathematical models of micropolar plates and shells are considered within the framework of the approximation approach. The governing equations of the theories are written in a thermodynamically consistent form of the conservation laws. This ensures hyperbolicity and correctness of the initial boundary value problems. For numerical solution, we propose parallel algorithms for supercomputers with graphics processing units. The algorithms are based on the splitting method with respect to spatial variables. We present the results of numerical computations of wave propagation in micropolar rectangular plates and cylindrical panels for media with different types of microstructure particles.     

An implicit three‐dimensional numerical model to simulate transport processes in coastal water bodies

A three‐dimensional baroclinic numerical model has been developed to compute water levels and water particle velocity distributions in coastal waters. The numerical model consists of hydrodynamic, transport and turbulence model components. In the hydrodynamic model component, the Navier–Stokes equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. The transport model component consists of the pollutant transport model and the water temperature and salinity transport models. In this component, the three‐dimensional convective diffusion equations are solved for each of the three quantities. In the turbulence model, a two‐equation k–ϵ formulation is solved to calculate the kinetic energy of the turbulence and its rate of dissipation, which provides the variable vertical turbulent eddy viscosity. Horizontal eddy viscosities can be simulated by the Smagorinsky algebraic sub grid scale turbulence model. The solution method is a composite finite difference–finite element method. In the horizontal plane, finite difference approximations, and in the vertical plane, finite element shape functions are used. The governing equations are solved implicitly in the Cartesian co‐ordinate system. The horizontal mesh sizes can be variable. To increase the vertical resolution, grid clustering can be applied. In the treatment of coastal land boundaries, the flooding and drying processes can be considered. The developed numerical model predictions are compared with the analytical solutions of the steady wind driven circulatory flow in a closed basin and of the uni‐nodal standing oscillation. Furthermore, model predictions are verified by the experiments performed on the wind driven turbulent flow of an homogeneous fluid and by the hydraulic model studies conducted on the forced flushing of marinas in enclosed seas. Copyright © 2000 John Wiley & Sons, Ltd.     

Non-isothermal flow of a Bingham fluid with pressure and temperature dependent viscosity

In this paper we investigate the non-isothermal flow of a Bingham fluid whose viscosity and yield stress depend on temperature and pressure. We consider two situations: in the first one we assume that the buoyancy effects are dominant and influence the development and evolution of the unyielded plug. In this case the governing equations are obtained via the Oberbeck–Boussinesq approximation which is derived using a perturbative approach. We show that within this approximation the heat generated by viscous friction can be safely neglected. In the second situation we assume that the frictional heating effects are dominant and influence the flow via the viscosity and yield stress that depend on temperature. For both situations we investigate the simple unidirectional flow between plates subjected to given thermal conditions. We derive the equations for the steady fully developed flow and we determine the exact position of the yield surfaces separating the yielded and the unyielded domain. We also show some plots to assess the effects due to the dependence of the rheological parameters on the temperature and pressure.     

Film boiling on a hot body moving in a fluid

We investigate the transient film boiling in the vicinity of a stagnation point on the frontal surface of a very hot blunt body which moves with a constant velocity in an incompressible viscous fluid in the presence of a vapor layer near the body surface. Within the unsteady two-phase boundary layer approximation, the equations of motion of the liquid and vapor phases are formulatedwith account of the conservation of mass, momentum, and energy on the a priori unknown phase interface. In the vicinity of the stagnation point on the body surface, the solution of the boundary layer equations is sought in the form of series in the longitudinal coordinate. For the leading terms of the series, a parabolic system of partial differential equations is obtained, which is solved numerically. The similarity parameters controlling the film boiling process are determined. On the basis of parametric numerical calculations, the dynamics of the vapor layer are investigated for the case of a plane hot body moving in water with the room pressure and temperature. In the space of governing parameters, the limits of the existence of steady and unsteady film boiling regimes are found.     

Heat transfer in non-Newtonian falling liquid film on a horizontal circular cylinder

 This study aims to investigate numerically the laminar flow and heat transfer in a pseudoplastic non-Newtonian falling liquid film on a horizontal cylinder for the constant heat flux and isothermal boundary conditions. The inertia terms are taken into account. An implicit finite difference method is carried out to solve the governing boundary layer equations. The effects of operational parameters on the hydrodynamic and heat transfer characteristics are examined and discussed in detail. The results presented show that the local and average Nusselt numbers varies significantly as a function of the concentration of aqueous carboxymethylcellulose (CMC) solutions and the cylinder diameter. Higher concentration of aqueous CMC solutions generate larger heat transfer coefficients. Finally, a comparison with the experimental and numerical results available in the literature for Newtonian fluids shows clearly that the present analysis is reasonably accurate. Received on 29 March 2001 / Published online: 29 November 2001     

Integral treatment of subcooled forced convection film boiling on a flat plate

Subcooled forced convection film boiling on a flat plate has been analysed by means of an integral method. Following the two phase boundary layer theory, the momentum and energy equations for both liquid and vapor layers are considered along with the compatibility conditions on the liquid-vapor interface. Subsequently, the governing equations are reduced to a set of algebraic equations which can readily be solved for given parameters. Comparison of the present solution with the Cess and Sparrow solution reveals an excellent performance of the present solution procedure. The effects of superheating, subcooling and liquid Prandtl number on the hydrodynamic and heat transfer characteristics are fully discussed. Furthermore, the asymptotic formulas are derived for the local Nusselt number and skin friction coefficient through a careful examination of the physical limiting conditions.     

Nonlinear waves in a liquid film on a near-horizontal surface

Nonlinear waves in a liquid film on a slightly inclined rigid plane are studied. A mathematical model is reduced to a system of two evolutionary equations for the layer thickness and the local fluid mass flow. In addition to viscous forces, gravity, and surface tension, the pressure difference over the layer thickness, induced by the gravity force projection on the normal to the underlying surface, is also taken into account. Spatially periodic solutions developing with time from small initial disturbances into regular nonlinear waves are considered. A spectral representation of the solution, the Galerkin method with respect to the uniform coordinate, and subsequent numerical calculation of the corresponding dynamic system on large time intervals are employed. Different variants in the space of the three governing parameters are calculated and some basic mechanisms of nonlinear dynamics of the two-dimensional waves are detected. The calculation results are compared with the existing experimental data. It is shown that the theoretical conclusions can be used to interpret and predict experiments.     

Leveling and thixotropic characteristics of concentrated zirconia inks for screen-printing

Screen-printing is a cost-effective method for the mass manufacture of zirconia-based solid oxide fuel cells (SOFCs) and oxygen separation membranes. The present work outlines an investigation into the leveling, thixotropic, and screen-printing characteristics of concentrated zirconia inks by a variety of rheological and imaging methods. A combination of viscosity, shear rate jump experiments, creep and recovery analysis, and yield stress measurements were used to assess ink thixotropy. Oscillatory rheometry and scanning electron microscopy/optical microscopy revealed a consistent effect of ethyl cellulose (binder) content upon the thixotropic and leveling characteristics of zirconia inks. While the yield stress (τ ₀), extent of recovery R(%), and rate of recovery (K) increase with increasing binder content, so did the surface roughness and thickness of the screen-printed films. Increasing the binder content not only increases the network strength of the thick films but also leads to increased leveling time. As a result, rheological modifiers are proposed to be necessary to improve the leveling characteristics of zirconia inks without losing the green strength of the thick films.     

Convective heat and mass transfer from a falling film to a laminar gas stream

A numerical study has been made of convective heat and mass transfer from a falling film to a laminar gas stream between vertical parallel plates. The effects of gas-liquid phase coupling, variable thermophysical properties, and film vaporization have been considered. Simultaneous mass, momentum and heat transfer between liquid film and gas stream is numerically studied by solving the respective governing equations for the liquid film and gas stream together. The influences of the inlet liquid temperature and liquid flowrate on the cooling of liquid film are examined for air-water and air-ethanol systems. Results show that the heat transfer from the gas-liquid interface to the gas stream is predominantly determined by the latent heat transfer connected with film evaporation. Additionally, better liquid film cooling is noticed for the system having a higher inlet liquid temperature or a lower liquid flowrate.     

可记录光盘染料旋涂工艺的流动机制研究

染料旋涂是可记录光盘复制中的核心工艺，依据实验给出旋涂工艺的简化模型，并运用旋转圆盘理论分析流体流动，考虑基片上的预刻槽对流动的影响，考虑基片表面和流体间的滑移效应，导出运动方程、边界方程和连续方程，最终导出染料旋涂的控制方程。     

Convective heat exchange of a viscoplastic fluid with temperature-dependent rheological characteristics during structural flow in a circular pipe

Heat exchange in a viscoplastic liquid moving in a circular pipe is investigated, taking into account the dependence of plastic viscosity and ultimate shear stress on temperature. A system of motion, energy, and continuity equations transformed under the assumption that the Pe and Pr numbers are much greater than 1 is solved on a computer by the method of finite differences using iterations. Results of the numerical solutions for the exponential form of the dependences of the rheological characteristics on temperature are analyzed in detail. A comparison of the numerical solutions with well-known theoretical solutions in particular cases and also with experimental data indicates their high precision.     

Crossflow of viscoplastic materials through tube bundles

The flow of viscoplastic materials through staggered arrays of tubes is analyzed. The mechanical behavior of the materials is assumed to obey the generalized Newtonian liquid (GNL) model, with a viscosity function given by the biviscosity law. The governing equations of this flow are solved numerically using a finite-volume method with a non-orthogonal mesh. For a representative range of the relevant parameters, results are presented in the form of velocity, pressure and viscosity fields. The pressure drop is also given as a function of rheological and geometric parameters.     

A rheological theory for liquid crystal thin films

A macroscopic rheological theory for compressible isothermal nematic liquid crystal films is developed and used to characterize the interfacial elastic, viscous, and viscoelastic material properties. The derived expression for the film stress tensor includes elastic and viscous components. The asymmetric film viscous stress tensor takes into account the nematic ordering and is given in terms of the film rate of deformation and the surface Jaumann derivative. The material function that describes the anisotropic viscoelasticity is the dynamic film tension, which includes the film tension and dilational viscosities. Viscous dissipation due to film compressibility is described by the anisotropic dilational viscosity. Three characteristic film shear viscosities are defined according to whether the nematic orientation is along the velocity direction, the velocity gradient, or the unit normal. In addition the dependence of the rheological functions on curvature and film thickness has been identified. The rheological theory provides a theoretical framework to future studies of thin liquid crystal film stability and hydrodynamics, and liquid crystal foam rheology. Received: 9 October 2000 Accepted: 6 April 2001     

进动圆筒内液体流动不稳定的非线性特征

在建立进动充液圆筒内液体偏差流动方程的基础上，结合液体惯性波和轴向二次流动线性解，通过对定常二次流动的线性稳定性分析，提出了函数空间表达的流动不稳定性非线性分岔分析方程. 对非惯性坐标系下液体流动的Navier-Stokes方程进行了数值求解，并对惯性波发生破裂(实验提供的3种主模态下得出的共振破裂现象)时的压力时间序列进行分析，得出了液体流动不稳定的基本非线性特征.     

Hydrodynamic resistance of a spheroidal particle with uniform internal heat release

The Stokes approximation is used to describe the stationary motion of a heated hydrosol spheroidal particle in a viscous incompressible liquid in which internal, uniformly distributed heat sources (sinks) of constant capacity act. It was assumed that the average particle surface temperature could differ significantly from the temperature of the ambient liquid. An analytical expression for the hydrodynamic force acting on the uniformly heated spheroidal particle was obtained by solving hydrodynamic equations with the temperature dependence of the viscosity represented as an exponential power series.     

Rheology of a dilute suspension of axisymmetric Brownian particles

Explicit results are presented for the complete rheological properties of dilute suspensions of rigid, axisymmetric Brownian particles possessing fore-aft symmetry, when suspended in a Newtonian liquid subjected to a general three-dimensional shearing flow, either steady or unsteady. It is demonstrated that these rheological properties can be expressed in terms of five fundamental material constants (exclusive of the solvent viscosity), which depend only upon the sizes and shapes of the suspended particles. Expressions are presented for these scalar constants for a number of solids of revolution, including spheroids, dumbbells of arbitrary aspect ratio and long slender bodies. These are employed to calculate rheological properties for a variety of different shear flows, including uniaxial and biaxial extensional flows, simple shear flows, and general two-dimensional shear flows. It is demonstrated that the rheological properties appropriate to a general two-dimensional shear flow can be deduced immediately from those for a simple shear flow. This observation greatly extends the utility of much of the prior Couette flow literature, especially the extensive numerical calculations of Scheraga et al. (1951, 1955).The commonality of many disparate results dispersed and diffused in earlier publications is emphasized, and presented from a unified hydrodynamic viewpoint.