Análise quantitativa de dispersão de vazamentos de substâncias inflamáveis e/ou tóxicas em ambientes com barreiras ou semi confinados.

Com o atual desenvolvimento industrial e tecnológico da sociedade, a presença de substâncias inflamáveis e/ou tóxicas aumentou significativamente em um grande número de atividades. A possível dispersão de gases perigosos em instalações de armazenamento ou em operações de transporte representam uma grande ameaça à saúde e ao meio ambiente. Portanto, a caracterização de uma nuvem inflamável e/ou tóxica é um ponto crítico na análise quantitativa de riscos. O objetivo principal desta tese foi fornecer novas perspectivas que pudessem auxiliar analistas de risco envolvidos na análise de dispersões em cenários complexos, por exemplo, cenários com barreiras ou semi-confinados. A revisão bibliográfica mostrou que, tradicionalmente, modelos empíricos e integrais são usados na análise de dispersão de substâncias tóxicas / inflamáveis, fornecendo estimativas rápidas e geralmente confiáveis ao descrever cenários simples (por exemplo, dispersão em ambientes sem obstruções sobre terreno plano). No entanto, recentemente, o uso de ferramentas de CFD para simular dispersões aumentou de forma significativa. Estas ferramentas permitem modelar cenários mais complexos, como os que ocorrem em espaços semi-confinados ou com a presença de barreiras físicas. Entre todas as ferramentas CFD disponíveis, consta na bibliografia que o software FLACS® tem bom desempenho na simulação destes cenários. Porém, como outras ferramentas similares, ainda precisa ser totalmente validado. Após a revisão bibliográfica sobre testes de campo já executados ao longo dos anos, alguns testes foram selecionados para realização de um exame preliminar de desempenho da ferramenta CFD utilizado neste estudo. Foram investigadas as possíveis fontes de incertezas em termos de capacidade de reprodutibilidade, de dependência de malha e análise de sensibilidade das variáveis de entrada e parâmetros de simulação. Os principais resultados desta fase foram moldados como princípios práticos a serem utilizados por analistas de risco ao realizar análise de dispersão com a presença de barreiras utilizando ferramentas CFD. Embora a revisão bibliográfica tenha mostrado alguns dados experimentais disponíveis na literatura, nenhuma das fontes encontradas incluem estudos detalhados sobre como realizar simulações de CFD precisas nem fornecem indicadores precisos de desempenho. Portanto, novos testes de campo foram realizados a fim de oferecer novos dados para estudos de validação mais abrangentes. Testes de campo de dispersão de nuvem de propano (com e sem a presença de barreiras obstruindo o fluxo) foram realizados no campo de treinamento da empresa Can Padró Segurança e Proteção (em Barcelona). Quatro testes foram realizados, consistindo em liberações de propano com vazões de até 0,5 kg/s, com duração de 40 segundos em uma área de descarga de 700 m2. Os testes de campo contribuíram para a reavaliação dos pontos críticos mapeados durante as primeiras fases deste estudo e forneceram dados experimentais para serem utilizados pela comunidade internacional no estudo de dispersão e validação de modelos. Simulações feitas utilizando-se a ferramenta CFD foram comparadas com os dados experimentais obtidos nos testes de campo. Em termos gerais, o simulador mostrou bom desempenho em relação às taxas de concentração da nuvem. O simulador reproduziu com sucesso a geometria complexa e seus efeitos sobre a dispersão da nuvem, mostrando claramente o efeito da barreira na distribuição das concentrações. No entanto, as simulações não foram capazes de representar toda a dinâmica da dispersão no que concerne aos efeitos da variação do vento, uma vez que as nuvens simuladas diluíram mais rapidamente do que nuvens experimentais.

Desenvolvimento da modelagem de turbulência e interação fluido-estrutura para as vibrações induzidas por vórtices de cilindro rígido.

