Search NASA⌕ Search

SEARCH · Search NASA

Results for “Edge Based”

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 73 records · Page 4

The design and implementation of a parallel unstructured Euler solver using software primitives

This paper is concerned with the implementation of a 3D unstructured-grid Euler-solver on massively parallel distributed-memory computer architectures. The goal is to minimize solution time by achieving high computational rates with a numerically efficient algorithm. An unstructured multigrid algorithm with an edge-based data-structure has been adopted, and a number of optimizations have been devised and implemented in order to accelerate the parallel computational rates. The implementation is carried out by creating a set of software tools, which ease the implementation of computational problems on parallel architecture machines by relieving the user of the low-level machine specific issues. The quantitative effect of the various optimizations are demonstrated, and we show that the combined effect of these optimizations leads to roughly a factor of three performance improvement. The overall solution efficiency is compared with that obtained on the CRAY-YMP vector supercomputer.

Das, R.↗

Seeing ghosts - Photometry of Saturn's G Ring

Saturn's faint and narrow G Ring is only visible to the eye in two Voyager images, each taken at a rather high solar phase angle of about 160 deg. In this paper we introduce a new photometric technique for averaging across multiple Voyager images, and use it to detect the G Ring at several additional viewing geometries. The resultant phase curve suggests that the G Ring is composed of dust particles obeying a very steep power-law size distribution. The dust is generally smaller than that seen in other rings, ranging down to 0.03 micron. The G Ring occupies the region between orbital radii 166,000 and 173,000 km, and has a peak somewhat closer to the inner edge. Based on these limits, we demonstrate that Voyager 2 passed through and directly sampled this ring during its 1981 encounter with Saturn. Combined analysis of additional data sets suggests that a population of larger bodies is also present in the G Ring; these bodies occupy a narrower band near the observed peak and are likely the source for the visible dust. Based on some preliminary dynamical models, we propose that these larger bodies represent leftover debris from the collisional breakup of a small moon in Saturn's distant past.

Showalter, Mark R.↗

A three dimensional multigrid Reynolds-averaged Navier-Stokes solver for unstructured meshes

A three-dimensional unstructured mesh Reynolds averaged Navier-Stokes solver is described. Turbulence is simulated using a single field-equation model. Computational overheads are minimized through the use of a single edge-based data-structure, and efficient multigrid solution technique, and the use of multi-tasking on shared memory multi-processors. The accuracy and efficiency of the code are evaluated by computing two-dimensional flows in three dimensions and comparing with results from a previously validated two-dimensional code which employs the same solution algorithm. The feasibility of computing three-dimensional flows on grids of several million points in less than two hours of wall clock time is demonstrated.

Mavriplis, D. J.↗

Vector quantization for efficient coding of upper subbands

This paper examines the application of vector quantization (VQ) to exploit both intra-band and inter-band redundancy in subband coding. The focus here is on the exploitation of inter-band dependency. It is shown that VQ is particularly suitable and effective for coding the upper subbands. Three subband decomposition-based VQ coding schemes are proposed here to exploit the inter-band dependency by making full use of the extra flexibility of VQ approach over scalar quantization. A quadtree-based variable rate VQ (VRVQ) scheme which takes full advantage of the intra-band and inter-band redundancy is first proposed. Then, a more easily implementable alternative based on an efficient block-based edge estimation technique is employed to overcome the implementational barriers of the first scheme. Finally, a predictive VQ scheme formulated in the context of finite state VQ is proposed to further exploit the dependency among different subbands. A VRVQ scheme proposed elsewhere is extended to provide an efficient bit allocation procedure. Simulation results show that these three hybrid techniques have advantages, in terms of peak signal-to-noise ratio (PSNR) and complexity, over other existing subband-VQ approaches.

Zeng, W. J.↗

User's and test case manual for FEMATS

The FEMATS program incorporates first-order edge-based finite elements and vector absorbing boundary conditions into the scattered field formulation for computation of the scattering from three-dimensional geometries. The code has been validated extensively for a large class of geometries containing inhomogeneities and satisfying transition conditions. For geometries that are too large for the workstation environment, the FEMATS code has been optimized to run on various supercomputers. Currently, FEMATS has been configured to run on the HP 9000 workstation, vectorized for the Cray Y-MP, and parallelized to run on the Kendall Square Research (KSR) architecture and the Intel Paragon.

Chatterjee, Arindam↗

A unified multigrid solver for the Navier-Stokes equations on mixed element meshes

