SEARCH · Engineering Papers

Results for “Planetary Boundary Height”

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.

SourceResearch connectionFrom yearThrough year

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

196 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗

NASA NTRS1992

Influence of aerodynamic roughness length on aeolian processes: Earth, Mars, Venus

The aerodynamic roughness length (z sub 0) is the height at which a wind profile assumes a zero velocity. The lower part of the atmospheric boundary layer will be impeded by friction with the surface. An increase in surface roughness will also increase the shear stress required to initiate particle entrainment by the wind. Bagnold (1941) estimated z sub 0 as being 1/30 of the mean particle size. In Nature, surface roughness is composed of nonerodible elements as well as sand-size erodible particles. To assess z sub 0 values as a function of terrain, field experiments were conducted to obtain wind profiles monitored over natural surfaces at 15 sites in the Mojave Desert, Death Valley, and Nye County, Nevada. These sites span a variety of arid-land terrains, including smooth playas, alluvial fans, and lava flows; z sub 0 values ranged from 0.0001 cm to 1 cm. These values were incorporated in a threshold model and a flux model to assess transport efficiency over such terrains in three planetary environments (Venus, Earth, and Mars), and for particle sizes ranging from 60-500 micron. Threshold and flux are a function of planetary environment, particle density and size (Dp), and z sub 0, and the shear velocity of 1.2 x U sub *t (for Dp = 250 micron and z sub 0 = 0.84). Results show that flux on Mars is approximately 14 g/(cm x s), on Earth it is approximately 3 g/(cm x s), and on Venus 0.5 g/(cm x s). Under all planetary environments, the results also show a dramatic decrease in the flux for particles greater than 200 microns when z sub 0 increases above 0.0085 cm (corresponding to sites consisting of alluvium). When z sub 0 approaches 0.03 cm (corresponding to a mantled pahoehoe lava), the flux diminishes.

Blumberg, Dan G.↗

DOE OSTI2021

Evaluating the performance of WRF urban schemes and PBL schemes over Dallas-Fort Worth during a dry summer and a wet summer

This study evaluated the Weather Research and Forecasting (WRF) model sensitivity to different planetary boundary layer (PBL) schemes (the YSU and MYJ schemes) and urban schemes including the bulk scheme (BULK), single-layer urban canopy model (UCM), multi-layer building environment parameterization (BEP) model, and multi-layer building energy model (BEM). Daily reinitialization simulations were conducted over Dallas-Fort Worth during a dry summer month (July 2011) and a wet summer month (July 2015) with weaker (stronger) daytime (nocturnal) UHI in 2011 than 2015. All urban schemes overestimated the urban daytime 2m temperature in both summers, but BEP and BEM still reproduced the daytime urban cool island in dry summer. All urban schemes reproduced the nocturnal urban heat island, with BEP producing the weakest one due to its unrealistic urban cooling. BULK and UCM overestimated the urban canopy wind speed, while BEP and BEM underestimated it. The urban schemes showed prominent impact on daytime PBL profiles. UCM+MYJ showed a superior performance than other configurations. The relatively large (small) aspect ratio between building height and road width in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy temperature. The relatively low (high) building height in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy wind speed. Improving urban schemes and providing realistic urban parameters were critical for improving urban canopy simulation.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1989

A PBL-radiation model for application to regional numerical weather prediction

Often in the short-range limited-area numerical weather prediction (NWP) of extratropical weather systems the effects of planetary boundary layer (PBL) processes are considered secondarily important. However, it may not be the case for the regional NWP of mesoscale convective systems over the arid and semi-arid highlands of the southwestern and south-central United States in late spring and summer. Over these dry regions, the PBL can grow quite high up into the lower middle troposphere (600 mb) due to very effective solar heating and hence a vigorous air-land thermal interaction can occur. The interaction representing a major heat source for regional dynamical systems can not be ignored. A one-dimensional PBL-radiation model was developed. The model PBL consists of a constant-flux surface layer superposed with a well-mixed (Ekman) layer. The vertical eddy mixing coefficients for heat and momentum in the surface layer are determined according to the surface similarity theory, while their vertical profiles in the Ekman layer are specified with a cubic polynomial. Prognostic equations are used for predicting the height of the nonneutral PBL. The atmospheric radiation is parameterized to define the surface heat source/sink for the growth and decay of the PBL. A series of real-data numerical experiments has been carried out to obtain a physical understanding how the model performs under various atmospheric and surface conditions. This one-dimensional model will eventually be incorporated into a mesoscale prediction system. The ultimate goal of this research is to improve the NWP of mesoscale convective storms over land.

