Search NASA⌕ Search

SEARCH · Search NASA

Results for “smoothers”

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 199 records · Page 11

Venus gravity

The anomalous gravity field of Venus shows high correlation with surface features revealed by radar. We extract gravity models from the Doppler tracking data from the Pioneer Venus Orbiter (PVO) by means of a two-step process. In the first step, we solve the nonlinear spacecraft state estimation problem using a Kalman filter-smoother. The Kalman filter was evaluated through simulations. This evaluation and some unusual features of the filter are discussed. In the second step, we perform a geophysical inversion using a linear Bayesian estimator. To allow an unbiased comparison between gravity and topography, we use a simulation technique to smooth and distort the radar topographic data so as to yield maps having the same characteristics as our gravity maps. The maps presented cover 2/3 of the surface of Venus and display the strong topography-gravity correlation previously reported. The topography-gravity scatter plots show two distinct trends.

Reasenberg, Robert D.↗

On the stability analysis of approximate factorization methods for 3D Euler and Navier-Stokes equations

The convergence characteristics of various approximate factorizations for the 3D Euler and Navier-Stokes equations are examined using the von-Neumann stability analysis method. Three upwind-difference based factorizations and several central-difference based factorizations are considered for the Euler equations. In the upwind factorizations both the flux-vector splitting methods of Steger and Warming and van Leer are considered. Analysis of the Navier-Stokes equations is performed only on the Beam and Warming central-difference scheme. The range of CFL numbers over which each factorization is stable is presented for one-, two-, and three-dimensional flow. Also presented for each factorization is the CFL number at which the maximum eigenvalue is minimized, for all Fourier components, as well as for the high frequency range only. The latter is useful for predicting the effectiveness of multigrid procedures with these schemes as smoothers. Further, local mode analysis is performed to test the suitability of using a uniform flow field in the stability analysis. Some inconsistencies in the results from previous analyses are resolved.

Demuren, A. O.↗

Preliminary results of polarization signatures for glacial moraines in the Mono Basin, Eastern Sierra Nevada

The valleys of the Mono Basin contain several sets of lateral and terminal moraines representing multiple stages of glaciation. The semi-arid climate with slow weathering rates preserved sequences of nested younger moraines within older ones. There is a well established relative chronology and recently exposure dating provided a new set of numerical dates. The moraines span the late Wisconsin (11-25 ka) to the Illinoian (130-190 ka) glaciations. The Mono Basin area was used as a 'calibration site' to establish remote dating techniques for eventual transfer to the more inaccessible but geomorphically and climatically similar moraines of the South American Andes Mountains. Planned polarimetric synthetic aperture radar (SAR) imagery acquired by JPL AIRSAR (South American Campaign) and SIR-C (Andes super-site) are analyzed to establish chronologies of previously undated moraine sequences in a study of Pleistocene climatic change in the Southern Hemisphere. The dry climate and sparse vegetation is also favorable for correlation of ground surface roughness with radar polarization signature. The slow weathering processes acting over thousands of years reduce the size, frequency, and angularity of surface boulders while increasing soil development on the moraines. Field observations based on this hypothesis result in relative ages consistent with those inferred from nested position within the valley. Younger moraines, therefore, will appear rougher than the older smoother moraines at scales measurable at AIRSAR wavelengths. Previously documented effects of ground surface roughness on polarization signatures suggest that analysis of moraine polarization signatures can be useful for relative dating. The technique may be extended to predict numerical ages. The data set reported were acquired on 8 Sep. 1989 with the JPL Airborne SAR (AIRSAR) collecting polarimetric imagery at C- (5.6 cm), L- (24 cm), and P-band (68 cm) with a flight-line parallel to the strike of the mountains. Phase calibration was performed on the analyzed scene by setting the co-phase of a smooth lake to zero as described. Absolute amplitude calibration was not possible because corner reflectors were not deployed.

Forster, Richard R.↗

Detecting surface roughness effects on the atmospheric boundary layer via AIRSAR data: A field experiment in Death Valley, California

The part of the troposphere influenced by the surface of the earth is termed the atmospheric boundary layer. Flow within this layer is influenced by the roughness of the surface; rougher surfaces induce more turbulence than smoother surfaces and, hence, higher atmospheric transfer rates across the surface. Roughness elements also shield erodible particles, thus decreasing the transport of windblown particles. Therefore, the aerodynamic roughness length (z(sub 0)) is an important parameter in aeolian and atmospheric boundary layer processes as it describes the aerodynamic properties of the underlying surface. z(sub 0) is assumed to be independent of wind velocity or height, and dependent only on the surface topography. It is determined using in situ measurements of the wind speed distribution as a function of height. For dry, unvegetated soils the intensity of the radar backscatter (sigma(sup 0)) is affected primarily by surface roughness at a scale comparable with the radar wavelength. Thus, both wind and radar respond to surface roughness variations on a scale of a few meters or less. Greeley showed the existence of a correlation between z(sub 0) and sigma(sup 0). This correlation was based on measurements over lava flows, alluvial fans, and playas in the southwest deserts of the United States. It is shown that the two parameters behave similarly also when there are small changes over a relatively homogeneous surface.

