Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 613 records · Page 34

Study of pre-storm environment by using rawinsonde and satellite observations

Rawinsonde and satellite remote sensing data were utilized to examine the prestorm environment and mechanisms for the initiation of four groups of severe storms. The storms in Altus, Oklahoma, Pampas, Texas, Bennett, Colorado, and Red River Valley, Oklahoma are described. The geographical distributions of the areas of high moisture concentration and variations of tropopause heights for the storm groups are analyzed. It is detected that in the area of a low-level high concentration of moisture, the local tropopause height is lowest at the time of the storm cloud formation and development, and the potential energy storage per unit areas for the overshootiong clouds penetrating above the tropopause is related to the intensity of the storms produced. Numerical cloud modeling was performed for the storms. The model data are compared with the satellite and rawinsonde observations, and it is noted that the data correlate well.

Hung, R. J.↗

Repeated buckling of composite shear panels

Failures in service of aerospace structures and research at the Technion Aircraft Structures Laboratory have revealed that repeatedly buckled stiffened shear panels might be susceptible to premature fatigue failures. Extensive experimental and analytical studies have been performed at Technion on repeated buckling, far in excess of initial buckling, for both metal and composite shear panels with focus on the influence of the surrounding structure. The core of the experimental investigation consisted of repeated buckling and postbuckling tests on Wagner beams in a three-point loading system under realistic test conditions. The effects of varying sizes of stiffeners, of the magnitude of initial buckling loads, of the panel aspect ratio and of the cyclic shearing force, V sub cyc, were studied. The cyclic to critical shear buckling ratios, (V sub cyc/V sub cr) were on the high side, as needed for efficient panel design, yet all within possible flight envelopes. The experiments were supplemented by analytical and numerical analyses. For the metal shear panels the test and numerical results were synthesized into prediction formulas, which relate the life of the metal shear panels to two cyclic load parameters. The composite shear panels studied were hybrid beams with graphite/epoxy webs bonded to aluminum alloy frames. The test results demonstrated that composite panels were less fatigue sensitive than comparable metal ones, and that repeated buckling, even when causing extensive damage, did not reduce the residual strength by more than 20 percent. All the composite panels sustained the specified fatigue life of 250,000 cycles. The effect of local unstiffened holes on the durability of repeatedly buckled shear panels was studied for one series of the metal panels. Tests on 2024 T3 aluminum panels with relatively small unstiffened holes in the center of the panels demonstrated premature fatigue failure, compared to panels without holes. Preliminary tests on two graphite epoxy shear panels with small holes in the center showed no similar fatigue life degradation and no shift in failure mode. Further tests on the effect of holes are in progress.

Singer, Josef↗

Theory of finite disturbances in a centrifugal compression system with a vaneless radial diffuser

A previous small perturbation analysis of circumferential waves in circumferential compression systems, assuming inviscid flow, is shown to be consistent with observations that narrow diffusers are more stable than wide ones, when boundary layer displacement effect is included. The Moore-Greitzer analysis for finite strength transients containing both surge and rotating stall in axial machines is adapted for a centrifugal compression system. Under certain assumptions, and except for a new second order swirl, the diffuser velocity field, including resonant singularities, can be carried over from the previous inviscid linear analysis. Nonlinear transient equations are derived and applied in a simple example to show that throttling through a resonant value of flow coefficient must occur in a sudden surge-like drop, accompanied by a transient rotating wave. This inner solution is superseded by an outer surge response on a longer time scale. Surge may occur purely as result of circumferential wave resonance. Numerical results are shown for various parametric choices relating to throttle schedule and the characteristic slope. A number of circumferential modes considered simultaneously is briefly discussed.

Moore, F. K.↗

Cosmic gamma-ray bursts

