DK-1D: a drift-kinetic simulation tool for modelling the shear Alfvén wave and its interaction with collisionless plasma

We present a highly accurate tool for the simulation of shear Alfven waves (SAW) in collisionless plasma. SAW are important in space plasma environments because for small perpendicular scale lengths they can support an electric field parallel to the ambient magnetic field. Electrons can be accelerated by the parallel electric field and these waves have been implicated as the source of vibrant auroral displays. However, the parallel electric field carried by SAW is small in comparison to the perpendicular electric field of the wave, making it difficult to measure directly in the laboratory, or by satellites in the near-Earth plasma environment. In this paper, we present a simulation code that provides a means to study in detail the SAW-particle interaction in both space and laboratory plasma. Using idealised, small-amplitude propagating waves with a single perpendicular wavenumber, the simulation code accurately reproduces the damping rates and parallel electric field amplitudes predicted by linear theory for varying temperatures and perpendicular scale lengths. We present a rigorous kinetic derivation of the parallel electric field strength for small-amplitude SAW and show that commonly-used inertial and kinetic approximations are valid except for where the ratio of thermal to Alfv\'{e}n speed is between 0.7 and 1.0. We also present nonlinear simulations of large-amplitude waves and show that in cases of strong damping, the damping rates and parallel electric field strength deviate from linear predictions when wave energies are greater than only a few percent of the plasma kinetic energy, a situation which is often observed in the magnetosphere. The drift-kinetic code provides reliable, testable predictions of the parallel electric field strength which can be investigated directly in the laboratory, and will help to bridge the gap between studies of SAW in man-made and naturally occuring plasma.

The mineralizing fluid in the Piaba gold deposit, São Luís cratonic fragment (NW-Maranhão, Brazil) based on fluid inclusion studies in quartz veins

O depósito aurífero de Piaba tornou-se a primeira mina em operação no fragmento cratônico São Luís, noroeste do Maranhão. Seu ambiente geológico compreende rochas metavulcanossedimentares do Grupo Aurizona e granitoides da Suíte Tromaí, entre outras unidades menores, formadas em ambiente de arcos de ilhas entre 2240 e 2150 Ma, juntamente com outras unidades menores. A mineralização é caracterizada por uma trama stockwork de veios e vênulas de quartzo com seus halos de alteração (clorita + muscovita + carbonato + pirita + calcopirita e ouro) hospedada em um granodiorito granofírico fino (Granófiro Piaba) e em rocha subvulcânica andesítica do Grupo Aurizona. O corpo mineralizado é espacialmente limitado à zona de cisalhamento ENE-WSW rúptil-dúctil (Falha Piaba). Estudos petrográficos, microtermométricos e por espectroscopia microRaman no quartzo definiram inclusões aquo-carbônicas bifásicas e trifásicas, produzidas por aprisionamento heterogêneo durante separação de fases, e fluidos aquosos tardios. A solução mineralizadora corresponde a um fluido aquo-carbônico composto por CO₂ (5 - 24 mol%, densidade de 0,96 - 0,99 g/cm³), H₂O (74 - 93 mol%), N₂ (< 1 mol%), CH₄ (<1mol%) e 5,5 % em peso NaCl equivalente. O minério depositou a 267 - 302ºC e 1,25 - 2,08 kbar, correspondendo a profundidades de 4 a 7 km, em consonância com o regime estrutural. A composição e o intervalo de P-T do fluido mineralizador, combinados com o caráter redutor (log ƒ_O2 -31,3 a -34,3) e a sulfetação da rocha hospedeira, sugerem que o ouro foi transportado como um complexo sulfetado. O minério foi depositado em consequência da separação de fase, redução da atividade de enxofre e da ƒO₂ pela interação fluido-rocha.

Evaluation of Wave energy transmission through a floating breakwater using the SPH method

[EN]Energy transmission through a box-shaped floating breakwater (FB) is examined, under simplified conditions, by using the smoothed particle hydrodynamics (SPH) method, a mesh-free particle numerical approach. The efficiency of the structure is assessed in terms of the coefficient of transm ission as a function of the wave period and the location of the floating breakwater relative to the zone to be protected. Preliminary results conceming wave energy transmission reveals a clear improvement of the efficiency as wave period decreases andan important role ofthe bathymetry.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Hydroelastic impacts of deformable wedges