Blumberg, Dan G.↗

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗

A multigrid method for steady Euler equations on unstructured adaptive grids

A flux-difference splitting type algorithm is formulated for the steady Euler equations on unstructured grids. The polynomial flux-difference splitting technique is used. A vertex-centered finite volume method is employed on a triangular mesh. The multigrid method is in defect-correction form. A relaxation procedure with a first order accurate inner iteration and a second-order correction performed only on the finest grid, is used. A multi-stage Jacobi relaxation method is employed as a smoother. Since the grid is unstructured a Jacobi type is chosen. The multi-staging is necessary to provide sufficient smoothing properties. The domain is discretized using a Delaunay triangular mesh generator. Three grids with more or less uniform distribution of nodes but with different resolution are generated by successive refinement of the coarsest grid. Nodes of coarser grids appear in the finer grids. The multigrid method is started on these grids. As soon as the residual drops below a threshold value, an adaptive refinement is started. The solution on the adaptively refined grid is accelerated by a multigrid procedure. The coarser multigrid grids are generated by successive coarsening through point removement. The adaption cycle is repeated a few times. Results are given for the transonic flow over a NACA-0012 airfoil.

Riemslagh, Kris↗

Multigrid solution of the Navier-Stokes equations on highly stretched grids with defect correction

Relaxation-based multigrid solvers for the steady incompressible Navier-Stokes equations are examined to determine their computational speed and robustness. Four relaxation methods with a common discretization have been used as smoothers in a single tailored multigrid procedure. The equations are discretized on a staggered grid with first order upwind used for convection in the relaxation process on all grids and defect correction to second order central on the fine grid introduced once per multigrid cycle. A fixed W(1,1) cycle with full weighting of residuals is used in the FAS multigrid process. The resulting solvers have been applied to three 2D flow problems, over a range of Reynolds numbers, on both uniform and highly stretched grids. In all cases the L(sub 2) norm of the velocity changes is reduced to 10(exp -6) in a few 10's of fine grid sweeps. The results from this study are used to draw conclusions on the strengths and weaknesses of the individual relaxation schemes as well as those of the overall multigrid procedure when used as a solver on highly stretched grids.

Sockol, Peter M.↗

Mapping robust parallel multigrid algorithms to scalable memory architectures

The convergence rate of standard multigrid algorithms degenerates on problems with stretched grids or anisotropic operators. The usual cure for this is the use of line or plane relaxation. However, multigrid algorithms based on line and plane relaxation have limited and awkward parallelism and are quite difficult to map effectively to highly parallel architectures. Newer multigrid algorithms that overcome anisotropy through the use of multiple coarse grids rather than relaxation are better suited to massively parallel architectures because they require only simple point-relaxation smoothers. In this paper, we look at the parallel implementation of a V-cycle multiple semicoarsened grid (MSG) algorithm on distributed-memory architectures such as the Intel iPSC/860 and Paragon computers. The MSG algorithms provide two levels of parallelism: parallelism within the relaxation or interpolation on each grid and across the grids on each multigrid level. Both levels of parallelism must be exploited to map these algorithms effectively to parallel architectures. This paper describes a mapping of an MSG algorithm to distributed-memory architectures that demonstrates how both levels of parallelism can be exploited. The result is a robust and effective multigrid algorithm for distributed-memory machines.

Overman, Andrea↗

Modelling TOVS radiances of synoptic systems