Esta tese apresenta o desenvolvimento e aplicação de modelos de turbulência, transição laminar-turbulenta e de interações fluido-estrutura ao escoamento externo em cilindro rígido estacionário e em vibrações induzidas por vórtices. Tais desenvolvimentos foram realizados no código ReFRESCO, baseado em técnicas de dinâmica de fluidos computacional (CFD). Realizou-se um estudo quanto ao desempenho do modelo k- SST em extensa faixa de números de Reynolds, segundo o qual se identificaram as deficiências de modelagem para este escoamento. A modelagem adaptativa das escalas (SAS) e o modelo de transição por correlações locais (LCTM), ambos combinados ao SST, melhoraram a aderência aos resultados experimentais para este escoamento, em uma contribuição original deste trabalho. A aplicação de técnicas de verificação e validação possibilitou a estimação de incertezas e erros para os modelos e números de Reynolds e também de identificada como outra contribuição deste trabalho. A combinação da modelagem em SST, SAS e LCTM com movimentos impostos de realizada para números de Reynolds moderados, diferentes frequências e amplitudes de vibração, algo que poucas publicações abordam em detalhes. Com relação aos movimentos livres, este trabalho traz contribuições com a aplicação dos modelos SST e SAS ao estudo de vibrações induzidas por vórtices em dois graus de liberdade, baixa razão de massa e números de Reynolds moderados, mais altos do que normalmente observados na literatura. Por fim, a investigação da importância relativa de efeitos da turbulência aos casos de movimentos livres e impostos, com relação ao caso de cilindro estacionário, comprovou a conjetura formulada na parte inicial deste trabalho, no que tange à escolha do modelo de turbulência em determinadas aplicações. Tal escolha mostrou-se menos decisiva no caso do cilindro em movimento imposto e ainda menos nos movimentos livres, em comparação ao caso estacionário, uma vez que a resposta em movimentos do corpo filtra grande parte dos efeitos turbulentos de ordem superior. Esta observação mostra-se relevante, uma vez que pode permitir simplificações na modelagem e aplicação de ferramentas de CFD em uma classe importante de projetos de engenharia.

Otimização acústica e análise numérica do escoamento ao redor de um conjunto cilindro-placa separadora.

RESUMO Simulações de aeroacústica computacional demandam uma quantidade considerável de tempo, o que torna complicada a realização de estudos paramétricos. O presente trabalho propõe uma metodologia viável para otimização aeroacústica. Através da análise numérica utilizando dinâmica dos fluidos computacional, foi estudada a aplicação de uma placa separadora desacoplada como método de controle passivo da esteira turbulenta de um cilindro e avaliou-se a irradiação de ruído causado pela interação do escoamento com ambos os corpos, empregando ferramentas de aeroacústica computacional baseadas no método de Ffowcs-Williams e Hawkings. Algumas abordagens distintas de metodologias de otimização de projeto foram aplicadas neste problema, com o objetivo de chegar a uma configuração otimizada que permita a redução do nível sonoro ao longe. Assim, utilizando uma ferramenta de otimização multidisciplinar, pode-se avaliar a capacidade de modelos heurísticos e a grande vantagem do emprego de algoritmos baseados em método de superfície de resposta quando aplicados em um problema não linear, pois requerem a avaliação de um menor número de alternativas para se obter um ponto ótimo. Além disso, foi possível identificar e agrupar os resultados em 5 clusters baseados em seus parâmetros geométricos, nível de pressão sonora global e o valor quadrático médio do coeficiente de arrasto, confirmando a eficiência da aplicação de placas separadoras longas desacopladas posicionadas próximas ao cilindro na estabilização da esteira turbulenta, enquanto que o posicionamento de placas acima de um espaçamento crítico aumentou o nível de pressão acústica irradiado devido à formação de vórtices no espaço entre o cilindro e a placa separadora.

Coupled Biomechanical Response of the Cornea Assessed by Non-Contact Tonometry. A Simulation Study

The mechanical response of the cornea subjected to a non-contact air-jet tonometry diagnostic test represents an interplay between its geometry, the corneal material behavior and the loading. The objective is to study this interplay to better understand and interpret the results obtained with a non-contact tonometry test. A patient-specific finite element model of a healthy eye, accounting for the load free configuration, was used. The corneal tissue was modeled as an anisotropic hyperelastic material with two preferential directions. Three different sets of parameters within the human experimental range obtained from inflation tests were considered. The influence of the IOP was studied by considering four pressure levels (10–28 mmHg) whereas the influence of corneal thickness was studied by inducing a uniform variation (300–600 microns). A Computer Fluid Dynamics (CFD) air-jet simulation determined pressure loading exerted on the anterior corneal surface. The maximum apex displacement showed a linear variation with IOP for all materials examined. On the contrary, the maximum apex displacement followed a cubic relation with corneal thickness. In addition, a significant sensitivity of the apical displacement to the corneal stiffness was also obtained. Explanation to this behavior was found in the fact that the cornea experiences bending when subjected to an air-puff loading, causing the anterior surface to work in compression whereas the posterior surface works in tension. Hence, collagen fibers located at the anterior surface do not contribute to load bearing. Non-contact tonometry devices give useful information that could be misleading since the corneal deformation is the result of the interaction between the mechanical properties, IOP, and geometry. Therefore, a non-contact tonometry test is not sufficient to evaluate their individual contribution and a complete in-vivo characterization would require more than one test to independently determine the membrane and bending corneal behavior.

