Results for “Computational geometry” — Page 25

NASA NTRS1987

A computational analysis of flow separation over five different airfoil geometries at high angles-of-attack

The flux splitting method of Van Leer (1982) and an implicit, upwind-biased, finite-volume scheme of Rumsey (1985) for solving the two-dimensional Reynolds-averaged Navier-Stokes equations were applied to analyze five airfoil geometries (the NACA 0012, 12-percent Joukowski, NACA 4412, NACA 65-015, and LS1-0417 sections), at high angles of attack. The effect of grid density and grid extent on the Strouhal number of the unsteady flow was computed, as well as the effect of variations in Mach number, Reynolds number, and angle of attack. It was found, that past an angle of attack of 20 deg, Strouhal numbers were independent of angle of attack (at a constant value of 0.115), agreeing well with the experimental findings of Tyler (1931).

Rumsey, Christopher L.↗

NASA NTRS1995

Experience with k-epsilon turbulence models for heat transfer computations in rotating

This viewgraph presentation discusses geometry and flow configuration, effect of y+ on heat transfer computations, standard and extended k-epsilon turbulence model results with wall function, low-Re model results (the Lam-Bremhorst model without wall function), a criterion for flow reversal in a radially rotating square duct, and a summary.

Tekriwal, Prabbat↗

DOE OSTI2015

A mesh partitioning algorithm for preserving spatial locality in arbitrary geometries

Highlights: •An algorithm for partitioning computational meshes is proposed. •The Morton order space-filling curve is modified to achieve improved locality. •A spatial locality metric is defined to compare results with existing approaches. •Results indicate improved performance of the algorithm in complex geometries. -- Abstract: A space-filling curve (SFC) is a proximity preserving linear mapping of any multi-dimensional space and is widely used as a clustering tool. Equi-sized partitioning of an SFC ignores the loss in clustering quality that occurs due to inaccuracies in the mapping. Often, this results in poor locality within partitions, especially for the conceptually simple, Morton order curves. We present a heuristic that improves partition locality in arbitrary geometries by slicing a Morton order curve at points where spatial locality is sacrificed. In addition, we develop algorithms that evenly distribute points to the extent possible while maintaining spatial locality. A metric is defined to estimate relative inter-partition contact as an indicator of communication in parallel computing architectures. Domain partitioning tests have been conducted on geometries relevant to turbulent reactive flow simulations. The results obtained highlight the performance of our method as an unsupervised and computationally inexpensive domain partitioning tool.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1978

A geometry package for generation of input data for a three-dimensional potential-flow program

The preparation of geometric data for input to three-dimensional potential flow programs was automated and simplified by a geometry package incorporated into the NASA Langley version of the 3-D lifting potential flow program. Input to the computer program for the geometry package consists of a very sparse set of coordinate data, often with an order of magnitude of fewer points than required for the actual potential flow calculations. Isolated components, such as wings, fuselages, etc. are paneled automatically, using one of several possible element distribution algorithms. Curves of intersection between components are calculated, using a hybrid curve-fit/surface-fit approach. Intersecting components are repaneled so that adjacent elements on either side of the intersection curves line up in a satisfactory manner for the potential-flow calculations. Many cases may be run completely (from input, through the geometry package, and through the flow calculations) without interruption. Use of the package significantly reduces the time and expense involved in making three-dimensional potential flow calculations.

Halsey, N. D.↗

NASA NTRS1979

An evaluation of four single element airfoil analytic methods

A comparison of four computer codes for the analysis of two-dimensional single element airfoil sections is presented for three classes of section geometries. Two of the computer codes utilize vortex singularities methods to obtain the potential flow solution. The other two codes solve the full inviscid potential flow equation using finite differencing techniques, allowing results to be obtained for transonic flow about an airfoil including weak shocks. Each program incorporates boundary layer routines for computing the boundary layer displacement thickness and boundary layer effects on aerodynamic coefficients. Computational results are given for a symmetrical section represented by an NACA 0012 profile, a conventional section illustrated by an NACA 65A413 profile, and a supercritical type section for general aviation applications typified by a NASA LS(1)-0413 section. The four codes are compared and contrasted in the areas of method of approach, range of applicability, agreement among each other and with experiment, individual advantages and disadvantages, computer run times and memory requirements, and operational idiosyncrasies.

Freuler, R. J.↗

NASA NTRS1999

Computational Analysis of Compressibility Effects on a High-Lift Wing

