Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Ecophysiologie de la maturation sexuelle et de la ponte de l'huître creuse Crassostrea gigas : réponses métaboliques (respiration) et alimentaires (filtration, absorption) en fonction des différents stades de maturation

The study of maturation and spawning of the oyster is part of a research program to investigate the summer mortalities of the oysters, Crassostrea gigas in Marennes-Oléron Bay. Four maturity stages were simultaneously obtained by diet and thermal conditioning (immature, low maturation, mature and post-spawning stages). Measurements of clearance, filtration, absorption and respiration rates allowed a calculation of the scope for growth and hence an estimation of the oyster's energetic budget at various maturity stages. Male and female oysters had similar physiological responses. The filtration rate ranged from 2.4 to 2.6 1.h(-1) at the early stages of maturation and decreased to 1.8 1.h.' during the maturity stage. Growth rate resulting from gonad development did not induce filtration rate changes. Mature 2.5 and 1.5-year-old oysters showed a negative energy budget reaching -15 and -90 J.h(-1) respectively. By contrast, non-ripe oysters had scope for growth in the range 110 to 170 J.h(-1). A negative energy budget during the high maturation stage resulted from a reduced absorption efficiency. A new allometric relationship for the respiration model of C. gigas was defined during vitellogenesis with a 0.574 coefficient value. Based on Our results, the oyster's physiological weakness during vitellogenesis should be considered as a part of explanation for spring and summer mortalities of cultured oysters in Marennes-Oléron Bay.

Ecophysiologie de la maturation sexuelle et de la ponte de l'huître creuse Crassostrea gigas : réponses métaboliques (respiration) et alimentaires (filtration, absorption) en fonction des différents stades de maturation

The study of sexual maturation and spawning in the Pacific oyster (Crassostrea gigas) is part of a vast research programme that endeavours to understand the causes of mortality that occur sporadically during the spring and summer seasons in the Marennes-Oléron Bay. Thermal and diet conditioning were used to obtain oysters at each stage of maturity simultaneously. Using the measured rates of clearance, consumption, absorption and respiration provided estimates of growth potential and gave the energetic budget of oysters at different stages of sexual maturity. Physiological responses were similar for males and females. Filtration decreased from 2.4 to 2.6 l.h (-1) to 1.8 l.h (-1) with increasing maturity. Weight gain was associated with gonad development and did not appear to have an effect on the clearance rate. Oysters 2.5 years old showed a negative energy budget (-15 J h (-1)) at later maturity stages. This deficit was confirmed (90 J.h (-1)) in oysters 1.5 years old at the same stage of maturity. On the contrary, immature oysters, in the early stages of maturity or post-spawning, had a growth potential of 110 to 170 J.h (-1). The energy deficit observed at later stages of maturity was primarily due to absorption, which decreased sharply during peak gametogenesis. Using measured respiration rates, an allometric relationship specific to gonad growth was determined with a coefficient of 0.574. Low physiological performance of oysters, observed at later stages of sexual maturity, must be taken into account in research on the factors responsible for spring and summer mortalities affecting oyster farms in Marennes-Oléron.

Thermal regime of the Southeast Indian Ridge between 88°E and 140°E: Remarks on the subsidence of the ridge flanks

The flanks of the Southeast Indian Ridge are characterized by anomalously low subsidence rates for the 0-25 Ma period: less than 300 m Ma(-1/2) between 101 degrees E and 120 degrees E and less than 260 m Ma(-1/2) within the Australian-Antarctic Discordance (AAD), between 120 degrees E and 128 degrees E. The expected along-axis variation in mantle temperature (similar to 50 degrees C) is too small to explain this observation, even when the temperature dependence of the mantle physical properties is accounted for. We successively analyze the effect on subsidence of different factors, such as variations in crustal thickness; the dynamic contribution of an old, detached slab supposedly present within the mantle below the AAD; and depletion in phi(m), a parameter here defined as the "ubiquitously distributed melt fraction" within the asthenosphere. These effects may all contribute to the observed, anomalously low subsidence rate of the ridge flanks, with the most significant contribution being probably related to the depletion in phi(m). However, these effects have a deep-seated origin that cannot explain the abruptness of the transition across the fracture zones that delineate the boundaries of the AAD, near 120 degrees E and near 128 degrees E, respectively.

