R410A condensation inside a commercial brazed plate heat exchanger

This paper presents the heat transfer coefficients and the pressure drop measured during HFC-410A condensation inside a commercial brazed plate heat exchanger: the effects of saturation temperature, refrigerant mass flux and vapour super-heating are investigated. The heat transfer coefficients show weak sensitivity to saturation temperature and great sensitivity to refrigerant mass flux and vapour super-heating. At low refrigerant mass flux (<20 kg/m² s) the saturated vapour condensation heat transfer coefficients are not dependent on mass flux and are well predicted by Nusselt [W. Nusselt, Die oberflachenkondensation des wasserdampfes, Energy 60 (1916) 541–546, 569–575] analysis for vertical surface: the condensation process is gravity controlled. For higher refrigerant mass flux (>20 kg/m²s) the saturated vapour condensation heat transfer coefficients depend on mass flux and are well predicted by Akers et al. [W.W. Akers, H.A. Deans, O.K. Crosser, Condensing heat transfer within horizontal tubes, Chem. Eng. Prog. Symp. Series 55 (1959) 171–176] equation: forced convection condensation occurs. In the forced convection condensation region the heat transfer coefficients show a 30% increase for a doubling of the refrigerant mass flux. The condensation heat transfer coefficients of super-heated vapour are 8–10% higher than those of saturated vapour and are well predicted by Webb [R.L. Webb, Convective condensation of superheated vapor, ASME J. Heat Transfer 120 (1998) 418–421] model. A simple linear equation based on the kinetic energy per unit volume of the refrigerant flow is proposed for the computation of the frictional pressure drop.     

Impactless biped walking on a slope

Walking without impacts has been considered in dynamics as a motion/force control problem. In order to avoid impacts, an approach for both the specified motion of the biped and its ground reaction forces was presented yielding a combined motion and force control problem. As an application, a walker on a horizontal plane has been considered. In this paper, it is shown how the control of the ground reaction forces and the energy consumption depend on the gradient of a slope. The biped dynamics and the constraints within the biped system and on the ground are discussed. A motion control synthesis is developed using the inverse dynamics principle proven to be most efficient for human walking research, too. The impactless walking with controlled legs is illustrated by a seven-link biped. The “flying” biped has nine degrees of freedom, with six control inputs. During locomotion, the standing leg has three scleronomic constraints, and the trunk has three rheonomic constraints. However, there are three rheonomic constraints for the prescribed leg motion or three scleronomic constraints for reaction forces of the trailing leg, respectively. The nominal control action for impactless walking can be precomputed and stored. The model proposed allows the investigation of several problems: uphill and downhill walking, optimization of step length, stiction of the feet on the slope and many more. All these findings are also of interest in biomechanics.     

A cryomechanics technique to measure dissipated energies of 10 nJ

A new cryomechanics measurement technique has been developed to measure fracture-induced dissipated energies as small as 10 nJ (10×10⁻⁹ J) at temperatures near 4.2 K. The technique, with much less stringent instrumentation requirements than those used for measurement of ∼10 nJ energies, was applied to an induced fracture experiment where dissipation energies were of the order of ∼100 μJ. Fracture of 0.5-mm diameter pencil leads of two different hardnesses gave rise to measured energies of 65 ∼ 110 μJ. A two-dimensional finite-element analysis was used to interpret the experimental measurements. Based on the analysis, approximately 50 μJ of 65 ∼ 110 μJ measured is estimated to be the dissipated energy associated with crack formation and propagation.     

Hybrid dynamics of two coupled oscillators that can impact a fixed stop

