Physical investigations of small particles : (I) aerosol particle charging and flux enhancement and (II) whispering gallery mode sensing

Part I

Particles are a key feature of planetary atmospheres. On Earth they represent the greatest source of uncertainty in the global energy budget. This uncertainty can be addressed by making more measurement, by improving the theoretical analysis of measurements, and by better modeling basic particle nucleation and initial particle growth within an atmosphere. This work will focus on the latter two methods of improvement.

Uncertainty in measurements is largely due to particle charging. Accurate descriptions of particle charging are challenging because one deals with particles in a gas as opposed to a vacuum, so different length scales come into play. Previous studies have considered the effects of transition between the continuum and kinetic regime and the effects of two and three body interactions within the kinetic regime. These studies, however, use questionable assumptions about the charging process which resulted in skewed observations, and bias in the proposed dynamics of aerosol particles. These assumptions affect both the ions and particles in the system. Ions are assumed to be point monopoles that have a single characteristic speed rather than follow a distribution. Particles are assumed to be perfect conductors that have up to five elementary charges on them. The effects of three body interaction, ion-molecule-particle, are also overestimated. By revising this theory so that the basic physical attributes of both ions and particles and their interactions are better represented, we are able to make more accurate predictions of particle charging in both the kinetic and continuum regimes.

The same revised theory that was used above to model ion charging can also be applied to the flux of neutral vapor phase molecules to a particle or initial cluster. Using these results we can model the vapor flux to a neutral or charged particle due to diffusion and electromagnetic interactions. In many classical theories currently applied to these models, the finite size of the molecule and the electromagnetic interaction between the molecule and particle, especially for the neutral particle case, are completely ignored, or, as is often the case for a permanent dipole vapor species, strongly underestimated. Comparing our model to these classical models we determine an “enhancement factor” to characterize how important the addition of these physical parameters and processes is to the understanding of particle nucleation and growth.

Part II

Whispering gallery mode (WGM) optical biosensors are capable of extraordinarily sensitive specific and non-specific detection of species suspended in a gas or fluid. Recent experimental results suggest that these devices may attain single-molecule sensitivity to protein solutions in the form of stepwise shifts in their resonance wavelength, \lambda_{R}, but present sensor models predict much smaller steps than were reported. This study examines the physical interaction between a WGM sensor and a molecule adsorbed to its surface, exploring assumptions made in previous efforts to model WGM sensor behavior, and describing computational schemes that model the experiments for which single protein sensitivity was reported. The resulting model is used to simulate sensor performance, within constraints imposed by the limited material property data. On this basis, we conclude that nonlinear optical effects would be needed to attain the reported sensitivity, and that, in the experiments for which extreme sensitivity was reported, a bound protein experiences optical energy fluxes too high for such effects to be ignored.

Phase-dependent properties for absorption and dispersion in a closed equispaced three-level ladder system

It is shown that in a closed equispaced three-level ladder system, by controlling the relative phase of two applied coherent fields, the conversion from absorption with inversion to lasing without inversion (LWI) can be realized; a large index of the refraction with zero absorption can be gotten; considerable increasing of the spectrum region and value of the LWI gain can be achieved. Our study also reveals that the incoherent pumping will produce a remarkable effect oil the phase-dependent properties of the system. Modifying value of the incoherent pumping can change the property of the system from absorption to amplification and enhance significantly LWI gain. If the incoherent pumping is absent, we cannot get any gain for any value of the relative phase. (c) 2007 Elsevier GmbH. All rights reserved.

Propagation of a few-cycle laser pulse in a V-type three-level system

Propagation of a few-cycle laser pulse in a V-type three-level system (fine structure levels of rubidium) is investigated numerically. The full three-level Maxwell-Bloch equations without the rotating wave approximation and the standing slowly varying envelope approximation are solved by using a finite-difference time-domain method. It is shown that, when the usual unequal oscillator strengths are considered, self-induced transparency cannot be recovered and higher spectral components can be produced even for small-area pulses. (c) 2005 Pleiades Publishing, Inc.

Test particle motion in a Lorentz gas

This is a two-part thesis concerning the motion of a test particle in a bath. In part one we use an expansion of the operator PLe^it(1-P)LLP to shape the Zwanzig equation into a generalized Fokker-Planck equation which involves a diffusion tensor depending on the test particle's momentum and the time.

