Exact solutions and transformation properties of nonlinear partial differential equations from general relativity

Various families of exact solutions to the Einstein and Einstein-Maxwell field equations of General Relativity are treated for situations of sufficient symmetry that only two independent variables arise. The mathematical problem then reduces to consideration of sets of two coupled nonlinear differential equations.

The physical situations in which such equations arise include: a) the external gravitational field of an axisymmetric, uncharged steadily rotating body, b) cylindrical gravitational waves with two degrees of freedom, c) colliding plane gravitational waves, d) the external gravitational and electromagnetic fields of a static, charged axisymmetric body, and e) colliding plane electromagnetic and gravitational waves. Through the introduction of suitable potentials and coordinate transformations, a formalism is presented which treats all these problems simultaneously. These transformations and potentials may be used to generate new solutions to the Einstein-Maxwell equations from solutions to the vacuum Einstein equations, and vice-versa.

The calculus of differential forms is used as a tool for generation of similarity solutions and generalized similarity solutions. It is further used to find the invariance group of the equations; this in turn leads to various finite transformations that give new, physically distinct solutions from old. Some of the above results are then generalized to the case of three independent variables.

Some problems in nonlinear elasticity

Two separate problems are discussed: axisymmetric equilibrium configurations of a circular membrane under pressure and subject to thrust along its edge, and the buckling of a circular cylindrical shell.

An ordinary differential equation governing the circular membrane is imbedded in a family of n-dimensional nonlinear equations. Phase plane methods are used to examine the number of solutions corresponding to a parameter which generalizes the thrust, as well as other parameters determining the shape of the nonlinearity and the undeformed shape of the membrane. It is found that in any number of dimensions there exists a value of the generalized thrust for which a countable infinity of solutions exist if some of the remaining parameters are made sufficiently large. Criteria describing the number of solutions in other cases are also given.

Donnell-type equations are used to model a circular cylindrical shell. The static problem of bifurcation of buckled modes from Poisson expansion is analyzed using an iteration scheme and pertubation methods. Analysis shows that although buckling loads are usually simple eigenvalues, they may have arbitrarily large but finite multiplicity when the ratio of the shell's length and circumference is rational. A numerical study of the critical buckling load for simple eigenvalues indicates that the number of waves along the axis of the deformed shell is roughly proportional to the length of the shell, suggesting the possibility of a "characteristic length." Further numerical work indicates that initial post-buckling curves are typically steep, although the load may increase or decrease. It is shown that either a sheet of solutions or two distinct branches bifurcate from a double eigenvalue. Furthermore, a shell may be subject to a uniform torque, even though one is not prescribed at the ends of the shell, through the interaction of two modes with the same number of circumferential waves. Finally, multiple time scale techniques are used to study the dynamic buckling of a rectangular plate as well as a circular cylindrical shell; transition to a new steady state amplitude determined by the nonlinearity is shown. The importance of damping in determining equilibrium configurations independent of initial conditions is illustrated.

Some approximate solutions of dynamic problems in the linear theory of thin elastic shells

Some aspects of wave propagation in thin elastic shells are considered. The governing equations are derived by a method which makes their relationship to the exact equations of linear elasticity quite clear. Finite wave propagation speeds are ensured by the inclusion of the appropriate physical effects.

The problem of a constant pressure front moving with constant velocity along a semi-infinite circular cylindrical shell is studied. The behavior of the solution immediately under the leading wave is found, as well as the short time solution behind the characteristic wavefronts. The main long time disturbance is found to travel with the velocity of very long longitudinal waves in a bar and an expression for this part of the solution is given.

When a constant moment is applied to the lip of an open spherical shell, there is an interesting effect due to the focusing of the waves. This phenomenon is studied and an expression is derived for the wavefront behavior for the first passage of the leading wave and its first reflection.

For the two problems mentioned, the method used involves reducing the governing partial differential equations to ordinary differential equations by means of a Laplace transform in time. The information sought is then extracted by doing the appropriate asymptotic expansion with the Laplace variable as parameter.

A biochemical and structural characterization of Drosophila neuroglian

