Investigation of chaotic mixing in laminar fluid systems by the use of computational fluid dynamics

This work presents significant development into chaotic mixing induced through periodic boundaries and twisting flows. Three-dimensional closed and throughput domains are shown to exhibit chaotic motion under both time periodic and time independent boundary motions, A property is developed originating from a signature of chaos, sensitive dependence to initial conditions, which successfully quantifies the degree of disorder withjn the mixing systems presented and enables comparisons of the disorder throughout ranges of operating parameters, This work omits physical experimental results but presents significant computational investigation into chaotic systems using commercial computational fluid dynamics techniques. Physical experiments with chaotic mixing systems are, by their very nature, difficult to extract information beyond the recognition that disorder does, does not of partially occurs. The initial aim of this work is to observe whether it is possible to accurately simulate previously published physical experimental results through using commercial CFD techniques. This is shown to be possible for simple two-dimensional systems with time periodic wall movements. From this, and subsequent macro and microscopic observations of flow regimes, a simple explanation is developed for how boundary operating parameters affect the system disorder. Consider the classic two-dimensional rectangular cavity with time periodic velocity of the upper and lower walls, causing two opposing streamline motions. The degree of disorder within the system is related to the magnitude of displacement of individual particles within these opposing streamlines. The rationale is then employed in this work to develop and investigate more complex three-dimensional mixing systems that exhibit throughputs and time independence and are therefore more realistic and a significant advance towards designing chaotic mixers for process industries. Domains inducing chaotic motion through twisting flows are also briefly considered. This work concludes by offering possible advancements to the property developed to quantify disorder and suggestions of domains and associated boundary conditions that are expected to produce chaotic mixing.

Investigating viscous fluid flow in an internal mixer using computational fluid dynamics

This thesis presents an effective methodology for the generation of a simulation which can be used to increase the understanding of viscous fluid processing equipment and aid in their development, design and optimisation. The Hampden RAPRA Torque Rheometer internal batch twin rotor mixer has been simulated with a view to establishing model accuracies, limitations, practicalities and uses. As this research progressed, via the analyses several 'snap-shot' analysis of several rotor configurations using the commercial code Polyflow, it was evident that the model was of some worth and its predictions are in good agreement with the validation experiments, however, several major restrictions were identified. These included poor element form, high man-hour requirements for the construction of each geometry and the absence of the transient term in these models. All, or at least some, of these limitations apply to the numerous attempts to model internal mixes by other researchers and it was clear that there was no generally accepted methodology to provide a practical three-dimensional model which has been adequately validated. This research, unlike others, presents a full complex three-dimensional, transient, non-isothermal, generalised non-Newtonian simulation with wall slip which overcomes these limitations using unmatched ridding and sliding mesh technology adapted from CFX codes. This method yields good element form and, since only one geometry has to be constructed to represent the entire rotor cycle, is extremely beneficial for detailed flow field analysis when used in conjunction with user defined programmes and automatic geometry parameterisation (AGP), and improves accuracy for investigating equipment design and operation conditions. Model validation has been identified as an area which has been neglected by other researchers in this field, especially for time dependent geometries, and has been rigorously pursued in terms of qualitative and quantitative velocity vector analysis of the isothermal, full fill mixing of generalised non-Newtonian fluids, as well as torque comparison, with a relatively high degree of success. This indicates that CFD models of this type can be accurate and perhaps have not been validated to this extent previously because of the inherent difficulties arising from most real processes.

The development of a novel Bio-MEMS filtration chip for the separation of specific cells in fluid suspension