Aceleración de corrientes eólicas de reducida velocidad para la generación eléctrica doméstica

En los últimos años, los sistemas que utilizan como fuente recursos renovables se han posicionado como una interesante alternativa para la producción de energía. Entre las fuentes disponibles, la energía eólica viene configurándose como una de las fuentes de energía renovable con mayor crecimiento en los últimos años. En este trabajo se propone el aprovechamiento de las corrientes eólicas circulantes sobre la superficie terrestre para la producción de electricidad con el fin de abastecer buena parte de la demanda en viviendas aisladas, pequeñas instalaciones agropecuarias, equipamiento de servicio ubicado en lugares remotos, etc. Por lo general estas brisas tienen una baja densidad energética, por ello proponemos una interfaz mecánica que concentre las masas de aire, acelerando su circulación y alcanzando importantes incrementos en la velocidad de impulsión. La primera parte se centra en la elaboración de un procedimiento de caracterización a partir de la metodología científica con el cual modelar una estructura concentradora de flujo eólico válida para un aerogenerador de eje vertical. Este método trata el diseño de un elemento acelerador capaz de optimizar el aprovechamiento de estas brisas con independencia de la dirección de éstas. Su diseño viene dado por la resolución de un conjunto de objetivos fundamentales, dotando al sistema de unas prestaciones particulares en relación a su arquitectura y operatividad. Estos objetivos son los siguientes: - Operatividad ante cualquier dirección eólica adoptada - Incremento del rendimiento potencial de la turbina de eje vertical - Minimizar el desarrollo de efectos turbulentos alrededor del sistema integrado - Capacidad resolutiva ante la presencia de fuertes vientos - Estabilidad estructural - Compatibilidad ante instalaciones propias del volumen arquitectónico - Control global del rendimiento del sistema. La segunda parte aborda el modelado del prototipo y el análisis de su comportamiento mediante simulaciones en el ámbito de la fluidodinámica computacional. El resultado es un prototipo caracterizado por una arquitectura capaz de sectorizar la entrada de viento en diferentes tramos inyectando el flujo eólico estratégicamente. La incorporación de la interfaz sobre el rotor aumenta la superficie de captación eólica, facilitando su entrada a través de las diferentes aberturas y llevando a cabo su concentración según avanza por los tramo de circulación. Una vez finalizado dicho avance, el flujo es inyectado en un rango angular de nido por la elevada fuerza de sustentación capaz de generarse gracias a la incidencia del flujo eólico, aprovechando las particularidades que ofrece la rotación de este tipo de rotores. El resultado de la inyección sectorizada es el desarrollo de una circulación interior vorticial que incide permanentemente en el rango de sustentación característico del perfil aerodinámico que define la geometría de la pala de rotación. Ello provoca que se alcance un funcionamiento nominal a velocidades reducidas. En este proceso se incluyen las acciones necesarias para dar una respuesta eficiente a cualquier tipo e solicitación eólica. En presencia de velocidades de relativa importancia, la interfaz concentradora adapta su arquitectura con el fin de regular la entrada de flujo, retrasando la activación de los dispositivos reguladores propios de la turbina eólica. En presencia de vientos importantes, la interfaz dispone de los mecanismos necesarios para proceder al cierre de las aberturas procediendo a la parada del rotor. La validación de los prototipos elaborados se ha llevado a cabo mediante simulación computacional CFD. Los resultados confirman la consecución de un funcionamiento nominal a velocidades más reducidas, una mayor superficie de captación y periodo de tiempo de funcionamiento efectivo en comparación a las turbinas convencionales. Para el caso práctico modelado los resultados mejoran en más de 2.5 veces la potencia generada y multiplican por cuatro la fuerza ejercida sobre las palas de la turbina. La elaboración de un método preciso para el diseño de este tipo de estructuras concentradoras posibilita que se pueda alcanzar un diseño en función a las necesidades del usuario, o las condiciones eólicas de un emplazamiento dado. A ello hay que unir su compatibilidad de uso y montaje con sistemas de captación solar, conformando un sistema híbrido capaz de aprovechar tanto la energía solar como eólica para el abastecimiento autónomo. Esta característica incrementa el ratio de zonas geográficas donde se puede llevar a cabo su implantación. Las últimas páginas están reservadas a esbozar las líneas futuras de desarrollo y evolución, tanto en términos de e ciencia productiva como en su incorporación a nuevos volúmenes arquitectónicos y estructuras civiles en general.