This work investigates the slamming phenomenon experienced during the water entry of deformable bodies. Wedges are chosen as reference geometry due to their similarity to a generic hull section. Hull slamming is a phenomenon occurring when a ship re-enters the water after having been partially or completely lifted out the water. While the analysis of rigid structures entering the water has been extensively studied in the past and there are analytical solutions capable of correctly predicting the hydrodynamic pressure distribution and the overall impact dynamics, the effect of the structural deformation on the structural force is still a challenging problem to be solved. In fact, in case of water impact of deformable bodies, the dynamic deflection could interact with the fluid flow, changing the hydrodynamic load. This work investigates the hull-slamming problem by experiments and numerical simulations of the water entry of elastic wedges impacting on an initially calm surface. The effect of asymmetry due to horizontal velocity component or initial tilt angle on the impact dynamics is also studied. The objective of this work is to determine an accurate model to predict the overall dynamics of the wedge and its deformations. More than 1200 experiments were conducted by varying wedge structural stiffness, deadrise angle, impact velocity and mass. On interest are the overall impact dynamics and the local structural deformation of the panels composing the wedge. Alongside with the experimental analysis, numerical simulations based on a coupled Smoothed Particle Hydrodynamics (SPH) and FEM method are developed. The experimental results provide evidence of the mutual interaction between hydrodynamic load and structural deformation. It is found a simple criterion for the onset of fluid structure interaction (FSI), giving reliable information on the cases where FSI should been taken into account.

Aeroelastic stability of structures: flutter analysis using Computational Fluid Dynamics

Thanks to the increasing slenderness and lightness allowed by new construction techniques and materials, the effects of wind on structures became in the last decades a research field of great importance in Civil Engineering. Thanks to the advances in computers power, the numerical simulation of wind tunnel tests has became a valid complementary activity and an attractive alternative for the future. Due to its flexibility, during the last years, the computational approach gained importance with respect to the traditional experimental investigation. However, still today, the computational approach to fluid-structure interaction problems is not as widely adopted as it could be expected. The main reason for this lies in the difficulties encountered in the numerical simulation of the turbulent, unsteady flow conditions generally encountered around bluff bodies. This thesis aims at providing a guide to the numerical simulation of bridge deck aerodynamic and aeroelastic behaviour describing in detail the simulation strategies and setting guidelines useful for the interpretation of the results.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

Numerical issues in Lagrangian tracking and topological evolution of fluid particles in wallbounded turbulent flows

The determination of the local Lagrangian evolution of the flow topology in wall-bounded turbulence, and of the Lagrangian evolution associated with entrainment across the turbulent / non-turbulent interface into a turbulent boundary layer, require accurate tracking of a fluid particle and its local velocity gradients. This paper addresses the implementation of fluid-particle tracking in both a turbulent boundary layer direct numerical simulation and in a fully developed channel flow simulation. Determination of the sub-grid particle velocity is performed using both cubic B-spline, four-point Hermite spline and higher-order Hermite spline interpolation. Both wall-bounded flows show similar oscillations in the Lagrangian tracers of both velocity and velocity gradients, corresponding to the movement of particles across the boundaries of computational cells. While these oscillation in the particle velocity are relatively small and have negligible effect on the particle trajectories for time-steps of the order of CFL = 0.1, they appear to be the cause of significant oscillations in the evolution of the invariants of the velocity gradient tensor.

Dynamic interactions between solids and viscous liquids with free surface = Interacciones dinámicas entre sólidos y líquidos viscosos con superficie libre

