Search NASA⌕ Search

SEARCH · Search NASA

Results for “global refinement”

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 19 records

Application of Parallel Adjoint-Based Error Estimation and Anisotropic Grid Adaptation for Three-Dimensional Aerospace Configurations

This paper demonstrates the extension of error estimation and adaptation methods to parallel computations enabling larger, more realistic aerospace applications and the quantification of discretization errors for complex 3-D solutions. Results were shown for an inviscid sonic-boom prediction about a double-cone configuration and a wing/body segmented leading edge (SLE) configuration where the output function of the adjoint was pressure integrated over a part of the cylinder in the near field. After multiple cycles of error estimation and surface/field adaptation, a significant improvement in the inviscid solution for the sonic boom signature of the double cone was observed. Although the double-cone adaptation was initiated from a very coarse mesh, the near-field pressure signature from the final adapted mesh compared very well with the wind-tunnel data which illustrates that the adjoint-based error estimation and adaptation process requires no a priori refinement of the mesh. Similarly, the near-field pressure signature for the SLE wing/body sonic boom configuration showed a significant improvement from the initial coarse mesh to the final adapted mesh in comparison with the wind tunnel results. Error estimation and field adaptation results were also presented for the viscous transonic drag prediction of the DLR-F6 wing/body configuration, and results were compared to a series of globally refined meshes. Two of these globally refined meshes were used as a starting point for the error estimation and field-adaptation process where the output function for the adjoint was the total drag. The field-adapted results showed an improvement in the prediction of the drag in comparison with the finest globally refined mesh and a reduction in the estimate of the remaining drag error. The adjoint-based adaptation parameter showed a need for increased resolution in the surface of the wing/body as well as a need for wake resolution downstream of the fuselage and wing trailing edge in order to achieve the requested drag tolerance. Although further adaptation was required to meet the requested tolerance, no further cycles were computed in order to avoid large discrepancies between the surface mesh spacing and the refined field spacing.

Lee-Rausch, E. M.↗

Empirical comparison of heuristic load distribution in point-to-point multicomputer networks

The study compared several load placement algorithms using instrumented programs and synthetic program models. Salient characteristics of these program traces (total computation time, total number of messages sent, and average message time) span two orders of magnitude. Load distribution algorithms determine the initial placement for processes, a precursor to the more general problem of load redistribution. It is found that desirable workload distribution strategies will place new processes globally, rather than locally, to spread processes rapidly, but that local information should be used to refine global placement.

Grunwald, Dirk C.↗

Global/local stress analysis of composite panels

A method for performing a global/local stress analysis is described, and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as boundary conditions for the local model. Then, the local model is analyzed independent of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Ransom, Jonathan B.↗

Global/local stress analysis of composite structures

A method for performing a global/local stress analysis is described and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as boundary conditions for the local model. Then, the local model is analyzed independent of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Ransom, Jonathan B.↗

Global/local stress analysis of composite panels

A method for performing a global/local stress analysis is described and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as 'boundary conditions' for the local model. Then, the local model is analyzed independently of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Ransom, J. B.↗

Dynamic Load Balancing for Adaptive Unstructured Grids

Dynamic mesh adaptation on unstructured grids is a powerful tool for computing unsteady three-dimensional problems that require grid modifications to efficiently resolve solution features. By locally refining and coarsening the mesh to capture phenomena of interest, such procedures make standard computational methods more cost effective. Highly refined meshes are required to accurately capture shock waves, contact discontinuities, vortices, and shear layers in fluid flow problems. Adaptive meshes have also proved to be useful in several other areas of computational science and engineering like computer vision and graphics, semiconductor device modeling, and structural mechanics. Local mesh adaptation provides the opportunity to obtain solutions that are comparable to those obtained on globally-refined grids but at a much lower cost. Additional information is contained in the original extended abstract.

Biswas, Rupak↗

TOPEX orbit determination and gravity recovery using Global Positioning System data from repeat orbits

