Whirl and Stall Flutter Simulation Using CFD

This paper presents recent research on numerical methods for whirl and stall flutter using computational fluid dynamics. The method involves coupling of the HMB3 CFD solver of the University of Glasgow and a NASTRAN derived structural model. Based upon a literature survey, a significant amount of research has been conducted on the numerical investigation of tiltrotors, with a focus on the XV-15 and V-22 aircraft. Within this paper, the coupling procedure is presented along with a steady CFD computation to highlight the accuracy of the high-fidelity method. In addition to this, a simple method is used to investigate the whirl flutter boundary of a standard propeller and the XV-15 blade

Inference of epidemiological parameters from household stratified data

We consider a continuous-time Markov chain model of SIR disease dynamics with two levels of mixing. For this so-called stochastic households model, we provide two methods for inferring the model parameters---governing within-household transmission, recovery, and between-household transmission---from data of the day upon which each individual became infectious and the household in which each infection occurred, as would be available from first few hundred studies. Each method is a form of Bayesian Markov Chain Monte Carlo that allows us to calculate a joint posterior distribution for all parameters and hence the household reproduction number and the early growth rate of the epidemic. The first method performs exact Bayesian inference using a standard data-augmentation approach; the second performs approximate Bayesian inference based on a likelihood approximation derived from branching processes. These methods are compared for computational efficiency and posteriors from each are compared. The branching process is shown to be an excellent approximation and remains computationally efficient as the amount of data is increased

Report of the sensor cooler technology panel

The Sensor Cooler Technology Panel identified three major areas in which technology development must be supported in order to meet the system performance requirements for the Astrotech 21 mission set science objectives. They are: long life vibration free refrigerators; mechanical refrigeration for 2 K to 5 K; and flight testing of emerging prototype refrigerators. A development strategy and schedule were recommended for each of the three areas

X-ray Reflection from Inhomogeneous Accretion Disks: II. Emission Line Variability and Implications for Reverberation Mapping

One of the principal scientific objectives of the upcoming Constellation-X mission is to attempt to map the inner regions of accretion disks around black holes in Seyfert galaxies by reverberation mapping of the Fe K fluorescence line. This area of the disk is likely radiation pressure dominated and subject to various dynamical instabilities. Here, we show that density inhomogeneities in the disk atmosphere resulting from the photon bubble instability (PBI) can cause rapid changes in the X-ray reflection features, even when the illuminating flux is constant. Using a simulation of the development of the PBI, we find that, for the disk parameters chosen, the Fe K and O VIII Ly\alpha lines vary on timescales as short as a few hundredths of an orbital time. In response to the changes in accretion disk structure, the Fe K equivalent width (EW) shows variations as large as ~100 eV. The magnitude and direction (positive or negative) of the changes depends on the ionization state of the atmosphere. The largest changes are found when the disk is moderately ionized. The O VIII EW varies by tens of eV, as well as exhibiting plenty of rapid, low-amplitude changes. This effect provides a natural explanation for some observed instances of short timescale Fe K variability which was uncorrelated with the continuum (e.g., Mrk 841). New predictions for Fe K reverberation mapping should be made which include the effects of this accretion disk driven line variability and a variable ionization state. Reflection spectra averaged over the evolution of the instability are well fit by constant density models in the 2-10 keV region.Comment: 20 pages, 3 figures. Accepted by Ap

Fast Mesh Refinement in Pseudospectral Optimal Control

Mesh refinement in pseudospectral (PS) optimal control is embarrassingly easy --- simply increase the order $N$ of the Lagrange interpolating polynomial and the mathematics of convergence automates the distribution of the grid points. Unfortunately, as $N$ increases, the condition number of the resulting linear algebra increases as $N^2$; hence, spectral efficiency and accuracy are lost in practice. In this paper, we advance Birkhoff interpolation concepts over an arbitrary grid to generate well-conditioned PS optimal control discretizations. We show that the condition number increases only as $\sqrt{N}$ in general, but is independent of $N$ for the special case of one of the boundary points being fixed. Hence, spectral accuracy and efficiency are maintained as $N$ increases. The effectiveness of the resulting fast mesh refinement strategy is demonstrated by using \underline{polynomials of over a thousandth order} to solve a low-thrust, long-duration orbit transfer problem.Comment: 27 pages, 12 figures, JGCD April 201

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