Desarrollo y Aplicación de Herramientas Teóricas, Numéricas y Códigos Computacionales en Dinámica de Gases para Ingeniería Aeronáutica

Objetivos generales: • Desarrollo de un código utilizando programación orientada objeto y Fortran 2003 para simulación de flujos de gases de alta entalpía incluyendo efectos químicos. • Mejora de técnicas numéricas con la finalidad de simular flujos gas dinámicos con alta entalpía y químicamente activos. • Mejora de la comprensión del fenómeno de intermitencia caótica. • Mejorar en la comprensión, implementación y utilización del código openFOAM para la simulación de flujos de interés en la industria aeroespacial. Objetivos específicos: • Poseer la capacidad de simular las ecuaciones de la dinámica de gases tridimensionales, químicamente activas y dependientes del tiempo usando volúmenes finitos mediante códigos propios. • Continuar con el desarrollo de la nueva teoría que permite una descripción más adecuada y precisa del fenómeno de intermitencia como una de las rutas hacia el caos. • Adaptar el código openFoam para la simulación de la dinámica del flujo de gases con elevada entalpía y/o reactivos. • Aplicación de los códigos desarrollados a la simulación de flujos de interés aeronáutico tales como chorros supersónicos, flujos supersónicos alrededor de objetos en vuelo atmosférico y procesos de combustión difusiva. Resultados esperados • Obtención de un código numérico en volúmenes finitos mediante programación orientada a objeto desarrollado en la UNC con la capacidad de simular flujos gas-dinámicos, tridimensionales, viscosos, químicamente activos y dependientes del tiempo. Este código será la base para un programa de multifísica en flujos gaseosos que en un futuro deberá incluir ionización y simulación de flujos magnetohidrodinámicos compresibles. • Desarrollo y/o mejora de técnicas numéricas para la solución de flujos gas-dinámicos químicamente activos. Se continuará con la investigación que usa distintas funciones limitadoras para distintas ondas. • Mejora del conocimiento sobre fenómenos que presentan intermitencia caótica. Se continuará con la teoría que está actualmente en desarrollo por integrantes del grupo en colaboración con investigadores de universidades de España. Principalmente se pondrá énfasis para obtener una formulación más general de la teoría mediante el desarrollo de una técnica de perturbación que permita aproximar cualquier función no lineal que gobierna el proceso de reinyección. • Implementación del código openFOAM para simular flujos supersónicos y procesos de combustión. Se generarán librerías específicas para utilizar este software en aplicaciones de flujos compresibles. • Se espera alcanza una mejor descripción numérica del fenómeno de chorros supersónicos sub-expandidos que impactan sobre paredes sólidas.

Análise de balística externa de um projétil de calibre 155 mm

O objetivo desta Dissertação é propor uma geometria para um projétil de combate a incêndios. É utilizado o código de volume finito Star - CCM+R para calcular forças e momentos que atuam no projétil e é criado um código no MathematicaR para definir as caraterísticas da trajetória. O PRODASV 3 R é utilizado de modo a verificar e validar os modelos implementados. Inicialmente, o procedimento adotado consistiu em criar uma metodologia completa para uma munição existente que permitisse ser verificada e validada com a literatura disponível. Foi realizada uma extensa caraterização aerodinâmica através de CFD, a fim de se estabelecer o comportamento das forças e momentos que atuam nos projéteis e assim fornecer os dados necessários para o modelo 6 DOF. As forças e momentos dependem do ângulo de ataque. Concluiu-se que a cinta de travamento, a cinta de guiamento e o efeito de Coriolis eram desprezáveis quando comparados com os resultados obtidos. Na fase seguinte foram testadas seis geometrias diferentes, variando-se o comprimento e o nariz. Definiu-se para este projétil um a velocidade inicial de 100m/s, 10kg, concluindo-se que o comprimento total seria de 697,1mm e perfil semiesférico, com capacidade para transportar 7,5 dm3 no interior. Foram realizados testes com vento para avaliar a estabilidade em voo, e o software e processo de análise desenvolvido neste trabalho estão disponíveis, sendo um contributo para o projeto FIREND.