A covariance analysis is presented for satellite tracking and gravity recovery with a differential Global Positioning System-based technique to be demonstrated on TOPEX in the early 1990s. The technique employs data from an ensemble of repeat ground tracks to recover a unique satellite epoch state for each track and a set of invariant positional parameters common to all tracks. The positional parameters represent the effect of mismodeled gravitational field on the satellite orbit. At an altitude of 1336 km, where gravity modeling is the dominant systematic error, averaging of random error over many arcs and adjustment of the gravity model reduce the final satellite position error. The positional parameters can then be used to produce a refined global gravity model. The analysis indicates that errors ranging from 5 to 8 cm in TOPEX altitude and 0.05 to 0.2 mGal for the gravity field can be achieved, depending on the number of repeat arcs used.

Wu, Jiun-Tsong↗

Multigrid solution of the Euler equations on unstructured and adaptive meshes

A multigrid algorithm has been developed for solving the steady-state Euler equations in two dimensions on unstructured triangular meshes. The method assumes the various coarse and fine grids of the multigrid sequence to be independent of one another, thus decoupling the grid generation procedure from the multigrid algorithm. The transfer of variables between the various meshes employs a tree-search algorithm which rapidly identifies regions of overlap between coarse and fine grid cells. Finer meshes are obtained either by regenerating new globally refined meshes, or by adaptively refining the previous coarser mesh. For both cases, the observed convergence rates are comparable to those obtained with structured multigrid Euler solvers. The adaptively generated meshes are shown to produce solutions of higher accuracy with fewer mesh points.

Mavriplis, Dimitri↗

Topex orbit determination and gravity recovery using GPS data from repeat orbits

A covariance analysis is presented for a differential GPS-based technique for satellite tracking and gravity recovery to be demonstrated on Topex in the early 1990's. The technique employs data from an ensemble of repeat ground tracks to recover a unique satellite epoch state for each track and a set of invariant local gravity parameters common to all tracks. For satellites above about 800 km altitude, where gravity modeling is the dominant systematic error, averaging of random error over many arcs and adjustment of the gravity model reduces the final satellite position error. The local gravity parameters can then be used to produce a refined global gravity model. The analysis indicates that errors ranging from 5 to 8 cm in Topex altitude and 0.05 to 0.2 mgal for the gravity field can be achieved, depending on the number of repeat arcs used.

Wu, Jiun-Tsong↗

Time-distance helioseismology

It is shown here that it is possible to extract time-distance information from temporal cross-correlations of the intensity fluctuation on the solar surface. This approach opens the way for seismic studies of local solar phenomena such subsurface inhomogeneities near sunspots and should help to refine global models of the internal velocity stratification in the sun.

Duvall, T. L., Jr.↗

Use of satellite ocean color observations to refine understanding of global geochemical cycles

In October 1978, the first satellite-borne color sensor, the Coastal Zone Color Scanner (CZCS), was launched aboard Nimbus-7 with four visible and two infrared bands, permitting a sensitivity about 60 times that of the Landsat-1 multispectral scanner. The CZCS radiance data can be utilized to estimate ocean chlorophyll concentrations by detecting shifts in sea color, particularly in oceanic waters. The obtained data can be used in studies regarding problems of overfishing, and, in addition, in investigations concerning the consequences of man's accelerated extraction of nitrogen from the atmosphere and addition of carbon to the atmosphere. The satellite data base is considered along with a simulation analysis, and ships providing ground-truth chlorophyll measurements in the ocean.

Walsh, J. J.↗

Grid Refinement Techniques for the 𝜸-Re𝜽𝒕 Transition Model in FUN3D

There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.

Grid Adaptation↗

Grid Refinement Techniques for the 𝜸-Re 𝜽 𝒕 Transition Model in FUN3D

There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.

Grid Adaptation↗

Antarctic Rebound and the Time-Dependence of the Earth's Shape