A number of cell-cell interactions in the nervous system are mediated by immunoglobulin gene superfamily members. For example, neuroglian, a homophilic neural cell adhesion molecule in Drosophila, has an extracellular portion comprising six C- 2 type immunoglobulin-like domains followed by five fibronectin type III (FnIII) repeats. Neuroglian shares this domain organization and significant sequence identity with Ll, a murine neural adhesion molecule that could be a functional homologue. Here I report the crystal structure of a proteolytic fragment containing the first two FnIII repeats of neuroglian (NgFn 1,2) at 2.0Å. The interpretation of photomicrographs of rotary shadowed Ng, the entire extracellular portion of neuroglian, and NgFnl-5, the five neuroglian Fn III domains, is also discussed.

The structure of NgFn 1,2 consists of two roughly cylindrical β-barrel structural motifs arranged in a head-to-tail fashion with the domains meeting at an angle of ~120, as defined by the cylinder axes. The folding topology of each domain is identical to that previously observed for single FnIII domains from tenascin and fibronectin. The domains of NgFn1,2 are related by an approximate two fold screw axis that is nearly parallel to the longest dimension of the fragment. Assuming this relative orientation is a general property of tandem FnIII repeats, the multiple tandem FnIII domains in neuroglian and other proteins are modeled as thin straight rods with two domain zig-zag repeats. When combined with the dimensions of pairs of tandem immunoglobulin-like domains from CD4 and CD2, this model suggests that neuroglian is a long narrow molecule (20 - 30 Å in diameter) that extends up to 370Å from the cell surface.

In photomicrographs, rotary shadowed Ng and NgFn1-5 appear to be highly flexible rod-like molecules. NgFn 1-5 is observed to bend in at least two positions and has a mean total length consistent with models generated from the NgFn1,2 structure. Ng molecules have up to four bends and a mean total length of 392 Å, consistent with a head-to-tail packing of neuroglian's C2-type domains.

Three-dimensional models of mantle convection : influence of phase transitions and temperature-dependent viscosity

Two of the most important questions in mantle dynamics are investigated in three separate studies: the influence of phase transitions (studies 1 and 2), and the influence of temperature-dependent viscosity (study 3).

(1) Numerical modeling of mantle convection in a three-dimensional spherical shell incorporating the two major mantle phase transitions reveals an inherently three-dimensional flow pattern characterized by accumulation of cold downwellings above the 670 km discontinuity, and cylindrical 'avalanches' of upper mantle material into the lower mantle. The exothermic phase transition at 400 km depth reduces the degree of layering. A region of strongly-depressed temperature occurs at the base of the mantle. The temperature field is strongly modulated by this partial layering, both locally and in globally-averaged diagnostics. Flow penetration is strongly wavelength-dependent, with easy penetration at long wavelengths but strong inhibition at short wavelengths. The amplitude of the geoid is not significantly affected.

(2) Using a simple criterion for the deflection of an upwelling or downwelling by an endothermic phase transition, the scaling of the critical phase buoyancy parameter with the important lengthscales is obtained. The derived trends match those observed in numerical simulations, i.e., deflection is enhanced by (a) shorter wavelengths, (b) narrower up/downwellings (c) internal heating and (d) narrower phase loops.

(3) A systematic investigation into the effects of temperature-dependent viscosity on mantle convection has been performed in three-dimensional Cartesian geometry, with a factor of 1000-2500 viscosity variation, and Rayleigh numbers of 10^5-10^7. Enormous differences in model behavior are found, depending on the details of rheology, heating mode, compressibility and boundary conditions. Stress-free boundaries, compressibility, and temperature-dependent viscosity all favor long-wavelength flows, even in internally heated cases. However, small cells are obtained with some parameter combinations. Downwelling plumes and upwelling sheets are possible when viscosity is dependent solely on temperature. Viscous dissipation becomes important with temperature-dependent viscosity.

The sensitivity of mantle flow and structure to these various complexities illustrates the importance of performing mantle convection calculations with rheological and thermodynamic properties matching as closely as possible those of the Earth.

