Theory of viscous and thermal attenuation of sound by small spheres

The problem is to calculate the attenuation of plane sound waves passing through a viscous, heat-conducting fluid containing small spherical inhomogeneities. The attenuation is calculated by evaluating the rate of increase of entropy caused by two irreversible processes: (1) the mechanical work done by the viscous stresses in the presence of velocity gradients, and (2) the flow of heat down the thermal gradients. The method is first applied to a homogeneous fluid with no spheres and shown to give the classical Stokes-Kirchhoff expressions. The method is then used to calculate the additional viscous and thermal attenuation when small spheres are present. The viscous attenuation agrees with Epstein's result obtained in 1941 for a non-heat-conducting fluid. The thermal attenuation is found to be similar in form to the viscous attenuation and, for gases, of comparable magnitude. The general results are applied to the case of water drops in air and air bubbles in water.

For water drops in air the viscous and thermal attenuations are camparable; the thermal losses occur almost entirely in the air, the thermal dissipation in the water being negligible. The theoretical values are compared with Knudsen's experimental data for fogs and found to agree in order of magnitude and dependence on frequency. For air bubbles in water the viscous losses are negligible and the calculated attenuation is almost completely due to thermal losses occurring in the air inside the bubbles, the thermal dissipation in the water being relatively small. (These results apply only to non-resonant bubbles whose radius changes but slightly during the acoustic cycle.)

On the thermal budget of Pahoehoe lava flows

In this thesis I investigate some aspects of the thermal budget of pahoehoe lava flows. This is done with a combination of general field observations, quantitative modeling, and specific field experiments. The results of this work apply to pahoehoe flows in general, even though the vast bulk of the work has been conducted on the lavas formed by the Pu'u 'O'o - Kupaianaha eruption of Kilauea Volcano on Hawai'i. The field observations rely heavily on discussions with the staff of the United States Geological Survey's Hawaiian Volcano Observatory (HVO), under whom I labored repeatedly in 1991-1993 for a period totaling about 10 months.

The quantitative models I have constructed are based on the physical processes observed by others and myself to be active on pahoehoe lava flows. By building up these models from the basic physical principles involved, this work avoids many of the pitfalls of earlier attempts to fit field observations with "intuitively appropriate" mathematical expressions. Unlike many earlier works, my model results can be analyzed in terms of the interactions between the different physical processes. I constructed models to: (1) describe the initial cooling of small pahoehoe flow lobes and (2) understand the thermal budget of lava tubes.

The field experiments were designed either to validate model results or to constrain key input parameters. In support of the cooling model for pahoehoe flow lobes, attempts were made to measure: (1) the cooling within the flow lobes, (2) the amount of heat transported away from the lava by wind, and (3) the growth of the crust on the lobes. Field data collected by Jones [1992], Hon et al. [1994b], and Denlinger [Keszthelyi and Denlinger, in prep.] were also particularly useful in constraining my cooling model for flow lobes. Most of the field observations I have used to constrain the thermal budget of lava tubes were collected by HVO (geological and geophysical monitoring) and the Jet Propulsion Laboratory (airborne infrared imagery [Realmuto et al., 1992]). I was able to assist HVO for part of their lava tube monitoring program and also to collect helicopterborne and ground-based IR video in collaboration with JPL [Keszthelyi et al., 1993].

The most significant results of this work are (1) the quantitative demonstration that the emplacement of pahoehoe and 'a'a flows are the fundamentally different, (2) confirmation that even the longest lava flows observed in our Solar System could have formed as low effusion rate, tube-fed pahoehoe flows, and (3) the recognition that the atmosphere plays a very important role throughout the cooling of history of pahoehoe lava flows. In addition to answering specific questions about the thermal budget of tube-fed pahoehoe lava flows, this thesis has led to some additional, more general, insights into the emplacement of these lava flows. This general understanding of the tube-fed pahoehoe lava flow as a system has suggested foci for future research in this part of physical volcanology.

Non-linear dispersive waves with a small dissipation

The general theory of Whitham for slowly-varying non-linear wavetrains is extended to the case where some of the defining partial differential equations cannot be put into conservation form. Typical examples are considered in plasma dynamics and water waves in which the lack of a conservation form is due to dissipation; an additional non-conservative element, the presence of an external force, is treated for the plasma dynamics example. Certain numerical solutions of the water waves problem (the Korteweg-de Vries equation with dissipation) are considered and compared with perturbation expansions about the linearized solution; it is found that the first correction term in the perturbation expansion is an excellent qualitative indicator of the deviation of the dissipative decay rate from linearity.

