Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid edge”

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 271 records · Page 15

Parallel implementation of an adaptive scheme for 3D unstructured grids on the SP2

Dynamic mesh adaption on unstructured grids is a powerful tool for computing unsteady flows that require local grid modifications to efficiently resolve solution features. For this work, we consider an edge-based adaption scheme that has shown good single-processor performance on the C90. We report on our experience parallelizing this code for the SP2. Results show a 47.0X speedup on 64 processors when 10 percent of the mesh is randomly refined. Performance deteriorates to 7.7X when the same number of edges are refined in a highly-localized region. This is because almost all the mesh adaption is confined to a single processor. However, this problem can be remedied by repartitioning the mesh immediately after targeting edges for refinement but before the actual adaption takes place. With this change, the speedup improves dramatically to 43.6X.

Strawn, Roger C.↗

Parallel Implementation of an Adaptive Scheme for 3D Unstructured Grids on the SP2

Dynamic mesh adaption on unstructured grids is a powerful tool for computing unsteady flows that require local grid modifications to efficiently resolve solution features. For this work, we consider an edge-based adaption scheme that has shown good single-processor performance on the C90. We report on our experience parallelizing this code for the SP2. Results show a 47.OX speedup on 64 processors when 10% of the mesh is randomly refined. Performance deteriorates to 7.7X when the same number of edges are refined in a highly-localized region. This is because almost all mesh adaption is confined to a single processor. However, this problem can be remedied by repartitioning the mesh immediately after targeting edges for refinement but before the actual adaption takes place. With this change, the speedup improves dramatically to 43.6X.

Oliker, Leonid↗

A Numerical Evaluation of Icing Effects on a Natural Laminar Flow Airfoil

As a part of CFD code validation efforts within the Icing Branch of NASA Glenn Research Center, computations were performed for natural laminar flow (NLF) airfoil, NLF-0414. with 6 and 22.5 minute ice accretions. Both 3-D ice castings and 2-D machine-generated ice shapes were used in wind tunnel tests to study the effects of natural ice is well as simulated ice. They were mounted in the test section of the Low Turbulence Pressure Tunnel (LTPT) at NASA Langley that the 2-dimensionality of the flow can be maintained. Aerodynamic properties predicted by computations were compared to data obtained through the experiment by the authors at the LTPT. Computations were performed only in 2-D and in the case of 3-D ice, the digitized ice shape obtained at one spanwise location was used. The comparisons were mainly concentrated on the lift characteristics over Reynolds numbers ranging from 3 to 10 million and Mach numbers ranging from 0.12 to 0.29. WIND code computations indicated that the predicted stall angles were in agreement with experiment within one or two degrees. The maximum lift values obtained by computations were in good agreement with those of the experiment for the 6 minute ice shapes and the minute 3-D ice, but were somewhat lower in the case of the 22.5 minute 2-D ice. In general, the Reynolds number variation did not cause much change in the lift values while the variation of Mach number showed more change in the lift. The Spalart-Allmaras (S-A) turbulence model was the best performing model for the airfoil with the 22.5 minute ice and the Shear Stress Turbulence (SST) turbulence model was the best for the airfoil with the 6 minute ice and also for the clean airfoil. The pressure distribution on the surface of the iced airfoil showed good agreement for the 6 minute ice. However, relatively poor agreement of the pressure distribution on the upper surface aft of the leading edge horn for the 22.5 minute ice suggests that improvements are needed in the grid or turbulence models.

Chung, James J.↗

A Unified Overset Grid Generation Graphical Interface and New Concepts on Automatic Gridding Around Surface Discontinuities

For many years, generation of overset grids for complex configurations has required the use of a number of different independently developed software utilities. Results created by each step were then visualized using a separate visualization tool before moving on to the next. A new software tool called OVERGRID was developed which allows the user to perform all the grid generation steps and visualization under one environment. OVERGRID provides grid diagnostic functions such as surface tangent and normal checks as well as grid manipulation functions such as extraction, extrapolation, concatenation, redistribution, smoothing, and projection. Moreover, it also contains hyperbolic surface and volume grid generation modules that are specifically suited for overset grid generation. It is the first time that such a unified interface existed for the creation of overset grids for complex geometries. New concepts on automatic overset surface grid generation around surface discontinuities will also be briefly presented. Special control curves on the surface such as intersection curves, sharp edges, open boundaries, are called seam curves. The seam curves are first automatically extracted from a multiple panel network description of the surface. Points where three or more seam curves meet are automatically identified and are called seam corners. Seam corner surface grids are automatically generated using a singular axis topology. Hyperbolic surface grids are then grown from the seam curves that are automatically trimmed away from the seam corners.

Chan, William M.↗