Thermal Tolerance of the Coffee Berry Borer Hypothenemus hampei: Predictions of Climate Change Impact on a Tropical Insect Pest

Coffee is predicted to be severely affected by climate change. We determined the thermal tolerance of the coffee berry borer, Hypothenemus hampei, the most devastating pest of coffee worldwide, and make inferences on the possible effects of climate change using climatic data from Colombia, Kenya, Tanzania, and Ethiopia. For this, the effect of eight temperature regimes (15, 20, 23, 25, 27, 30, 33 and 35 degrees C) on the bionomics of H. hampei was studied. Successful egg to adult development occurred between 20-30 degrees C. Using linear regression and a modified Logan model, the lower and upper thresholds for development were estimated at 14.9 and 32 degrees C, respectively. In Kenya and Colombia, the number of pest generations per year was considerably and positively correlated with the warming tolerance. Analysing 32 years of climatic data from Jimma (Ethiopia) revealed that before 1984 it was too cold for H. hampei to complete even one generation per year, but thereafter, because of rising temperatures in the area, 1-2 generations per year/coffee season could be completed. Calculated data on warming tolerance and thermal safety margins of H. hampei for the three East African locations showed considerably high variability compared to the Colombian site. The model indicates that for every 1 degrees C rise in thermal optimum (T(opt)), the maximum intrinsic rate of increase (r(max)) will increase by an average of 8.5%. The effects of climate change on the further range of H. hampei distribution and possible adaption strategies are discussed. Abstracts in Spanish and French are provided as supplementary material Abstract S1 and Abstract S2.

FLOW REGIME DRIVEN THERMAL ENHANCEMENT IN INTERNALLY-GROOVED TUBES

Internally-grooved refrigeration tubes maximize tube-side evaporative heat transfer rates and have been identified as a most promising technology for integration into compact cold plates. Unfortunately, the absence of phenomenological insights and physical models hinders the extrapolation of grooved-tube performance to new applications. The success of regime-based heat transfer correlations for smooth tubes has motivated the current effort to explore the relationship between flow regimes and enhanced heat transfer in internally-grooved tubes. In this thesis, a detailed analysis of smooth and internally-grooved tube data reveals that performance improvement in internally-grooved tubes at low-to-intermediate mass flux is a result of early flow regime transition. Based on this analysis, a new flow regime map and corresponding heat transfer coefficient correlation, which account for the increased wetted angle, turbulence, and Gregorig effects unique to internally-grooved tubes, were developed. A two-phase test facility was designed and fabricated to validate the newly-developed flow regime map and regime-based heat transfer coefficient correlation. As part of this setup, a non-intrusive optical technique was developed to study the dynamic nature of two-phase flows. It was found that different flow regimes result in unique temporally varying film thickness profiles. Using these profiles, quantitative flow regime identification measures were developed, including the ability to explain and quantify the more subtle transitions that exist between dominant flow regimes. Flow regime data, based on the newly-developed method, and heat transfer coefficient data, using infrared thermography, were collected for two-phase HFE-7100 flow in horizontal 2.62mm - 8.84mm diameter smooth and internally-grooved tubes with mass fluxes from 25-300 kg/m²s, heat fluxes from 4-56 kW/m², and vapor qualities approaching 1. In total, over 6500 combined data points for the adiabatic and diabatic smooth and internally-grooved tubes were acquired. Based on results from the experiments and a reinterpretation of data from independent researchers, it was established that heat transfer enhancement in internally-grooved tubes at low-to-intermediate mass flux is primarily due to early flow regime transition to Annular flow. The regime-based heat transfer coefficient outperformed empirical correlations from the literature, with mean and absolute deviations of 4.0% and 32% for the full range of data collected.

