Oblique-Shock-Wave/Boundary-Layer Interaction Using Near-Wall Reynolds-Stress Models

A synthesis is presented of the predictive capability of a family of near-wall wall-normal free Reynolds stress models (which are completely independent of wall topology, i.e., of the distance fromthe wall and the normal-to-thewall orientation) for oblique-shock-wave/turbulent-boundary-layer interactions. For the purpose of comparison, results are also presented using a standard low turbulence Reynolds number k–ε closure and a Reynolds stress model that uses geometric wall normals and wall distances. Studied shock-wave Mach numbers are in the range MSW = 2.85–2.9 and incoming boundary-layer-thickness Reynolds numbers are in the range Reδ0 = 1–2×106. Computations were carefully checked for grid convergence. Comparison with measurements shows satisfactory agreement, improving on results obtained using a k–ε model, and highlights the relative importance of redistribution and diffusion closures, indicating directions for future modeling work.

Biotribological assessment for artificial synovial joints : the role of boundary lubrication

Biotribology, the study of lubrication, wear and friction within the body, has become a topic of high importance in recent times as we continue to encounter debilitating diseases and trauma that destroy function of the joints. A highly successful surgical procedure to replace the joint with an artificial equivalent alleviates dysfunction and pain. However, the wear of the bearing surfaces in prosthetic joints is a significant clinical problem and more patients are surviving longer than the life expectancy of the joint replacement. Revision surgery is associated with increased morbidity and mortality and has a far less successful outcome than primary joint replacement. As such, it is essential to ensure that everything possible is done to limit the rate of revision surgery. Past experience indicates that the survival rate of the implant will be influenced by many parameters, of primary importance, the material properties of the implant, the composition of the synovial fluid and the method of lubrication. In prosthetic joints, effective boundary lubrication is known to take place. The interaction of the boundary lubricant and the bearing material is of utmost importance. The identity of the vital active ingredient within synovial fluid (SF) to which we owe the near frictionless performance of our articulating joints has been the quest of researchers for many years. Once identified, tribo tests can determine what materials and more importantly what surfaces this fraction of SF can function most optimally with. Surface-Active Phospholipids (SAPL) have been implicated as the body’s natural load bearing lubricant. Studies in this thesis are the first to fully characterise the adsorbed SAPL detected on the surface of retrieved prostheses and the first to verify the presence of SAPL on knee prostheses. Rinsings from the bearing surfaces of both hip and knee prostheses removed from revision operations were analysed using High Performance Liquid Chromatography (HPLC) to determine the presence and profile of SAPL. Several common prosthetic materials along with a novel biomaterial were investigated to determine their tribological interaction with various SAPLs. A pin-on-flat tribometer was used to make comparative friction measurements between the various tribo-pairs. A novel material, Pyrolytic Carbon (PyC) was screened as a potential candidate as a load bearing prosthetic material. Friction measurements were also performed on explanted prostheses. SAPL was detected on all retrieved implant bearing surfaces. As a result of the study eight different species of phosphatidylcholines were identified. The relative concentrations of each species were also determined indicating that the unsaturated species are dominant. Initial tribo tests employed a saturated phosphatidylcholine (SPC) and the subsequent tests adopted the addition of the newly identified major constituents of SAPL, unsaturated phosphatidylcholine (USPC), as the test lubricant. All tribo tests showed a dramatic reduction in friction when synthetic SAPL was used as the lubricant under boundary lubrication conditions. Some tribopairs showed more of an affinity to SAPL than others. PyC performed superior to the other prosthetic materials. Friction measurements with explanted prostheses verified the presence and performance of SAPL. SAPL, in particular phosphatidylcholine, plays an essential role in the lubrication of prosthetic joints. Of particular interest was the ability of SAPLs to reduce friction and ultimately wear of the bearing materials. The identification and knowledge of the lubricating constituents of SF is invaluable for not only the future development of artificial joints but also in developing effective cures for several disease processes where lubrication may play a role. The tribological interaction of the various tribo-pairs and SAPL is extremely favourable in the context of reducing friction at the bearing interface. PyC is highly recommended as a future candidate material for use in load bearing prosthetic joints considering its impressive tribological performance.

Vibration characteristics of plate elements subjected in-plane loads (axial loads)

