General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

I. Configurational stability and redistribution equilibria in organomagnesium compounds. II. Redistribution equilibria in organocadmium compounds

I. CONFIGURATIONAL STABILITY AND REDISTRIBUTION EQUILIBRIA IN ORGANOMAGNESIUM COMPOUNDS

The dependence of the rate of inversion of a dialkylmagnesium compound on the solvent has been studied.

Examination of the temperature dependence of the nuclear magnetic resonance spectrum of 1-phenyl-2-propylmagnesium bromide in diethyl ether solution indicates that inversion of configuration at the methylene group of this Grignard reagent occurs with an approximate rate of 2 sec^-1 at room temperature. This is the first example of a rapid inversion rate in a secondary Grignard reagent.

The rates of exchange of alkyl groups between dineopentylmagnesium and di-s-butylmagnesium, bis-(2-methylbutyl)-magnesium and bis-(4, 4-dimethyl-2-pentyl)-magnesium respectively in diethyl ether solution were found to be fast on the nmr time scale. However, the alkyl group exchange rate was found to be slow in a diethyl ether solution of dineopentylmagnesium and bis-(2-methylbutyl)-magnesium containing N, N, N', N'-tetramethylethylenediamine. The unsymmetrical species neopentyl-2-methylbutyl-magnesium was observed at room temperature in the nmr spectrum of the solution containing the diamine.

II. REDISTRIBUTION EQUILIBRIA IN ORGANOCADMIUM COMPOUNDS

The exchange of methyl groups in dimethylcadmium has been studied by nuclear magnetic resonance spectroscopy. Activation parameters for the methyl group exchange have been measured for a neat sample and for a solution in tetrahydrofuran. The exchange is faster in the basic solvent tetrahydrofuran relative to the neat sample and in tetrahydrofuran solution is retarded by the solvating agent N, N, N’, N’-tetramethylethylenediamine and greatly increased by cadmium bromide. The addition of methanol to a solution of dimethylcadmium in tetrahydrofuran appears to have very little effect on the rate of exchange. The exchange was found to proceed with retention of configuration. The rate-limiting step for the exchange of methyl groups in a basic solvent appears to be the dissociation of coordinating solvent from dimethylcadmium.

The equilibrium between methylcadmium bromide, dimethylcadmium and cadmium bromide in tetrahydrofuran solution has also been studied. At room temperature the interconversion of the species is very fast on the nmr time scale but at -100° distinct absorptions for methylcadmium bromide and imethylcadmium are observed.

The species ethylmethylcadmium has been observed in the nmr spectrum.

The rate of exchange of vinyl groups in a solution of divinylcadmium in tetrahydrofuran has been found to be fast on the nmr time scale.

The stability of traveling wave solutions of parabolic equations

A method for determining by inspection the stability or instability of any solution u(t,x) = ɸ(x-ct) of any smooth equation of the form u_t = f(u_(xx),u_x,u where ∂/∂a f(a,b,c) > 0 for all arguments a,b,c, is developed. The connection between the mean wavespeed of solutions u(t,x) and their initial conditions u(0,x) is also explored. The mean wavespeed results and some of the stability results are then extended to include equations which contain integrals and also to include some special systems of equations. The results are applied to several physical examples.

Slender-body hypervelocity boundary-layer instability

With novel application of optical techniques, the slender-body hypervelocity boundary-layer instability is characterized in the previously unexplored regime where thermo-chemical effects are important. Narrowband disturbances (500-3000~kHz) are measured in boundary layers with edge velocities of up to 5~km/s at two points along the generator of a 5 degree half angle cone. Experimental amplification factor spectra are presented. Linear stability and PSE analysis is performed, with fair prediction of the frequency content of the disturbances; however, the analysis over-predicts the amplification of disturbances. The results of this work have two key implications: 1) the acoustic instability is present and may be studied in a large-scale hypervelocity reflected-shock tunnel, and 2) the new data set provides a new basis on which the instability can be studied.

Periodic solutions of integro-differential equations which arise in population dynamics

The problem of the existence and stability of periodic solutions of infinite-lag integra-differential equations is considered. Specifically, the integrals involved are of the convolution type with the dependent variable being integrated over the range (- ∞,t), as occur in models of population growth. It is shown that Hopf bifurcation of periodic solutions from a steady state can occur, when a pair of eigenvalues crosses the imaginary axis. Also considered is the existence of traveling wave solutions of a model population equation allowing spatial diffusion in addition to the usual temporal variation. Lastly, the stability of the periodic solutions resulting from Hopf bifurcation is determined with aid of a Floquet theory.