The objective of this study was to investigate compressibility effects on a high-lift flowfield by simulating the flow about a three-dimensional multi-element wing. The computations were performed by solving both the incompressible and compressible Navier-Stokes equations (using the INS3D and OVERFLOW codes) on structured, overset grids. Turbulence was modeled via the one-equation, fully turbulent Spalart-Allmaras model. The computational results were validated with surface pressure measurements acquired at the NASA Ames 7- by 10-Foot Wind Tunnel. The geometry used for all computations consisted of an unswept wing in a landing configuration with a half-span flap and a three-quarter-span slat mounted inside a rectangular duct approximating the wind tunnel walls. The solutions were carefully examined to account for effects due to differences in algorithms. Compressibility effects were demonstrated by comparing surface particle traces, sectional pressure coefficient and boundary layer profile plots. It was found that small regions of compressibility near the slat and main-element leading edge can largely impact the flow. Even small compressibility regions can have significant global effects on the circulation and separation of each of the high-lift elements.

Baker, M. David↗

NASA NTRS1987

Conservative zonal schemes for patched grids in 2 and 3 dimensions

The computation of flow over complex geometries, such as realistic aircraft configurations, poses difficult grid generation problems for computational aerodynamicists. The creation of a traditional, single-module grid of acceptable quality about an entire configuration may be impossible even with the most sophisticated of grid generation techniques. A zonal approach, wherein the flow field is partitioned into several regions within which grids are independently generated, is a practical alternative for treating complicated geometries. This technique not only alleviates the problems of discretizing a complex region, but also facilitates a block processing approach to computation thereby circumventing computer memory limitations. The use of such a zonal scheme, however, requires the development of an interfacing procedure that ensures a stable, accurate, and conservative calculation for the transfer of information across the zonal borders.

Hessenius, Kristin A.↗

DOE OSTI2020

Investigation of the nacelle blockage effect for a downwind turbine

As downwind turbines garner increasing research interest, nacelle blockage becomes an important consideration. This paper examines nacelle blockage effects under a variety of Reynolds numbers, turbulence intensities, and nacelle geometries and proposes a computationally inexpensive nacelle blockage engineering model. Results show minimal nacelle blockage impacts on rotor loading and performance, except for tall (low-aspect-ratio) nacelles, which can cause an increase in rotor C p of up to 0.6%. In all cases, nacelle blockage increases rotor loading and performance metrics, however little. Finally, (more expensive) geometry-resolved and (less expensive) body-force computational fluid dynamics modeling techniques are compared, and body-force models are found to need improvement to adequately model nacelle blockage.

17 WIND ENERGY↗

DOE OSTI2024

Nonlinear elasticity with the Shifted Boundary Method

Here, we propose a new unfitted/immersed computational framework for nonlinear solid mechanics, which bypasses the complexities associated with the generation of CAD representations and subsequent body-fitted meshing. This approach allows to speed up the cycle of design and analysis in complex geometry and requires relatively simple computer graphics representations of the surface geometries to be simulated, such as the Standard Tessellation Language (STL format). Complex data structures and integration on cut elements are avoided by means of an approximate boundary representation and a modification (shifting) of the boundary conditions to maintain optimal accuracy. An extensive set of computational experiments in two and three dimensions is included.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1992

Radiative transfer in the electrode region of an arc heater

Radiation from hot gases in the electrode region of an arc heater are discussed. Major original contributions are made: (1) an analytical expression for the radiative heat flux, for an axisymmetric geometry and a grey reflecting wall, is derived, (2) the expression is rendered analytically differentiable/integrable by using appropriate techniques to remove the singularities, and (3), finally a two band model is devised to compute the radiative heat load at the walls in NASA Ames' jet facility. The advantage of using these techniques is that the time taken for the numerical integration is considerably reduced and computation of a 39 x 39 axisymmetric geometry becomes manageable. Fluid equations, representative of a high current arc heater, are numerically solved for an axisymmetric geometry. The computed average heat load of 400 Watts/sq cm using the two band model agrees well with the experimental estimates made at NASA Ames' facility.

Durgapal, P.↗

NASA NTRSDate not supplied

Mixed-element USM3D Contributions to the 4th AIAA High-Lift Prediction Workshop

This paper discusses results of the mixed-element USM3D (USM3D-ME) simulations performed for the 4th AIAA High-Lift Prediction Workshop. The workshop was separated into six Technical Focus Groups to investigate the impact of geometry modeling, grid, and computational methods for predicting high lift flows. This work was performed under the Fixed Grid RANS Technical Focus Group. The primary geometry selected for the workshop was the High-Lift Common Research Model. The performed simulations included a flap deflection study and both a grid refinement study and pitch sweep for the nominal flap deflection configuration. The results show that USM3D-ME RANS solutions, generally, tends to underpredict the lift coefficient and to predict a less negative pitching moment relative to the experimental data. The predicted drag coefficient values agree better with experiment for smaller angles of attack but were observed to be larger than experiment for the largest angle of attack simulated. The results of the grid refinement study demonstrated a lack of grid convergence for the provided grid family. The results of the grid refinement study are consistent with the submissions to the Fixed Grid and Mesh Adaptation Technical Focus Groups. Grid convergence for the provided grid family remains elusive for the international community. A 2D Multielement Airfoil configuration was included to enable a turbulence model verification study, which illustrated favorable agreement between USM3D-ME and the solutions provided by other flow solvers.

