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 379 records · Page 21

Applications of Fractal Analytical Techniques in the Estimation of Operational Scale

The observational scale and the resolution of remotely sensed imagery are essential considerations in the interpretation process. Many atmospheric, hydrologic, and other natural and human-influenced spatial phenomena are inherently scale dependent and are governed by different physical processes at different spatial domains. This spatial and operational heterogeneity constrains the ability to compare interpretations of phenomena and processes observed in higher spatial resolution imagery to similar interpretations obtained from lower resolution imagery. This is a particularly acute problem, since longterm global change investigations will require high spatial resolution Earth Observing System (EOS), Landsat 7, or commercial satellite data to be combined with lower resolution imagery from older sensors such as Landsat TM and MSS. Fractal analysis is a useful technique for identifying the effects of scale changes on remotely sensed imagery. The fractal dimension of an image is a non-integer value between two and three which indicates the degree of complexity in the texture and shapes depicted in the image. A true fractal surface exhibits self-similarity, a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution, and the slope of the fractal dimension-resolution relationship would be zero. Most geographical phenomena, however, are not self-similar at all scales, but they can be modeled by a stochastic fractal in which the scaling properties of the image exhibit patterns that can be described by statistics such as area-perimeter ratios and autocovariances. Stochastic fractal sets relax the self-similarity assumption and measure many scales and resolutions to represent the varying form of a phenomenon as the pixel size is increased in a convolution process. We have observed that for images of homogeneous land covers, the fractal dimension varies linearly with changes in resolution or pixel size over the range of past, current, and planned space-borne sensors. This relationship differs significantly in images of agricultural, urban, and forest land covers, with urban areas retaining the same level of complexity, forested areas growing smoother, and agricultural areas growing more complex as small pixels are aggregated into larger, mixed pixels. Images of scenes having a mixture of land covers have fractal dimensions that exhibit a non-linear, complex relationship to pixel size. Measuring the fractal dimension of a difference image derived from two images of the same area obtained on different dates showed that the fractal dimension increased steadily, then exhibited a sharp decrease at increasing levels of pixel aggregation. This breakpoint of the fractal dimension/resolution plot is related to the spatial domain or operational scale of the phenomenon exhibiting the predominant visible difference between the two images (in this case, mountain snow cover). The degree to which an image departs from a theoretical ideal fractal surface provides clues as to how much information is altered or lost in the processes of rescaling and rectification. The measured fractal dimension of complex, composite land covers such as urban areas also provides a useful textural index that can assist image classification of complex scenes.

Emerson, Charles W.↗

Satellite Data Assimilation

Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.

Fukumori, Ichiro↗

Current Background Noise Sources and Levels in the NASA Ames 40- by 80-Foot Wind Tunnel: A Status Report

Background noise measurements were made of the acoustic environment in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel (40x80) at NASA Ames Research Center. The measurements were acquired subsequent to the 40x80 Aeroacoustic Modernization Project, which was undertaken to improve the anechoic characteristics of the 40x80's closed test section as well as reduce the levels of background noise in the facility. The resulting 40x80 anechoic environment was described by Soderman et. al., and the current paper describes the resulting 40x80 background noise, discusses the sources of the noise, and draws comparisons to previous 40x80 background noise levels measurements. At low wind speeds or low frequencies, the 40x80 background noise is dominated by the fan drive system. To obtain the lowest fan drive noise for a given tunnel condition, it is possible in the 40x80 to reduce the fans' rotational speed and adjust the fans' blade pitch, as described by Schmidtz et. al. This idea is not new, but has now been operationally implemented with modifications for increased power at low rotational speeds. At low to mid-frequencies and at higher wind speeds, the dominant noise mechanism was thought to be caused by the surface interface of the previous test section floor acoustic lining. In order to reduce this noise mechanism, the new test section floor lining was designed to resist the pumping of flow in and out of the space between the grating slats required to support heavy equipment. In addition, the lining/flow interface over the entire test section was designed to be smoother and quieter than the previous design. At high wind speeds or high frequencies, the dominant source of background noise in the 40x80 is believed to be caused by the response of the in-flow microphone probes (required by the nature of the closed test section) to the fluctuations in the freestream flow. The resulting background noise levels are also different for probes of various diameters and types. The inflow microphone support strut is also a source of background noise but this source's impact may be minimized by careful design of the strut. In the present paper, the mechanisms mentioned above are discussed in detail. Their frequency and velocity ranges of dominance are defined and the differences between past and current facility background noise levels are presented. This paper gives valuable information for those wishing to make acoustic measurements in the 40x80. With this report and an estimate of the noise levels produced by the noise source of interest, it should be possible to determine the signal-to-noise ratios and measurement locations to successfully perform aeroacoustic testing in the NASA Ames Research Center's 40- by 80-Foot Wind Tunnel.

