Natural convection from a plane vertical porous surface in non-isothermal surroundings

In this paper, laminar natural convection flow from a permeable and isothermal vertical surface placed in non-isothermal surroundings is considered. Introducing appropriate transformations into the boundary layer equations governing the flow derives non-similar boundary layer equations. Results of both the analytical and numerical solutions are then presented in the form of skin-friction and Nusselt number. Numerical solutions of the transformed non-similar boundary layer equations are obtained by three distinct solution methods, (i) the perturbation solutions for small � (ii) the asymptotic solution for large � (iii) the implicit finite difference method for all � where � is the transpiration parameter. Perturbation solutions for small and large values of � are compared with the finite difference solutions for different values of pertinent parameters, namely, the Prandtl number Pr, and the ambient temperature gradient n.

Scaling of thermo-magnetic convection

The scaling to characterize unsteady boundary layer development for thermo-magnetic convection of paramagnetic fluids with the Prandtl number greater than one is developed. Under the consideration is a square cavity with initially quiescent isothermal fluid placed in microgravity condition (g = 0) and subject to a uniform, vertical gradient magnetic field. A distinct magnetic thermal-boundary layer is produced by sudden imposing of a higher temperature on the vertical sidewall and as an effect of magnetic body force generated on paramagnetic fluid. The transient flow behavior of the resulting boundary layer is shown to be described by three stages: the start-up stage, the transitional stage and the steady state. The scaling is verified by numerical simulations with the magnetic momentum parameter m variation and the parameter γRa variation.

Thermo-magnetic convection of paramagnetic fluids with non-instantaneous heating

The unsteady boundary-layer development for thermomagnetic convection of paramagnetic fluids inside a square cavity has been considered in this study. The cavity is placed in a microgravity condition (no gravitation acceleration) and under a uniform magnetic field which acts vertically. A ramp temperature boundary condition is applied on left vertical side wall of the cavity where the temperature initially increases with time up to some specific time and maintain constant thereafter. A distinct magnetic convection boundary layer is developed adjacent to the left vertical wall due to the effect of the magnetic body force generated on the paramagnetic fluid. An improved scaling analysis has been performed using triple-layer integral method and verified by numerical simulations. The Prandtl number has been chosen greater than unity varied over 5-100. Moreover, the effect of various values of the magnetic parameter and magnetic Rayleigh number on the fluid flow and heat transfer has been shown.

The effects of cold air (-110°C) and water (8°C) cryotherapy on knee skin temperature

Exposure to cold air, whole body cryotherapy (WBC), is a novel treatment employed by athletes. In WBC individuals dressed in minimal clothing are exposed to a temperature below -100°C for 2-4 min. The use of WBC has been advocated as a treatment for various knee injuries. PURPOSE: To compare the effects of two modalities of cryotherapy, -110°C WBC and 8°C cold water immersion (CWI) on knee skin temperature (Tsk). METHODS: With ethical approval and written informed consent 10 healthy active male participants (26.5±4.9 yr, 183.5±6.0 cm, 90.7±19.9 kg, 26.8±5.0 kg/m2, 23.0±9.3% body fat (measured by DXA), 7.6 ± 2.0 mm patellar skin fold; mean±SD) were exposed to 4 min of CWI and WBC. The treatment order was randomised in a controlled crossover design, with a minimum of 7 days between treatments. During WBC participants stood in a chamber (-60±3°C) for 20 s before entering the main chamber (-110°C±3°C) where they remained for 3 min and 40 s. For CWI participants were seated in a tank filled with cold water (8±0.3°C) and immersed to the level of the sternum for 4 min. Right knee Tsk was assessed via non-contact, infrared thermal imaging. A quadrilateral region of interest was created using inert markers placed 5 cm above and below the most superior and inferior aspect of the patella. Tsk within this quadrilateral was recorded pre, immediately post and every 10 min thereafter for 60 min. Tsk changes were examined using a two-way (treatment x time) repeated measures analyses of variance. In addition, a paired sample t-test was used to compare baseline Tsk before both treatments. RESULTS: Knee Tsk was similar before treatment (WBC: 29.9±0.7°C, CWI: 29.6±0.9°C, p>0.05). There was a significant main effect for treatment (p<0.05) and time (p<0.001). Compared to baseline, Tsk was significantly reduced (p<0.05) immediately post and at 10, 20, 30, 40, 50 and 60 min after both cooling modalities. Knee Tsk was lower (p<0.05) immediately after WBC (19.0±0.9°C) compared to CWI (20.5±0.6°C). However, from 10 to 60 min post, knee Tsk was lower (p<0.05) following the CWI treatment. CONCLUSION: WBC elicited a greater decrease in knee Tsk compared to CWI immediately after treatment. However, both modalities display different recovery patterns and Tsk after CWI was significantly lower than WBC at 10, 20, 30, 40, 50 and 60 min after treatment.

