Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fixed grid”

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 55 records · Page 3

Image Navigation and Registration (INR) Performance Assessment Tool Set (IPATS) for the GOES-R Advanced Baseline Imager and Geostationary Lightning Mapper

The GOES-R Flight Project has developed an Image Navigation and Registration (INR) Performance Assessment Tool Set (IPATS) for measuring Advanced Baseline Imager (ABI) and Geostationary Lightning Mapper (GLM) INR performance metrics in the post-launch period for performance evaluation and long term monitoring. For ABI, these metrics are the 3-sigma errors in navigation (NAV), channel-to-channel registration (CCR), frame-to-frame registration (FFR), swath-to-swath registration (SSR), and within frame registration (WIFR) for the Level 1B image products. For GLM, the single metric of interest is the 3-sigma error in the navigation of background images (GLM NAV) used by the system to navigate lightning strikes. 3-sigma errors are estimates of the 99.73rd percentile of the errors accumulated over a 24 hour data collection period. IPATS utilizes a modular algorithmic design to allow user selection of data processing sequences optimized for generation of each INR metric. This novel modular approach minimizes duplication of common processing elements, thereby maximizing code efficiency and speed. Fast processing is essential given the large number of sub-image registrations required to generate INR metrics for the many images produced over a 24 hour evaluation period. Another aspect of the IPATS design that vastly reduces execution time is the off-line propagation of Landsat based truth images to the fixed grid coordinates system for each of the three GOES-R satellite locations, operational East and West and initial checkout locations. This paper describes the algorithmic design and implementation of IPATS and provides preliminary test results.

Image Registration↗

Assessment of the SMAP Passive Soil Moisture Product

The National Aeronautics and Space Administration (NASA) Soil Moisture Active Passive (SMAP) satellite mission was launched on January 31, 2015. The observatory was developed to provide global mapping of high-resolution soil moisture and freeze-thaw state every two to three days using an L-band (active) radar and an L-band (passive) radiometer. After an irrecoverable hardware failure of the radar on July 7, 2015, the radiometer-only soil moisture product became the only operational Level 2 soil moisture product for SMAP. The product provides soil moisture estimates posted on a 36 kilometer Earth-fixed grid produced using brightness temperature observations from descending passes. Within months after the commissioning of the SMAP radiometer, the product was assessed to have attained preliminary (beta) science quality, and data were released to the public for evaluation in September 2015. The product is available from the NASA Distributed Active Archive Center at the National Snow and Ice Data Center. This paper provides a summary of the Level 2 Passive Soil Moisture Product (L2_SM_P) and its validation against in situ ground measurements collected from different data sources. Initial in situ comparisons conducted between March 31, 2015 and October 26, 2015, at a limited number of core validation sites (CVSs) and several hundred sparse network points, indicate that the V-pol Single Channel Algorithm (SCA-V) currently delivers the best performance among algorithms considered for L2_SM_P, based on several metrics. The accuracy of the soil moisture retrievals averaged over the CVSs was 0.038 cubic meter per cubic meter unbiased root-mean-square difference (ubRMSD), which approaches the SMAP mission requirement of 0.040 cubic meter per cubic meter.

passive↗

Geometry Modeling for Unstructured Mesh Adaptation

The quantification and control of discretization error is critical to obtaining reliable simulation results. Adaptive mesh techniques have the potential to automate discretization error control, but have made limited impact on production analysis workflow. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic mesh adaptation mechanics. However, the poor integration of initial mesh generation and adaptive mesh mechanics to typical sources of geometry has hindered adoption of adaptive mesh techniques, where these geometries are often created in Mechanical Computer- Aided Design (MCAD) systems. The difficulty of this coupling is compounded by two factors: the inherent complexity of the model (e.g., large range of scales, bodies in proximity, details not required for analysis) and unintended geometry construction artifacts (e.g., translation, uneven parameterization, degeneracy, self-intersection, sliver faces, gaps, large tolerances be- tween topological elements, local high curvature to enforce continuity). Manual preparation of geometry is commonly employed to enable fixed-grid and adaptive-grid workflows by reducing the severity and negative impacts of these construction artifacts, but manual process interaction inhibits workflow automation. Techniques to permit the use of complex geometry models and reduce the impact of geometry construction artifacts on unstructured grid workflows are models from the AIAA Sonic Boom and High Lift Prediction are shown to demonstrate the utility of the current approach.