In the clinical/microbiological laboratory there are currently several ways of separating specific cells from a fluid suspension. Conventionally cells can be separated based on size, density, electrical charge, light-scattering properties, and antigenic surface properties. Separating cells using these parameters can require complex technologies and specialist equipment. This paper proposes new Bio-MEMS (microelectromechanical systems) filtration chips manufactured using deep reactive ion etching (DRIE) technology that, when used in conjunction with an optical microscope and a syringe, can filter and grade cells for size without the requirement for additional expensive equipment. These chips also offer great versatility in terms of design and their low cost allows them to be disposable, eliminating sample contamination. The pumping mechanism, unlike many other current filtration techniques, leaves samples mechanically and chemically undamaged. In this paper the principles behind harnessing passive pumping are explored, modelled, and validated against empirical data, and their integration into a microfluidic device to separate cells from a mixed population suspension is described. The design, means of manufacture, and results from preliminary tests are also presented. © IMechE 2007.

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis.

Quantum transition state theory and solving a first order master equation for complex reactions numerically

The field of chemical kinetics is an exciting and active field. The prevailing theories make a number of simplifying assumptions that do not always hold in actual cases. Another current problem concerns a development of efficient numerical algorithms for solving the master equations that arise in the description of complex reactions. The objective of the present work is to furnish a completely general and exact theory of reaction rates, in a form reminiscent of transition state theory, valid for all fluid phases and also to develop a computer program that can solve complex reactions by finding the concentrations of all participating substances as a function of time. To do so, the full quantum scattering theory is used for deriving the exact rate law, and then the resulting cumulative reaction probability is put into several equivalent forms that take into account all relativistic effects if applicable, including one that is strongly reminiscent of transition state theory, but includes corrections from scattering theory. Then two programs, one for solving complex reactions, the other for solving first order linear kinetic master equations to solve them, have been developed and tested for simple applications.

Computational fluid dynamics (CFD) based approach to consequence assessment of accidental release of hydrocarbon during storage and transportation

This thesis investigated the risk of accidental release of hydrocarbons during transportation and storage. Transportation of hydrocarbons from an offshore platform to processing units through subsea pipelines involves risk of release due to pipeline leakage resulting from corrosion, plastic deformation caused by seabed shakedown or damaged by contact with drifting iceberg. The environmental impacts of hydrocarbon dispersion can be severe. Overall safety and economic concerns of pipeline leakage at subsea environment are immense. A large leak can be detected by employing conventional technology such as, radar, intelligent pigging or chemical tracer but in a remote location like subsea or arctic, a small chronic leak may be undetected for a period of time. In case of storage, an accidental release of hydrocarbon from the storage tank could lead pool fire; further it could escalate to domino effects. This chain of accidents may lead to extremely severe consequences. Analyzing past accident scenarios it is observed that more than half of the industrial domino accidents involved fire as a primary event, and some other factors for instance, wind speed and direction, fuel type and engulfment of the compound. In this thesis, a computational fluid dynamics (CFD) approach is taken to model the subsea pipeline leak and the pool fire from a storage tank. A commercial software package ANSYS FLUENT Workbench 15 is used to model the subsea pipeline leakage. The CFD simulation results of four different types of fluids showed that the static pressure and pressure gradient along the axial length of the pipeline have a sharp signature variation near the leak orifice at steady state condition. Transient simulation is performed to obtain the acoustic signature of the pipe near leak orifice. The power spectral density (PSD) of acoustic signal is strong near the leak orifice and it dissipates as the distance and orientation from the leak orifice increase. The high-pressure fluid flow generates more noise than the low-pressure fluid flow. In order to model the pool fire from the storage tank, ANSYS CFX Workbench 14 is used. The CFD results show that the wind speed has significant contribution on the behavior of pool fire and its domino effects. The radiation contours are also obtained from CFD post processing, which can be applied for risk analysis. The outcome of this study will be helpful for better understanding of the domino effects of pool fire in complex geometrical settings of process industries. The attempt to reduce and prevent risks is discussed based on the results obtained from the numerical simulations of the numerical models.

Geochemistry of sheeted dike complex minerals of ODP Hole 140-504B