Chang, Chia-Bo↗

DOE OSTI2026

Observed Temperature and Moisture Advection Characteristics in the Southern Great Plains

Until recently, advection profiles in the planetary boundary layer have been difficult to quantify from observations. As a result, not much is known about its basic characteristics, like magnitude and diurnal cycles, or how these differ as a function of large-scale environments. Here, to provide insight into advective characteristics, a Green’s Theorem–based method for calculating profiles of advection in the lowest 3 km from an array of ground-based thermodynamic and kinematic profiling instruments has been applied to 2 years of observations from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observing site. Since advection largely depends on synoptic-scale forcing, a self-organizing map (SOM) was used to classify days based on mean sea level pressure analyses; from that, the diurnal evolution of vertical profiles of advection for different synoptic environments was quantified. The overall mean magnitude of the potential temperature advection is approximately ±3 K h −1 , while the moisture advection is approximately ±1.5 g kg −1 h −1 , with substantial variability in the sign and magnitude of the advective tendencies at different heights and throughout the diurnal cycle. Advection magnitude is strongly connected to the strength of synoptic-scale forcing, which varies depending on the time of year. Variability in advection is larger for potential temperature than it is for moisture, and strongly forced environments, like midlatitude cyclones, feature greater variability in the advection magnitude than weakly forced or quiescent environments.

advection↗

DOE OSTI2022

High-Resolution Wind Resource Data Set of the Greater Puerto Rico Region

In February 2022, the U.S. Department of Energy and six national laboratories launched the Puerto Rico Grid Resilience and Transitions to 100% Renewable Energy Study (PR100). PR100 aims to provide a comprehensive analysis of possible pathways for Puerto Rico's energy future, with a goal of 100% renewable energy by 2050. As a part of the renewable energy potential assessment in this project, we developed 20 years (2001-2020) of data using a numerical weather prediction (NWP) model for onshore and offshore wind resource assessment for the Puerto Rico. The research steps in developing the long-term wind resource data sets based on the NWP model were: 1. Model wind resource based on the Weather Research and Forecasting (WRF) model. 2. Develop WRF model configurations for Puerto Rico. 3. Test WRF with 11 different physics parameterizations for planetary boundary layer (PBL). 4. Assess WRF output from the different PBL schemes against observations. 5. Select a final model configuration which can produce the modeled wind speed with sufficient accuracy. 6. Produce 20 years wind resource data sets for Puerto Rico region. In the first stage of our framework for developing wind resource data, we developed a WRF model configuration using two nested domains (9 km and 3 km) to cover Puerto Rico and U.S. Virgin Islands and downscale the ERA5 reanalysis data (0.25 degrees x 0.25 degrees; hourly interval) to a 3-km domain. For the second stage, we implemented one-year simulations focused on using 11 different PBL physics parameterizations to find a combination of WRF physics parameterizations that could provide accurately modeled wind speed for Puerto Rico. We also analyzed the sensitivity of the modeled wind speed to PBL schemes for onshore and offshore locations. The WRF output resulting from the 11 WRF experiments using different PBL parameterizations were evaluated against observations obtained from the National Data Buoy Center (NDBC) as well as at hub height for a location for which measurements were available. A final model setup selected through the validation with observational data was used to produce 20 years of data with 3-km spatial and 5-minute temporal resolution. The WRF model output was post-processed to include wind profiles and basic atmospheric variables in a format that can be easily used for downstream modeling. The 20 years of wind resource data will be made available through NREL and support the estimation of wind energy development costs for the PR100 study.

17 WIND ENERGY↗

NASA NTRS2015

Intercomparison of Martian Lower Atmosphere Simulated Using Different Planetary Boundary Layer Parameterization Schemes

