The effect of thermal prehistory on turbulent separated flow at sudden tube expansion

Results of numerical investigation of the effect of heat boundary layer thickness in front of a sudden expansion of a round tube on turbulent transfer in the zone of flow separation, attachment, and relaxation are presented. Before separation the flow was hydrodynamically stable, and the heat layer in front of expansion could change its thickness in maximally possible limits: from zero to a half of tube diameter. The Reynolds number varied from 6.7·10³ to 1.33·10⁵. It was found that the growth of heat layer thickness leads to reduction of heat transfer intensity in the separation area and moving away of the coordinate of maximal heat transfer from the place of tube expansion. Generalizing dependence for the maximum Nusselt number is given for variation of the heat layer thickness. Comparison with experimental data of [1] proved the main behavior tendencies of heat and mass transfer processes in separation flows behind a backward-facing step with different thermal prehistory.     

Scaling of heat transfer in gas-gas injector combustor

The scaling of heat transfer in gas-gas injector combustor is investigated theoretically, numerically and experimentally based on the previous study on the scaling of gas-gas combustion flowfield. The similarity condition of the gas-gas injector combustor heat transfer is obtained by conducting a formulation analysis of the boundary layer Navier-Stokes equations and a dimensional analysis of the corresponding heat transfer phenomenon. Then, a practicable engineering scaling criterion of the gas-gas injector combustor heat transfer is put forward. The criterion implies that when the similarity conditions of inner flowfield are satisfied, the size and the pressure of gas-gas combustion chamber can be changed, while the heat transfer can still be qualitatively similar to the distribution trend and quantitatively correlates well with the size and pressure as q ∝ pc0 .8d t-0.2. Based on the criterion, single-element injector chambers with different geometric sizes and at different chamber pressures ranging from 1 MPa to 20 MPa are numerically simulated. A single-element injector chamber is designed and hot-fire tested at seven chamber pressures from 0.92 MPa to 6.1 MPa. The inner wall heat flux are obtained and analysed. The numerical and experimental results both verified the scaling criterion in gas-gas injector combustion chambers under different chamber pressures and geometries.     

Experimental study and mathematical modelling of heat transfer processes in heat accumulating media

Experimental results on reversing non-stationary heat transfer are presented for filtration of an air flow through an immobile heat accumulating medium consisting of lead (D = 2.0, 3.5, and 4.5 mm) and glass (D = 3.2 mm) balls. The studied device imitated the cyclic modes of heat regeneration in the ventilation system for domestic and office rooms. Dependency between the time of flow switching and Re number was measured. The mathematical model describing heat transfer between a gas flow and an immobile layer of balls was developed. Good correspondence between the experimental data and calculation results is observed for high Reynolds numbers. For low Re numbers the effect of heat losses is considerable, and experimental time of flow switching is shorter than the calculation one. The work was financially supported by the President of the Russian Federation (Grant No. NSh 6526.2006.3), Russian Foundation for Basic Research (Grant No. 06-08-00982), Foundation “Global energy” and Program “Energy saving of SB RAS”.     

Heat transfer and crisis phenomena with intense boiling in the falling wave liquid films

Experimental data on heat transfer with intense evaporation in the falling films of liquid nitrogen were analysed. According to data generalization, heat transfer at evaporation becomes more intense under the precrisis modes at high heat fluxes for two studied boundary conditions on the heat-releasing surface: T _w ≈ const and q _w ≈ const. The relative contributions of conductive and convective components of heat transfer for different heat fluxes were estimated due to statistical treatment of the wave characteristics carried out by the capacitance probes for measurement of the local liquid film thickness. It was found out that heat transfer intensification is mainly caused by a drastic decrease in thermal resistance of the local zones with intensely evaporating residual layer between large waves. At that, the convective component of heat transfer related to wave perturbations on a free surface of a liquid film decreases significantly with a rise of heat fluxes. New data on pulsations of the local temperature of the heat-releasing surface were obtained at different points along the flow with the modes of “dry spot” formation. The work was financially supported by the Russian Foundation for Basic Research (grant No. 03-02-04027-NNIO-a) and DFG (Deutsche Forschungsgemeinschaft, project Nos. Re-463-37-1 and SFB-540).     