COMPREHENSIVE ELECTRICAL/OPTICAL/THERMAL CHARACTERIZATIONS OF HIGH POWER LIGHT EMITTING DIODES AND LASER DIODES

Thermal characterizations of high power light emitting diodes (LEDs) and laser diodes (LDs) are one of the most critical issues to achieve optimal performance such as center wavelength, spectrum, power efficiency, and reliability. Unique electrical/optical/thermal characterizations are proposed to analyze the complex thermal issues of high power LEDs and LDs. First, an advanced inverse approach, based on the transient junction temperature behavior, is proposed and implemented to quantify the resistance of the die-attach thermal interface (DTI) in high power LEDs. A hybrid analytical/numerical model is utilized to determine an approximate transient junction temperature behavior, which is governed predominantly by the resistance of the DTI. Then, an accurate value of the resistance of the DTI is determined inversely from the experimental data over the predetermined transient time domain using numerical modeling. Secondly, the effect of junction temperature on heat dissipation of high power LEDs is investigated. The theoretical aspect of junction temperature dependency of two major parameters – the forward voltage and the radiant flux – on heat dissipation is reviewed. Actual measurements of the heat dissipation over a wide range of junction temperatures are followed to quantify the effect of the parameters using commercially available LEDs. An empirical model of heat dissipation is proposed for applications in practice. Finally, a hybrid experimental/numerical method is proposed to predict the junction temperature distribution of a high power LD bar. A commercial water-cooled LD bar is used to present the proposed method. A unique experimental setup is developed and implemented to measure the average junction temperatures of the LD bar. After measuring the heat dissipation of the LD bar, the effective heat transfer coefficient of the cooling system is determined inversely. The characterized properties are used to predict the junction temperature distribution over the LD bar under high operating currents. The results are presented in conjunction with the wall-plug efficiency and the center wavelength shift.

Effect of cooling rate on properties of plasma nitrided AISI 1010 steel

In this work, AISI 1010 steel samples were plasma nitrided into 20% N 2 100 Pa and 400 Pa for N 2 and H 2 , respectively), temperatures of 500 and 580 °C, during 2 h. Three different procedures for cooling were accomplished after nitriding. In the first procedure the cooling occurred naturally, that is, the sample was kept on substrate holder. In the second one the sample was pulled off and cooling in a cold surface. Finally, in the third cooling process the sample was pulled off the substrate holder down into special reservoir filled with oil held at ambient temperature. The properties of the AISI 1010 steel samples were characterized by optical and electron microscopy, X-ray diffraction, Mössbauer spectroscopy and microhardness tests. Thermal gradient inside the sample kept on substrate holder during cooling process was measured by three inserted thermocouples at different depths. When samples were cooled rapidly the transformation of ϵ-Fe 2 − 3 N to γ′-Fe 4 N was inhibited. Such effect is indicated by the high concentration of ϵ-Fe compound zone. To get solid state solution of nitrogen in the diffusion zone, instead of precipitates of nitride phases, the cooling rate should be higher than a critical value of about 0.95 °C/s. When this value is reached at any depth of the diffusion zone, two distinct diffusion zones will appear. Temperature gradients were measured inside the samples as a consequence of the plasma treatment. It's suggested the need for standardization of the term “treatment temperature” for plasma treatment because different nitrided layer properties could be reported for the same “treatment temperature”.

Self-sealing of thermal fatigue and mechanical damage in fiber-reinforced composite materials

