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

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

Phase change interactions and singular fronts

The balances of mass, linear momentum and energy for a continuum provide jump relations between values of the physical variables on the two sides of a singular surface, either a boundary of the medium or an interior surface. In the case of a mixture, an overlap of interacting continua, there are jump relations for each constituent. While an elementary phase change front across which one phase of a constituent is transformed completely to a different phase can be treated as a single constituent, more general situations have co-existing phases on one side of the front, each with their own density, velocity, stress and internal energy fields, which must be treated as separate constituents. The phase change is now a mass transfer between constituents which becomes a surface production term in the mass balance jump relation for each constituent. In turn this implies surface production contributions to the momentum and energy relations associated with the surface mass transfer, including interaction body force and energy transfer contributions as well as the direct transfer terms. The general jump relations with such surface production contributions are formulated, and are illustrated for a number of situations arising in polythermal ice sheets and wet snow packs.     

Energy release,friction, and supplemental relations at phase interfaces

An Eulerian expression for the dissipation due to a moving surface of discontinuity in mass density, velocity, stress, energy, and heat flux is obtained. This leads to Eulerian measures forenergy release andfriction — measures that are work conjugate tomass flux andvelocity slip, respectively. Constitutive equations involving these quantities are proposed as a means to determine transition kinetics in materials that change phase.     

Propagation of Thermo-Elastic Waves at Several Typical Interfaces Based on the Theory of Dipolar Gradient Elasticity

The reflection and transmission properties of thermo-elastic waves at five possible interfaces between two different strain gradient thermo-elastic solids are investigated based on the generalized thermo-elastic theory without energy dissipation(the GN theory). First, the function of free energy density is postulated and the constitutive relations are defined. Then,the temperature field and the displacement field are obtained from the motion equation in the form of displacement and the thermal transport equation without energy dissipation in the strain gradient thermo-elastic solid. Finally, the five types of thermo-elastic interfacial conditions are used to calculate the amplitude ratios of the reflection and transmission waves with respect to the incident wave. Further, the reflection and transmission coefficients in terms of energy flux ratio are calculated and the numerical results are validated by the energy conservation along the normal direction. It is found that there are five types of dispersive waves, namely the coupled longitudinal wave(the CP wave), the coupled thermal wave(the CT wave), the shear wave, and two evanescent waves(the coupled SP wave and SS wave), that become the surface waves at an interface. The mechanical interfacial conditions mainly influence the coupled CP waves, SV waves, and surface waves, while the thermal interfacial conditions mainly influence the coupled CT waves.     

Viscous dissipation effects on fully developed combined free and forced non-Newtonian convection in a vertical tube

An investigation which includes the simultaneous effects of viscous dissipation and combined free and forced laminar non-Newtonian convection is presented. The problem under consideration is that of fully developed upflow in a vertical, circular tube which is heated with a constant wall heat flux. All properties are assumed to be constants in the analysis except for a temperature dependent density in the body force term which generates the free convection effects. The coupled continuity, momentum, and energy equations are solved using a finite difference technique. Numerical solutions are presented as a function of the parameters of the problem-flow behavior index n, Grashof number over Reynolds number ratio Gr/Re, and the Eckert number-Prandtl number product E Pr. The results show that heating due to viscous dissipation distorts the velocity profile, increases the friction factor, and decreases the Nusselt number.     

Direct simulation of multiphase flows with modeling of dynamic interface contact angle

We describe a modeling technique for dynamic contact angle between a phase interface and a solid wall using a generalized Navier boundary condition in the context of a front-tracking-based multiphase method. The contact line motion is determined by the generalized Navier slip boundary condition in order to eliminate the infinite shear stress at the contact line. Applying this slip boundary condition only to the interface movement with various slip ratios shows good agreement with experimental results compared to allowing full fluid slip along the solid surface. The interface slip model performs well on grid convergence tests using both the slip ratio and slip length models. A detailed energy analysis was performed to identify changes in kinetic, surface, and potential energies as well as viscous and contact line dissipation with time. A friction coefficient for contact line dissipation was obtained based on the other computed energy terms. Each energy term and the friction coefficient were compared for different grid resolutions. The effect of varying the slip ratio as well as the contact angle distribution versus contact line speed was analyzed. The behavior of drop impact on a solid wall with different advancing and receding angles was investigated. Finally, the proposed dynamic contact model was extended to three dimensions for large-scale parallel calculations. The impact of a droplet on a solid cylinder was simulated to demonstrate the capabilities of the proposing formulation on general solid structures. Widely different contact angles were tested and showed distinctive characteristic behavior clearly.     

