A finite volume technique to simulate the flow of a viscoelastic fluid

Hydrodynamically developing flow of Oldroyd B fluid in the planar die entrance region has been investigated numerically using SIMPLER algorithm in a non-uniform staggered grid system. It has been shown that for constant values of the Reynolds number, the entrance length increases as the Weissenberg number increases. For small Reynolds number flows the center line velocity distribution exhibit overshoot near the inlet, which seems to be related to the occurrence of numerical breakdown at small values of the limiting Weissenberg number than those for large Reynolds number flows. The distributions of the first normal stress difference display clearly the development of the flow characteristics from extensional flow to shear flow.List of symbols D rate of strain tensor - L slit halfheight - P pressure, indeterminate part of the Cauchy stress tensor - R the Reynolds number - t time - U average velocity in the slit - u velocity vector - u,v velocity components - W the Weissenberg number based on the difference between stress relaxation time and retardation time - W ₁ the Weissenberg number based on stress relaxation time - x,y rectangular Cartesian coordinates - ratio of retardation time to stress relaxation time - zero-shear-rate viscosity, ₁ + ₂ - ₁ non-Newtonian contribution to - ₂ Newtonian contribution to - ₁ stress relaxation time - ₂ retardation time - density - (, , ) xx, yy and xy components of ₁, respectively - determinate part of the Cauchy stress tensor - ₁ non-Newtonian contribution to - ₂ Newtonian contribution to      

Thermodynamic properties of antiferromagnetic superconductors containing nonmagnetic,magnetic, and spin-orbit impurities

Effects of all types of impurities (nonmagnetic, magnetic, and spin-orbit) on an antiferromagnetic superconductor (AFSC) have been investigated by studying the transition temperatureT _c and the specific heat jump. We have assumed a one-dimensional electron band. The impurity scattering is treated within the self-consistent Born approximation. We find that: (a) the molecular fieldH _Q and the magnetic impurities depress superconductivity of AFSC and their pair-breaking effect is additive; (b) the effect of spin-orbit impurities is the same as that of nonmagnetic impurities—these enhance superconductivity by screening the molecular field; and (c) in the extreme dirty limit, the AFSC is described in terms of an effective pair-breaking parameter given by 1/_eff=1/₂+H _Q ² where 1/=1/₁+2/3_so(1/₁, 1/₂, and 1/_so, respectively, are the scattering rates from nonmagnetic, magnetic and spin-orbit impurities).     

Effect of the elastic factor on the hydrodynamic stability of a structurally viscous medium

