Development of secondary flows in viscoelastic curved ducts and the influence of outlet region

This paper presents numerical simulations of Newtonian and viscoelastic flows through a 180° curved duct of square cross section with a long straight outlet region. A particular attention is paid to the development of the flow in the output rectangular region after the curved part. The viscoelastic fluid is modeled using the constitutive equation proposed by Phan–Thien–Tanner (PTT). The numerical results, obtained with a finite-volume method, are shown for three different Dean numbers (125,137,150)$(125,137,150)$ and for three Deborah numbers (0.1,0.2,0.3)$(0.1,0.2,0.3)$. The necessary outlet length to impose boundary conditions is presented and discussed for these cases. Streamlines and vortex formation are shown to illustrate and analyze the evolution of the secondary flow in this region.     

具有不规则自由表面的粘弹性介质中声波模拟的有限差分算法

声波在具有不规则自由表面的粘弹性介质传播的数值模拟工作在地震勘探,地震预测中非常重要,特别是当模型有起伏的自由界面和较强的衰减特性时更是如此。基于镜像方法和直接方法,本文发展了一种二维有限差分算法,可以模拟不规则自由表面引起的声波散射问题。该方法将自由表面条件与速度应力方程结合求解粘弹性波动方程,在垂直和水平自由段及拐点处设置相应的边界条件。该方法假设自由表面穿过剪切应力和相应参数的网格点。为了提高计算精度,分别计算了水平和垂直方向上应力镜像值,而对质点速度,采用先水平方向后垂直方向分析进行镜像计算和更新。将粘弹性水平自由表面镜像方法和不同倾斜度平滑自由表面的改进方法的数值结果进行对比,并验证了算法的精度。     

反冲质子磁分析技术用于氘氚中子能谱测量研究

介绍了一种基于反冲质子法和磁分析技术的氘氚聚变诊断方法, 适用于稳态及脉冲条件下的等离子体温度、燃料密度和中子产额的精确诊断. 设计了小型的原理性装置, 磁分析器使用高性能钕铁硼二极永磁铁, 焦平面上使用CR-39固体径迹探测器或PIN探测器测量质子位置分布. 使用²³⁹Pu α 源对磁分析器进行了实验标定, 建立了配套的模拟程序. 利用蒙特卡罗方法模拟分析了装置整体性能, 并在K-400加速器上进行了中子实验研究.     

Deterministic and stochastic equations of motion arising in Oldroyd fluids of order one: existence,uniqueness, exponential stability and invariant measures

Abstract

In this work, we consider the two-dimensional viscoelastic fluid flow equations, arising from the Oldroyd model for the non-Newtonian fluid flows. We investigate the well-posedness of such models in two-dimensional bounded and unbounded (Poincaré domains) domains, both in deterministic and stochastic settings. The existence and uniqueness of weak solution in the deterministic case is proved via a local monotonicity property of the linear and nonlinear operators and a localized version of the Minty-Browder technique. Some results on the exponential stability of stationary solutions are also established. The global solvability results for the stochastic counterpart are obtained by a stochastic generalization of the Minty-Browder technique. The exponential stability results in the mean square as well as in the pathwise (almost sure) sense are also discussed. Using the exponential stability results, we finally prove the existence of a unique invariant measure, which is ergodic and strongly mixing.     

Dynamic toughness in elastic nonlinear viscous solids

This work addresses the interrelationship among dissipative mechanisms—material separation in the fracture process zone (FPZ), nonelastic deformation in the surrounding background material and kinetic energy—and how they affect the macroscopic dynamic fracture toughness as well as the limiting crack speed in strain rate sensitive materials. To this end, a micromechanics-based model for void growth in a nonlinear viscous solid is incorporated into a microporous strip of cell elements that forms the FPZ. The latter is surrounded by background material described by conventional constitutive relations. In the first part of the paper, the background material is assumed to be purely elastic. Here, the computed dynamic fracture toughness is a convex function of crack velocity. In the second part, the background material as well as the FPZ are described by similar rate-sensitivity parameters. Voids grow in the strain rate strengthened FPZ as the crack velocity increases. Consequently, the work of separation increases. By contrast, the inelastic dissipation in the background material appears to be a concave function of crack velocity. At the lower crack velocity regime, where dissipation in the background material is dominant, toughness decreases as crack velocity increases. At high crack velocities, inelastic deformation enhanced by the inertial force can cause a sharp increase in toughness. Here, the computed toughness increases rapidly with crack velocity. There exist regimes where the toughness is a non-monotonic function of the crack velocity. Two length scales—the width of the FPZ and the ratio of the shear wave speed to the reference strain rate—can be shown to strongly affect the dynamic fracture toughness. Our computational simulations can predict experimental data for fracture toughness vs. crack velocity reported in several studies for amorphous polymeric materials.     