Microstructure-performance relations of ultra-high-performance concrete accounting for effect of alpha-quartz-to-coesite silica phase transformation

The effect of the α-quartz-to-coesite silica phase transformation on the load-carrying and energy-dissipation capacities of ultra-high-performance concrete (UHPC) under dynamic loading with hydrostatic pressures of up to 10 GPa is evaluated. The model resolves essential deformation and failure mechanisms and provides a phenomenological account of the transformation. Four modes of energy dissipated are tracked, including inelastic deformation, distributed cracking, interfacial friction, and the energy dissipated through transformation of the quartz aggregate. Simulations are carried out over a range of volume fractions of the constituent phases. Results show that the phase transformation has a significant effect on the energy-dissipation capacity of UHPC for the conditions studied. Although transformation accounts for less than 2% of the total energy dissipation, the transformation leads to a twofold increase in the crack density and yields almost an 18% increase in the overall energy dissipation. Structure-response relations that can be used for materials design are established.     

Diffuse interface model for high speed cavitating underwater systems

High speed underwater systems involve many modelling and simulation difficulties related to shocks, expansion waves and evaporation fronts. Modern propulsion systems like underwater missiles also involve extra difficulties related to non-condensable high speed gas flows. Such flows involve many continuous and discontinuous waves or fronts and the difficulty is to model and compute correctly jump conditions across them, particularly in unsteady regime and in multi-dimensions. To this end a new theory has been built that considers the various transformation fronts as ‘diffuse interfaces’. Inside these diffuse interfaces relaxation effects are solved in order to reproduce the correct jump conditions. For example, an interface separating a compressible non-condensable gas and compressible water is solved as a multiphase mixture where stiff mechanical relaxation effects are solved in order to match the jump conditions of equal pressure and equal normal velocities. When an interface separates a metastable liquid and its vapor, the situation becomes more complex as jump conditions involve pressure, velocity, temperature and entropy jumps. However, the same type of multiphase mixture can be considered in the diffuse interface and stiff velocity, pressure, temperature and Gibbs free energy relaxation are used to reproduce the dynamics of such fronts and corresponding jump conditions. A general model, based on multiphase flow theory is thus built. It involves mixture energy and mixture momentum equations together with mass and volume fraction equations for each phase or constituent. For example, in high velocity flows around underwater missiles, three phases (or constituents) have to be considered: liquid, vapor and propulsion gas products. It results in a flow model with 8 partial differential equations. The model is strictly hyperbolic and involves waves speeds that vary under the degree of metastability. When none of the phase is metastable, the non-monotonic sound speed is recovered. When phase transition occurs, the sound speed decreases and phase transition fronts become expansion waves of the equilibrium system. The model is built on the basis of asymptotic analysis of a hyperbolic total non-equilibrium multiphase flow model, in the limit of stiff mechanical relaxation. Closure relations regarding heat and mass transfer are built under the examination of entropy production. The mixture equation of state (EOS) is based on energy conservation and mechanical equilibrium of the mixture. Pure phases EOS are used in the mixture EOS instead of cubic one in order to prevent loss of hyperbolicity in the spinodal zone of the phase diagram. The corresponding model is able to deal with metastable states without using Van der Waals representation.     

Impact comminution of solids due to local kinetic energy of high shear strain rate: I. Continuum theory and turbulence analogy