Structure of the earth's core and lowermost mantle from seismic PKP waves

This thesis addresses the fine structure, both radial and lateral, of compressional wave velocity and attenuation of the Earth's core and the lowermost mantle using waveforms, differential travel times and amplitudes of PKP waves, which penetrate the Earth's core.

The structure near the inner core boundary (ICB) is studied by analyzing waveforms of a regional sample. The waveform modeling approach is demonstrated to be an effective tool for constrainning the ICB structure. The best model features a sharp velocity jump of 0.78km/s at the ICB and a low velocity gradient at the lowermost outer core (indicating possible inhomogeneity) and high attenuation at the top of the inner core.

A spherically symmetric P-wave model of the core, is proposed from PKP differential times, waveforms and amplitudes. The ICB remains sharp with a velocity jump of 0. 78km/ s. A very low velocity gradient at the base of the fluid core is demonstrated to be a robust feature, indicating inhomogeneity is practically inevitable. The model also indicates that the attenuation in the inner core decreases with depth. The velocity at D" is smaller than PREM.

The inner core is confirmed to be very anisotropic, possessing a cylindrical symmetry around the Earth spin axis with the N-S direction 3% faster than the E-W direction. All of the N-S rays through the inner core were found to be faster than the E-W rays by 1.5 to 3.5s. Exhaustive data selection and efforts in insolating contributions from the region above ensure that this is an inner core feature.

The anisotropy at the very top of the inner core is found to be distinctly different from the deeper part. The top 60km of the inner core is not anisotropic. From 60km to 150km, there appears to be a transition from isotropy to anisotropy.

PKP differential travel times are used to study the P velocity structure in D". Systematic regional variations of up to 2s in AB-DF times were observed, attributed primarily to heterogeneities in the lower 500km of the mantle. However, direct comparisons with tomographic models are not successful.

Coupled effects of mechanics, geometry, and chemistry on bio-membrane behavior

Lipid bilayer membranes are models for cell membranes--the structure that helps regulate cell function. Cell membranes are heterogeneous, and the coupling between composition and shape gives rise to complex behaviors that are important to regulation. This thesis seeks to systematically build and analyze complete models to understand the behavior of multi-component membranes.

We propose a model and use it to derive the equilibrium and stability conditions for a general class of closed multi-component biological membranes. Our analysis shows that the critical modes of these membranes have high frequencies, unlike single-component vesicles, and their stability depends on system size, unlike in systems undergoing spinodal decomposition in flat space. An important implication is that small perturbations may nucleate localized but very large deformations. We compare these results with experimental observations.

We also study open membranes to gain insight into long tubular membranes that arise for example in nerve cells. We derive a complete system of equations for open membranes by using the principle of virtual work. Our linear stability analysis predicts that the tubular membranes tend to have coiling shapes if the tension is small, cylindrical shapes if the tension is moderate, and beading shapes if the tension is large. This is consistent with experimental observations reported in the literature in nerve fibers. Further, we provide numerical solutions to the fully nonlinear equilibrium equations in some problems, and show that the observed mode shapes are consistent with those suggested by linear stability. Our work also proves that beadings of nerve fibers can appear purely as a mechanical response of the membrane.

Oxidative chemisorption of SO_2 on γ - Al_2O_3 and deposition of H_2-permselective SiO_2 films

The interaction of SO_2 with γ - Al_2O_3 and the deposition of H_2 permselective SiO_2 films have been investigated. The adsorption and oxidative adsorption of SO_2 on γ - Al_2O_3 have been examined at temperatures 500-700°C by Fourier transform infrared spectroscopy (FTIR) and thermogravimetric analysis (TGA). At temperatures above 500°C most of SO_2 adsorbed on the strong sites on alumina. The adsorbed SO_2 species was characterized by an IR band at 1065 cm^(-1). The equilibrium coverage and initial rate of adsorption decreased with temperature suggesting a two-step adsorption. When γ - Al_2O_3 was contacted with a mixture of SO_2 and O_2, adsorption of SO_2 and oxidation of the adsorbed SO_2 to a surface sulfate characterized by broad IR bands at 1070 cm^(-1), 1390 cm^(-1) took place. The results of a series of TGA experiments under different atmospheres strongly suggest that surface SO_2 and surface sulfate involve the same active sites such that SO_2 adsorption is inhibited by already formed sulfate. The results also indicate a broad range of site strengths.

