Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A coronal magnetic field model with horizontal volume and sheet currents

When globally mapping the observed photospheric magnetic field into the corona, the interaction of the solar wind and magnetic field has been treated either by imposing source surface boundary conditions that tacitly require volume currents outside the source surface or by limiting the interaction to thin current sheets between oppositely directed field regions. Yet observations and numerical Magnetohydrodynamic (MHD) calculations suggest the presence of non-force-free volume currents throughout the corona as well as thin current sheets in the neighborhoods of the interfaces between closed and open field lines or between oppositely directed open field lines surrounding coronal helmet-streamer structures. This work presents a model including both horizontal volume currents and streamer sheet currents. The present model builds on the magnetostatic equilibria developed by Bogdan and Low and the current-sheet modeling technique developed by Schatten. The calculation uses synoptic charts of the line-of-sight component of the photospheric magnetic field measured at the Wilcox Solar Observatory. Comparison of an MHD model with the calculated model results for the case of a dipole field and comparison of eclipse observations with calculations for CR 1647 (near solar minimum) show that this horizontal current-current-sheet model reproduces polar plumes and axes of corona streamers better than the source-surface model and reproduces polar plumes and axes of corona streamers better than the source-surface model and reproduces coro nal helmet structures better than the current-sheet model.

Zhao, Xuepu↗
NASA NTRS1994

A two-dimensional model study of the QBO signal in SAGE II NO2 and O3

Calculations of the quasi biennial oscillation (QBO) signal in Stratospheric Aerosol and Gas Experiment (SAGE) II O3 and NO2 data between 1984 and 1991 are presented and have been investigated by using a two-dimensional model. The isentropic 2D model is a fully interactive radiative-dynamical-chemical model in which the eddy fluxes of chemical species are calculated in a consistent manner. The QBO in the model has been forced by relaxing the equatorial zonal wind toward the observations at Singapore allowing the comparison of the model with observations from specific years. The model reproduces the observed vertical structure of the equatorial ozone anomaly with the well-known transition from dynamical to photochemical control at around 28km. The model also reproduces the observed vertical structure of the SAGE II observed NO2 anomaly. The model studies have shown that it is the QBO modulation of NO2 which the main cause of QBO signal in O3 above 30 km. The model also reproduces the observed latitudinal structure of the QBO signals in O3 and NO2. Due to the differing horizontal distribution of O3 and NO(y) the ozone signal shows a distinct phase change in the subtropics whereas the NO2 anomaly gives a broader signal.

Chipperfield, M. P.↗
NASA NTRS1994

Mirages and the nature of Pluto's atmosphere

We present model occultation lightcurves demonstrating that a strong thermal inversion layer at the base of Pluto's stratosphere can reproduce the minimum flux measured by the Kuiper Airborne Observatory (KAO) during the 1988 occultation of a star by Pluto. The inversion layer also forms the occultation equivalent of a mirage at a radius of 1198 km, which is capable of hiding tropospheres of significant depth. Pluto's surface lies below 1198 km, its radius depending on the depth of the troposphere. We begin by computing plausible temperature structures for Pluto's lower atmosphere, constrained by a calculation of the temperature of the atmosphere near the surface. We then trace rays from the occulted star through the model atmosphere, computing the resultant bending of the ray. Model light curves are obtained by summing the contribution of individual rays within the shadow of Pluto on Earth. We find that we can reproduce the KAO lightcurve using model atmospheres with a temperature inversion and no haze. We have explored models with tropospheres as deep as 40 km (implying a Pluto radius of 1158 km) that reproduce the suite of occultation data. Deeper tropospheres can be fitted to the data, but the mutual event radius of 1150 km probably provides a lower bound. If Pluto has a shallow or nonexistent troposphere, its density is consistent with formation in the solar nebula with modest water loss due to impact ejection. If the troposhere is relatively deep, implying a smaller radius and larger density, significant amounts of water loss are required.

Stansberry, J. A.↗
NASA NTRS1997

Development of an Antimicrobial Susceptibility Testing Method Suitable for Performing During Space Flight