The quasiharmonic approximation (QHA), in its simplest form also called the statically constrained (SC) QHA, has been shown to be a straightforward method to compute thermoelastic properties of crystals. Recently we showed that for non-cubic solids SC-QHA calculations develop deviatoric thermal stresses at high temperatures. Relaxation of these stresses leads to a series of corrections to the free energy that may be taken to any desired order, up to self-consistency. Here we show how to correct the elastic constants obtained using the SC-QHA. We exemplify the procedure by correcting to first order the elastic constants of MgSiO$_3$-perovskite and MgSiO$_3$-post-perovskite, the major phases of the Earth's lower mantle. We show that this first order correction is quite satisfactory for obtaining the aggregated elastic averages of these minerals and their velocities in the lower mantle. This type of correction is also shown to be applicable to experimental measurements of elastic constants in situations where deviatoric stresses can develop, such as in diamond anvil cells.Comment: 4 figures, 1 table, submitted to Phys. Rev. B, July 200

Spin Fluctuations and the Pseudogap in Organic Superconductors

We show that there are strong similarities in the spin lattice relaxation of non-magnetic organic charge transfer salts, and that these similarities can be understood in terms of spin fluctuations. Further, we show that, in all of the kappa-phase organic superconductors for which there is nuclear magnetic resonance data, the energy scale for the spin fluctuations coincides with the energy scale for the pseudogap. This suggests that the pseudogap is caused by short-range spin correlations. In the weakly frustrated metals k-(BEDT-TTF)_2Cu[N(CN)_2]Br, k-(BEDT-TTF)_2Cu(NCS)_2, and k-(BEDT-TTF)_2Cu[N(CN)_2]Cl (under pressure) the pseudogap opens at the same temperature as coherence emerges in the (intralayer) transport. We argue that this is because the spin correlations are cut off by the loss of intralayer coherence at high temperatures. We discuss what might happen to these two energy scales at high pressures, where the electronic correlations are weaker. In these weakly frustrated materials the data is well described by the chemical pressure hypothesis (that anion substitution is equivalent to hydrostatic pressure). However, we find important differences in the metallic state of k-(BEDT-TTF)_2Cu_2(CN)_3, which is highly frustrated and displays a spin liquid insulating phase. We also show that the characteristic temperature scale of the spin fluctuations in (TMTSF)_2ClO_4 is the same as superconducting critical temperature, which may be evidence that spin fluctuations mediate the superconductivity in the Bechgaard salts.Comment: 7 pages, 4 figures; to appear in PR

Comparison of methods for estimating continuous distributions of relaxation times

The nonparametric estimation of the distribution of relaxation times approach is not as frequently used in the analysis of dispersed response of dielectric or conductive materials as are other immittance data analysis methods based on parametric curve fitting techniques. Nevertheless, such distributions can yield important information about the physical processes present in measured material. In this letter, we apply two quite different numerical inversion methods to estimate the distribution of relaxation times for glassy \lila\ dielectric frequency-response data at 225 \kelvin. Both methods yield unique distributions that agree very closely with the actual exact one accurately calculated from the corrected bulk-dispersion Kohlrausch model established independently by means of parametric data fit using the corrected modulus formalism method. The obtained distributions are also greatly superior to those estimated using approximate functions equations given in the literature.Comment: 4 pages and 4 figure

Oxygen Isotopic Imaging of Refractory Inclusions from the Miller Range (MIL) 090019 CO3 Chondrite: A Perovskite Perspective

Calcium-Aluminum-rich Inclusions (CAIs) in primitive meteorites are the first solids to condense in the Solar System. The oxygen isotopic compositions recorded in various mineral components of CAIs provide clues about their origins and post-formation histories, recording processes such as condensation, melting, nebular alteration, and fluidrock reactions on the parent body. MIL 090019 is similar to some rare carbonaceous chondrites such as Acfer 094, DOM 08004/6 and ALH 77303 that contain high abundances of a variety of refractory inclusions. This provides an opportunity to study the oxygen isotopic record of different types of refractory inclusions within the same meteorite. We analyzed CAIs specifically targeting primary minerals that are direct nebular condensates, such as corundum and perovskite, with the goal of gaining insights into the O isotopic composition of the nebular gas(es) from which these CAIs condensed. As MIL 090019 is a classified as CO3.1, it shows some signs of thermal metamorphism, compared to the more primitive CO3 meteorites (e.g., DOM 08004/06). A second goal of this study is to search for evidence of nebular processes in phases such as perovskite and melilite that are susceptible to parent body alteration to varying degrees. We analyzed the oxygen isotopic compositions of various CAIs from the MIL 090019 CO3 carbonaceous chondrite by ion imaging using the NanoSIMS 50L (Nano Secondary Ion Mass Spectrometer) at JSC following methods described in. An advantage of ion imaging over traditional spot analyses is that it provides spatial context to the oxygen isotopic data. This work builds on previously reported oxygen isotopic composition of two other CAIs (CAI-44 and CAI-E2) from the same meteorite thin section