An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Multi-scale micro-structure generation strategy for up-scaling transport in clays

Clays and claystones are used as backfill and barrier materials in the design of waste repositories, because they act as hydraulic barriers and retain contaminants. Transport through such barriers occurs mainly by molecular diffusion. There is thus an interest to relate the diffusion properties of clays to their structural properties. In previous work, we have developed a concept for up-scaling pore-scale molecular diffusion coefficients using a grid-based model for the sample pore structure. Here we present an operational algorithm which can generate such model pore structures of polymineral materials. The obtained pore maps match the rock’s mineralogical components and its macroscopic properties such as porosity, grain and pore size distributions. Representative ensembles of grains in 2D or 3D are created by a lattice Monte Carlo (MC) method, which minimizes the interfacial energy of grains starting from an initial grain distribution. Pores are generated at grain boundaries and/or within grains. The method is general and allows to generate anisotropic structures with grains of approximately predetermined shapes, or with mixtures of different grain types. A specific focus of this study was on the simulation of clay-like materials. The generated clay pore maps were then used to derive upscaled effective diffusion coefficients for non-sorbing tracers using a homogenization technique. The large number of generated maps allowed to check the relations between micro-structural features of clays and their effective transport parameters, as is required to explain and extrapolate experimental diffusion results. As examples, we present a set of 2D and 3D simulations and investigated the effects of nanopores within particles (interlayer pores) and micropores between particles. Archie’s simple power law is followed in systems with only micropores. When nanopores are present, additional parameters are required; the data reveal that effective diffusion coefficients could be described by a sum of two power functions, related to the micro- and nanoporosity. We further used the model to investigate the relationships between particle orientation and effective transport properties of the sample.

Dielectrophoresis-based analyte separation and analysis

Dielectrophoresis—the tendency of a material of high dielectric permittivity to migrate in an electrical field gradient to a region of maximum field strength—provides an ideal motive force for manipulating small volumes of biological analytes in microfluidic microsystems. The work described in this thesis was based on the hypothesis that dielectrophoresis could be exploited to provide high-resolution cell separations in microsystems as well as a means for the electrically-controllable manipulation of solid supports for molecular analysis. To this end, a dielectrophoretic/gravitational field-flow-fractionation (DEP/G-FFF) system was developed and the separation performance evaluated using various types and sizes of polystyrene microspheres as model particles. It was shown that separation of the polystyrene beads was based on the differences in their effective dielectrophoretic properties. The ability of an improved DEP/G-FFF system to separate genetically identical, but phenotypically dissimilar cell types was demonstrated using mixtures of 6m2 mutant rat kidney cells grown under transforming and non-transforming culture conditions. Additionally, a panel of engineered dielectric microspheres was designed with specific, predetermined dielectrophoretic properties such that their dielectrophoretic behaviors would be controllable and predictable. The fabrication method involved the use of gold-coated polystyrene microsphere cores coated with a self-assembled monolayer of alkanethiol and, optionally, a self-assembled monolayer of phospholipid to form a thin-insulating-shell-over-conductive-interior structure. The successful development of the DEP/G-FFF separation system and the dielectrically engineered microspheres provides proof-of-principle demonstrations of enabling dielectrophoresis-based microsystem technology that should provide powerful new methods for the manipulation, separation and identification of analytes in many diverse fields. ^

Mejora del comportamiento térmico de los morteros de cal aditivados y su empleo en la rehabilitación de inmuebles = Improvement of thermal behaviour of additived lime mortars and their use in rehabilitation of buildings

