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 775 records · Page 43

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↗

Stochastic Microgrid Scheduling With Chance‐Constrained Resilience Consideration

Traditionally, it is assumed that microgrids transition seamlessly from grid‐connected operation to islanded mode in the event of sudden main grid outages. In reality, the islanding process, especially unintentional islanding, is rarely seamless. Instead, it is subject to voltage and frequency fluctuations caused by the instantaneous disconnection of the point of common coupling (PCC) switch, variations in loads and renewable generation output and even the protection tripping of distributed energy resources (DERs). To mitigate these fluctuations and facilitate a smooth islanding process, we propose a stochastic microgrid scheduling model that incorporates chance‐constrained resilience measures. Specifically, the resilience measure is defined as the probability of successful islanding (PSI), that is, the probability that a microgrid can mitigate the generation‐demand imbalance caused by the disconnection of the PCC switch, variations in load and renewable generation and DER tripping. This measure is modelled using chance constraints. Unlike existing reliability and resilience indices, which typically neglect the possibility of microgrid/DER failure under extreme events and assume their survival while primarily focussing on reducing impact duration or magnitude, the proposed PSI‐based framework explicitly addresses microgrid and DER survival during the islanding transition. The formulated nonlinear chance constraints are approximated using a multiinterval approach and equivalently represented as a mixed‐integer linear programming (MILP) formulation. Case study results validate the proposed method, showing that the PSI estimation error is reduced to less than 8%, compared to approximately 28% with existing methods. Various sensitivity analyses on the DER tripping rate and PSI settings were performed to validate the robustness of the proposed method. In particular, the necessity of accounting for DER tripping in the PSI calculation was demonstrated.

chance constrained optimization↗

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↗

Hybrid Curve Fitting For Reducing Motion Commands In Object Construction

Existing slicing software for additive manufacturing typically requires a triangulated mesh as its input. Triangulated meshes are approximate representations of exact CAD models. Despite the loss in dimensional accuracy, triangulated meshes are used because they are computationally easier to cross-section and offset than the exact geometry in CAD format. When a triangulated object is prepared, the resulting machine instructions include only linear motion commands. Numerous modern motion controllers can move in arc and spline motions; however, the absence of slicing software that supports curvature prevents these commands from being leveraged. To address this limitation, this paper presents a method for the hybrid reconstruction of arcs and splines as a post-processing step to traditional slicing. This method can greatly reduce the number of motion commands required to construct an object by printing smooth curved surfaces. This concise representation of tool-pathing allows for more even extrusion and is computed without a major impact on slicing time.

Wade, Charles↗

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.↗

A Spatially Coherent Global Soil Moisture Product with Improved Temporal Resolution

Global soil moisture products that are completely independent of any type of ancillary data and solely rely on satellite observations are presented. Additionally, we further develop an existing downscaling technique that enhances the spatial resolution of such products to approximately 11 km. These products are based on internal modules of the Land Parameter Retrieval Model (LPRM), an algorithm that uses the radiative transfer equation to link soil moisture, vegetation optical depth and land surface temperature to observed brightness temperatures. The soil moisture product that is independent of any type of ancillary data uses the internally calculated dielectric constant as a soil moisture proxy. This data product is not influenced by errors associated with coarse-scale global soil property maps or by any other type of forcing (e.g. re-analysis) data and is therefore solely based on satellite microwave observations. The second step builds upon recent developments to increase the spatial resolution of the LPRM retrievals using a smoothing filter downscaling method. With this method we can attain a spatial resolution that can be more useful at the scale of local and regional hydrological studies as well. The steps presented in this paper were applied to observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). The newly derived data sets were validated using ground-based observations from the International Soil Moisture Network (ISMN). The internally calculated dielectric constant product results in significantly more days with valid retrievals than the original soil moisture data products, in particular over arid regions. The dielectric constant product resulted in similar correlations with in situ data as the original soil moisture data product. Together, these findings demonstrate the usefulness of this new dielectric constant product for the hydrological modeling community and climate studies. A case study on the Australian Fitzroy catchment demonstrated that the downscaled data product has a more detailed spatial description of soil moisture, especially during wet and dry conditions with more pronounced dry and wet regions within the catchment. The increased resolution data products could therefore improve runoff predictions and this study demonstrated the potential added value of a transitioning from a spatial resolution of 56 km toward a higher resolution of 11 km. The hydrological implications of these newly developed data records are not only linked to AMSR-E satellite data, but also to the next generation Soil Moisture Active and Passive (SMAP) mission where a 9 km spatial resolution is the target resolution for satellite soil moisture products. The new data products will not replace the current LPRM products, but will be added to the existing array of data products and will become publicly available through our data portals.