Axial deformations resulting from in-plane loads (axial forces) of plate elements impact significantly on their vibration characteristics. Although, numerous methods have been developed to quantify axial forces and hence deformations of individual plate elements with different boundary conditions based on their natural frequencies, these methods are unable to apply to the plate elements in a structural system. This is because the natural frequency is a global parameter for the entire structure. Thus, this paper proposes a comprehensive vibration based procedure to quantify axial deformations of plate elements in a structural framing system. Unique capabilities of the proposed method present through illustrative examples. Keywords- Plate Elements, Dynamic Stiffness Matrix, Finite Element Method, Vibration Characteristics, Axial Deformation

Stability boundary analysis of the dynamic voltage restorer in weak systems with dynamic loads

In this contribution, a stability analysis for a dynamic voltage restorer (DVR) connected to a weak ac system containing a dynamic load is presented using continuation techniques and bifurcation theory. The system dynamics are explored through the continuation of periodic solutions of the associated dynamic equations. The switching process in the DVR converter is taken into account to trace the stability regions through a suitable mathematical representation of the DVR converter. The stability regions in the Thevenin equivalent plane are computed. In addition, the stability regions in the control gains space, as well as the contour lines for different Floquet multipliers, are computed. Besides, the DVR converter model employed in this contribution avoids the necessity of developing very complicated iterative map approaches as in the conventional bifurcation analysis of converters. The continuation method and the DVR model can take into account dynamics and nonlinear loads and any network topology since the analysis is carried out directly from the state space equations. The bifurcation approach is shown to be both computationally efficient and robust, since it eliminates the need for numerically critical and long-lasting transient simulations.

Radiation effects on natural convection laminar flow from a horizontal circular cylinder

The effect of radiation on natural convection flow from an isothermal circular cylinder has been investigated numerically in this study. The governing boundary layer equations of motion are transformed into a non-dimensional form and the resulting nonlinear systems of partial differential equations are reduced to convenient boundary layer equations, which are then solved numerically by two distinct efficient methods namely: (i) implicit finite differencemethod or the Keller-Box Method (KBM) and (ii) Straight Forward Finite Difference Method (SFFD). Numerical results are presented by velocity and temperature distribution of the fluid as well as heat transfer characteristics, namely the shearing stress and the local heat transfer rate in terms of the local skin-friction coefficient and the local Nusselt number for a wide range of surface heating parameter and radiation-conduction parameter. Due to the effects of the radiation the skin-friction coefficients as well as the rate of heat transfer increased and consequently the momentum and thermal boundary layer thickness enhanced.

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.

Quantifying in-plane deformation of plate elements using vibration characteristics

Plate elements are used in many engineering applications. In-plane loads and deformations have significant influence on the vibration characteristics of plate elements. Numerous methods have been developed to quantify the effects of in-plane loads and deformations of individual plate elements with different boundary conditions based on their natural frequencies. However, these developments cannot be applied to the plate elements in a structural system as the natural frequency is a global parameter for the entire structure. This highlights the need for a method to quantify in-plane deformations of plate elements in structural framing systems. Motivated by this gap in knowledge, this research has developed a comprehensive vibration based procedure to quantify in-plane deformation of plate elements in a structural framing system. This procedure with its unique capabilities to capture the influence of load migration, boundary conditions and different tributary areas is presented herein and illustrated through examples.

MHD natural convection flow from an isothermal horizontal circular cylinder under consideration of temperature dependent viscosity

Magnetohydrodynamic (MHD) natural convection laminar flow from an iso-thermal horizontal circular cylinder immersed in a fluid with viscosity proportional to a linear function of temperature will be discussed with numerical simulations. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear system of partial differential equa-tions are reduced to convenient form, which are solved numerically by two very efficient methods, namely, (i) Implicit finite difference method together with Keller box scheme and (ii) Direct numerical scheme. Numerical results are presented by velocity and temperature distributions of the fluid as well as heat transfer characteristics, namely the shearing stress and the local heat transfer rate in terms of the local skin-friction coefficient and the local Nusselt number for a wide range of magnetohydrodynamic parameter, viscosity-variation parameter and viscous dissipation parameter. MHD flow in this geometry with temperature dependent viscosity is absent in the literature. The results obtained from the numerical simulations have been veri-fied by two methodologies.

