SEARCH · Engineering Papers
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.
Full CI calculations have been performed on the CH3 radical. The full CI results are compared to those obtained using CASSCF/multireference CI and coupled-pair functional methods, both at the equilibrium CH distance and at geometries with the three CH bonds extended. In general, the performance of the approximate methods is similar to that observed in calculations on other molecules in which one or two bonds were stretched.
RE J1034+396 is one of the few active galactic nuclei (AGNs) with a significant quasiperiodic oscillation (QPO). The QPO has been observed in over 1 Ms of XMM-Newton observations spanning over a decade. We investigate the power spectral density function (PSD) of seven long (∼90 ks) XMM-Newton observations of the AGN RE J1034+396 in two energy bands. The soft (0.3–0.5 keV) band targets emission from the disk, while the hard (2–7 keV) band isolates the primary X-ray continuum emission from the corona. The QPO is significantly detected in the hard band of five of the seven observations. The best-fitting models indicate that the QPO detection in both bands is entirely attributable to the coronal emission with no additional contribution from the disk. This explains the strong coherence between the hard and soft bands at the QPO frequency. The covariance spectrum is consistent with this picture as the variability at QPO frequencies is attributed solely to fluctuations in the hot corona. The time lag as a function of energy is well described by a ∼2000 s intrinsic soft lag, resulting from the disk responding to emission from the corona, that undergoes phase wrapping at approximately the QPO frequency. By demonstrating that in this system the QPO arises in the corona, we provide new insights into the mechanisms generating QPOs.
Rock 71597, a coarse-grained plagioclase-poikilitic basalt that contains 15.8% MgO and 19.3% modal olivine, is described. Olivine and ilmenite crystals occur in two size generations, and large subequant olivine crystals, small 'matrix' olivines, and large ilmenite crystals are characterized. Since large olivine and ilmenite crystals with skeletal habits suggestive of crystallization under relatively rapid cooling are found in a rock whose texture is consonant with slower cooling, it is suggested that the large skeletal crystals formed toward the more rapidly cooled margin of a flow and accumulated by gravitational crystal settling toward the flow's interior. In support of this hypothesis, calculations indicate that the observed olivine compositions are approximately 10 mol% too Fe-rich to have crystallized from a melt of composition corresponding to that of 71597.
This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.
In 1952, Lighthill developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley revisited Proudman's work and applied the Lighthill relationship to evaluate directly the radiated acoustic power from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale based on Reynolds numbers 2.0 x 10(exp 3) less than or equal to R(sub lambda) less than or equal to 12.7 x 10(exp 3). Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 - 3.3 and the second order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68 - 3.68.
Nuclear energy is a highly desirable source of energy, and He-3 is the most prized of the fusion reactants. As the Wisconsin Group has emphasized, He-3 may be the only true economic ore on the Moon. The lack of a shielding atmosphere on the Moon permits solar-wind alpha particles to impinge upon the lunar regolith and become implanted into the various solid components. In particular, large quantities of helium (5 to 50 ppm) are presented. The measured parameter of I(sub s)/FeO, a direct indicator of maturity and exposure age, can be used as a first approximation to predict the abundances of many solar-wind components in the soils. However, because ilmenite has a much higher retentivity for helium than the other phases, the TiO2 contents of the soils are better indicators of helium contents (Taylor, Space 90). High-Ti mare bassalt regions, such as at the Apollo 17 locale, appear to be the best areas for He mining (15 to 50 ppm He(sub T)), versus 3 to 9 ppm in the Highlands. However, the relationships between I(sub s)/FeO, TiO2 and He-3 contents are complicated - e.g., many of the most He-rich soils are immature to submature. The amount of He-3 in the regolith of the moon is estimated at 220,000 tons in the outer 2 m of the Maria.
Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.
The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.
The paper concentrates on a linear approximation method for predicting the changes occurring in steady-state numerical solutions of the Euler equations as a consequence of small changes in the independent variables which control the problem. The importance of proper boundary-condition treatment and other issues concerning the problem are covered along with the importance of proper algorithm selection for a fully supersonic inviscid flow. The method is applied to a subsonic nozzle involving variation of the pressure on the outflow boundary and to a supersonic inlet involving variation of the inflow Mach number. In the subsonic test case, the comparisons between the predicted and conventional numerical solutions are shown to be good, while in the supersonic test case, the agreement between the approximation method and conventional numerical solution starts out well but rapidly degenerates at some point in the flowfield as the perturbation of the boundary conditions is increased.
We have studied the thermal consequences of very big impacts on Titan [1]. Titan's thick atmosphere and volatile-rich surface cause it to respond to big impacts in a somewhat Earth-like manner. Here we construct a simple globally-averaged model that tracks the flow of energy through the environment in the weeks, years, and millenia after a big comet strikes Titan. The model Titan is endowed with 1.4 bars of N2 and 0.07 bars of CH4, methane lakes, a water ice crust, and enough methane underground to saturate the regolith to the surface. We assume that half of the impact energy is immediately available to the atmosphere and surface while the other half is buried at the site of the crater and is unavailable on time scales of interest. The atmosphere and surface are treated as isothermal. We make the simplifying assumptions that the crust is everywhere as methane saturated as it was at the Huygens landing site, that the concentration of methane in the regolith is the same as it is at the surface, and that the crust is made of water ice. Heat flow into and out of the crust is approximated by step-functions. If the impact is great enough, ice melts. The meltwater oceans cool to the atmosphere conductively through an ice lid while at the base melting their way into the interior, driven down in part through Rayleigh-Taylor instabilities between the dense water and the warm ice. Topography, CO2, and hydrocarbons other than methane are ignored. Methane and ethane clathrate hydrates are discussed quantitatively but not fully incorporated into the model.
A number of systems were used to study relationships between lightning and other storm parameters. Instrumentation includes a system for locating cloud to ground lightning strike points, an acoustic lightning mapper and very high frequencies (VHF) (30 to 80 MHz) lightning mapper to study the structure and location of lightning, and S band radar to provide other storm data. Several characteristics of the reconstructed lightning have been relatively constant: (1) lightning structure is usually predominately horizontal; (2) lightning tends to occur preferrentially in certain altitude ranges that are only a few kilometers thick and above the O C isotherm; (3) it tends to occur near, but not necessarily inside, cores of high reflectivity (approximately or 45 dBz), and it often connects with regions of relatively weak reflectivity downstream from stronger reflectivity and updraft cores. Other relationships are also discussed.
The importance of blade design parameters in rotor vibratory response and the design of a minimum vibration blade based upon this understanding are examined. Various design approaches are examined for a 4 bladed articulated rotor operating at a high speed flight condition. Blade modal shaping, frequency placement, structural and aerodynamic coupling, and intermodal cancellation are investigated to systematically identify and evaluate blade design parameters that influence blade airloads, blade modal response, hub loads, and fuselage vibration. The relative contributions of the various components of blade force excitation and response to the vibratory hub loads transmitted to the fuselage are determined in order to isolate primary candidates for vibration alleviation. A blade design is achieved which reduces the predicted fuselage vibration from the baseline blade by approximately one half. Blade designs are developed that offer significant reductions in vibration (and fatigue stresses) without resorting to special vibration alleviation devices, radical blade geometries, or weight penalties.
Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.
Cyclotron radiation emission spectroscopy (CRES) is a modern technique for high-precision energy spectroscopy, in which the energy of a charged particle in a magnetic field is measured via the frequency of the emitted cyclotron radiation. The He6-CRES collaboration aims to use CRES to probe beyond the standard model physics at the TeV scale by performing high-resolution and low-background beta-decay spectroscopy of 6 He and 19 Ne. Having demonstrated the first observation of individual, high-energy (0.1–2.5 MeV) positrons and electrons via their cyclotron radiation, the experiment provides a novel window into the radiation of relativistic charged particles in a waveguide via the time-derivative (slope) of the cyclotron radiation frequency, df c /dt. We show that analytic predictions for the total cyclotron radiation power emitted by a charged particle in circular and rectangular waveguides are approximately consistent with the Larmor formula, each scaling with the Lorentz factor of the underlying e ± as γ 4 . This hypothesis is corroborated with experimental CRES slope data.
To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.
The surface of the planet Mars is characterized by two different hemispheres: old densely cratered structures in the Southern Uplands, and sparsely cratered younger topographically lower regions covering approximately the northern third of the planet's surface. The model for explaining these global surface structures is characterized by the following features: (1) homogeneous accretion of Mars; (2) formation of a metal melt layer; (3) northward translation of the undifferentiated core due to a Rayleigh-Taylor instability; (4) vigorous convection in the southern parts and formation of the Southern Uplands' crust; (5) fragmentation of the primordial core, slowly dissolving rockbergs, beginning of Tharsis uplift and volcanism; (6) formation of the Northern Lowlands' crust from only weakly differentiated silicatic material; and (7) reaching of the present state with symmetrically placed core and further thermal evolution.
Statistical estimates of wind shear in the planetary boundary layer are important in the design of V/STOL aircraft, and for the design of the Space Shuttle. The data analyzed in this study consist of eleven sets of longitudinal turbulent velocity fluctuation time histories digitized at 0.2 sec intervals with approximately 18,000 data points per time history. The longitudinal velocity fluctuations were calculated with horizontal wind and direction data collected at the 18-, 30-, 60-, 90-, 120-, and 150-m levels. The data obtained confirm the result that Eulerian time spectra transformed to wave-number spectra with Taylor's frozen eddy hypothesis possess inertial-like behavior at wave-numbers well out of the inertial subrange.
Spectra of Jupiter in the 5-micron window region have been acquired. From these it is deduced that (1) the presence of CO in Jupiter is verified; (2) the column abundance is approximately 0.016 cm amagat (identical to 4.3 by 10 to the 17th power mol/sq cm); and (3) the CO is probably nonuniformly mixed in the atmosphere, being concentrated into the stratosphere.