A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Non-Linear Scale Interactions in a Forced Turbulent Boundary Layer

This thesis explores the dynamics of scale interactions in a turbulent boundary layer through a forcing-response type experimental study. An emphasis is placed on the analysis of triadic wavenumber interactions since the governing Navier-Stokes equations for the flow necessitate a direct coupling between triadically consist scales. Two sets of experiments were performed in which deterministic disturbances were introduced into the flow using a spatially-impulsive dynamic wall perturbation. Hotwire anemometry was employed to measure the downstream turbulent velocity and study the flow response to the external forcing. In the first set of experiments, which were based on a recent investigation of dynamic forcing effects in a turbulent boundary layer, a 2D (spanwise constant) spatio-temporal normal mode was excited in the flow; the streamwise length and time scales of the synthetic mode roughly correspond to the very-large-scale-motions (VLSM) found naturally in canonical flows. Correlation studies between the large- and small-scale velocity signals reveal an alteration of the natural phase relations between scales by the synthetic mode. In particular, a strong phase-locking or organizing effect is seen on directly coupled small-scales through triadic interactions. Having characterized the bulk influence of a single energetic mode on the flow dynamics, a second set of experiments aimed at isolating specific triadic interactions was performed. Two distinct 2D large-scale normal modes were excited in the flow, and the response at the corresponding sum and difference wavenumbers was isolated from the turbulent signals. Results from this experiment serve as an unique demonstration of direct non-linear interactions in a fully turbulent wall-bounded flow, and allow for examination of phase relationships involving specific interacting scales. A direct connection is also made to the Navier-Stokes resolvent operator framework developed in recent literature. Results and analysis from the present work offer insights into the dynamical structure of wall turbulence, and have interesting implications for design of practical turbulence manipulation or control strategies.

Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

A New Criterion Number for the Boundary Conditions at the Solid/Liquid Interface in Nanoscale

A simple, but important three-atom model was proposed at the solid/liquid interface, leading to a new criterion number, lambda, governing the boundary conditions (BCs) in nanoscale. The solid wall is considered as the face-centered-cubic (fcc) structure. The fluid is the liquid argon with the well-known LJ potential. Based on the concept, the two micro-systems have the same BCs if they have The same criterion number. The degree of the locking BCs is enhanced when lambda equals to 0.757. Such critical criterion number results in the substantial epitaxial ordering and one, two, or even three liquid layers are locked by the solid wall, depending on the coupling energy scale ratio of the solid and liquid atoms. With deviation from the critical criterion number, the flow approaches the slip BCs and there are little ordering structures within the liquid. Always at the same criterion number, the degree of the slip is decreased or the locking is enhanced with increasing the coupling energy scale ratio of the solid and liquid atoms. The above analysis is well confirmed by the molecular dynamics (MD) simulation. The slip length is well correlated in terms of the new criterion number. The future work is suggested to extend the present theory for other microstructures of the solid wall atoms and quasi-LJ potentials.

Numerical approximation of a fractional-in-space diffusion equation (II) - with nonhomogeneous boundary conditions

In this paper, a space fractional di®usion equation (SFDE) with non- homogeneous boundary conditions on a bounded domain is considered. A new matrix transfer technique (MTT) for solving the SFDE is proposed. The method is based on a matrix representation of the fractional-in-space operator and the novelty of this approach is that a standard discretisation of the operator leads to a system of linear ODEs with the matrix raised to the same fractional power. Analytic solutions of the SFDE are derived. Finally, some numerical results are given to demonstrate that the MTT is a computationally e±cient and accurate method for solving SFDE.

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.

Nanoparticle melting as a Stefan moving boundary problem

The melting of spherical nanoparticles is considered from the perspective of heat flow in a pure material and as a moving boundary (Stefan) problem. The dependence of the melting temperature on both the size of the particle and the interfacial tension is described by the Gibbs-Thomson effect, and the resulting two-phase model is solved numerically using a front-fixing method. Results show that interfacial tension increases the speed of the melting process, and furthermore, the temperature distribution within the solid core of the particle exhibits behaviour that is qualitatively different to that predicted by the classical models without interfacial tension.

The moderating role of self-efficacy on the relationship between emotional exhaustion and job strain for boundary-spanning service employees

Appropriate behaviours toward customers often requires employees to suppress some genuine emotions and/or express other emotions; genuine or contrived. Managing one's emotions in this way gives rise to emotional exhaustion. This can have consequences for psychological ill health, in the form of work place strain, and ultimately employee's desire to leave. This student examines the relationships between emotional management, emotional exhaustion and turnover intentions amongst diversional therapy professionals. We find that some forms of emotional management have a significant impact on emotional exhaustion and that this predicts workplace strain. Furthermore, the deleterious effects of emotional exhaustion are mitigated somewhat for employees who have strong beliefs in their ability to provide good service, compared to employees with lower self efficacy beliefs.