In part two the resultant equation is studied in some detail for the case of test particle motion in a weakly coupled Lorentz Gas. The diffusion tensor for this system is considered. Some of its properties are calculated; it is computed explicitly for the case of a Gaussian potential of interaction.

The equation for the test particle distribution function can be put into the form of an inhomogeneous Schroedinger equation. The term corresponding to the potential energy in the Schroedinger equation is considered. Its structure is studied, and some of its simplest features are used to find the Green's function in the limiting situations of low density and long time.

Energy spectrum of ground state and excitation spectrum of quasi-particle for hard-core boson in optical lattices

We investigate the energy spectrum of ground state and quasi-particle excitation spectrum of hard-core bosons, which behave very much like spinless noninteracting fermions, in optical lattices by means of the perturbation expansion and Bogoliubov approach. The results show that the energy spectrum has a single band structure, and the energy is lower near zero momentum; the excitation spectrum gives corresponding energy gap, and the system is in Mott-insulating state at Tonks limit. The analytic result of energy spectrum is in good agreement with that calculated in terms of Green's function at strong correlation limit.

Dynamic characterization of micro-particle systems

Ordered granular systems have been a subject of active research for decades. Due to their rich dynamic response and nonlinearity, ordered granular systems have been suggested for several applications, such as solitary wave focusing, acoustic signals manipulation, and vibration absorption. Most of the fundamental research performed on ordered granular systems has focused on macro-scale examples. However, most engineering applications require these systems to operate at much smaller scales. Very little is known about the response of micro-scale granular systems, primarily because of the difficulties in realizing reliable and quantitative experiments, which originate from the discrete nature of granular materials and their highly nonlinear inter-particle contact forces.

In this work, we investigate the physics of ordered micro-granular systems by designing an innovative experimental platform that allows us to assemble, excite, and characterize ordered micro-granular systems. This new experimental platform employs a laser system to deliver impulses with controlled momentum and incorporates non-contact measurement apparatuses to detect the particles’ displacement and velocity. We demonstrated the capability of the laser system to excite systems of dry (stainless steel particles of radius 150 micrometers) and wet (silica particles of radius 3.69 micrometers, immersed in fluid) micro-particles, after which we analyzed the stress propagation through these systems.

We derived the equations of motion governing the dynamic response of dry and wet particles on a substrate, which we then validated in experiments. We then measured the losses in these systems and characterized the collision and friction between two micro-particles. We studied wave propagation in one-dimensional dry chains of micro-particles as well as in two-dimensional colloidal systems immersed in fluid. We investigated the influence of defects to wave propagation in the one-dimensional systems. Finally, we characterized the wave-attenuation and its relation to the viscosity of the surrounding fluid and performed computer simulations to establish a model that captures the observed response.

The findings of the study offer the first systematic experimental and numerical analysis of wave propagation through ordered systems of micro-particles. The experimental system designed in this work provides the necessary tools for further fundamental studies of wave propagation in both granular and colloidal systems.

On the transport properties of fluid-particle flow

The hydrodynamic forces acting on a solid particle in a viscous, incompressible fluid medium at low Reynolds number flow is investigated mathematically as a prerequisite to the understanding of transport processes in two-phase flow involving solid particles and fluid. Viscous interaction between a small number of spherical particles and continuous solid boundaries as well as fluid interface are analyzed under a “point-force” approximation. Non-spherical and elastic spherical particles in a simple shear flow area are then considered. Non-steady motion of a spherical particle is briefly touched upon to illustrate the transient effect of particle motion.

A macroscopic continuum description of particle-fluid flow is formulated in terms of spatial averages yielding a set of particle continuum and bulk fluid equations. Phenomenological formulas describing the transport processes in a fluid medium are extended to cases where the volume concentration of solid particles is sufficiently high to exert an important influence. Hydrodynamic forces acting on a spherical solid particle in such a system, e.g. drag, torque, rotational coupling force, and viscous collision force between streams of different sized particles moving relative to each other are obtained. Phenomenological constants, such as the shear viscosity coefficient, and the diffusion coefficient of the bulk fluid, are found as a function of the material properties of the constituents of the two-phase system and the volume concentration of solid. For transient heat conduction phenomena, it is found that the introduction of a complex conductivity for the bulk fluid permits a simple mathematical description of this otherwise complicated process. The rate of heat transfer between particle continuum and bulk fluid is also investigated by means of an Oseen-type approximation to the energy equation.