We use the mesoscale modeling capability of Mars Weather Research and Forecasting (MarsWRF) model to study the sensitivity of the simulated Martian lower atmosphere to differences in the parameterization of the planetary boundary layer (PBL). Characterization of the Martian atmosphere and realistic representation of processes such as mixing of tracers like dust depend on how well the model reproduces the evolution of the PBL structure. MarsWRF is based on the NCAR WRF model and it retains some of the PBL schemes available in the earth version. Published studies have examined the performance of different PBL schemes in NCAR WRF with the help of observations. Currently such assessments are not feasible for Martian atmospheric models due to lack of observations. It is of interest though to study the sensitivity of the model to PBL parameterization. Typically, for standard Martian atmospheric simulations, we have used the Medium Range Forecast (MRF) PBL scheme, which considers a correction term to the vertical gradients to incorporate nonlocal effects. For this study, we have also used two other parameterizations, a non-local closure scheme called Yonsei University (YSU) PBL scheme and a turbulent kinetic energy closure scheme called Mellor- Yamada-Janjic (MYJ) PBL scheme. We will present intercomparisons of the near surface temperature profiles, boundary layer heights, and wind obtained from the different simulations. We plan to use available temperature observations from Mini TES instrument onboard the rovers Spirit and Opportunity in evaluating the model results.

Natarajan, Murali↗

NASA NTRSDate not supplied

Differential absorption lidar (DIAL) for water vapor and aerosol profiling from airborne and space-based platforms

Many current challenges in weather and climate prediction involve cloud and aerosol processes, which are intimately linked to processes that initiate or take place within the planetary boundary layer (PBL). The current space-based observing system provides a wealth of information on the spatial distribution of key atmospheric variables such as profiles of water vapor (WV) and temperature, but lacks the sensitivity, accuracy, and vertical resolution in the PBL to overcome the observational gaps that inhibit the parametrization and advancement of current climate and weather models. To overcome some of these observational gaps and improve understanding and monitoring of these key processes, mature remote sensing techniques like Differential Absorption Lidar (DIAL) can be employed from various platforms to capture processes on a wide range of spatial scales. This presentation will provide an overview of NASA’s water vapor DIAL capabilities, present the strengths, weaknesses, and synergies of the DIAL technique relative to traditional methods, and provide representative data examples from the HALO airborne water vapor DIAL which was designed as a technology testbed for future space-based DIAL missions as well as a workhorse instrument in support of suborbital process studies. HALO is an airborne multifunction lidar system, simultaneously profiling WV with the DIAL technique and aerosol and cloud optical properties and PBL heights via the high spectral resolution lidar (HSRL) technique. HALO has a compact design that allows it to fly on most research aircraft, enabling nadir profiling from the lower stratosphere down through the PBL, which is critical for studying processes such as convection and cloud formation that can link the PBL and free troposphere. We will end with an outlook on the benefits and challenges of transitioning the mature DIAL technique to space to compliment other active and passive sounders and enable new PBL observations for the future Earth Observing System. The Atmospheric Boundary layer Lidar pathfindEr (ABLE), a program aimed at advancing technologies to enable the first demonstration of DIAL from space with cross-cutting observables of PBL heights, aerosol and cloud profiling, and methane columns from space will also be presented.

Brian J. Carroll↗

NASA NTRS2020

Ancient Martian aeolian sand dune deposits recorded in the stratigraphy of Valles Marineris and implications for past climates