spatial↗

16-38 micron spectroscopy of NGC 7027

Medium-resolution spectroscopy of the planetary nebula NGC 7027 over the wavelength range from 16 to 38 microns is reported which was performed from NASA's Kuiper Airborne Observatory. It is found that within the errors, the spectrum is smooth and has the shape of a dilute 130-K blackbody but falls below that curve at the longer wavelengths. The absence of two theoretically predicted features in the 16-40-micron band, a forbidden fine-structure line of O IV at 25.91 microns and a broad emission feature in the 22-35-micron region due to carbonate materials, is discussed; a lower limit of 25,000 per cu cm is placed on the electron density in the ionized region. A simple model involving small graphite grains at a temperature around 90 K is proposed to explain the observed featureless continuum. This model is shown to require a dust mass of about 0.025 solar mass and a gas mass of approximately 1 solar mass, which imply that most of this mass would be associated with the molecular cloud surrounding the ionized region.

Mccarthy, J. F.↗

Electronic Structure Correspondence of Singlet-Triplet Scale Separation in Strained Sr2RuO4

At a temperature of roughly 1 K, Sr 2 RuO 4 undergoes a transition from a normal Fermi liquid to a superconducting phase. Even while the former is relatively simple and well understood, the superconducting state has not even been understood after 25 years of study. More recently, it has been found that critical temperatures can be enhanced by the application of uniaxial strain, up to a critical strain, after which it falls off. In this work, we take an “instability” approach and seek divergences in susceptibilities. This provides an unbiased way to distinguish tendencies to competing ground states. We show that in the unstrained compound, the singlet and triplet instabilities of the normal Fermi liquid phase are closely spaced. Under uniaxial strain, electrons residing on all orbitals contributing to the Fermiology become more coherent, while the electrons of the Ru-dxy character become heavier, and the electrons of the Ru d x z , y z characters become lighter. In the process, Im χ(q,ω) increases rapidly around q = ( 0.3 , 0.3 , 0 ) 2 π / a and q = ( 0.5 , 0.25 , 0 ) 2 π / a , while it gets suppressed at all other commensurate vectors, in particular at q = 0, which is essential for spin-triplet superconductivity. We observe that the magnetic anisotropy under strain drops smoothly, which is concomitant with the increment in singlet instability. Thus, the triplet superconducting instability remains the lagging instability of the system, and the singlet instability enhances under strain, leading to a large energy-scale separation between these competing instabilities. However, since this happens even without spin-orbit coupling, we believe it is primarily the enhancement in the spin fluctuation glue around quasi-anti-ferromagnetic vectors that drives the Cooper pairing instead of the magnetic anisotropy. At large strain, an instability to a spin density wave overtakes the superconducting one. The analysis relies on a high-fidelity, ab initio description of the one-particle properties and two-particle susceptibilities, based on the quasiparticle self-consistent GW approximation augmented by dynamical mean field theory. This approach is described and its high fidelity confirmed by comparing to observed one- and two-particle properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Influence of Tire Tread Pattern and Runway Surface Condition on Braking Friction and Rolling Resistance of a Modern Aircraft Tire