The modeling of high velocity impact into brittle or quasibrittle solids is hampered by the unavailability of a constitutive model capturing the effects of material comminution into very fine particles. The present objective is to develop such a model, usable in finite element programs. The comminution at very high strain rates can dissipate a large portion of the kinetic energy of an impacting missile. The spatial derivative of the energy dissipated by comminution gives a force resisting the penetration, which is superposed on the nodal forces obtained from the static constitutive model in a finite element program. The present theory is inspired partly by Grady's model for expansive comminution due to explosion inside a hollow sphere, and partly by analogy with turbulence. In high velocity turbulent flow, the energy dissipation rate gets enhanced by the formation of micro-vortices (eddies) which dissipate energy by viscous shear stress. Similarly, here it is assumed that the energy dissipation at fast deformation of a confined solid gets enhanced by the release of kinetic energy of the motion associated with a high-rate shear strain of forming particles. For simplicity, the shape of these particles in the plane of maximum shear rate is considered to be regular hexagons. The particle sizes are assumed to be distributed according to the Schuhmann power law. The condition that the rate of release of the local kinetic energy must be equal to the interface fracture energy yields a relation between the particle size, the shear strain rate, the fracture energy and the mass density. As one experimental justification, the present theory agrees with Grady's empirical observation that, in impact events, the average particle size is proportional to the (−2/3) power of the shear strain rate. The main characteristic of the comminution process is a dimensionless number B_a (Eq. (37)) representing the ratio of the local kinetic energy of shear strain rate to the maximum possible strain energy that can be stored in the same volume of material. It is shown that the kinetic energy release is proportional to the (2/3)-power of the shear strain rate, and that the dynamic comminution creates an apparent material viscosity inversely proportional to the (1/3)-power of that rate. After comminution, the interface fracture energy takes the role of interface friction, and it is pointed out that if the friction depends on the slip rate the aforementioned exponents would change. The effect of dynamic comminution can simply be taken into account by introducing the apparent viscosity into the material constitutive model, which is what is implemented in the paper that follows.     

Fluid penetration effects in porous media contact

The contact boundary conditions at the interface between two fluid-saturated porous bodies are derived. The general derivation is performed within the well-founded framework of the Theory of Porous Media (TPM) based on the constituent balance relations of mass, momentum, and energy accounting for finite discontinuities at the contact surface. Particular attention is drawn to the effects associated with the interstitial fluid flux across the interface. The derived contact conditions include two kinematic continuity conditions for the solid velocity and the fluid seepage velocity as well as two jump conditions for the effective solid stress and the pore-fluid pressure. As an application, the common case of biphasic porous media contact proceeding from materially incompressible constituents and inviscid fluid properties is discussed in detail.      

Some Thermodynamical Aspects of the Solidification of Two-phase Flows

The solidification of a two-phase flow is analyzed. It is shown how fluid-to-fluid phase transformations and kinetic stresses due to relative diffusion between fluid components influence dissipation of energy, friction and heat flux at the interface solid–fluid.     

Transient free convection with mass transfer MHD micropolar fluid in a porous plate by the network method

The transient problem of coupled heat and mass transfer of a micropolar fluid in magneto‐hydrodynamic free convection from a vertical infinite porous plate with an exponentially decaying heat generating considering the viscous dissipation and ohmic heating effects is studied. Joule heating must be considered when the viscous dissipation and the Prandtl number are large. The non‐dimensional equations for the conservation of mass, momentum, energy and concentration are solved by means a numerical technique based on electric analogy (network simulation method). This method provides the numerical response of the system by running the network in circuit resolution software with the solution to both transient and steady‐state problems at the same time, and its programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The effects of the material parameters, viscous dissipation, internal generation and Joule heating on velocity, angular momentum and temperature fields across the boundary layer are investigated. In addition, the skin‐friction coefficient, couple stress coefficient, Nusselt number and Sherwood number are shown in tabular form. The numerical results for velocity and temperature distributions of micropolar fluids are compared with the corresponding flow problems for a Newtonian fluid. Copyright © 2007 John Wiley & Sons, Ltd.     

微液滴在PDMS软基体表面的动态摩擦行为研究