En esta tesis se investiga la interacción entre un fluido viscoso y un cuerpo sólido en presencia de una superficie libre. El problema se expresa teóricamente poniendo especial atención a los aspectos de conservación de energía y de la interacción del fluido con el cuerpo. El problema se considera 2D y monofásico, y un desarrollo matemático permite una descomposición de los términos disipativos en términos relacionados con la superficie libre y términos relacionados con la enstrofía. El modelo numérico utilizado en la tesis se basa en el método sin malla Smoothed Particle Hydrodynamics (SPH). De manera análoga a lo que se hace a nivel continuo, las propiedades de conservación se estudian en la tesis con el sistema discreto de partículas. Se tratan también las condiciones de contorno de un cuerpo que se mueve en un flujo viscoso, implementadas con el método ghost-fluid. Se ha desarrollado un algoritmo explícito de interacción fluido / cuerpo. Se han documentado algunos casos de modo detallado con el objetivo de comprobar la capacidad del modelo para reproducir correctamente la disipación de energía y el movimiento del cuerpo. En particular se ha investigado la atenuación de una onda estacionaria, comparando la simulación numérica con predicciones teóricas. Se han realizado otras pruebas para monitorizar la disipación de energía para flujos más violentos que implican la fragmentación de la superficie libre. La cantidad de energía disipada con los diferentes términos se ha evaluado en los casos estudiados con el modelo numérico. Se han realizado otras pruebas numéricas para verificar la técnica de modelización de la interacción fluido / cuerpo, concretamente las fuerzas ejercidas por las olas en cuerpos con formas simples, y el equilibrio de un cuerpo flotante con una forma compleja. Una vez que el modelo numérico ha sido validado, se han realizado simulaciones numéricas para obtener una comprensión más completa de la física implicada en casos (casi) realistas sobre los había aspectos que no se conocían suficientemente. En primer lugar se ha estudiado el el flujo alrededor de un cilindro bajo la superficie libre. El estudio se ha realizado con un número de Reynolds moderado, para un rango de inmersiones del cilindro y números de Froude. La solución numérica permite una investigación de los patrones complejos que se producen. La estela del cilindro interactúa con la superficie libre. Se han identificado algunos inestabilidades características. El segundo estudio se ha realizado sobre el problema de sloshing, tanto experimentalmente como numéricamente. El análisis se restringe a aguas poco profundas y con oscilación horizontal, pero se ha estudiado un gran número de condiciones, lo que lleva a una comprensión bastante completa de los sistemas de onda involucradas. La última parte de la tesis trata también sobre un problema de sloshing pero esta vez el tanque está oscilando con rotación y hay acoplamiento con un sistema mecánico. El sistema se llama pendulum-TLD (Tuned Liquid Damper - con líquido amortiguador). Este tipo de sistema se utiliza normalmente para la amortiguación de las estructuras civiles. El análisis se ha realizado analíticamente, numéricamente y experimentalmente utilizando líquidos con viscosidades diferentes, centrándose en características no lineales y mecanismos de disipación. ABSTRA C T The subject of the present thesis is the interaction between a viscous fluid and a solid body in the presence of a free surface. The problem is expressed first theoretically with a particular focus on the energy conservation and the fluid-body interaction. The problem is considered 2D and monophasic, and some mathematical development allows for a decomposition of the energy dissipation into terms related to the Free Surface and others related to the enstrophy. The numerical model used on the thesis is based on Smoothed Particle Hydrodynamics (SPH): a computational method that works by dividing the fluid into particles. Analogously to what is done at continuum level, the conservation properties are studied on the discrete system of particles. Additionally the boundary conditions for a moving body in a viscous flow are treated and discussed using the ghost-fluid method. An explicit algorithm for handling fluid-body coupling is also developed. Following these theoretical developments on the numerical model, some test cases are devised in order to test the ability of the model to correctly reproduce the energy dissipation and the motion of the body. The attenuation of a standing wave is used to compare what is numerically simulated to what is theoretically predicted. Further tests are done in order to monitor the energy dissipation in case of more violent flows involving the fragmentation of the free-surface. The amount of energy dissipated with the different terms is assessed with the numerical model. Other numerical tests are performed in order to test the fluid/body interaction method: forces exerted by waves on simple shapes, and equilibrium of a floating body with a complex shape. Once the numerical model has been validated, numerical tests are performed in order to get a more complete understanding of the physics involved in (almost) realistic cases. First a study is performed on the flow passing a cylinder under the free surface. The study is performed at moderate Reynolds numbers, for various cylinder submergences, and various Froude numbers. The capacity of the numerical solver allows for an investigation of the complex patterns which occur. The wake from the cylinder interacts with the free surface, and some characteristical flow mechanisms are identified. The second study is done on the sloshing problem, both experimentally and numerically. The analysis is restrained to shallow water and horizontal excitation, but a large number of conditions are studied, leading to quite a complete understanding of the wave systems involved. The last part of the thesis still involves a sloshing problem but this time the tank is rolling and there is coupling with a mechanical system. The system is named pendulum-TLD (Tuned Liquid Damper). This kind of system is normally used for damping of civil structures. The analysis is then performed analytically, numerically and experimentally for using liquids with different viscosities, focusing on non-linear features and dissipation mechanisms.

Pulsatile blood flow interaction with arterial walls of aorta : autoregulation and impedance pressure boundary condition and its biomedical applications