A review of the cosmic gamma-ray burst phenomenon is presented. Both the light curves and the energy spectra of these short transient events display a great diversity. However, rapid rise times and periodicities sometimes observed in the light curves suggest a compact object origin. Similarly, absorption and emission features in the energy spectra argue strongly in favor of this interpretation. Counterparts to gamma-bursters in other energy ranges, such as optical and sort x-ray, have still not been identified, however, leading to a large uncertainty in the distances to bursters. Although gamma-ray burst sources have not yet been observed to repeat, numerous bursts from three objects which may be related to the gamma-bursters, called Soft Gamma Repeaters, have been recorded; there is weak evidence that they may be relatively distant on a galactic scale. Future missions, particularly those emphasizing high energy, time, and/or spatial resolution, as well as a multiwavelength approach, are likely to advance our understanding of this enigmatic phenomenon.

Hurley, K.↗

Ordinal judgments of numerical symbols by macaques (Macaca mulatta)

Two rhesus monkeys (Macaca mulatta) learned that the arabic numerals 0 through 9 represented corresponding quantities of food pellets. By manipulating a joystick, the monkeys were able to make a selection of paired numerals presented on a computer screen. Although the monkeys received a corresponding number of pellets even if the lesser of the two numerals was selected, they learned generally to choose the numeral of greatest value even when pellet delivery was made arrhythmic. In subsequent tests, they chose the numerals of greater value when presented in novel combinations or in random arrays of up to five numerals. Thus, the monkeys made ordinal judgments of numerical symbols in accordance with their absolute or relative values.

Washburn, David A.↗

The relative importance of mass and wind data in the FGGE observing system

There are two theoretical arguments that have been used to discuss the relative importance of mass and wind data in numerical weather prediction (NWP). We will analyze these arguments in this section as clearly as possible in order to draw conclusions which may help to interpret experimental results on four-dimensional data assimilation, simulations of future observing systems, as well as give guidance on how to improve the efficiency with which we use the present observing system.

Kalnay, E.↗

A Two Element Laminar Flow Airfoil Optimized for Cruise

Numerical and experimental results are presented for a new two-element, fixed-geometry natural laminar flow airfoil optimized for cruise Reynolds numbers on the order of three million. The airfoil design consists of a primary element and an independent secondary element with a primary to secondary chord ratio of three to one. The airfoil was designed to improve the cruise lift-to-drag ratio while maintaining an appropriate landing capability when compared to conventional airfoils. The airfoil was numerically developed utilizing the NASA Langley Multi-Component Airfoil Analysis computer code running on a personal computer. Numerical results show a nearly 11.75 percent decrease in overall wing drag with no increase in stall speed at sailplane cruise conditions when compared to a wing based on an efficient single element airfoil. Section surface pressure, wake survey, transition location, and flow visualization results were obtained in the Texas A&M University Low Speed Wind Tunnel. Comparisons between the numerical and experimental data, the effects of the relative position and angle of the two elements, and Reynolds number variations from 8 x 10(exp 5) to 3 x 10(exp 6) for the optimum geometry case are presented.

Steen, Gregory Glen↗

[Activities of the Department of Electrical Engineering, Howard University]

Theoretical derivations, computer analysis and test data are provided to demonstrate that the cavity model is a feasible one to analyze thin-substrate, rectangular-patch microstrip antennas. Seven separate antennas were tested. Most of the antennas were designed to resonate at L-band frequencies (1-2 GHz). One antenna was designed to resonate at an S-band (2-4 GHz) frequency of 2.025 GHz. All dielectric substrates were made of Duroid, and were of varying thicknesses and relative dielectric constant values. Theoretical derivations to calculate radiated free space electromagnetic fields and antenna input impedance were performed. MATHEMATICA 2.2 software was used to generate Smith Chart input impedance plots, normalized relative power radiation plots and to perform other numerical manipulations. Network Analyzer tests were used to verify the data from the computer programming (such as input impedance and VSWR). Finally, tests were performed in an anechoic chamber to measure receive-mode polar power patterns in the E and H planes. Agreement between computer analysis and test data is presented. The antenna with the thickest substrate (e(sub r) = 2.33,62 mils thick) showed the worst match to theoretical impedance data. This is anticipated due to the fact that the cavity model generally loses accuracy when the dielectric substrate thickness exceeds 5% of the antenna's free space wavelength. A method of reducing computer execution time for impedance calculations is also presented.