Experimental study of the evolution of periodic controlled disturbances in a hypersonic viscous shock layer on a flat plate

Characteristics of the fields of mean density and density fluctuations measured with introduction of periodic disturbances into a hypersonic viscous boundary layer on a flat plate are presented. The experiments are performed for a flow Mach number M_∞ = 21, Reynolds number per meter Re_1∞ = 6·10⁵ m⁻¹, and temperature factor of the surface T _w /T ₀ = 0.26. The disturbances are introduced into the shock layer by an oblique gasdynamic whistle. The work was financially supported by the Russian Foundation for Basic Research (Grants Nos. 04-01-00474 and 05-08-33436).     

An Experimental Investigation of Liquid Jet Impingement and Single-Phase Spray Cooling Using Polyalphaolefin

Experiments on triangular and rectangular array jet impingement and single-phase spray cooling have been performed to determine the effect of both cooling techniques on heat transfer coefficient (h) and the coolant mass flux required for a given cooling load. Experiments were performed with circular orifices and nozzles for different H/D values from 1.5 to 26 and Reynolds number range of 219 to 837, which is quite lower than the ranges employed in widely used correlations. The coolant used was polyalphaolefin. The experiments simulated the boundary condition produced at the surface of the stator of a high power low-density generator or motor. For the custom fabricated orifices, commercial nozzles, and conditions used in this study, both cooling configurations showed enhancement of heat transfer coefficient as H/D increases to a certain limit after which it starts to decrease. The heat transfer coefficient always increases with Reynolds number. In keeping with previous studies, single-phase spray cooling technique can provide the same heat transfer coefficient as jets at a slightly lower mass flux, but with much higher-pressure head. Special Nu_d correlations that account for the range of parameters and coolant studied in this work are derived.     

HEAT AND MASS TRANSFER AND WALL SHEAR STRESS IN VERTICAL GAS-LIQUID FLOW

Results of an experimental investigation of heat and mass transfer and wall shear stress at gas-liquid flow in a vertical tube are presented. Local wall shear stress and mass transfer coefficients were measured by an electrochemical method. Experiments were performed in the range of Reynolds number variation with respect to liquid Rc_i, = 8.5 × 10³-5.4 × 10⁴, gas Re_g = 3 × 10³-1.4 × 10⁵, pressure 0.1-1 MPa. The relationship between heat and mass transfer and wall shear at gas-liquid flows is shown to exist. The results of measuring heat and mass transfer coefficients are generalized by formulas applied to calculate heat and mass transfer in single-phase turbulent flow.     

Infrared-induced temperature distributions of solid D2 ices

A specific-wavelength infrared (IR) light (λ=3140 nm) was irradiated into a solid D₂ ice prepared in a cylinder target cell. The temperature in the solid D₂ ice oscillated periodically with a high amplitude when irradiated by the IR light. The temperature oscillation has been well explained based on the two-dimensional heat transfer theory plus the IR-irradiation effect. The transmission optical imaging reveals that such a temperature oscillation is favorable to recrystallize the solid D₂ ice from multicrystal to quasi single crystal. This suggests an efficient method to layer the solid hydrogen-isotope ice for the inertial-confinement-fusion (ICF) experiments.     

Implementation Aspects of 3D Lattice-BGK: Boundaries, Accuracy, and a New Fast Relaxation Method