通过固液界面摩擦力测试装置研究了微液滴在PDMS软基体表面运动时的动态摩擦学行为,并对微液滴体积、滑动速度及软基体力学性能对固液界面动态摩擦行为的影响进行了分析. 结果表明:微液滴在软基体表面运动时表现出最大静摩擦力和动态摩擦力. 最大静摩擦力与微液滴黏度和速度梯度呈正比,动态摩擦力与微液滴体积、滑动速度和基体力学性能有关. 随着微液滴体积的增加,三相接触线长度增加,动态摩擦力增加;随着相对滑动速度增加,三相接触线长度及接触角滞后增加,动态摩擦力增加;随着软基体弹性模量降低,固液界面黏附力增加,固液界面运动能量耗散增加,动态摩擦力增加. 研究结果可为PDMS软基体表面微液滴的精确驱动和运动参数优化提供理论指导,也可进一步丰富固液界面摩擦理论.      

Slip weakening in rocks and analog materials at co-seismic slip rates

Determination of co-seismic slip resistance in earth faults is critical for understanding the magnitude of shear-stress reduction and hence the near-fault acceleration that can occur during earthquakes. Also, knowledge of shear resistance dependency on slip velocity, slip distance, normal stress, and surface roughness is fundamental information for understanding earthquake physics and the energy released during such events. In the present study, plate-impact pressure-shear friction experiments are employed to investigate the frictional resistance in fine-grained Arkansas novaculite rock and an analog material (soda-lime glass), under relevant interfacial conditions. The results of the experiments indicate that a wide range of frictional slip conditions can exist during a single slip event. These conditions range from initial no-slip followed by slip weakening, strengthening, and seizure at the frictional interface. For the case of glass vs. glass experiments, the first slip-weakening, with μ in the range of 0.4-0.2, is understood to be most likely due to thermal-induced flash heating and incipient melting at asperity junctions, while the slip-strengthening, with μ in the range of 0.1-0.4, is understood to be a result of coalescence and solidification of local melt patches. For the case of the fine grain Arkansas Novaculite rocks a similar range of slip conditions is observed.     

Interface propagation and microstructure evolution in phase field models of stress-induced martensitic phase transformations

Analytical solutions for diffuse interface propagation are found for two recently developed Landau potentials that account for the phenomenology of stress-induced martensitic phase transformations. The solutions include the interface profile and velocity as a function of temperature and stress tensor. An instability in the interface propagation near lattice instability conditions is studied numerically. The effect of material inertia is approximately included. Two methods for introducing an athermal interface friction in phase field models are discussed. In the first method an analytic expression defines the location of the diffuse interface, and the rate of change of the order parameters is required to vanish if the driving force is below a threshold. As an alternative and more physical approach, we demonstrate that the introduction of spatially oscillatory stress fields due to crystal defects and the Peierls barrier, or to a jump in chemical energy, reproduces the effect of an athermal threshold. Finite element simulations of microstructure evolution with and without an athermal threshold are performed. In the presence of spatially oscillatory fields the evolution self-arrests in realistic stationary microstructures, thus the system does not converge to an unphysical single-phase final state, and rate-independent temperature- and stress-induced phase transformation hysteresis are exhibited.     

Active and passive controls of nanoparticles in Maxwell stagnation point flow over a slipped stretched surface

A steady stagnation-point flow of an incompressible Maxwell fluid towards a linearly stretching sheet with active and passive controls of nanoparticles is studied numerically. The momentum equation of the Maxwell nanofluid is inserted with an external velocity term as a result of the flow approaches the stagnation point. Conventional energy equation is modified by incorporation of nanofluid Brownian and thermophoresis effects. The condition of zero normal flux of nanoparticles at the stretching surface is defined to impulse the particles away from the surface in combination with nonzero normal flux condition. A hydrodynamic slip velocity is also added to the initial condition as a component of the entrenched stretching velocity. The governing partial differential equations are then reduced into a system of ordinary differential equations by using similarity transformation. A classical shooting method is applied to solve the nonlinear coupled differential equations. The velocity, temperature and nanoparticle volume fraction profiles together with the reduced skin friction coefficient, Nusselt number and Sherwood number are graphically presented to visualize the effects of particular parameters. Temperature distributions in passive control model are consistently lower than in the active control model. The magnitude of the reduced skin friction coefficient, Nusselt number and Sherwood number decrease as the hydrodynamic slip parameter increases while the Brownian parameter has negligible effect on the reduced heat transfer rate when nanoparticles are passively controlled at the surface. It is also found that the stagnation parameter contributes better heat transfer performance of the nanofluid under both active and passive controls of normal mass flux.     