A method for deriving necessary and sufficient conditions for the existence of a general uniform wavetrain solution is presented and illustrated in the plasma dynamics problem. Peaking of the plasma wave is demonstrated, and it is shown that the necessary and sufficient existence conditions are essentially equivalent to the statement that no wave may have an amplitude larger than the peaked wave.

A new type of fully non-linear stability criterion is developed for the plasma uniform wavetrain. It is shown explicitly that this wavetrain is stable in the near-linear limit. The nature of this new type of stability is discussed.

Steady shock solutions are also considered. By a quite general method, it is demonstrated that the plasma equations studied here have no steady shock solutions whatsoever. A special type of steady shock is proposed, in which a uniform wavetrain joins across a jump discontinuity to a constant state. Such shocks may indeed exist for the Korteweg-de Vries equation, but are barred from the plasma problem because entropy would decrease across the shock front.

Finally, a way of including the Landau damping mechanism in the plasma equations is given. It involves putting in a dissipation term of convolution integral form, and parallels a similar approach of Whitham in water wave theory. An important application of this would be towards resolving long-standing difficulties about the "collisionless" shock.

Thin metastructures with engineered thermal expansion

The geometry and constituent materials of metastructures can be used to engineer the thermal expansion coefficient. In this thesis, we design, fabricate, and test thin thermally stable metastructures consisting of bi-metallic unit cells and show how the coefficient of thermal expansion (CTE) of these metastructures can be finely and coarsely tuned by varying the CTE of the constituent materials and the unit cell geometry. Planar and three-dimensional finite element method modeling is used to drive the design and inform experiments, and predict the response of these metastructures. We demonstrate computationally the significance of out-of-plane effects in the metastructure response. We develop an experimental setup using digital image correlation and an infrared camera to experimentally measure full displacement and temperature fields during testing and accurately measure the metastructures’ CTE. We experimentally demonstrate high aspect ratio metastructures of Ti/Al and Kovar/Al which exhibit near-zero and negative CTE, respectively. We demonstrate robust fabrication procedures for thermally stable samples with high aspect ratios in thin foil and thin film scales. We investigate the lattice structure and mechanical properties of thin films comprising a near-zero CTE metastructure. The mechanics developed in this work can be used to engineer metastructures of arbitrary CTE and can be extended to three dimensions.

Fracture of materials undergoing solid-solid phase transformation

A large number of technologically important materials undergo solid-solid phase transformations. Examples range from ferroelectrics (transducers and memory devices), zirconia (Thermal Barrier Coatings) to nickel superalloys and (lithium) iron phosphate (Li-ion batteries). These transformations involve a change in the crystal structure either through diffusion of species or local rearrangement of atoms. This change of crystal structure leads to a macroscopic change of shape or volume or both and results in internal stresses during the transformation. In certain situations this stress field gives rise to cracks (tin, iron phosphate etc.) which continue to propagate as the transformation front traverses the material. In other materials the transformation modifies the stress field around cracks and effects crack growth behavior (zirconia, ferroelectrics). These observations serve as our motivation to study cracks in solids undergoing phase transformations. Understanding these effects will help in improving the mechanical reliability of the devices employing these materials.

In this thesis we present work on two problems concerning the interplay between cracks and phase transformations. First, we consider the directional growth of a set of parallel edge cracks due to a solid-solid transformation. We conclude from our analysis that phase transformations can lead to formation of parallel edge cracks when the transformation strain satisfies certain conditions and the resulting cracks grow all the way till their tips cross over the phase boundary. Moreover the cracks continue to grow as the phase boundary traverses into the interior of the body at a uniform spacing without any instabilities. There exists an optimal value for the spacing between the cracks. We ascertain these conclusion by performing numerical simulations using finite elements.