The desorption of adsorbed SO_2 and the reductive desorption of oxidatively adsorbed SO_2 have been investigated by microreactor experiments and thermogravimetric analysis (TGA). Temperature programmed reduction (TPR) of adsorbed SO_2 showed that SO_2 was desorbed without significant reaction with H_2 when H_2 concentration was low while considerable reaction occurred when 100% H_2 was used. SO_2 adsorbed on the strong sites on alumina was reduced to sulfur and H_2S. The isothermal reduction experiments of oxidatively adsorbed SO_2 reveal that the rate of reduction is very slow below 550°C even with 100% H_2. The reduction product is mainly composed of SO_2. TPR experiments of oxidatively adsorbed SO_2 showed that H_2S arose from a sulfate strongly chemisorbed on the surface.

Films of amorphous SiO_2 were deposited within the walls of porous Vycor tubes by SiH_4 oxidation in an opposing reactants geometry : SiH_4 was passed inside the tube while O_2 was passed outside the tube. The two reactants diffused opposite to each other and reacted within a narrow front inside the tube wall to form a thin SiO_2 film. Once the pores were plugged the reactants could not reach each other and the reaction stopped. At 450°C and 0.1 and 0.33 atm of SiH_4 and O_2, the reaction was complete within 15 minutes. The thickness of the SiO_2 film was estimated to be about 0.1 µm. Measurements of H_2 and N_2 permeation rates showed that the SiO_2 film was highly selective to H_2 permeation. The H_2:N_2 flux at 450°C varied between 2000-3000.

Thin SiO_2 films were heat treated in different gas mixtures to determine their stability in functioning as high-temperature hydrogen-permselective membranes. The films were heat-treated at 450-700°C in dry N_2, dry O_2, N_2-H_2O, and O_2-H_2O mixtures. The permeation rates of H_2 and N_2 changed depending on the original conditions of film formation as well as on the heat treatment. Heating in dry N_2 slowly reduced the permeation rates of both H_2 and N_2. Heating in a N_2-H_2O atmosphere led to a steeper decline of H_2 permeability. But the permeation rate of N_2 increased or decreased according to whether the film deposition had been carried out in the absence or presence of H_2O vapor, respectively. Thermal treatment in O_2 caused rapid decline of the permeation rates of H_2 and N_2 in films that were deposited under dry conditions. The decline was moderate in films deposited under wet conditions.

Self-assembly of brush polymers

The unique structure and properties of brush polymers have led to increased interest in them within the scientific community. This thesis describes studies on the self-assembly of these brush polymers.

Chapter 2 describes a study on the rapid self-assembly of brush block copolymers into nanostructures with photonic bandgaps spanning the entire visible spectrum, from ultraviolet to near infrared. Linear relationships are observed between the peak wavelengths of reflection and polymer molecular weights. This work enables "bottom-up" fabrication of photonic crystals with application-tailored bandgaps, through synthetic control of the polymer molecular weight and the method of self-assembly.

Chapter 3 details the analysis of the self-assembly of symmetrical brush block copolymers in bulk and thin films. Highly ordered lamellae with domain spacing ranging from 20 to 240 nm are obtained by varying molecular weight of the backbone. The relationship between degree of polymerization and the domain spacing is reported, and evidence is provided for how rapidly the brush block copolymers self-assemble and achieve thermodynamic equilibrium.

Chapter 4 describes investigations into where morphology transitions take place as the volume fraction of each block is varied in asymmetrical brush block copolymers. Imaging techniques are used to observe a transition from lamellar to a cylindrical morphology as the volume fraction of one of the blocks exceeds 70%. It is also shown that the asymmetric brush block copolymers can be kinetically trapped into undulating lamellar structures by drop casting the samples.