Prospective decrease in progesterone concentrations in female lightweight rowers during the competition season compared with the off-season: A controlled study examining weight loss and intensive exercise

The purpose of this study was to monitor ovarian hormone function response to intense exercise and body weight changes in female athletes. Ovarian hormone function was evaluated in 12 female lightweight rowers and 10 age-height-weight matched sedentary controls. Ovarian hormone function was assessed during consecutive competition season and off season, by measurement of peak and average alternative day overnight urinary oestrone glucuronide (E1G) and pregnanediol glucuronide (PdG) excretion. Competition season was associated with a 5.8 kg (9.3%) body weight loss in the lightweight rowers. Significantly lower competition season peak and average urinary excretion of PdG were found in the lightweight rowers compared with the controls. Lower competition season peak and average urinary excretion of E1G were also found in the lightweight rowers compared with the controls, but the difference did not reach significance. The number of rowing training hours was a significant determinant of peak PdG excretion in the rowers (R2 = 0.40; p<0.02). The seasonal suppression of PdG excretion was associated with degree of weight loss (R2 = 0.46; p<0.01). The competition related decrease in E1G and PdG excretion for the lightweight rowers was predominantly restored during the off season when exercise intensity and duration were decreased and body weight increased. These results showed a significant (p<0.05) reduction in progesterone metabolite excretion and a non-significant decrease in oestrone metabolite excretion associated with intensive competition season training loads and body weight reduction in female lightweight rowers.

Computation of Thermal Stress Intensity Factors for Bimaterial Interface Cracks Using Domain Integral Method

The concept of domain integral used extensively for J integral has been applied in this work for the formulation of J(2) integral for linear elastic bimaterial body containing a crack at the interface and subjected to thermal loading. It is shown that, in the presence of thermal stresses, the J(k) domain integral over a closed path, which does not enclose singularities, is a function of temperature and body force. A method is proposed to compute the stress intensity factors for bimaterial interface crack subjected to thermal loading by combining this domain integral with the J(k) integral. The proposed method is validated by solving standard problems with known solutions.

Adhesion-Induced Instability in Asperities

Adhesive forces between two approaching asperities will deform the asperities, and under certain conditions this will result in a sudden runaway deformations leading to a jump-to-contact instability. We present finite element-based numerical studies on adhesion-induced deformation and instability in asperities. We consider the adhesive force acting on an asperity, when it is brought near a rigid half-space, due to van der Waals interaction between the asperity and the half-space. The adhesive force is considered to be distributed over the volume of the asperity (body force), thus resulting in more realistic simulations for the length scales considered. Iteration scheme based on a ``residual stress update'' algorithm is used to capture the effect of deformation on the adhesion force, and thereby the equilibrium configuration and the corresponding force. The numerical results are compared with the previous approximate analytical solutions for adhesion force, deformation of the asperity and adhesion-induced mechanical instability (jump-to-contact). It is observed that the instability can occur at separations much higher,and could possibly explain the higher value of instability separation observed in experiments. The stresses in asperities, particularly in case of small ones, are found to be high enough to cause yielding before jump -to-contact. The effect of roughness is considered by modeling a spherical protrusion on the hemispherical asperity.This small-scale roughness at the tip of the asperities is found to control the deformation behavior at small separations, and hence are important in determining the friction and wear due to the jump-to-contact instability.

