Velocity of second sound and superfluid density near the superfluid transition in3He-4He mixtures

Using superleak condenser transducers, the velocity of second soundU ₂ has been measured near the superfluid transition temperature T in³He-⁴He mixtures with molar concentrationsX of³He of 0.0, 0.038, 0.122, 0.297, and 0.440. We have obtained the superfluid density _s/ fromU ₂ on the basis of linearized two-fluid hydrodynamics. The results for _s/ are consistent with those obtained from the oscillating disk method, as expected from two-fluid hydrodynamics. The value of _s/ at eachX could be expressed by a single power law, _s/=k, where =1-T/R, with the experimental uncertainty. It is found that the exponent is independent of concentration forX0.44 within the experimental uncertainty. This concentration independence of is in agreement with the universality concept. From the conclusion that the values of are universal forX0.44, the concentration dependence of the superfluid component _s is expressed by an empirical equation _s(X, )=²_s(0, ). It is found that corresponds to the volume fraction of⁴He in the superfluid³He-⁴He mixture. The value of is in agreement with that obtained from the measurement of the molar volume by others.This paper is based on a thesis submitted to Tokyo University of Education in partial fulfillment of the requirements for the Ph.D. degree.     

Sound propagation in 3He-4He mixtures near the superfluid transition

Measurements of the acoustic attenuation and dispersion in liquid ³He-⁴He mixtures near the superfluid transition T (x) are reported. The frequency range is /2gp=1–45 MHz and the ³He mole fraction X of the mixtures is 0.007, 0.05, 0.15, and 0.36. Comparisons are made with the measurements of Buchal and Pobell for similar mixtures obtained in the kHz region, and on the whole, the consistency between the two experiments is very satisfactory. An analysis is then performed using both the kHz and MHz data. In the normal phase, where the energy dissipation is caused by order parameter fluctuations having a lifetime _F , the attenuation data can all be scaled according to the expression = (T )f(_F. Here (T )^1+y, with y being a function of the mole fraction X and _F(T–T )^–x, with x increasing weakly with X. In the superfluid phase, we attempt a similar scaling representation, which is found to be fairly successful, but where x(T\s-T ) is roughly 15% larger than x(T>T ). In the superfluid phase we also analyze the attenuation data, assuming the additivity of relaxation and fluctuation-dissipation mechanism, and discuss the relaxation times so derived. In contrast to the attenuation, the dispersion data cannot be brought satisfactorily into a scaling representation. However, at T , we find U()-U(0)^y as predicted by Kawasaki, where y is in good agreement with the values from attenuation experiments.Supported by a grant from the National Science Foundation.     

Universality and the scaling laws in the superfluid phase transition of3He-4He mixtures

From the second-sound velocityU ₂ near the superfluid transition point, the superfluid densities in³He-⁴He mixtures, _s (X) and _s (), were deduced along the paths of constant³He concentrationX and of constant chemical potential difference of³He and⁴He. The following critical exponents of _s are determined: (a) =XX for _s (X) in the(X, T) plane,(b) X for _s (X) in the(, T) plane, and(c) for _s () in the(, T) plane. It is found that and _X change by about 4–6% relative to with increasing³He concentration up toX=0.4 and by 8–10% up toX=0.53. It seems that, belowX=0.53, universality hold for . Values of have been found to be in good agreement with the critical exponent of _s in pure⁴He under constant pressure. The values of and _X forX0.53 are also found to be consistent with the scaling relations in the (,T) plane of³He-⁴He mixture.Work performed in part while at the Electrotechnical Laboratory.     

Shear viscosity of liquid4He and3He-4He mixtures,especially near the superfluid transition

High-resolution measurements of are reported for liquid⁴He and³He-⁴He mixtures at saturated vapor pressures between 1.2 and 4.2 K with particular emphasis on the superfluid transition. Here is the mass density, the shear viscosity, and in the superfluid phase both and are the contributions from the normal component of the fluid ( _n and _n ). The experiments were performed with a torsional oscillator operating at 151 Hz. The mole fraction X of³He in the mixtures ranged from 0.03 to 0.65. New data for the total density and data for _n by various authors led to the calculation of . For⁴He, the results for are compared with published ones, both in the normal and superfluid phases, and also with predictions in the normal phase both over a broad range and close to T. The behavior of and of in mixtures if presented. The sloped/dT near T and its change at the superfluid transition are found to decrease with increasing³He concentration. Measurements at one temperature of versus pressure indicate a decreasing dependence of on molar volume asX(³He) increases. Comparison of at T, the minimum of _n in the superfluid phase and the temperature of this minimum is made with previous measurements. Thermal conductivity measurements in the mixtures, carried out simultaneously with those of , revealed no difference in the recorded superfluid transition, contrary to earlier work. In the appendices, we present data from new measurements of the total density for the same mixtures used in viscosity experiments. Furthermore, we discuss the data for _n determined for⁴He and for³He-⁴He mixtures, and which are used in the analysis of the data.     

Relaxation times and mass diffusion in superfluid dilute3He-4He mixtures

Relaxation times are reported from the transients observed during thermal conductivity _eff and thermal diffusionk _T ^* measurements in superfluid mixtures of³He in⁴He with a layer thickness of 1.81 mm. The experiments extend from 1.7 K toT and over a³He concentration range 10^–6X<5×10^–2. The agreement between the measured and the predicted from the two-fluid thermohydrodynamic equations is satisfactory forX>10^–3 but deteriorates for smaller³He concentrations. This situation is similar to that for _eff andk _T ^* results and indicates that the transport properties in very dilute mixtures with layers of finite thickness are not well understood. ForX>10^–3, the mass diffusion coefficientD _iso for isolated³He in⁴He has been determined from and from _eff measurements. There is an inconsistency by a constant numerical factor between these determinations. This problem might be related to the observations that in the superfluid phase, the relaxation times for different cell heightsh do not scale withh ². FromD _iso derived via the _eff data, the³He impurity-roton scattering cross section is determined. Comparisons with previous work are made.     

Direct and inverse heat-conduction problems for a semiinfinite rod for a partial outflow of heat through the surface

Analytical solutions of the direct and inverse problems of nonstationary heat conduction in a thin semiinfinite rod are given for the case of radiative heat fluxes at the lateral surfaces and a partial outflow of heat by convection and radiation through the end of the rod.Notation thermal diffusivity - x₁ coordinate along the length of the rod - t₁ time - t=t₁/d² dimensionless time (Fourier number) - x=X₁/d relative coordinate - T_o initial temperature - Boltzmann constant - Sk=aT_c ³d/ Stark number - Bi=d/ reduced Biot number - emissivity Translated from Inzhenero-Fizicheskii Zhurnal, Vol. 47, No. 1, pp. 148–153, July, 1984.     

Natural convection in a rectangular enclosure in the presence of a magnetic field with uniform heat flux from the side walls

Summary A numerical study is presented for magnetohydrodynamic free convection of an electrically conducting fluid in a two-dimensional rectangular enclosure in which two side walls are maintained at uniform heat flux condition. The horizontal top and bottom walls are thermally insulated. A finite difference scheme comprising of modified ADI (Alternating Direction Implicit) method and SOR (Successive-Over-Relaxation) method is used to solve the governing equations. Computations are carried out over a wide range of Grashof number, Gr and Hartmann number, Ha for an enclosure of aspect ratio 1 and 2. The influences of these parameters on the flow pattern and the associated heat transfer characteristics are discussed. Numerical results show that with the application of an external magnetic field, the temperature and velocity fields are significantly modified. When the Grashof number is low and Hartmann number is high, the central streamlines are elongated and the isotherms are almost parallel representing a conduction state. For sufficiently large magnetic field strength the convection is suppressed for all values of Gr. The average Nusselt number decreases with an increase of Hartmann number and hence a magnetic field can be used as an effective mechanism to control the convection in an enclosure.List of symbols Ar aspect ratio,H/L - B ₀ induction magnetic field - H ₀ magnetic field,H ₀=B ₀/ _m - g gravitational acceleration - Gr Grashof number,gq(L/k)L ³/v ² - H height of the enclosure - Ha Hartmann number, - k thermal conductivity - Nu local Nusselt number - average Nusselt number - p pressure - Pr Prandtl number, / - q heat flux - t time - T dimensionless temperature, (–₀)/q(L/k) - u vertical velocity - U dimensionless vertical velocity,uL/ - v horizontal velocity - V dimensionless horizontal velocity,vL/ - x vertical coordinate - X dimensionless vertical coordinate,x/L - y horizontal coordinate - Y dimensionless horizontal coordinate,y/L - thermal diffusivity - thermal expansion coefficient - temperature - ₀ reference temperature - density - kinematic viscosity - viscosity - _m magnetic permeability - electrical conductivity - stream function - dimensionless stream function, / - dimensionless time,t/L ² - vorticity - dimensionless vorticity, L ²/ - X grid spacing inX-direction - Y grid spacing inY-direction - time increment - ² Laplacian operator     

Thermal conductivity of methane

This paper reports thermal conductivity data for methane measured in the temperature range 120–400 K and pressure range 25–700 bar with a maximum uncertainty of ± 1%. A simple correlation of these data accurate to within about 3% is obtained and used to prepare a table of recommended values.Nomenclature a _k ,b _ij ,b _k Parameters of the regression model, k= 0 to n; i =0 to m; j =0 to n - P Pressure (MPa or bar) - Q _kl Heat flux per unit length (mW · m^–1) - t time (s) - T Temperature (K) - T _cr Critical temperature (K) - T _r reduced temperature (= T/T _cr) - T _w Temperature rise of wire between times t ₁ and t ₂ (deg K) - T ^* Reduced temperature difference (TT _cr)/T _cr - Thermal conductivity (mW · m^–1 · K^–1) - ₁ Thermal conductivity at 1 bar (mW · m^–1 · K^–1) - _bg Background thermal conductivity (mW · m^–1 · K^–1) - _cr Anomalous thermal conductivity (mW · m^–1 · K^–1) - _e Excess thermal conductivity (mW · m^–1 · K^–1) - Density (g · cm^–3) - _cr Critical density (g · cm^–3) - _r Reduced density (= / _cr) - ^* Reduced density difference (– _cr )/ _cr      

Experimental study of the boundary between single-phase convection and boiling in underheated cryogenic liquids

The effect of pressure and underheating on the position of the boundary between heat-transfer regimes in liquid helium and hydrogen is investigated.Notation q heat flux - p pressure - =T_s–T underheating - T_s saturation temperature - T temperature of liquid - T=T_wa – T T_s=T_wa – T_s - T_wa temperature of heat-emitting surface - A,a, B, b, C constants - m, n indices - Nu Nusselt number - Ra Rayleigh number - thermal conductivity - coefficient of cubical expansion - kinematic viscosity - g acceleration - standard deviation Indices 01 conditions of convection-boiling transition - 02 conditions of boiling-convection transition Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 42, No. 1, pp. 5–11, January, 1982.     

The behavior of the4He viscosity near the superfluid transition

Measurements of the shear viscosity at saturated vapor pressure through the lambda transition indicate a singular behavior of the form |1 – (/)|=A ^x , (where =|1–(T/T )|, with equal values for the critical exponent on both sides of the transition.Work sponsored by Consiglio Nazionale delle Ricerche, Rome (Italy).     

Surface tension of helium at the superfluid critical end-point

The singularity in the vapor-liquid interfacial tension, (T), of helium at the transition to superfluidity is analyzed theoretically. The universal amplitude ratio R ⁺ =K₊( ₀ ⁺ )^d–1/k _B T , where K₊ and K_– are the amplitudes of the |T–T|^µ singularity in , with =1.34 ₃ , is known from recent work to first order in =4–d for the general n-vector model in d dimensions. Extrapolation to d=3 for n=2 indicates R ⁺ =0.05–0.08, which is shown to be consistent with the experimental data. Further analysis of the experiments establishes that the universal ratio Q=K₊/K _– exceeds 0.35, and is consistent with the recent prediction Q0.9; this demonstrates the inadequacies of earlier theoretical treatments. The existence in the observed surface tension of an anomalous, negative contribution of unknown origin at a few millikelvin beneath T is stressed.     

Nonlinear asymptotic theory of hypersonic flow past a circular cone

Summary The hypersonic small-disturbance theory is reexamined in this study. A systematic and rigorous approach is proposed to obtain the nonlinear asymptotic equation from the Taylor-Maccoll equation for hypersonic flow past a circular cone. Using this approach, consideration is made of a general asymptotic expansion of the unified supersonic-hypersonic similarity parameter together with the stretched coordinate. Moreover, the successive approximate solutions of the nonlinear hypersonic smalldisturbance equation are solved by iteration. Both of these approximations provide a closed-form solution, which is suitable for the analysis of various related flow problems. Besides the velocity components, the shock location and other thermodynamic properties are presented. Comparisons are also made of the zeroth-order with first-order approximations for shock location and pressure coefficient on the cone surface, respectively. The latter (including the nonlinear effects) demonstrates better correlation with exact solution than the zeroth-order approximation. This approach offers further insight into the fundamental features of hypersonic small-disturbance theory.Notation a speed of sound - H unified supersonic-hypersonic similarity parameter, - K hypersonic similarity parameter, M - M freestream Mach number - P pressure - T temperature - S entropy - u, v radial, polar velocities - V freestream velocity - shock angle - cone angle - density - density ratio, /() - ratio of specific heats - polar angle - stretched polar angle, / - (), (), () gage functions     

Tests on a flow cryostat with series cooling

Measurements and calculations on a flow cryostat with serial cooling have given equivalent thermal schemes that have been tested for adequacy and consequent simple working formulas.Notation T_c, T_i, T_w, T_f temperatures of case, body i, tube wall, and flowing coolant in K - T₀ and T_e coolant temperatures at inlet and exit for heat exchanger and pipes in K - T_wi mean pipe wall temperature at points of attachment of bundles from body i in K - T_wn pipe wall temperature at point of attachment for bundle n in K - (_i)_n and _i thermal conductivities of bundle n and all bundles from body i in W/K - _ij thermal conductivity between bodies i and j in W/K - _ci, , _cw thermal conductivities from case to body i and total and radiative conductivities from case to pipe in W/K - _c convective heat-transfer coefficient between pipe and coolant in W/m²·K - _r radiative heat-transfer coefficient between case and pipe in W/m²·K - pipe material thermal conductivity in W/m·K - c specific heat of helium at constant pressure in J/kg·K - q and q_r correspondingly densities of the total heat flux and radiative flux to the pipe in W/m² - P_r heat flux along bundle r in W - M coolant mass flow rate in kg/sec - F tube cross section area in m² - S_i and S_o inside and outside surface areas of pipe in m² - L pipe length in m - ¯x=x/L relative coordinate along pipe axis - ¯x_r relative coordinate for bundle r attachment - R total number of bundles - N_i number of bundles cooling body i - J_i number of bodies linked by heat bridges to body i - _i relative error in calculating the temperature of body i by comparison with numerical result in % - _w mean relative error in heat exchanger temperature calculated numerically by comparison with temperature from (4) taken at ten equally separated points in % - (¯x-¯x_r) Dirac function Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 56, No. 5, pp. 760–767, May, 1989.     

Shear viscosity of the B phase of superfluid3He

The shear viscosity (T) in the Balian-Werthamer (BW) state of superfluid ³ He is calculated variationally throughout the region 0t 1(t=T/T _c) from the transport equation for Bogoliubov quasiparticles. Coherence factors are treated exactly in the calculation of the collision integral. The numerical result for =_s= _s(T)/_n(T_c) agree very well with experiment in the range 0.8t1.0. Analytic expressions = 0.577 (1–1.0008t) and =1–(23/64) [=(T)/k _B T] are obtained in the low-temperature region and in the vicinity ofT _c, respectively. From the numerical analysis it is shown that the latter equation is valid only in the temperature range 0.9997t1.0.Supported by the Research Institute for Fundamental Physics, Kyoto University.     

On the Bose-Einstein condensation of nonzero-momentum cooper pairs. A second-order phase transition

Based on the BCS Hamiltonian, the normal-to-super phase transition is investigated, approaching the critical temperatureT _c from the high-temperature side. Nonzero-momentum Cooper pairs, that is, pairs of electrons (holes) with antiparallel spins and nearly opposite momenta aboveT _c in the bulk limit, are shown to move like independent bosons with the energy vs. momentump relation=1/2v_F , where _F represents the Fermi velocity (1/2m* _F ² _FFermi energy). The system of free Cooper pairs undergoes a phase transition of the second order with the critical temperatureT _c given byk _B T _c=1/2(²³ _F ³ n/1.20257)^1/3 wheren is the number density of Cooper pairs. The ratio of the jump of the heat capacity, C, to the maximum heat capacity,C _s, is a universal constant: C/C _s=0.60874; this number is close to the universal constant 0.588 obtained by the finite-temperature BCS theory. The physical significance of these results is discussed, referring to the well-known BCS theory, which treats the many-Cooper-pair ground state exactly and the thermodynamic state belowT _c approximately. An explanation is proposed on the question why sodium should remain normal down to 0 K, based on the band structures with the hypothesis that the supercondensate composed of zero-momentum electron and hole Cooper pairs is electrically neutral.     

Energy loss distribution of charged particles in thin absorbers

Conclusions 1. For thin absorbers, the Landau theory incorporating certain above-listed corrections is a general-purpose one and is in good agreement with experiment both for heavy and light charged particles 0.01. 2. For intermediate layers (0.01 1 and 1), exact solutions are provided by the Vavilov theory [11]. 3. The numerical-analytic method of plotting the energy loss distribution function proposed in [13] is suitable for very thin ( 0.01) absorbers.Translated from Izmeritel'naya Tekhnika, No. 3, pp. 60–62, March, 1970.     

Vibrational properties of wood along the grain

The dynamic Young's modulus (E_L) and loss tangent (tan _L ) along the grain, dynamic shear modulus (G_L) and loss tangent (tan _S) in the vertical section, and density () of a hundred spruce wood specimens used for the soundboards of musical instruments were determined. The relative acoustic conversion efficiency ( ) and a ratio reflecting the anisotropy of wood (, (E_L/G_L)(tan _S/tan _L)) were defined in order to evaluate the acoustic quality of wood along the grain. There was a positive correlation between and , and the variation in was larger than that in . It seemed logical to evaluate the acoustic quality of spruce wood by a measure of . By using a cell wall model, those acoustic factors were expressed with the physical properties of the cell wall constituents. This model predicted that the essential requirement for an excellent soundboard is smaller fibril angle of the cell wall, which yields higher and higher . On the other hand, the effects of chemical treatments on the and of wood were clarified experimentally and analyzed theoretically. It was suggested that the and of wood cannot be improved at the same time by chemical treatment.     

Thermal Expansion and Isothermal Compressibility of TlIn1 – xNdxSe2 (0 ≤ x≤ 0.08) Crystals

The thermal expansion coefficient () and isothermal compressibility ( _T ) of TlIn_{1 – x} Nd _x Se₂(0 x 0.08) crystals were measured between 77 and 400 K. In the range 77–160 K, both and _T increase with temperature, the increase in being much steeper. At higher temperatures, and _T change very little. The observed composition dependences of and _T are interpreted in terms of energy-band structure.     

Solution of thermal conductivity problems with boundary condition of the third kind and arbitrary coordinate and time dependence of the biot number

An analytical solution of the thermal conductivity problem with boundary conditions of the third kind and arbitrary coordinate and time dependence of the Biot number is found in the form of a converging series of quadratures.Notation , z dimensionless coordinates - dimensionless temperature - Q dimensionless volume heat-liberation density per unit time - Fo=/² Fourier number - Bi₁(, Fo)=(, Fo) · / Biot number - thermal diffusivity coefficient - plate thickness - time - (, Fo) heat-liberation coefficient - thermal conductivity coefficient - i summation index - Jo zero order Bessel function of the first kind Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 41, No. 3, pp. 536–540, September, 1981.     

Nonlinear axisymmetric response of orthotropic thin truncated conical and spherical caps

Summary An approximate analytical solution of the large deflection axisymmetric response of polar orthotropic thin truncated conical and spherical shallow caps is presented. Donnell type equations are employed. The deflection is approximated by a one term mode shape satisfying the boundary conditions. The Galerkin's method is used to get the governing equation for the deflection at the hole. Nonlinear free vibration response and the response under uniformly distributed static and step function loads are obtained. The effect of various parameters is investigated.Notations A, A ^* Inward and outward amplitudes - a, b, h Base radius, inner radius and thickness of the cap - D M h ³/[12(–v ² )] - E ,E Young's moduli - H ^*,H Apex height, dimensionless apex heght:H ^*/h - N _, Stress resultants - p ^1/2 - q Uniformly distributed load - Q,Q₀ Dimensionless load: , dimensionless step load - Q, Q ₀ Dimensionless load: , step load - t, Time, dimensionless time: t - T _A Ratio of nonlinear periodT for inward amplitudeA and the linear periodT _L - w ^* Normal displacement at middle surface - w Dimensionless displacement:w ^*/h - ₁ Linear parameter of static response - Orthotropic Parameter:E /E - Mass density - ₂,₃ Quadratic and cubic nonlinearity parameters - b/a - v ,v Poisson's ratios - Dimensionless radius:r/a - ^*, Stress function, dimensionless stress function: - ₀ ^* ,₀ Linnear frequency, dimensionless frequency: With 7 Figures