Aeolian sediment transport, deposition, and erosion have been ongoing throughout Mars’ history. This record of widespread aeolian processes is preserved in landforms and geologic units that retain important clues about past environmental conditions including wind patterns. In this study we describe landforms within Melas Chasma, Valles Marineris, that occur in distinct groups with linear to crescentic shapes, arranged with a characteristic wavelength; some possess slope profiles analogous to modern sand dunes, yet show evidence for lithification. Based on the features’ dimensions, asymmetry, and spatial patterns relative to modern equivalents, we interpret these landforms to be two classes of aeolian bedforms: decameter-scale mega-ripples and sand dunes. The presence of superposed erosional features and depositional units indicates these landforms were cemented and likely ancient. Melas paleo-dunes are found atop Hesperian-aged layered deposits, but we estimate them to be younger, likely lithified in the Amazonian period. Although a range of degradation was observed, some paleo-dunes are >10 m tall and maintain steep lee sides (>25°), an uncommon scenario for terrestrial examples as other geologic processes lead to dune obliteration. The preserved paleo-bedform geometries are largely consistent with those of modern aeolian indicators, suggesting no major shifts in wind regime or contributing boundary conditions. Finally, we propose their appearance and context require sequential periods of dune migration, stabilization following catastrophic burial, cementation, differential erosion, exposure, and burial. The presence of wholly preserved duneforms appears to be more common on Mars compared to the Earth and may signal something important about Martian landscape evolution. Plain Language Summary Wind-driven sand dunes are common on modern Mars and the presence of certain sedimentary rock layers indicates these landforms occurred there in the past. Here, we explore features in the canyons of Valles Marineris which show characteristics frequently attributed to dunes, yet evidence for their lithification and burial is clear. Their horizontal extent, height, shape, slopes, and collective spatial patterns led to our interpretation of them as lithified dune fields. The heavily eroded appearance along with superposed elements such as craters and boulders led to the conclusion that these are relatively ancient landforms. Despite this erosion This level of preservation is rare for terrestrial sand dunes due to erosion and tectonics, thus providing an opportunity to reconstruct the various factors contributing to their history. Understanding the circumstances that led to the preservation of these ancient duneforms will yield crucial information regarding planetary sedimentary processes and the geologic history of the region.

Matthew Chojnacki↗

NASA NTRS1994

Numerical modeling studies of wake vortex transport and evolution within the planetary boundary layer

The proposed research involves four tasks. The first of these is to simulate accurately the turbulent processes in the atmospheric boundary layer. TASS was originally developed to study meso-gamma scale phenomena, such as tornadic storms, microbursts and windshear effects in terminal areas. Simulation of wake vortex evolution, however, will rely on appropriate representation of the physical processes in the surface layer and mixed layer. This involves two parts. First, a specified heat flux boundary condition must be implemented at the surface. Using this boundary condition, simulation results will be compared to experimental data and to other model results for validation. At this point, any necessary changes to the model will be implemented. Next, a surface energy budget parameterization will be added to the model. This will enable calculation of the surface fluxes by accounting for the radiative heat transfer to and from the ground and heat loss to the soil rather than simple specification of the fluxes. The second task involves running TASS with prescribed wake vortices in the initial condition. The vortex models will be supplied by NASA Langley Research Center. Sensitivity tests will be performed on different meteorological environments in the atmospheric boundary layer, which include stable, neutral, and unstable stratifications, calm and severe wind conditions, and dry and wet conditions. Vortex strength may be varied as well. Relevant non-dimensional parameters will include the following: Richardson number or Froude number, Bowen ratio, and height to length scale ratios. The model output will be analyzed and visualized to better understand the transport, decay, and growth rates of the wake vortices. The third task involves running simulations using observed data. MIT Lincoln Labs is currently planning field experiments at the Memphis airport to measure both meteorological conditions and wake vortex characteristics. Once this data becomes available, it can be used to validate the model for vortex behavior under different atmospheric conditions. The fourth task will be to simulate the wake in a more realistic environment covering a wider area. This will involve grid nesting, since high resolution will be required in the wake region but a larger total domain will be used. During the first allocation year, most of the first task will be accomplished.

Lin, Yuh-Lang↗

DOE OSTI2024

Aerosol Vertical Turbulent Mass Flux Retrievals Through Novel Remote Sensing Algorithm

Abstract Integrated measurements of aerosol, radiation, cloud, and turbulent transport in the planetary boundary layer (PBL) are essential for understanding and modeling climate and air quality. Here, we developed a new technique for the identification of convective turbulent regions and deriving the vertical distribution of aerosol turbulent mass fluxes within PBL. The algorithm uses retrievals from coherent Doppler lidars and a high spectral resolution lidar. The technique was applied to study particle mass fluxes over 2 months (November–December 2020) during the campaign conducted at the DOE Atmospheric Radiation Measurement Southern Great Plains (SGP) site in Lamont, Oklahoma. The algorithm developed here is capable of continuously deriving vertically resolved (curtains) aerosol mass fluxes. Our data analysis shows that at the site, the 30‐min averaged fluxes at 135 m above the surface were mainly positive (upward) at ∼1 μg m −2 s −1 , suggesting that the surface is the primary source of the particle mass supplied to the boundary layer at the SGP site. Analyses of the individual case studies have revealed that not all the derived fluxes can be linked to surface emissions. Both positive and negative values in a range of ±5 μg m −2 s −1 can be caused by convective thermals interacting between the residual layer and the mixed layer and by rotation of the horizontal wind with the height. Large erroneous negative fluxes can also be caused by drizzling/precipitating clouds. We anticipate that the application of the current technique will lead to a more realistic representation of aerosol mass budgets and bidirectional mixing rates.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