Allen, Christopher S.↗

The DAO Retrospective Data Assimilation System

The fixed-lag Kalman smoother of Cohn et al. (1994), or a version of a more computationally feasible approximation of it developed by Todling et al. (1998), is under implementation at the Data Assimilation Office (DAO) to build a retrospective data assimilation system (RDAS). The initial version of the RDAS uses the physical-space statistical analysis system and a modification of the incremental analysis update (IAU) procedure of Bloom et al. (1996); a follow up version of the RDAS Nill require the use of the adjoint of the DAO general circulation model (developed by Y. Yang and M Navon). The retrospective procedure is designed to produce improved analyses as well as improved assimilated fields consequently providing an improved climate representation through data assimilation.

Todling, Richardo↗

TPS Materials and Costs for Future Reusable Launch Vehicles

There is considerable interest in developing new reusable launch vehicles (RLVs) for reducing the cost of transporting payload to and from orbit. This work reviews thirteen candidate thermal protection system (TPS) options currently available for RLVs. It is useful to begin with the current Shuttle TPS layout as a reference. The nose cap and wing leading edge , which reach the highest temperatures, are made of reinforced carbon-carbon (RCC) that is protected from oxidation by an external coating (about 0.020" thick) of silicon-carbide. Most of the windward surface is 9 lb/cubic ft ceramic tiles (LI-900) with a thin (about 0.012") coating of Reaction Cured Glass (RCG). The leeward side of the vehicle is covered largely by AFRSI, a quilted ceramic blanket, and FRSI, a polyamide felt. These four materials can be considered first generation reusable TPS. Since the time of the Shuttle design, considerable progress has been made advancing TPS technologies in terms of thermal performance, robustness, and cost. For each of the major systems, a second generation ceramic TPS has been developed, tested, and characterized. Metallic-based systems have also been developed. For applications requiring RCC in the past, advanced carbon-carbon (ACC) is now available. This material has better mechanical properties, somewhat higher temperature capability to 2900F and greatly increased oxidation resistance. New carbon fiber reinforced silicon-carbide matrix composites (C/SiCs) have shown additional improvement in properties over ACC with use temperatures to 3000F and above. For rigid tiles, NASA Ames has made two significant advancements. The first is a tile substrate called Alumina Enhanced Thermal Barrier, or AETB, that incorporates alumina fibers for improved dimensional stability at high temperatures, to 2600F and above. This material can be made to densities as low as 8 lb/cubic ft. The second is a coating preparation called Toughened Uni-piece Fibrous Insulation, or TUFT, that penetrates about 0.1 in. into the tile substrate. The resulting composite, with a functionally gradient density near the surface, provides orders of magnitude increased damage resistance compared with RCG coated LI-900, with only a small weight increase. The TPS that combines these two developments is called AETB-8/TUFI and has been adopted for high damage areas on the Shuttles. Two notable developments have occur-red in flexible ceramic blanket technology. The first is aluminoborosilicate-based fibers with use temperatures of 2200F and above," in comparison to quartz and silica fiber used in AFRSI which have multi-use temperature limits of 1200 to 1400F. Blankets incorporating these new high temperature fibers are referred to as AFRSI-HT. The second is an integral weaving techniques that produces a fluted core blanket with a smoother surface and greater resistance to aero-acoustic noise, to levels as high as 170 dB. This Ames innovation is called Tailorable Advanced Blanket Insulation, or TABI. Finally, for felt-based TPS, Boeing is developing Polybenzimidazole Blanket Insulation, or PBI, with a multi-use temperature limit of 1000F and above, in contrast to Shuttle FRSI which has a multi-use temperature limit of about 700F. 1.6 NASA Langley and BF Goodrich (formerly Rohr Corp.) have led the development of metallic-based TPS. This activity uses essentially three approaches: metallic tiles which encase a fibrous ceramic batting in a box fabricated largely from metallic honeycombs, typically Nickel based alloys; metallic honeycomb sheets, made of Nickel-based alloys, incorporating a fibrous back-side insulation encapsulated in a metallic foil bag, providing reduced weight; and metallic multi-wall, which is comprised of dimpled Titanium metal sheets, which are stacked and then diffusion bonded at contact points to form the TPS. The Nickel-based systems can be used up to temperatures of about 1800F, and the Titanium system to about 100F. These thirteen TPS materials have pros and cons to their usage in terms of temperature capability, weight, initial cost, and maintenance. Carbon-carbon and C/SiC systems have the highest temperature capability but are relatively expensive and heavy, requiring significant time, expertise, and costly facilities and tools for design and fabrication. Second generation ceramic tiles are relatively light, durable, simple to fabricate and easy to install; however, waterproofing is a concern. Blankets and felts are light, simple, inexpensive, and easy to install over curved vehicle surfaces, but durability and waterproofing are concerns. Metallics are robust and appear to have eliminated waterproofing, but they tend to be heavy and relatively expensive, requiring costly facilities and tools. If thin metal sheets are used to reduce weights, then issues arise from possible metal fatigue and corrosion caused by thermal cycling, pressure oscillations, and environmental exposure. For application to future RLVs, system analyses show that a significant component of the vehicle life cycle cost is from the TPS; however, it is difficult to quantify and to compare the potential savings of advanced systems without performing full vehicle designs using each of the different options. Because this entails a considerable effort and also tends to submerge TPS cost impacts under unrelated vehicle design assumptions, there is a clear need for a simpler quantitative method to evaluate the cost impact of different TPS options. To this end, this work introduces a TPS life-cycle cost parameter which is easily computed and applicable to generic RLVs.