Effects of Elevated Free-Stream Turbulence on Active Control of a Separation Bubble

The Effects of elevated free-stream turbulence (FST) on the natural and periodically excited separation bubbles were studied experimentally, due to the relevance of this flow to low-pressure turbine blades at low Reynolds numbers. A bubble was formed at the leading edge of a flat plate and the FST level was altered by placing a grid across the flow at different locations upstream of the plate. The mixing across the separated shear-layer, forming the free boundary of the bubble, increased due to the elevated FST and due to nominally two-dimensional periodic excitation, both flattening and shortening the bubble. Periodic excitation at frequencies that were at least an order of magnitude lower than those associated with the initial shear-layer instability, were very effective at low FST, because the amplitudes of the excitation frequency and its harmonic were amplified over the bubble. High frequency excitation (F+ 3, based on the length of the baseline low FST bubble) had a major effect close to the separation location, while farther downstream the excited fluctuations rapidly decayed in the reattachment region. Low frequency excitation, that generated waves comparable to the length of the unperturbed bubble (F+ 1) were less effective and their magnitude decayed at a slower rate downstream of reattachment. An increase in the level of the FST reduced the net effect of the periodic excitation on the mixing enhancement and subsequent reattachment process, probably due to a destructive interference between the nominally 2D excitation and the random (in space and time) FST, reducing the spanwise coherence and therefore the effectiveness of the current control strategy. However, even at the reduced effectiveness of 2D periodic excitation at elevated FST, it accelerated the reattachment process and the recovery rate of the reattached boundary layer, enhancing the boundary layer resistance to repeat separation and reducing its momentum loss further downstream.

Ashpis, D.↗

Fabrication of a Kilopixel Array of Superconducting Microcalorimeters with Microstripline Wiring

A document describes the fabrication of a two-dimensional microcalorimeter array that uses microstrip wiring and integrated heat sinking to enable use of high-performance pixel designs at kilopixel scales (32 X 32). Each pixel is the high-resolution design employed in small-array test devices, which consist of a Mo/Au TES (transition edge sensor) on a silicon nitride membrane and an electroplated Bi/Au absorber. The pixel pitch within the array is 300 microns, where absorbers 290 microns on a side are cantilevered over a silicon support grid with 100-micron-wide beams. The high-density wiring and heat sinking are both carried by the silicon beams to the edge of the array. All pixels are wired out to the array edge. ECR (electron cyclotron resonance) oxide underlayer is deposited underneath the sensor layer. The sensor (TES) layer consists of a superconducting underlayer and a normal metal top layer. If the sensor is deposited at high temperature, the ECR oxide can be vacuum annealed to improve film smoothness and etch characteristics. This process is designed to recover high-resolution, single-pixel x-ray microcalorimeter performance within arrays of arbitrarily large format. The critical current limiting parts of the circuit are designed to have simple interfaces that can be independently verified. The lead-to-TES interface is entirely determined in a single layer that has multiple points of interface to maximize critical current. The lead rails that overlap the TES sensor element contact both the superconducting underlayer and the TES normal metal

Chervenak, James↗

Shock-vortex interaction over a 65-degree delta wing in transonic flow

Transonic flow over a 65-deg swept-back, sharp-edged, cropped delta wing is investigated computationally using the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with an implicit, upwind, flux-difference splitting, finite-volume scheme. Coarse and fine O-H grids are used to obtain the solution. The grid consists of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively. The results are presented for an angle of attack of 20 deg Mach number of 0.85 and Reynolds number of 3.23 x 10 exp 6. With the fine grid, the results show that a system of shocks has been captured over the upper wing surface and that the leading-edge vortex core experiences an unsteady supersonic vortex breakdown after passing through a spanwise shock near the wing trailing edge. The computed results at a certain time are in good agreement with the experimental data. Topological aspects of the vortex breakdown flowfield are also presented and discussed.

Kandil, Osama A.↗

An efficient advancing front algorithm for Delaunay triangulation

There has been some recent interest in fluid dynamics calculations on unstructured meshes. One method of unstructured mesh generation involves Delaunay triangulation. This method has certain advantages but it can be expensive to implement. Furthermore, there can be problems with crossing grid lines near boundaries. A method shown here avoids many of the robustness and efficiency problems previously associated with Delaunay triangulation. As an added feature, a simple algorithm is shown which allows removal of diagonal edges from cells that are nearly rectangular. This can result in significant savings in the cost per iteration of a flow solver using this grid.

Merriam, Marshal L.↗

Polar microwave brightness temperatures from Nimbus-7 SMMR: Time series of daily and monthly maps from 1978 to 1987