Park, Michael A.↗

TPSAS-NF1676L-10980-DND

We plan to perform the following sets of computations on unadapted (fixed) grids: 1) Structured RANS set 1 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), 2) Structured RANS set 2 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 3) Structured RANS set 3 (time permitting) (Code: CFL3D, Grid: Str-OnetoOne-B-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 4) Unstructured RANS set 1 (Code: FUN3D, Grid: Unst-Mixed-FromTet-Nodecentered-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), and 5) Unstructured RANS set 2 (time permitting) (Code: FUN3D, Grid: Unst-Hex-FromOnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras),. Optional case 3 is not being computed. CFL3D is a structured upwind-biased cell-centered RANS code,1 and FUN3D is an unstructured upwind-biased node-centered RANS code

Elizabeth M Lee-Rausch↗

TPSAS-NF1676L-10454-DND

We plan to perform the following sets of computations on unadapted (fixed) grids: 1) Structured RANS set 1 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), 2) Structured RANS set 2 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 3) Structured RANS set 3 (time permitting) (Code: CFL3D, Grid: Str-OnetoOne-B-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 4) Unstructured RANS set 1 (Code: FUN3D, Grid: Unst-Mixed-FromTet-Nodecentered-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), and 5) Unstructured RANS set 2 (time permitting) (Code: FUN3D, Grid: Unst-Hex-FromOnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras),. Optional case 3 is not being computed. CFL3D is a structured upwind-biased cell-centered RANS code,1 and FUN3D is an unstructured upwind-biased node-centered RANS code

C L Rumsey↗

HLPW-4/GMGW-3: Overview and Workshop Summary

The Fourth AIAA CFD High Lift Prediction Workshop and the Third Geometry and Mesh Generation Workshop were held collaboratively with the common goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this joint endeavor was the use of Technology Focus Groups, an innovative new approach for workshops involving close collaboration between participants. These groups, which included both mesh generation and flow solver experts, worked to accelerate advancements for their particular methodologies by addressing key questions of importance {\em prior} to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop. Measured experimental wind tunnel data were available for comparison. The workshop also included a two-dimensional turbulence model verification exercise based on the CRM-HL wing shape. Altogether, 44 participants submitted a total of 184 data sets of CFD results. This paper provides a high-level summary of the results and conclusions from the workshop. Like at past workshops, fixed-grid Reynolds-averaged Navier-Stokes continued to be inaccurate and inconsistent for high lift. However, mesh adaptation definitively brought more consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics.

Christopher L Rumsey↗

HLPW-5: Overview and Workshop Summary

The Fifth AIAA CFD High-Lift Prediction Workshop was held with the goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this endeavor was the use of Technology Focus Groups (TFG), which included both mesh generation and flow-solver experts working together to accelerate advancements for their particular CFD methodology, by addressing key questions of importance prior to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop, and was used for three unique test cases. Wind-tunnel data were available for comparison for one of the test cases. Altogether, 365 datasets of CFD results were submitted by 47 teams, with 41 teams contributing to the multiple configurations of Case 1, 40 to Case 2, and 18 to Case 3. This paper provides a high-level summary of the results and conclusions from the workshop. As concluded from past workshops, fixed-grid Reynolds-averaged Navier-Stokes methods continued to be inaccurate and inconsistent for high-lift flows near maximum lift. However, application of mesh-adaptation technology helped to achieve improved consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics, particularly at maximum lift. Best practices for these methods were refined over the course of the workshop and new challenges were identified.

Adam M Clark↗

High-Lift Prediction Workshop 5: Summary of Reynolds-Averaged Navier–Stokes Technology Focus Group Summary