La envolvente de la edificación es la responsable de equilibrar el intercambio energético entre el interior y el exterior, por lo tanto cualquier actuación encaminada a la reducción del consumo energético ha de establecer, como uno de sus objetivos prioritarios, la mejora del comportamiento de la misma. Las edificaciones anteriores a 1940 constituyen la mayor parte de las existentes en áreas rurales y centros urbanos. En ellas, la repercusión de la fachada sobre las transmitancias globales pone de manifiesto la necesidad de intervención. Sin embargo, su elevada inercia térmica y los importantes saltos térmicos característicos de gran parte de España plantean la importancia de que aquélla se efectúe por el exterior. A tal respecto, la falta de disponibilidad de espesor suficiente para implantar sistemas tipo SATE deriva en que, frecuentemente, la única solución viable sea la de aislar por el interior perdiendo con ello la capacidad de acumulación térmica del muro y con el asociado riesgo de condensaciones. La amplia tradición en el empleo de revestimientos, especialmente en base de cal, permiten que éstos sean utilizados no sólo como elemento estético o de protección de la obra de fábrica antigua sino también para la mejora del comportamiento térmico del soporte, si se aprovecha el mecanismo de transmisión térmica por radiación. Éste es el objetivo de la presente Tesis Doctoral en la que se estudia la modificación de las propiedades radiantes de los morteros de revestimiento para la mejora de la eficiencia energética de las construcciones históricas, principalmente las constituidas por muros monolíticos, aunque podría ser de aplicación a otro tipo de construcciones compuestas por diversas capas. Como punto de partida, se estudió y revisó la documentación disponible sobre las investigaciones de las tres áreas científico-tecnológicas que convergen en la Tesis Doctoral: rehabilitación, material y comportamiento térmico, a partir de lo cual se comprobó la inexistencia de estudios similares al objeto de la presente Tesis Doctoral. Complementariamente, se analizaron los revestimientos en lo concerniente a los materiales que los constituyen, la composición de las mezclas y características de cada una de las capas así como al enfoque que, desde el punto de vista térmico, se estimaba más adecuado para la obtención de los objetivos propuestos. Basándonos en dichos análisis, se preseleccionaron ochenta materiales que fueron ensayados en términos de reflectancia y emisividad para elegir cuatro con los que se continuó la investigación. Éstos, junto con la cal elegida para la investigación y el árido marmóreo característico de la última capa de revestimiento, fueron caracterizados térmicamente, de forma pormenorizada, así como química y físicamente. Los fundamentos teóricos y los estudios preliminares desarrollados con distintos materiales, en estado fresco y endurecido, fueron empleados en la dosificación de componentes de las mezclas, en dos proporciones distintas, para el estudio del efecto del agregado. Éstas se ensayaron en estado fresco, para comprobar su adecuación de puesta en obra y prever su VI adherencia al soporte, así como en estado endurecido a 28 y 90 días de curado, para conocer las propiedades que permitieran prever su compatibilidad con aquél y estimar el alcance de la reducción de transferencias térmicas lograda. Además, se estudiaron las características generales de las mezclas que sirvieron para establecer correlaciones entre distintas propiedades y entender los resultados mecánicos, físicos (comportamiento frente al agua) y energéticos. Del estudio conjunto de las distintas propiedades analizadas se propusieron dos mezclas, una blanca y otra coloreada, cuyas características permiten garantizar la compatibilidad con la obra de fábrica antigua, desde el punto de vista físico y mecánico, y preservar la autenticidad de los revestimientos, en cuanto a la técnica de aplicación tradicional en sistemas multicapa. El comportamiento térmico de las mismas, sobre una obra de fábrica de 40 cm de espesor, se estimó, en estado estacionario y pseudo-transitorio, demostrándose reducciones del flujo térmico entre 16-48%, en condiciones de verano, y entre el 6-11%, en invierno, en función de la coloración y de la rugosidad de la superficie, en comparación con el empleo de la mezcla tradicional. Por lo que, se constata la viabilidad de los materiales compuestos propuestos y su adecuación al objetivo de la investigación. VII ABSTRACT The envelope is responsible for balancing the energy exchange between the inside and outside in buildings. For this reason, any action aimed at reducing energy consumption must establish, as one of its key priorities, its improvement. In rural areas and urban centers, most of the constructions were built before 1940. In them, the impact of the façade on the global transmittance highlights the need for intervention. However, its high thermal inertia and fluctuation of temperatures in the majority of Spain bring up that it should be placed outside the insulation. In this regard, the lack of availability of enough thickness to implement systems such as ETICS results in that often the only viable solution is to isolate the interior, losing thereby the wall’s heat storage capacity with the associated risk of condensation. The tradition in the use of renders, especially lime-based, allows us to use them not only as an aesthetic element or to protect the ancient masonry, but also for improved thermal performance of the support by taking advantage of the heat transfer mechanism by radiation. This is the aim of this Doctoral Thesis in which the modification of the radiative properties of lime mortars for renders to improve the energy efficiency of historic buildings, mainly composed of monolithic walls, is studied, although it could be applied to other structures composed of several layers. As a starting point, the available literature in the three scientific-technological areas that converge at the Doctoral Thesis: rehabilitation, material and thermal behaviour, was reviewed, and confirmed the absence of researches similar to this Doctoral Thesis. Additionally, the renders were studied with regard to the materials that constitute them, the composition of the mixtures and the characteristics of each layer, as well as to the approach which, from a thermal point of view, was deemed the most suitable for achieving the objectives sets. Based on thre aforementioned analysis, eighty materials tested in terms of reflectance and emissivity were pre-selected, to choose four with which the research was continued. Common marble sand, used in the last layer of the renders, together with the appointed materials and hydrated lime were characterized thermally, in detail, as well as chemically and physically. The theoretical bases and preliminary studies with different materials, in fresh and hardened state, were used in the dosage of the composition of the mixtures. In order to study their effect they were used in two different proportions, that is, ten mixtures in total. These were tested in their fresh state to evaluate their setting-up suitability and foresee their adhesion to the support, as well as in their hardened state, at 28 and 90 days of curing, to establish the properties which enabled us to anticipate their compatibility with the old masonry walls and estimate the scope of the reduction of heat transfers achieved. In addition, the general characteristics of the mixtures used to establish correlations and to understand the mechanical, physical and energy results were studied. Two mixtures, one white and one colored, were proposed as the result of the different properties analysed, whose characteristics allow the guarantee of mechanical and physical compatibility VIII with the old masonry and preserve the authenticity of the renders. The thermal behavior of both, applied on a masonry wall 40 cm thick, was estimated at a steady and pseudo-transient state, with heat flow reductions between 16-48% during summertime and 6-11% during wintertime, depending on the color and surface roughness, compared to the use of the traditional mixture. So, the viability of the proposed composite materials and their fitness to the aim of the research are established.