Chapter 5 explores the capability of macromolecules to interdigitate into densely grafted molecular brush copolymers using stereocomplex formation as a driving force. The stereocomplex formation between complementary linear polymers and brush copolymers is demonstrated, while the stereocomplex formation between complementary brush copolymers is shown to be restricted.

A multi-scale approach to shaping carbon nanotube structures for hollow microneedles

The concept of a carbon nanotube microneedle array is explored in this thesis from multiple perspectives including microneedle fabrication, physical aspects of transdermal delivery, and in vivo transdermal drug delivery experiments. Starting with standard techniques in carbon nanotube (CNT) fabrication, including catalyst patterning and chemical vapor deposition, vertically-aligned carbon nanotubes are utilized as a scaffold to define the shape of the hollow microneedle. Passive, scalable techniques based on capillary action and unique photolithographic methods are utilized to produce a CNT-polymer composite microneedle. Specific examples of CNT-polyimide and CNT-epoxy microneedles are investigated. Further analysis of the transport properties of polymer resins reveals general requirements for applying arbitrary polymers to the fabrication process.

The bottom-up fabrication approach embodied by vertically-aligned carbon nanotubes allows for more direct construction of complex high-aspect ratio features than standard top-down fabrication approaches, making microneedles an ideal application for CNTs. However, current vertically-aligned CNT fabrication techniques only allow for the production of extruded geometries with a constant cross-sectional area, such as cylinders. To rectify this limitation, isotropic oxygen etching is introduced as a novel fabrication technique to create true 3D CNT geometry. Oxygen etching is utilized to create a conical geometry from a cylindrical CNT structure as well as create complex shape transformations in other CNT geometries.

CNT-polymer composite microneedles are anchored onto a common polymer base less than 50 µm thick, which allows for the microneedles to be incorporated into multiple drug delivery platforms, including modified hypodermic syringes and silicone skin patches. Cylindrical microneedles are fabricated with 100 µm outer diameter and height of 200-250 µm with a central cavity, or lumen, diameter of 30 µm to facilitate liquid drug flow. In vitro delivery experiments in swine skin demonstrate the ability of the microneedles to successfully penetrate the skin and deliver aqueous solutions.

An in vivo study was performed to assess the ability of the CNT-polymer microneedles to deliver drugs transdermally. CNT-polymer microneedles are attached to a hand actuated silicone skin patch that holds a liquid reservoir of drugs. Fentanyl, a potent analgesic, was administered to New Zealand White Rabbits through 3 routes of delivery: topical patch, CNT-polymer microneedles, and subcutaneous hypodermic injection. Results demonstrate that the CNT-polymer microneedles have a similar onset of action as the topical patch. CNT-polymer microneedles were also vetted as a painless delivery approach compared to hypodermic injection. Comparative analysis with contemporary microneedle designs demonstrates that the delivery achieved through CNT-polymer microneedles is akin to current hollow microneedle architectures. The inherent advantage of applying a bottom-up fabrication approach alongside similar delivery performance to contemporary microneedle designs demonstrates that the CNT-polymer composite microneedle is a viable architecture in the emerging field of painless transdermal delivery.

Radiation from a short electric dipole antenna in a hot uniaxial plasma

The effects of electron temperature on the radiation fields and the resistance of a short dipole antenna embedded in a uniaxial plasma have been studied. It is found that for ω < ω_p the antenna excites two waves, a slow wave and a fast wave. These waves propagate only within a cone whose axis is parallel to the biasing magnetostatic field B_o and whose semicone angle is slightly less than sin ^(-1) (ω/ω_p). In the case of ω > ω_p the antenna excites two separate modes of radiation. One of the modes is the electromagnetic mode, while the other mode is of hot plasma origin. A characteristic interference structure is noted in the angular distribution of the field. The far fields are evaluated by asymptotic methods, while the near fields are calculated numerically. The effects of antenna length ℓ, electron thermal speed, collisional and Landau damping on the near field patterns have been studied.

The input and the radiation resistances are calculated and are shown to remain finite for nonzero electron thermal velocities. The effect of Landau damping and the antenna length on the input and radiation resistances has been considered.