The first chapter is devoted to linear integro-differential equations with constant coefficients utilizing the method of semi-groups of operators. The second chapter analyzes the Hopf bifurcation providing an existence theorem. Also, the two-timing perturbation procedure is applied to construct the periodic solutions. The third chapter uses two-timing to obtain traveling wave solutions of the diffusive model, as well as providing an existence theorem. The fourth chapter develops a Floquet theory for linear integro-differential equations with periodic coefficients again using the semi-group approach. The fifth chapter gives sufficient conditions for the stability or instability of a periodic solution in terms of the linearization of the equations. These results are then applied to the Hopf bifurcation problem and to a certain population equation modeling periodically fluctuating environments to deduce the stability of the corresponding periodic solutions.

Limits and tradeoffs in the control of autocatalytic systems

Despite the complexity of biological networks, we find that certain common architectures govern network structures. These architectures impose fundamental constraints on system performance and create tradeoffs that the system must balance in the face of uncertainty in the environment. This means that while a system may be optimized for a specific function through evolution, the optimal achievable state must follow these constraints. One such constraining architecture is autocatalysis, as seen in many biological networks including glycolysis and ribosomal protein synthesis. Using a minimal model, we show that ATP autocatalysis in glycolysis imposes stability and performance constraints and that the experimentally well-studied glycolytic oscillations are in fact a consequence of a tradeoff between error minimization and stability. We also show that additional complexity in the network results in increased robustness. Ribosome synthesis is also autocatalytic where ribosomes must be used to make more ribosomal proteins. When ribosomes have higher protein content, the autocatalysis is increased. We show that this autocatalysis destabilizes the system, slows down response, and also constrains the system’s performance. On a larger scale, transcriptional regulation of whole organisms also follows architectural constraints and this can be seen in the differences between bacterial and yeast transcription networks. We show that the degree distributions of bacterial transcription network follow a power law distribution while the yeast network follows an exponential distribution. We then explored the evolutionary models that have previously been proposed and show that neither the preferential linking model nor the duplication-divergence model of network evolution generates the power-law, hierarchical structure found in bacteria. However, in real biological systems, the generation of new nodes occurs through both duplication and horizontal gene transfers, and we show that a biologically reasonable combination of the two mechanisms generates the desired network.

Mechanistic and reactivity studies of bent metallocene complexes of niobium and tantalum

A variety of olefin hydride complexes of niobium and tantalum has been prepared in order to study their reactivity and to gain insight into organometallic reaction mechanisms. Examination of a series of ethylene and propylene complexes of niobocene (CP_2Nb; Cp = η^5-C_5H_5), permethylniobocene (Cp*_2Nb; Cp* = η^5-C_5(CH_3)_5), tantalocene, and permethyltantalocene has indicated that there are both large electronic and steric effects deriving from the metal (and its ancillary ligands) in the olefin insertion (β-migratory insertion) process. Furthermore, a thermodynamic and kinetic analysis has been completed for a series of substituted styrene complexes of niobocene in order to better understand the important electronic properties of the olefin. The results are in accord with a concerted four-center process with only moderate charge development.

The special case of β-migratory insertion of a hydride ligand into coordinated benzyne has also been studied for the permethyltantalocene system. The coordinatively unsaturated (sixteen electron) phenyl tautomer, which is made accessible by the facile benzyne hydride insertion reaction, readily reacts with a variety of ligands, L, to afford Cp*_2 Ta(C_6H_5)L complexes (L = CO, O_2, NC≡R, :CH_2, H_2, etc.). This family of compounds exhibits interesting reactivity (a-migratory insertion, O_2 activation, and reductive elimination) which is discussed in some detail.

Finally a series of paramagnetic seventeen electron Cp*_2 TaX_2 (X = halide, alkyl, hydride) complexes, and the corresponding cationic and anionic species, have been prepared and studied. The odd electron neutral complexes exhibit surprising thermal stability and undergo very little reactivity. While the chemistry of the anionic compounds is almost completely dominated by their potent reducing power, that of the cations is quite diverse and amenable for study. Therefore the syntheses and reactivity (1 ,2-eliminations, ligand insertions, and deprotonation reactions) of these coordinatively unsaturated sixteen electron species are presented.

Synthesis and characterization of high-spin organic materials: prototypes for the polaronic ferromagnet