Modelización de las propiedades radiativas e interacción de láseres ultraintensos con la materia

Desde el año 2004 el código ARWEN ha sido utilizado con éxito para simular y diseñar experimentos relacionados con blancos para fusión por confinamiento inercial [146], astrofísica de laboratorio [145], plasmas como amplificadores de láseres de rayos X [107] o plasmas creados por láser para la medición de espectros de transmisión. Para la realización de estas simulaciones es necesario, además de métodos de alto orden precisos y que presenten buenas propiedades conservativas, conocer ciertas propiedades de los plasmas. En el caso de la fluidodinámica y la conducción electrónica necesitaremos conocer la ecuación de estado [94, 49, 36], y para el transporte de la radiación será preciso disponer de los datos de absorción y emisión [104, 95, 40]. Hasta el año 2009 ARWEN dependía de códigos externos para la generación de estas tablas de opacidad, careciendo de control sobre los métodos empleados para su generación. Además estos códigos asumían equilibrio local termodinámico (LTE), limitando su validez a rangos de alta densidad y baja temperatura. En el marco de esta tesis se ha desarrollado el código BIGBART para la generación de tablas detalladas de opacidad y emisividad para su uso en el módulo de transporte de radiación. De esta forma el grupo dispondrá de su propia herramienta de generación de propiedades radiativas. El código desarrollado es capaz de tratar plasmas en estado fuera de equilibrio (non-LTE) mediante el modelo colisional-radiativo, extendiendo así el rango de validez de las tablas generadas. El trabajo desarrollado para implementar un código LTE/non-LTE estacionario es el siguiente Cálculo de estructura y datos atómicos. Se ha acoplado en código FAC a BIGBART, incorporando la capacidad para generar potenciales atómicos para una configuración y el cálculo de funciones de onda de electrones en orbitales ligados y libres. Aproximaciones y métodos para la obtención de tasas y secciones eficaces de procesos. Se han incluido y programado los modelos implementados en FAC para el cálculo de secciones eficaces de fotoionización, y tasas de decaimiento de emisión espontánea y autoionización. Además se ha incluido el modelo Plane-Wave Born (PWBA) para el cálculo de las secciones eficaces de ionización y excitación colisional. Modelos para la obtención de la distribución de estados iónicos dentro del plasma. Se ha programado un solver LTE basado en la ecuación de Saha-Boltzmann con efectos de ionización por presión debida a los iones adyacentes. También se ha implementado un modelo non-LTE colisionalradiativo para la resolución del sistema de ecuaciones que nos permite obtener la densidad de estados iónicos fuera de equilibrio. Modelo non-LTE RADIOM. Se ha implementado el modelo RADIOM para aproximar efectos de no-equilibrio mediante cálculos LTE a una temperatura equivalente, menor o igual que la temperatura electrónica real. Cálculo de las propiedades espectrales de absorción y emisión. Se han implementado los modelos para el cálculo de los perfiles espectrales de absorción y emisión para procesos entre niveles ligados, ligado-libre y librelibre. Aprovechando el trabajo realizado en este sentido, durante el transcurso de esta tesis se amplió el código BIGBART para tratar problemas con dependencia temporal. La extensión para tratar este tipo de problemas se orientó a la simulación numérica de la interacción de láseres ultra intensos en el rango XUV/rayos X. Para ello, además de adaptar el modelo non-LTE colisionalradiativo se incluyeron procesos adicionales asociados a la interacción de la materia con fotones altamente energéticos. También se han incluido modelos para el cálculo de las propiedades ópticas, y por ende las propiedades dieléctricas de la materia irradiada, de gran interés en algunas aplicaciones novedosas de estos láseres intensos. Debido a la naturaleza fuertemente fuera de equilibrio en la interacción de fotones de alta energía con la materia, se incluyó el tratamiento de la distribución de electrones libres fuera de equilibrio en la aproximación de Fokker-Planck, tanto para condiciones degeneradas como no degeneradas. El trabajo desarrollado en el código non-LTE con dependencia temporal es el siguiente Procesos asociados a láseres intensos XUV/rayos X. Se ha implementado el cálculo de procesos radiativos estimulados de absorción y emisión por el láser. También se han incluido procesos asociados a la creación de vacantes en capas internas electrónicas (Shake), además de doble autoionización y doble fotoionización. Cálculo de propiedades ópticas y dieléctricas en blancos sólidos. Se ha implementado un modelo para la absorción por bremsstrahlung inverso en blancos en estado sólido. Con el coeficiente de extinción debido a procesos de fotoabsorción resonante, fotoionización y bremsstrahlung inverso se obtiene el ´ındice de refracción mediante la relación de Kronig-Kramers. Electrones fuera de equilibrio. Se ha tratado la evolución de la distribución de electrones, cuando no está justificado asumir que es Maxwelliana o de Fermi-Dirac, mediante la aproximación de Fokker-Planck para la colisión entre electrones libres. En la resolución de la ecuación de Fokker-Planck se han incluido los procesos inelásticos por colisiones con iones y términos fuente por interacción con el láser y otros procesos. ABSTRACT Since 2004 the ARWEN code has been successfully used to simulate and design targets for inertial confinement fusion experiments [146], laboratory astrophysics [145], plasmas as X-ray lasers amplifiers [107] or laser created plasmas for measuring transmission spectra. To perform these simulations it is necessary, in addition to high order precise methods with good conservative properties, to know certain properties of plasmas. For fluid dynamic and electronic conduction we need to know the equation of state [94, 49, 36], and for radiation transport it will be necessary to have the data of the absorption and emission [104, 95, 40]. Until 2009 ARWEN depended on external codes to generate these opacity tables, lacking of control over the methods used for their generation. Besides, these codes assumed local thermodynamic equilibrium (LTE), limiting their validity ranges to high densities and low temperatures. As part of this thesis it has been developed the BIGBART code for generating detailed opacity and emissivity tables for use in the radiation transport module. This group will have its own tool for the generation of radiative properties. The developed code is capable of treating plasmas out of equilibrium (non-LTE) by means of a collisional-radiative model, extending the range of validity of the generated tables. The work to implement an LTE/non-LTE steady-state code is as follows Calculation of structure and atomic data. the FAC code was coupled to BIGBART, incorporating the ability to generate atomic potentials for calculating configuration wave functions for bound and free electrons. Approaches and methods for obtaining cross sections and processes rates. We have included and reprogrammed in Fortran the models implemented in FAC for calculation of photoionization cross sections and decay rates of spontaneous emission and autoionization. We also included the Plane- Wave Born (PWBA) model to calculate the cross sections of ionization and collisional excitation. Models for the obtention of the distribution of ionic states within the plasma. We programmed a LTE solver based on the Saha-Boltzmann equation with pressure ionization effects due to adjacent ions. It has also been implemented a non-LTE collisional-radiative model for solving the system of equations that allows us to obtain the density of ionic states out of equilibrium. Non-LTE RADIOM model. We have implemented the non-LTE RADIOM model to approximate non-equilibrium effects with LTE data at an equivalent temperature, lower or equal to the actual electronic temperature. Calculation of the spectral absorption and emission properties. Models have been implemented for the calculation of the spectral profiles of absorption and emission processes between bound levels, free-bound and free-free. Taking advantage of the work done in this direction throughout the course of this thesis the code BIGBART was extended to treat time-dependent problems. The extension to treat such problems is oriented to the numerical simulation of the interaction of ultra intense lasers in the XUV/X-ray range. For this range, in addition to adapting the non-LTE collisional-radiative model, additional processes associated with the interaction of matter with high energy photons. We also included models for calculation of the optical properties, and therefore the dielectric properties of the irradiated material, of great interest in some novel applications of these intense lasers. Due to the strong non-equilibrium nature of the interaction of high energy photons with matter, we included the treatment of the distribution of free electrons out of equilibrium in the Fokker-Planck approximation for both degenerate and non-degenerate conditions. The work in the non-LTE time-dependent code is as follows Processes associated with intense XUV/X-ray lasers. We have implemented the calculation of stimulated radiative processes in absorption and emission. Also we included processes associated with the creation of electronic vacancies in inner shells (Shake), double autoionization and double photoionization. Calculation of optical and dielectric properties in solid targets. We have implemented a model for inverse bremsstrahlung absorption in solid targets. With the extinction coefficient from resonant photoabsorption, photoionization and inverse bremsstrahlung the refractive index is obtained by the Kramers-Kronig relation. Electrons out of equilibrium. We treat the evolution of the electron distribution, when it is not justified to assume a Maxwellian or Fermi-Dirac distribution, by the Fokker-Planck approximation for collisions between electrons. When solving the Fokker-Planck equation we included inelastic collision processes with ions and source terms by interaction with the laser and other processes.

