Evaluation of the drag force by integrating the energy dissipation rate in stokes flow for 2D domains using the FEM

The total drag force on the surface of a body, which is the sum of the form drag and the skin friction drag in a 2D domain, is numerically evaluated by integrating the energy dissipation rate in the whole domain for an incompressible Stokes fluid. The finite element method is used to calculate both the energy dissipation rate in the whole domain as well as the drag on the boundary of the body. The evaluation of the drag and the energy dissipation rate are post-processing operations which are carried out after the velocity field and the pressure field for the flow over a particular profile have been obtained. The results obtained for the flow over three different but constant area profiles—a circle, an ellipse and a cross-section of a prolate spheroid—with uniform inlet velocity are presented and it is shown that the total drag force times the velocity is equal to the total energy dissipation rate in the entire finite flow domain. Hence, by calculating the energy dissipation rate in the domain with unit velocity specified at the far-field boundary enclosing the domain, the drag force on the boundary of the body can be obtained.     

Frictional energy dissipation in contact of nominally flat rough surfaces under harmonically varying loads

If the nominal contact tractions at an interface are everywhere below the Coulomb friction limit throughout a cycle of oscillatory loading, the introduction of surface roughness will generally cause local microslip between the contacting asperities and hence some frictional dissipation. This dissipation is important both as a source of structural damping and as an indicator of potential fretting damage. Here we use a strategy based on the Ciavarella-Jäger superposition and a recent solution of the general problem of the contact of two half spaces under oscillatory loading to derive expressions for the dissipation per cycle which depend only on the normal incremental stiffness of the contact, the external forces and the local coefficient of friction. The results show that the dissipation depends significantly on the relative phase between the oscillations in normal and tangential load—a factor which has been largely ignored in previous investigations. In particular, for given load amplitudes, the dissipation is significantly larger when the loads are out of phase. We also establish that for small amplitudes the dissipation varies with the cube of the load amplitude and is linearly proportional to the second derivative of the elastic compliance function for all contact geometries, including those involving surface roughness. It follows that experimental observations of less than cubic dependence on load amplitude cannot be explained by reference to roughness alone, or by any other geometric effect in the contact of half spaces.     

Dissipation in a dilute suspension of spheres in a second-order fluid

We examine the effective medium properties of a dilute suspension of spheres in a second-order fluid under linear shear. Since the second-order fluid is the first step toward the general viscoelastic fluid, the results obtained may provide a qualitative feel for the problem in which the suspending fluid obeys a more complicated (and realistic) constitutive relation.The dissipation in the medium is calculated by determining the rate of working by surface forces; this is compared to the dissipation in a homogeneous fluid to give the effective properties. The results show that the term linear in volume fraction increases the corresponding rheological coefficient, just as in the Newtonian case. It is to be noted that the second-order dissipation is zero for simple shear and other weak flows, whereas for strong flows the small correction may increase or decrease the overall dissipation.     

The Development of the Asymptotic Viscous Dissipation Profile in a Vertical Free Convective Boundary Layer Flow in a Porous Medium

The effect of viscous dissipation on the development of the boundary layer flow from a cold vertical surface embedded in a Darcian porous medium is investigated. It is found that the flow evolves gradually from the classical Cheng–Minkowycz form to the recently discovered asymptotic dissipation profile which is a parallel flow.     

ENERGY DISSIPATION IN TAPPING MODE ATOMIC FORCE MICROSCOPY

原子力显微镜有多种成像模式,其中轻敲模式是最为常用的扫描方式.轻敲模式能获取样品表面形貌的高度信息和相位信息,其中相位信息具有更多的价值,如能反映样品的表面能、弹性、亲疏水性等.依据振动力学理论,相位与振动系统的能量耗散有关.探针样品间的能量耗散对于理解轻敲模式下原子力显微镜的成像机理至关重要,样品特性和测量环境会影响能量耗散.本文在不考虑毛细力影响下,基于JKR接触模型,给出了探针样品相互作用下的加卸载曲线,结合原子力显微镜力曲线实验,给出了探针-样品分离失稳点的位置,从而计算一个完整接触分离过程的能量耗散,进而讨论考虑表面粗糙度对能量耗散的影响.在轻敲模式下考虑毛细力影响,通过特征时间对比,证明挤出效应是液桥生成的主导因素,在等容条件下,用数值方法计算了不同相对湿度对能量耗散的影响.通过一维振子模型,简要说明原子力显微镜相位像与样品表面能、杨氏模量、表面粗糙度、相对湿度之间的关系.分析表明,表面粗糙度和环境湿度均会引起相位的变化,进而认为它们是引起赝像的因素.     

Steady mixed convection flow of Maxwell fluid over an exponentially stretching vertical surface with magnetic field and viscous dissipation