We consider two linearly coupled masses, where one mass can have inelastic impacts with a fixed, rigid stop. This leads to the study of a two degree of freedom, piecewise linear, frictionless, unforced, constrained mechanical system. The system is governed by three types of dynamics: coupled harmonic oscillation, simple harmonic motion and discrete rebounds. Energy is dissipated discontinuously in discrete amounts, through impacts with the stop. We prove the existence of a non-zero measure set of orbits that lead to infinite impacts with the stop in a finite time. We show how to modify the mathematical model so that forward existence and uniqueness of solutions for all time is guaranteed. Existence of hybrid periodic orbits is shown. A geometrical interpretation of the dynamics based on action coordinates is used to visualize numerical simulation results for the asymptotic dynamics.     

In situ observations of unusual drop deformation and wobbling in simple shear flow

Deformation and wobbling of a liquid drop immersed in a liquid matrix were studied under mild shear conditions for various viscosity ratios. In situ visualization experiments were conducted on a homemade transparent Couette cell incorporated to the Paar Physica MCR500 shear rheometer. The effect of drop or matrix elasticity was examined and was found to play a major role in both deformation and wobbling processes. Experimental results were compared to Jackson and Tucker (J Rheol 47:659–682, 2003), Maffettone and Minale (J Non-Newton Fluid Mech 78:227–241, 1998) and Yu and Bousmina (J Rheol 47:1011–1039, 2003) ellipsoidal models. It was found that the agreement between the Newtonian models and the experimental results required an increase in the drop viscosity. Such increment in viscosity was found to scale with the first normal stress difference.     

An experimental study of the unsteady vortex structures in the wake of a root-fixed flapping wing

An experimental study was conducted to characterize the evolution of the unsteady vortex structures in the wake of a root-fixed flapping wing with the wing size, stroke amplitude, and flapping frequency within the range of insect characteristics for the development of novel insect-sized nano-air-vehicles (NAVs). The experiments were conducted in a low-speed wing tunnel with a miniaturized piezoelectric wing (i.e., chord length, C = 12.7 mm) flapping at a frequency of 60 Hz (i.e., f = 60 Hz). The non-dimensional parameters of the flapping wing are chord Reynolds number of Re = 1,200, reduced frequency of k = 3.5, and non-dimensional flapping amplitude at wingtip h = A/C = 1.35. The corresponding Strouhal number (Str) is 0.33, which is well within the optimal range of 0.2 < Str < 0.4 used by flying insects and birds and swimming fishes for locomotion. A digital particle image velocimetry (PIV) system was used to achieve phased-locked and time-averaged flow field measurements to quantify the transient behavior of the wake vortices in relation to the positions of the flapping wing during the upstroke and down stroke flapping cycles. The characteristics of the wake vortex structures in the chordwise cross planes at different wingspan locations were compared quantitatively to elucidate underlying physics for a better understanding of the unsteady aerodynamics of flapping flight and to explore/optimize design paradigms for the development of novel insect-sized, flapping-wing-based NAVs.     

Two-phase co-current flow simulations using periodic boundary conditions in horizontal, 4, 10 and 90° inclined eccentric annulus,flow prediction using a modified interFoam solver and comparison with experimental results

Two-phase oil and gas flow were simulated in an entirely eccentric annulus and compared with experimental data at horizontal, 4, 10, and 90° inclination. The gas-phase was sulphur hexafluoride and the liquid phase a mixture of Exxsol D60 and Marcol 82 for the inclined cases (5–16), and pure Exxsol D60 for the horizontal cases (1–4). The diameter of the outer and inner cylinders was 0.1 and 0.04 m, respectively, for the inclined domains and 0.1 and 0.05 m for the horizontal domain. The cases studied consist of liquid phase fractions between 0.3 and 0.65 and mixture velocities from 1.2 to 4.25 m/s. The mean pressure gradient is within 33% of the expected experimental behavior for all inclined cases. In contrast, the low-velocity horizontal domains exhibit significant deviation, with a drastic over-prediction of the mean pressure gradient by as much as 200–335% for cases 1 and 2. The two remaining horizontal cases (3 and 4) are within 22% of the expected mean pressure gradient. Cases 13–16 are a replication of cases 5–8 at an increased inclination; the mean pressure gradient is within 6.5% of the expected increase due to the increase in hydrostatic pressure. By comparing cases 1–4 to previous published simulations at a lower eccentricity, we found a decrease of the mean pressure gradient by 30–40%, which is in line with existing literature, although for single-phase flows. The simulated and experimental liquid holdup profiles are in good agreement when comparing the fractional data; wave and slug frequencies match to within 0.5 Hz; however, at closer inspection, it is apparent that there is a decrease in the amount of phase-mixing of the simulations compared to the experiments. When increasing the mesh density from 115 k cells/m to 2 million cells/m, the simulations exhibit significantly more phase mixing, but are still unable to produce conventional slugs. In a simplified case, conventional slugs are observed at grid sizing of 1 × 1 × 1 mm, whereas the cells of the 2 million cells/m mesh are roughly 1.5 × 1.5 × 1.5 mm.     