Sensing Water Vapon via Spacecraft Radio Occultation Observations

The radio occultation technique has been used to characterize planetary atmospheres since the 1960's spanning atmospheric pressures from 16 microbars to several bars. In 1988, the use of GPS signals to make occultation observations of Earth's atmosphere was realized by Tom Yunck and Gunnar Lindal at JPL. In the GPS to low-Earth-orbiter limb- viewing occultation geometry, Fresnel diffraction yield a unique combination of high vertical resolution of 100 m to 1 km at long wavelengths (approx. 20 cm) insensitive to particulate scattering which allows routine limb sounding from the lower mesosphere through the troposphere. A single orbiting GPS/GLONASS receiver can observe - 1000 to 1400 daily occultations providing as many daily, high vertical resolution soundings as the present global radiosonde network, but with far more evenly distributed, global coverage. The occultations yield profiles of refractivity as a function of height. In the cold, dry conditions of the upper troposphere and above (T less than 240 K), profiles of density, pressure (geopotential), and temperature can be derived. Given additional temperature information, water vapor can be derived in the midddle and lower troposphere with a unique combination of vertical resolution, global distribution and insensitivity to clouds and precipitation to an accuracy of approx. 0.2 g/kg. At low latitudes, moisture profiles will be accurate to 1-5% within the convective boundary layer and better than 20% below 6 to 7 km. Accuracies of climatological averages should be approx. 0. 1 g/kg limited by the biases in the temperature estimates. To use refractivity to constrain water vapor, knowledge of temperature is required. The simplest approach is to use the temperature field from an analysis such as the 6 hour ECMWF global analysis interpolated to the locations of each occultation. A better approach is to combine the temperature and moisture fields from such an analysis with the occultation refractivity in a weighting scheme based on the errors in each data field. A ID variational combinational approach has been developed at the UKMO. We win present results from both approaches from GPS/MET data taken in June and July 1995 and compare them with the ECMWF global 6 hour moisture analyses which are derived largely from TOVS and radiosonde data. Overall, the atmosphere below the 500 mb level appears somewhat drier in general than the ECNIWF humidity field. A 2-D (latitude vs. height) climatological snapshot derived from a 2-week span of GPS/MET data will be compared to the humidity climatology of Peixoto and Oort derived from radiosonde data from 1963-1973. Differences between the GPS results and Peixoto and Oort may be the signature of a climate trend over the past 30 years.

Kursinski, E. Robert↗

NASA NTRS2018

Development of a Mars Lidar (MARLI) for Measuring Wind and Aerosol Profiles from Orbit

Our understanding of the Mars atmosphere and the coupled atmospheric processes that drive its seasonal cycles is limited by a lack of observation data, particularly measurements that capture diurnal and seasonal variations on a global scale. As outlined in the 2011 Planetary Science Decadal Survey and the recent Mars Exploration Program Analysis Group(MEPAG) Goals Document, near-polar-orbital measurements of height-resolved aerosol backscatter and wind profiles area high-priority for the scientific community and would be valuable science products as part of a next-generation orbital science package. To address these needs, we have designed and tested a breadboard version of a direct detection atmospheric wind lidar for Mars orbit. It uses a single-frequency, seeded Nd:YAG laser ring oscillator operating at 1064nm (4 kHz repetition rate), with a 30-ns pulse duration amplified to 4 mJ pulse energy. The receiver uses a Fabry-Perotetalon as part of a dual-edge optical discrimination technique to isolate the Doppler-induced frequency shift of the back scattered photons. To detect weak aerosol backscatter profiles, the instrument uses a 4x4 photon-counting HgCdTeAPD detector with a 7 MHz bandwidth and < 0.4 fW/Hz(exp 1/2) noise equivalent power. With the MARLI lidar breadboard instrument, we were able to measure Doppler shifts continuously between 1 and 30 m/s by using a rotating chopper wheel to impart a Doppler shift to incident laser pulses. We then coupled the transmitter and receiver systems to a laser ranging telescope at the Goddard Geophysical and Astronomical Observatory (GGAO) to measure backscatter and Doppler wind profiles in the atmosphere from the ground. We measured a 5.3 ± 0.8 m/s wind speed from clouds in the planetary boundary layer at a range of 4 to 6 km. This measurement was confirmed with a range-over-time measurement to the same clouds as well as compared to EMC meteorological models. Here we describe the lidar approach and the breadboard instrument, and report some early results from ongoing field experiments.