I. The attenuation and dispersion of sound in a condensing medium. II. The flow of a gas-particle mixture downstream of a normal shock in a nozzle

I. The attenuation of sound due to particles suspended in a gas was first calculated by Sewell and later by Epstein in their classical works on the propagation of sound in a two-phase medium. In their work, and in more recent works which include calculations of sound dispersion, the calculations were made for systems in which there was no mass transfer between the two phases. In the present work, mass transfer between phases is included in the calculations.

The attenuation and dispersion of sound in a two-phase condensing medium are calculated as functions of frequency. The medium in which the sound propagates consists of a gaseous phase, a mixture of inert gas and condensable vapor, which contains condensable liquid droplets. The droplets, which interact with the gaseous phase through the interchange of momentum, energy, and mass (through evaporation and condensation), are treated from the continuum viewpoint. Limiting cases, for flow either frozen or in equilibrium with respect to the various exchange processes, help demonstrate the effects of mass transfer between phases. Included in the calculation is the effect of thermal relaxation within droplets. Pressure relaxation between the two phases is examined, but is not included as a contributing factor because it is of interest only at much higher frequencies than the other relaxation processes. The results for a system typical of sodium droplets in sodium vapor are compared to calculations in which there is no mass exchange between phases. It is found that the maximum attenuation is about 25 per cent greater and occurs at about one-half the frequency for the case which includes mass transfer, and that the dispersion at low frequencies is about 35 per cent greater. Results for different values of latent heat are compared.

II. In the flow of a gas-particle mixture through a nozzle, a normal shock may exist in the diverging section of the nozzle. In Marble’s calculation for a shock in a constant area duct, the shock was described as a usual gas-dynamic shock followed by a relaxation zone in which the gas and particles return to equilibrium. The thickness of this zone, which is the total shock thickness in the gas-particle mixture, is of the order of the relaxation distance for a particle in the gas. In a nozzle, the area may change significantly over this relaxation zone so that the solution for a constant area duct is no longer adequate to describe the flow. In the present work, an asymptotic solution, which accounts for the area change, is obtained for the flow of a gas-particle mixture downstream of the shock in a nozzle, under the assumption of small slip between the particles and gas. This amounts to the assumption that the shock thickness is small compared with the length of the nozzle. The shock solution, valid in the region near the shock, is matched to the well known small-slip solution, which is valid in the flow downstream of the shock, to obtain a composite solution valid for the entire flow region. The solution is applied to a conical nozzle. A discussion of methods of finding the location of a shock in a nozzle is included.

I. Quantum mechanics of one-dimensional two-particle models. II. A quantum mechanical treatment of inelastic collisions

Part I

Solutions of Schrödinger’s equation for system of two particles bound in various stationary one-dimensional potential wells and repelling each other with a Coulomb force are obtained by the method of finite differences. The general properties of such systems are worked out in detail for the case of two electrons in an infinite square well. For small well widths (1-10 a.u.) the energy levels lie above those of the noninteresting particle model by as much as a factor of 4, although excitation energies are only half again as great. The analytical form of the solutions is obtained and it is shown that every eigenstate is doubly degenerate due to the “pathological” nature of the one-dimensional Coulomb potential. This degeneracy is verified numerically by the finite-difference method. The properties of the square-well system are compared with those of the free-electron and hard-sphere models; perturbation and variational treatments are also carried out using the hard-sphere Hamiltonian as a zeroth-order approximation. The lowest several finite-difference eigenvalues converge from below with decreasing mesh size to energies below those of the “best” linear variational function consisting of hard-sphere eigenfunctions. The finite-difference solutions in general yield expectation values and matrix elements as accurate as those obtained using the “best” variational function.

The system of two electrons in a parabolic well is also treated by finite differences. In this system it is possible to separate the center-of-mass motion and hence to effect a considerable numerical simplification. It is shown that the pathological one-dimensional Coulomb potential gives rise to doubly degenerate eigenstates for the parabolic well in exactly the same manner as for the infinite square well.

Part II