Performance of an oscillating subsoiler in breaking a hardpan

A single shank tractor mounted oscillating subsoiler was developed to break hardpan, common in sugarcane (Saccharum officinarum) farms especially after harvest when heavy trucks transport the cut canes from the field to the sugar factory. Field experiments were conducted to determine the optimum combination of performance parameters of the subsoiler. Field tests were conducted at frequencies of oscillation of 3.7, 5.67, 7.85, 9.48 and 11.45 Hz; amplitudes of 18, 21, 23.5, 34 and 36.5 mm; and forward speeds of 1.85, 2.20 and 3.42 km h⁻¹ at moisture contents close to the lower plastic limit of the clay soil. A reduction in draft but an increase in total power requirement was found for oscillating compared to non-oscillating subsoiler. The draft and power ratios were significantly affected by the forward speed, frequency and amplitude. Their combined interaction, expressed in terms of the velocity ratio (the ratio of peak tool velocity to forward speed), however, had the strongest influence. At the same velocity ratio, the draft reduction and power increase were less at higher amplitude of oscillation. For the field conditions tested, the optimum operation for least energy expenditure was obtained at an amplitude of 36.5 mm, frequency of 9.48 Hz and speed of 2.20 km h⁻¹ with a draft ratio of 0.33 and power ratio of only 1.24. It could be concluded that the oscillating subsoiler reduces draft for breaking hardpan, reduces soil compaction and promotes the use of lighter tractors by utilizing tractor power-take-off (p.t.o.) power to achieve higher efficiency of power transmission. ©     

Decaying homogeneous turbulence under strong density stratification at low-Prandtl number

The characteristics of decaying homogeneous turbulence under strong density stratification have been studied using direct numerical simulations. While our previous study dealt with rotating stratified turbulence, here we investigate the detailed flow structure of stratified turbulence without rotation especially at low-Prandtl number. By assuming a low-Prandtl-number fluid, e.g. liquid sodium: Pr ≈ 0.01, gallium: Pr ≈ 0.025, internal gravity waves are markedly attenuated due to the large thermal conductivity, and turbulence soon reaches a two-component state, where vertical energy, coupled with potential energy, significantly decays, and becomes negligible as observed experimentally (Praud et al. in J Fluid Mech 522:1–33, 2005). In the horizontal plane, there appear large-scale vortices with vertical vorticity, and those with the same sign of vorticity increase their horizontal length scale by merging with each other. In the vertical plane, highly sheared regions represented by horizontal vorticity also tend to horizontally increase their length scale and become layered structures by the combined effects of vortex coalescence and energy cascade into higher vertical wavenumbers.      

Measurement method of complex viscoelastic material properties