The steady mixed convection flow and heat transfer from an exponentially stretching vertical surface in a quiescent Maxwell fluid in the presence of magnetic field, viscous dissipation and Joule heating have been studied. The stretching velocity, surface temperature and magnetic field are assumed to have specific exponential function forms for the existence of the local similarity solution. The coupled nonlinear ordinary differential equations governing the local similarity flow and heat transfer have been solved numerically by Chebyshev finite difference method. The influence of the buoyancy parameter, viscous dissipation, relaxation parameter of Maxwell fluid, magnetic field and Prandtl number on the flow and heat transfer has been considered in detail. The Nusselt number increases significantly with the Prandtl number, but the skin friction coefficient decreases. The Nusselt number slightly decreases with increasing viscous dissipation parameter, but the skin friction coefficient slightly increases. Maxwell fluid reduces both skin friction coefficient and Nusselt number, whereas buoyancy force enhances them.     

Conjugate heat transfer in thermally developing laminar flow with viscous dissipation effects

In this study, thermally developing laminar forced convection in a pipe including viscous dissipation and wall conductance is investigated numerically. The constant heat flux is assumed to be imposed at the outer surface of the pipe wall. The finite volume method is used. The distributions for the developing temperature and local Nusselt number in the entrance region are obtained. The dependence of the results on the Brinkman number and the dimensionless thermal conductivity are shown. The viscous heating effect on the wall is shown. Significant viscous dissipation effects have been observed for large Br.     

Measurements of the dissipation coefficient at the wavemaker in the process of generation of the resonant standing waves in a tank

The complex dissipation at the wavemaker used for direct excitation of nonlinear standing waves in the vicinity of the cut-off frequency, is measured experimentally. The results indicate that the absolute value of this dissipation coefficient exceeds significantly the estimate based on purely viscous dissipation. This is attributed to the turbulent effects resulting from the vortex shedding at the discontinuities of the wavemaker surface. It is shown that incorporation of this dissipation in the boundary condition at the wavemaker as a nonlinear term in a form generally used to describe dissipation in a turbulent boundary layer, is more appropriate allowing to remove hysteresis in the value of the wavemaker dissipation coefficient. Such hysteresis is observed when a linear dissipation model is applied.     

弹塑性微凸体侧向接触相互作用能耗

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设,忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用,这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面,从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力,获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理,从而获取加、卸载法向合力与总变形之间的关系,由于法向分力产生的塑性变形及切向分力产生的摩擦,导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分,获得一个周期的微凸体接触能耗,包括应变能耗及摩擦能耗.仿真分析表明:微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大,能耗越大,且在弹性阶段最为明显.在弹性阶段,仅存在侧向的摩擦能耗,故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段,由于微凸体接触能耗为应变能耗,且接触角对其能耗影响甚微,故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型,本文提出的模型与Bartier的实验结果更吻合.     

波浪与新型双排开孔圆筒防波堤相互作用三维数值模拟

双排开孔圆筒防波堤是基于圆筒、板式结构的一种复合式新型结构型式;基于不可压缩两相流模型建立三维数值波浪水槽,通过RNG k-ε湍流模型进行湍流封闭,并采用TruVOF方法捕捉自由液面,开展波浪与双排开孔圆筒防波堤相互作用数值模拟,探究相对排间距、开孔率对新型双排开孔圆筒防波堤消浪性能的影响,分析了后排开孔圆筒防波堤附近的复杂水动力现象和流动特性.结果表明,在本文研究工况范围内,沿程平均波高随相对排间距的增大先增大后减小,随开孔率的增大而增大,周期对沿程平均波高的影响没有明显规律;当B/D=9, e=23.11%时,新型双排开孔圆筒防波堤消浪效果最优,反射系数在0.4～0.46之间,透射系数在0.3～0.35之间,耗散系数在0.8～0.85之间;自由液面破碎、水气掺混、环状涡运动演化是新型双排开孔圆筒防波堤紊动耗能消波的主要原因;相对排间距会引起后排防波堤附近涡量分布以及剪切层形态的变化,从而导致不同的紊动特性,影响双排开孔圆筒防波堤消浪特性.研究结果可以为新型双排开孔圆筒防波堤工程设计与消浪机理研究提供理论支撑.     

Strain-energy failure criterion for rigid-plastic bodies

In the framework of the theory of strain-hardening rigid-plastic bodies, we suggest a strain-energy failure criterion. The strain states of a rigid-plastic body are depicted as points in the principal strain space and form a third-order hyperbolic surface for incompressible bodies. The strain processes are depicted by lines on this surface. The loading surface related to the strain surface lines is introduced, which permits calculating the energy dissipation as a parameter determining the failure.     

Towards the numerical simulation of the horizontal slug front