A general method of treating inelastic collisions quantum mechanically is developed and applied to several one-dimensional models. The formalism is first developed for nonreactive “vibrational” excitations of a bound system by an incident free particle. It is then extended to treat simple exchange reactions of the form A + BC →AB + C. The method consists essentially of finding a set of linearly independent solutions of the Schrödinger equation such that each solution of the set satisfies a distinct, yet arbitrary boundary condition specified in the asymptotic region. These linearly independent solutions are then combined to form a total scattering wavefunction having the correct asymptotic form. The method of finite differences is used to determine the linearly independent functions.

The theory is applied to the impulsive collision of a free particle with a particle bound in (1) an infinite square well and (2) a parabolic well. Calculated transition probabilities agree well with previously obtained values.

Several models for the exchange reaction involving three identical particles are also treated: (1) infinite-square-well potential surface, in which all three particles interact as hard spheres and each two-particle subsystem (i.e. BC and AB) is bound by an attractive infinite-square-well potential; (2) truncated parabolic potential surface, in which the two-particle subsystems are bound by a harmonic oscillator potential which becomes infinite for interparticle separations greater than a certain value; (3) parabolic (untruncated) surface. Although there are no published values with which to compare our reaction probabilities, several independent checks on internal consistency indicate that the results are reliable.

Evaluation of the In Situ Ichthyoplankton Imaging System (ISIIS): comparison with the traditional (bongo net) sampler

Plankton and larval fish sampling programs often are limited by a balance between sampling frequency (for precision) and costs. Advancements in sampling techniques hold the potential to add considerable efficiency and, therefore, add sampling frequency to improve precision. We compare a newly developed plankton imaging system, In Situ Ichthyoplankton Imaging System (ISIIS), with a bongo sampler, which is a traditional plankton sampling gear developed in the 1960s. Comparative sampling was conducted along 2 transects ~30–40 km long. Over 2 days, we completed 36 ISIIS tow-yo undulations and 11 bongo oblique tows, each from the surface to within 10 m of the seafloor. Overall, the 2 gears detected comparable numbers of larval fishes, representing similar taxonomic compositions, although larvae captured with the bongo were capable of being identified to lower taxonomic levels, especially larvae in the small (<5 mm), preflexion stages. Size distributions of the sampled larval fishes differed considerably between these 2 sampling methods, with the size range and mean size of larval fishes larger with ISIIS than with the bongo sampler. The high frequency and fine spatial scale of ISIIS allow it to add considerable sampling precision (i.e., more vertical sections) to plankton surveys. Improvements in the ISIIS technology (including greater depth of field and image resolution) should also increase taxonomic resolution and decrease processing time. When coupled with appropriate net sampling (for the purpose of collecting and verifying the identification of biological samples), the use of ISIIS could improve overall survey design and simultaneously provide detailed, process-oriented information for fisheries scientists and oceanographers.

Mechanisms of erosion of unfilled elastomers by solid particle impact

The mechanisms of material removal were studied during the erosion of two unfilled elastomers (natural rubber and epoxidised natural rubber). The effects of impact velocity and of lubrication by silicone oil were investigated. The development of surface features due to single impacts and during the early stages of erosion was followed by scanning electron microscopy. The basic material removal mechanism at impact angles of both 30° and 90° involves the formation and growth of fine fatigue cracks under the tensile surface stresses caused by impact. No damage was observed after single impacts; it was found that many successive impacts are necessary for material removal. It was found that the erosion rate has a very strong dependance on impact velocity above about 50 ms-1.

Effects of supplemental feed and fertilizer on growth and survival of Macrobrachium rosenbergii (de Man) post larvae in pond nursery system

The study was conducted to compare the performance of different nursing practices of giant freshwater prawn (Macrobrachium rosenbergii) post-larvae (PL). Three treatments such as only fertilizers (T1), fertilizers with 5% supplementary feed (local feed) (T2), and 10% commercial feed (T3) were applied in the nursing system of prawn PLs in earthen pond. An earthen pond (315m²) was divided into nine equal small ponds by fine meshed nylon nets. Feeds were used once daily on a tray placed near the pond bottom. There was a significant difference (p<0.05) in some water quality parameters like pH and total alkalinity, but all measured water quality parameters viz. water temperature, transparency, dissolved oxygen and ammonia-nitrogen were within the acceptable range for nursing of prawn PL. The results showed that the mean final lengths of prawn post-larvae were 6.3±0.07 cm, 7.12±0.22 cm and 8.17±0.16 cm in T1, T2 and T3, respectively. There were significant difference (p<0.05) in mean final length of prawn PL among the treatments. Significantly higher (p<0.05) average daily weight gain was observed in T3 (0.071 ±0.007 g) than in T2 (0.052±0.006 g) and T1 (0.031 ±0.002 g). The specific growth rate (SGR) of T3 (8.81±0.26) was found significantly higher (p<0.05) than T2 (8.35±0.22) and T1 (7.42±0.11). Survival rate (%) was also significantly higher (p<0.05) in T3 (66.24±1.58) than in T2 (60.52±1.64) and T1 (53.86±2.71). Therefore, it may be concluded that the growth and survival in prawn nursery was better in commercial feed than only fertilizers and fertilizers with local feeds.