A measurement technique of viscoelastic properties of polymers is proposed to investigate complex Poisson’s ratio as a function of frequency. The forced vibration responses for the samples under normal and shear deformation are measured with varying load masses. To obtain modulus of elasticity and shear modulus, the present method requires only knowledge of the load mass, geometrical characteristics of a sample, as well as both the amplitude ratio and phase lag of the forcing and response oscillations. The measured data were used to obtain the viscoelastic properties of the material based on a 2D numerical deformation model of the sample. The 2D model enabled us to exclude data correction by the empirical form factor used in 1D model. Standard composition (90% PDMS polymer + 10% catalyst) of silicone RTV rubber (Silastic^® S2) were used for preparing three samples for axial stress deformation and three samples for shear deformation. Comprehensive measurements of modulus of elasticity, shear modulus, loss factor, and both real and imaginary parts of Poisson’s ratio were determined for frequencies from 50 to 320 Hz in the linear deformation regime (at relative deformations 10^?6 to 10^?4) at temperature 25 °C. In order to improve measurement accuracy, an extrapolation of the obtained results to zero load mass was suggested. For this purpose measurements with several masses need to be done. An empirical requirement for the sample height-to-radius ratio to be more than 4 was found for stress measurements. Different combinations of the samples with different sizes for the shear and stress measurements exhibited similar results. The proposed method allows one to measure imaginary part of the Poisson’s ratio, which appeared to be about 0.04–0.06 for the material of the present study.     

Damping Characterization and Viscoelastic Behavior of Laminated Polymer Matrix Composites Using a Modified Classical Lamination Theory

The damping property of materials can be defined as the ratio of dissipated energy over the total strain energy during the loading–unloading process, called the specific damping capacity (SDC). In this study, in order to characterize the damping properties of materials, a test plan in designed to extract the SDC of a single layer composite from hysteresis data. The theory of linear viscoelasticity incorporates a varying Young’s Modulus by using a complex stiffness modulus. Considering different lay-ups, the modified classical lamination plate theory is modified to represent both stiffness and SDC of laminates. The results are compared with experimental results for symmetric laminated specimen. This evaluation shows a very good agreement between theoretical and experimental results in the range of low frequency loading from 0.2 to 4 Hz. The complex compliance matrix changes the governing equation in to a complex form which contains both stiffness and damping properties.     

Zonal fracturing mechanism in deep crack-weakened rock masses

The mechanical behaviors of deep crack-weakened rock masses are different from those of shallow crack-weakened rock masses. The surrounding rock in shallow crack-weakened rock mass engineering is classified into loose zone, plastic zone and elastic zone, while the surrounding rock in deep crack-weakened rock mass engineering is classified into fractured zone and non-fractured zone, which occur alternatively. It is assumed that the deep rock masses contain one joint set, in which the probability density function describing the distribution of sizes is assumed to follow the Rayleigh distribution, and the probability density function describing the distribution of spacing is assumed to follow the Weibull distribution. On the basis of strength criterion of deep rock mass, the near-field stress redistribution around circular opening induced by excavation is determined. The strong interaction among cracks is investigated by using the dislocation model. The nucleation, growth, interaction and coalescence of cracks were analyzed based on the strain energy density factor theory. When cracks coalesce, failure of deep crack-weakened rock masses occurs, fractured zone is formed. Then, size and quantity of fractured zone and non-fractured zone are given out. The size and quantity of fractured zone increase with decreasing strength of rock mass. The size and quantity of fractured zone increase with increasing in situ stress. Zonal fracturing phenomenon occurs once value of in situ stress is larger than the unaxial compressive strength of rock masses. The size and quantity of fractured zone decrease with increasing λ when p₂ > p₁. The size and quantity of fractured zone increase with increasing λ when p₂ < p₁.     

Very high Reynolds number boundary layers over 3D sparse roughness and obstacles: the mean flow

