SEARCH · Search NASA

Results for “Navier Stokes solver”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗

DOE OSTI2023

Moment Representation of Regularized Lattice Boltzmann Methods on NVIDIA and AMD GPUs

The lattice Boltzmann method is a highly scalable Navier-Stokes solver that has been applied to flow problems in a wide array of domains. However, the method is bandwidth-bound on modern GPU accelerators and has a large memory footprint. In this paper, we present new 2D and 3D GPU implementations of two different regularized lattice Boltzmann methods, which are not only able to achieve an acceleration of ∼ 1.4 × w.r.t. reference lattice Boltzmann implementations but also reduce the memory requirements by up to 35% and 47% in 2D and 3D simulations respectively. These new approaches are evaluated on NVIDIA and AMD GPU architectures.

Valero Lara, Pedro↗

DOE OSTI2021

Scalable Interpolation on GPUs for Thermal Fluids Applications

The need for efficient interpolation arises in many use cases for spectral element simulations. Towards this end, we present a fast and robust GPU-enable interpolation utility for NekRS, a spectral-element, Navier-Stokes solver. This utility is based on the findpts and findpts eval routines from the gslib library. To demonstrate the functionality, we added support for particle tracking. The GPU implementation of the interpolation allows for tracking 100 times more particles for the same computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Aerodynamic Sensitivities over Separable Shape Tensors

Here, we present a comprehensive aerodynamic sensitivity analysis of airfoil parameterization informed by separable shape tensors. This parameterization approach uniquely benefits the design process by isolating various well-studied shape characteristics, such as airfoil thickness, and providing a well-regulated low-dimensional parameter domain for aerodynamic designs. Exploring the aerodynamic sensitivities of this novel parameterization can provide valuable insights for more robust designs and future manufacturing efforts. We construct a data-driven parameter space of airfoils using principal geodesic analysis of separable shape tensors informed by a curated database containing almost 20,000 suitable engineering airfoils. Analyzing the shape reconstruction error and the maximum mean discrepancy between joint distributions of aerodynamic quantities, we study the dimensionality of the learned parameter space. This simple numerical experiment demonstrates a dramatic dimension reduction that retains design effectiveness and promotes regularity of the shape representations. Finally, we generate new airfoils and use the HAM2D Reynolds-averaged Navier–Stokes solver to predict lift, drag, and moment coefficients. We compute multiple sensitivity metrics to quantify and assert the consistency of parameter influence on the aerodynamic quantities. We also explore low-dimensional polynomial ridge approximations to motivate physical intuitions and offer explanations of the approximated sensitivities.

17 WIND ENERGY↗

NASA NTRS1977

Numerical aerodynamic simulation facility preliminary study, volume 1

A technology forecast was established for the 1980-1985 time frame and the appropriateness of various logic and memory technologies for the design of the numerical aerodynamic simulation facility was assessed. Flow models and their characteristics were analyzed and matched against candidate processor architecture. Metrics were established for the total facility, and housing and support requirements of the facility were identified. An overview of the system is presented, with emphasis on the hardware of the Navier-Stokes solver, which is the key element of the system. Software elements of the system are also discussed.

Source record↗

NASA NTRS1978

Preliminary study for a Numerical Aerodynamic Simulation Facility

The state of the art of relevant technologies, of systems and processor architectures, and the measurable computational requirements of the two existing Navier-Stokes solution programs were assessed by Control Data Corporation to determine the best approach for designing a system for aerodynamic simulation. Standard parts and components were used to identify the support processing system, which is composed of commercially available equipment and software. Technological achievements in large scale integration technology and system organization of subcomponents borrowed from the STAR-100C project were used in the design of the Navier-Stokes solver.

Lincoln, N. R.↗

NASA NTRS1982

Numerical solution of Space Shuttle Orbiter flow field

The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.

Venkatapathy, E.↗

NASA NTRS1983

Theoretical investigations of high lift aerodynamics