Fiber reinforced composite tanks provide a promising method of storage for liquid oxygen and hydrogen for aerospace applications. The inherent thermal fatigue of these vessels leads to the formation of microcracks, which allow gas phase leakage across the tank walls. In this dissertation, self-healing functionality is imparted to a structural composite to effectively seal microcracks induced by both mechanical and thermal loading cycles. Two different microencapsulated healing chemistries are investigated in woven glass fiber/epoxy and uni-weave carbon fiber/epoxy composites. Self-healing of mechanically induced damage was first studied in a room temperature cured plain weave E-glass/epoxy composite with encapsulated dicyclopentadiene (DCPD) monomer and wax protected Grubbs' catalyst healing components. A controlled amount of microcracking was introduced through cyclic indentation of opposing surfaces of the composite. The resulting damage zone was proportional to the indentation load. Healing was assessed through the use of a pressure cell apparatus to detect nitrogen flow through the thickness direction of the damaged composite. Successful healing resulted in a perfect seal, with no measurable gas flow. The effect of DCPD microcapsule size (51 um and 18 um) and concentration (0 - 12.2 wt%) on the self-sealing ability was investigated. Composite specimens with 6.5 wt% 51 um capsules sealed 67% of the time, compared to 13% for the control panels without healing components. A thermally stable, dual microcapsule healing chemistry comprised of silanol terminated poly(dimethyl siloxane) plus a crosslinking agent and a tin catalyst was employed to allow higher composite processing temperatures. The microcapsules were incorporated into a satin weave E-glass fiber/epoxy composite processed at 120C to yield a glass transition temperature of 127C. Self-sealing ability after mechanical damage was assessed for different microcapsule sizes (25 um and 42 um) and concentrations (0 - 11 vol%). Incorporating 9 vol% 42 um capsules or 11 vol% 25 um capsules into the composite matrix leads to 100% of the samples sealing. The effect of microcapsule concentration on the short beam strength, storage modulus, and glass transition temperature of the composite specimens was also investigated. The thermally stable tin catalyzed poly(dimethyl siloxane) healing chemistry was then integrated into a [0/90]s uniweave carbon fiber/epoxy composite. Thermal cycling (-196C to 35C) of these specimens lead to the formation of microcracks, over time, formed a percolating crack network from one side of the composite to the other, resulting in a gas permeable specimen. Crack damage accumulation and sample permeability was monitored with number of cycles for both self-healing and traditional non-healing composites. Crack accumulation occurred at a similar rate for all sample types tested. A 63% increase in lifetime extension was achieved for the self-healing specimens over traditional non-healing composites.

Effect of different feed-rate in bone drilling: experimental and numerical study

The behaviour of bone tissue during drilling has been subject of recent studies due to its great importance. Because of thermal nature of the bone drilling, high temperatures and thermal mechanical stresses are developed during drilling that affect the process quality. However, there is still a lack information with regard to the distribution of mechanical and thermal stresses during bone drilling. The present paper describes a sequentially coupled thermal-stress analysis to assess the mechanical and thermal stress distribution during bone drilling. A three-dimensional thermo-mechanical model was developed using the ANSYS/LSDYNA finite element code under different drilling conditions. The model incorporates the dynamic characteristics of drilling process, as well as the thermo-mechanical properties of the involved materials. Experimental tests with polyurethane foam materials were also carried out. It was concluded that the use of higher feed-rates lead to a decrease of normal stresses and strains in the foam materials. The experimental and numerical results were compared and showed good agreement. The proposed numerical model could be used to predict the better drilling parameters and minimize the bone injuries.

Effect of different feed-rate in bone drilling: experimental and numerical study

The behaviour of bone tissue during drilling has been subject of recent studies due to its great importance. Because of thermal nature of the bone drilling, high temperatures and thermal mechanical stresses are developed during drilling that affect the process quality. However, there is still a lack information with regard to the distribution of mechanical and thermal stresses during bone drilling.

Thermal analysis in the drilling of solid rigid foam materials and ex-vivo bovine bones

Bone is a dynamic, highly vascularized tissue with a unique capacity to heal and regenerate without scarring. However, drilling remains a concern in several clinical procedures due to thermal damage of the bone and surrounding tissue. The success of this surgeries is dependent of many factors and also in temperature generation during the drilling bone. When an excessive heat is produced during the drilling, thermal necrosis can occur and the bone suffers injuries. Studies have shown that the increased temperature is directly related with the drilling parameters, particularly, the drill speed, feed-rate, applied force, the depth of cut, the geometry of the drill bit, the use or not of a cooling system and also the type of bone.