Two operational analysis models, produced by ECMWF and NMC, are compared against observed satellite radiances by converting the analyzed variables into synthetic radiances using a radiative transfer model (RTM). Observed TOVS satellite radiances are used as a ground truth, since they are the only source of continuous, synoptic-scale resolution observations available over the data-sparse tropical Pacific Ocean. It is assumed that if the analysis models correctly simulate atmospheric conditions, the observed and modelled radiance fields will have similar patterns and magnitudes. In each of seven TOVS channels examined, the mean temperature difference between the models was small, but the differences between the models and the observations were large. The mean temperature difference within each channel was larger than the standard deviation in most cases. The moisture channels of the RTM were sensitive to changes in moisture content, but none of the channels were significantly affected by changes in temperature. Single-day modelled radiance analyses contain synoptic information similar to the satellite observations. ECMWF is often too noisy, producing features not verified by the satellite; NMC is often too smooth. NMC more closely resembles the TOVS contours. Modelled water vapor channel analyses resemble observed monthly and composited observations well, but do not accurately reproduce the highly-variable single-day fields. ECMWF generally has more accurate gradients. Neither model simulates the tropics well. The water vapor channels, especially channel eleven (weighted at 700 mb), most closely resembled the tropical plume composite described by McGuirk and Ulsh. The only plume evidence found in the thermal or microwave channels, in either the observations or the models, was a composite trough situated in the correct location, extending from the mid-latitudes into the subtropics. This feature was not readily discernable until the definition and mature stages. Additional information besides the modelled temperatures would be required to produce an accurate forecast for the tropics. Model patterns more closely resembled the observations in the quiescent atmosphere, though the magnitudes were not any closer. ECMWF simulates the TOVS quiescent state better than NMC. The moisture channels continue to present the most realistic patterns. The models contained more detail in the tropics than was verified by the satellite observations, apparently trying to force an active ITCZ. Lag correlations showed TOVS to be more persistent than the models. TOVS also showed a smoother lag trend, but all products appeared to have adequate ability to react to changes in the atmosphere over time.

Coe, Thomas Eddy↗

A comparison of TSS and TRASYS in form factor calculation

As the workstation and personal computer become more popular than a centralized mainframe to perform thermal analysis, the methods for space vehicle thermal analysis will change. Already, many thermal analysis codes are now available for workstations, which were not in existence just five years ago. As these changes occur, some organizations will adopt the new codes and analysis techniques, while others will not. This might lead to misunderstandings between thermal shops in different organizations. If thermal analysts make an effort to understand the major differences between the new and old methods, a smoother transition to a more efficient and more versatile thermal analysis environment will be realized.

Golliher, Eric↗

Improved solution accuracy for TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajectory Determination System (GTDS) and an extended Kalman filter estimation system to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. GTDS is the operational orbit determination system used by the FDD in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. The extended Kalman filter was implemented in an orbit determination analysis prototype system, closely related to the Real-Time Orbit Determination System/Enhanced (RTOD/E) system. In addition, the Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generated an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the geodynamics (GEODYN) orbit determination system with laser ranging and Doppler Orbitography and Radiopositioning integrated by satellite (DORIS) tracking measurements. The TOPEX/Poseidon trajectories were estimated for November 7 through November 11, 1992, the timeframe under study. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch-least-squares solutions were assessed based on the solution residuals, while the sequential solutions were assessed based on primarily the estimated covariances. The batch-least-squares and sequential orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 2 meters for the batch-least-squares and less than 13 meters for the sequential estimation solutions. After the sequential estimation solutions were processed with a smoother algorithm, position differences with POD orbit solutions of less than 7 meters were obtained. The differences among the POD, GTDS, and filter/smoother solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C. E.↗

Improved solution accuracy for Landsat-4 (TDRSS-user) orbit determination

This paper presents the results of a study to compare the orbit determination accuracy for a Tracking and Data Relay Satellite System (TDRSS) user spacecraft, Landsat-4, obtained using a Prototype Filter Smoother (PFS), with the accuracy of an established batch-least-squares system, the Goddard Trajectory Determination System (GTDS). The results of Landsat-4 orbit determination will provide useful experience for the Earth Observing System (EOS) series of satellites. The Landsat-4 ephemerides were estimated for the January 17-23, 1991, timeframe, during which intensive TDRSS tracking data for Landsat-4 were available. Independent assessments were made of the consistencies (overlap comparisons for the batch case and convariances for the sequential case) of solutions produced by the batch and sequential methods. The filtered and smoothed PFS orbit solutions were compared with the definitive GTDS orbit solutions for Landsat-4; the solution differences were generally less than 15 meters.

Oza, D. H.↗

System concept definition of the Grumman superconducting Electromagnetic Suspension (EMS) Maglev design