Metodologie di calcolo automatico di schiere palettate

La progettazione e la modellazione delle geometrie di corpi complessi come le schiere palettate delle turbomacchine, da sempre impegna il tecnico che, dapprima su carta, e successivamente in forma digitale, deve scontrarsi con le difficoltà sia analitiche di risoluzione di sistemi di equazioni differenziali, che geometriche a causa della doppia curvatura dei profili stessi. L’avvento dei calcolatori ha inevitabilmente giocato un ruolo fondamentale nella rapida evoluzione di tecniche di modellazione, calcolo e rappresentazione, per aiutare il progettista a risolvere completamente il problema, o almeno riscontrare risultati approssimativamente corretti, al fine di ridurre i tempi di realizzazione e i costi dell’impresa. Si vuole dunque cercare di descrivere le fasi che la progettazione oggi richiede, sfruttando quello che i software moderni mettono a disposizione, con l’obiettivo di mostrare uno dei molteplici percorsi che il progettista oggi può seguire per riuscire nel suo scopo.

Produzione di energia elettrica mediante energia marina

La tesi affronta il tema della conversione dell’energia marina in energia elettrica. Tale argomento viene inquadrato nell’ambito della produzione di energia elettrica da fonti rinnovabili, di cui si riportano le stime dei costi relativi, e confronti con la produzione da fonti fossili. La tesi affronta il tema della conversione dell’energia marina discutendone: la definizione, l’analisi economica, i progetti principali in cui sono coinvolte le maggiori imprese del settore, cenni alle stime del potenziale energetico dell’energia marina lungo le coste italiane e agli impianti in corso di realizzazione nel suolo nazionale. Vengono descritte inoltre le caratteristiche principali di un impianto per la conversione di energia marina in energia elettrica, ed in particolare: il funzionamento delle turbine marine e la loro classificazione, il principio di funzionamento di un convertitore di energia ondosa ed in particolare del tipo a galleggiamento. La tesi discute le prospettive di sviluppo della tecnologia in oggetto e le possibili applicazioni future. Infine, è presente un’analisi dei principali benefici e svantaggi dell’energia marina come fonte per la produzione di energia elettrica.

Karel Skreta, maljř : wlastenská směšnohra w trogjm děgstwu /

With: Mašinkář / Jan Nep. Nürnberger. W Praze : Tisk Jar. Pospíšila, 1846 -- Weselý šwec, aneb Proměněné ženy / Jan Nep. Štěpánek. W Praze : Nákl. Wáclawa Špinky, 1841.

Axisymmetric shock wave interaction with a cone: a benchmark test

Results of the benchmark test are presented of comparing numerical schemes solving shock wave of M-s = 2.38 in nitrogen and argon interacting with a 43 degrees semi-apex angle cone and corresponding experiments. The benchmark test was announced in Shock Waves Vol. 12, No. 4, in which we tried to clarify the effects of viscosity and heat conductivity on shock reflection in conical flows. This paper summarizes results of ten numerical and two experimental applications. State of the art in studies regarding the shock/cone interaction is clarified.

A parametric study of droplet deformation through a microfluidic contraction: Low viscosity Newtonian droplets

Numerical simulations are conducted to investigate how a droplet of Newtonian liquid. entrained in a higher viscosity Newtonian liquid, behaves when passing through an axisymmetric microfluidic contraction. Simulations are performed using a transient Volume of Fluid finite volume algorithm, and cover ranges of Reynolds and Weber numbers relevant to microfluidic flows. Results are presented for a droplet to surrounding fluid viscosity ratio of 0.001. In contrast to behaviour at higher viscosity ratios obtained previously by the authors, shear and interfacial tension driven instabilities often develop along the droplet Surface. leading to complex shape development, and in some instances, droplet breakup. (c) 2006 Elsevier Ltd. All rights reserved.