Diabases were recovered during Legs 137 and 140 at Hole 504B from depths between 1621.5 and 2000.4 meters below seafloor in the lower sheeted dike complex. The samples contain multiple generations of millimetric to centimetric veins. The orientation of the measured veins suggests that two main vein sets exist: one characterized by shallow dipping and the other by random trend. Thermal contraction during rock cooling is considered the main mechanism responsible for fracture formation. Vein infill is related to the circulation of hydrothermal fluids near the spreading axis. Some veins are surrounded by millimeter-sized alteration halos due to fluid percolation from the fractures through the host rock. Vein-filling minerals are essentially amphibole, chlorite, and zeolites. Amphibole composition is controlled by the microstructural site of the rock. Actinolite is the main amphibole occurring in the veins and also in the groundmass away from the halos. In the alteration halos, amphibole shows composition of actinolitic hornblende and Mg-hornblende. Late-stage tension gashes and interstitial spaces in some amphibole-bearing veins are filled with zeolites, suggesting that the veins likely suffered multiple opening stages that record the cooling history of the circulating fluids. Evidence of deformation recorded by the recovered samples seems to be restricted to veins that clearly represent elements of weakness of the rock. On the basis of vein geometry and microstructure we infer structural interpretations for the formation mechanism and for deformation of veins.

Polymer Fluid Dynamics: Continuum and Molecular Appoaches

To solve problems in polymer fluid dynamics, one needs the equation of continuity, motion, and energy. The last two equations contain the stress tensor and the heat-flux vector for the material. There are two ways to formulate the stress tensor: (1) one can write a continuum expression for the stress tensor in terms of kinematic tensors, or (2) one can select a molecular model that represents the polymer molecule, and then develop an expression for the stress tensor from kinetic theory. The advantage of the kinetic theory approach is that one gets information about the relation between the molecular structure of the polymers and the rheological properties. In this review, we restrict the discussion primarily to the simplest stress tensor expressions or “constitutive equations” containing from two to four adjustable parameters, although we do indicate how these formulations may be extended to give more complicated expressions. We also explore how these simplest expressions are recovered as special cases of a more general framework, the Oldroyd 8-constant model. The virtue of studying the simplest models is that we can discover some general notions as to which types of empiricisms or which types of molecular models seem to be worth investigating further. We also explore equivalences between continuum and molecular approaches. We restrict the discussion to several types of simple flows, such as shearing flows and extensional flows. These are the flows that are of greatest importance in industrial operations. Furthermore, if these simple flows cannot be well described by continuum or molecular models, then it is not necessary to lavish time and energy to apply them to more complex flow problems.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

Induced Damping on Vibrating Circular Plates Submerged in still Fluid

When a structure vibrates immersed in a fluid it is known that the dynamic properties of the system are modified. The surrounding fluid will, in general, contribute to the inertia, the rigidity and the damping coefficient of the coupled fluid-structure system. For light structures, like spacecraft antennas, even when the fluid is air the contribution to the dynamic properties can be important. For not so light structures the ratio of the equivalent fluid/structure mass and rigidity can be very small and the fluid contribution could be neglected. For the ratio of equivalent fluid/structure damping both terms are of the same order and therefore the fluid contribution must be studied. The working life of the spacecraft structure would be on space and so without any surrounding fluid. The response of a spacecraft structure on its operational life would be attenuated by the structural damping alone but when the structure is dynamically tested on the earth the dynamic modal test is performed with the fluid surrounding it. The results thus are contaminated by the effects of the fluid. If the damping added by the fluid is of the same order as the structural damping the response of the structure in space can be quite different to the response predicted on earth. It is therefore desirable to have a method able to determine the amount of damping induced by the fluid and that should be subtracted of the total damping measured on the modal vibration test. In this work, a method for the determination of the effect of the surrounding fluid on the dynamic characteristics of a circular plate has been developed. The plate is assumed to vibrate harmonically with the vacuum modes and the generalized forces matrix due to the fluid is thus computed. For a compressible fluid this matrix is formed by complex numbers including terms of inertia, rigidity and damping. The matrix due to the fluid loading is determined by a boundary element method (BEM). The BEM used is of circular rings on the plate surface so the number of elements to obtain an accurate result is very low. The natural frequencies of the system are computed by an iteration procedure one by one and also the damping fluid contribution. Comparisons of the present method with various experimental data and other theories show the efficiency and accuracy of the method for any support condition of the plate.