Rasky, Dan J.↗

Multigrid approaches to non-linear diffusion problems on unstructured meshes

The efficiency of three multigrid methods for solving highly non-linear diffusion problems on two-dimensional unstructured meshes is examined. The three multigrid methods differ mainly in the manner in which the nonlinearities of the governing equations are handled. These comprise a non-linear full approximation storage (FAS) multigrid method which is used to solve the non-linear equations directly, a linear multigrid method which is used to solve the linear system arising from a Newton linearization of the non-linear system, and a hybrid scheme which is based on a non-linear FAS multigrid scheme, but employs a linear solver on each level as a smoother. Results indicate that all methods are equally effective at converging the non-linear residual in a given number of grid sweeps, but that the linear solver is more efficient in cpu time due to the lower cost of linear versus non-linear grid sweeps.

Mavriplis, Dimitri J.↗

Cloud Structure and Cloud-Radiation Interaction: From Smallest Observable Scales (Less than 10 m) to GCM Grid Scales (More than 100 km)

In this presentation we review the fractal nature of internal cloud structure from cm- to km-scales as captured by in-situ probes during long horizontal penetrations by aircraft. We uncover the non-Poissonian spatial distribution of cloud droplets at submeter scales and confirm scale-invariant behavior for large scales. Based on these structural characteristics, we generate simple fractal cloud models that reproduce statistical scaling properties of real clouds. These stochastic models represent a link between nonlinear science, in general, and cloud-radiation interaction, in particular. Next we run three-dimensional radiative transfer computations on these synthetic fractal clouds and compare the structure of the resulting radiation fields with the known structure of the cloud model and with satellite images of real clouds. The different behaviors observed for small and large-scale variabilities will be discussed in detail. We find that while the large-scale fluctuations of the resulting radiation fields resemble those in the original scale-invariant cloud structure, the radiation at small scales is much smoother than its cloud liquid water counterpart. This violates scale-invariance and produces a scale-break at 0.2-0.5 km that is clearly observed in high-resolution satellite data such as from Landsat. Finally, we show how radiative transfer Green function theory in the photon diffusion limit explains (and predicts) the above phenomena of "radiative smoothing."

Marshak, Alexander↗

Navigation Algorithms for the SeaWiFS Mission