An Analysis of Turbulent Motions In and Around a Differentially Forced Density Interface

The Rapid-Distortion-Theory-based analysis proposed by Fernando and Hunt [1] is extended to study the nature of turbulence in and around a density interface sandwiched between turbulent layers with dissimilar properties. It is shown that interfacial motions consist of low-frequency, resonantly excited, nonlinear internal waves and high-frequency, linear internal waves driven by background turbulence. Based on the assumptions that (i) all resonant waves and some nonresonant waves having frequencies close to the resonant frequencies grow rapidly, break, and cause interfacial mixing, (ii) the spectral amplitude of the vertical velocity in the wave-breaking regime is constant, and (iii) kinetic energy is equipartitioned between linear and nonlinear breaking wave regimes, the r.m.s. vertical velocity at the interface and the turbulent kinetic energy flux into the interface are calculated. The migration velocity of the interface is calculated using the additional assumption that the buoyancy flux into a given turbulent layer is a fixed fraction of the turbulent kinetic energy flux supplied to the interface by the same layer. The calculations are found to be in good agreement with the entrainment data obtained in previous laboratory experiments in the parameter regime where the interface is dominated by internal wave dynamics. Received 23 July 1997 and accepted 8 January 1999     

A compressible two‐fluid multiphase model for CO2 leakage through a wellbore

This paper introduces an effectively mesh‐independent and computationally efficient model for CO₂ leakage through wellbores. A one‐dimensional compressible two‐fluid domain, representing a homogeneous air gas and a multiphase CO₂ with a jump at the interface between them, is modeled. The physical domain is modeled using the drift‐flux model, and the governing equations are solved using a mixed finite‐element discretization scheme. The standard Galerkin FEM, the partition of unity method, and the level‐set method are integrated to solve the problem. All important physical phenomena and processes occurring along the wellbore path, including fluid dynamics, buoyancy, phase change, compressibility, thermal interaction, wall friction, and slip between phases, together with the jump in density and enthalpy between air and CO₂, are considered. Two numerical examples illustrating the computational capability and efficiency of the model are presented. Copyright © 2015 John Wiley & Sons, Ltd.     

Photoelastic studies on progress of separation in interference fits

Interference fits are used extensively in aircraft structural joints because of their improved fatigue performance. Recent advances in analysis of these joints have increased understanding of the nonlinear load-contact and load-interfacial slip variations in these joints. Experimental work in these problems is lacking due to difficulties in determining partial contact and partial slip along the pin-hole interface. In this paper, an experimental procedure is enumerated for determining load-contact relations in interference/clearance fits, using photoelastic models and applying a technique for detecting progress of separation/contact up to predetermined locations. The study incorporates a detailed procedure for model making, controlling interference, locating break of contact up to known locations around the interface, estimating optically the degree of interference, determining interfacial friction and evaluating stresses in the sheet. Experiments, simulating joints in large sheets, were carried out under both pin and plate loads. The present studies provide load-separation behavior in interference joint with finite interfacial friction.     

Thermomechanics of the interface between a body and its environment

We formulate integral statements of force balance, energy balance, and entropy imbalance for an interface between a body and its environment. These statements account for interfacial energy, entropy, and stress but neglect the inertia of the interface. Our final results consist of boundary conditions describing thermomechanical interactions between the body and its environment. In their most general forms, these results are partial differential equations that account for dissipation and encompass as special cases Navier’s slip law, Newton’s law of cooling, and Kirchhoff’s law of radiation. When dissipation is neglected, our results reduce to the well-known zero-slip, free-surface, zero-shear, prescribed temperature, and flux-free conditions. Dedicated to James K. Knowles: teacher, colleague, friend