A time series of daily brightness temperature gridded maps (October 25, 1978 through August 15, 1987) were generated from all ten channels of the Nimbus-7 Scanning Multichannel Microwave Radiometer orbital data. This unique data set can be utilized in a wide range of applications including heat flux, ocean circulation, ice edge productivity, and climate studies. Two sets of data in polar stereographic format are created for the Arctic region: one with a grid size of about 30 km on a 293 by 293 array similar to that previously utilized for the Nimbus-5 Electrically Scanning Microwave Radiometer, while the other has a grid size of about 25 km on a 448 by 304 array identical to what is now being used for the DMSP Scanning Multichannel Microwave Imager. Data generated for the Antaractic region are mapped using the 293 by 293 grid only. The general technique for mapping, and a quality assessment of the data set are presented. Monthly and yearly averages are also generated from the daily data and sample geophysical ice images and products derived from the data are given. Contour plots of monthly ice concentrations derived from the data for October 1978 through August 1987 are presented to demonstrate spatial and temporal detail which this data set can offer, and to show potential research applications.

Comiso, Josefino C.↗

An adaptive remeshing scheme for vortex dominated flows using three-dimensional unstructured grids

An adaptive remeshing procedure for vortex dominated flows is described, which uses three-dimensional unstructured grids. Surface grid adaptation is achieved using the static pressure as an adaptation parameter, while entropy is used in the field to accurately identify high vorticity regions. An emphasis has been placed in making the scheme as automatic as possible so that a minimum user interaction is required between remeshing cycles. Adapted flow solutions are obtained on two sharp-edged configurations at low speed, high angle-of-attack flow conditions. The results thus obtained are compared with fine grid CFD solutions and experimental data, and conclusions are drawn as to the efficiency of the adaptive procedure.

Parikh, Paresh↗

Quasi-optical grids with thin rectangular patch/aperture elements

Theoretical analysis is presented for an efficient and accurate performance evaluation of quasi-optical grids comprised of thin rectangular patch/aperture elements with/without a dielectric substrate/superstrate. The convergence rate of this efficient technique is improved by an order of magnitude with the approximate edge conditions incorporated in the basis functions of the integral equation solution. Also presented are the interesting applications of this efficient analytical technique to the design and performance evaluation of the coupling grids and beam splitters in the optical systems as well as thermal protection sunshields used in the communication systems of satellites and spacecrafts.

Wu, Te-Kao↗

Navier-Stokes analysis of blunt trailing edge airfoils

The flow around blunt trailing edge airfoils was studied by solving the Reynolds-averaged Navier-Stokes equations. The solution procedure combines a grid around the airfoil with a second grid for the wake so that the time advancement over the domain is fully implicit. This is not only very efficient for the algorithm but also allows implicit solutions of a one equation turbulence model appropriate for both boundary layers and wakes. An algebraic and two one-equation turbulence models are tested for a blunt RAE 2822 airfoil section and detailed comparisons with experimental data are presented in the trailing edge region.

Stanaway, Sharon↗

Navier-Stokes analysis of solid propellant rocket motor internal flows

A multidimensional implicit Navier-Stokes analysis that uses numerical solution of the ensemble-averaged Navier-Stokes equations in a nonorthogonal, body-fitted, cylindrical coordinate system has been applied to the simulation of the steady mean flow in solid propellant rocket motor chambers. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. The code was validated by comparing computed results with the experimental data obtained in cylindrical-port cold-flow tests. The agreement between the computed and experimentally measured mean axial velocities is excellent. The axial location of transition to turbulent flow predicted by the two-equation (k-epsilon) turbulence model used in the computations also agrees well with the experimental data. Computations performed to simulate the axisymmetric flowfield in the vicinity of the aft field joint in the Space Shuttle solid rocket motor using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot.

Sabnis, J. S.↗

Navier-Stokes calculation of transonic flow past the NTF 65 deg delta wing

Models of four delta wings were built and tested in the 8 x 8 ft Transonic Wind Tunnel at LaRC. The wings are identical in planform shape with a swept-back angle of 65 degrees, but bear different leading edge profiles. The models were tested under pressurized and cryogenic conditions to simulate true flight Reynolds numbers. Data on the aerodynamic forces and pressure distributions at various locations on the wings were taken at various flight Mach numbers and angles of attack. Effects of high Reynolds number and leading edge radius on the aerodynamic characteristics of the wings are being accessed. To thoroughly understand the turbulent, vortical flow around the wings, an effort to perform computational aerodynamic analysis is being made. The objective of the analysis is to supplement and validate the experimental data and explain the high Reynolds number and leading edge effects. GRIDGEN, a software developed by General Dynamics, is being used to generate the grid topology for the flow field around the wings. The flow solver to be used is CFL3D, a computation fluid dynamics code developed at LaRC. Based on the geometric description of the wings, a FORTRAN program called WINGSURF was written to generate the databases defining the surface geometry of the wings. To match the true geometry of the models in the wind tunnel for realistic comparison with experimental results, databases defining the sting support for the wing models were also created by two other FORTRAN programs, STINJOIN and STINREAR. Listing of the programs are attached and the geometry of a typical wing model with the sting support is shown. Other aspects of the investigation are discussed.