Cremons, Daniel R.↗

NASA NTRS1993

Aspects of modelling the tectonics of large volcanoes on the terrestrial planets

Analytic solutions for the responses of planetary lithospheres to volcanic loads have been used to model faulting and infer elastic plate thicknesses. Predictions of the distribution of faulting around volcanic loads, based on the application of Anderson's criteria for faulting to the results of the models, do not agree well with observations. Such models do not give the stress state in the load itself, but only suggest a state of horizontal compressive stress there. Further, these models have considered only the effect of an instantaneously emplaced load. They do not address the time evolution of stresses, nor do they consider the effect of a load which grows. A finite element approach allows us to assign elements to the load itself, and thus permits calculation of the stress state and stress history within the edifice. The effects of episodic load growth can also be treated. When these effects are included, models give much better agreement with observations. We use the finite element code TECTON to construct axisymmetric models of volcanoes resting on an elastic lithospheric plate overlying a viscoelastic asthenosphere. We have implemented time-dependent material properties in order to simulate incremental volcano growth. The viscoelastic layer was taken to extend to a sufficient depth so that a rigid lower boundary has no significant influence on the results. The code first calculates elastic deformations and stresses and then determines the time-dependent viscous deformations and stresses. Time in the model scales as the Maxwell time tau(m) in the asthenosphere. We consider a volcano 25 km in height and 200 km in radius on an elastic lithosphere 40 km thick (parameters approximately appropriate to Ascraeus Mons). The volcano consists of three load increments applied at intervals of 1000 tau(m). Contours of maximum deviatoric stress in the fully-grown edifice at the conclusion of flexure (t = 3000 tau(m)) are shown.

Mcgovern, Patrick J.↗

NASA NTRS1997

Experimental and Computational Studies of the Flow Over a Sting Mounted Planetary Probe Configuration

This paper summarizes the results of a series of experimental studies in the LENS shock tunnel and computations with DSMC and Navier Stokes codes which have been made to examine the aerothermal and flowfield characteristics of the flow over a sting-supported planetary probe configuration in hypervelocity air and nitrogen flows. The experimental program was conducted in the LENS hypervelocity shock tunnel at total enthalpies of 5and 10 MJkg for a range of reservoir pressure conditions from 70 to 500 bars. Heat transfer and pressure measurements were made on the front and rear face of the probe and along the supporting sting. High-speed and single shot schlieren photography were also employed to examine the flow over the model and the time to establish the flow in the base recirculation region. Predictions of the flowfield characteristics and the distributions of heat transfer and pressure were made with DSMC codes for rarefied flow conditions and with the Navier-Stokes solvers for the higher pressure conditions where the flows were assumed to be laminar. Analysis of the time history records from the heat transfer and pressure instrumentation on the face of the probe and in the base region indicated that the base flow was fully established in under 4 milliseconds from flow initiation or between 35 and 50 flow lengths based on base height. The measurements made in three different tunnel entries with two models of identical geometries but with different instrumentation packages, one prepared by NASA Langley and the second prepared by CUBRC, demonstrated good agreement between heat transfer measurements made with two different types of thin film and coaxial gage instrumentation. The measurements of heat transfer and pressure to the front face of the probe were in good agreement with theoretical predictions from both the DSMC and Navier Stokes codes. For the measurements made in low density flows, computations with the DSMC code were found to compare well with the pressure and heat transfer measurements on the sting, although the computed heat transfer rates in the recirculation region did not exhibit the same characteristics as the measurements. For the 10MJkg and 500 bar reservoir match point condition, the measurements and heat transfer along the sting from the first group of studies were in agreement with the Navier Stokes solutions for laminar conditions. A similar set of measurements made in later tests where the model was moved to a slightly different position in the test section indicated that the boundary layer in the reattachment compression region was close to transition or transitional where small changes in the test environment can result in larger than laminar heating rates. The maximum heating coefficients on the sting observed in the present studies was a small fraction of similar measurements obtained at nominally the same conditions in the HEG shock tunnel, where it is possible for transition to occur in the base flow, and in the low enthalpy studies conducted in the NASA Langley high Reynolds number Mach 10 tunnel where the base flow was shown to be turbulent. While the hybrid Navier- StokedDMSC calculations by Gochberg et al. (Reference 1) suggested that employing the Navier- Stokes calculations for the entire flowfield could be seriously in error in the base region for the 10 MJkg, 500 bar test case, similar calculations performed by Cornell, presented here, do not.

