Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 613 records · Page 34

Effects of the LBV Primary's Mass-loss Rate on the 3D Hydrodynamics of eta Carinae's Colliding Winds

At the heart of eta Carinae's spectacular "Homunculus" nebula lies an extremely luminous (L(sub Total) greater than approximately 5 × 10(exp 6) solar luminosity) colliding wind binary with a highly eccentric (e approximately 0.9), 5.54-year orbit (Figure 1). The primary of the system, a Luminous Blue Variable (LBV), is our closest (D approximately 2.3 kpc) and best example of a pre-hypernova or pre-gamma ray burst environment. The remarkably consistent and periodic RXTE X-ray light curve surprisingly showed a major change during the system's last periastron in 2009, with the X-ray minimum being approximately 50% shorter than the minima of the previous two cycles1. Between 1998 and 2011, the strengths of various broad stellar wind emission lines (e.g. Halpha, Fe II) in line-of-sight (l.o.s.) also decreased by factors of 1.5 - 3 relative to the continuum2. The current interpretation for these changes is that they are due to a gradual factor of 2 - 4 drop in the primary's mass-loss rate over the last approximately 15 years1, 2. However, while a secular change is seen for a direct view of the central source, little to no change is seen in profiles at high stellar latitudes or reflected off of the dense, circumbinary material known as the "Weigelt blobs"2, 3. Moreover, model spectra generated with CMFGEN predict that a factor of 2 - 4 drop in the primary's mass-loss rate should lead to huge changes in the observed spectrum, which thus far have not been seen. Here we present results from large- (plus or minus 1620 AU) and small- (plus or minus 162 AU) domain, full 3D smoothed particle hydrodynamics (SPH) simulations of eta Car's massive binary colliding winds for three different primary-star mass-loss rates (2.4, 4.8, and 8.5 × 10(exp -4) solar mass/yr). The goal is to investigate how the mass-loss rate affects the 3D geometry and dynamics of eta Car's optically-thick wind and spatially-extended wind-wind collision (WWC) regions, both of which are known sources of observed X-ray, optical, UV, and near-IR emission and absorption. We use two domain sizes in order to better understand how the primary's mass-loss rate influences the various observables that form at different length scales. The 3D simulations provide information important for helping constrain η Car's recent mass-loss history and future state.

lavestigate↗

Progress Toward Evaluating Roughness-Resolved Aerothermodynamic Simulations for the Mars Sample Return Earth Entry System

This work pursues the efficient Navier-Stokes simulation, using NASA’s standard LAURA and DPLR simulation tools, of the flowfield and convective heating resulting from the patterned-roughness surface of the Mars Sample Return (MSR) Earth Entry System (EES). This type of simulation has been infeasible for past NASA missions. However, for EES, the combination of a patterned roughness heatshield approximated with a sinusoidal pattern, the zero degrees angle of attack nominal flow, and the supersonic boundary layer edge Mach number, all combine to make these roughness resolved grid simulations feasible. The patterned roughness and zero degrees angle of attack nominal flow allow the simulation domain to be reduced to a thin axial flowfield slice that captures a minimum of single period of the roughness pattern. With a periodic boundary condition applied along the boundaries of this slice, simulation on the reduced domain is identical to the full domain but with significant cost saving. The supersonic boundary layer edge Mach number prevents the upstream influence, which allows the smooth stagnation region to be frozen. This freezing of the stagnation region flow is required because of the numerical instability resulting from the tight spacing of the grid in the axial direction as the stagnation point is approached. Applying the developed analysis approach to a peak heating EES case results in a roughness augmentation factor for the convective heating of 1.32, which agrees well with correlations that predict values between 1.30 and 1.37. Analysis approach developed in this work provides the framework for future studies of turbulence modeling impacts and the influence of an ablating surface.

Earth Entry System↗

Gas kinematics and ionization along the extended sight line to HD 116852