Analysis of conical-cylinder shells fundamental characteristics for structural health diagnostics

This paper describes the formulation for the free vibration of joined conical-cylindrical shells with uniform thickness using the transfer of influence coefficient for identification of structural characteristics. These characteristics are importance for structural health monitoring to develop model. This method was developed based on successive transmission of dynamic influence coefficients, which were defined as the relationships between the displacement and the force vectors at arbitrary nodal circles of the system. The two edges of the shell having arbitrary boundary conditions are supported by several elastic springs with meridional/axial, circumferential, radial and rotational stiffness, respectively. The governing equations of vibration of a conical shell, including a cylindrical shell, are written as a coupled set of first order differential equations by using the transfer matrix of the shell. Once the transfer matrix of a single component has been determined, the entire structure matrix is obtained by the product of each component matrix and the joining matrix. The natural frequencies and the modes of vibration were calculated numerically for joined conical-cylindrical shells. The validity of the present method is demonstrated through simple numerical examples, and through comparison with the results of previous researchers.

Corneal changes following near work in myopic anisometropia

PURPOSE: To examine the symmetry of corneal changes following near work in the fellow eyes of non-amblyopic myopic anisometropes. METHODS: Thirty-four non-amblyopic myopic anisometropes (minimum 1 D spherical equivalent anisometropia) were recruited. Corneal topography was measured with the Medmont E300 Videokeratoscope before and after a controlled near task. Subjects were positioned to minimise head movements and read continuous text on a computer monitor for 10 minutes at an angle of 25 degrees downward gaze and an accommodation demand of 2.5 D. Measures of palpebral aperture morphology during primary and downward gaze were also obtained using digital photography and analysed with customised software. RESULTS: Significant changes in corneal topography were observed after ten minutes of reading. Localised superior regions of corneal topographical change (a hyperopic shift in corneal power) were typically exhibited in both eyes following the near task. The mean change in the corneal sphero-cylinder was +0.02/-0.11 x 113 and +0.02/-0.06 x 68 for the more and less myopic eyes respectively for a 6 mm corneal diameter. A significantly greater change in corneal astigmatism power vector J0 (a larger increase in against the rule astigmatism) was observed in the more myopic eyes (p < 0.01 for a 6 mm diameter). The more and less myopic eyes exhibited a high degree of interocular symmetry for measures of palpebral aperture morphology during both primary and downward gaze. Changes in corneal power vectors following reading were associated with eyelid position during downward gaze. CONCLUSIONS: Changes in corneal topography observed following a controlled reading task were highly symmetrical between the fellow eyes of myopic anisometropes due to the interocular symmetry of the palpebral aperture. However, the more myopic eye did exhibit a small but significantly greater magnitude of change in corneal astigmatism compared to the less myopic eye following reading. These findings may have implications for understanding the mechanism of development of non-amblyopic myopic anisometropia.

Shear thinning and shear thickening non-Newtonian confined fluid flow over rotating cylinder

In present work, numerical solution is performed to study the confined flow of power-law non Newtonian fluids over a rotating cylinder. The main purpose is to evaluate drag and thermal coefficients as functions of the related governing dimensionless parameters, namely, power-law index (0.5 ≤ n ≤ 1.4), dimensionless rotational velocity (0 ≤ α ≤ 6) and the Reynolds number (100 ≤ Re ≤ 500). Over the range of Reynolds number, the flow is known to be steady. Results denoted that the increment of power law index and rotational velocity increases the drag coefficient due to momentum diffusivity improvement which is responsible for low rate of heat transfer, because the thicker the boundary layer, the lower the heat transfer is implemented.

Non-Newtonian natural convection flow along an isothermal horizontal circular cylinder using modified power-law model

Laminar two-dimensional natural convection boundary-layer flow of non-Newtonian fluids along an isothermal horizontal circular cylinder has been studied using a modified power-law viscosity model. In this model, there are no unrealistic limits of zero or infinite viscosity. Therefore, the boundary-layer equations can be solved numerically by using marching order implicit finite difference method with double sweep technique. Numerical results are presented for the case of shear-thinning as well as shear thickening fluids in terms of the fluid velocity and temperature distributions, shear stresses and rate of heat transfer in terms of the local skin-friction and local Nusselt number respectively.