A program which generates a coordinate system for a two element airfoil with the mesh points concentrated in areas of significant vorticity, i.e., boundary layer and wake is operational. The 'imbedded' grid method developed allows a transition from the scale of the main airfoil to the scale of the flap. This requirement is essential for the modeling of viscous flows over the flap and slat of a multielement airfoil. An airfoil mounted in a 2-D wind tunnel was formulated. The program is ready for a fine grid and a large number of planes to explore the characteristics of a Navier-Stokes solver in a quasi-3D case. The program was converted to a form suitable for the STAR computer. Runs were made to map a three dimensional flow field for a wall airfoil intersection with and without lift.

Bennett, G.↗

NASA NTRS1984

Comparison of secondary flows predicted by a viscous code and an inviscid code with experimental data for a turning duct

A comparison of the secondary flows computed by the viscous Kreskovsky-Briley-McDonald code and the inviscid Denton code with benchmark experimental data for turning duct is presented. The viscous code is a fully parabolized space-marching Navier-Stokes solver while the inviscid code is a time-marching Euler solver. The experimental data were collected by Taylor, Whitelaw, and Yianneskis with a laser Doppler velocimeter system in a 90 deg turning duct of square cross-section. The agreement between the viscous and inviscid computations was generally very good for the streamwise primary velocity and the radial secondary velocity, except at the walls, where slip conditions were specified for the inviscid code. The agreement between both the computations and the experimental data was not as close, especially at the 60.0 deg and 77.5 deg angular positions within the duct. This disagreement was attributed to incomplete modelling of the vortex development near the suction surface.

Schwab, J. R.↗

NASA NTRS1984

Numerical simulation of the tip vortex off a low-aspect-ratio wing at transonic speed

The viscous transonic flow around a low aspect ratio wing was computed by an implicit, three dimensional, thin-layer Navier-Stokes solver. The grid around the geometry of interest is obtained numerically as a solution to a Dirichlet problem for the cube. A low aspect ratio wing with large sweep, twist, taper, and camber is the chosen geometry. The topology chosen to wrap the mesh around the wing with good tip resolution is a C-O type mesh. The flow around the wing was computed for a free stream Mach number of 0.82 at an angle of attack of 5 deg. At this Mach number, an oblique shock forms on the upper surface of the wing, and a tip vortex and three dimensional flow separation off the wind surface are observed. Particle path lines indicate that the three dimensional flow separation on the wing surface is part of the roots of the tip vortex formation. The lifting of the tip vortex before the wing trailing edge is observed by following the trajectory of particles release around the wing tip.

Mansour, N. N.↗

NASA NTRS1984

Numerical simulation of the tip vortex off a low-aspect-ratio wing at transonic speed

The viscous transonic flow around a low-aspect-ratio wing has been computed using an implicit, three-dimensional, 'thin-layer' Navier-Stokes solver. The grid around the geometry of interest is obtained numerically as a solution to a Dirichlet problem for the cube. The geometry chosen for this study is a low-aspect-ratio wing with large sweep, twist, taper, and camber. The topology chosen to wrap the mesh around the wing with good tip resolution is a C-O type mesh. Using this grid, the flow around the wing was computed for a free-stream Mach number of 0.82 at an angle of attack of 5 deg. At this Mach number, an oblique shock forms on the upper surface of the wing, and a tip vortex and three-dimensional flow separation off the wing surface are observed. Particle path lines indicate that the three-dimensional flow separation on the wing surface is part of the roots of the tip-vortex formation. The lifting of the tip vortex before the wing trailing edge is clearly observed by following the trajectory of particles released around the wing tip.

Mansour, N. N.↗

NASA NTRS1983

Comparison of secondary flows predicted by a viscous code and an inviscid code with experimental data for a turning duct

