Results for “unstructured” — Page 47

NASA NTRSDate not supplied

Verification of Viscous Goal-Based Anisotropic Mesh Adaptation

Adaptive unstructured mesh techniques have a limited, but growing impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic mesh adaptation mechanics. Anisotropic metric construction methods are evaluated with analytically defined primal and adjoint fields. This allows the comparison of different metric formulations and different implementations of the same formulation without the complications of a flow and adjoint solution method. Unstructured mesh adaptation tools are verified by comparison on analytic primal and dual field before verification on benchmark aerodynamics cases. The documentation of these verification exercises helps to prepare these goal-based methods for routine use in production simulation workflows.

Mesh adaptation↗

NASA NTRS2020

Systems Development, Datamining and Knowledge Discovery Abstract

This summer, 2020, during my NASA internship I worked with my mentor, Ali Shaykhian, as well as a group of four other interns: Javel Gramling, Janelisse Morales, Tristian Running Crane, and Zulmarie Jiménez. Our research and projects all differ but work together to solve datamining unstructured data into an easy to read format. Turning unstructured data into something easier to follow is important for helping quickly pull data out of larger documents, so that one doesn’t have to go through multiple pages to find certain data points. By being able to structure data pulled from a document, it can be usedto collect data from mass amounts of forms and arrange it in an easy to glance at table instead of multiple forms. I chose to focus mainly on creating form templates with both Microsoft Word and Excel and getting used to the types of data that can be collected; as well as learning where both programs differed. After I was familiar with what could be gathered, I worked towards taking data collected by a Word form and importing it into an Excel spreadsheet. By being able to transfer data from a Word document to an Excel document, there is an added layer of functionality to the datamining. Moving data around between Excel sheets isn’t that complex of a process, but when you try to import from a Word document a lot of formatting and readability can be lost. The purpose of my research is to reduce that loss by using Visual Basic scripts to clean and arrange imported data.

Makayla Amber Renfro↗

NASA NTRSDate not supplied

Verification of Viscous Goal-Based Anisotropic Mesh Adaptation

Adaptive unstructured mesh techniques have a limited, but growing impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, anisotropic metric construction methods, and anisotropic mesh adaptation mechanics. A key ingredient for the broader acceptance of unstructured mesh adaptation is the verification of these implementations. Anisotropic metric construction methods are evaluated with analytically defined primal fields and the corresponding entropy variables as adjoint fields. This allows the comparison of different metric formulations and different implementations of the same formulation without the complications of a flow and adjoint solution method. The convergence of the output associated with the entropy variable adjoint is studied for mesh adaptation to these fields and a manufactured solution. Mesh adapted drag output is studied for two simple wings in compressible laminar flow to show fine-mesh convergence of multiple metric construction methods to less than a single drag count. The documentation of these verification exercises helps to prepare these goal-based methods for routine use in more complex simulations for production workflows.

mesh adaptation↗

NASA NTRSDate not supplied

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the acceleration and separation regions of the flow. An unstructured finite-volume solver is used along with an equilibrium wall model. The Mach number of the oncoming flow selected for the simulation is 0.176, and two Reynolds numbers (Re) are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Spanwise periodic results for the lower Re case are compared with the available DNS data, while the full three-dimensional simulation results for the higher Re case are compared to available experimental data. Sensitivity of WMLES results are assessed for a number of factors including grid resolution, wall model exchange location, type of wall model, and unstructured grid topology. The WMLES results for the lower Re case agree well with available DNS data in terms of the wall pressure variation, velocity and turbulent stress profile comparisons. The skin-friction predictions show a reduced tendency to separate. This appears to be accentuated for the higher Re case, for the medium resolution grid used in this study, where WMLES does not predict any flow separation in the mean, which is in contrast to the large separation bubble observed in the experiments.

turbulence↗

NASA NTRS2021

Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities such as line generation for general unstructured grids and optimization of the iterative solver. The efficiency of the new developments is assessed in a massively parallel environment. The assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Three benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere cylinder configuration, a transonic flow around the ONERA M6 wing, and a supersonic flow through a square duct are considered.

Mohagna J Pandya↗

NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗

NASA NTRSDate not supplied

The Influence of Adaptive Mesh Refinement on the Prediction of Vortex Interactions about a Generic Missile Airframe

The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.

slender body aerodynamics↗

NASA NTRSDate not supplied

Derivation of Effective Properties Based on Porous Scale Simulations Using Filtering Techniques