We present Goddard High-Resolution Spectrograph intermediate observations of the interstellar medium toward HD 116852, a low halo star at a distance of 4.8 kpc (z = -1.3 kpc) in the direction l = 304. deg 9, b = 16.deg 1. The small science aperture observations have signal-to-noise ratios ranging from 30 to 90 and resolutions of 11 to 18 km/s (FWHM). We confirm the optical MK classification of this star through an analysis of its ultraviolet photosperic and stellar wind profiles. We detect interstellar lines of Al III, Si IV, C IV, and N V together with lines of C I, C I*, C I**, Si II, Ge II, P II, and Ni II. We convert the Mg II, P II, S II, Al III, Si IV, C IV, and N V profiles into measure of apparent column density as a function of LSR velocity. Gas scale height, velocity dispersion, and differential Galactic rotation effects govern the profile shapes. A simple computer model of the expected sight line column density profiles for the low and high ion species indicates that the gas velocity dispersions and scale heights increase as the ionization level of the gas increases. We find scale heigts H greater than or = 1 kpc for the high ions, which are comparable to the z-distance of the star, whereas we find H approximatley = 0.6 to 0.7 kpc for A1 III and H approximatley 0.1 kpc for P II and Ge II. An enhancement in the Al II profile near -15 km/s accounts for approximately 25% of the A1 III column along the sight line and probably arises within gas located approximately 500 pc below the Sagittarius-Carina spiral arm link. Portions of the broad underlying A1 III distribution are associated with the higher ionization lines, perhaps in conductive interfaces. The presence of N V and the column density ratios of Si IV, C IV, and N V favor the interpretation that much of the high ion absorption is produced by collisional ionization in gas with T = 1-3 x 10(exp 5) K. An enhancement near -35 km/s in both the Si IV and C IV profiles may be due to an outflow from the Norma spiral arm at a z-distance of about -1 kpc. The smooth decrease of the N(C IV)/N(Si IV) ratio at negative velocities may be due to an ISM ionization structure for hot gas that changes distance from the Galactic plane or to a two phase gas distribution in which the relative contribution from each phase changes with distance from the Galactic plane.

Sembach, Kenneth R.↗

Preliminary Stratigraphic Basis for Geologic Mapping of Venus