Holden, Michael S.↗

NASA NTRS1999

Waves in Radial Gravity Using Magnetic Fluid

Terrestrial laboratory experiments studying various fluid dynamical processes are constrained, by being in an Earth laboratory, to have a gravitational body force which is uniform and unidirectional. Therefore fluid free-surfaces are horizontal and flat. Such free surfaces must have a vertical solid boundary to keep the fluid from spreading horizontally along a gravitational potential surface. In atmospheric, oceanic, or stellar fluid flows that have a horizontal scale of about one-tenth the body radius or larger, sphericity is important in the dynamics. Further, fluids in spherical geometry can cover an entire domain without any sidewall effects, i.e. have truly periodic boundary conditions. We describe spherical body-force laboratory experiments using ferrofluid. Ferrofluids are dilute suspensions of magnetic dipoles, for example magnetite particles of order 10 nm diameter, suspended in a carrier fluid. Ferrofluids are subject to an additional body force in the presence of an applied magnetic field gradient. We use this body force to conduct laboratory experiments in spherical geometry. The present study is a laboratory technique improvement. The apparatus is cylindrically axisymmetric. A cylindrical ceramic magnet is embedded in a smooth, solid, spherical PVC ball. The geopotential field and its gradient, the body force, were made nearly spherical by careful choice of magnet height-to-diameter ratio and magnet size relative to the PVC ball size. Terrestrial gravity is eliminated from the dynamics by immersing the "planet" and its ferrofluid "ocean" in an immiscible silicone oil/freon mixture of the same density. Thus the earth gravity is removed from the dynamics of the ferrofluid/oil interface and the only dynamically active force there is the radial magnetic gravity. The entire apparatus can rotate, and waves are forced on the ferrofluid surface by exterior magnets. The biggest improvement in technique is in the wave visualization. Fluorescing dye is added to the oil/freon mixture and an argon ion laser generates a horizontal light that can be scanned vertically. Viewed from above, the experiment is a black circle with wave deformations surrounded by a light background. A contour of the image intensity at any light sheet position gives the surface of the ferrofluid "ocean" at that "latitude". Radial displacements of the waves as a function of longitude are obtained by subtracting the contour line positions from a no-motion contour at that laser sheet latitude. The experiments are run by traversing the forcing magnet with the laser sheet height fixed and images are frame grabbed to obtain a time-series at one latitude. The experiment is then re-run with another laser-sheet height to generate a full picture of the three-dimensional wave structure in the upper hemisphere of the ball as a function of time. We concentrate here on results of laboratory studies of waves that are important in Earth's atmosphere and especially the ocean. To get oceanic scaling in the laboratory, the experiment must rotate rapidly (4-second rotation period) so that the wave speed is slow compared to the planetary rotation speed as in the ocean. In the Pacific Ocean, eastward propagating Kelvin waves eventually run into the South American coast. Theory predicts that some of the wave energy should scatter into coastal-trapped Kelvin waves that propagate north and south along the coast. Some of this coastal wave energy might then scatter into mid-latitude Rossby waves that propagate back westward. Satellite observations of the Pacific Ocean sea-surface temperature and height seem to show signatures of westward propagating mid-latitude Rossby waves, 5 to 10 years after the 1982-83 El Nino. The observational data is difficult to interpret unambiguously owing to the large range of motions that fill the ocean at shorter timescales. This series of reflections giving eastward, north- ward, and then westward traveling waves is observed cleanly in the laboratory experiments, confirming the theoretical expectations.

Ohlsen, D. R.↗