The fifth High-Lift Prediction Workshop (HLPW-5), which involved the high-lift version of the NASA common research model in several configurations, assessed various computational fluid dynamics methods, including Reynolds-averaged Navier-Stokes (RANS) and hybrid large-eddy simulations. This paper summarizes RANS solutions computed on fixed grids. Case 1, a verification case, considered a simple wing-body configuration and focused on grid convergence of lift, drag, and pitching moment coefficients. Case 2 is a configuration buildup case that focused on predicting the effects of increasing geometric complexity. For buildup configurations with slats, flaps, and a nacelle/pylon (configurations 2.2, 2.3, and 2.4), wind-tunnel data were provided by ONERA. Case 3 focused on the reference landing configuration at four Reynolds-number conditions. For Case 1, grid-converged RANS solutions were achieved using the standard Spalart-Allmaras (SA) turbulence model and the SA model with a quadratic constitutive relation and a rotation correction. Agreement between RANS solutions was observed for the simplest configuration 2.1 of Case 2 with the standard SA model. For other configurations, agreement between RANS solutions was hampered by insufficient iterative and grid convergence, especially at high angles of attack. In comparison with the experiment, RANS solutions qualitatively showed the correct configuration buildup trend but underpredicted lift and overpredicted both drag and pitching moment at high angles of attack.

Boris Diskin↗

A ground track control algorithm for the Topographic Mapping Laser Altimeter (TMLA)

The results of an analysis of an algorithm that will provide autonomous onboard orbit control using orbits determined with Global Positioning System (GPS) data. The algorithm uses the GPS data to (1) compute the ground track error relative to a fixed longitude grid, and (2) determine the altitude adjustment required to correct the longitude error. A program was written on a personal computer (PC) to test the concept for numerous altitudes and values of solar flux using a simplified orbit model including only the J sub 2 zonal harmonic and simple orbit decay computations. The algorithm was then implemented in a precision orbit propagation program having a full range of perturbations. The analysis showed that, even with all perturbations (including actual time histories of solar flux variation), the algorithm could effectively control the spacecraft ground track and yield more than 99 percent Earth coverage in the time required to complete one coverage cycle on the fixed grid (220 to 230 days depending on altitude and overlap allowance).

Blaes, V.↗

Aerodynamic and Aerothermal Simulations of Mars Concept Vehicles Using Overset DPLR

NASA's Descent System Studies (DSS) Program is studying various concept vehicles to enable landing of heavy payloads on the surface of Mars. While it is desirable to run high-fidelity CFD simulations to accurately assess the aerodynamic and aerothermal effects of various design changes during EDL, it is usually difficult to quickly generate high-quality grids suitable for such analyses. One approach to address this bottleneck in mesh generation is through the use oversetting grids. Although the overset approach is efficient and powerful in solving partial differential equations on complex geometries, new users often find it challenging to apply overset concepts for their simulations. For example, generating hyperbolic grids with sufficient overlap; priority in hole-cutting on multiple overlapping grids; and fixes to assemble overlapping viscous grids at the body surface. The objective of this presentation is to introduce a simple process that combines the advantages of near-body, point-matched, structured grids with oversetting background grids suitable for grid alignment. This approach allows for grids that can be sequenced, reclustering of mesh spacing at the wall, and grid alignment with the bow shock. The current methodology is tested on a Mid-L/D configuration using the overset DPLR code.

Overset Grids↗

Efficient Algorithm for Rectangular Spiral Search

An algorithm generates grid coordinates for a computationally efficient spiral search pattern covering an uncertain rectangular area spanned by a coordinate grid. The algorithm does not require that the grid be fixed; the algorithm can search indefinitely, expanding the grid and spiral, as needed, until the target of the search is found. The algorithm also does not require memory of coordinates of previous points on the spiral to generate the current point on the spiral.

Brugarolas, Paul↗

A Survey of Composite Grid Generation for General Three-dimensional Sections