Great strides have been made during the past 30 years in refining models of the last global glaciation. The refinements draw upon a vastly expanded relative sea level and sedimentary core record. Furthermore, we now possess a sharpened understanding of the mechanisms that drive climate changes associated with deglaciation. Some 15 years ago, using only 5.5 years of ranging data, analyses of the drift in LAGEOS I node acceleration was used to infer that postglacial rebound was responsible for a secular change in the Earth's ellipsoidal shape (Yoder et al., .1983]. Today there exists a wealth of geodynamics satellite orbit data that constrain the secular time-dependence of the Earth's shape and low order gravity field, which includes mass redistribution from present-day glacier and great ice sheet imbalance and from postglacial rebound. We have shown that an unambiguous determination of the secular variation in the Earth's pear shaped harmonic (l = 3, m = 0) might provide information that bears on the present-day mass balance of Antarctica. This issue is revisited in light of new constraints on glacial loading during the late-Pleistocene and Holocene. An especially critical issue for the interpretation of secular odd degree zonal harmonics, l = 3 to 7, is the timing and magnitude of the deglaciation of Antarctica from Last Glacial Maximum. We explore ways in which the recovery of secular variation in both zonal and non-zonal harmonics for l = 2 through 7 can improve constraints on both rebound and present-day ice sheet balance.

Ivins, Erik R.↗

Proposal for the design of a zero gravity tool storage device

Astronauts frequently use a variety of hand tools during space missions, especially on repair missions. A toolbox is needed to allow storage and retrieval of tools with minimal difficulties. The toolbox must contain tools during launch, landing, and on-orbit operations. The toolbox will be used in the Shuttle Bay and therefore must withstand the hazardous space environment. The three main functions of the toolbox in space are: to protect the tools from the space environment and from damaging one another, to allow for quick, one-handed access to the tools; and to minimize the heat transfer between the astronaut's hand and the tools. This proposal explores the primary design issues associated with the design of the toolbox. Included are the customer and design specifications, global and refined function structures, possible solution principles, concept variants, and finally design recommendations.

Stuckwisch, Sue↗

[Analysis of Multiplatform CO (Carbon Monoxide) Measurements During Trace-P Mission]

Carbon monoxide is considered mission critical (TRACE-P NRA) because it is one of the gases involved in controlling the oxidizing power of the atmosphere and, as a tracer gas, is valuable in interpreting mission data sets. Carbon monoxide exhibits interannual differences, suggesting relatively short-term imbalances in sources and sinks. Sources of CO are dominated by fossil fuel combustion, biomass burning, and the photochemical oxidation of CH4 and nonmethane hydrocarbons while reaction with OH is believed to be the major sink for atmospheric CO, with additional losses due to soil uptake. Uncertainties in the magnitude and distribution of both sources and sinks remain fairly large however, and additional data are required to refine the global budget. Seasonal changes and a northern hemispheric latitudinal gradient have been described for a variety of Pacific basin sites through long-term monitoring of surface background levels. Latitudinal variations have also recently been described at upper tropospheric altitudes over a multi-year period by. TRACE-P will provide an aircraft survey of CO over the northern Pacific in the northern spring when CO concentrations are at their seasonal maximum in the northern hemisphere (NH) and at their seasonal minimum in the southern hemisphere (SH). Previous GTE missions, Le., PEM West-B and PEM Tropics-B, ground-based, and satellite observations (MAPS, April 1994) give us a general picture of the distribution of CO over the northern Pacific during this season. Based on these measurements, background CO levels over remote ocean areas are anticipated to be in the range of 110 - 180 ppbv, while those closer to the Asian continent may rise as high as 600 ppbv. These measurements also reveal high spatial variability (both horizontal and vertical) as well as temporal variations in CO over the area planned for the TRACE-P mission. This variability is a result of multiple CO sources, the meteorological complexity of transport processes, and the photochemical aging of air masses. The influence of biomass burning in the southern Pacific should be relatively small since the mission coincides with the southern tropical wet season when agricultural burning is at its seasonal low. The proposed CO measurements taken during TRACE-P should therefore largely be a function of the impact of various NH sources, primarily Asian and predominantly fossil fuel combustion and biomass burning. These processes are also major sources of many other atmospheric pollutants, consequently making accurate and precise CO measurements is one of the highest TRACE-P priorities [TRACE-P NRA]. The TRACE-P mission emphasizes the dual objectives of assessing the magnitude of the transport of chemically and radiatively important gases such as CO from Asia to the western Pacific, and determining how emissions change and are modified during this transport.

