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.
Magnetic models, guided by results of gravity-topography admittance studies, suggest that the anomaly pattern in the central southern highlands of Mars results from large blocks of coherently magnetized crust separated by 'non-magnetic' areas. Additional information is contained in the original extended abstract.
The origin of the Martian dichotomy, which divides highlands from lowlands, is unknown. We examine a section of the dichotomy ( 50 - 90E) defined by steep scarps and normal faults. Stratigraphy and age relationships preclude the formation of the 2.5 km high boundary via erosion. The abrupt disappearance of topographic knobs similar to 300 - 500 km to the northeast is interpreted as a buried fault. Alignment of the buried fault with grabens, stratigraphy, and age determinations using crater counts indicate that the lowland bench is down faulted highlands crust. The estimated local strain (3.5%) and fault pattern are broadly consistent with gravitational relaxation of a plateau boundary. Magnetic and gravity anomalies occur on either side of the buried fault. Admittance analysis indicates isostatic compensation. Although nonunique, a model with a 10 km thick intracrustal block under the lowland bench, a 20 km thick block under the plains, and an excess density of 200 kg/m(3) provides a good fit to the isostatic anomaly. A good fit to a profile of the magnetic field perpendicular to the dichotomy is produced using uniformly polarized intracrustal blocks 10 - 20 km thick, an intensity of 6 Am/m, a field inclination of -30 degrees, and gaps aligned with the isostatic anomalies. One interpretation is that high-density intrusions demagnetized the crust after dynamo cessation and that low-lying magnetized areas could be down faulted highlands crust. Another model (inclination of 30 degrees) has magnetized crust beneath the isostatic anomalies, separated by gaps. The gaps could result from hydrothermal alteration of the crust along fault zones.
Unlike Earth, on Venus long wavelength geoid anomalies correlate well with topography. Venus's admittance curve between harmonic degrees 3 and 18 is inconsistent with Airy isostasy but is consistent with dynamic support from convection being the dominant mechanism of compensation on Venus. We model dynamic compensation on Venus using simple flow models which assume a spherically symmetric Newtonian mantle viscosity profile. Preliminary models parameterize the viscosity variation with depth as a 2 layer model with a boundary at 720 km depth. A model in which viscosity in the lower mantle is a factor of 10 lower than in the upper mantle can explain Venus's observed admittance curve for degrees 3 through 18. Dynamic models which include a chemical boundary between the upper and lower mantle do not successfully explain the observed admittance curve, indicating that Venus does not have a chemically layered mantle.
Electric utilities are facing the need for better monitoring, analysis, and control of their distribution systems. An accurate mathematical model is a key to both the development of cutting-edge, scalable model-based algorithms and the assessment of emerging technologies such as distributed energy resources (DER) for grid planning and operation. However, the constantly evolving nature of power distribution systems poses challenges to maintaining accurate models. In this paper, we propose a novel load flow based approach to validate power distribution models. Networked equipment models described according to the Common Information Model (CIM) standard and a measurement model are used to formulate the distribution load flow problem. First, a system admittance matrix (Ybus) is derived from device-level CIM parameters. Next, the operational parameters (dynamic Ybus and nodal injections) are extracted from the measurement model using sensor configuration and equipment state. An iterative power flow method is then used to compute nodal voltages and branch flows that are compared against the measurement data to find any inconsistencies in the networked equipment model. This approach is implemented within GridAPPS-D, an open-source standards-based platform for advanced distribution management system (ADMS) application development, and demonstrated on the IEEE 13-bus, 123-bus, and 8500-node test feeders.
Measurements of reactive admittances that describe the capabilities of the combustion processes associated with coaxial gaseous fuel injectors to amplify combustor oscillations are presented. These admittances are needed in (1) stability analyses of rocket motors and gas turbine combustors; (2) for the evaluation of 'driving' capabilities of injectors; and (3) for checking the applicability of a theoretical model. The modified standing-wave technique was used to determine the admittances of the combustion processes in coaxial injectors utilizing air-acetylene mixtures. The measurements indicated that (1) coaxial injectors can sustain combustion instabilities over wide frequency ranges; (2) the maximum driving capability of an injector decreases with an increased equivalence ratio; (3) the frequency at which maximum driving is observed decreases with the increased equivalence ratio; and (4) the characteristic combustion time of an injector design decreases with increased frequency.
The time evolving strain field contains a wealth of information that can be used to interpret subsurface behavior. For example, injecting or removing fluids from reservoirs or aquifers causes deformation that can be used as a diagnostic signal in some cases, while it can interfere with geodetic interpretations in other cases. We've previously demonstrated the feasibility of measuring the strain tensor at a depth of 30m caused by injection into a reservoir at 530m. The observed strain signals were interpreted using four independent analytic and numerical methods that resulted in estimates of the poroelastic properties and geometry of the reservoir that was consistent with data from well logs. However, studies like these are only possible if these deformations can be reliably measured. Years of lab and field work has culminated in the development of a novel borehole strainmeter capable of resolving multiple components of strain using embedded optical fibers configured as Michelson interferometers. It features four horizontal gauges separated by 45° to resolve the horizontal strain tensor as well as a vertical strain gauge and a sixth null component for state-of-health monitoring. The downhole sensing package also includes an open pipe through its center for grout circulation during single-trip deployments and a fully welded stainless steel exterior for robustness and longevity. These instruments have a resolution of 2x10-13 strain that can easily measure the solid earth tides. Preliminary data are available from four strainmeters in shale at our Oklahoma site and four in compacted sand and gravel in Utah. These are deployed from 40-60m, except one of the strainmeters in Oklahoma is deployed at 500m. The data include strains from the initial grout curing, comparisons to predicted earth tide models and in-situ calibration results, barometric pressure admittances and spectral analyses as well as signals from underground injections and surface waves from teleseismic events. Preliminary analyses indicate behavior consistent with other strainmeter deployments, and comparison to data from a Gladwin strainmeter at the Oklahoma site validate the performance of the new design. Analyses from a suite of six well tests at the Oklahoma site show for the first time how the strain tensor field varies with location during well testing.
Gravity information released by the Soviet Union allows the quantitative assessment of how the geopotential research mission (GRM) mission might effect the ability to use global gravity data for continental tectonic interpretation. The information is of an isostatic response spectra for eight individual tectonic units in the USSR. The regions examined include the Caroathians, Caucasus, Urals, Pamirs, Tien-Shan, Altal, Chersky Ridge, and East Siberian Platform. The 1 deg x 1 deg gravity data are used to calculate the admittances are used in two different sorts of tectonic studies of mountain belts in the USSR: (1) interpretation of isostatic responses in terms of plate models of compensation for mountainous terrain. Using geologic information concerning time of the orogeny, lithospheric plates involved, and polarity of subduction in collision zones, they convert the best-fitting flexural rigidity to an elastic plate thickness for the lithospheric plate inferred to underlie the mountains; the isostatic admittance functions is an attempt to directly model gravity and topography data for a few select regions in the Soviet Union. By knowing the value of the expected correlation between topography and gravity from the admittances, the Artemjev's map in mountainous areas can be calibrated, and the maps are converted back to Bouguer gravity. This procedure is applied to the Caucasus and southern Urals.
Finite element theory is used to calculate the acoustic field of a propeller in a soft walled circular wind tunnel and to compare the radiation patterns to the same propeller in free space. Parametric solutions are present for a "Gutin" propeller for a variety of flow Mach numbers, admittance values at the wall, microphone position locations, and propeller to duct radius ratios. Wind tunnel boundary layer is not included in this analysis. For wall admittance nearly equal to the characteristic value of free space, the free field and ducted propeller models agree in pressure level and directionality. In addition, the need for experimentally mapping the acoustic field is discussed.
Finite element theory is used to calculate the acoustic field of a propeller in a soft walled circular wind tunnel and to compare the radiation patterns to the same propeller in free space. Parametric solutions are present for a 'Gutin' propeller for a variety of flow Mach numbers, admittance values at the wall, microphone position locations, and propeller to duct radius ratios. Wind tunnel boundary layer is not included in this analysis. For wall admittance nearly equal to the characteristic value of free space, the free field and ducted propeller models agree in pressure level and directionality. In addition, the need for experimentally mapping the acoustic field is discussed.
Magellan gravity data confirm that several of the large, tectonically deformed, plateau-like highlands on venus are shallowly compensated, most likely by crustal thickness variations. Apparent depths of isostatic compensation, computed in the spatial domain, range from 30 to 50 km for Alpha, Tellus, Ovda, and Thetis Regiones. Using a two-layer model for isostatic compensation, Alpha, Tellus, and Ovda are best represented as nearly completely compensated in crust that is regionally 20-40 km thick around these highlands, with little contribution from deeper mantle sources. In contrast to these three areas, a stronger regional gravity high associated with Thetis requires a significant upper mantle component to compensation. This is evident in the spectral admittance as a pronounced deep, long-wavelength anomaly. In the two-layer isostatic model, a broad, deeply compensated upland underlies a shallowly compensated central block of Thetis. If this deep component is interpreted as a thermal anomaly, the loci of maximum upwelling agree well with sites of recent extension. The plateau highlands are thus physiographically and isostatically equivalent to terrestrial continents, though probably not compositionally. They also share the record of a long tectonic history. The large regional gravity anomaly of Thetis indicates that active mantle proceses continue even beneath some areas (tessera) thought to be a relic of a former geological regime. The excellent agreement of modeled crustal thicknesses around Alpha, Tellus, and Ovda Regiones suggests that 20-40 km is a representative global value for the plains. Such a crust is thicker than previously estimated and about twice as thick as the expected thickness of crust produced at venusian spreading centers
A new analysis of the Doppler tracking data from the Lunar Prospector mission in 1999 revealed a number of previously-unseen gravity anomalies at spatial scales as small as 27 km over the nearside. The tracking data at low altitudes (50 km or below) were better analyzed to resolve the nearside features without dampening from a power law constraint, by partitioning the gravity parameters concentrated on either the nearside or farside. The resulting model presents gravity anomalies correlated with topography with a correlation coefficient of 0.7 or higher from degree 50 to 150, the widest bandwidth yet. The gravity-topography admittance of approx. 70 mGal/km is found from numerous craters of which diameters are 60 km or less. In addition, the new model produces orbits that fit to independent radio tracking data from the Lunar Reconnaissance Orbiter and Kaguya (SELENE) better than previous gravity models. This high-resolution model can be of immediate use to geophysical analysis of small craters. Our technique could be applied to an upcoming mission, the Gravity Recovery And Interior Laboratory and useful to extract short wavelength signals from the MESSENGER Doppler data.
We present a method with which we determined the local lunar gravity field model over the South Pole-Aitken (SPA) basin on the farside of the Moon by estimating adjustments to a global lunar gravity field model using SELENE tracking data. Our adjustments are expressed in localized functions concentrated over the SPA region in a spherical cap with a radius of 45deg centered at (191.1 deg E, 53.2 deg S), and the resolution is equivalent to a 150th degree and order spherical harmonics expansion. The new solution over SPA was used in several applications of geophysical analysis. It shows an increased correlation with high-resolution lunar topography in the frequency band l = 40-70, and admittance values are slightly different and more leveled when compared to other, global gravity field models using the same data. The adjustments expressed in free-air anomalies and differences in Bouguer anomalies between the local solution and the a priori global solution correlate with topographic surface features. The Moho structure beneath the SPA basin is slightly modified in our solution, most notably at the southern rim of the Apollo basin and around the Zeeman crater
[i] From gravity and topography data collected by the Mars Global Surveyor spacecraft we calculate gravity/topography admittances and correlations in the spectral domain and compare them to those predicted from models of lithospheric flexure. On the basis of these comparisons we estimate the thickness of the Martian elastic lithosphere (T(sub e)) required to support the observed topographic load since the time of loading. We convert T(sub e) to estimates of heat flux and thermal gradient in the lithosphere through a consideration of the response of an elastic/plastic shell. In regions of high topography on Mars (e.g., the Tharsis rise and associated shield volcanoes), the mass-sheet (small-amplitude) approximation for the calculation of gravity from topography is inadequate. A correction that accounts for finite-amplitude topography tends to increase the amplitude of the predicted gravity signal at spacecraft altitudes. Proper implementation of this correction requires the use of radii from the center of mass (collectively known as the planetary shape ) in lieu of topography referenced to a gravitational equipotential. Anomalously dense surface layers or buried excess masses are not required to explain the observed admittances for the Tharsis Montes or Olympus Mons volcanoes when this correction is applied. Derived T, values generally decrease with increasing age of the lithospheric load, in a manner consistent with a rapid decline of mantle heat flux during the Noachian and more modest rates of decline during subsequent epochs.
The present paper focuses on resonance and energetics of the daily tides, especially in the southern ocean, the distribution of gravitational power input of daily and half-daily tides, and comparison with other estimates of global dissipation rates. The present global tidal maps, derived from Geosat altimetry, compare favorably with ground truth data at about the same rms level as the models of Schwiderski (1983), and are slightly better in lunar than in solar tides. Diurnal admittances clearly show Kelvin wave structure in the southern ocean and confirm the resonant mode of Platzman (1984) at 28.5 + or - 0.1 hr with an apparent Q of about 4. Driving energy is found to enter dominantly in the North Pacific for the daily tides and is strongly peaked in the tropical oceans for the half-daily tides. Global rates of working on all major tide constituents except S2 agree well with independent results from analyses of gravity through satellite tracking. Comparison at S2 is improved by allowing for the air tide in gravitational results but suggests deficiencies in all solar tide models.
A procedure which can accurately characterize injector designs for large thrust (0.5 to 1.5 million pounds), high pressure (500 to 3000 psia) LOX/hydrocarbon engines is currently under development. In this procedure, a rectangular cross-sectional combustion chamber is to be used to simulate the lower traverse frequency modes of the large scale chamber. The chamber will be sized so that the first width mode of the rectangular chamber corresponds to the first tangential mode of the full-scale chamber. Test data to be obtained from the rectangular chamber will be used to assess the full scale engine stability. This requires the development of combustion stability models for rectangular chambers. As part of the combustion stability model development, a computer code, NOAD based on existing theory was developed to calculate the nozzle admittances for both rectangular and axisymmetric nozzles. This code is detailed.
The results presented represent the first successful attempt at the measurement of the driving capabilities of coaxial gaseous propellant rocket injectors. The required data have been obtained by employing the modified impedance tube technique with compressed air as the oxidizer and acetylene gas as the fuel. The data describe the frequency dependence of the injector admittances, from which the frequency dependence of the injection response factors can be calculated. The measured injector admittances have been compared with the predictions of the Feiler and Heidmann (1967) analytical model assuming different values for the characteristic combustion time. The values of combustion time which result in a best fit between the measured and predicted data are indicated for different equivalence ratios. It is shown that for the coaxial injector system investigated in this study the characteristic combustion times vary between .7 and 1.2 msec for equivalence ratios in the range of .57 to 1.31. The experimental data clearly show that the tested injector system could indeed drive combustion instabilities over a frequency range that is in qualitative agreement with the predictions of the Feiler and Heidmann model.
This study reports the electronic properties of the indium tin oxide (ITO)/β-Ga2O3 Schottky diode through the temperature (50–320 K)-dependent current density–voltage, capacitance–voltage, and admittance spectroscopy experiments. The room-temperature turn-on voltage of the ITO/β-Ga2O3 diode was observed to be 0.83 V, and it exhibited a slow increase with decreasing temperature. The ITO/β-Ga2O3 diode's ideality factor (n) varied from 3.92 to 1.05, and the zero-bias Schottky barrier height (φbo) varied from 0.31 to 1.28 eV. The temperature-dependent n and φbo indicate the spatial inhomogeneities of the potential barrier at the Ga2O3-ITO interface. We exploit the Gaussian distribution model to explain the Schottky barrier inhomogeneities, which could be associated with the defects observed by the admittance spectroscopy method.
The geoid and topographic fields of the central Pacific were delineated and shown to correlate closely at intermediate wavelengths (500 to 2500 km). The associated admittance shows that anomalies having wavelengths less than about 1000 km are probably supported by the elastic strength of the lithosphere. Larger wavelength anomalies are due to dynamic effects in the sublithosphere. Direct modeling of small scale convection in the asthenosphere shows that the amplitudes of observed geoid and topographic anomalies can be independently matched, but that the observed admittance cannot. Only by imposing an initial regional variation in the thermal regime is it possible to match the admittance. It is proposed that this variation may be due to differences in the onset time of convection beneath the lithosphere of different ages. That is, convection beneath thickening lithosphere is strongly dependent on the rate of thickening (V) relative to the rise time for convection. The critical Rayleigh number contains the length scale K/V, where K is thermal diffusivity. Young, fast growing lithosphere stabilizes the underlying asthenosphere unless it has an unusually low viscosity. Lithosphere of different age, separated by fracture zones, will go unstable at different times, producing regional horizontal temperature gradient that may strongly influence convection. Laboratory and numerical experiments are proposed to study this form of convection and its influence on the geoid.