Yalamanchili, Raj C.↗

The Uncertain Nature of Cometary Motions

The number of active short- and long-periodic comets crossing the Earth's orbit each year is less than 10 percent of the corresponding number of asteroids crossing the Earth's orbit. However, the higher relative velocities of comets with respect to the Earth and the uncertainties associated with accurately computing their future trajectories can cause considerable problems when assessing the risks of Earth-crossing objects. Unlike asteroids, the motions of active comets are often affected by so-called nongravitational (outgassing) forces that are imperfectly modeled. In addition, the astrometric optical observations that are used to refine a comet's orbit are often imprecise because a comet's center of mass can be hidden by atmospheric gas and dust. For long-period comets, there is the additional problem of having to base orbital solutions on relatively short observational data intervals. Long-term numerical integrations extending two centuries into the future have been carried out to investigate upcoming Earth-close approaches by known periodic comets. Error analyses and impact probabilities have been computed for those comets that will pass closest to the Earth. Although there are no known comets that will make dangerously close Earth approaches in the next two centuries, there are a few objects that warrant future monitoring.

Yeomans, Donald K.↗

Interannual Fluctuations in Atmospheric Angular Momentum Simulated by the National Centers for Environmental Prediction Medium Range Forecast Model

An earlier study established the existence of globally coherent interannual fluctuations in atmospheric angular momentum (AAM), associated with the El Nino-Southern Oscillation (ENSO) cycle. In this paper, we pursue the origin and the structure of these fluctuations using an ensemble of experiments generated by the National Centers for Environmental Prediction, medium range forecast model version 9. In the control experiments, where the observed sea surface temperatures (SSTs) were used as the lower boundary conditions, the model captures the characteristic V-like structure in time-latitude plots of zonally averaged AAM, while experiments with climatological SSTs and those with either perpetual warm or cold ENSO conditions superimposed on the climatological SSTs failed to reproduce this structure. The numerical results indicate that these AAM structures are related to SST variations associated with transitions between different phases of the ENSO cycle and have both propagating and standing components. The largest zonal wind contribution from the levels studied (850, 500, and 200 hPa) is at 200 hPa, where the tropical convective outflow is the strongest. Composites of zonal wind and geopotential height show a clear relationship between the stages of the global AAM oscillation and the ENSO cycle. The strong similarity between the simulated and observed AAM series attests to the model's ability to realistically simulate the interannual response of the atmosphere to ENSO SST anomalies.

Mo, Kingtse C.↗

Exchange Forces in Dispersion Relations Investigated Using Circuit Relations

We propose a novel method to compute in an exact manner the left-hand cut discontinuity of the electron-Atom partial wave scattering amplitude in the complex energy plane within the static stage approximation. Zero energy dispersion relations for electron-Hydrogen scattering are computed numerically for illustration.

Vrinceanu, D.↗

Hyperbolic Methods for Surface and Field Grid Generation

This chapter describes the use of hyperbolic partial differential equation methods for structured surface grid generation and field grid generation. While the surface grid generation equations are inherently three dimensional, the field grid generation equations can be formulated in two or three dimensions. The governing equations are derived from orthogonality relations and cell area/volume constraints; and are solved numerically by marching from an initial curve or surface. The marching step size and marching distance can be prescribedly the user. Exact specifications of the side and outer boundaries are not possible with a one sweep marching scheme but limited control is achievable. Excellent orthogonality and grid clustering characteristics are provided by hyperbolic methods with one to two orders of magnitude savings in time over typical elliptic methods. Since hyperbolic grid generation methods do not require the exact specifications of the side and outer boundaries of a grid, these methods are particularly well suited for the overlapping grid approach for solving problems on complex configurations. Grid generation software based on hyperbolic methods and their applications on several complex configurations will be described.

Chan, William M.↗

Self-Scheduling Parallel Methods for Multiple Serial Codes with Application to WOPWOP