On the possibility of coherently stimulated recombination and cosmological structure generation: recombination instability.

Possible instabilities during cosmological recombination may produce an epoch of nonlinear density growth and fractal-like structural patterns out to the horizon scale at that epoch (approximately 200 Mpc today). With this motivation, we examine the consequences of the change in effective radiative recombination reaction rate coefficients produced by intense stimulated emission. The proton-electron recombination is considered as a natural laser, leading to the formation of spatially nonuniform distributions of neutral matter earlier than the recombination epoch.

Physicochemical characteristics of chlorofluorohydrocarbon based inhalation aerosols

The aim of this work was to gain a better understanding of the physiochemical factors which affect the formulation of suspension inhalation aerosols. This has been attempted by applying the principles of colloid science to aerosol formulation. Both a drug system and a model colloid system have been used. The adsorption of six nonionic and cationic surfactants onto Spherisorb has been investigated. The results were analysed by calculating the area occupied by one adsorbed molecule at the surface and by comparing these values for each surfactant. The amount of each surfactant adsorbed was correlated with the number of sites on that surfactant molecule which could interact with the surface. The stability of suspensions, produced by both the model colloid Spherisorb, and by the drug isoprenaline sulphate, after adsorption of the surfactants, has been assessed by measuring settling times and rising times. The most stable suspensions were found to be those which had the greatest amounts of long chain fatty acid surfactant adsorbed on their surface. A comparison was made between the effective stabilising properties of Span 85 and oleic acid on various drug suspensions. It was found that Span 85 gave the most stable suspensions. Inhalation aerosol suspensions of isoprenaline sulphate were manufactured using the same surfactants used in the adsorption and suspension stability studies and were analysed by measuring the particle size distributions of the suspension and the emitted doses. The results were found to correlate with the adsorption and suspension stability studies and it was concluded that a deflocculated suspension was preferable to a flocculated suspension in inhalation aerosols provided that the drug density was less than the propellant density. The application of this work to preformulation studies was also discussed.

