Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Experimental investigation on TBAB clathrate hydrate slurry flows in a horizontal tube: Forced convective heat transfer behaviors

Tetra-n-butyl-ammonium bromide (TBAB) clathrate hydrate slurry (CHS) is one kind of secondary refrigerants, which is promising to be applied into air-conditioning or latent-heat transportation systems as a thermal storage or cold carrying medium for energy saving. It is a solid-liquid two phase mixture which is easy to produce and has high latent heat and good fluidity. In this paper, the heat transfer characteristics of TBAB slurry were investigated in a horizontal stainless steel tube under different solid mass fractions and flow velocities with constant heat flux. One velocity region of weakened heat transfer was found. Moreover, TBAB CHS was treated as a kind of Bingham fluids, and the influences of the solid particles, flow velocity and types of flow on the forced convective heat transfer coefficients of TBAB CHS were investigated. At last, criterial correlations of Nusselt number for laminar and turbulent flows in the form of power function were summarized, and the error with experimental results was within 20%.

Fully developed laminar flow and heat transfer in smooth trapezoidal microchannel

Numerical analysis of fully developed laminar slip flow and heat transfer in trapezoidal micro-channels has been studied with uniform wall heat flux boundary conditions. Through coordinate transformation, the governing equations are transformed from physical plane to computational domain, and the resulting equations are solved by a finite-difference scheme. The influences of velocity slip and temperature jump on friction coefficient and Nusselt number are investigated in detail. The calculation also shows that the aspect ratio and base angle have significant effect on flow and heat transfer in trapezoidal micro-channel. (c) 2005 Elsevier Ltd. All rights reserved.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

Generalization of Response Number for Dynamic Plastic Response of Shells Subjected to impulsive Loading

A dimensionless number, termed as response number in Zhao [Archive of Applied Mechanics 68 (1998) 524], has been suggested for the dynamic plastic response of beams and plates made up of rigidly perfect plastic materials subjected to dynamic loading. Many theoretical and experimental results can be reformulated into new concise forms with the response number. The concept of a new dimensionless number, response number, termed as Rn(n), is generalized in Zhao [Forschung im Ingenieurwesen 65 (1999) 107] to study the elastic, plastic, dynamic elastic as well as dynamic plastic buckling problems of columns, plates as well as shells. The response number Rn(n) is generalized to the dynamic behaviour of shells of various shapes in the present paper.

Application of response number for dynamic plastic response of plates subjected to impulsive loading

A dimensionless number, termed response number, is applied to the dynamic plastic response of plates subjected to dynamic loading. Many theoretical and experimental results presented by different researchers are reformulated into new concise forms with the response number. The advantage of the new forms is twofold: (1) they are more physically meaningful, and (2) they are independent of the choice of units, thus, they have wider range of applications.

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Suggestion of a new dimensionless number for dynamic plastic response of beams and plates

A dimensionless number, termed response number in the present paper, is suggested for the dynamic plastic response of beams and plates made of rigid-perfectly plastic materials subjected to dynamic loading. It is obtained at dimensional reduction of the basic governing equations of beams and plates. The number is defined as the product of the Johnson's damage number and the square of the half of the slenderness ratio for a beam; the product of the damage number and the square of the half of the aspect ratio for a plate or membrane loaded dynamically. Response number can also be considered as the ratio of the inertia force at the impulsive loading to the plastic limit load of the structure. Three aspects are reflected in this dimensionless number: the inertia of the applied dynamic loading, the resistance ability of the material to the deformation caused by the loading and the geometrical influence of the structure on the dynamic response. For an impulsively loaded beam or plate, the final dimensionless deflection is solely dependent upon the response number. When the secondary effects of finite deflections, strain-rate sensitivity or transverse shear are taken into account, the response number is as useful as in the case of simple bending theory. Finally, the number is not only suitable to idealized dynamic loads but also applicable to dynamic loads of general shape.

Towards an Understanding of the Influence of Sedimentation on Colloidal Aggregation by Peclet Number

The Peclet number is a useful index to estimate the importance of sedimentation as compared to the Brownian motion. However, how to choose the characteristic length scale for the Peclet number evaluation is rather critical because the diffusion length increases as the square root of the time whereas the drifting length is linearly related to time. Our Brownian dynamics simulation shows that the degree of sedimentation influence on the coagulation decreases when the dispersion volume fraction increases. Therefore using a fixed length, such as the diameter of particle, as the characteristic length scale for Peclet number evaluation is not a good choice when dealing with the influence of sedimentation on coagulation. The simulations demonstrated that environmental factors in the coagulation process, such as dispersion volume fraction and size distribution, should be taken into account for more reasonable evaluation of the sedimentation influence.

Instability from steady and axisymmetric to steady and asymmetric floating half zone convection in a fat liquid bridge of larger Prandtl number

The linear instability analysis of the present paper shows that the thermocapillary convection in a half floating zone of larger Prandtl number has a steady instability mode w(i) = 0 and m = 1 for a fat liquid bridge V = 1.2 with small geometrical aspect ratio A = 0.6. This conclusion is different from the usual idea of hydrothermal instability, and implies that the instability of the system may excite a steady and axial asymmetric state before the onset of oscillation in the ease of large Prandtl number.

Prediction of structural dynamic plastic shear failure by Johnson's damage number

It is proved that Johnson's damage number is the sole similarity parameter for dynamic plastic shear failure of structures loaded impulsively, therefore, dynamic plastic shear failure can be understood when damage number reaches a critical value. It is suggested that the damage number be generally used to predict the dynamic plastic shear failure of structures under various kinds of dynamic loads (impulsive loading, rectangular pressure pulse, exponential pressure pulse, etc.,). One of the advantages for using the damage number to predict such kind of failure is that it is conveniently used for dissimilar material modeling.

Two Bifurcation Transition Processes in Floating Half Zone Convection of Larger Prandtl Number Fluid

Processes of the onset oscillation in the thermocapillaxy convection under the Earth's gravity are investigated by the numerical simulation and experiments in a floating half zone of large Prandtl number with different volume ratio. Both computational and experimental results show that the steady and axisymmetric convection turns to the oscillatory convection of m=1 for the slender liquid bridge, and to the oscillatory convection before a steady and 3D asymmetric state for the case of a fat liquid bridge. It implies that, there are two critical Marangoni numbers related, respectively, to these two bifurcation transitions for the fat liquid bridge. The computational results agree with the results of ground-based experiments.