Wu, Chivey C.↗

Sea ice lead dynamics from ERS-1 SAR

The use of ERS-1 SAR (Synthetic Aperture Radar) images to measure ice lead changes and ice deformation, and thus to determine their relationships and compare them to theory, is reported. ERS-1 SAR data collected at the Alaska (U.S.) SAR facility are processed into low resolution geocoded images which are used as input to the Geophysical Processor System (GPS). The GPS uses cross correlation and edge detection on sequential pairs of images to track a regular 5 km array of grid points. The result for each sequential image pair is a field of ice displacement that contains, on average, more than 200 vectors. The GPS also classifies all the pixels of the images into one of four ice types based on the normalized backscatter coefficients and a seasonal look-up table. This allows all pixels to be labeled as either multiyear ice or as the first year ice and open water that comprise leads.

Stern, Harry L.↗

New multigrid approach for three-dimensional unstructured, adaptive grids

A new multigrid method with adaptive unstructured grids is presented. The three-dimensional Euler equations are solved on tetrahedral grids that are adaptively refined or coarsened locally. The multigrid method is employed to propagate the fine grid corrections more rapidly by redistributing the changes-in-time of the solution from the fine grid to the coarser grids to accelerate convergence. A new approach is employed that uses the parent cells of the fine grid cells in an adapted mesh to generate successively coaser levels of multigrid. This obviates the need for the generation of a sequence of independent, nonoverlapping grids as well as the relatively complicated operations that need to be performed to interpolate the solution and the residuals between the independent grids. The solver is an explicit, vertex-based, finite volume scheme that employs edge-based data structures and operations. Spatial discretization is of central-differencing type combined with a special upwind-like smoothing operators. Application cases include adaptive solutions obtained with multigrid acceleration for supersonic and subsonic flow over a bump in a channel, as well as transonic flow around the ONERA M6 wing. Two levels of multigrid resulted in reduction in the number of iterations by a factor of 5.

Parthasarathy, Vijayan↗

CFD Analysis of Modular Thrusters Performance

The effective performance of modular thrusters in an aerospike configuration is difficult to determine. Standard analytical tools are applicable to conventional nozzle shapes, but are limited when applied to an aerospike nozzle (An aerospike nozzle is an altitude compensating external nozzle). Three baseline nozzle shapes are derived using standard analytical procedures. The baseline nozzle sizes are restricted to fill a volume envelope. The three shapes are an axi-symmetric round nozzle, a two dimensional planar square exit nozzle, and a super elliptic round to nearly square nozzle. The integrated (thruster/aerospike) performance of the three nozzles is determined through the use of three dimensional viscous computational fluid dynamic (CFD) calculation where complex features of the flow field can be accurately captured. The resulting installed performance is then used to evaluate the efficiency of these nozzle shapes for aerospike applications. The determination of effective performance of a thruster nozzle integrated into an aerospike nozzle require the solution of the three dimensional turbulent Navier-Stokes equations. The model used in this study consisted of two zones; one of the upstream thruster cowl surface so freestream conditions can be accurately predicted, and two, the aerospike surface beginning with the thruster outflow and extending to the end of the aerospike surface. The numerical grid consisted of over 120,000 nodes and used symmetry on the thruster centerline and edge. A two species non-reacting chemistry model was used to capture the variation of fluid properties between the hot plume base and freestream air. From the results of the three baseline nozzle aerospike calculations, the effictive performance of the nozzle was determined. The flow fields of these calculations do show some variation between the cases. Recirculation zones on the cowl surface is predicted for the two dimensional planar nozzle and a smaller one for the super elliptic nozzle. The recirculation is caused by the strong pressure gradient between the plume and freestream flows. The axi-symmetric nozzle results indicate recirculation zones on the thruster face. These recirculation zones smooth the pressure gradient between the plume and freestream flow limiting the formation of recirculation on the cowl surface. Thruster to thruster interaction is evident for the axi-symmetric and supper elliptic calculation while the two dimensional planar nozzle did not have any lateral expansion in the nozzle, so thruster to thruster interaction is limited. The integrated performance results, at the altitude choosen, show very little variation between the three thruster shapes. This result allows for nozzle shape determination based on additional considerations (thermal, structural, weight) besides performance.

Ungewitter, Ronald J.↗