The design, synthesis and magnetic characterization of thiophene-based models for the polaronic ferromagnet are described. Synthetic strategies employing Wittig and Suzuki coupling were employed to produce polymers with extended π-systems. Oxidative doping using AsF_5 or I_2 produces radical cations (polarons) that are stable at room temperature. Magnetic characterization of the doped polymers, using SQUID-based magnetometry, indicates that in several instances ferromagnetic coupling of polarons occurs along the polymer chain. An investigation of the influence of polaron stability and delocalization on the magnitude of ferromagnetic coupling is pursued. A lower limit for mild, solution phase I_2 doping is established. A comparison of the variable temperature data of various polymers reveals that deleterious antiferromagnetic interactions are relatively insensitive to spin concentration, doping protocols or spin state. Comparison of the various polymers reveals useful design principles and suggests new directions for the development of magnetic organic materials. Novel strategies for solubilizing neutral polymeric materials in polar solvents are investigated.

The incorporation of stable bipyridinium spin-containing units into a polymeric high-spin array is explored. Preliminary results suggest that substituted diquat derivatives may serve as stable spin-containing units for the polaronic ferromagnet and are amenable to electrochemical doping. Synthetic efforts to prepare high-spin polymeric materials using viologens as a spin source have been unsuccessful.

A model for energy and morphology of crystalline grain boundaries with arbitrary geometric character

It has been well-established that interfaces in crystalline materials are key players in the mechanics of a variety of mesoscopic processes such as solidification, recrystallization, grain boundary migration, and severe plastic deformation. In particular, interfaces with complex morphologies have been observed to play a crucial role in many micromechanical phenomena such as grain boundary migration, stability, and twinning. Interfaces are a unique type of material defect in that they demonstrate a breadth of behavior and characteristics eluding simplified descriptions. Indeed, modeling the complex and diverse behavior of interfaces is still an active area of research, and to the author's knowledge there are as yet no predictive models for the energy and morphology of interfaces with arbitrary character. The aim of this thesis is to develop a novel model for interface energy and morphology that i) provides accurate results (especially regarding "energy cusp" locations) for interfaces with arbitrary character, ii) depends on a small set of material parameters, and iii) is fast enough to incorporate into large scale simulations.

In the first half of the work, a model for planar, immiscible grain boundary is formulated. By building on the assumption that anisotropic grain boundary energetics are dominated by geometry and crystallography, a construction on lattice density functions (referred to as "covariance") is introduced that provides a geometric measure of the order of an interface. Covariance forms the basis for a fully general model of the energy of a planar interface, and it is demonstrated by comparison with a wide selection of molecular dynamics energy data for FCC and BCC tilt and twist boundaries that the model accurately reproduces the energy landscape using only three material parameters. It is observed that the planar constraint on the model is, in some cases, over-restrictive; this motivates an extension of the model.

In the second half of the work, the theory of faceting in interfaces is developed and applied to the planar interface model for grain boundaries. Building on previous work in mathematics and materials science, an algorithm is formulated that returns the minimal possible energy attainable by relaxation and the corresponding relaxed morphology for a given planar energy model. It is shown that the relaxation significantly improves the energy results of the planar covariance model for FCC and BCC tilt and twist boundaries. The ability of the model to accurately predict faceting patterns is demonstrated by comparison to molecular dynamics energy data and experimental morphological observation for asymmetric tilt grain boundaries. It is also demonstrated that by varying the temperature in the planar covariance model, it is possible to reproduce a priori the experimentally observed effects of temperature on facet formation.

Finally, the range and scope of the covariance and relaxation models, having been demonstrated by means of extensive MD and experimental comparison, future applications and implementations of the model are explored.

Stability of parametrically excited differential equations

Sufficient stability criteria for classes of parametrically excited differential equations are developed and applied to example problems of a dynamical nature.

Stability requirements are presented in terms of 1) the modulus of the amplitude of the parametric terms, 2) the modulus of the integral of the parametric terms and 3) the modulus of the derivative of the parametric terms.

The methods employed to show stability are Liapunov’s Direct Method and the Gronwall Lemma. The type of stability is generally referred to as asymptotic stability in the sense of Liapunov.

The results indicate that if the equation of the system with the parametric terms set equal to zero exhibits stability and possesses bounded operators, then the system will be stable under sufficiently small modulus of the parametric terms or sufficiently small modulus of the integral of the parametric terms (high frequency). On the other hand, if the equation of the system exhibits individual stability for all values that the parameter assumes in the time interval, then the actual system will be stable under sufficiently small modulus of the derivative of the parametric terms (slowly varying).