A series of taxiing tests was conducted at the Langley landing­ loads track with both braked and unbraked (freely rolling) single and tandem wheels equipped with 32x8.8 type VII aircraft tires of dif­ferent tread designs to obtain data on tire and braking characteristics during operation on dry and on contaminated concrete and asphalt run­ ways. Contaminants used were water, slush, JP-4 jet fuel, and organic and detergent fire-extinguishing foams. Forward velocities for the tests ranged from approximately 13 to 104 knots. Vertical loads of approximately 9,000 to 22,000 pounds and tire inflation pressures of 85 to 350 pounds per square inch were used. Results indicated that the unbraked tire rolling resistance increased with increasing forward velocity on dry and on contaminated runway surfaces. Peak tire-ground friction coefficients developed during wheel braking decreased rapidly with increasing velocity on contaminated runways but remained relatively unchanged on dry runways as the forward velocity was increased. Dry-runway friction coeffi­cients were found to be relatively insensitive to tire tread pattern. However, the magnitude of the friction coefficients developed by tires on contaminated runways was extremely sensitive to the tire tread pat­ tern used, with circumferential-groove treads developing the highest values of friction coefficient, and smooth and dimple treads the lowest values for the tread patterns and runway conditions investigated.

Horne, Walter B.↗

Vertical temperature and density patterns in the Arctic mesosphere analyzed as gravity waves

Three series of rocket soundings including pitot soundings, grenade soundings, and paired pitot-grenade soundings, were conducted from high latitude sites during winter. Temperature and wind profiles and one density profile were observed independently to obtain the thermodynamic structure, the wind structure, and thus their interdependence in the mesosphere. Temperature profiles from all soundings in each series were averaged, and a smooth curve (or series of smooth curves) drawn through the points. A hydrostatic atmosphere based on the average, measured temperature profile was computed, and deviations from the mean atmosphere were analyzed in terms of gravity wave theory. The vertical wavelengths of the deviations were 10-20 km, and the wave amplitudes slowly increased with height. The experimental data were matched by calculated gravity waves having a period ranging between 15 and 80 minutes and horizontal wavelengths of 60 to 280 km. The interpretation is generally consistent with the results of others who have studied gravity-acoustic waves in the atmosphere. The wind measurements are consistent with the thermodynamic measurements. The results also suggest that gravity waves traveled from east to west with a horizontal phase velocity of approximately 60 m/sec.

Eberstein, I. J.↗

Atmospheric Infrared Sounder

The microwave 'first-guess' algorithm was run on the cloudy test simulations. Eight datasets were considered in the cloudy test, comprising approximately 360 retrievals, of which one failed to converge. Retrievals were done on the AMSU-A grid. Examination of the true profiles (provided for the A, C and D tracks) showed numerous cases of very pronounced temperature inversion layers in the troposphere which the retrieval does not have enough vertical resolution to reproduce. A typical example with an inversion layer near 700 mbar is shown in Figure 1. The inversion layers also exhibit strong vertical gradients of water vapor which are not resolved in the retrieval. (Water vapor volume density is given in g/sq cm per layer. Layer thickness is 20 mb from 200 to 400 mb, and 25 mb from 400 to 1,000 mb.) The retrievals do reproduce the overall smoothed shape of the profiles, and therefore as a first guess should be within the range of linear methods for IR retrievals using AIRS.

Rosenkranz, Philip W.↗

ENO-Osher schemes for Euler equations

The combination of the Osher approximate Riemann solver for the Euler equations and various ENO schemes is discussed for one-dimensional flow. The three basic approaches, viz., the ENO scheme using primitive variable reconstruction, either with the Cauchy-Kowalewski procedure for time integration or the TVD Runge-Kutta scheme, and the flux-ENO method are tested on different shock tube cases. The ENO-Osher scheme using the Cauchy-Kowalewski procedure for time integration is found to be the most accurate and robust compared with the other methods and is also computationally efficient. The tests showed that the ENO schemes perform reasonably well, but have problems in cases where two discontinuites are close together. In that case there are not enough points in the smooth part of the flow to create a nonoscillatory interpolation.

Van Der Vegt, Jacobus J.↗