The generation and use of composite grids for general three-dimensional physical boundary configurations is discussed, and the availability of several codes or procedures is noted. With the composite framework, the physical region is segmented into sub-regions, each bounded by six curved sides, and a grid is generated in each sub-region. These grids may be joined at the interfaces between the sub-regions with various degrees of continuity. This structure allows codes to be constructed to operate on rectangular blocks in computational space, so that existing solution procedures can be readily incorporated in the construction of codes for general configurations. Numerical grid generation is an integral part of the numerical solution of partial differential equations and is one of the pacing items in the development of codes for general configurations. The numerically generated grid frees the computational simulation from restriction to certain boundary shapes and allows general codes to be written in which the boundary shape is specified simply by input. The numerically generated grid allows all computation to be done on a fixed square grid in the computational space, which is always rectangular by construction.

Thompson, J. F.↗

Ionospheric calibration for single frequency altimeter measurements

This report investigates the potential of using Global Positioning System (GPS) data and a model of the ionosphere to supply a measure of the sub-satellite Total Electron Current (TEC) of the required accuracy (10 TECU rms) for the purpose of calibrating single frequency radar altimeter measurements. Since climatological (monthly mean) models are known to be in error by as much as 50 percent, this work focused on the Parameterized Real-Time Ionospheric Specification Model (PRISM) which has the capability to improve model accuracy by ingesting (adjusting to) in situ ionospheric measurements. A set of globally distributed TEC measurements were generated using GPS data and were used as input to improve the accuracy of the PRISM model. The adjusted PRISM TEC values were compared to TOPEX dual frequency TEC measurements (which are considered truth) for a number of TOPEX sub-satellite tracks. The adjusted PRISM values generally compared to the TOPEX measurements within the 10 TECU accuracy requirements when the sub-satellite track passed within 300 to 400 km of the GPS TEC data or when the track passed through a night time ionosphere. However, when the sub-satellite points were greater than 300 to 400 km away from the GPS TEC data or when a local noon ionosphere was sampled, the adjusted PRISM values generally differed by greater than 10 TECU rms with data excursions from the TOPEX TEC measurements of as much as 40 TECU (an 8 cm path delay error at K band). Therefore, it can be concluded from this analysis that an unrealistically large number of GPS stations would be needed to predict sub-satellite TEC at the 10 TECU level in the day time ionosphere using a model such as PRISM. However, a technique currently being studied at the Jet Propulsion Laboratory (JPL) may provide a means of supplying adequate TEC data to meet the 10 TECU ionospheric correction accuracy when using a realistic number of ionospheric stations. This method involves using global GPS TEC data to estimate a global grid of vertical ionospheric TEC as a function of time (i.e. every one half hour) in a sun-fixed longitude frame. Working in a sun-fixed longitude frame, one is not limited by the spatial decorrelation distance of the ionosphere, but instead is limited more by the temporal correlations of the ionosphere in the sun-fixed frame which are a smaller effect. It is the opinion of the authors that using the global sun-fixed TEC grid data, in particular, ingesting it into PRISM, offers the best possibility of meeting the 10 TECU ionospheric correction accuracy requirement, and should be the subject of further study.

Schreiner, William S.↗

An Assessment of the Differences Between Spatial Resolution and Grid Size for the SMAP Enhanced Soil Moisture Product over Homogeneous Sites

Satellite-based passive microwave remote sensing typically involves a scanning antenna that makes measurements at irregularly spaced locations. These locations can change on a day to day basis. Soil moisture products derived from satellite-based passive microwave remote sensing are usually resampled to a fixed Earth grid that facilitates their use in applications. In many cases the grid size is finer than the actual spatial resolution of the observation, and often this difference is not well understood by the user. Here, this issue was examined for the Soil Moisture Active Passive (SMAP) enhanced version of the passive-based soil moisture product, which has a grid size of 9-km and a nominal spatial resolution of 33-km. In situ observations from core validation sites were used to compute comparison metrics. For sites that satisfied the established reliability and scaling criteria, the impact of validating the 9-km grid product with in situ data collected over a 9-km versus a 33-km domain was very small for the sites studied (0.039 cu. m/cu. m unbiased root mean square difference for the 9-km case versus 0.037 cu. m/cu. m for the 33-km case). This result does not mean that the resolution of the product is 9-km but that for the conditions studied here the soil moisture estimated from in situ observations over 9-km is a close approximation of the soil moisture estimated from in situ observations over the 33-km resolution. The implication is that using the enhanced SMAP product at its grid resolution of 9-km should not introduce large errors in most applications.