A unified multigrid solution technique is presented for solving the Euler and Reynolds-averaged Navier-Stokes equations on unstructured meshes using mixed elements consisting of triangles and quadrilaterals in two dimensions, and of hexahedra, pyramids, prisms, and tetrahedra in three dimensions. While the use of mixed elements is by no means a novel idea, the contribution of the paper lies in the formulation of a complete solution technique which can handle structured grids, block structured grids, and unstructured grids of tetrahedra or mixed elements without any modification. This is achieved by discretizing the full Navier-Stokes equations on tetrahedral elements, and the thin layer version of these equations on other types of elements, while using a single edge-based data-structure to construct the discretization over all element types. An agglomeration multigrid algorithm, which naturally handles meshes of any types of elements, is employed to accelerate convergence. An automatic algorithm which reduces the complexity of a given triangular or tetrahedral mesh by merging candidate triangular or tetrahedral elements into quadrilateral or prismatic elements is also described. The gains in computational efficiency afforded by the use of non-simplicial meshes over fully tetrahedral meshes are demonstrated through several examples.

Mavriplis, D. J.↗

Design and implementation of a parallel unstructured Euler solver using software primitives

This paper is concerned with the implementation of a three-dimensional unstructured-grid Euler solver on massively parallel distributed-memory computer architectures. The goal is to minimize solution time by achieving high computational rates with a numerically efficient algorithm. An unstructured multigrid algorithm with an edge-based data structure has been adopted, and a number of optimizations have been devised and implemented to accelerate the parallel computational rates. The implementation is carried out by creating a set of software tools, which provide an interface between the parallelization issues and the sequential code, while providing a basis for future automatic run-time compilation support. Large practical unstructured grid problems are solved on the Intel iPSC/860 hypercube and Intel Touchstone Delta machine. The quantitative effects of the various optimizations are demonstrated, and we show that the combined effect of these optimizations leads to roughly a factor of 3 performance improvement. The overall solution efficiency is compared with that obtained on the Cray Y-MP vector supercomputer.

Das, R.↗

Simulation of Conformal Spiral Slot Antennas on Composite Platforms

During the course of the grant, we wrote and distributed about 12 reports and an equal number of journal papers supported fully or in part by this grant. The list of reports (title & abstract) and papers are given in Appendices A and B. This grant has indeed been instrumental in developing a robust hybrid finite element method for the analysis of complex broadband antennas on doubly curved platforms. Previous to the grant, our capability was limited to simple printed patch antennas on mostly planar platforms. More specifically: (1) mixed element formulations were developed and new edge-based prisms were introduced; (2) these elements were important in permitting flexibility in geometry gridding for most antennas of interest; (3) new perfectly matched absorbers were introduced for mesh truncations associated with highly curved surfaces; (4) fast integral algorithms were introduced for boundary integral truncations reducing CPU time from O(N-2) down to O(N-1.5) or less; (5) frequency extrapolation schemes were developed for efficient broadband performance evaluations. This activity has been successfully continued by NASA researchers; (6) computer codes were developed and extensively tested for several broadband configurations. These include FEMA-CYL, FEMA-PRISM and FEMA-TETRA written by L. Kempel, T. Ozdemir and J. Gong, respectively; (7) a new infinite balun feed was designed nearly constant impedance over the 800-3000 MHz operational band; (8) a complete slot spiral antenna was developed, fabricated and tested at NASA Langley. This new design is a culmination of the projects goals and integrates the computational and experimental efforts. this antenna design resulted in a U.S. patent and was revised three times to achieve the desired bandwidth and gain requirements from 800-3000 MHz.

Volakis, J. L.↗

Hypersonic Flow Computations on Unstructured Meshes