Development of a Coupling Model for Fluid-Structure Interaction using the Mesh-free Finite Element Method and the Lattice Boltzmann Method

In the presented thesis work, the meshfree method with distance fields was coupled with the lattice Boltzmann method to obtain solutions of fluid-structure interaction problems. The thesis work involved development and implementation of numerical algorithms, data structure, and software. Numerical and computational properties of the coupling algorithm combining the meshfree method with distance fields and the lattice Boltzmann method were investigated. Convergence and accuracy of the methodology was validated by analytical solutions. The research was focused on fluid-structure interaction solutions in complex, mesh-resistant domains as both the lattice Boltzmann method and the meshfree method with distance fields are particularly adept in these situations. Furthermore, the fluid solution provided by the lattice Boltzmann method is massively scalable, allowing extensive use of cutting edge parallel computing resources to accelerate this phase of the solution process. The meshfree method with distance fields allows for exact satisfaction of boundary conditions making it possible to exactly capture the effects of the fluid field on the solid structure.

Decoding complex flow phenomena via Dynamic Mode Decomposition

In this thesis, the viability of the Dynamic Mode Decomposition (DMD) as a technique to analyze and model complex dynamic real-world systems is presented. This method derives, directly from data, computationally efficient reduced-order models (ROMs) which can replace too onerous or unavailable high-fidelity physics-based models. Optimizations and extensions to the standard implementation of the methodology are proposed, investigating diverse case studies related to the decoding of complex flow phenomena. The flexibility of this data-driven technique allows its application to high-fidelity fluid dynamics simulations, as well as time series of real systems observations. The resulting ROMs are tested against two tasks: (i) reduction of the storage requirements of high-fidelity simulations or observations; (ii) interpolation and extrapolation of missing data. The capabilities of DMD can also be exploited to alleviate the cost of onerous studies that require many simulations, such as uncertainty quantification analysis, especially when dealing with complex high-dimensional systems. In this context, a novel approach to address parameter variability issues when modeling systems with space and time-variant response is proposed. Specifically, DMD is merged with another model-reduction technique, namely the Polynomial Chaos Expansion, for uncertainty quantification purposes. Useful guidelines for DMD deployment result from the study, together with the demonstration of its potential to ease diagnosis and scenario analysis when complex flow processes are involved.

Mineral paragenesis and fluid inclusion constraints on the Torre di Rio skarn, Island of Elba

The historical iron ore deposits of eastern Elba held great importance for the region and were its primary source of iron. The Torre di Rio skarn, despite its easily accessible outcrop and vicinity to the larger Rio Marina deposit, was never properly characterized. The results of petrographic and microthermometric study presented in this work provide new constraints on the Torre di Rio skarn. Mineral assemblage of ilvaite, calcite, quartz, iron oxides and sulphides combined with textural evidence indicate that Torre di Rio skarn does not fit into classical skarn model. The complex paragenetic sequence and overlapping of skarn and ore mineralogy is result of fast formation at relatively low temperatures evidenced by the silicon enrichment and pervasive nature of limonite alteration. Hematite-magnetite textural relationship points to boundary conditions of the ore fluid in terms of oxygen fugacity. Eutectic temperatures range from -16 to -33 °C indicating complex fluids. Calculated salinities range from 1.4 to 17.4 wt% NaCleq suggesting multiple fluids of different compositions. Total homogenization temperatures vary from 330 °C to 150 °C with both homogeneously and heterogeneously trapped FIAs. Ore deposition is concentrated where skarn formation was controlled primarily by phase separation during boiling. Calculated fluid pressure at boiling suggest shallow formation depth of a few hundred meters and constrains maximum temperature of ore deposition to c. 260 °C. This work suggest that relatively low salinities of fluid inclusions could indicate dominant marine origin of the hydrothermal fluids that were activated by the Porto Azzurro pluton emplacement and that scavenged Fe from sedimentary host rocks. During boiling at shallow depths and decreasing iron solubility, these fluids started precipitating Fe-minerals at Torre di Rio mineralization. Mixing with batches of more saline fluids at around 236 °C increased salinity abruptly and marked the end of ore deposition.