Very little is known regarding the affects of the microgravity environment of space flight upon the action of antimicrobial agents on bacterial pathogens. This study was undertaken to develop a simple method for conducting antibacterial susceptibility tests during a Space Shuttle mission. Specially prepared susceptibility test research cards (bioMerieux Vitek, Hazelwood, MO) were designed to include 6-11 serial two-fold dilutions of 14 antimicrobial agents, including penicillins, cephalosporins, a Beta-lactamase inhibitor, vancomycin, erythromycin, tetracycline, gentamicin, ciprofloxacin, and trimethoprim/sulfamethoxazole. Minimal inhibitory concentrations (MICS) of the drugs were determined by visual reading of color endpoints in the Vitek research cards made possible by incorporation of a colorimetric growth indicator (alamarBlue(Trademark), Accumed International, Westlake, OH). This study has demonstrated reproducible susceptibility results when testing isolates of Staphylococcus aurezis, Group A Streptococcus, Enterococcusfaecalis, Escherichia coli (beta-lactamase positive and negative strains), Klebsiella pneumoniae, Enterobacter cloacae, and Pseudomoiias aeruginosa. In some instances, the MICs were comparable to those determined using a standard broth microdilution method, while in some cases the unique test media and format yielded slightly different values, that were themselves reproducible. The proposed in-flight experiment will include inoculation of the Vitek cards on the ground prior to launch of the Space Shuttle, storage of inoculated cards at refrigeration temperature aboard the Space Shuttle until experiment initiation, then incubation of the cards for 18-48 h prior to visual interpretation of MICs by the mission's astronauts. Ground-based studies have shown reproducible MICs following storage of inoculated cards for 7 days at 4-8 C to accommodate the mission's time schedule and the astronauts' activities. For comparison, ground-based control (normal gravity) MIC values will be generated by simultaneous inoculation and incubation of a second set of test cards in a laboratory at the launch site. This procedure can provide a safe and compact experiment that should yield new information on the affects of microgravity on the biological activities of various classes of antibiotics.

Jorgensen, James H.↗
NASA NTRS1999

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗
NASA NTRS1997

Comparison of TOPEX/Poseidon Sea Level and Linear Model Results forced by Various Wind Products for the Tropical Pacific

The goal of this paper is to compare TOPEX/Posaidon (T/P) sea level with sea level results from linear ocean model experiments forced by several different wind products for the tropical Pacific. During the period of this study (October 1992 - October 1995), available wind products include satellite winds from the ERS-1 scatterometer product of [HALP 97] and the passive microwave analysis of SSMI winds produced using the variational analysis method (VAM) of [ATLA 91]. In addition, atmospheric GCM winds from the NCEP reanalysis [KALN 96], ECMWF analysis [ECMW94], and the Goddard EOS-1 (GEOS-1) reanalysis experiment [SCHU 93] are available for comparison. The observed ship wind analysis of FSU [STRI 92] is also included in this study. The linear model of [CANE 84] is used as a transfer function to test the quality of each of these wind products for the tropical Pacific. The various wind products are judged by comparing the wind-forced model sea level results against the T/P sea level anomalies. Correlation and RMS difference maps show how well each wind product does in reproducing the T/P sea level signal. These results are summarized in a table showing area average correlations and RMS differences. The large-scale low-frequency temporal signal is reproduced by all of the wind products, However, significant differences exist in both amplitude and phase on regional scales. In general, the model results forced by satellite winds do a better job reproducing the T/P signal (i.e. have a higher average correlation and lower RMS difference) than the results forced by atmospheric model winds.

Hackert, Eric C.↗
NASA NTRS2000

Modeling the Diurnal Tides in the MLT Region with the Doppler Spread Parameterization of Gravity Waves

Special Session: SA03 The mesosphere/lower thermosphere region: Structure, dynamics, composition, and emission. Ground based and satellite observations in the upper mesosphere and lower thermosphere (MLT) reveal large seasonal variations in the horizontal wind fields of the diurnal and semidiurnal tides. To provide an understanding of the observations, we discuss results obtained with our Numerical Spectral Model (NMS) that incorporates the gravity wave Doppler Spread Parameterization (DSP) of Hines. Our model reproduces many of the salient features observed, and we discuss numerical experiments that delineate the important processes involved. Gravity wave momentum deposition and the seasonal variations in the tidal excitation contribute primarily to produce the large equinoctial amplitude maxima in the diurnal tide. Gravity wave induced variations in eddy viscosity, not accounted for in the model, have been shown by Akmaev to be important too. For the semidiurnal tide, with amplitude maximum observed during winter solstice, these processes also contribute, but filtering by the mean zonal circulation is more important. A deficiency of our model is that it cannot reproduce the observed seasonal variations in the phase of the semidiurnal tide, and numerical experiments are being carried out to diagnose the cause and to alleviate this problem. The dynamical components of the upper mesosphere are tightly coupled through non-linear processes and wave filtering, and this may constrain the model and require it to reproduce in detail the observed phenomenology.

Mayr, H. G.↗
NASA NTRS2002