Numerical simulation of magnetite segregation in a dense medium cyclone

Numerical simulations of turbulent driven flow in a dense medium cyclone with magnetite medium have been conducted using Fluent. The predicted air core shape and diameter were found to be close to the experimental results measured by gamma ray tomography. It is possible that the Large eddy simulation (LES) turbulence model with Mixture multi-phase model can be used to predict the air/slurry interface accurately although the LES may need a finer grid. Multi-phase simulations (air/water/medium) are showing appropriate medium segregation effects but are over-predicting the level of segregation compared to that measured by gamma-ray tomography in particular with over prediction of medium concentrations near the wall. Further, investigated the accurate prediction of axial segregation of magnetite using the LES turbulence model together with the multi-phase mixture model and viscosity corrections according to the feed particle loading factor. Addition of lift forces and viscosity correction improved the predictions especially near the wall. Predicted density profiles are very close to gamma ray tomography data showing a clear density drop near the wall. The effect of size distribution of the magnetite has been fully studied. It is interesting to note that the ultra-fine magnetite sizes (i.e. 2 and 7 mu m) are distributed uniformly throughout the cyclone. As the size of magnetite increases, more segregation of magnetite occurs close to the wall. The cut-density (d(50)) of the magnetite segregation is 32 gm, which is expected with superfine magnetite feed size distribution. At higher feed densities the agreement between the [Dungilson, 1999; Wood, J.C., 1990. A performance model for coal-washing dense medium cyclones, Ph.D. Thesis, JKMRC, University of Queensland] correlations and the CFD are reasonably good, but the overflow density is lower than the model predictions. It is believed that the excessive underflow volumetric flow rates are responsible for under prediction of the overflow density. (c) 2006 Elsevier Ltd. All rights reserved.

Study of flow behaviour in a three-product cyclone using computational fluid dynamics

Simplicity in design and minimal floor space requirements render the hydrocyclone the preferred classifier in mineral processing plants. Empirical models have been developed for design and process optimisation but due to the complexity of the flow behaviour in the hydrocyclone these do not provide information on the internal separation mechanisms. To study the interaction of design variables, the flow behaviour needs to be considered, especially when modelling the new three-product cyclone. Computational fluid dynamics (CFD) was used to model the three-product cyclone, in particular the influence of the dual vortex finder arrangement on flow behaviour. From experimental work performed on the UG2 platinum ore, significant differences in the classification performance of the three-product cyclone were noticed with variations in the inner vortex finder length. Because of this simulations were performed for a range of inner vortex finder lengths. Simulations were also conducted on a conventional hydrocyclone of the same size to enable a direct comparison of the flow behaviour between the two cyclone designs. Significantly, high velocities were observed for the three-product cyclone with an inner vortex finder extended deep into the conical section of the cyclone. CFD studies revealed that in the three-product cyclone, a cylindrical shaped air-core is observed similar to conventional hydrocyclones. A constant diameter air-core was observed throughout the inner vortex finder length, while no air-core was present in the annulus. (c) 2006 Elsevier Ltd. All rights reserved.

Numerical analysis of gas and micro-particle interactions in a hand-held shock-tube device

A unique hand-held gene gun is employed for ballistically delivering biomolecules to key cells in the skin and mucosa in the treatment of the major diseases. One of these types of devices, called the Contoured Shock Tube (CST), delivers powdered micro-particles to the skin with a narrow and highly controllable velocity distribution and a nominally uniform spatial distribution. In this paper, we apply a numerical approach to gain new insights in to the behavior of the CST prototype device. The drag correlations proposed by Henderson (1976), Igra and Takayama (1993) and Kurian and Das (1997) were applied to predict the micro-particle transport in a numerically simulated gas flow. Simulated pressure histories agree well with the corresponding static and Pitot pressure measurements, validating the CFD approach. The calculated velocity distributions show a good agreement, with the best prediction from Igra & Takayama correlation (maximum discrepancy of 5%). Key features of the gas dynamics and gas-particle interaction are discussed. Statistic analyses show a tight free-jet particle velocity distribution is achieved (570 +/- 14.7 m/s) for polystyrene particles (39 +/- 1 mu m), representative of a drug payload.