This paper presents a scheme for efficiently running a large number of serial jobs on parallel computers. Two examples are given of computer programs that run relatively quickly, but often they must be run numerous times to obtain all the results needed. It is very common in science and engineering to have codes that are not massive computing challenges in themselves, but due to the number of instances that must be run, they do become large-scale computing problems. The two examples given here represent common problems in aerospace engineering: aerodynamic panel methods and aeroacoustic integral methods. The first example simply solves many systems of linear equations. This is representative of an aerodynamic panel code where someone would like to solve for numerous angles of attack. The complete code for this first example is included in the appendix so that it can be readily used by others as a template. The second example is an aeroacoustics code (WOPWOP) that solves the Ffowcs Williams Hawkings equation to predict the far-field sound due to rotating blades. In this example, one quite often needs to compute the sound at numerous observer locations, hence parallelization is utilized to automate the noise computation for a large number of observers.

Long, Lyle N.↗

Numerical Study of Wake Vortex Behavior in Turbulent Domains with Ambient Stratification

A three-dimensional large eddy simulation model is used to investigate the sensitivity of ambient stratification with turbulence on the behavior of aircraft wake vortices. Modeled ambient turbulence levels range from very weak to moderate, and stratification levels range from strongly stable to unstable. The results of profound significance from this study are: 1) very little sensitivity between vortex linking time and the level of stratification, 2) the mean vortex separation remained nearly constant regardless of stratification and turbulence (at least prior to linking), 3) the wake vortices did not rise regardless of the level of stratification, and 4) for very strong stratification, the vortex stopped descending and quickly dissipated even before vortex linking could occur. These results are supported by experimental data and are contrary to conclusions from other numerical studies that assume laminar flow and/or relatively-low Reynolds numbers.

Switzer, George F.↗

EARLY IMPACT MELTING AND SPACE EXPOSURE HISTORY OF THE PAT91501 LCHONDRITE

Collisions probably occurred frequently in the early history of the asteroid belt. Their effects, which should be recorded in meteorites, must have included heating and melting along with shock alteration of mineral textures. Some non-chondritic meteorite types e.g., eucrites and IIE and IAB irons - do indeed give evidence of extensive impact heating more than 3.4 Gyr ago. The ordinary chondrites, in contrast, show little evidence of early impact heating. The Ar-Ar and Rb-Sr ages of ordinary chondrites that experienced intense shock are for the most part relatively young, many less than 1.5 Gyr. The numerous L-chondrites with Ar- Ar ages clustering near 0.5 Gy are a well-known example. One of them, the 105-kg Chico Lchondrite, shows the effects of unusually intense heating. It is approximately 60% impact melt and likely formed as a dyke beneath a large crater when the L-chondrite parent body underwent a very large impact approximately 0.5 Gyr ago. In rare instances, older shock dates are indicated for ordinary chondrites. Dixon et al show early impact resetting of Ar-Ar ages of a few LL-chondrites including MIL 99301 at 4.23 0.03 Gyr, but in none of these stones did shock lead to extensive melting. As of 2003, searches for chondritic melts attributable to early shock had turned up only the Shaw L-chondrite, which has an Ar-Ar age of approximately 4.42 Gyr. PAT91501 is an 8.55-kg L-chondrite containing vesicles and metal-troilite nodules. It is a unique, near-total impact melt, unshocked, depleted in siderophile and chalcophile elements, and contains only approximately 10% relic chondritic material. The authors conclude that PAT91501 crystallized rapidly and from a much more homogeneous melt than did Shaw. They suggest that PAT resembles Chico and likely formed as an impact melt vein within an impact crater. To define the history of PAT, we have determined its Ar-39-Ar-40 age and measured several radioactive and stable nuclides produced during its space exposure to cosmic rays.

Bogard, Donald D.↗

Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997-2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the 'truth' signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.

Combination↗

Shear Thinning Near the Critical Point of Xenon