This study presents a method for derivation of effective properties at the interface and in-depth of porous materials. The method defines a Representative Elementary Volume (REV) and applies filtering techniques to computer effective properties such as porosity and flow quantities, such as velocity and pressure. The script, developed to process the data was tested on the VTK type files that contain the mesh information and the flow solution. The method allows to choose between two types of filters, such as cellular and top-hat and define the size of the REV and number of samples along the domain. Extraction of the REV from the domain is performed to exact boundaries requested for the user. This is done using a triangulation technique and cutting through the cells to comply to the requested boundaries of the volume. The method can be applied to both structured and unstructured meshes. Filtering the material porosity and flow quantities involves integration of the numerical data. The algorithm provides three integration methods, such as Riemann sum, Monte Carlo and Quadrature rule to perform the integration. The Monte-Carlo technique permits the use of either uniform or linearly spaced distribution of points. The Quadrature rule is currently applicable to tetrahedral element types. The Monte Carlo and Quadrature rule methods require interpolation of the flow quantities at the sample points. For interpolation, two methods were tested and are readily available, Gaussian interpolation and re-sampling. It has been shown that re-sampling method has better consistency and acceptable accuracy in interpolation of the data. The algorithm was written in Python language and uses a number of modules. The major module besides numpy is PyVista. It is used to process the computational domain, clip the REV and interpolate the data. Quadrature rule integration was performed using a quadpy module. ParaView software was used externally to convert the flow solution to the VTK (or more specifically VTU) format. Integration of ParaView in the same environment with PyVista encountered problems and could not be implemented in this work. The developed algorithm is expected to be applicable to unstructured meshes and more complex porous structures as soon as the data can be passed in VTK type format. With the report is provided Python script for filtering the solution and a Matlab script for simple generation and processing of 2-D and 3-D porous channel geometries. The two scripts don't communicate.

Alexsander Zibitsker↗

NASA NTRSDate not supplied

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M.V. Fischels↗

NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗

NASA NTRSDate not supplied

Wall-modeled LES of the Three-dimensional Speed Bump Experiment

We evaluate the performance of wall-modeled large eddy simulation (WMLES) in predicting smooth-body turbulent flow separation for a three-dimensional Gaussian-shaped speed bump geometry. The Reynolds number based on the bump length is 2 million, and an unstructured compressible finite-volume solver is used with an equilibrium wall model and dynamic subgrid scale model. Spanwise periodic simulations of the centerline two-dimensional bump were used to assess grid resolution requirements, and resolving the thin internal layer in the accelerating region was found to be necessary to correctly capture the downstream separated flow region. Based on these insights, an optimized grid that was smaller by a factor of two provided comparable accuracy to a finer grid that has been used by us and other researchers in past studies. Using insights gained from the spanwise periodic simulation, an unstructured polyhderal grid with about 250 million cells was used for the three-dimensional (3D) configuration with inviscid tunnel side and top walls. Detailed comparisons of wall skin-friction coefficient, wall pressure, velocity and turbulent stresses with available experimental data indicated excellent agreement. While the low-Reynolds Number Spalart-Allmaras RANS model with rotation/curvature correction gave qualitatively good agreement with experiments, WMLES showed significantly more accurate quantitative predictions of the flowfield. The 3D WMLES showed good agreement with experiments in the separated region, and the centerplane results indicated a different separation topology compared to the spanwise periodic simulation.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗

NASA NTRSDate not supplied

Time-Accurate Computational Fluid Dynamics Analysis of a Supersonic Turbine to Assess Structural Mode Excitation

Gas turbines can experience significant unsteady loading during operation which can cause damage such as blade cracking or bearing wear. While computational fluid dynamics (CFD) modeling of full three-dimensional gas turbine flow to verify meanline design performance is standard practice within the industry, the use of time-accurate CFD makes it possible to assess the unsteady loading which can be especially problematic at operating conditions where unsteady fluid loading aligns with structural modes. The Fluid Dynamics branch (ER42) at NASA Marshall Space Flight Center has developed a time-accurate CFD analysis methodology using the unstructured, density-based, Loci-CHEM solver which is capable of simulating the relative motion required for turbomachinery applications through sliding interface meshes. The highly-parallel unstructured solver allows for full three-dimensional, time-accurate simulation of turbines including upstream and downstream components such as manifolds, nozzles, and stators, any of which could excite a turbine structural mode. Recently, Loci-CHEM was used to analyze the FASTRAC supersonic turbine to assess the loading on the rotor. Overall performance and time-averaged loading on the rotor was computed, including the axial force over the entire rotor disk which is required for force balance calculations on the turbopump. The FASTRAC turbine has supersonic nozzles upstream of the rotor that result in significant unsteady loading on the rotor blades at multiples of the nozzle pass frequency. The unsteady blade pressure loading was decomposed into the frequency domain for direct use in structural forced response analysis. Utilizing this type of CFD analysis coupled with structural analysis has proven to be a valuable tool in identifying structural modes that are being excited by unsteady fluid loading.

