Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate control variates”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Exposure of Highly Reflective Far-Ultraviolet Coatings to LEO Environment for TRL Advancement for the Habitable Worlds Observatory

Over the past several years, our team at the NASA Goddard Space Flight Center (GSFC), in collaboration with the Jet Propulsion Laboratory (JPL) and the Naval Research Laboratory (NRL), have developed several new protected Al coating technologies with improved Far-Ultraviolet (FUV) reflectance performance and more robust environmental stability. These include Al protected with lithium-based fluoride (LiF) coatings such as hot deposited LiF (eLiF), XeF 2 -passivated LiF coatings (XeLiF), Li 3 AlF 6 overcoats, AlF 3 (plasma passivated), and XeF 2 passivated MgF 2 coatings (XeMgF 2 ). However, despite these developments, relatively little experimental information exists regarding the long-term stability of these coatings after simultaneous exposure to the Low Earth Orbit (LEO) environment. To evaluate their environmental durability and advance their Technology Readiness Level (TRL), representative samples of these coating technologies, Al+XeLiF, Al+eLiF, Al+Li 3 AlF 6 , Al+AlF 3 , Al+XeMgF 2 , and bare Al as reference samples, were flown onboard the International Space Station (ISS) as part of the Materials International Space Station Experiment 20 (MISSE-20). The MISSE platform is operated by Aegis Aerospace for relatively long-duration (e.g. 6 months) external exposure experiments in a space environment. The coatings were deposited on ULE and Zerodur substrates, and witness samples deposited on the same coatings runs, but on glass slides were left on Earth as control. A subset of the flight samples were fitted with MgF 2 windows to reduce direct exposure to atomic oxygen (AO) and charged particles, and still let through UV and other types of radiation on these samples. After approximately 6 months of exposure in the forward-facing (RAM) direction, the samples were returned to Earth and characterized by performing FUV/NUV/VIS/NIR reflectance measurements to quantify degradation, spectroscopic ellipsometry to evaluate coating thickness and optical constants, atomic force microscopy (AFM) to quantify surface roughness evolution during the flight in the ISS. The MgF 2 protective windows were also analyzed through FUV transmission measurements. Comparison with the witness samples left on the ground and with the flown bare Al samples (with just the naturally occurring Al 2 O 3 layer) provided a quantitative assessment of coating thickness variations, optical performance degradation, and morphological changes induced by prolonged exposure to the LEO environment.

Far Ultraviolet (FUV)↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Boundary-layer and wake measurements on a swept, circulation-control wing

Wind-tunnel measurements of boundary-layer and wake velocity profiles and surface static-pressure distributions are presented for a swept, circulation-control wing. The model is an aspect-ratio-four semispan wing mounted on the tunnel side wall at a sweep angle of 45 deg. A full-span, tangential, rearward-blowing, circulation-control slot is located ahead of the trailing edge on the upper surface. Flow surveys were obtained at mid-semispan at freestream Mach numbers of 0.425 and 0.70. Boundary-layer profiles measured on the forward portions of the wing are approximately streamwise and two-dimensional. The flow in the vicinity of the jet exit and in the near wake is highly three-dimensional. The jet flow near the slot on the Coanda surface is directed normal to the slot. Near-wake surveys show large outboard flows at the center of the wake. At Mach 0.425 and a 5-deg angle of attack, a range of jet blowing rates was found for which an abrupt transition from incipient separation to attached flow occurs in the boundary layer upstream of the slot. The variation in the lower-surface separation location with blowing rate was determined from boundary-layer measurements at Mach 0.425.

Spaid, Frank W.↗
NASA NTRS1943

A theoretical investigation of the lateral oscillations of an airplane with free rudder with special reference to the effect of friction

Charts showing the variation in dynamic stability with the rudder hinge-moment characteristics are presented. A stabilizing rudder floating tendency combined with a high degree of aerodynamic balance is shown to lead to oscillations of increasing amplitude. This dynamic instability is increased by viscous-friction in the rudder control system. The presence of solid friction in the rudder control system will cause steady oscillations of constant amplitude if the floating angle of the rudder per unit angle of sideslip is stabilizing and greater than a certain critical value that depends on other airplane parameters, such as vertical-tail area and airplane moment of inertia about the vertical axis. The amplitude of the steady oscillation is proportional to the amount of friction and is generally quite small but increases as the condition of dynamic instability is approached. An approximate method of calculating the amplitudes of the steady oscillation is explained and is illustrated by a numerical example. A more accurate step-by-step calculation of the motion is also made and it is shown that the agreement with the approximate method is good.