A method for computing inviscid hypersonic flow over complex configurations using unstructured meshes is presented. The unstructured grid solver uses an edge{based finite{volume formulation. Fluxes are computed using a flux vector splitting scheme that is capable of representing constant enthalpy solutions. Second{order accuracy in smooth flow regions is obtained by linearly reconstructing the solution, and stability near discontinuities is maintained by locally forcing the scheme to reduce to first-order accuracy. The implementation of the algorithm to parallel computers is described. Computations using the proposed method are presented for a sphere-cone configuration at Mach numbers of 5.25 and 10.6, and a complex hypersonic re-entry vehicle at Mach numbers of 4.5 and 9.8. Results are compared to experimental data and computations made with established structured grid methods. The use of the solver as a screening tool for rapid aerodynamic assessment of proposed vehicles is described.

Bibb, K. L.↗

An Accuracy Evaluation of Unstructured Node-Centred Finite Volume Methods

Node-centred edge-based finite volume approximations are very common in computational fluid dynamics since they are assumed to run on structured, unstructured and even on mixed grids. We analyse the accuracy properties of both first and second derivative approximations and conclude that these schemes can not be used on arbitrary grids as is often assumed. For the Euler equations first-order accuracy can be obtained if care is taken when constructing the grid. For the Navier-Stokes equations, the grid restrictions are so severe that these finite volume schemes have little advantage over structured finite difference schemes. Our theoretical results are verified through extensive computations.

Svard, Magnus↗

Image Processing for Planetary Limb/Terminator Extraction

A novel image segmentation technique for extracting limb and terminator of planetary bodies is proposed. Conventional edge- based histogramming approaches are used to trace object boundaries. The limb and terminator bifurcation is achieved by locating the harmonized segment in the two equations representing the 2-D parameterized boundary curve. Real planetary images from Voyager 1 and 2 served as representative test cases to verify the proposed methodology.

image processing planetary limb planetary terminat↗

Orbiter Entry Aeroheating Working Group Viscous CFD Boundary Layer Transition Trailblazer Solutions

Boundary layer transition correlations for the Shuttle Orbiter have been previously developed utilizing a two-layer boundary layer prediction technique. The particular two-layer technique that was used is limited to Mach numbers less than 20. To allow assessments at Mach numbers greater than 20, it is proposed to use viscous CFD to the predict boundary layer properties. This report addresses if the existing Orbiter entry aeroheating viscous CFD solutions, which were originally intended to be used for heat transfer rate predictions, adequately resolve boundary layer edge properties and if the existing two-layer results could be leveraged to reduce the number of needed CFD solutions. The boundary layer edge parameters from viscous CFD solutions are extracted along the wind side centerline of the Space Shuttle Orbiter at reentry conditions, and are compared with results from the two-layer boundary layer prediction technique. The differences between the viscous CFD and two-layer prediction techniques vary between Mach 6 and 18 flight conditions and Mach 6 wind tunnel conditions, and there is not a straightforward scaling between the viscous CFD and two-layer values. Therefore: it is not possible to leverage the existing two-layer Orbiter flight boundary layer data set as a substitute for a viscous CFD data set; but viscous CFD solutions at the current grid resolution are sufficient to produce a boundary layer data set suitable for applying edge-based boundary layer transition correlations.

Wood, William A.↗

Three-Dimensional Navier-Stokes Calculations Using the Modified Space-Time CESE Method

The space-time conservation element solution element (CESE) method is modified to address the robustness issues of high-aspect-ratio, viscous, near-wall meshes. In this new approach, the dependent variable gradients are evaluated using element edges and the corresponding neighboring solution elements while keeping the original flux integration procedure intact. As such, the excellent flux conservation property is retained and the new edge-based gradients evaluation significantly improves the robustness for high-aspect ratio meshes frequently encountered in three-dimensional, Navier-Stokes calculations. The order of accuracy of the proposed method is demonstrated for oblique acoustic wave propagation, shock-wave interaction, and hypersonic flows over a blunt body. The confirmed second-order convergence along with the enhanced robustness in handling hypersonic blunt body flow calculations makes the proposed approach a very competitive CFD framework for 3D Navier-Stokes simulations.

Chang, Chau-lyan↗

Blunt-Body Entry Vehicle Aerothermodynamics: Transition and Turbulence on the CEV and MSL Configurations

Recent, current, and planned NASA missions that employ blunt-body entry vehicles pose aerothermodynamic problems that challenge the state-of-the art of experimental and computational methods. The issues of boundary-layer transition and turbulent heating on the heat shield have become important in the designs of both the Mars Science Laboratory and Crew Exploration Vehicle. While considerable experience in these general areas exists, that experience is mainly derived from simple geometries; e.g. sharp-cones and flat-plates, or from lifting bodies such as the Space Shuttle Orbiter. For blunt-body vehicles, application of existing data, correlations, and comparisons is questionable because an all, or mostly, subsonic flow field is produced behind the bow shock, as compared to the supersonic (or even hypersonic) flow of other configurations. Because of the need for design and validation data for projects such as MSL and CEV, many new experimental studies have been conducted in the last decade to obtain detailed boundary-layer transition and turbulent heating data on this class of vehicle. In this paper, details of several of the test programs are reviewed. The laminar and turbulent data from these various test are shown to correlate in terms of edge-based Stanton and Reynolds number functions. Correlations are developed from the data for transition onset and turbulent heating augmentation as functions of momentum thickness Reynolds number. These correlation can be employed as engineering-level design and analysis tools.

Hollis, Brian R.↗

Effects of Mesh Irregularities on Accuracy of Finite-Volume Discretization Schemes

The effects of mesh irregularities on accuracy of unstructured node-centered finite-volume discretizations are considered. The focus is on an edge-based approach that uses unweighted least-squares gradient reconstruction with a quadratic fit. For inviscid fluxes, the discretization is nominally third order accurate on general triangular meshes. For viscous fluxes, the scheme is an average-least-squares formulation that is nominally second order accurate and contrasted with a common Green-Gauss discretization scheme. Gradient errors, truncation errors, and discretization errors are separately studied according to a previously introduced comprehensive methodology. The methodology considers three classes of grids: isotropic grids in a rectangular geometry, anisotropic grids typical of adapted grids, and anisotropic grids over a curved surface typical of advancing layer grids. The meshes within the classes range from regular to extremely irregular including meshes with random perturbation of nodes. Recommendations are made concerning the discretization schemes that are expected to be least sensitive to mesh irregularities in applications to turbulent flows in complex geometries.

Diskin, Boris↗

Autonomous Onboard Science Data Analysis for Comet Missions

Coming years will bring several comet rendezvous missions. The Rosetta spacecraft arrives at Comet 67P/Churyumov-Gerasimenko in 2014. Subsequent rendezvous might include a mission such as the proposed Comet Hopper with multiple surface landings, as well as Comet Nucleus Sample Return (CNSR) and Coma Rendezvous and Sample Return (CRSR). These encounters will begin to shed light on a population that, despite several previous flybys, remains mysterious and poorly understood. Scientists still have little direct knowledge of interactions between the nucleus and coma, their variation across different comets or their evolution over time. Activity may change on short timescales so it is challenging to characterize with scripted data acquisition. Here we investigate automatic onboard image analysis that could act faster than round-trip light time to capture unexpected outbursts and plume activity. We describe one edge-based method for detect comet nuclei and plumes, and test the approach on an existing catalog of comet images. Finally, we quantify benefits to specific measurement objectives by simulating a basic plume monitoring campaign.

comets↗

Evaluation of the Minifilament-Eruption Scenario for Solar Coronal Jets in Polar Coronal Holes

Solar coronal jets are suspected to result from magnetic reconnection low in the Sun's atmosphere. Sterling et al. (2015) looked as 20 jets in polar coronal holes, using X-ray images from the Hinode/X-Ray Telescope (XRT) and EUV images from the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA). They suggested that each jet was driven by the eruption of twisted closed magnetic field carrying a small-scale filament, which they call a 'minifilament', and that the jet was produced by reconnection of the erupting field with surrounding open field. In this study, we carry out a more extensive examination of polar coronal jets. From 180 hours of XRT polar coronal hole observations spread over two years (2014-2016), we identified 130 clearly-identifiable X-ray jet events and thus determined an event rate of over 17 jets per day per in the Hinode/XRT field of view. From the broader set, we selected 25 of the largest and brightest events for further study in AIA 171, 193, 211, and 304 Angstrom images. We find that at least the majority of the jets follow the minifilament-eruption scenario, although for some cases the evolution of the minifilament in the onset of its eruption is more complex than presented in the simplified schematic of Sterling et al. (2015). For all cases in which we could make a clear determination, the spire of the X-ray jet drifted laterally away from the jet-base-edge bright point; this spire drift away from the bright point is consistent with expectations of the minifilament-eruption scenario for coronal-jet production. This work was supported with funding from the NASA/MSFC Hinode Project Office, and from the NASA HGI program.

holes↗

Economical Third-Order Methods for Accurate Surface Heating Predictions on Simplex Element Meshes

A node-centered edge-based finite volume discretization of the compressible Navier-Stokes equations is presented with the heat flux reformulated as a first order system. A dissipation vector is derived for the reformulated system, such that the heat flux can be upgraded to O(h^3) on simplex element meshes in the same fashion as the inviscid fluxes. The method of manufactured solutions is used to demonstrate this uniform order property in heat flux. This new system is shown to efficiently produce accurate surface heating predictions on hypersonic hemisphere flow using an anisotropic simplex element mesh, achieving O(h^3) accuracy at relatively low computational cost compared to similar methods.

Computational Fluid Dynamics↗