Albumin Is Synthesized In Epididymis And Aggregates In A High Molecular Mass Glycoprotein Complex Involved In Sperm-egg Fertilization.

The epididymis has an important role in the maturation of sperm for fertilization, but little is known about the epididymal molecules involved in sperm modifications during this process. We have previously described the expression pattern for an antigen in epididymal epithelial cells that reacts with the monoclonal antibody (mAb) TRA 54. Immunohistochemical and immunoblotting analyses suggest that the epitope of the epididymal antigen probably involves a sugar moiety that is released into the epididymal lumen in an androgen-dependent manner and subsequently binds to luminal sperm. Using column chromatography, SDS-PAGE with in situ digestion and mass spectrometry, we have identified the protein recognized by mAb TRA 54 in mouse epididymal epithelial cells. The ∼65 kDa protein is part of a high molecular mass complex (∼260 kDa) that is also present in the sperm acrosomal vesicle and is completely released after the acrosomal reaction. The amino acid sequence of the protein corresponded to that of albumin. Immunoprecipitates with anti-albumin antibody contained the antigen recognized by mAb TRA 54, indicating that the epididymal molecule recognized by mAb TRA 54 is albumin. RT-PCR detected albumin mRNA in the epididymis and fertilization assays in vitro showed that the glycoprotein complex containing albumin was involved in the ability of sperm to recognize and penetrate the egg zona pellucida. Together, these results indicate that epididymal-derived albumin participates in the formation of a high molecular mass glycoprotein complex that has an important role in egg fertilization.

Effects Of Partial Inhibition Of Respiratory Complex I On H2o 2 Production By Isolated Brain Mitochondria In Different Respiratory States.

The aim of this work was to characterize the effects of partial inhibition of respiratory complex I by rotenone on H2O2 production by isolated rat brain mitochondria in different respiratory states. Flow cytometric analysis of membrane potential in isolated mitochondria indicated that rotenone leads to uniform respiratory inhibition when added to a suspension of mitochondria. When mitochondria were incubated in the presence of a low concentration of rotenone (10 nm) and NADH-linked substrates, oxygen consumption was reduced from 45.9 ± 1.0 to 26.4 ± 2.6 nmol O2 mg(-1) min(-1) and from 7.8 ± 0.3 to 6.3 ± 0.3 nmol O2 mg(-1) min(-1) in respiratory states 3 (ADP-stimulated respiration) and 4 (resting respiration), respectively. Under these conditions, mitochondrial H2O2 production was stimulated from 12.2 ± 1.1 to 21.0 ± 1.2 pmol H2O2 mg(-1) min(-1) and 56.5 ± 4.7 to 95.0 ± 11.1 pmol H2O2 mg(-1) min(-1) in respiratory states 3 and 4, respectively. Similar results were observed when comparing mitochondrial preparations enriched with synaptic or nonsynaptic mitochondria or when 1-methyl-4-phenylpyridinium ion (MPP(+)) was used as a respiratory complex I inhibitor. Rotenone-stimulated H2O2 production in respiratory states 3 and 4 was associated with a high reduction state of endogenous nicotinamide nucleotides. In succinate-supported mitochondrial respiration, where most of the mitochondrial H2O2 production relies on electron backflow from complex II to complex I, low rotenone concentrations inhibited H2O2 production. Rotenone had no effect on mitochondrial elimination of micromolar concentrations of H2O2. The present results support the conclusion that partial complex I inhibition may result in mitochondrial energy crisis and oxidative stress, the former being predominant under oxidative phosphorylation and the latter under resting respiration conditions.