On the burst of branched polymer melts during inflation

Two molten low-density polyethylene melts, shaped as plates, have been inflated into a circular cylinder during isothermal conditions. Lowering the inflation rates allow the plates to be inflated into a larger volume of the cylinder before bursting. Numerical simulations of the inflations have been performed, using a time-strain separable constitutive K-BKZ equation based on the potential function from the Doi–Edwards theory. The material parameters in the constitutive model are based on liner viscoelastic and time dependent uniaxial elongational viscosities. The numerical calculations show quantitative agreement with the experiments, including the appearance of the burst, for a wide range of experimental conditions. This strongly suggests that the initiation of the burst in the polymer melts is a hydrodynamic phenomenon.     

The study of the mechanical impedance of foods and biomaterials to characterize their linear viscoelastic behavior at high frequencies

A new instrumental setup, which is based on the oscillatory squeezing flow model, was developed to characterize the viscoelastic properties of foods and biopolymers. Analysis of the data was performed by two different approaches. The first approach employed principles of vibration to determine viscous damping and elasticity of the sample harmonically compressed between two plates. The second approach involved the use of additional calculations based on the squeezing flow model which were linked to concepts of vibration analysis, such as mechanical impedance, to determine fundamental rheological parameters like complex viscosity and related (elastic and viscous) moduli. The experimental setup for the method is simple to use and could be attached to existing commercial instruments such as texture analyzers and universal testing machines. The use of the proposed method with this type of instruments would provide them with the additional capability of performing dynamic rheological testing. The dynamic mass of the instrument was significantly low when compared to that of other instruments that use similar principles. This low dynamic mass enabled the use of relatively higher frequencies for the testing of the samples. Comparison of the viscous and elastic moduli obtained with the proposed method and from conventional rheometers for a variety of foods and biopolymers showed good agreements.     

Large frequency bands viscoelastic properties of honey

Honey is a viscoelastic material which presents a crystallization phase transition at low temperatures. This phase transition limits the studies of dynamic behavior at high frequencies from the classical rheometry using time–temperature superposition (TTS). In order to characterize the viscoelastic properties of honey at high frequencies, we have developed a multiple ultrasonic reflection device (MUR). The viscoelastic properties of honey were measured by MUR and classical rheometry at high and low frequencies, respectively, between 13.1 and 31.3 °$°$C. Matching both results, we built the master curve over a frequency range covering nine decades, from which we determined the main rheological parameters of honey. Finally, from an inverse approach, we extracted from this master curve, the acoustical parameters for frequencies inaccessible by the ultrasonic methods.     

Onset of electroconvective instability of Oldroydian viscoelastic liquid layer in Brinkman porous medium

The problem of the onset of electrohydrodynamic instability in a horizontal layer of Oldroydian viscoelastic dielectric liquid through Brinkman porous medium under the simultaneous action of a certical ac electric field and a vertical temperature gradient is analyzed. Applying linear stability theory, we derive an equation of eight order. Under somewhat suitable boundary conditions, this equation can be solved exactly to yield the required eigenvalue relationship from which various critical values are determined in detail. Both the cases of stationary and oscillatory instabilities are discussed if the liquid layer is heated from below or above. The effects of the porosity of porous medium, the medium permeability, the Prandtl number, the ratio of retardation time to relaxation time, the elastic number, in the presence or absence of Rayleigh number are shown graphically for both cases. Some of the known results are derived as special cases. The electrical force has been shown to be the sole agency causing instability of the considered system since it is much more important than the buoyancy force even if the medium is porous.     

Entropy analysis for viscoelastic magnetohydrodynamic flow over a stretching surface

This paper presents the application of the second law analysis of thermodynamics to viscoelastic magnetohydrodynamic flow over a stretching surface. The velocity and temperature profiles are obtained analytically using the Kummer's functions and used to compute the entropy generation number. The effects of the magnetic parameter, the Prandtl number, the heat source/heat sink parameter and the surface temperature parameter on velocity and temperature profiles are presented. The influences of the same parameters, the Hartmann number, the dimensionless group parameter and the Reynolds number on the entropy generation are also discussed.