Gravity-driven fingering simulations for a thin liquid film flowing down the outside of a vertical cylinder

A numerical study is presented to examine the fingering instability of a gravity-driven thin liquid film flowing down the outer wall of a vertical cylinder. The lubrication approximation is employed to derive an evolution equation for the height of the film, which is dependent on a single parameter, the dimensionless cylinder radius. This equation is identified as a special case of that which describes thin film flow down an inclined plane. Fully three-dimensional simulations of the film depict a fingering pattern at the advancing contact line. We find the number of fingers observed in our simulations to be in excellent agreement with experimental observations and a linear stability analysis reported recently by Smolka & SeGall (Phys Fluids 23, 092103 (2011)). As the radius of the cylinder decreases, the modes of perturbation have an increased growth rate, thus increasing cylinder curvature partially acts to encourage the contact line instability. In direct competition with this behaviour, a decrease in cylinder radius means that fewer fingers are able to form around the circumference of the cylinder. Indeed, for a sufficiently small radius, a transition is observed, at which point the contact line is stable to transverse perturbations of all wavenumbers. In this regime, free surface instabilities lead to the development of wave patterns in the axial direction, and the flow features become perfectly analogous to the two-dimensional flow of a thin film down an inverted plane as studied by Lin & Kondic (Phys Fluids 22, 052105 (2010)). Finally, we simulate the flow of a single drop down the outside of the cylinder. Our results show that for drops with low volume, the cylinder curvature has the effect of increasing drop speed and hence promoting the phenomenon of pearling. On the other hand, drops with much larger volume evolve to form single long rivulets with a similar shape to a finger formed in the aforementioned simulations.

Numerical solutions for thin film flow down the outside and inside of a vertical cylinder

We consider a model for thin film flow down the outside and inside of a vertical cylinder. Our focus is to study the effect that the curvature of the cylinder has on the gravity-driven instability of the advancing contact line and to simulate the resulting fingering patterns that form due to this instability. The governing partial differential equation is fourth order with a nonlinear degenerate diffusion term that represents the stabilising effect of surface tension. We present numerical solutions obtained by implementing an efficient alternating direction implicit scheme. When compared to the problem of flow down a vertical plane, we find that increasing substrate curvature tends to increase the fingering instability for flow down the outside of the cylinder, whereas flow down the inside of the cylinder substrate curvature has the opposite effect. Further, we demonstrate the existence of nontrivial travelling wave solutions which describe fingering patterns that propagate down the inside of a cylinder at constant speed without changing form. These solutions are perfectly analogous to those found previously for thin film flow down an inclined plane.

Raman, infrared and near-infrared spectroscopic characterization of the herderite-hydroxylherderite mineral series

Natural single-crystal specimens of the herderite-hydroxylherderite series from Brazil, with general formula CaBePO4(F,OH), were investigated by electron microprobe, Raman, infrared and near-infrared spectroscopies. The minerals occur as secondary products in granitic pegmatites. Herderite and hydroxylherderite minerals show extensive solid solution formation. The Raman spectra of hydroxylherderite are characterized by bands at around 985 and 998 cm-1, assigned to ν1 symmetric stretching mode of the HOPO33- and PO43- units. Raman bands at around 1085, 1128 and 1138 cm-1 are attributed to both the HOP and PO antisymmetric stretching vibrations. The set of Raman bands observed at 563, 568, 577, 598, 616 and 633 cm-1 are assigned to the ν4 out of plane bending modes of the PO4 and H2PO4 units. The OH Raman stretching vibrations of hydroxylherderite were observed ranging from 3626 cm-1 to 3609 cm-1. The infrared stretching vibrations of hydroxylherderites were observed between 3606 cm-1 and 3599 cm-1. By using a Libowitzky type function, hydrogen bond distances based upon the OH stretching bands were calculated. Characteristic NIR bands at around 6961 and 7054 cm-1 were assigned to the first overtone of the fundamental, whilst NIR bands at 10194 and 10329 cm-1 are assigned to the second overtone of the fundamental OH stretching vibration. Insight into the structure of the herderite-hydroxylherderite series is assessed by vibrational spectroscopy.