Para las decisiones urgentes sobre intervenciones quirúrgicas en el sistema cardiovascular se necesitan simulaciones computacionales con resultados fiables y que consuman un tiempo de cálculo razonable. Durante años los investigadores han trabajado en diversos métodos numéricos de cálculo que resulten atractivos para los cirujanos. Estos métodos, precisos pero costosos desde el punto de vista del coste computacional, crean un desajuste entre la oferta de los ingenieros que realizan las simulaciones y los médicos que operan en el quirófano. Por otra parte, los métodos de cálculo más simplificados reducen el tiempo de cálculo pero pueden proporcionar resultados no realistas. El objetivo de esta tesis es combinar los conceptos de autorregulación e impedancia del sistema circulatorio, la interacción flujo sanguíneo-pared arterial y modelos geométricos idealizados tridimensionales de las arterias pero sin pérdida de realismo, con objeto de proponer una metodología de simulación que proporcione resultados correctos y completos, con tiempos de cálculo moderados. En las simulaciones numéricas, las condiciones de contorno basadas en historias de presión presentan inconvenientes por ser difícil conocerlas con detalle, y porque los resultados son muy sensibles ante pequeñas variaciones de dichas historias. La metodología propuesta se basa en los conceptos de autorregulación, para imponer la demanda de flujo aguas abajo del modelo en el ciclo cardiaco, y la impedancia, para representar el efecto que ejerce el flujo en el resto del sistema circulatorio sobre las arterias modeladas. De este modo las historias de presión en el contorno son resultados del cálculo, que se obtienen de manera iterativa. El método propuesto se aplica en una geometría idealizada del arco aórtico sin patologías y en otra geometría correspondiente a una disección Stanford de tipo A, considerando la interacción del flujo pulsátil con las paredes arteriales. El efecto de los tejidos circundantes también se incorpora en los modelos. También se hacen aplicaciones considerando la interacción en una geometría especifica de un paciente anciano que proviene de una tomografía computarizada. Finalmente se analiza una disección Stanford tipo B con tres modelos que incluyen la fenestración del saco. Clinicians demand fast and reliable numerical results of cardiovascular biomechanic simulations for their urgent pre-surgery decissions. Researchers during many years have work on different numerical methods in order to attract the clinicians' confidence to their colorful contours. Though precise but expensive and time-consuming methodologies create a gap between numerical biomechanics and hospital personnel. On the other hand, simulation simplifications with the aim of reduction in computational time may cause in production of unrealistic outcomes. The main objective of the current investigation is to combine ideas such as autoregulation, impedance, fluid-solid interaction and idealized geometries in order to propose a computationally cheap methodology without excessive or unrealistic simplifications. The pressure boundary conditions are critical and polemic in numerical simulations of cardiovascular system, in which a specific arterial site is of interest and the rest of the netwrok is neglected but represented by a boundary condition. The proposed methodology is a pressure boundary condition which takes advantage of numerical simplicity of application of an imposed pressure boundary condition on outlets, while it includes more sophisticated concepts such as autoregulation and impedance to gain more realistic results. Incorporation of autoregulation and impedance converts the pressure boundary conditions to an active and dynamic boundary conditions, receiving feedback from the results during the numerical calculations and comparing them with the physiological requirements. On the other hand, the impedance boundary condition defines the shapes of the pressure history curves applied at outlets. The applications of the proposed method are seen on idealized geometry of the healthy arotic arch as well as idealized Stanford type A dissection, considering the interaction of the arterial walls with the pulsatile blood flow. The effect of surrounding tissues is incorporated and studied in the models. The simulations continue with FSI analysis of a patient-specific CT scanned geometry of an old individual. Finally, inspiring of the statistic results of mortality rates in Stanford type B dissection, three models of fenestrated dissection sac is studied and discussed. Applying the developed boundary condition, an alternative hypothesis is proposed by the author with respect to the decrease in mortality rates in patients with fenestrations.

Estudo experimental do fenômeno de vibração induzida por vórtices em cilindro rígido livre para oscilar com dois graus de liberdade.

O fenômeno de vibração induzida por vórtices (VIV) é um problema fundamental dentro da Mecânica dos Fluidos e um exemplo importante de interação fluido-estrutura. Esta tese investiga fenômeno de VIV quando um cilindro rígido, submetido a escoamento uniforme, está livre para oscilar na direção transversal e alinhada com a corrente incidente. A tese foi estruturada ao redor de sete perguntas relacionadas ao fenômeno de VIV: 1) O fenômeno e os resultados experimentais são repetitivos? 2) Como ocorre a transição entre ramos de resposta? 3) Qual é o papel da inércia da estrutura oscilante? 4) Qual é o papel de sua rigidez? 5) Quais são as frequências naturais mais importantes da estrutura? 6) Quais padrões de esteira se desenvolvem para VIV com dois graus de liberdade? 7) Quais são os efeitos do movimento na direção alinhada com a corrente no processo de formação e desprendimento de vórtices? O fenômeno de VIV é estudado de maneira experimental em uma base elástica pendular capaz de oscilar com o mesmo momento de inércia e frequência natural nas duas direções. Os experimentos de VIV foram realizados em canal de água recirculante e com diferentes condições de inércia e rigidez. A técnica de velocimetria por imagem de partículas foi usada e permitiu identificar diferentes padrões de esteira de vórtices. Verificou-se que o VIV é repetitivo a nível de amplitudes médias e frequências dominantes. A transição dos ramos pode ocorrer de maneira intermitente ou com histerese. Os parâmetros de inércia e rigidez da estrutura são capazes de mudar o regime de oscilação e, para algumas condições, suprimir as vibrações alinhadas com a corrente. Dentre os padrões de esteira observados, um deles não havia sido relatado na literatura e é definido nesta tese. O novo modo de emissão apresenta dois vórtices com circulação oposta e elevada intensidade emitidos por ciclo. A influência da direção alinhada com o escoamento está relacionada a dois efeitos: a velocidade relativa entre o cilindro e o fluido, responsável pelo aumento da circulação dos vórtices na esteira, e o ângulo de fase do movimento nas direções alinhada e transversal, capaz de mudar o processo de formação dos vórtices.