Two-dimensional viscoelastic discrete triangular system with negative-stiffness components

Potential applications of high-damping and high-stiffness composites have motivated extensive research on the effects of negative-stiffness inclusions on the overall properties of composites. Recent theoretical advances have been based on the Hashin-Shtrikman composite models, one-dimensional discrete viscoelastic systems and a two-dimensional nested triangular viscoelastic network. In this paper, we further analyze the two-dimensional triangular structure containing pre-selected negative-stiffness components to study its underlying deformation mechanisms and stability. Major new findings are structure-deformation evolution with respect to the magnitude of negative stiffness under shear loading and the phenomena related to dissipation-induced destabilization and inertia-induced stabilization, according to Lyapunov stability analysis. The evolution shows strong correlations between stiffness anomalies and deformation modes. Our stability results reveal that stable damping peaks, i.e. stably extreme effective damping properties, are achievable under hydrostatic loading when the inertia is greater than a critical value. Moreover, destabilization induced by elemental damping is observed with the critical inertia. Regardless of elemental damping, when the inertia is less than the critical value, a weaker system instability is identified.

Toy model to describe the effect of positional blocklike disorder in metamaterials composites

We study theoretically the effect of a new type of blocklike positional disorder on the effective electromagnetic properties of one-dimensional chains of resonant, high-permittivity dielectric particles, where particles are arranged into perfectly well-ordered blocks whose relative position is a random variable. This creates a finite order correlation length that mimics the situation encountered in metamaterials fabricated through self-assembled techniques, whose structures often display short-range order between near neighbors but long-range disorder, due to stacking defects. Using a spectral theory approach combined with a principal component statistical analysis, we study, in the long-wavelength regime, the evolution of the electromagnetic response when the composite filling fraction and the block size are changed. Modifications in key features of the resonant response (amplitude, width, etc.) are investigated, showing a regime transition for a filling fraction around 50%.