Greenberg, Harry↗
NASA NTRS2013

File Specification for GEOS-5 FP-IT (Forward Processing for Instrument Teams)

The GEOS-5 FP-IT Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP-IT incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP-IT uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The analysis is performed at a horizontal resolution of 0.625-degree longitude by 0.5-degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational GEOS-5 FP-IT processing in support of the EOS instrument teams. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary.

3D Var↗
NASA NTRS1986

Stresses and deformations in cross-ply composite tubes subjected to a uniform temperature change

This study investigates the effects of a uniform temperature change on the stresses and deformations of composite tubes and determines the accuracy of an approximate solution based on the principle of complementary virtual work. Interest centers on tube response away from the ends and so a planar elasticity approach is used. For the approximate solution a piecewise linear variation of stresses with the radial coordinate is assumed. The results from the approximate solution are compared with the elasticity solution. The stress predictions agree well, particularly peak interlaminar stresses. Surprisingly, the axial deformations also agree well, despite the fact that the deformations predicted by the approximate solution do not satisfy the interface displacement continuity conditions required by the elasticity solution. The study shows that the axial thermal expansion coefficient of tubes with a specific number of axial and circumferential layers depends on the stacking sequence. This is in contrast to classical lamination theory, which predicts that the expansion will be independent of the stacking arrangement. As expected, the sign and magnitude of the peak interlaminar stresses depend on stacking sequence. For tubes with a specific number of axial and circumferential layers, thermally induced interlaminar stresses can be controlled by altering stacking arrangement.

Hyer, M. W.↗
NASA NTRS1987

Experience in distributed parameter modeling of the Spacecraft Control Laboratory Experiment (SCOLE) structure

The Spacecraft Control Laboratory Experiment (SCOLE) configuration is used to compare exact and approximate solutions of the partial differential equations which define its structural dynamics. The need for a proof model for evaluating competing control laws demands that solutions be generated which not only exhibit accurate modal characteristics, but precise static deflections as well. Because precise pointing is required, the motion of the end bodies of the Shuttle-attached antenna must be known with great accuracy. Modal models are attractive because of their stable solutions but require hundreds of modes to obtain a static deflection accuracy of only one percent. Although proportional damping in bending agrees well with experimental results using the SCOLE experimental apparatus, modes which involve both torsion and bending differ significantly from proportional damping. A lumped mass model is used to generate exact static deflections, but only approximate modal characteristics. Asymptotic solutions to the distributed parameter system approximate very accurately the modal characteristics at high mode numbers. Ways are examined for refining the approximate solutions by applying a first-order variation and by employing singular perturbation techniques which are usually limited to ordinary differential equations. The most accurate solutions of the distributed parameter model of SCOLE are obtained by combining exact and asymptotic solutions.

Taylor, L. W.↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Process Model for Friction Stir Welding

Friction stir welding (FSW) is a relatively new process being applied for joining of metal alloys. The process was initially developed by The Welding Institute (TWI) in Cambridge, UK. The FSW process is being investigated at NASA/MSEC as a repair/initial weld procedure for fabrication of the super-light-weight aluminum-lithium shuttle external tank. The FSW investigations at MSFC were conducted on a horizontal mill to produce butt welds of flat plate material. The weldment plates are butted together and fixed to a backing plate on the mill bed. A pin tool is placed into the tool holder of the mill spindle and rotated at approximately 400 rpm. The pin tool is then plunged into the plates such that the center of the probe lies at, one end of the line of contact, between the plates and the shoulder of the pin tool penetrates the top surface of the weldment. The weld is produced by traversing the tool along the line of contact between the plates. A lead angle allows the leading edge of the shoulder to remain above the top surface of the plate. The work presented here is the first attempt at modeling a complex phenomenon. The mechanical aspects of conducting the weld process are easily defined and the process itself is controlled by relatively few input parameters. However, in the region of the weld, plasticizing and forging of the parent material occurs. These are difficult processes to model. The model presented here addresses only variations in the radial dimension outward from the pin tool axis. Examinations of the grain structure of the weld reveal that a considerable amount of material deformation also occurs in the direction parallel to the pin tool axis of rotation, through the material thickness. In addition, measurements of the axial load on the pin tool demonstrate that the forging affect of the pin tool shoulder is an important process phenomenon. Therefore, the model needs to be expanded to account for the deformations through the material thickness and the forging affect of the shoulder. The energy balance at the boundary of the plastic region with the environment required that energy flow away from the boundary in both radial directions. One resolution to this problem may be to introduce a time dependency into the process model, allowing the energy flow to oscillate across this boundary. Finally, experimental measurements are needed to verify the concepts used here and to aid in improving the model.