Hot-wire and oil-film interferometry measurements are taken for 3D rough wall boundary layers at very high Reynolds numbers (61,000 < Re θ < 120,000) with low blockage ratios, 10 < δ/H < 135, and high roughness, 100 < H ⁺ < 4,900. The results cover flows over both rough walls and over obstacles and are compared with and provide extension to recent lower Reynolds number results. The validity of the Townsend ‘wall similarity hypothesis’ in relation to consistently increasing 3D roughness is interrogated. In agreement with recent work, Schultz and Flack (J Fluid Mech 580:381–405, 2007) and Castro (J Fluid Mech 585:469–485, 2007) found that, for relatively low roughness, Townsend’s hypothesis holds for the mean velocity field. With increasing roughness, the equilibrium layer diminishes and gradually vanishes. The viscous component of the wall shear stress decreases, while the turbulent component increases as the roughness effects extend across the boundary layer.     

Plasma jets produced by low energy laser pulse interaction with planar and cratered targets

Laser experiments of the plasma jet formation using nanosecond laser pulses with low energy, i.e., <20 J, are presented. Planar and cratered gadolinium and aluminum targets are irradiated with laser intensities of several 10¹⁴ W/cm². Spatially-resolved time-integrated X-ray spectra were recorded in the spectral range from 7 to 10 Å. A jet-like structure is obtained from aluminum targets with a preformed crater, which is not seen in planar target irradiation. For gadolinium, a jet is observed from both planar and preformed cratered targets, suggesting that the collimation is dominated by radiative cooling. A radiation-hydrodynamics code coupled to a non-LTE ionization code was used to model the plasma. The calculated plasma emission was found to be consistent with the experimental results.     

The load distribution among three legs on the wall: model predictions for cockroaches

In a former study on terrestrial locomotion of cockroaches in the sagittal plane, it was hypothesised that the ground reaction force distribution among three legs synchronously in contact with a substrate is predominantly explained by joint torque minimisation within all three legs. We verified this hypothesis with a simple mechanical model in two dimensions, consisting of one body and three mass-less legs. Hereto, we calculated force distributions resulting from different optimisation criteria for varying slope angles of the substrate. We compared these distributions to each other and the few experimental findings available. We found that, for any slope angle, the force distribution rather seems to be derived from the fundamental “table” solution, i.e. equalised vertical and vanishing horizontal components (equivalent to pure force minimisation at zero slope), than from pure torque minimisation. For cockroaches, the “table” solution is likely to be modified by torque minimisation within the leading and the trailing leg. We demonstrate that the minimisation of leg forces and of interaction forces is fully equivalent. Moreover, our model predicts the force distribution for arbitrary slope angles. Based on our model calculations, we speculate that in terrestrial locomotion, some animals may rely on spring-mass model dynamics whatever slope angle to be overcome. This might only become evident when focusing movement analyses strictly on a gravity rather than on a substrate-based coordinate system.     

Visualization of natural convection heat transfer on horizontal cylinders

Natural convection heat transfer phenomena on horizontal cylinders were investigated experimentally in order to explore the applicability of analogy experimental method using the copper electroplating system and to visualize the local heat transfer depending on the angular position and the diameter of the horizontal cylinder. The diameters of the cylinders are varied from 0.01 to 0.15 m, which correspond to the Rayleigh numbers of 1.73 × 10⁷–5.69 × 10¹¹. The measured mass transfer coefficients show good agreements with the existing heat transfer correlations. The patterns of copper plated on the aluminum cathodes for various Rayleigh numbers reveal and visualize the local heat transfer depending on the angular position and show good agreement with the works of Kitamura et al. The hydrogen bubbles produced at higher applied potential visualize the plumes appeared on top region of the cylinders.     

Mass and isotopic concentrations of water-insoluble refractory carbon in total suspended particulates at Mt. Waliguan Observatory (China)