Modeling and analysis of thermoelectric energy conversion efficiency in nanostructures.

116 p.

The Processing of Aluminum Gasarites via Thermal Decomposition of Interstitial Hydrides

Gasarite structures are a unique type of metallic foam containing tubular pores. The original methods for their production limited them to laboratory study despite appealing foam properties. Thermal decomposition processing of gasarites holds the potential to increase the application of gasarite foams in engineering design by removing several barriers to their industrial scale production. The following study characterized thermal decomposition gasarite processing both experimentally and theoretically. It was found that significant variation was inherent to this process therefore several modifications were necessary to produce gasarites using this method. Conventional means to increase porosity and enhance pore morphology were studied. Pore morphology was determined to be more easily replicated if pores were stabilized by alumina additions and powders were dispersed evenly. In order to better characterize processing, high temperature and high ramp rate thermal decomposition data were gathered. It was found that the high ramp rate thermal decomposition behavior of several hydrides was more rapid than hydride kinetics at low ramp rates. This data was then used to estimate the contribution of several pore formation mechanisms to the development of pore structure. It was found that gas-metal eutectic growth can only be a viable pore formation mode if non-equilibrium conditions persist. Bubble capture cannot be a dominant pore growth mode due to high bubble terminal velocities. Direct gas evolution appears to be the most likely pore formation mode due to high gas evolution rate from the decomposing particulate and microstructural pore growth trends. The overall process was evaluated for its economic viability. It was found that thermal decomposition has potential for industrialization, but further refinements are necessary in order for the process to be viable.

Development of measurement and modeling techniques to quantify atmospheric deposition of persistent, bioaccumulative and toxic chemicals in the Great Lakes

Measurement and modeling techniques were developed to improve over-water gaseous air-water exchange measurements for persistent bioaccumulative and toxic chemicals (PBTs). Analytical methods were applied to atmospheric measurements of hexachlorobenzene (HCB), polychlorinated biphenyls (PCBs), and polybrominated diphenyl ethers (PBDEs). Additionally, the sampling and analytical methods are well suited to study semivolatile organic compounds (SOCs) in air with applications related to secondary organic aerosol formation, urban, and indoor air quality. A novel gas-phase cleanup method is described for use with thermal desorption methods for analysis of atmospheric SOCs using multicapillary denuders. The cleanup selectively removed hydrogen-bonding chemicals from samples, including much of the background matrix of oxidized organic compounds in ambient air, and thereby improved precision and method detection limits for nonpolar analytes. A model is presented that predicts gas collection efficiency and particle collection artifact for SOCs in multicapillary denuders using polydimethylsiloxane (PDMS) sorbent. An approach is presented to estimate the equilibrium PDMS-gas partition coefficient (Kpdms) from an Abraham solvation parameter model for any SOC. A high flow rate (300 L min-1) multicapillary denuder was designed for measurement of trace atmospheric SOCs. Overall method precision and detection limits were determined using field duplicates and compared to the conventional high-volume sampler method. The high-flow denuder is an alternative to high-volume or passive samplers when separation of gas and particle-associated SOCs upstream of a filter and short sample collection time are advantageous. A Lagrangian internal boundary layer transport exchange (IBLTE) Model is described. The model predicts the near-surface variation in several quantities with fetch in coastal, offshore flow: 1) modification in potential temperature and gas mixing ratio, 2) surface fluxes of sensible heat, water vapor, and trace gases using the NOAA COARE Bulk Algorithm and Gas Transfer Model, 3) vertical gradients in potential temperature and mixing ratio. The model was applied to interpret micrometeorological measurements of air-water exchange flux of HCB and several PCB congeners in Lake Superior. The IBLTE Model can be applied to any scalar, including water vapor, carbon dioxide, dimethyl sulfide, and other scalar quantities of interest with respect to hydrology, climate, and ecosystem science.