Soil moisture↗

Streamwise Upwind, Moving-Grid Flow Algorithm

Extension to moving grids enables computation of transonic flows about moving bodies. Algorithm computes unsteady transonic flow on basis of nondimensionalized thin-layer Navier-Stokes equations in conservation-law form. Solves equations by use of computational grid based on curvilinear coordinates conforming to, and moving with, surface(s) of solid body or bodies in flow field. Simulates such complicated phenomena as transonic flow (including shock waves) about oscillating wing. Algorithm developed by extending prior streamwise upwind algorithm solving equations on fixed curvilinear grid described in "Streamwise Algorithm for Simulation of Flow" (ARC-12718).

Goorjian, Peter M.↗

Sensitivity of Spaceborne and Ground Radar Comparison Results to Data Analysis Methods and Constraints

With the availability of active weather radar observations from space from the Precipitation Radar (PR) on board the Tropical Rainfall Measuring Mission (TR.MM) satellite, numerous studies have been performed comparing PR reflectivity and derived rain rates to similar observations from ground-based weather radars (GR). These studies have used a variety of algorithms to compute matching PR and GR volumes for comparison. Most studies have used a fixed 3-dimensional Cartesian grid centered on the ground radar, onto which the PR and GR data are interpolated using a proprietary approach and/or commonly available GR analysis software (e.g., SPRINT, REORDER). Other studies have focused on the intersection of the PR and GR viewing geometries either explicitly or using a hybrid of the fixed grid and PR/GR common fields of view. For the Dual-Frequency Precipitation Radar (DPR) of the upcoming Global Precipitation Measurement (GPM) mission, a prototype DPR/GR comparison algorithm based on similar TRMM PR data has been developed that defines the common volumes in terms of the geometric intersection of PR and GR rays, where smoothing of the PR and GR data are minimized and no interpolation is performed. The PR and GR volume-averaged reflectivity values of each sample volume are accompanied by descriptive metadata, for attributes including the variability and maximum of the reflectivity within the sample volume, and the fraction of range gates in the sample average having reflectivity values above an adjustable detection threshold (typically taken to be 18 dBZ for the PR). Sample volumes are further characterized by rain type (Stratiform or Convective), proximity to the melting layer, underlying surface (land/water/mixed), and the time difference between the PR and GR observations. The mean reflectivity differences between the PR and GR can differ between data sets produced by the different analysis methods; and for the GPM prototype, by the type of constraints and categorization applied to the data. In this paper, we will show results comparing the 3-D gridded analysis "black box" approach to the GPM prototype volume-matching approach, using matching TRMM PR and WSR-88D ground radar data. The affects of applying data constraints and data categorizations on the volume-matched data to the results will be shown, and explanations of the differences in terms of data and analysis algorithm characteristics will be presented. Implications of the differences to the determination of PR/DPR calibration differences and use of ground radar data to evaluate the PR and DPR attenuation correction algorithms will be discussed.

Morris, Kenneth R.↗

Numerical solution of three-dimensional free turbulent shear flows

Numerical results are presented for three-dimensional compressible turbulent jet and wake flows. An alternating direction implicit numerical procedure is used to solve the finite-difference form of the parabolic elliptic Navier-Stokes equations. A coordinate transformation maps the boundaries at infinity into a finite computational domain in order to properly specify infinity boundary conditions as well as contain the downstream growth of the viscous flow field in a fixed computational grid. Turbulence closure is achieved through an algebraic mixing length eddy viscosity model. Numerical results for supersonic flow are presented for an axisymmetric jet, an elliptical jet, an elliptical wake, and two interacting rectangular jets. Experimental data were not available for comparison with the numerical results. However, the results compare well with empirical results for free shear flows.

Oh, Y. H.↗