In many realistic fluid-dynamical simulations the specification of the boundary conditions, the error sources, and the number of time steps to reach a steady state are important practical considerations. In this paper we study these issues in the case of the lattice-BGK model. The objective is to present a comprehensive overview of some pitfalls and shortcomings of the lattice-BGK method and to introduce some new ideas useful in practical simulations. We begin with an evaluation of the widely used bounce-back boundary condition in staircase geometries by simulating flow in an inclined tube. It is shown that the bounce-back scheme is first-order accurate in space when the location of the non-slip wall is assumed to be at the boundary nodes. Moreover, for a specific inclination angle of 45 degrees, the scheme is found to be second-order accurate when the location of the non-slip velocity is fitted halfway between the last fluid nodes and the first solid nodes. The error as a function of the relaxation parameter is in that case qualitatively similar to that of flat walls. Next, a comparison of simulations of fluid flow by means of pressure boundaries and by means of body force is presented. A good agreement between these two boundary conditions has been found in the creeping-flow regime. For higher Reynolds numbers differences have been found that are probably caused by problems associated with the pressure boundaries. Furthermore, two widely used 3D models, namelyD₃Q₁₅andD₃Q₁₉, are analysed. It is shown that theD₃Q₁₅model may induce artificial checkerboard invariants due to the connectivity of the lattice. Finally, a new iterative method, which significantly reduces the saturation time, is presented and validated on different benchmark problems.     

Thermodiffusion of interstitials in a film

A general expression is obtained for the thermodiffusion coefficientD _T of a film, also for its dependence on temperature, film thickness and type of concentration distribution of the interstitials, according to its profile. In neglecting interaction of interstitials with each other in a model of a film with single-plane boundaries with the vacuum and the substrate, it is shown thatD _T>D _T ^(∞) is possible, whereD _T ^(∞) is the thermodiffusion coefficient of a bulk specimen. Ural State Technical University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 9, pp. 43–46, September, 1996.     

Numerical simulation of turbulent natural convection of oxide heat generating melt in a core catcher at NPP with VVER-1000

The problem statement and simulation results are presented concerning turbulent natural convection in a vertical cylindrical molten pool with internal heat generation and other parameters (inner Rayleigh number Ra _i ∼ 10¹⁶–10¹⁷) corresponding to oxide core melt in a core catcher for NPP with VVER-1000. Commercial code FLUENT 6.3 was used for CFD calculations. The results on heat transfer are approximated by power law correlations for mean Nusselt numbers vs. Rayleigh number and pool height, describing the heat transfer at upper, lateral, and total boundaries of the cylinder. The influence of volumetric heat generation and material properties is studied. Spatial distribution of wall heat transfer is analyzed for different pool heights possible in the real core catcher. Along with serial calculations with isothermal boundary conditions, the cases with heat radiation conditions are considered. The results may be used for estimations of heat transfer and melt overheating in a VVER core catcher and for coefficient identification of simplified models of integrated system severe-accident codes.     

The effects of GaN capping layer thickness on two-dimensional electron mobility in GaN/AlGaN/GaN heterostructures

In this paper we present a study of the effect of GaN capping layer thickness on the two-dimensional (2D)-electron mobility and the two-dimensional electron gas (2DEG) sheet density which is formed near the AlGaN barrier/buffer GaN layer. This study is undertaken using a fully numerical calculation for GaN/Al_xGa_1−xN/GaN heterostructures with different Al mole fraction in the Al_xGa_1−xN barrier, and for various values of barrier layer thickness. The results of our analysis clearly indicate that increasing the GaN capping layer thickness leads to a decrease in the 2DEG density. Furthermore, it is found that the room-temperature 2D-electron mobility reaches a maximum value of approximately 1.8×10³ cm² /Vs⁻¹ for GaN capping layer thickness grater than 100 Å with an Al_0.32Ga_0.68N barrier layer of 200 Å thick. In contrast, for same structure, the 2DEG density decreases monotonically with GaN capping layer thickness, and eventually saturates at approximately 6×10¹² cm⁻² for capping layer thickness greater than 500 Å. A comparison between our calculated results with published experimental data is shown to be in good agreement for GaN capping layers up to 500 Å thickness.     

Non-stationary conjugate heat exchange and phase transitions at the high-energy surface processing. Part 2. Simulation of the technological processes