Investigating Marine Boundary Layer Aerosol Budgets and Variability

Thesis (Master's)--University of Washington, 2016-08

Inhomogeneous cloud coverage through the Coulomb explosion of dust in substellar atmospheres

Context. Recent observations of brown dwarf spectroscopic variability in the infrared infer the presence of patchy cloud cover. Aims. This paper proposes a mechanism for producing inhomogeneous cloud coverage due to the depletion of cloud particles through the Coulomb explosion of dust in atmospheric plasma regions. Charged dust grains Coulomb-explode when the electrostatic stress of the grain exceeds its mechanical tensile stress, which results in grains below a critical radius a < a Coul crit being broken up. Methods. This work outlines the criteria required for the Coulomb explosion of dust clouds in substellar atmospheres, the effect on the dust particle size distribution function, and the resulting radiative properties of the atmospheric regions. Results. Our results show that for an atmospheric plasma region with an electron temperature of Te = 10 eV (≈105 K), the critical grain radius varies from 10−7 to 10−4 cm, depending on the grains’ tensile strength. Higher critical radii up to 10−3 cm are attainable for higher electron temperatures. We find that the process produces a bimodal particle size distribution composed of stable nanoscale seed particles and dust particles with a ≥ a Coul crit , with the intervening particle sizes defining a region devoid of dust. As a result, the dust population is depleted, and the clouds become optically thin in the wavelength range 0.1–10 μm, with a characteristic peak that shifts to higher wavelengths as more sub-micrometer particles are destroyed. Conclusions. In an atmosphere populated with a distribution of plasma volumes, this will yield regions of contrasting radiative properties, thereby giving a source of inhomogeneous cloud coverage. The results presented here may also be relevant for dust in supernova remnants and protoplanetary disks.

Measurements of small ice crystals in arctic cirrus: observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC)

In-situ observations on the size and shape of particles in arctic cirrus are less common than those in mid-latitude and tropical cirrus with considerable uncertainty about the contributions of small ice crystals (maximum dimension D<50 µm) to the mass and radiative properties that impact radiative forcing. In situ measurements of small ice crystals in arctic cirrus were made during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) in April 2008 during transits of the National Research Council of Canada Convair-580 between Fairbanks and Barrow, Alaska and during Mixed Phase Arctic Cloud Experiment (MPACE) in October 2004 with the University of North Dakota (UND) Citation over Barrow, Alaska. Concentrations of small ice crystals with D < 50 μm from a Cloud and Aerosol Spectrometer (CAS), a Cloud Droplet Probe (CDP), a Forward Scattering Spectrometer Probe (FSSP), and a two-dimensional stereo probe (2DS) were compared as functions of the concentrations of crystals with D > 100 μm measured by a Cloud Imaging Probe (CIP) and two-dimensional stereo probe (2DS) in order to assess whether the shattering of large ice crystals on protruding components of different probes artificially amplified measurements of small ice crystal concentrations. The dependence of the probe comparison on other variables as CIP N>100 (number concentrations greater than diameter D>100 μm),temperature, relative humidity respect to ice (RHice), dominant habit from the Cloud Particle Imager (CPI), aircraft roll, pitch, true air speed and angle of attack was examined to understand potential causes of discrepancies between probe concentrations. Data collected by these probes were also compared against the data collected by a CAS, CDP and CIP during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) and by a CAS and 2DS during the Tropical Composition, Cloud and Climate Coupling (TC4) missions. During ISDAC, the CAS and FSSP both overestimated measurements of small ice crystals compared to both the CDP and 2DS by 1-2 orders of magnitude. Further, the amount of overestimation increased with the concentrations from the CIP2 (N>100 > 0.1 L-1). There was an unexplained discrepancy in concentrations of small crystals between the CDP and 2DS during ISDAC. In addition, there was a strong dependence on RHice of the average ratios of the N3-50, CAS/N3-50,CDP, N3-50, FSSP096/N3-50,CDP, N3-50, CAS/N3-50,FSSP096, N10-50, CDP/N3-50,2DS, N10-50, FSSP096/N10-50,2DS. Continued studies are needed to understand the discrepancy of these probes.