The navigation algorithms for the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) were designed to meet the requirement of 1-pixel accuracy-a standard deviation (sigma) of 2. The objective has been to extract the best possible accuracy from the spacecraft telemetry and avoid the need for costly manual renavigation or geometric rectification. The requirement is addressed by postprocessing of both the Global Positioning System (GPS) receiver and Attitude Control System (ACS) data in the spacecraft telemetry stream. The navigation algorithms described are separated into four areas: orbit processing, attitude sensor processing, attitude determination, and final navigation processing. There has been substantial modification during the mission of the attitude determination and attitude sensor processing algorithms. For the former, the basic approach was completely changed during the first year of the mission, from a single-frame deterministic method to a Kalman smoother. This was done for several reasons: a) to improve the overall accuracy of the attitude determination, particularly near the sub-solar point; b) to reduce discontinuities; c) to support the single-ACS-string spacecraft operation that was started after the first mission year, which causes gaps in attitude sensor coverage; and d) to handle data quality problems (which became evident after launch) in the direct-broadcast data. The changes to the attitude sensor processing algorithms primarily involved the development of a model for the Earth horizon height, also needed for single-string operation; the incorporation of improved sensor calibration data; and improved data quality checking and smoothing to handle the data quality issues. The attitude sensor alignments have also been revised multiple times, generally in conjunction with the other changes. The orbit and final navigation processing algorithms have remained largely unchanged during the mission, aside from refinements to data quality checking. Although further improvements are certainly possible, future evolution of the algorithms is expected to be limited to refinements of the methods presented here, and no substantial changes are anticipated.

Hooker, Stanford B.↗

Progress in the Simulation of Steady and Time-Dependent Flows with 3D Parallel Unstructured Cartesian Methods

The proposed paper will present recent extensions in the development of an efficient Euler solver for adaptively-refined Cartesian meshes with embedded boundaries. The paper will focus on extensions of the basic method to include solution adaptation, time-dependent flow simulation, and arbitrary rigid domain motion. The parallel multilevel method makes use of on-the-fly parallel domain decomposition to achieve extremely good scalability on large numbers of processors, and is coupled with an automatic coarse mesh generation algorithm for efficient processing by a multigrid smoother. Numerical results are presented demonstrating parallel speed-ups of up to 435 on 512 processors. Solution-based adaptation may be keyed off truncation error estimates using tau-extrapolation or a variety of feature detection based refinement parameters. The multigrid method is extended to for time-dependent flows through the use of a dual-time approach. The extension to rigid domain motion uses an Arbitrary Lagrangian-Eulerlarian (ALE) formulation, and results will be presented for a variety of two- and three-dimensional example problems with both simple and complex geometry.

Aftosmis, M. J.↗

Use of TRMM Data to Test an Improved Parameterization of Stratiform Precipitation

We successfully implemented our EAUcup parameterization of convection in our most recent geodesic-grid version of the Colorado State University General Circulation Model (CSU GCM). Tendencies due to convection and large-scale cloud microphysics are computed using the same time-step, allowing smoother interactions between convective and large-scale moist processes. Figures illustrate our most recent progress on the simulation of convective and stratiform precipitations for January. Comparisons against TRMM precipitation radar (PR) data would reveal a systematic overestimation of the simulated monthly-mean convective rainfall over well-known areas of deep convection. One major deficiency in the global distribution of stratiform precipitation simulated in the CSU GCM is the lack of stratiform precipitation over land in the tropics. This result strongly disagree with Fowler and Randall who show that convective and stratiform precipitations always coexist and that their ratio to the total precipitation is about 0.5. We believe that stratiform precipitation over land evaporates before reaching the ground because the lower troposphere remains too dry in response to the lack of convective detrainment occurring immediately above the Planetary Boundary Layer (PBL). This issue is further investigated at present. We made significant progress in the development of our parameterization of fractional cloudiness. As an example of our ongoing research effort, a figure shows the global distribution of the cloud ice mixing ratio and corresponding fractional area of ice clouds at 250 hPa. These results were obtained using our latitude-longitude grid version of the CSU GCM. There is a nice correspondence between areas of high (low) values of the cloud ice mixing ratio and the cloud fraction. A nice correspondence between the global distribution of the cloud water mixing ratio and fraction of water clouds in the lower troposphere can also be found. We are pursuing this effort, but in the geodesic version of the GCM. We hope to be able to make comparison of the simulated precipitation against TRMM precipitation data in the near future.