Steps towards the numerical simulation of the flow behind the slug front in horizontal slug flow performed with a streamfunction-vorticity representation of the mean flow and an energy dissipation model for the turbulence are discussed. The flow field consists of two vortices, one saddle point and four stagnation regions. Attention is focused on the following boundary conditions: moving wall jet, moving wall, free jet velocity discontinuity and vertical liquid-gas open surface. A dissipation flux boundary condition is suggested to simulate the interaction of the turbulent eddies with the open surface. A method to assess the necessity to use a transport model equation for the dissipation rather than a geometric specification of a length is suggested. Three different ways to characterize the mixing zone length are proposed.     

Geometrically nonlinear continuum thermomechanics with surface energies coupled to diffusion

Surfaces can have a significant influence on the overall response of a continuum body but are often neglected or accounted for in an ad hoc manner. This work is concerned with a nonlinear continuum thermomechanics formulation which accounts for surface structures and includes the effects of diffusion and viscoelasticity. The formulation is presented within a thermodynamically consistent framework and elucidates the nature of the coupling between the various fields, and the surface and the bulk. Conservation principles are used to determine the form of the constitutive relations and the evolution equations. Restrictions on the jump in the temperature and the chemical potential between the surface and the bulk are not a priori assumptions, rather they arise from the reduced dissipation inequality on the surface and are shown to be satisfiable without imposing the standard assumptions of thermal and chemical slavery. The nature of the constitutive relations is made clear via an example wherein the form of the Helmholtz energy is explicitly given.     

Heat transfer effects in extrudate swell of elastic liquids

Numerical solutions are given which reveal the effects of viscous dissipation, shear-thinning, and elasticity on the extrudate swell of Nylon-6. Results for a cold isothermal free surface are compared with previous solutions for an adiabatic free surface.     

青藏铁路冻土区环境问题对工程安全可靠性影响

当大涡模拟用于研究化学反应流动时，传统的滤波方法会导致化学反应项不封闭. 为克服这 个困难，发展了条件滤波大涡模拟方法. 在选择适当的条件变量后，条件滤波的化学反 应项可以表达为一个封闭项. 但同时也带来了新的问题：条件滤波耗散或条件滤波扩散项的 不封闭. 为解决这一问题，采用了直接数值模拟方法研究了它们在大小尺度上的统计特 性. 研究结果表明：条件滤波耗散和扩散对于大尺度的依赖主要体现在大尺度标量场中扩散 层结构的影响，同时小尺度脉动的变化几乎与条件滤波扩散无关，而它对条件滤波耗散却显 现出明显的作用. 在构造条件滤波耗散的亚格子模型时，小尺度脉动的作用不容忽视.     

可燃颗粒悬浮流中的正激波

本文讨论了激波后可燃颗粒点火、燃烧的机制和波后流场的松弛结构,结果表明,点火阶段颗粒外层浓度梯度极大,燃烧速率由非均相化学反应速率决定,而在燃烧阶段,外层的浓度梯度极大,燃烧速率由气体的扩散速率所决定。在考虑管壁效应时,对守恒方程中的能量耗散项作了说明。     

Reflection and transmission of regular waves at a surface-pitching slotted barrier

The interactions between regular surface waves and a surface-pitching slotted barrier are investigated both analytically and experimentally.A quasi-linear theory is developed using the eigenfunction expansion method.The energy dissipation within the barriers is modeled by a quadratic friction factor,and an equivalent linear dissipation coefficient,which is depth-varying,wave-height dependent,is introduced to linearize the matching condition at the surface-pitching barrier.By comparing the theoretical results with laboratory experiments,it is shown that the present method can satisfactorily predict the variation of the reflection and transmission coefficients with wave height.     

Dynamic fracture instabilities in brittle crystals generated by thermal phonon emission: Experiments and atomistic calculations

Dynamic cleavage fracture experiments of brittle single crystal silicon revealed several length scales of surface and path instabilities: macroscale path selection, mesoscale crack deflection, and nanoscale surface ridges. These phenomena cannot be predicted or explained by any of the continuum mechanics based equations of motion of dynamic cracks, as presumably critical energy dissipation mechanisms are not fully accounted for in the theories. Experimentally measured maximum crack speed, always lower than the theoretical limit, is another phenomenon that is as yet not well understood.We suggest that these phenomena depend on velocity dependent and anisotropic material property that resists crack propagation. The basic approach is that the bond breaking mechanisms during dynamic crack propagation vibrate the atoms at the crack front to generate thermal phonon emission, or heat, which provides additional energy dissipation mechanisms. This energy dissipation mechanism is a material property that resists crack propagation. To evaluate this property, we combined the continuum based elastodynamic Freund equation of motion with molecular dynamics atomistic computer “experiments”.We analyzed the above experimental dynamic fracture instabilities in silicon with the obtained velocity dependent and anisotropic material property and show its importance in cleavage of brittle crystals.