The age relations between geologic formations have been studied at 36 1000x1000 km areas centered at the dark paraboloid craters. The geologic setting in all these sites could be characterized using only 16 types of features and terrains (units). These units form a basic stratigraphic sequence (from older to younger: (1) Tessera (Tt); (2-3) Densely fractured terrains associated with coronae (COdf) and in the form of remnants among plains (Pdf); (4) Fractured and ridged plains (Pfr); (5) Plains with wrinkle ridges (Pwr); (6-7) Smooth and lobate plains (Ps/Pl); and (8) Rift-associated fractures (Fra). The stratigraphic position of the other units is determined by their relation with the units of the basic sequence: (9) Ridge bells (RB), contemporary with Pfr; (10-11) Ridges of coronae and arachnoids annuli (COar/Aar), contemporary with wrinkle ridges of Pwr; (12) Fractures of coronae annuli (COaf) disrupt Pwr and Ps/Pl; (13) Fractures (F) disrupt Pwr or younger units; (14) Craters with associated dark paraboloids (Cdp), which are on top of all volcanic and tectonic units except the youngest episodes of rift-associated fracturing and volcanism; (15-16) Surficial streaks (Ss) and surficial patches (Sp) are approximately contemporary with Cdp. These units may be used as a tentative basis for the geologic mapping of Venus including VMAP. This mapping should test the stratigraphy and answer the question of whether this stratigraphic sequence corresponds to geologic events which were generally synchronous all around the planet or whether the sequence is simply a typical sequence of events which occurred in different places at diffferent times.

Basilevsky, A. T.↗

The Acoustic Influence of Cell Depth on the Rotordynamic Characteristics of Smooth-Rotor/Honeycomb-Stator Annular Gas Seals

A two-control volume is employed for honeycomb-stator/smooth-rotor seals, with a conventional control-volume used for the through flow and a 'capacitance accumulator' model for the honeycomb cells. The control volume for the honeycomb cells is shown to cause a dramatic reduction in the effective acoustic velocity of the main flow, dropping the lowest acoustic frequency into the frequency range of interest for rotordynamics. In these circumstances, the impedance functions for the seals can not be modeled with conventional (frequency-independent) stiffness, damping, and mass coefficients. More general transfer functions are required to account for the reaction forces, and calculated here as a lead-lag term for the direct force function and a lag term for the cross-coupled function. These first order functions are simple compared to transfer functions for magnetic bearings or foundations, For synchronous response to imbalance, they can be approximated by running-speed-dependent stiffness and damping coefficients in conventional rotordynamic codes. Correct predictions for stability and transient response will require more general algorithms, pressumably using a state-space format.

Childs, Dara W.↗

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

Calibration of the Soft X-ray Telescopes (SXT) Onboard the ASTRO-H Satellite

ASTRO-H is an astrophysics satellite dedicated for non-dispersive X-ray spectroscopic study on selective celestial X-ray sources. Among the onboard instruments there are four Wolter-I X-ray mirrors of their reflectors' figure in conical approximation. Two of the four are soft X-ray mirrors, of which the energy range is from a few hundred eV to 15 keV. The focal point instruments will be a calorimeter (SXS) and a CCD camera (SXI), respectively. The mirrors were in quadrant configuration with photons being reflected consecutively in the primary and secondary stage before landing on the focal plane of 5.6 m away from the interface between the two stages. The reflectors of the mirror are made of heat-formed aluminum substrate of the thickness gauged of 152 m, 229 m, and 305 m of the alloy 5052 H-19, followed by epoxy replication on gold-sputtered smooth Pyrex cylindrical mandrels to acquire the X-ray reflective surface. The epoxy layer is 10 m nominal and surface gold layer of 0.2 m. Improvements on angular response over its predecessors, e.g. Astro-E1/Suzaku mirrors, come from error reduction on the figure, the roundness, and the grazing angle/radius mismatching of the reflecting surface, and tighter specs and mechanical strength on supporting structure to reduce the reflector positioning and the assembly errors. Each soft x-ray telescope (SXT), FM1 or FM2, were integrated from four independent quadrants of mirrors. The stray-light baffles, in quadrant configuration, were mounted onto the integrated mirror. Thermal control units were attached to the perimeter of the integrated mirror to keep the mirror within operating temperature in space. The completed instrument went through a series of optical alignment, thus made the quadrant images confocal and their optical axes in parallel to achieve highest throughput possible. Environmental tests were carried out, and optical quality of the telescopes has been confirmed. The optical and x-ray calibrations also include: angular resolution, effective area in the energy range of ~ 0.4 - 12keV, off-axis response, etc. Some of those are being carried out by our counterpart at JAXA/ISAS, Japan. We report the calibration results of the FM1 and FM2 that were obtained at Goddard Space Flight Center.

ASTRO-H Satellite↗

Fading Coronal Structure and the Onset of Turbulence in the Young Solar Wind

Above the top of the solar corona, the young, slow solar wind transitions from low-beta, magnetically structured flow dominated by radial structures to high-beta, less structured flow dominated by hydrodynamics. This transition, long inferred via theory, is readily apparent in the sky region close to 10deg from the Sun in processed, background-subtracted solar wind images. We present image sequences collected by the inner Heliospheric Imager instrument on board the Solar-Terrestrial Relations Observatory (STEREO/HI1) in 2008 December, covering apparent distances from approximately 4deg to 24deg from the center of the Sun and spanning this transition in the large-scale morphology of the wind. We describe the observation and novel techniques to extract evolving image structure from the images, and we use those data and techniques to present and quantify the clear textural shift in the apparent structure of the corona and solar wind in this altitude range. We demonstrate that the change in apparent texture is due both to anomalous fading of the radial striae that characterize the corona and to anomalous relative brightening of locally dense puffs of solar wind that we term "flocculae." We show that these phenomena are inconsistent with smooth radial flow, but consistent with the onset of hydrodynamic or magnetohydrodynamic instabilities leading to a turbulent cascade in the young solar wind.

DeForest, C. E.↗

Mineralogy and Thermal Properties of V-Type Asteroid 956 Elisa: Evidence for Diogenitic Material from the Spitzer IRS (5-35 Micrometers) Spectrum

We present the thermal infrared (5-35 micrometer) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph ("IRS"; Houck,1.R. et .11. [20041. Astrophys, 1. SuppL 154, 18-24) together with new ground-based lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 +/- 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (Hv) at that rotational phase to be 12.58 +/- 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 +/- 0.4 km with a visible albedo pv = 0.142+/- 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 +/- 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 +/- 2.8 K and beaming parameter eta = 1.16 +/- 0.05. Thermophysical modeling places a lower limit of 20 J m(exp -2)K(exp -1)s(exp -1/2) on the thermal inertia of the asteroid's surface layer (if the surface is very smooth) but more likely values fall between 30 and 150 J m(exp -2)K(exp -1)s(exp -1/2) depending on the sense of rotation. The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 micrometer reststrahlen band, the 15-16.5 micrometer Si-O-Si stretching region, and the 16-25 micrometer reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range W0(sub 2+/-1)En(sub 74+/-2)Fs(sub 24+/-1). Spectral deconvolution of the 9-12 micrometer reststrahlen features indicates that up to approximately 20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component. Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J,M., Pieters, C.M., Pratt, S.F.[1990]. J. Geophys. Res, 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite- like composition, which has very rarely been observed among asteroids.

Lim, Lucy F.↗

AR-40 AR-39 Age of an Impact-Melt Lithology in DHOFAR 961

The South Pole-Aitken (SPA) basin is the stratigraphically oldest identifiable lunar basin and is therefore one of the most important targets for absolute age-dating to help understand whether ancient lunar bombardment history smoothly declined or was punctuated by a cataclysm. The SPA basin also has another convenient property, a geochemically distinct interior, unobscured by extensive mare basalt fill. A case has been made for the possible origin of the Dhofar 961 lunar meteorite in the South Pole-Aitken (SPA) basin, based on comparing its composition with Lunar Prospector gamma-ray data for the interior of the SPA basin. Dhofar 961 contains several different impact-melt (IM) lithologies. Jolliff et al. described two classes of mafic impact-melt lithologies, one dominated by olivine (Lithology A) and the other by plagioclase (An 95-96.5) (Lithology B). Broad-beam analyses of these lithologies yielded (is) approximately 14.0 wt% FeO, 11.7 wt% MgO, and 15.4 wt% Al2O3. Lithologies A and B differ by approximately 2.5% Al2O3, 1.5% FeO and 1.5% MgO, consistent with the occurrence of olivine phenocrysts in A and plagioclase clasts in B. Both lithologies are considerably more mafic than the Apollo mafic impact-melt breccias, corresponding to olivine gabbronorite. Joy et al. used U-Pb dating to investigate phosphate fragments in the Dhofar 961 matrix and impact-melt clasts. Matrix phosphates have 4.34 to 4 Ga ages, consistent with ancient KREEP-driven magmatic episodes and Pre-Nectarian ((is) greater than 3.92 Ga). Phosphates found within Dhofar 961 crystalline impact melt breccia clasts range from 4.26 to 3.89 Ga, potentially recording events throughout the basin forming epoch of lunar history. The youngest reset ages in the Dhofar 961 sample represent an upper limit for the time of formation of the meteorite. Joy et al suggested this age represents the final impact that mixed and consolidated several generations of precursor rocks into the Dhofar meteorite group, although they note that further age dating of all the stones is required to test this hypothesis. We received a split of Dhofar 961 from R. Zeigler consisting of a large clast of IM Lithology B, with some light-colored, friable matrix clinging to the external margins of the impact-melt clast. This lithology was not present in the samples investigated by Joy et al. and thus does not have corresponding U-Pb ages on it. We created multiple subsplits of both the IM and matrix lithologies, each weighing several tens of micrograms. We conducted Ar-40 Ar-39 dating of this candidate SPA material by high-resolution step heating and comparing it with the regolith that surrounds it.

Lithology↗

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

Real-time orbit estimation for ATS-6 from redundant attitude sensors

A program installed in the ATSOCC on-line computer operates with attitude sensor data to produce a smoothed real-time orbit estimate. This estimate is obtained from a Kalman filter which enables the estimate to be maintained in the absence of T/M data. The results are described of analytical and numerical investigations into the sensitivity of Control Center output to the position errors resulting from the real-time estimation. The results of the numerical investigation, which used several segments of ATS-6 data gathered during the Sensor Data Acquisition run on August 19, 1974, show that the implemented system can achieve absolute position determination with an error of about 100 km, implying pointing errors of less than 0.2 deg in latitude and longitude. This compares very favorably with ATS-6 specifications of approximately 0.5 deg in latitude-longitude.

Englar, T. S., Jr.↗

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties↗

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM↗

The roughness of the Martian surface: A scale dependent model

In the coming decade, several lander missions to Mars are planned (e.g., MESUR Pathfinder, MESUR). One of the dangers facing planners of these missions is the rough topography observed at both Viking Lander sites. Both landing sites are ubiquitously covered with meter-scale boulders. Objects of this size pose obvious threats to soft landers, especially at Mars where the distance from Earth causes prohibitive time lags between the transmission of commands and feedback from the spacecraft. An obvious solution is to scout for a 'smooth' site prior to the landing. However, the best resolutions realizable on current and future missions (i.e., Mars Observer) are on the order of several meters. Even at this scale, boulders of 1-2 meters in size are unresolvable. Additionally, the amount of time and spacecraft resources required to search even a small area of the planet are unrealistic given other mission objectives. An alternative is to determine the 'roughness' of the surface at a subpixel scale using bidirectional reflectance observations. Much larger areas of the planet can be searched, and much of the search can easily be automated. The morphology of the martian plains observed by the Viking Landers is physically simple. The surface is covered with a layer (approximately flat lying) of aeolian sediment from which numerous outcrops of bedrock and boulders protrude. This morphology, while simple, will be difficult to characterize from orbit using traditional bidirectional reflectance models for two reasons. First, modeling the surface as facets with Gaussian or exponential slope distributions is not realistic given the morphology described above. Second, the roughness parameter is an 'average' of the roughness at scales ranging from the wavelength of light being scattered to the pixel size of the observation. Thus, there is no definite scale of roughness that can be extracted from the Hapke roughness parameter. Using the concepts of geometric and boolean models developed by several workers, we have developed a model for the bidirectional reflectance of a surface morphology comparable to that observed at the Viking Lander sites.

Shepard, M. K.↗

The Structure of the Statocyst of the Freshwater Snail Biomphalaria Glabrata (Pulmonata, Basommatophora)

The structure of the statocyst of the freshwater snail Biomphalaria glabrata has been examined by light and electron microscopy. The two statocysts are located on the dorsal-lateral side of the left and right pedal ganglion. The statocysts are spherical, fluid-filled capsules with a diameter of approximately 60 microns for young and 110 microns for adult snails. The wall of the cyst is composed of large receptor cells and many smaller supporting cells. The receptor cells bear cilia which are evenly distributed on the apical surface. The cilia have the typical 9+2 internal tubule configuration. Striate rootlets originate from the base of the basal body and run downward into the cytoplasm. Side-roots arise from one side of the basal body and a basal foot from the other. For each receptor cell, the basal foot always points to the periphery of the surface, indicating that the receptor cell is non-polarized. The receptor cells contain cytoplasmic organelles such as mitochondria, ribosomes, rough and smooth endoplasmic reticulum, compact Golgi bodies and multivesicular bodies. Supporting cells bearing microvilli are interposed between the receptor cells. The junction complex between the supporting cells and the receptor cells is composed of adherens and septate junctions, while between supporting cells only the adherens junctions are present. The static nerve arises from the lateral side of the cyst and contains axons in which parallel neurotubules and mitochondria are found. The axons arise directly from the base of the receptor cells without synapse. In the cyst lumen there are unattached statoconia. The statoconia have a plate-like or concentric membranous ring structure. Based on the morphology, the function of the statocyst in Biomphalaria is discussed.

Gao, Wenyuan↗

Hydrostatic Stress Effect on the Yield Behavior of Inconel 100

Classical metal plasticity theory assumes that hydrostatic stress has negligible effect on the yield and postyield behavior of metals. Recent reexaminations of classical theory have revealed a significant effect of hydrostatic stress on the yield behavior of various geometries. Fatigue tests and nonlinear finite element analyses (FEA) of Inconel 100 (IN100) equal-arm bend specimens and new monotonic tests and nonlinear finite element analyses of IN100 smooth tension, smooth compression, and double-edge notch tension (DENT) test specimens have revealed the effect of internal hydrostatic tensile stresses on yielding. Nonlinear FEA using the von Mises (yielding is independent of hydrostatic stress) and the Drucker-Prager (yielding is linearly dependent on hydrostatic stress) yield functions were performed. A new FEA constitutive model was developed that incorporates a pressure-dependent yield function with combined multilinear kinematic and multilinear isotropic hardening using the ABAQUS user subroutine (UMAT) utility. In all monotonic tensile test cases, the von Mises constitutive model, overestimated the load for a given displacement or strain. Considering the failure displacements or strains for the DENT specimen, the Drucker-Prager FEM s predicted loads that were approximately 3% lower than the von Mises values. For the failure loads, the Drucker Prager FEM s predicted strains that were up to 35% greater than the von Mises values. Both the Drucker-Prager model and the von Mises model performed equally-well in simulating the equal-arm bend fatigue test.

Allen, Phillip A.↗