A comparison of the secondary flows computed by the viscous Kreskovsky-Briley-McDonald code and the inviscid Denton code with benchmark experimental data for turning duct is presented. The viscous code is a fully parabolized space-marching Navier-Stokes solver while the inviscid code is a time-marching Euler solver. The experimental data were collected by Taylor, Whitelaw, and Yianneskis with a laser Doppler velocimeter system in a 90 deg turning duct of square cross-section. The agreement between the viscous and inviscid computations was generally very good for the streamwise primary velocity and the radial secondary velocity, except at the walls, where slip conditions were specified for the inviscid code. The agreement between both the computations and the experimental data was not as close, especially at the 60.0 deg and 77.5 deg angular positions within the duct. This disagreement was attributed to incomplete modeling of the vortex development near the suction surface.

Schwab, J. R.↗

NASA NTRS1985

A parabolized Navier-Stokes algorithm for separated supersonic internal flows

A stable global-iteration procedure is applied to a fully implicit noniterative parabolized Navier-Stokes algorithm to improve the solutions of the single-sweep method as well as to resolve small separation bubbles due to shock-wave/boundary-layer interactions. A modified form of the FLARE approximation is employed in the separated regions and some requirements for stable iterations are discussed. The results of the global-iteration scheme for two test cases of supersonic internal flows are compared against the results of the single-sweep method as well as the results of the full Navier-Stokes solver.

Chitsomboon, T.↗

NASA NTRS1985

Numerical study of Reynolds stress in compressible flows

A second order closure has been implemented in an implicit Navier-Stokes solver to study the behavior of the Reynolds stresses under the influence of severe pressure gradients. In the boundary layer zone, the strongly sheared character of the mean flow dominates the turbulence generation mechanisms. However, the pressure gradients also play a very important role for these processes, but at different locations within the boundary layer.

Vandromme, D.↗

NASA NTRS1985

The behaviour of turbulence anisotropy through shock waves and expansions

A second order closure has been implemented in an implicit Navier-Stokes solver to study the behavior of the Reynolds stresses under the influence of severe pressure gradients. In the boundary layer zone, the strongly sheared character of the mean flow dominates the turbulence generation mechanisms. However, the pressure gradients play also a very important role for these processes, but at different locations within the boundary layer. This aspect may be emphasized by the analysis of turbulence anisotropy through shock waves and expansions.

Minh, H. H.↗

NASA NTRS1986

Real gas effects on the numerical simulation of a hypersonic inlet

The present use of a parabolized Navier-Stokes solver to accurately simulate the flowfield in a supersonic inlet yields good agreement between numerical analysis and experiment for a Mach 7.4 inlet under cruise conditions, with an internal compression ratio of 8. The significance of real gas effects on the performance calculation of a hypersonic inlet is demonstrated, with small changes in the ratio of specific heats resulting in a substantial change in the calculated pitot pressure ratio.

Ng, W.-F.↗

NASA NTRS1986

The calculation of supersonic viscous flows using the parabolized Navier-Stokes equations

Solution of the parabolic Navier-Stokes (PNS) equations for supersonic flows is discussed, and compatibility of the PNS method with the triple-deck theory of Stewartson (1974) is demonstrated. Characteristic and stability analyses show that use of an appropriate filter on the pressure term in the x-momentum equation can suppress the departure solutions, giving the usual desired weak interaction solution. An Alternating Direction Explicit procedure, with minimal computer storage requirements compared to the full Navier-Stokes solvers, is proposed to calculate strongly interacting flows using a global iteration procedure for the PNS equations. The PNS equations are used to solve the hypersonic viscous interaction problem, and good agreement is found with experimental results.

Davis, R. T.↗

NASA NTRS1987

Aerodynamic design modification of a hypersonic wind tunnel nozzle by CSCM with high order accuracy

An improved version of the CSCM implicit Navier-Stokes solver with a flexible data structure has been applied, in conjunction with a sophisticated patched grid system, to the aerodynamic analysis of hypersonic axisymmetric contoured nozzles of the NASA-Ames 3.5-foot hypersonic wind tunnel. The results obtained are found to be in good agreement with experiment. Based on the analysis, a newly designed throat contour for the Mach 14 nozzle has been recommended.

Yang, J. Y.↗