UV durable colour pigment doped SmA liquid crystal composites for outdoor trans-reflective bi-stable displays

High brightness trans-reflective bi-stable displays based on smectic A (SmA) liquid crystals (LCs) can have nearly perfect transparency in the clear state and very high reflection in the scattered state. Because the LC material in use is stable under UV radiation, this kind of displays can stand for strong day-light and therefore be ideal for outdoor applications from e-books to public signage and advertisement. However, the colour application has been limited because the traditional colourants in use are conventional dyes which are lack of UV stability and that their colours are easily photo bleached. Here we present a colour SmA display demonstrator using pigments as colourant. Mixing pigments with SmA LCs and maintain the desirable optical switching performance is not straightforward. We show here how it can be done, including how to obtain fine sized pigment nano-particles, the effects of particle size and size distribution on the display performance. Our optimized pigments/SmA compositions can be driven by a low frequency waveform (∼101Hz) to a scattered state to exhibit colour while by a high frequency waveform (∼103Hz) to a cleared state showing no colour. Finally, we will present its excellent UV life-time (at least >7.2 years) in comparison with that of dye composition (∼2.4 years). The complex interaction of pigment nano-particles with LC molecules and the resulting effects on the LC electro-optical performances are still to be fully understood. We hope this work will not only demonstrate a new and practical approach for outdoor reflective colour displays but also provide a new material system for fundamental liquid crystal colloid research work. © 2012 SPIE.

An integrated 2-D model of a lithium ion battery: The effect of material parameters and morphology on storage particle stress

An integrated 2-D model of a lithium ion battery is developed to study the mechanical stress in storage particles as a function of material properties. A previously developed coupled stress-diffusion model for storage particles is implemented in 2-D and integrated into a complete battery system. The effect of morphology on the stress and lithium concentration is studied for the case of extraction of lithium in terms of previously developed non-dimensional parameters. These non-dimensional parameters include the material properties of the storage particles in the system, among other variables. We examine particles functioning in isolation as well as in closely-packed systems. Our results show that the particle distance from the separator, in combination with the material properties of the particle, is critical in predicting the stress generated within the particle. © 2012 Springer-Verlag.

Field programmable gate array based predictive control system for spacecraft rendezvous in elliptical orbits

Copyright © 2014 John Wiley & Sons, Ltd. Copyright © 2014 John Wiley & Sons, Ltd. Summary A field programmable gate array (FPGA) based model predictive controller for two phases of spacecraft rendezvous is presented. Linear time-varying prediction models are used to accommodate elliptical orbits, and a variable prediction horizon is used to facilitate finite time completion of the longer range manoeuvres, whilst a fixed and receding prediction horizon is used for fine-grained tracking at close range. The resulting constrained optimisation problems are solved using a primal-dual interior point algorithm. The majority of the computational demand is in solving a system of simultaneous linear equations at each iteration of this algorithm. To accelerate these operations, a custom circuit is implemented, using a combination of Mathworks HDL Coder and Xilinx System Generator for DSP, and used as a peripheral to a MicroBlaze soft-core processor on the FPGA, on which the remainder of the system is implemented. Certain logic that can be hard-coded for fixed sized problems is implemented to be configurable online, in order to accommodate the varying problem sizes associated with the variable prediction horizon. The system is demonstrated in closed-loop by linking the FPGA with a simulation of the spacecraft dynamics running in Simulink on a PC, using Ethernet. Timing comparisons indicate that the custom implementation is substantially faster than pure embedded software-based interior point methods running on the same MicroBlaze and could be competitive with a pure custom hardware implementation.