Who or What? Self-Replication and Function-Reproduction in the Origin of Life

In this presentation, we will present results on the fundamental properties of two classes of replicating systems: autocatalytic replicators that reproduce exact copies of a template molecule, and function reproducers that maintain a set of essential functions without replicating the identities of the functional moieties. We will describe the stability and behavior in-the-large of autocatalytic replicators. Most importantly, we have found no sharp distinction between an autocatalytic and a non-autocatalytic domain. We will also present a new derivation of von Kiedrowski's square-root rate law. Function - reproducers are proposed as an important component of protocells and we will present theoretical results on a simple model system that incorporates known peptide biophysics. For a wide range of parameters, we have shown that this type of system can improve its overall performance, even in the absence of any method for information storage. This type of system improvement is defined to be non-genomic evolution.

New, Michael H.↗
NASA NTRS2004

Decelerating Flows in the TeV Blazars

TeV emission from a class of BL Lacertae (BL) objects is commonly modeled as radiation from relativistically moving homogeneous plasma blobs. In the context of these models, the blob Lorentz factors needed to reproduce the (corrected for absorption by the IR background) TeV emission are large ($\delta \gtrsim 50$) are required to reproduce via Synchrotron-Self Compton (SSC) the observed TeV emission. The main reason for this is that stronger beaming eases the problem of the lack of $\sim$ IR-UV synchrotron seed photons needed to produce the de-absorbed $\sim$ few TeV peak of the spectral energy distribution (SED). However, such high Doppler factors are in strong disagreement with the unified scheme, according to which BLs are FR I radio galaxies with their jets closely aligned to the line of sight. Here, motivated by the detection of sub-luminal velocities in the sub-pc scale jets of the best studied TeV blazars, MKN 421 and MKN 501. we examine the possibility that the relativistic flow in the TeV BLs is longitudinally decelerating. In this case, the problem of the missing seed photons is solved due to Upstream Compton (UC) scattering, a process in which the upstream energetic electrons from the fast base of the flow 'see' the synchrotron seed photons produced in the slow part of the flow relativistically beamed. Modest Lorentz factors ($\Gamma kim 15s). decelerating down to values compatible with the recent radio interferometric observations, reproduce the $\sim$ few TeV peak energy of these sources. Furthermore, such decelerating flows are shown to be in agreement with the BL - FR I unification.

Kazanas, Demosthenes↗
NASA NTRS2005

Construction of Covariance Functions with Variable Length Fields

This article focuses on construction, directly in physical space, of three-dimensional covariance functions parametrized by a tunable length field, and on an application of this theory to reproduce the Quasi-Biennial Oscillation (QBO) in the Goddard Earth Observing System, Version 4 (GEOS-4) data assimilation system. These Covariance models are referred to as multi-level or nonseparable, to associate them with the application where a multi-level covariance with a large troposphere to stratosphere length field gradient is used to reproduce the QBO from sparse radiosonde observations in the tropical lower stratosphere. The multi-level covariance functions extend well-known single level covariance functions depending only on a length scale. Generalizations of the first- and third-order autoregressive covariances in three dimensions are given, providing multi-level covariances with zero and three derivatives at zero separation, respectively. Multi-level piecewise rational covariances with two continuous derivatives at zero separation are also provided. Multi-level powerlaw covariances are constructed with continuous derivatives of all orders. Additional multi-level covariance functions are constructed using the Schur product of single and multi-level covariance functions. A multi-level powerlaw covariance used to reproduce the QBO in GEOS-4 is described along with details of the assimilation experiments. The new covariance model is shown to represent the vertical wind shear associated with the QBO much more effectively than in the baseline GEOS-4 system.

Gaspari, Gregory↗
NASA NTRS2005

Decelerating Flows in TeV Blazars: A Resolution to the BL Lac -- FR I Unification Problem