Seismic Bearing Capacity of Foundations on Cohesionless Slopes

By applying the lower bound theorem of limit analysis in conjunction with finite elements and nonlinear optimization, the bearing capacity factor N has been computed for a rough strip footing by incorporating pseudostatic horizontal seismic body forces. As compared with different existing approaches, the present analysis is more rigorous, because it does not require an assumption of either the failure mechanism or the variation of the ratio of the shear to the normal stress along the footing-soil interface. The magnitude of N decreases considerably with an increase in the horizontal seismic acceleration coefficient (kh). With an increase in kh, a continuous spread in the extent of the plastic zone toward the direction of the horizontal seismic body force is noted. The results obtained from this paper have been found to compare well with the solutions reported in the literature. (C) 2013 American Society of Civil Engineers.

A finite element method for the Saint-Venant torsion and bending problems for prismatic beams

In this work, we present a finite element formulation for the Saint-Venant torsion and bending problems for prismatic beams. The torsion problem formulation is based on the warping function, and can handle multiply-connected regions (including thin-walled structures), compound and anisotropic bars. Similarly, the bending formulation, which is based on linearized elasticity theory, can handle multiply-connected domains including thin-walled sections. The torsional rigidity and shear centers can be found as special cases of these formulations. Numerical results are presented to show the good coarse-mesh accuracy of both the formulations for both the displacement and stress fields. The stiffness matrices and load vectors (which are similar to those for a variable body force in a conventional structural mechanics problem) in both formulations involve only domain integrals, which makes them simple to implement and computationally efficient. (C) 2014 Elsevier Ltd. All rights reserved.

Current kinematics and dynamics of Africa and the East African Rift System

Although the East African Rift System (EARS) is an archetype continental rift, the forces driving its evolution remain debated. Some contend buoyancy forces arising from gravitational potential energy (GPE) gradients within the lithosphere drive rifting. Others argue for a major role of the diverging mantle flow associated with the African Superplume. Here we quantify the forces driving present-day continental rifting in East Africa by (1) solving the depth averaged 3-D force balance equations for 3-D deviatoric stress associated with GPE, (2) inverting for a stress field boundary condition that we interpret as originating from large-scale mantle tractions, (3) calculating dynamic velocities due to lithospheric buoyancy forces, lateral viscosity variations, and velocity boundary conditions, and (4) calculating dynamic velocities that result from the stress response of horizontal mantle tractions acting on a viscous lithosphere in Africa and surroundings. We find deviatoric stress associated with lithospheric GPE gradients are similar to 8-20 MPa in EARS, and the minimum deviatoric stress resulting from basal shear is similar to 1.6 MPa along the EARS. Our dynamic velocity calculations confirm that a force contribution from GPE gradients alone is sufficient to drive Nubia-Somalia divergence and that additional forcing from horizontal mantle tractions overestimates surface kinematics. Stresses from GPE gradients appear sufficient to sustain present-day rifting in East Africa; however, they are lower than the vertically integrated strength of the lithosphere along most of the EARS. This indicates additional processes are required to initiate rupture of continental lithosphere, but once it is initiated, lithospheric buoyancy forces are enough to maintain rifting.

The Uniaxial Tensile Deformation of Ni Nanowire: Atomic-Scale Computer Simulations

The mechanical deformations of nickel nanowire subjected to uniaxial tensile strain at 300 K are simulated by using molecular dynamics with the quantum corrected Sutten-Chen many-body force field. We have used common neighbor analysis method to investigate the structural evolution of Ni nanowire during the elongation process. For the strain rate of 0.1%/ps, the elastic limit is up to about 11% strain with the yield stress of 8.6 GPa. At the elastic stage, the deformation is carried mainly through the uniform elongation of the distances between the layers (perpendicular to the Z-axis) while the atomic structure remains basically unchanged. With further strain, the slips in the {111} planes start to take place in order to accommodate the applied strain to carry the deformation partially, and subsequently the neck forms. The atomic rearrangements in the neck region result in a zigzag change in the stress-strain curve; the atomic structures beyond the region, however, have no significant changes. With the strain close to the point of the breaking, we observe the formation of a one-atom thick necklace in Ni nanowire. The strain rates have no significant effect on the deformation mechanism, but have some influence on the yield stress, the elastic limit, and the fracture strain of the nanowire.

Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