The radiation condition for solving Maxwell's equations is discussed and the phase and group velocities for propagation given. It is found that for ω < ω_p in the radial direction (cylindrical coordinates) the power flow is in the opposite direction to that of the phase propagation. For ω > ω_p the hot plasma mode has similar characteristics.

Effect of pressure on mechanical behavior of filled elastomers

The Young's modulus, stress-strain curves, and failure properties of glass bead-filled EPDM vulcanizates were studied under superposed hydrostatic pressure. The glass bead-filled EPDM was employed as a representation of composite systems, and the hydrostatic pressure controls the filler-elastomer separation under deformation. This separation shows up as a volume change of the system, and its infuence is reflected in the mechanical behavior as a reinforcing effect of variable degree.

The strain energy stored in the composite system in simple tension was calculated by introducing a model which is described as a cylindrical block of elastomer with two half spheres of filler on each end with their centers on the axis of the cylinder. In the derivation of the strain energy, assumptions were made to obtain the strain distribution in the model, and strain energy-strain relation for the elastomer was also assumed. The derivation was carried out for the case of no filler-elastomer separation and was modified to include the case of filler-elastomer separation.

The resulting strain energy, as a function of stretch ratio and volume of the system, was used to obtain stress-strain curves and volume change-strain curves of composite systems under superposed hydrostatic pressure.

Changes in the force and the lateral dimension of a ring specimen were measured as it was stretched axially under a superposed hydrostatic pressure in order to calculate the mechanical properties mentioned above. A tensile tester was used which is capable of sealing the whole system to carry out a measurement under pressure. A thickness measuring device, based on the Hall effect, was built for the measurement of changes in the lateral dimension of a specimen.

The theoretical and experimental results of Young's modulus and stress-strain curves were compared and showed fairly good agreement.

The failure data were discussed in terms of failure surfaces, and it was concluded that a failure surface of the glass-bead-filled EPDM consists of two cones.

Imperfection Insensitive Thin Shells

The buckling of axially compressed cylindrical shells and externally pressurized spherical shells is extremely sensitive to even very small geometric imperfections. In practice this issue is addressed by either using overly conservative knockdown factors, while keeping perfect axial or spherical symmetry, or adding closely and equally spaced stiffeners on shell surface. The influence of imperfection-sensitivity is mitigated, but the shells designed from these approaches are either too heavy or very expensive and are still sensitive to imperfections. Despite their drawbacks, these approaches have been used for more than half a century.

This thesis proposes a novel method to design imperfection-insensitive cylindrical shells subject to axial compression. Instead of following the classical paths, focused on axially symmetric or high-order rotationally symmetric cross-sections, the method in this thesis adopts optimal symmetry-breaking wavy cross-sections (wavy shells). The avoidance of imperfection sensitivity is achieved by searching with an evolutionary algorithm for smooth cross-sectional shapes that maximize the minimum among the buckling loads of geometrically perfect and imperfect wavy shells. It is found that the shells designed through this approach can achieve higher critical stresses and knockdown factors than any previously known monocoque cylindrical shells. It is also found that these shells have superior mass efficiency to almost all previously reported stiffened shells.

Experimental studies on a design of composite wavy shell obtained through the proposed method are presented in this thesis. A method of making composite wavy shells and a photogrametry technique of measuring full-field geometric imperfections have been developed. Numerical predictions based on the measured geometric imperfections match remarkably well with the experiments. Experimental results confirm that the wavy shells are not sensitive to imperfections and can carry axial compression with superior mass efficiency.

An efficient computational method for the buckling analysis of corrugated and stiffened cylindrical shells subject to axial compression has been developed in this thesis. This method modifies the traditional Bloch wave method based on the stiffness matrix method of rotationally periodic structures. A highly efficient algorithm has been developed to implement the modified Bloch wave method. This method is applied in buckling analyses of a series of corrugated composite cylindrical shells and a large-scale orthogonally stiffened aluminum cylindrical shell. Numerical examples show that the modified Bloch wave method can achieve very high accuracy and require much less computational time than linear and nonlinear analyses of detailed full finite element models.