M V Fischels↗

NASA NTRSDate not supplied

USM3D-ME Revised SA-neg Implementation and Assessment for Seventh Drag Prediction Workshop Solutions

Simulations performed by an unstructured-grid, cell-centered, finite-volume, mixed-element flow solver USM3D-MEwere presented at the Seventh Drag Prediction Workshop (DPW-7). Reynolds-averaged Navier-Stokes (RANS) solutions for the NASA high-speed common research model wing-body configuration were computed using a new implementation of the negative variant of the Spalart-Allmaras turbulence model (SA-neg) with hierarchical adaptive nonlinear iteration method (HANIM).ThisUSM3D-ME enhancement improved efficiency and robustness of RANS solutions resulting in fast iterative convergence to machine zero residuals on mixed-element unstructured grids provided by the DPW-7 organizing committee. This paper presents theoretical foundations for the new implementation of the SA-neg model and analyzesUSM3D-ME solutions for flow conditions corresponding to the required DPW-7 cases.

Boris Diskin↗

NASA NTRSDate not supplied

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

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

Grid Adaptation↗

NASA NTRSDate not supplied

High-Fidelity Simulations of Lift+Cruise VTOL Urban Air Mobility Concept Aircraft in Hover

This paper presents progress of high-fidelity multidisciplinary simulations for the NASA lift+cruise vertical takeoff and landing (VTOL) urban air mobility concept aircraft in hover. The study focuses on the unsteady flow solutions featuring strong wake interactions between rotors and the airframe. The simulations couple a high-fidelity aerodynamic model with a comprehensive rotorcraft aeromechanics tool with rotor rotation-speed trim. The aerodynamic model is based on the Reynolds-averaged Navier-Stokes equations using the one-equation Spalart-Allmaras turbulence model with rotation correction. Unsteady aerodynamic flows are computed on a dynamic, deformable, unstructured, overset grid system. An integrated overset-grid assembler is used to construct the composite grid from 23 component grids and facilitate communications between individual component grids. Rotor performance and airframe forces and moments are computed and compared on two unstructured grids. The study demonstrates capabilities of high-fidelity multidisciplinary analysis tools to capture the strong unsteady flowfield around this multi-rotor aircraft in hover.

CFD↗

NASA NTRSDate not supplied

Validation of the Hypersolve CFD Solver for Entry Descent and Landing Applications

The functional equivalence of the HyperSolve unstructured edge-based, finite-volume computational fluid dynamics code to the Langley Aerothermodynamic Upwind Relaxation Algorithm multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing community. A suite of cases using a range of thermochemical gas models on relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flowfield quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. The functional equivalence of the HyperSolve unstructured edge-based finite-volume computational fluid dynamics (CFD) code to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. A suite of cases using a range of thermochemical gas models on EDL-relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flow field quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. In general, HyperSolve predictions for surface pressure and surface heat flux are in close agreement with those predicted by LAURA.

hypersolve↗

NASA NTRSDate not supplied

Assessment of UMUSCL Scheme for DNS of Turbulent Flows

Direct numerical simulation (DNS) is performed using an unstructured node-based finite-volume methodology for canonical turbulent flow configurations to assess the effect of the upwinding parameter (κ) in the unstructured monotonic upstream scheme for conservation laws (UMUSCL) reconstruction scheme. The NASA FUN3D solver is used, and four turbulent validation cases are considered: periodic Taylor-Green vortex, (2) flow through a plane channel, (3) flow through a channel with a constriction, and (4) flow over a flat plate. The results are presented for κ = 0.5, 0.9 and 1 on hexahedral grids, with the last value corresponding to a zero-numerical dissipation scheme. Detailed comparisons with available reference data indicate that for the smallest length-scales simulated, κ = 0.5 is too dissipative while κ = 1 could lead to an energy build-up. Overall, κ =0.9 yields satisfactory results in terms of accuracy and robustness for all the cases considered in this study.

Direct Numerical Simulation↗