Pougatchev, Nikita S.↗

What it Takes to Compute Highly Accurate Rovibrational Line Lists for Use in Astrochemistry

We review the Best Theory + Reliable High-resolution Experiment (BTRHE) strategy for obtaining highly accurate molecular rovibrational line lists with InfraRed (IR) intensities. The need for highly accurate molecular rovibrational line lists is two-fold: a) assignment of the many rovibrational lines for common stable molecules especially those that exhibit a large amplitude motion, such as NH3, or have a high density of states such as SO2; and b) characterization of the atmospheres of exoplanets which will be one of the main areas of research in astronomy in the coming decades. The first motivation arises due to the need to eliminate lines due to common molecules in an astronomical observation in order to identify lines from new molecules, while the second motivation arises due to the need to obtain accurate molecular opacities in order to characterize the atmosphere of an exoplanet. The BTRHE strategy first consists of using high quality ab initio quantum chemical methods to obtain a global potential energy surface (PES) and dipole moment surface (DMS) that contains the proper physics. The global PES is then refined using a subset of the reliable high-resolution experimental data. The refined PES then gives energy level predictions to an accuracy similar to the reproduction accuracy of the experimental data used in the refinement step in the interpolation region (i.e., within the range of the experimental data used in the refinement step). The accuracy of the energy levels will slowly degrade as they are extrapolated to spectral regions beyond the high-resolution experimental data used in the refinement step. However, because the degradation is slow, the predicted energy levels can be used to assign new high-resolution experiments, and the data from these can then be used in a subsequent refinement step. In this way, the global PES eventually can yield highly accurate energy levels for all desired spectral regions including to very high energies and high J values. We show that IR intensities computed with the BTRHE rovibrational wavefunctions and the DMS can be very accurate provided one has minimized the fitting error of the DMS and tested the completeness of the DMS. Some examples of our work on NH3, CO2, and SO2 are given to highlight the usefulness of the BTRHE strategy and to provide ideas on how to further improve its predictive power in the future. In particular, it is shown how successive refinement steps, once new high-resolution data is available, can lead to PESs that yield highly accurate transition energies to larger spectral regions. The importance of including non-adiabatic corrections to reduce the J-dependence of errors for H-containing molecules is shown with work on NH3. Another very important aspect of the BTRHE approach is the consistency across isotopologues, which allows for highly accurate line lists for any isotopologue once one is obtained for the main isotopologue (which has more high-resolution data available for refinement).

Xinchuan Huang↗

Vortex-dominated conical-flow computations using unstructured adaptively-refined meshes

A conical Euler/Navier-Stokes algorithm is presented for the computation of vortex-dominated flows. The flow solver involves a multistage Runge-Kutta time stepping scheme which uses a finite-volume spatial discretization on an unstructured grid made up of triangles. The algorithm also employs an adaptive mesh refinement procedure which enriches the mesh locally to more accurately resolve the vortical flow features. Results are presented for several highly-swept delta wing and circular cone cases at high angles of attack and at supersonic freestream flow conditions. Accurate solutions were obtained more efficiently when adaptive mesh refinement was used in contrast with refining the grid globally. The paper presents descriptions of the conical Euler/Navier-Stokes flow solver and adaptive mesh refinement procedures along with results which demonstrate the capability.

Batina, John T.↗