Grumman, under contract to the Army Corps of Engineers, completed a System Concept Definition (SCD) study to design a high-speed 134 m/s (300 m.p.h.) magnetically levitated (Maglev) transportation system. The primary development goals were to design a Maglev that is safe, reliable, environmentally acceptable, and low-cost. The cost issue was a predominant one, since previous studies have shown that an economically viable Maglev system (one that is attractive to investors for future models of passenger and/or freight transportation) requires a cost that is about $12.4 M/km ($20 Million per mile). The design is based on the electromagnetic suspension (EMS) system using superconducting iron-core magnets mounted along both sides of the vehicle. The EMS system has several advantages compared to the electrodynamic suspension (EDS) Maglev systems such as low stray magnetic fields in the passenger cabin and the surrounding areas, uniform load distribution along the full length of the vehicle, and small pole pitch for smoother propulsion and ride comfort. It is also levitated at all speeds and incorporates a wrap-around design of safer operation. The Grumman design has all the advantages of an EMS system identified above, while eliminating (or significantly improving) drawbacks associated with normal magnet powered EMS systems. Improvements include larger gap clearance, lighter weight, lower number of control servos, and higher off line switching speeds. The design also incorporates vehicle tilt (plus or minus 9 deg) for higher coordinated turn and turn out speed capability.

Proise, M.↗

A Hubble Space Telescope ultraviolet spectrum of SN 1993J

We obtained a Hubble Space Telescope (HST) ultraviolet spectrum of the Type II supernova SN 1993J in M81 on 1993 April 15. The approximately 1650-2900 A region is smoother than observed for SN 1987A and SN 1992A and lacks strong P Cygni absorptions caused by iron peak element lines. Synthetic spectra calculated using a parameterized Local Thermodynamic Equilibrium (LTE) procedure and a simple model atmosphere do not fit the UV region. Radio observations suggest that SN 1993J is embedded in a thick circumstellar envelope. The UV spectra of other supernovae that are believed to have thick circumstellar envelopes also have approximately 1650-2900 A regions lacking in strong P Cygni absorptions. Interaction of supernova ejecta and circumstellar matter may cause the smooth UV spectrum. If so, UV observations of supernovae will provide insight into the circumstellar environment of the supernova progenitors.

Jeffery, David J.↗

Bell Regio, Venus: Integration of remote sensing data and terrestrial analogs for geologic analysis

The geology and surface morphology of Bell Regio (18-42 deg N, 32-58 deg E) are investigated using a combination of Magellan, Venera, and analogous terrestrial data. The properties of surface units are compared to either direct terrestrial analog measurements or to the behaviors predicted by theoretical models. Five major volcanic sources are identified from geologic mapping (Tepev Mons, Nefertiti corona, a large shield volcano east of Tepev, and two small edifices southeast of Tepev). The volcano Api Mons lies northeast of the main Bell uplift. The oldest volcanic units are associated with an extensive low shield volcano east of Tepev Mons and a small edifice southeast of Tepev. The annular flow apron of Tepev Mons formed next, with volcanism at a second small edifice on the southeast flank of Tepev Mons producing the youngest flow units. Comparisons between Magellan data, terrestrial radar images, and field topography profiles suggest that only three units resemble terrestrial a'a flows; the remainder are consistent with smoother pahoehoe-type surfaces. This suggests that most of the flow units were erupted at relatively low volume effusion rates (less than 100 cu m/sec) over long periods of time or had very low viscosities

Campbell, Bruce A.↗

Space and time aliasing structure is monthly mean polar-orbiting satellite data

Monthly mean wind fields from the European Remote Sensing Satellite (ERS1) scatterometer are presented. A banded structure which resembles the satellite subtrack is clearly and consistently apparent in the isotachs as well as the u and v components of the routinely produced fields. The structure also appears in the means of data from other polar-orbiting satellites and instruments. An experiment is designed to trace the cause of the banded structure. The European Centre for Medium-Range Weather Forecast (ECMWF) gridded surface wind analyses are used as a control set. These analyses are also sampled with the ERS1 temporal-spatial samplig pattern to form a simulated scatterometer wind set. Both sets are used to create monthly averages. The banded structures appear in the monthly mean simulated data but do not appear in the control set. It is concluded that the source of the banded structure lies in the spatial and temporal sampling of the polar-orbiting satellite which results in undersampling. The problem involves multiple timescales and space scales, oversampling and under-sampling in space, aliasing in the time and space domains, and preferentially sampled variability. It is shown that commonly used spatial smoothers (or filters), while producing visually pleasing results, also significantly bias the true mean. A three-dimensional spatial-temporal interpolator is designed and used to determine the mean field. It is found to produce satisfactory monthly means from both simulated and real ERS1 data. The implications to climate studies involving polar-orbiting satellite data are discussed.

Zeng, Lixin↗

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

Controlling Flows Of Two Ingredients For Spraying

Closed-loop servo control subsystem incorporated, as modification, into system controlling flows of two ingredients mixed and sprayed to form thermally insulating foams on large tanks. Provides steady flows at specified rates. Foams produced smoother and of higher quality. Continued use of system results in substantial reduction in cost stemming from close control of application of foam and consequent reduced use of material.

Chandler, Huel H.↗