We measured shear thinning, a viscosity decrease ordinarily associated with complex liquids, near the critical point of xenon. The data span a wide range of reduced shear rate: 10(exp -3) < gamma-dot tau < 700, where gamma-dot tau is the shear rate scaled by the relaxation time tau of critical fluctuations. The measurements had a temperature resolution of 0.01 mK and were conducted in microgravity aboard the Space Shuttle Columbia to avoid the density stratification caused by Earth's gravity. The viscometer measured the drag on a delicate nickel screen as it oscillated in the xenon at amplitudes 3 mu,m < chi (sub 0) >430 mu, and frequencies 1 Hz < omega/2 pi < 5 Hz. To separate shear thinning from other nonlinearities, we computed the ratio of the viscous force on the screen at gamma-dot tau to the force at gamma-dot tau approximates 0: C(sub gamma) is identical with F(chi(sub 0), omega tau, gamma-dot tau )/F)(chi(sub 0, omega tau, 0). At low frequencies, (omega tau)(exp 2) < gamma-dot tau, C(sub gamma) depends only on gamma-dot tau, as predicted by dynamic critical scaling. At high frequencies, (omega tau)(exp 2) > gamma-dot tau, C(sub gamma) depends also on both x(sub 0) and omega. The data were compared with numerical calculations based on the Carreau-Yasuda relation for complex fluids: eta(gamma-dot)/eta(0)=[1+A(sub gamma)|gamma-dot tau|](exp - chi(sub eta)/3+chi(sub eta)), where chi(sub eta) =0.069 is the critical exponent for viscosity and mode-coupling theory predicts A(sub gamma) =0.121. For xenon we find A(sub gamma) =0.137 +/- 0.029, in agreement with the mode coupling value. Remarkably, the xenon data close to the critical temperature T(sub c) were independent of the cooling rate (both above and below T(sub c) and these data were symmetric about T(sub c) to within a temperature scale factor. The scale factors for the magnitude of the oscillator s response differed from those for the oscillator's phase; this suggests that the surface tension of the two-phase domains affected the drag on the screen below T(sub c).

Zimmerli, Gregory A.↗

Examining the Impacts of High-Resolution Land Surface Initialization on Model Predictions of Convection in the Southeastern U.S.

One of the most challenging weather forecast problems in the southeastern U.S. is daily summertime pulse convection. During the summer, atmospheric flow and forcing are generally weak in this region; thus, convection typically initiates in response to local forcing along sea/lake breezes, and other discontinuities often related to horizontal gradients in surface heating rates. Numerical simulations of pulse convection usually have low skill, even in local predictions at high resolution, due to the inherent chaotic nature of these precipitation systems. Forecast errors can arise from assumptions within physics parameterizations, model resolution limitations, as well as uncertainties in both the initial state of the atmosphere and land surface variables such as soil moisture and temperature. For this study, it is hypothesized that high-resolution, consistent representations of surface properties such as soil moisture and temperature, ground fluxes, and vegetation are necessary to better simulate the interactions between the land surface and atmosphere, and ultimately improve predictions of local circulations and summertime pulse convection. The NASA Short-term Prediction Research and Transition (SPORT) Center has been conducting studies to examine the impacts of high-resolution land surface initialization data generated by offline simulations of the NASA Land Informatiot~ System (LIS) on subsequent numerical forecasts using the Weather Research and Forecasting (WRF) model (Case et al. 2008, to appear in the Journal of Hydrometeorology). Case et al. presents improvements to simulated sea breezes and surface verification statistics over Florida by initializing WRF with land surface variables from an offline LIS spin-up run, conducted on the exact WRF domain and resolution. The current project extends the previous work over Florida, focusing on selected case studies of typical pulse convection over the southeastern U.S., with an emphasis on improving local short-term WRF simulations over the Mobile, AL and Miami, FL NWS county warning areas. Future efforts may involve examining the impacts of assimilating remotely-sensed soil moisture data, and/or introducing weekly greenness vegetation fraction composites (as opposed to monthly climatologies) into ol'fline NASA LIS runs. Based on positive impacts, the offline LIS runs could be transitioned into an operational mode, providing land surface initialization data to NWS forecast offices in real time.

Case, Jonathan L.↗