Adams, Glynn↗
NASA NTRS1994

Annealing temperature and O2 partial pressure dependence of T(sub c) in HgBa2CuO(4+delta)

Samples of HgBa2CuO(4+delta) (Hg-1201) were annealed under various conditions. After carefully controlling annealing time, annealing temperature (T(sub a)) and O2 partial pressure (P(sub 0)), we were able to find the reversible annealing conditions for Hg-1201. Under 1 atm O2 at 260 C less than or equal to T(sub a) less than or equal to 400 C, the obtained T(sub c) is nearly the same (approximately 97 K). However, it decreases quickly with T(sub a) greater than 300 C in high vacuum (P(sub 0) approximately 10(exp -8) atm), and reaches zero at T(sub a) = 400 C. On the other hand, T(sub c) decreases with the decrease of T(sub a) in high-pressure O2 (approximately 500 atm) and reaches approximately 20 K at about 240 C. In the entire annealing region, the oxygen surplus varies significantly from 0.03 to 0.4, and a wide range of T(sub c) variation (0 goes to 97 K goes to 20 K) was obtained with anion doping alone.

Xiong, Q.↗
NASA NTRS2015

A High-Resolution Merged Wind Dataset for DYNAMO: Progress and Future Plans

In order to support research on optimal data assimilation methods for the Cyclone Global Navigation Satellite System (CYGNSS), launching in 2016, work has been ongoing to produce a high‐resolution merged wind dataset for the Dynamics of the Madden Julian Oscillation (DYNAMO) field campaign, which took place during late 2011/early 2012. The winds are produced by assimilating DYNAMO observations into the Weather Research and Forecasting (WRF) three‐dimensional variational (3DVAR) system. Data sources from the DYNAMO campaign include the upper‐air sounding network, radial velocities from the radar network, vector winds from the Advanced Scatterometer (ASCAT) and Oceansat‐2 Scatterometer (OSCAT) satellite instruments, the NOAA High Resolution Doppler Lidar (HRDL), and several others. In order the prep them for 3DVAR, significant additional quality control work is being done for the currently available TOGA and SMART‐R radar datasets, including automatically dealiasing radial velocities and correcting for intermittent TOGA antenna azimuth angle errors. The assimilated winds are being made available as model output fields from WRF on two separate grids with different horizontal resolutions ‐ a 3‐km grid focusing on the main DYNAMO quadrilateral (i.e., Gan Island, the R/V Revelle, the R/V Mirai, and Diego Garcia), and a 1‐km grid focusing on the Revelle. The wind dataset is focused on three separate approximately 2‐week periods during the Madden Julian Oscillation (MJO) onsets that occurred in October, November, and December 2011. Work is ongoing to convert the 10‐m surface winds from these model fields to simulated CYGNSS observations using the CYGNSS End‐To‐End Simulator (E2ES), and these simulated satellite observations are being compared to radar observations of DYNAMO precipitation systems to document the anticipated ability of CYGNSS to provide information on the relationships between surface winds and oceanic precipitation at the mesoscale level. This research will improve our understanding of the future utility of CYGNSS for documenting key MJO processes.

Lang, Timothy J.↗
NASA NTRS2019

Seasonal- and Beta-Angle-Dependent Latitude Bias Variations in Natural Decays

Prior work has demonstrated pronounced statistical clustering of natural decays of medium-to-high-inclination orbital objects peaking approximately 30 degrees in Argument of Latitude ahead of nodal crossings. This effect is caused by the physical bulge in the Earth and the overlying atmosphere, that cyclically modifies effective altitude (and therefore density) faster than the trajectory's decay itself. While prior work has averaged seasonal and RAAN effects over all non-uniform atmosphere possibilities to support long-term characterization of the clustering of final entries in generating a pre-mission Expectation of Casualty, the current study characterizes seasonal and beta angle effects as potential influences on the near-term statistical risks of specific tactical decay scenarios, relative to the average. Such effects on the density profile along an orbit may be important considerations in any scenario where small orbital adjustments are used to optimize the timing and location of final entry trajectories. I.E., two identical spacecraft entering in different seasons and/or beta angles may have different minimum-risk scenarios for identical control capabilities and space weather conditions. Further, the early heating history of shallow trajectories is explored, examining the influence of dramatically different density profiles over the final orbit as the spacecraft either skims over or dives into the atmosphere.