The effect of relaxation phenomena on the hydrodynamic stability of the plane gradient flow of a structurally viscous medium is investigated using linear theory.Notation _ij stress tensor deviator - U_i components of the velocity vector - x_i coordinates - t time - P pressure - =₀L/^*V plasticity parameter - _o limiting shear stress - andc dimensionless wave number and the perturbation frequency - Re=VL/* Reynolds number - density - F_ij deformation rate tensor Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 35, No. 5, pp. 868–871, November, 1978.     

The Interplay of Electronic and Nuclear Magnetism in PtFe x at Milli-, Micro-, and Nanokelvin Temperatures

We have measured ac susceptibility, nuclear magnetic resonance, and nuclear heat capacity of two PtFe _x samples with concentrations of magnetic impurities x = 11 ppm and 41 ppm at magnetic fields (0 ± 0.05) mTB248 mT. The susceptibility data have been measured at temperatures of 0.3 KT100 mK, no hint for nuclear magnetic ordering could be detected to a temperature of 0.3 K. The nuclear heat capacity data taken at 1.4 KT10 mK show enhanced values which scale with x at low polarization. This effect is described by a model assuming an internal magnetic field caused by the impurities. No indication for nuclear magnetic ordering could be detected to 1.4 K. The nuclear magnetic resonance experiments have been performed on these samples at 0.8 KT0.5 mK and 2.5 mTB22.8 mT as well as on three other samples with x = 5, 10, 31 ppm in a different setup at 40 KT0.5 mK and at 5.4 mTB200 mT. Spin-lattice and effective spin-spin relaxation times ₁and ₂ ^* of ¹⁹⁵ Pt strongly depend on x and on the external magnetic field. No temperature dependence of ₁and ₂ ^* could be detected and the NMR data, too, give no hint for nuclear magnetic ordering to 0.8 K.     

Molar momentum and heat transfer

Universal relations governing the molar transfer of momentum and heat are derived on the basis of a hypothesis about the dependence of the boundaries of the molar transfer region on the flow structure and with the use of a special mathematical transformation.Notation u average longitudinal velocity, m/sec - T average temperature, °K - T_w wall temperature, °K - kinematic viscosity coefficient, m²/sec - density, kg/m³ - c_p specific heat, J/kg·K - tangential stress, N/m² - t_w tangential stress at wall, N/m² - q_w specific heat flux at wall, W/m² - u_*=_w/ dynamic velocity, m/sec - *=q_w/c_pu_* characteristic temperature, °K - thickness of boundary layer, m - ₀ thickness of laminar sublayer, m - l = /u transverse space scale of average mole at wall, m - y⁺ = y/l _22C6; dimensionless coordinate - u⁺=u/u* dimensionless velocity - ⁺=(T_w – T)/* dimensionless temperature - ⁺=/_w dimensionless tangential stress - R=In (y⁺/ _o ⁺ )/In (⁺/ _o ⁺ ) generalized dimensionless co-ordinate - U = (u⁺ - u _o ⁺ )/(u _o ⁺ - u _o ⁺ ) generalized dimensionless velocity - Pr Prandtl number Indices * flow parameters evaluated at y⁺=1 - parameters at y⁺=⁺ - 0 parameters at y⁺= _o ⁺ - w parameters at wall Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 41, No. 3, pp. 441–448, September, 1981.     

Signal rise time of the magnetic bolometer

We present measurements of the temperature dependent signal rise time _S and discuss a model for calculating _S(T). The bolometer consists of a paramagnetic sample and an absorber. The lattice is heated up by absorbing - particles, and the relaxation of the magnetization is measured with a SQUID. With decreasing temperatures _S first increases as ₁1/T, but then decreases strongly. At 30 mK it is reduced by orders of magnitude compared with ₁. This result is in agreement with a theoretical model which takes into account the heat capacities of the lattice, the resonant phonons, the spins, and thermal resistances between these capacities. Under the condition of the bottleneck effect _S is found to be proportional to T³. At low temperatures the lowest values of _S of 2 ms may already be limited by the Kapitza resistance. These are the first measurements of the spin-lattice relaxation times with the energy being transferred from the lattice to the spins.     

The critical curve in an antiferromagnetic superconductor with nonmagnetic impurities

The critical curve of a transition of the second kind in an antiferromagnetic superconductor (AFS) with nonmagnetic impurities has been studied. The AFS is described by using the mean-field model given by Nass, Levin, and Grest and assuming a one-dimensional electron band. We find that the points on the critical curve satisfy the thermodynamic stability condition for 0₁/₀5.04 and 0.49H_Q/₀1.64.Here ₁ is the inverse lifetime of a conduction electron for nonmagnetic impurity scattering,H _Q is the antiferromagnetic molecular field, ₀ is the zero-temperature order parameter of a superconductor in the absence ofH _Q and impurities. Further, ₁ and H_Q denote the values of these quantities for points on the critical curve. For ₁/₀>5.04 and H_Q/₀>1.64, the phase transition from the superconducting to the normal state is always of the second kind. Some thermal properties of the system near the critical curve have also been investigated and we find that these depends dramatically on the impurity concentration.     

Oscillatory behavior of photogenerated carrier diffusion in γ-Mo4O11 under magnetic field

The pulsed-laser induced transient thermoelectric effect (TTE) under electric and magnetic fields has been measured for the CDW material -Mo ₄O₁₁. The photo-induced TTE voltages decay with the characteristic relaxation times _i (i=2 and 3) for thermal diffusion of two types of holes. In the CDW state [Tc (=100 K)], the relaxation times _i oscillate periodically with the magnetic field, and the corresponding relaxation amplitudes a_i change stepwise at the minima in the _i-Bcurves. These oscillation periods B_i (=0.56–0.68 T) are independent of the electric field E, and these increase appreciably with temperature above 40 K, while the initial phase differences _i increase linearly with E.Possible mechanisms for such anomalous oscillations of the relaxation times are discussed from classical and conventional quantum mechanical viewpoints; however, its origin is uncertain at the moment.     

Theory of periodically driven,current-carrying,superconducting filaments. I. Spatially homogeneous states

Starting from the nonequilibrium theory of dirty superconductors in the Ginzburg-Landau regime, spatially homogeneous states with an applied currentI=I ₀+I ₁ cos (t) are considered. Expressions for the linear response (I₁ small) valid up to high frequencies (k _BT_c) are derived and evaluated analytically for the experimentally important case of smallI ₀ and ₀(T). Then the nonlinear response is treated for frequencies with _E1. Interesting new behavior is found for frequencies ₀ 1, where ₀ is essentially the GL relaxation time.     

Transient heat conduction in hollow spheres with a moving inner boundary

The finite integral transform method is used to obtain the solution of unsteady heat conduction problems for a hollow sphere with a moving internal boundary and various boundary conditions at the outer surface. For the solution of the problems of interest integral transform formulas are presented with kernels (16), (20), and (24) and the corresponding inversion formulas (18), (22), (26), (29) and characteristic equations (17), (21), (25), (28), (31), (33).Nomenclature a, thermal diffusivity and conductivity - t temperature of phase transformation - density - heat transfer coefficient - Q total quantity of heat passing through inner boundary - F latent heat of phase transformation - Fo(1,)=a/R ₁ ² , Fo(i,)=/r _i ² , Fo(i, _i)=a _i/r _i ² Fourier numbers - Bi₂=R₂/ Biot number     

Numerical calculation of a generalized thermal model of radio-electronic apparatus

A method is proposed for numerical calculation of the temperature field of a generalized model of electronic equipment with high component density.Notation x,y,z,x,y spatial coordinates, m - time, sec - L_x, L_v, L_z dimensions of heated zone, m - _x, _y, _z effective thermal-conductivity coefficients of heated zone, W/m·deg - ₂ thermal conductivity of chassis, W/m·deg - a _z thermal diffusivity of heated zone along z axis, m²/sec - c₁ effective specific heat of heated zone, J/kg·deg - ₁ effective density of heated zone, kg/m³ - c₃, ₃, c₂, ₂ thermophysical characteristics of cooling agent and chassis, J/kg·deg·kg/m³ - q_v(x, ), q(x, y) volume heat-source distribution, W/m³ - q_s (x) surface heat-source distribution, W/m² - p number of cooling agent channels - Fo Fourier number - Bi Biot number - Uⁱ coolant velocity in i-th channel, m/sec - T₁(x, ), T₂(x, ), T₃(x, ) temperature distribution of heated zone, chassis, and coolant, °K - T₃₀, T₁₀(x), T₂₀(x) initial temperatures, °K - T_3in coolant temperature at input to channel, °K - T_T(x) effective temperature distribution of heat loss elements, °K - T_C temperature of external medium, °K - dimensionless heated zone temperature - _v(x) local volume heat exchange coefficient, W/m³·deg - ₁₂(x), _1C(x), _1T(x) heat liberation coefficients - W/m²·sec; ₂₁(x, y), _2c(x, y), _2T(x, y) volume heat-exchange coefficients of chassis with heated zone, medium, and cooling elements, W/m³·deg Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 40, No. 5, pp. 876–882, May, 1981.     

Low-Frequency Relaxation Rate in the Superconducting YBa2Cu3O6.7

The submillimeter-wave 3 cm ^–1 < < 40 cm ^–1 complex conductivity of the reduced YBa ₂ Cu ₃ O _6.7 film, (T _C =56.5 K) was investigated for temperatures 4 K < T < 300 K. The frequency dependence of the effective quasiparticle scattering rate 1/*() was extracted from the spectra. 1/* is shown to be frequency independent at low frequencies and high temperatures. On decreasing temperature the scattering rate increases with increasing frequencies. Finally, at 6 K it follows a power-law, 1/* ^1.75±0.3.     

Plasticity of lead single crystals in the superconducting and normal states at 4.2K

The flow stress increment _SN associated with a superconducting-normal phase transition has been studied as a function of deformation and strain rate for several lead single crystals at 4.2 K. In the linear work hardening range, the total applied stress _S determines the size of _SN. The dislocation inertial model recently proposed by Granato for the enhanced plasticity of the superconducting state is evaluated in detail with respect to the dependence of _SN on the distance between obstacles to dislocation movement. Reasonable agreement between theory and experimental results is obtained for the linear work hardening range, if it is assumed that the obstacles are forest dislocations. For the normal-state electronic drag coefficient a valueB _eN1.1×10^–4 dyn · sec/cm² is estimated.The experimental part of this work was performed at Materials Science Division, Argonne National Laboratory, Argonne, Illinois, under the auspices of the U.S. Atomic Energy Commission.     

Hydrodynamics of a vortex ring

We considered the kinematics and dynamics of a vortex ring in an incompressible fluid in toroidal coordinates. We obtained the change in the pressure difference along the boundary between two flow regions in the case of a moving torus.Notation , , toroidal coordinates - (V ;V ;V ) velocity of a fluid particle and its projections in toroidal coordinates - g ,g ,g metric tensor components - the Jacobian of transition to curvilinear coordinates - V ₀ velocity at the center of a vortex ring on its symmetry axis - x, y, z Cartesian coordinates - z, y, cylindrical coordinates - a distance from the axis of a torus (V=0) to its axis of symmetry (Oz) - angle between the Oy axis and the line that connects a fluid particle on the streamline =const, which represents a circle [16], with the center of this circle - U _z,U _y velocities in the cylindrical system of coordinates - ₀ stream function of a stationary vortex ring - velocity circulation - U V ₁, velocity of a rectilinear flow at infinity - ₁ stream function of a rectilinear flow - = ₀ + ₁ superposition of two flows - n=k ⁴=V ₁/V ₀ velocity ratio coefficient - R radius of a vortical region - U velocity of fluid particles at the boundary in polar coordinates (r, ) with the center at the coordinate origin (point 0) - fluid density - p ₀,p pressure at infinity and at a certain point of flow - pressure difference Polotsk State University, Polotsk, Belarus. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 68, No. 4, pp. 531–536, July–August, 1995.     

Use of adjoint functions in investigations of heat conduction and transfer processes

An examination is made of the use of adjoint functions in heat conduction and convection theory. Formulas of perturbation theory are obtained for steady and unsteady cases, an interpretation of the physical meaning of adjoint temperature is given, and some applications of the theory are discussed.Notation (r,) thermal conductivity - t(r,) temperature - t *(r,) adjoint temperature - q_V(r,) density of heat release sources - p(r,) a parameter of adjoint equation - r generalized coordinate - time - (r_s, ) heat transfer coefficient - I linear functional of temperature - (r,;r₀,₀) and *(r,; r₀,₀) Green's function for t(r, ) and t *(r, ) - C(r,) volume specific heat - W(r, ) vector distribution of flow velocities - V, S volume and surface areas of body - R radius of HRE - r, radial and angular coordinates - F_in, F_out inlet and outlet flow areas of channel     

Motion of a thermal wave front in a nonlinear medium with absorption

We study the evolution of a thermal perturbation in a nonlinear medium whose thermal conductivity depends on the temperature and the temperature gradient according to a power law.Notation u temperature - k coefficient of thermal conductivity - t time - x spatial variable - x₊ a point on the thermal wave front - a ² generalized coefficient of thermal diffusivity - , , , and s parameters of the process - (x^s) Dirac delta-function - B[, ] a beta function - v(, x), (t) auxiliary functions - A, C, T_o, T_m, T*, R, r, p, and _m constants and parameters Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 39, No. 4, pp. 728–731, October, 1980.     

Interaction between a dislocation and an impurity in KCl: Mg2+ single crystals

The interaction between a dislocation and the impurity in KCl: Mg²⁺ (0.035 mol% in the melt) was investigated at 77–178 K with respect to the two models: one is the Fleischer's model and the other the Fleischer's model taking account of the Friedel relation. The latter is termed the F-F. The dependence of strain-rate sensitivity due to the impurities on temperature for the specimen was appropriate to the Fleischer's model than the F-F. Furthermore, the activation enthalpy, H, for the Fleischer's model appeared to be nearly proportional to the temperature in comparison with the F-F. The Friedel relation between effective stress and average length of the dislocation segments is exact for most weak obstacles to dislocation motion. However, above-mentioned results mean that the Friedel relation is not suitable for the interaction between a dislocation and the impurity in the specimen. Then, the value of H(T _c) at the Fleischer's model was found to be 0.61 eV. H(T _c) corresponds to the activation enthalpy for overcoming of the strain field around the impurity by a dislocation at 0 K. In addition, the Gibbs free energy, G ₀, concerning the dislocation motion was determined to be between 0.42 and 0.48 eV on the basis of the following equation ln / = G ₀/(kT_p0)1 – (T/T _c)^{1/2 –1}(T/T _c)^1/2 + ln ₀/where k is the Boltzmann's constant, T the temperature, T _c the critical temperature at which the effective stress due to the impurities is zero, _p0 the effective shear stress without thermal activation, and ₀ the frequency factor.     

Electrical conductivity,thermoelectric power and dielectric constant of NiWO4

The a.c. electrical conductivity ( _ac), thermoelectric power () and dielectric constant () of antiferromagnetic NiWO₄ are presented. _ac and have been measured in the temperature range 300 to 1000 K and in the temperature range 600 to 1000 K. Conductivity data are interpreted in the light of band theory of solids. The compound obeys the exponential law of conductivity = ₀ exp (–W/kT). Activation energy has been estimated as 0.75eV. The conductivity result is summarized in the following equation =2.86 exp (–0.75 eV/kT)^–1 cm^–1 in the intrinsic region. The material is p-type below 660 K and above 950 K, and is n-type between 660 and 950 K.     

Numerical study of turbulent heat transfer and flow characteristics of hot flow over a sudden-expansion with base mass injection

Summary This study presents the numerical calculations of the fluid flow and turbulent heat transfer characteristics of hot flow over a sudden-expansion with cold air base mass injection. The turbulent governing equations are solved by a control-volume-based finite-difference method with power-law scheme, the well knownk- model, and its associate wall function to describe the turbulent behavior. The velocity and pressure terms of momentum equations are solved by the SIMPLE (Semi-Implicid Method for Pressure-Linked Equation) method. In this study non uniform staggered grids are used. The parameters interested include the inlet Reynolds number (Re), inlet temperature (T₀), and the injection flow rate (Q). The numerical results show that the reattachment lengths are reasonably predicted with a maximum discrepancy within 9.1%. It also shows that the base mass injection suppresses the horizontal velocity and turbulence intensity. In these high temperature heat transfer characteristics, the heat transfer coefficient increased with increasing inlet temperature and inlet Reynolds number, but decreased with increasing injection flow rate of the cooling air.Nomenclature C ₁,C ₂,C turbulent constant - E constant - G generation rate of turbulent kinetic energy - H channel height at inlet - i turbulence intensity - k turbulent kinetic energy - Nu local Nusselt number - q _w heat flux - Re Reynolds number - S source term - T temperature - T ₀ inlet temperature - TI turbulent intensity - U ₀ inlet velocity - U friction velocity - U,V x, y component velocity - Reynolds shear stress - X reattachment length - y ⁺ dimensionless distance from the wall - dependent variables - diffusion coefficient of equation - thermal diffusivity of fluid - density - von Kármán constant - turbulent Prandtl number - dynamic viscosity - kinematic viscosity - _w wall shear stress - turbulent energy dissaption rate - length scale constant     

A model for the thermodynamic properties of metallic rare earth systems with an unstable valence

Rare earth compounds and alloys were regarded for many years as classic examples of ionic magnetism in a solid. Due to the limited spatial extent of the 4f wave functions, the 4f electrons do not significantly participate in the chemical bonding, and hence retain most of their ionic character during the transition from a gas to a solid. However, recent experimental evidence on several metallic rare earth systems indicated an apparent loss of the ionic magnetic moment. Extensive measurements on these systems of the magnetic susceptibility, lattice constant, Mössbauer isomer shift, x-ray photoelectron spectrum, and heat capacity can be qualitatively understood if one postulates (as was first done by Maple and Wohlleben using a model due to Hirst) that the rare earth ion fluctuates between two ionic configurations (valence states) which differ in occupation number by one electron (4f ⁿ , 4f ^n–1 5d ¹). We propose to quantify this simple idea by assigning each valence configuration a finite lifetime atT=0. The two lifetimes _n and _n–1 are converted to bands with widthsh/ _n andh/ _n–1 . The states in each band, however, are forced to retain the ionic properties as determined from Hund's rules. The temperature-dependent contribution of the 4f shell to the magnetic susceptibility, heat capacity, and thermal expansion coefficient is then numerically computed using Fermi-Dirac statistics. By varying _n and _n–1 it is possible to quantitatively describe different types of magnetic behavior: integral valence ( _n =, _n–1 =h/), configuration crossover ( _n = _n–1 ), and Kondo phenomena ( _n > _n–1 ). The results of the model are compared to three well-studied rate earth systems with unstable valence: YbAl ₃ , CePd ₃ , and (LaCe)Al ₂ .Work performed within the research program of the Sonderforschungbereich 125—Aachen/Jülich/Köln