Mass concentration and isotopic values δ¹³C and ¹⁴C are presented for the water-insoluble refractory carbon (WIRC) component of total suspended particulates (TSP), collected weekly during 2003, as well as from October 2005 to May 2006 at the WMO-GAW Mt. Waliguan (WLG) site. The overall average WIRC mass concentration was (1183 ± 120) ng/m³ (n = 79), while seasonal averages were 2081 ± 1707 (spring), 454 ± 205 (summer), 650 ± 411 (autumn), and 1019 ± 703 (winter) ng/m³. Seasonal variations in WIRC mass concentrations were consistent with black carbon measurements from an aethalometer, although WIRC concentrations were typically higher, especially in winter and spring. The δ¹³C PDB value (−25.3 ± 0.8)‰ determined for WIRC suggests that its sources are C₃ biomass or fossil fuel combustion. No seasonal change in δ¹³C PDB was evident. The average percent Modern Carbon (pMC) for ¹⁴C in WIRC for winter and spring was (67.2 ± 7.7)% (n = 29). Lower pMC values were associated with air masses transported from the area east of WLG, while higher pMC values were associated with air masses from the Tibetan Plateau, southwest of WLG. Elevated pMC values with abnormally high mass concentrations of TSP and WIRC were measured during a dust storm event.     

A numerical study on the absorption of water vapor into a film of aqueous LiBr falling along a vertical plate

Absorber is an important component in absorption machines and its characteristics have significant effects on the overall efficiency of absorption machines. This article reports a model of simultaneous heat and mass transfer process in absorption of refrigerant vapor into a lithium bromide solution of water––cooled vertical plate absorber in the Reynolds number range of 5 < Re < 150. The boundary layer assumptions were used for the transport of mass, momentum and energy equations and the fully implicit finite difference method was employed to solve the governing equations in the film flow. Dependence of lithium bromide aqueous properties to the temperature and concentration and film thickness to vapor absorption was employed. This model can predict temperature, concentration and properties of aqueous profiles as well as the absorption heat and mass fluxes, heat and mass transfer coefficients, Nusslet and Sherwood number of absorber. An analysis for linear distribution of wall temperature condition carries out to investigation the reliability of the present numerical method through comparing with previous investigation.     

Transport properties in mixtures involving carbon dioxide at low and moderate density: test of several intermolecular potential energies and comparison with experiment

It is the purpose of this paper to extract unlike intermolecular potential energies of five carbon dioxide-based binary gas mixtures including CO₂–He, CO₂–Ne, CO₂–Ar, CO₂–Kr, and CO₂–Xe from viscosity data and compare the calculated potentials with other models potential energy reported in literature. Then, dilute transport properties consisting of viscosity, diffusion coefficient, thermal diffusion factor, and thermal conductivity of aforementioned mixtures are calculated from the calculated potential energies and compared with literature data. Rather accurate correlations for the viscosity coefficient of afore-cited mixtures embracing the temperature range 200 K < T < 3273.15 K is reproduced from the present unlike intermolecular potentials energy. Our estimated accuracies for the viscosity are to within ±2%. In addition, the calculated potential energies are used to present smooth correlations for other transport properties. The accuracies of the binary diffusion coefficients are of the order of ±3%. Finally, the unlike interaction energy and the calculated low density viscosity have been employed to calculate high density viscosities using Vesovic–Wakeham method.     

Numerical modeling of incline plate LiBr absorber

Among major components of LiBr–H₂O absorption chillers is the absorber, which has a direct effect on the chillier size and whose characteristics have significant effects on the overall efficiency of absorption machines. In this article, heat and mass transfer process in absorption of refrigerant vapor into a lithium bromide solution of water-cooled incline plate absorber in the Reynolds number range of 5 < Re < 150 is performed numerically. The boundary layer assumptions are used for the mass, momentum and energy transport equations and the fully implicit finite difference method is employed to solve the governing equations. Dependence of lithium bromide aqueous properties to the temperature and concentration is employed as well as dependence of film thickness to vapor absorption. An analysis for linear distribution of wall temperature condition carries out to investigate the reliability of the present numerical method through comparing with previous investigation. The effect of plate angle on heat and mass transfer parameters is investigated and the results show that absorption mass flux and heat and mass transfer coefficient increase as the angle of the plate increase. The main parameters of absorber design, namely Nusselt and Sherwood numbers, are correlated as a function of Reynolds Number and the plate angle.