Bacon, John B.↗
NASA NTRS2020

A Physical Model of Moulin Evolution on the Greenland Ice Sheet

Nearly all proglacial water discharge from the Greenland Ice Sheet is routed englacially via moulins. Identification of these moulins in high-resolution imagery is a frequent topic of study, but the processes controlling how and where moulins form remain poorly understood. Because moulins may reasonably compose approximately 10-15% of the englacial-subglacial hydrologic system, the evolution and shape of moulins can alter both the timing and variability of meltwater inputs to the bed. This evolution can impact both the form of the subglacial hydrologic system and associated response of ice motion. Here, we develop a physical model of moulin formation and evolution to constrain the role of englacial processes in shaping the form and structure of the subglacial hydrologic system. Within this model, moulin geometry is controlled by a balance of viscous and elastic deformation and is dependent on that deformation, refreezing, and the dissipation of turbulent and sensible heat energy. All of which are dependent on the characteristics of the available supraglacial meltwater and the surrounding ice. We find moulin geometry is responsive to changes in these parameters over the course of hours to days, indicating that diurnal and multi-day variations in melt can substantially alter the geometry of a moulin and, consequently, the pressure-discharge relationship at the bed of the ice sheet. Therefore, there is no single moulin shape that can appropriately represent englacial storage across the Greenland Ice Sheet.

Andrews, Lauren C.↗
NASA NTRS1983

Laboratory evaluation of an airborne ozone instrument that compensates for altitude/sensitivity effects

One problem encountered in the use of air-quality instrumentation on aircraft is the variation of instrument sensitivity with pressure as the result of altitude changes of the aircraft. Many instruments experience sensitivity changes of as much as a factor of 2 at altitudes of 6 km. Discussed are recent modifications to a chemiluminescent (ethylene) ozone detector that allow the instrument to automatically compensate for pressure/sensitivity effects. The modification provides automated mass flow rate control for both the sample and ethylene gas flows. The flow control systems maintain flow rate to within 15 percent for a 100-torr instantaneous pressure change, and flow rates are returned to the desired set points within 10 s after the pressure change. During simulated altitude changes (300 m/min from mean sea level to 3-km altitude), flow rates were controlled to within 3 percent of the set point. Laboratory data are summarized verifying the operation of the instrument for a pressure range of 760 torr (sea level) to 350 torr (approximately 20,000 ft) and an ozone concentration range from 20 to approximately 700 ppb.

Gregory, G. L.↗
NASA NTRS1994

The in-flight calibration of the Hubble Space Telescope fine guidance sensors, 2 (a success story)

The Hubble Space Telescope's fine guidance sensors (FGS's) are unique in the performance levels being attempted; spacecraft control and astrometric research with accuracies better than 3 milli-arcseconds (mas) are the ultimate goals. This paper presents a review of the in-flight calibration of the sensors, describing both the algorithms used and the results achieved to date. The work was done primarily in support of engineering operations related to spacecraft pointing and control and secondarily in support of the astrometric science calibration effort led by the Space Telescope Astrometry Team. Calibration items of principal interest are distortion, sensor magnification, and relative alignment. An initial in-flight calibration of the FGS's was performed in December 1990; this calibration has been used operationally over the past few years. Followup work demonstrated that significant, unexpected temporal variations in the calibration parameters are occurring; provided good characterization of the variation; and set the stage for a distortion calibration designed to achieve the full design accuracy for one of the FGS's. This full distortion calibration, using data acquired in January 1993, resulted in a solution having single-axis residuals with a standard deviation of 2.5 mas. Scale and alignment calibration results for all of the FGS's have been achieved commensurate with the best ground-based astrometric catalogs (root-mean-square error approximately 25 mas). A calibration monitoring program has been established to allow regular updates of the calibration parameters as needed.

Welter, Gary L.↗
NASA NTRS2001

Titan Airship Surveyor