The developed general physical-mathematical model, FEM-based calculation procedure as well as the software were in practical use to simulate the processes of the non-stationary conjugate heat exchange and phase transformations during the processing of the surface with a high-concentrated energy fluxes, with a stationary, pulsed, and movable heating sources (the processing, including the surface fusing with a quasilaminar plasma jet, transfer electric arc and impulse electron beam; cleaning of the metal substrate surfaces from an oxide layer with the aid of a cathode vacuum arc, etc). The processes of practical importance with considerably different spatial and temporal scales featuring the density of the heat fluxes power q ∈ [10⁷; 10¹⁴] W/m² have been studied. The work was financially supported by the Russian Foundation for Basic Research (Grant No. 07-08-00209).     

Sub-Doppler fluorescence on the atomic <Emphasis Type="Italic">D</Emphasis><Subscript>2</Subscript> line of a submicron rubidium-vapor layer

An extremely thin cell (ETC) with the thickness of a Rb atomic vapor layer in the range of 100–300 nm was fabricated. It is demonstrated that a simple laser-diode technique with a single resonant light beam is sufficient to observe separately all of the atomic hyperfine transitions of the D ₂ line of Rb (780 nm) and also allows us to measure the relative transition probabilities of the hyperfine transitions. The onset of collisional self-broadening of the hyperfine transitions as the number density of atoms increases was studied. The detrimental role of the atoms with slow longitudinal velocity in the sub-Doppler response of the Rb ETC is demonstrated by studies in which the cell is tilted from normal incidence of the laser beam. It is also shown that using an ETC allows us to resolve in a moderate external magnetic field the Zeeman splitting of the hyperfine transitions of the ⁸⁷Rb D ₁ transition F _g=1F _e=1,2. Received: 19 February 2003 / Revised version: 4 April 2003 / Published online: 2 June 2003 RID="*" ID="*"Corresponding author. Fax: +374/32-31172, E-mail: david@ipr.sci.am     

Filler-matrix thermal boundary resistance of diamond-copper composite with high thermal conductivity

A composite material with a high thermal conductivity is obtained by capillary infiltration of copper into a bed of diamond particles of 400 μm size, the particles having been pre-coated with tungsten. The measured thermal conductivity of the composite decreases from 910 to 480 W m⁻¹ K⁻¹ when the coating thickness is increased from 110 to 470 nm. Calculations of the filler/matrix thermal boundary resistance R and the thermal conductivity of the coating layer λ _i using differential effective medium, Lichtenecker’s and Hashin’s models give similar numerical values of R and λ _i ≈ 1.5 W m⁻¹ K⁻¹. The minimal thickness of the coating h ∼ 100 nm necessary for ensuring production of a composite while maximizing its thermal conductivity, is of the same order as the free path of the heat carriers in diamond (phonons) and in copper (electrons). The heat conductance of the diamond/tungsten carbide coating/copper interface when h is of this thickness is estimated as (0.8–1) × 10⁸ W m⁻² K⁻¹ and is at the upper level of values characteristic for perfect dielectric/metal boundaries.     

Structural phase changes in a titanium-silicon system modified by high-current electron beams and compression plasma flows

Structural phase changes in a titanium-silicon system treated by low-energy high-current electron beams (HCEBs) and compression plasma flows (CPFs) with the duration 100 μs and the energy density 12–15 J/cm² are studied. Scanning electron microscopy, X-ray diffraction and electron microprobe analysis are used in this work. The formation of a titanium-doped silicon layer 10–25 μm thick, titanium silicides (TiSi₂ under HCEBs and Ti₅Si₃ under CPF treatment), silicon dendrites, and needle-like eutectics (typical size of precipitates is about 50 nm) is revealed. It is shown via the results of numerical simulation that the thickness of the metal-doped layer is mainly controlled by the power density value and the surface nonuniformity of the heat flow over the target surface. The thermodynamic regularities of phase formation are discussed, taking into account heat transfer between the silicide nuclei and solid silicon.     

Carmeli's Cosmology Fits Data for an Accelerating and Decelerating Universe Without Dark Matter or Dark Energy