Comparison of Lattice-Fluid Binary Parameters For Mixtures and Block Copolymers

The aim of this report is to discuss the method of determination of lattice-fluid binary interaction parameters by comparing well characterized immiscible blends and block copolymers of poly(methyl methacrylate) (PMMA) and poly(ϵ−caprolactone) (PCL). Experimental pressure-volume-temperature (PVT) data in the liquid state were correlated with the Sanchez—Lacombe (SL) equation of state with the scaling parameters for mixtures and copolymers obtained through combination rules of the characteristic parameters for the pure homopolymers. The lattice-fluid binary parameters for energy and volume were higher than those of block copolymers implying that the copolymers were more compatible due to the chemical links between the blocks. Therefore, a common parameter cannot account for both homopolymer blend and block copolymer phase behaviors based on current theory. As we were able to adjust all data of the mixtures with a single set of lattice-binary parameters and all data of the block copolymers with another single set we can conclude that both parameters did not depend on the composition for this system. This characteristic, plus the fact that the additivity law of specific volumes can be suitably applied for this system, allowed us to model the behavior of the immiscible blend with the SL equation of state. In addition, a discussion on the relationship between lattice-fluid binary parameters and the Flory–Huggins interaction parameter obtained from Leibler's theory is presented.

Adsorption of Lennard-Jones fluids in carbon slit pores of a finite length. A computer simulation study

The adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length was studied with Canonical Ensemble (NVT) and Gibbs Ensemble Monte Carlo Simulations (GEMC). The Canonical Ensemble was a collection of cubic simulation boxes in which a finite pore resides, while the Gibbs Ensemble was that of the pore space of the finite pore. Argon was used as a model for Lennard-Jones fluids, while the adsorbent was modelled as a finite carbon slit pore whose two walls were composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. The Lennard-Jones (LJ) 12-6 potential model was used to compute the interaction energy between two fluid particles, and also between a fluid particle and a carbon atom. Argon adsorption isotherms were obtained at 87.3 K for pore widths of 1.0, 1.5 and 2.0 nm using both Canonical and Gibbs Ensembles. These results were compared with isotherms obtained with corresponding infinite pores using Grand Canonical Ensembles. The effects of the number of cycles necessary to reach equilibrium, the initial allocation of particles, the displacement step and the simulation box size were particularly investigated in the Monte Carlo simulation with Canonical Ensembles. Of these parameters, the displacement step had the most significant effect on the performance of the Monte Carlo simulation. The simulation box size was also important, especially at low pressures at which the size must be sufficiently large to have a statistically acceptable number of particles in the bulk phase. Finally, it was found that the Canonical Ensemble and the Gibbs Ensemble both yielded the same isotherm (within statistical error); however, the computation time for GEMC was shorter than that for canonical ensemble simulation. However, the latter method described the proper interface between the reservoir and the adsorbed phase (and hence the meniscus).

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.

Theoretical modeling and experimental validation of fluid-loaded micro rectangular plates

This paper presents a theoretical model on the vibration analysis of micro scale fluid-loaded rectangular isotropic plates, based on the Lamb's assumption of fluid-structure interaction and the Rayleigh-Ritz energy method. An analytical solution for this model is proposed, which can be applied to most cases of boundary conditions. The dynamical experimental data of a series of microfabricated silicon plates are obtained using a base-excitation dynamic testing facility. The natural frequencies and mode shapes in the experimental results are in good agreement with the theoretical simulations for the lower order modes. The presented theoretical and experimental investigations on the vibration characteristics of the micro scale plates are of particular interest in the design of microplate based biosensing devices. Copyright © 2009 by ASME.