Randall, David A.↗

A Powerful New Engine

Through Small Business Innovation Research (SBIR) funding from NASA's Glenn Research Center, Moller International created a new coating for rotary engines, which significantly improves the fuel consumption of a vehicle while reducing emissions. The new coatings are offered in the new Rotapower(R) engine, which is produced and distributed by Moller subsidiary, Freedom Motors, Inc. The coating allows the Rotapower engine to function smoother than other models, reducing wear and protecting the engine. The Rotapower engine has the ability to operate on a variety of fuels, including gasoline, natural gas, diesel, alcohol, and kerosene. A small and lightweight engine, it is projected to replace many of today's bulkier versions. The 10 horsepower model fits in the palm of one's hand, while the 160 horsepower model fits into a 5-gallon bucket. The clean running Rotapower engine is environmentally appealing, because it eliminates over 98 percent of the total emissions given off by traditional piston engines. Fewer pollutants are spewed into the air, making it especially attractive in areas where air pollution is a major problem. Due to the clean-burning nature of the engine, it meets the stringent standards set by the California Air Resources Board. The engine also has numerous commercial benefits in several types of recreational, industrial, and transportation applications, including personal watercraft, snowmobiles, portable generators. and pumps.

Source record↗

The Estimation Theory Framework of Data Assimilation

Lecture 1. The Estimation Theory Framework of Data Assimilation: 1. The basic framework: dynamical and observation models; 2. Assumptions and approximations; 3. The filtering, smoothing, and prediction problems; 4. Discrete Kalman filter and smoother algorithms; and 5. Example: A retrospective data assimilation system

Cohn, S.↗

Impact of Climatic Variability on Atmospheric Mass Distribution and GRACE-Derived Gravity Fields

During the period we calculated the atmospheric data sets related to its mass and angular momentum distribution. For mass, we determined the various harmonics from the NCEP-NCAR reanalysis, especially the low-order harmonics that are useful in studying the gravitation distribution as will be determined from the GRACE mission. Atmospheric mass is also related to the atmospheric loading on the solid Earth; we cooperated with scientists who needed the atmospheric mass information for understanding its contributions to the overall loading, necessary for vertical and horizontal coordinate estimation. We calculated atmospheric angular momentum from the NCEP-NCAR reanalyses and 4 operational meteorological centers, based on the motion (wind) terms and the mass (surface pressure) terms. These are associated with motions of the planet, including its axial component causing changes in the length of day, more related to the winds, and the equatorial component related to motions of the pole, more related to the mass. Tasks related to the ocean mass and angular momentum were added to the project as well. For these we have noted the ocean impact on motions of the pole as well as the torque mechanisms that relate the transfer of angular momentum between oceans and solid earth. The activities of the project may be summarized in the following first manuscript written in December 2002, for a symposium that Dr. Salstein attended on Geodynamics. We have continued to assess ocean angular momentum (OAM) quantities derived from bottom pressure and velocity fields estimated with our finite-difference barotropic (single layer) model. Three years of output (1993-95) from a run without any data constraints was compared to output from a corresponding run that was constrained by altimeter data using a Kalman filter and smoother scheme. Respective OAM time series were combined with corresponding atmospheric series and compared to observed polar motion. The constrained OAM series provided slightly better variance reduction than the unconstrained series. Analysis provided a check on the estimation scheme and pointed to further work to improve the determination of OAM using this method. A significant effort was also devoted to quantifying effects of uncertainties in high frequency winds on the mean and seasonal momentum exchange between atmosphere and oceans.

Salstein, David A.↗

Production of Near-Mirror Surface Quality by Precision Grinding