Second, we model the effect of the semiconducting nature and dopants on cracks in ferroelectric perovskite materials, particularly barium titanate. Traditional approaches to model fracture in these materials have treated them as insulators. In reality, they are wide bandgap semiconductors with oxygen vacancies and trace impurities acting as dopants. We incorporate the space charge arising due the semiconducting effect and dopant ionization in a phase field model for the ferroelectric. We derive the governing equations by invoking the dissipation inequality over a ferroelectric domain containing a crack. This approach also yields the driving force acting on the crack. Our phase field simulations of polarization domain evolution around a crack show the accumulation of electronic charge on the crack surface making it more permeable than was previously believed so, as seen in recent experiments. We also discuss the effect the space charge has on domain formation and the crack driving force.

High fidelity probe and mitigation of mirror thermal fluctuations

Thermal noise arising from mechanical loss in high reflective dielectric coatings is a significant source of noise in precision optical measurements. In particular, Advanced LIGO, a large scale interferometer aiming to observed gravitational wave, is expected to be limited by coating thermal noise in the most sensitive region around 30–300 Hz. Various theoretical calculations for predicting coating Brownian noise have been proposed. However, due to the relatively limited knowledge of the coating material properties, an accurate approximation of the noise cannot be achieved. A testbed that can directly observed coating thermal noise close to Advanced LIGO band will serve as an indispensable tool to verify the calculations, study material properties of the coating, and estimate the detector’s performance.

This dissertation reports a setup that has sensitivity to observe wide band (10Hz to 1kHz) thermal noise from fused silica/tantala coating at room temperature from fixed-spacer Fabry–Perot cavities. Important fundamental noises and technical noises associated with the setup are discussed. The coating loss obtained from the measurement agrees with results reported in the literature. The setup serves as a testbed to study thermal noise in high reflective mirrors from different materials. One example is a heterostructure of Al_xGa_1−xAs (AlGaAs). An optimized design to minimize thermo–optic noise in the coating is proposed and discussed in this work.

Investigations of noise and of quantum interference in proximity effect bridges

This work reports investigations upon weakly superconducting proximity effect bridges. These bridges, which exhibit the Josephson effects, are produced by bisecting a superconductor with a short (<1µ) region of material whose superconducting transition temperature is below that of the adjacent superconductors. These bridges are fabricated from layered refractory metal thin films whose transition temperature will depend upon the thickness ratio of the materials involved. The thickness ratio is changed in the area of the bridge to lower its transition temperature. This is done through novel photolithographic techniques described in the text, Chapter 2.

If two such proximity effect bridges are connected in parallel, they form a quantum interferometer. The maximum zero voltage current through this circuit is periodically modulated by the magnetic flux through the circuit. At a constant bias current, the modulation of the critical current produces a modulation in the dc voltage across the bridge. This change in dc voltage has been found to be the result of a change in the internal dissipation in the device. A simple model using lumped circuit theory and treating the bridges as quantum oscillators of frequency ω = 2eV/h, where V is the time average voltage across the device, has been found to adequately describe the observed voltage modulation.

The quantum interferometers have been converted to a galvanometer through the inclusion of an integral thin film current path which couples magnetic flux through the interferometer. Thus a change in signal current produces a change in the voltage across the interferometer at a constant bias current. This work is described in Chapter 3 of the text.

The sensitivity of any device incorporating proximity effect bridges will ultimately be determined by the fluctuations in their electrical parameters. He have measured the spectral power density of the voltage fluctuations in proximity effect bridges using a room temperature electronics and a liquid helium temperature transformer to match the very low (~ 0.1 Ω) impedances characteristic of these devices.

We find the voltage noise to agree quite well with that predicted by phonon noise in the normal conduction through the bridge plus a contribution from the superconducting pair current through the bridge which is proportional to the ratios of this current to the time average voltage across the bridge. The total voltage fluctuations are given by <V^2(f ) > = 4kTR^2_d I/V where R_d is the dynamic resistance, I the total current, and V the voltage across the bridge . An additional noise source appears with a strong 1/f^(n) dependence , 1.5 < n < 2, if the bridges are fabricated upon a glass substrate. This excess noise, attributed to thermodynamic temperature fluctuations in the volume of the bridge, increases dramatically on a glass substrate due to the greatly diminished thermal diffusivity of the glass as compared to sapphire.

Hydrothermal experiments on the thermal stability of amino substances in sediments

Aspartic acid, threonine, serine and other thermally unstable amino acids have been found in fine-grained elastic sediments of advanced geologic age. The presence of these compounds in ancient sediments conflicts with experimental data determined for their simple thermal decomposition.