The structure of liquid argon as determined by x-ray diffraction

X-ray diffraction measurements and subsequent data analyses have been carried out on liquid argon at five states in the density range of 0.91 to 1.135 gm/cc and temperature range of 127 to 143°K. Duplicate measurements were made on all states. These data yielded radial distribution and direct correlation functions which were then used to compute the pair potential using the Percus-Yevick equation. The potential minima are in the range of -105 to -120°K and appear to substantiate current theoretical estimates of the effective pair potential in the presence of a weak three-body force.

The data analysis procedure used was new and does not distinguish between the coherent and incoherent absorption factors for the cell scattering which were essentially equal. With this simplification, the argon scattering estimate was compared to the gas scattering estimate on the laboratory frame of reference and the two estimates coincided, indicating the data normalized. The argon scattering on the laboratory frame of reference was examined for the existence of the peaks in the structure factor and the existence of an observable third peak was considered doubtful.

Numerical studies of the effect of truncation, normalization, the subsidiary peak phenomenon in the radial distribution function, uncertainties in the low angle data relative to errors in the direct correlation function and the distortion phenomenon are presented.

The distortion phenomenon for this experiment explains why the Mikolaj-Pings argon data yielded pair potential well depths from the Percus-Yevick equation that were too shallow and an apparent slope with respect to density that was too steep compared to theoretical estimates.

The data presented for each measurement are: empty cell and cell plus argon intensity, absorption factors, argon intensity, smoothed argon intensity, smoothed argon intensity corrected for distortion, structure factor, radial distribution function, direct correlation function and the pair potential from the Percus-Yevick equation.

Effect of crude oil and some petroleum products on Clarias Gariepinus fingerlings (Catfish: Claridae)

Ninety (90) hatchery bred fingerlings of Clarias gariepinus (mean weight: 0.96 ± 0.1g) were randomly placed in 15 plastic baths (25 litres each) at the Research laboratory and were exposed to different concentrations of oil products to determine their effects on the fish, to facilitate inferential deductions that will enhance effective aquatic environmental management. Three (3) replicate basins of 5 experimental treatments (crude oil, petrol oil, kerosene oil, engine oil and control) were used at a concentration of 1.25ml. L-1. The control experiment was devoid of oil treatment. Six (6) fingerlings were placed in each replicate basin, flooded with 20 litres of clean tap water and fed with nutrafin cichilid food, 2 times daily at 3% body weight. The results showed that the feeding behaviour and swimming performances of fish were reduced after 24 hours of the addition of the various oil pollutants. Mortality of fingerlings in the oiled basins increased as the hours of exposure increased (i.e. 24, 48, 72 and 96 hours). Recovery was not immediate in the treated basin while surviving fingerlings in the control basins grew up to post-fingerlings after 90 days (3 months). There were significant differences (P<0.01 and P<0.05) in the effect of crude oil and the petroleum products on the mortality rate of C. gariepinus when exposed to oil pollutants at 1.25ml. L-1 concentration

The performance of ejectors driven by sinusoidally unsteady jets

The operation of ejectors driven by a low-speed, sinusoidally unsteady jet has been studied. The thrust augmentation is shown to be highly dependent on the non-dimensional frequency of the driver jet, but independent of its Mach Number. Convective rather than acoustically propagated phenomena dominate the ejector flowfield. Unsteady pressure measurements on the internal surfaces of the ejector have enabled convecting ring vortices to be identified. The impingement of a ring vortex on the leading edge of the ejector causes the peak unsteady body force. The non-dimensional diameter of the ejector is shown to be the only geometric variable that affects the optimum non-dimensional frequency for thrust augmentation. An experimentally optimised geometry is presented. An expression relating the mechanical efficiency and thrust augmentation of the ejector is developed, and shown to be crucially dependent on the degree of unsteadiness in the ejector exit plane.