Mechanical components such as gears and bearings operate with the working surfaces in intimate contact with a mating part. The performance of such components will be influenced by the quality of the working surface. In general, a smoother surface will perform better than a rougher surface since the lubrication conditions are improved. For example, surfaces with a special near-mirror quality finish of low roughness performed better than ground surfaces when tested using a block-on-ring arrangement. Bearings with near-mirror quality have been tested and analyzed; lower running torques were measured and improved fatigue life was anticipated. Experiments have been done to evaluate the performance of gears with improved, low roughness surface finishing. The measured performance improvements include an increased scuffing (scoring) load capacity by a factor of 1.6, a 30-percent reduction of gear tooth running friction, and longer fatigue lives by a factor of about four. One can also anticipate that near-mirror quality surface finishing could improve the performance of other mechanical components such as mechanical seals and heavily loaded journal bearings. Given these demonstrated benefits, capable and economical methods for the production of mechanical components with near-mirror quality surfaces are desired. One could propose the production of near-mirror quality surfaces by several methods such as abrasive polishing, chemical assisted polishing, or grinding. Production of the surfaces by grinding offers the possibility to control the macro-geometry (form), waviness, and surface texture with one process. The present study was carried out to investigate the possibility of producing near-mirror quality surfaces by grinding. The present study makes use of a specially designed grinding machine spindle to improve the surface quality relative to the quality produced when using a spindle of conventional design.

Dimofte, Florin↗

The NASA/Ames Mars General Circulation Model: Model Improvements and Comparison with Observations

For many years, the NASA/Ames Mars General Circulation Model (GCM) has been built around the UCLA B-grid dynamical core. An attached tracer transport scheme based on the aerosol microphysical model of Toon et al. (1988) provided a tool for studying dust storm transport and feedbacks (Murphy et al., 1995). While we still use a B-grid version of the model, the Ames group is now transitioning to the ARIES/GEOS Goddard C-grid dynamical core (Suarez and Takacs, 1995). The C-grid produces smoother fields when the model top is raised above 50 km, and has a built in transport scheme for an arbitrary number of tracers. All of our transport simulations are now carried out with the C-grid. We have also been updating our physics package. Several years ago we replaced our bulk boundary layer scheme with a level 2 type diffusive scheme, and added a multi-level soil model (Haberle et al., 2000). More recently we replaced our radiation code with a more generalized two-stream code that accounts for aerosol multiple scattering and gaseous absorption. This code gives us much more flexibility in choosing aerosol optical properties and radiatively active gases.

Haberle, R. M.↗

Multigrid Solution of the Navier-Stokes Equations at Low Speeds with Large Temperature Variations

Multigrid methods for the Navier-Stokes equations at low speeds and large temperature variations are investigated. The compressible equations with time-derivative preconditioning and preconditioned flux-difference splitting of the inviscid terms are used. Three implicit smoothers have been incorporated into a common multigrid procedure. Both full coarsening and semi-coarsening with directional fine-grid defect correction have been studied. The resulting methods have been tested on four 2D laminar problems over a range of Reynolds numbers on both uniform and highly stretched grids. Two of the three methods show efficient and robust performance over the entire range of conditions. In addition none of the methods have any difficulty with the large temperature variations.

Sockol, Peter M.↗

NASA's Autonomous Formation Flying Technology Demonstration, Earth Observing-1(EO-1)

NASA's first autonomous formation flying mission, the New Millennium Program's (NMP) Earth Observing-1 (EO-1) spacecraft, recently completed its principal goal of demonstrating advanced formation control technology. This paper provides an overview of the evolution of an onboard system that was developed originally as a ground mission planning and operations tool. We discuss the Goddard Space Flight Center s formation flying algorithm, the onboard flight design and its implementation, the interface and functionality of the onboard system, and the implementation of a Kalman filter based GPS data smoother. A number of safeguards that allow the incremental phasing in of autonomy and alleviate the potential for mission-impacting anomalies from the on- board autonomous system are discussed. A comparison of the maneuvers planned onboard using the EO-1 autonomous control system to those from the operational ground-based maneuver planning system is presented to quantify our success. The maneuvers discussed encompass reactionary and routine formation maintenance. Definitive orbital data is presented that verifies all formation flying requirements.

Folta, David↗

Applications of Space-Filling-Curves to Cartesian Methods for CFD

The proposed paper presents a variety novel uses of Space-Filling-Curves (SFCs) for Cartesian mesh methods in 0. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, most are applicable on general body-fitted meshes -both structured and unstructured. We demonstrate the use of single O(N log N) SFC-based reordering to produce single-pass (O(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations. Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 512 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 10% of ideal even with only around 50,000 cells in each subdomain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with O(max(M,N)) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for finite-difference-based gradient design methods.

Aftosmis, Michael J.↗