CFD↗

NASA NTRSDate not supplied

Mixed-Element USM3D Contributions to the 4th AIAA High-Lift Prediction Workshop

This paper discusses results of the mixed-element USM3D (USM3D-ME) simulations performed for the 4th AIAA High-Lift Prediction Workshop. The workshop was separated into six Technical Focus Groups to investigate the impact of geometry modeling, grid, and computational methods for predicting high lift flows. This work was performed under the Fixed Grid RANS Technical Focus Group. The primary geometry selected for the workshop was the High-Lift Common Research Model. The performed simulations included a flap deflection study and both a grid refinement study and pitch sweep for the nominal flap deflection configuration. The results show that USM3D-ME RANS solutions, generally, tends to underpredict the lift coefficient and to predict a less negative pitching moment relative to the experimental data. The predicted drag coefficient values agree better with experiment for smaller angles of attack but were observed to be larger than experiment for the largest angle of attack simulated. The results of the grid refinement study demonstrated a lack of grid convergence for the provided grid family. The results of the grid refinement study are consistent with the submissions to the Fixed Grid and Mesh Adaptation Technical Focus Groups. Grid convergence for the provided grid family remains elusive for the international community. A 2D Multielement Airfoil configuration was included to enable a turbulence model verification study, which illustrated favorable agreement between USM3D-ME and the solutions provided by other flow solvers.

CFD↗

DOE OSTI2019

Investigating Machine Learning Based X-Ray Computed Tomography Reconstruction Methods to Enhance the Accuracy of CT Scans

Development of new materials and predictive capabilities of component performance hinges on the ability to accurately digitize "as-built" geometries. X-ray computed tomography (CT) offers a non-destructive method of capturing these details but current methodologies are unable to produce the required fidelity for critical component certification. This project focused on discovering the limitations of existing CT reconstruction algorithms and exploring machine learning (ML) methodologies to overcome these limitations. We found that existing CT reconstruction methods are insufficient for Sandia's critical component certification process and that ML algorithms are a viable path forward to improving the quality of CT images.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) for engineering flow measurements

Increasingly, fully 3D experimental measurements of flow in complex engineering geometries are required to validate computational fluid dynamics models that support and inform reactor design and licensing. One barrier to such measurements is the complexity of typical reactor components and subsequent lack of optical access in these systems. To overcome this, the deployment of coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) is explored for flow measurement in reactor thermal hydraulic components and model (scaled) systems. Through this methodology, fully 3D flow information (via PEPT) and detailed internal geometry (via CT) are captured in opaque systems such as pipes, rod bundles, packed beds, etc. The reconstructed flow field and geometry can then be overlain to reveal detailed flow features around internal structures within a given test section. This is enabled through the use of a combined preclinical PET/CT scanner with overlapping PET and CT fields of view. Such measurements are useful for characterizing flow inside such intricate nuclear thermal hydraulic components as core geometries and heat exchangers, among others, and providing valuable 3D validation data for CFD models. In this work, basic tests of this 3D flow/geometry mapping are presented, and the implications of such measurements are discussed. Further, preliminary measurements are made with both point sources and flow in a simple pipe flow geometry to evaluate the capabilities of this technique. PEPT and CT features are coregistered with up to 0.1 mm precision, and pipe flow mean velocity and Reynolds stresses are reconstructed with similar accuracy to previous PEPT demonstrations. The utility of PEPT/CT is shown herein, and suggestions for future measurements are made.

3D flow measurement↗

NASA NTRSDate not supplied

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗

NASA NTRSDate not supplied

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗

NASA NTRS1992

Analysis of a pneumatic forebody flow control concept about a full aircraft geometry

A full aircraft geometry is used to computationally analyze the effectiveness of a pneumatic forebody flow control concept. An overset grid technique is employed to model the aircraft and slot geometry. Steady-state solutions for both isolated forebody and full aircraft configurations are carried out using a thin-layer Navier-Stokes flow solver. A solution obtained using the full aircraft geometry and a flight sideslip condition investigates the effect of sideslip on the leading edge extention vortex burst point. A no-sideslip blowing solution using the isolated forebody at full-scale wind tunnel test conditions is compared with experimental data to determine the accuracy of the numerical method. A solution employing the full geometry and slot blowing at flight conditions is obtained.

Gee, Ken↗

DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