A new relation for the density parameter Ω is derived as a function of expansion velocity υ based on Carmeli's cosmology. This density function is used in the luminosity distance relation D _L. A heretofore neglected source luminosity correction factor (1 − (υ/c)²)^−1/2 is now included in D _L. These relations are used to fit type Ia supernovae (SNe Ia) data, giving consistent, well-behaved fits over a broad range of redshift 0.1 < z < 2. The best fit to the data for the local density parameter is Ω_m = 0.0401 ± 0.0199. Because Ω_m is within the baryonic budget there is no need for any dark matter to account for the SNe Ia redshift luminosity data. From this local density it is determined that the redshift where the universe expansion transitions from deceleration to acceleration is z _t = 1.095^+0.264 _−0.155. Because the fitted data covers the range of the predicted transition redshift z _t, there is no need for any dark energy to account for the expansion rate transition. We conclude that the expansion is now accelerating and that the transition from a closed to an open universe occurred about 8.54 Gyr ago.     

Capacitance-voltage behaviour inp-type ɛ-GaSe single crystal ∥c at different temperatures

Summary  Capacitance-voltage measurements have been carried out onp-type ɛ-GaSe single crystal ∥c in the temperature range 300 to 360 K, with applied voltages of -1, 0 and +1 V. TheC-V measurements in this temperature range have shown a shift in capacitanceC and conductanceG to the higher values with an increase in temperature. The depletion layer widthW, the Debye length L_D and the doping densityN _α have been worked out and plots ofN _α vs. W have shown a decrease inW with an increase in temperature. The plots of L_D vs. N _α vary as 1/N_α^1/2, which gives N_αL_D ⋍ 3.3 × 10¹¹ charges/m² for doping density of 10¹⁶m⁻³. The values ofG at different temperatures have been used to obtain the activation energies, which are found to be ΔE ⋍ 0.11 eV for -1 and +1 V applied voltages, and ΔE ⋍ 0.06 eV for zero volt. The authors of this paper have agreed to not receive the proofs for correction.     

Influence of europium doping on the sensitization of luminescence in GaAs/AlGaAs and InGaN/GaN quantum-well structures

A correlation between the photoluminescence spectra and structural parameters of Eu-doped quantum- well nanostructures InGaN/GaN and GaAs/AlGaAs is established. It is shown that the incorporation of rare-earth ions initiates lattice (as a rule, compressive) strains. The excitation migration in structures of high perfection stimulates transfer of nonequilibrium carriers to the ⁵ D ₂-⁵ D ₀ atomic levels of the Eu ion. In less perfect structures, the insertion of a rare-earth ion leads to the formation of isovalent traps in GaN layers capable of effectively capturing nonequilibrium carriers, which increases the intensity of photoluminescence of the structure by one order of magnitude.     

Pressure-induced phase transition in silicon nitride material

The equilibrium crystal structures,lattice parameters,elastic constants,and elastic moduli of the polymorphs α-,β-,and γ-Si3N4,have been calculated by first-principles method.β-Si3N4 is ductile in nature and has an ionic bonding.γSi3N4 is found to be a brittle material and has covalent chemical bonds,especially at high pressures.The phase boundary of the β→γ transition is obtained and a positive slope is found.This indicates that at higher temperatures it requires higher pressures to synthesize γ-Si3N4.On the other hand,the α→γ phase boundary can be described as P = 14.37198+ 3.27 × 10?3T-7.83911 × 10?7T2-3.13552 × 10?10T3.The phase transition from α-to γ-Si3N4 occurs at 16.1 GPa and 1700 K.Then,the dependencies of bulk modulus,heat capacity,and thermal expansion on the pressure P are obtained in the ranges of 0 GPa-30 GPa and 0 K-2000 K.Significant features in these properties are observed at high temperatures.It turns out that the thermal expansion of γ-Si3N4 is larger than that of α-Si3N4 over wide pressure and temperature ranges.The evolutions of the heat capacity with temperature for the Si3N4 polymorphs are close to each other,which are important for possible applications of Si3N4.