Saturn's moon Titan is considered to be one of the prime candidates for studying prebiotic materials - the substances that precede the formation of life but have disappeared from the Earth as a result of the evolution of life. A unique combination of a dense, predominantly nitrogen, atmosphere (more than four times that of the Earth), low gravity (six times less than on the Earth) and small temperature variations makes Titan the almost ideal planet for studies with lighter-than-air aerial platforms (aerobots). Moreover, since methane clouds and photochemical haze obscure the surface, low-altitude aerial platforms are the only practical means that can provide global mapping of the Titan surface at visible and infrared wavelengths. One major challenge in Titan exploration is the extremely cold atmosphere (approx. 90 K). However, current material technology the capability to operate aerobots at these very low temperatures. A second challenge is the remoteness from the Sun (10 AU) that makes the nuclear (radioisotopic) energy the only practical source of power. A third challenge is remoteness from the Earth (approx. 10 AU, two-way light-time approx. 160 min) which imposes restrictions on data rates and makes impractical any meaningful real-time control. A small-size airship (approx. 25 cu m) can carry a payload approximately 100 kg. A Stirling engine coupled to a radioisotope heat source would be the prime choice for producing both mechanical and electrical power for sensing, control, and communications. The cold atmospheric temperature makes Stirling machines especially effective. With the radioisotope power source the airship may fly with speed approximately 5 m/s for a year or more providing an excellent platform for in situ atmosphere measurements and a high-resolution remote sensing with unlimited access on a global scale. In a station-keeping mode the airship can be used for in situ studies on the surface by winching down an instrument package. Floating above the surface allows relatively simple means for flight control. Mission requirements and possible methods of navigation, control, data acquisition, and communications are discussed. The presentation describes also the state-of-the art and current progress in aerial deployed aerobots.

Kerzhanovich, V.↗
NASA NTRS1950

Free-Flight Investigation of the Static and Dynamic Longitudinal Stability Characteristics of 1/3.7-Scale Rocket-Powered Models of the Bell MX-776A

An investigation of the static and dynamic longitudinal stability characteristics of 1/3.7 scale rocket-powered model of the Bell MX-776A has been made for a Mach number range from 0.8 to 1.6. Two models were tested with all control surfaces at 0 degree deflection and centers of gravity located 1/4 and 1/2 body diameters, respectively, ahead of the equivalent design location. Both models were stable about the trim conditions but did not trim at 0 degree angle of attack because of slight constructional asymmetries. The results indicated that the variation of lift and pitching moment was not linear with angle of attack. Both lift-curve slope and pitching-moment-curve slope were of the smallest magnitude near 0 degree angle of attack. In general, an increase in angle of attack was accompanied by a rearward movement of the aerodynamic center as the rear wing moved out of the downwash from the forward surfaces. This characteristic was more pronounced in the transonic region. The dynamic stability in the form of total damping factor varied with normal-force coefficient but was greatest for both models at a Mach number of approximately 1.25. The damping factor was greater at the lower trim normal-force coefficients except at a Mach number of 1.0. At that speed the damping factor was of about the same magnitude for both models. The drag coefficient increased with trim normal-force coefficient and was largest in the transonic region.

Michal, David H.↗
NASA NTRS2001

Constraining the Mean Crustal Thickness on Mercury

The topography of Mercury is poorly known, with only limited radar and stereo coverage available. However, radar profiles reveal topographic contrasts of several kilometers over wavelengths of approximately 1000 km. The bulk of Mercury's geologic activity took place within the first 1 Ga of the planet's history), and it is therefore likely that these topographic features derive from this period. On Earth, long wavelength topographic features are supported either convectively, or through some combination of isostasy and flexure. Photographic images show no evidence for plume-like features, nor for plate tectonics; I therefore assume that neither convective support nor Pratt isostasy are operating. The composition and structure of the crust of Mercury are almost unknown. The reflectance spectrum of the surface of Mercury is similar to that of the lunar highlands, which are predominantly plagioclase. Anderson et al. used the observed center-of-mass center-of-figure offset together with an assumption of Airy isostasy to infer a crustal thickness of 100-300 km. Based on tidal despinning arguments, the early elastic thickness (T(sub e)) of the (unfractured) lithosphere was approximately equal to or less than 100 km. Thrust faults with lengths of up to 500 km and ages of about 4 Ga B.P. are known to exist on Mercury. Assuming a semicircular slip distribution and a typical thrust fault angle of 10 degrees, the likely vertical depth to the base of these faults is about 45 km. More sophisticated modelling gives similar or slightly smaller answers. The depth to the base of faulting and the elastic layer are usually similar on Earth, and both are thought to be thermally controlled. Assuming that the characteristic temperature is about 750 K, the observed fault depth implies that the heat flux at 4 Ga B.P. is unlikely to be less than 20 mW m(exp -2) for a linear temperature gradient. For an elastic thickness of 45 km, topography at 1000 km wavelength is likely to be about 60% compensated. There are thus likely to be considerable lateral variations in crustal thickness. Additional information is contained in the original extended abstract.

Nimmo, F.↗