Global oceanic precipitation: A joint view by TOPEX and the TOPEX microwave radiometer

The TOPEX/POSEIDON mission offers the first opportunity to observe rain cells over the ocean by a dual-frequency radar altimeter (TOPEX) and simultaneously observe their natural radiative properties by a three-frequency radiometer (TOPEX microwave radiometer (TMR)). This work is a feasibility study aimed at understanding the capability and potential of the active/passive TOPEX/TMR system for oceanic rainfall detection. On the basis of past experiences in rain flagging, a joint TOPEX/TMR rain probability index is proposed. This index integrates several advantages of the two sensors and provides a more reliable rain estimate than the radiometer alone. One year's TOPEX/TMR TMR data are used to test the performance of the index. The resulting rain frequency statistics show quantitative agreement with those obtained from the Comprehensive Ocean-Atmosphere Data Set (COADS) in the Intertropical Convergence Zone (ITCZ), while qualitative agreement is found for other regions of the world ocean. A recent finding that the latitudinal frequency of precipitation over the Southern Ocean increases steadily toward the Antarctic continent is confirmed by our result. Annual and seasonal precipitation maps are derived from the index. Notable features revealed include an overall similarity in rainfall pattern from the Pacific, the Atlantic, and the Indian Oceans and a general phase reversal between the two hemispheres, as well as a number of regional anomalies in terms of rain intensity. Comparisons with simultaneous Global Precipitation Climatology Project (GPCP) multisatellite precipitation rate and COADS rain climatology suggest that systematic differences also exist. One example is that the maximum rainfall in the ITCZ of the Indian Ocean appears to be more intensive and concentrated in our result compared to that of the GPCP. Another example is that the annual precipitation produced by TOPEX/TMR is constantly higher than those from GPCP and COADS in the extratropical regions of the northern hemisphere, especially in the northwest Pacific Ocean. Analyses of the seasonal variations of prominent rainy and dry zones in the tropics and subtropics show various behaviors such as systematic migration, expansion and contraction, merging and breakup, and pure intensity variations, The seasonality of regional features is largely influenced by local atmospheric events such as monsoon, storm, or snow activities. The results of this study suggest that TOPEX and its follow-on may serve as a complementary sensor to the special sensor microwave/imager in observing global oceanic precipitation.

Analytical formulae for electromechanical effective properties of 3-1 longitudinally porous piezoelectric materials

A unidirectional fiber composite is considered here, the fibers of which are empty cylindrical holes periodically distributed in a transversely isotropic piezoelectric matrix, The empty-fiber cross-section is circular and the periodicity is the same in two directions at an angle pi/2 or pi/3. Closed-form formulae for all electromechanical effective properties of these 3-1 longitudinally periodic porous piezoelectric materials are presented. The derivation of such expressions is based on the asymptotic homogenization method as a limit of the effective properties of two-phase transversely isotropic parallel fiber-reinforced composites when the fibers properties tend to zero. The plane effective coefficients satisfy the corresponding Schulgasser-Benveniste-Dvorak universal type of relations, A new relation among the antiplane effective constants from the solutions of two antiplane strains and potential local problems is found. This relation is valid for arbitrary shapes of the empty-fiber cross-sections. Based on such a relation, and using recent numerical results for isotropic conductive composites, the antiplane effective properties are computed for different geometrical shapes of the empty-fiber cross-section. Comparisons with other analytical and numerical theories are presented. (c) 2008 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.