This thesis presents parametric studies on a series of externally pressurized pseudo-spherical shells, i.e., polyhedral shells, including icosahedron, geodesic shells, and triambic icosahedra. Several optimization methods have been developed to further improve the performance of pseudo-spherical shells under external pressure. It has been shown that the buckling pressures of the shell designs obtained from the optimizations are much higher than the spherical shells and not sensitive to imperfections.

Experimental study on inertial effects in liquid-solid flows

This thesis presents experimental measurements of the rheological behavior of liquid-solid mixtures at moderate Reynolds (defined by the shear rate and particle diameter) and Stokes numbers, ranging from 3 ≤ Re ≤ 1.6 × 10³ and 0.4 ≤ St ≤ 195. The experiments use a specifically designed Couette cylindrical rheometer that allows for probing the transition from transporting a pure liquid to transporting a dense suspension of particles. Measurements of the shear stress are presented for a wide range of particle concentration (10 to 60% in volume) and for particle to fluid density ratio between 1 and 1.05. The effective relative viscosity exhibits a strong dependence on the solid fraction for all density ratios tested. For density ratio of 1 the effective viscosity increases with Stokes number (St) for volume fractions (φ) lower than 40% and becomes constant for higher φ. When the particles are denser than the liquid, the effective viscosity shows a stronger dependance on St. An analysis of the particle resuspension for the case with a density ratio of 1.05 is presented and used to predict the local volume fraction where the shear stress measurements take place. When the local volume fraction is considered, the effective viscosity for settling and no settling particles is consistent, indicating that the effective viscosity is independent of differences in density between the solid and liquid phase. Shear stress measurements of pure fluids (no particles) were performed using the same rheometer, and a deviation from laminar behavior is observed for gap Reynolds numbers above 4× 10³, indicating the presence of hydrodynamic instabilities associated with the rotation of the outer cylinder. The increase on the effective viscosity with Stokes numbers observed for mixtures with φ ≤ 30% appears to be affected by such hydrodynamic instabilities. The effective viscosity for the current experiments is considerably higher than the one reported in non-inertial suspensions.

Experimental, Numerical and Analytical Studies of the MHD-driven plasma jet, instabilities and waves

This thesis describes a series of experimental, numerical, and analytical studies involving the Caltech magnetohydrodynamically (MHD)-driven plasma jet experiment. The plasma jet is created via a capacitor discharge that powers a magnetized coaxial planar electrodes system. The jet is collimated and accelerated by the MHD forces.

We present three-dimensional ideal MHD finite-volume simulations of the plasma jet experiment using an astrophysical magnetic tower as the baseline model. A compact magnetic energy/helicity injection is exploited in the simulation analogous to both the experiment and to astrophysical situations. Detailed analysis provides a comprehensive description of the interplay of magnetic force, pressure, and flow effects. We delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms.

When the experimental jet is sufficiently long, it undergoes a global kink instability and then a secondary local Rayleigh-Taylor instability caused by lateral acceleration of the kink instability. We present an MHD theory of the Rayleigh-Taylor instability on the cylindrical surface of a plasma flux rope in the presence of a lateral external gravity. The Rayleigh-Taylor instability is found to couple to the classic current-driven instability, resulting in a new type of hybrid instability. The coupled instability, produced by combination of helical magnetic field, curvature of the cylindrical geometry, and lateral gravity, is fundamentally different from the classic magnetic Rayleigh-Taylor instability occurring at a two-dimensional planar interface.

In the experiment, this instability cascade from macro-scale to micro-scale eventually leads to the failure of MHD. When the Rayleigh-Taylor instability becomes nonlinear, it compresses and pinches the plasma jet to a scale smaller than the ion skin depth and triggers a fast magnetic reconnection. We built a specially designed high-speed 3D magnetic probe and successfully detected the high frequency magnetic fluctuations of broadband whistler waves associated with the fast reconnection. The magnetic fluctuations exhibit power-law spectra. The magnetic components of single-frequency whistler waves are found to be circularly polarized regardless of the angle between the wave propagation direction and the background magnetic field.