TeV emission from BL Lacertae (BL) objects i s commonly modeled as Synchrotron-Self Compton (SSC) radiation from relativistically moving homogeneous plasma blobs. In the context of these models, the blob Lorentz factors needed to reproduce the corrected for absorption by the diffuse IR background (DIRB) TeV emission are large ($\delta \gtrsim 50$). The main reason for this is that stronger beaming eases the problem of the lack of $\sim$ IR-UV synchrotron seed photons needed to produce the de-absorbed $\sim $ few TeV peak of the spectral energy distribution (SED). However, such high Doppler factors are in strong disagreement with the unified scheme, according to which BLs are FR I radio galaxies with their jets closely aligned to the line of sight. Here, motivated by the detection of sub-luminal velocities in the sub-pc scale jets of the best studied TeV blazars, MKN 421 and MKN 501, we examine the possibility that the relativistic flows in the TeV BLs decelerate. In this case, the problem of the missing seed photons is solved because of Upstream Compton (UC) scattering, a process in which the upstream energetic electrons from the fast base of the flow 'see' the synchrotron seed photons produced in the slow part of the flow relativistically beamed. Modest Lorentz factors ($\Gamma \sim 15$), decelerating down to values compatible with the recent radio interferometric observations, reproduce the $\sim $ few TeV peak energy of these sources. Furthermore, such decelerating flows are shown to be in agreement with the BL - FR I unification, naturally reproducing the observed BL/FR I broadband luminosity ratios.

Kazanas, Demosthenes↗
NASA NTRS2000

Creep Properties of NiAl-1Hf Single Crystals Re-Investigated

NiAl-1Hf single crystals have been shown to be quite strong at 1027 C, with strength levels approaching those of advanced Ni-based superalloys. Initial testing, however, indicated that the properties might not be reproducible. Study of the 1027 C creep behavior of four different NiAl-1Hf single-crystal ingots subjected to several different heat treatments indicated that strength lies in a narrow band. Thus, we concluded that the mechanical properties are reproducible. Recent investigations of the intermetallic NiAl have confirmed that minor alloying additions combined with single-crystal growth technology can produce elevated temperature strength levels approaching those of Ni-based superalloys. For example, General Electric alloy AFN 12 {Ni-48.5(at.%) Al-0.5Hf-1Ti-0.05Ga} has a creep rupture strength equivalent to Rene 80 combined with a approximately 30-percent lower density, a fourfold improvement in thermal conductivity, and the ability to form a self-protective alumina scale in aggressive environments. Although the compositions of strong NiAl single crystals are relatively simple, the microstructures are complex and vary with the heat treatment and with small ingot-toingot variations in the alloy chemistry. In addition, initial testing suggested a strong dependence between microstructure and creep strength. If these observations were true, the ability to utilize NiAl single-crystal rotating components in turbine machinery could be severely limited. To investigate the possible limitations in the creep response of high-strength NiAl single crystals, the NASA Glenn Research Center at Lewis Field initiated an in depth investigation of the effect of heat treatment on the microstructure and subsequent 1027 C creep behavior of [001]-oriented NiAl-1Hf with a nominal chemistry of Ni-47.5Al-1Hf-0.5Si. This alloy was selected since four ingots, grown over a number of years and possessing slightly different compositions, were available for study. Specimens taken from the ingots were subjected to several heat treatment schedules, examined by transmission electron microscopy, and tested in both compression and tension. An example of the microstructure found in a [001]-oriented NiAl-1Hf specimen after a solution treatment at 1317 C for 50 hr followed by air cooling is illustrated in the image on the left, where the NiAl matrix contains a uniform distribution of nanometer-scale Gphase (Ni16Hf6Si7) precipitates. Other heat treating schedules produced microstructures with nanometer-sized G-phase cubes and plates or, in an extreme case, produced a microstructure with all the G-phase converted to Heusler (Ni2AlHf) particles. The results of 1027 C creep strength and strain rate testing are illustrated which summarizes data from tensile and compressive testing of samples cut from all four NiAl-1Hf ingots and subjected to a variety of heat treatment schedules. With one exception, all the strength values lie in a narrow band that spans six orders of magnitude in strain rate. The only factor that produced results outside of this band was the heat treatment schedule that dissolved all the G-phase and replaced it with Heusler precipitates. The results portrayed in this figure lead to the important practical conclusion that the elevated-temperature creep properties of NiAl-1Hf single crystals are reproducible and are not affected by small variations in alloy chemistry from ingot to ingot or by different initial distributions of G-phase in the heat-treated alloy. The only variable in this study that produced a significant and delerious effect on mechanical strength was a post-solution heat treatment that lead to the complete disappearance of the G-phase in favor of Heusler precipitates.

Whittenberger, J. Daniel↗
NASA NTRS2005

Simulations of Dynamics and Transport during the September 2002 Antarctic Major Warming