Recent and Late Miocene sediments and their humic acid extracts, known to contain essentially complete suites of amino acids, were heated with H₂O in a bomb at temperatures up to 500°C in order to compare the thermal decomposition characteristics of the sedimentary amino compounds.

Most of the amino acids found in protein hydrolyzates are obtained from the Miocene rock in amounts 10 to 100 times less than from the Recent sediment. The two unheated humic acids are rather similar despite their great age difference. The Miocene rock appears uncontaminated by Recent carbon.

Yields of amino acids generally decline in the heated Recent sediment. Some amino compounds apparently increase with heating time in the Miocene rock.

Relative thermal stabilities of the amino acids in sediments are generally similar to those determined using pure aqueous solutions. The relative thermal stabilities of glutamic acid, glycine, and phenylalanine vary in the Recent sediment but are uniform in the Miocene rock.

Amino acids may occur in both proteins and humic complexes in the Recent sediment, while they are probably only present in stabilized organic substances in the Miocene rock. Thermal decomposition of protein amino acids may be affected by surface catalysis in the Recent sediment. The apparent activation energy for the decomposition of alanine in this sediment is 8400 calories per mole. Yields of amino compounds from the heated sediments are not affected by thermal decomposition only.

Amino acids in sediments may only be useful for geothermometry in a very general way.

A better picture of the amino acid content of older sedimentary rocks may be obtained if these sediments are heated in a bomb with H₂O at temperatures around 150°C prior to HCl hydrolysis.

Leucine-isoleucine ratios may prove to be useful as indicators of amino acid sources or for evaluating the fractionation of these substances during diagenesis. Leucine-isoleucine ratios of the Recent and Miocene sediments and humic acids are identical. The humic acids may have a continental source.

The carbon-nitrogen and carbon-hydrogen ratios of sediments and humic acids increase with heating time and temperature. Ratios comparable to those in some kerogens are found in the severely heated Miocene sediment and humic acid.

Bicyclo[3.2.0]hepta-1,4,6-triene: synthesis, thermal rearrangement, and anion formation

In order to determine the properties of the bicycloheptatrienyl anion (Ia) (predicted to be conjugatively stabilized by Hückel Molecular Orbital Theory) the neutral precursor, bicyclo[3. 2. 0] hepta-1, 4, 6-triene (I) was prepared by the following route.

Reaction of I with potassium-t-butoxide, potassium, or lithium dicyclohexylamide gave anion Ia in very low yield. Reprotonation of I was found to occur solely at the 1 or 5 position to give triene II, isolated as to its dimers.

A study of the acidity of I and of other conjugated hydrocarbons by means of ion cyclotron resonance spectroscopy resulted in determination of the following order of relative acidities:

H₂S ˃ C₅H₆ ˃ CH₃NO₂ ˃ 1, 4- C₅H₈ ˃ I ˃ C₂H₅OH ˃ H₂O; cyclo-C₇H₈ ˃ C₂ H₅OH; фCH₃ ˃ CH₃OH

In addition, limits for the proton affinities of the conjugate bases were determined:

350 kcal/mole ˂ PA(C₅ H₅^-) ˂ 360 kcal/mole

362 kcal/mole ˂ PA(C₅H₇^-, Ia, cyclo-C₇H₇^-) ˂ 377 kcal/mole PA(фCH₂^-) ˂ 385 kcal/mole

Gas phase kinetics of the trans-XVIII to I transformation gave the following activation parameters: E_a = 43.0 kcal/mole, log A = 15.53 and ∆S^ǂ (220°) = 9.6 cu. The results were interpreted as indicating initial 1,2 bond cleavage to give the 1,3-diradical which closed to I. Similar studies on cis-XVIII gave results consistent with a surface component to the reaction (E_a = 22.7 kcal/mole; log A = 9.23, ∆S^ǂ (119°) = -18.9 eu).

The low pressure (0.01 to 1 torr) pyrolysis of trans-XVIII gave in addition to I, fulvenallene (LV), ethynylcyclopentadiene (LVI) and heptafulvalene (LVII). The relative ratios of the C₇H₆ isomers were found to be dependent upon temperature and pressure, higher relative pressure and lower temperatures favoring formation of I. The results were found to be consistent with the intermediacy of vibrationally excited I and subsequent reaction to give LV and LVI.