A mechanistic model simulation initialized on 14 September 2002, forced by 100-hPa geopotential heights from Met Office analyses, reproduced the dynamical features of the 2002 Antarctic major warming. The vortex split on approx.25 September; recovery after the warming, westward and equatorward tilting vortices, and strong baroclinic zones in temperature associated with a dipole pattern of upward and downward vertical velocities were all captured in the simulation. Model results and analyses show a pattern of strong upward wave propagation throughout the warming, with zonal wind deceleration throughout the stratosphere at high latitudes before the vortex split, continuing in the middle and upper stratosphere and spreading to lower latitudes after the split. Three-dimensional Eliassen-Palm fluxes show the largest upward and poleward wave propagation in the 0(deg)-90(deg)E sector prior to the vortex split (coincident with the location of strongest cyclogenesis at the model's lower boundary), with an additional region of strong upward propagation developing near 180(deg)-270(deg)E. These characteristics are similar to those of Arctic wave-2 major warmings, except that during this warming, the vortex did not split below approx.600 K. The effects of poleward transport and mixing dominate modeled trace gas evolution through most of the mid- to high-latitude stratosphere, with a core region in the lower-stratospheric vortex where enhanced descent dominates and the vortex remains isolated. Strongly tilted vortices led to low-latitude air overlying vortex air, resulting in highly unusual trace gas profiles. Simulations driven with several meteorological datasets reproduced the major warming, but in others, stronger latitudinal gradients at high latitudes at the model boundary resulted in simulations without a complete vortex split in the midstratosphere. Numerous tests indicate very high sensitivity to the boundary fields, especially the wave-2 amplitude. Major warmings occurred for initial fields with stronger winds and larger vortices, but not smaller vortices, consistent with the initiation of wind-deceleration by upward-propagating waves near the poleward edge of the region where wave 2 can propagate above the jet core. Thus, given the observed 100-hPa boundary forcing, stratospheric preconditioning is not needed to reproduce a major warming similar to that observed. The anomalously strong forcing in the lower stratosphere can be viewed as the primary direct cause of the major warming.

winds↗
NASA NTRS2008

Quantitative Studies of the Optical and UV Spectra of Galactic Early B Supergiants

We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic B0-B5 supergiants of luminosity classes Ia, Ib, Iab, and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model predictions to observed UV spectral features is made. Methods. Fundamental parameters (e.g., T(sub eff), log L(sub *), mass-loss rates and CNO abundances) are derived for individual stars using CMFGEN, a nLTE, line-blanketed model atmosphere code. The impact of these newly derived parameters on the Galactic B supergiant Ten scale, mass discrepancy, and wind-momentum luminosity relation is examined. Results. The B supergiant temperature scale derived here shows a reduction of about 1000-3000 K compared to previous results using unblanketed codes. Mass-loss rate estimates are in good agreement with predicted theoretical values, and all of the 20 BO-B5 supergiants analysed show evidence of CNO processing. A mass discrepancy still exists between spectroscopic and evolutionary masses, with the largest discrepancy occuring at log (L/(solar)L approx. 5.4. The observed WLR values calculated for B0-B0.7 supergiants are higher than predicted values, whereas the reverse is true for B1-B5 supergiants. This means that the discrepancy between observed and theoretical values cannot be resolved by adopting clumped (i.e., lower) mass-loss rates as for O stars. The most surprising result is that, although CMFGEN succeeds in reproducing the optical stellar spectrum accurately, it fails to precisely reproduce key UV diagnostics, such as the N v and C IV P Cygni profiles. This problem arises because the models are not ionised enough and fail to reproduce the full extent of the observed absorption trough of the P Cygni profiles. Conclusions. Newly-derived fundamental parameters for early B supergiants are in good agreement with similar work in the field. The most significant discovery, however, is the failure of CMFGEN to predict the correct ionisation fraction for some ions. Such findings add further support to revising the current standard model of massive star winds, as our understanding of these winds is incomplete without a precise knowledge of the ionisation structure and distribution of clumping in the wind. Key words. techniques: spectroscopic - stars: mass-loss - stars: supergiants - stars: abundances - stars: atmospheres - stars: fundamental parameters

Searle, S. C.↗
NASA NTRS2007

Anatomy of the Binary Black Hole Recoil: A Multipolar Analysis

We present a multipolar analysis of the recoil velocity computed in recent numerical simulations of binary black hole coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including 1 = 4 are sufficient to accurately reproduce the final recoil velocity (= 98%) and that only a few dominant modes contribute significantly to it (2 95%). We describe how the relative amplitude, and more importantly, the relative phase, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ring-down phases. We also find that the numerical results can be reproduced, to a high level of accuracy, by an effective Newtonian formula for the multipole moments obtained by replacing in the Newtonian formula the radial separation with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes. Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of '.anti-kick" for each of the simulations. Lastly, we apply this multipolar analysis to understand the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Schnittman, Jeremy↗