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 217 records · Page 12

Influence of numerical dissipation in computing supersonic vortex-dominated flows

Steady supersonic vortex-dominated flows are solved using the unsteady Euler equations for conical and three-dimensional flows around sharp- and round-edged delta wings. The computational method is a finite-volume scheme which uses a four-stage Runge-Kutta time stepping with explicit second- and fourth-order dissipation terms. The grid is generated by a modified Joukowski transformation. The steady flow solution is obtained through time-stepping with initial conditions corresponding to the freestream conditions, and the bow shock is captured as a part of the solution. The scheme is applied to flat-plate and elliptic-section wings with a leading edge sweep of 70 deg at an angle of attack of 10 deg and a freestream Mach number of 2.0. Three grid sizes of 29 x 39, 65 x 65 and 100 x 100 have been used. The results for sharp-edged wings show that they are consistent with all grid sizes and variation of the artificial viscosity coefficients. The results for round-edged wings show that separated and attached flow solutions can be obtained by varying the artificial viscosity coefficients. They also show that the solutions are independent of the way time stepping is done. Local time-stepping and global minimum time-steeping produce same solutions.

Kandil, O. A.↗

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗

Sensitivity of 30-cm mercury bombardment ion thruster characteristics to accelerator grid design

The design of ion optics for bombardment thrusters strongly influences overall performance and lifetime. The operation of a 30 cm thruster with accelerator grid open area fractions ranging from 43 to 24 percent, was evaluated and compared with experimental and theoretical results. Ion optics properties measured included the beam current extraction capability, the minimum accelerator grid voltage to prevent backstreaming, ion beamlet diameter as a function of radial position on the grid and accelerator grid hole diameter, and the high energy, high angle ion beam edge location. Discharge chamber properties evaluated were propellant utilization efficiency, minimum discharge power per beam amp, and minimum discharge voltage.

Rawlin, V. K.↗

Sensitivity of 30-cm mercury bombardment ion thruster characteristics to accelerator grid design

The design of ion optics for bombardment thrusters strongly influences overall performance and lifetime. The operation of a 30-cm thruster with accelerator grid open area fractions ranging from 43 to 24 percent, was evaluated and compared with previously published experimental and theoretical results. Ion optics properties measured included the beam current extraction capability, the minimum accelerator grid voltage to prevent backstreaming, ion beamlet diameter as a function of radial position on the grid and accelerator grid hole diameter, and the high energy, high angle ion beam edge location. Discharge chamber properties evaluated were propellant utilization efficiency, minimum discharge power per beam amp, and minimum discharge voltage.

Rawlin, V. K.↗

Finite element analysis of transonic flows in highly staggered cascades

Finite element method is employed in the solution of full-potential equation for analyzing steady, transonic flows through a series of highly staggered cascades. Convergence characteristics and efficiency of the developed numerical procedure are discussed. Accuracy of the mathematical model is investigated in terms of the employed computational grid. In particular, the importance of accurate representation of the trailing edge together with the shock on the obtained flow results are studied. Numerical results include choked flows, for a cascade of NACA 0012 airfoils with pitch/chord length = 1, and stagger angle of 45 deg.

Akay, H. U.↗

Euler calculations for wings using Cartesian grids

A method is presented for the calculation of transonic flows past wings using Cartesian grids. The calculations are based on a finite volume formulation of the Euler equations. Results are presented for a rectangular wing with a flat tip and the ONERA M6 wing. In general, the results are in good agreement with other computations and available experiment. However, Cartesian grids require a greater number of points than body fitted grids in order to resolve the flow properties near the leading edge of a swept wing.

Gaffney, R. L., Jr.↗

Computational validation of a parabolized Navier-Stokes solver on a sharp-nose cone at hypersonic speeds

Perfect gas computational results from a newly-developed upwind, parabolized Navier-Stokes (PNS) solver are compared with an existing set of experimental laminar results for a 10-deg half-angle circular cone at freestream Mach number of 7.95. Comparisons were performed with surface pressure and heat transfer data, as well as with flowfield pitot measurements. The PNS code predicted the surface quantities accurately up through 20-deg angle-of-attack, including crossflow separation, and correctly defined the location of the bow shock and the edge of the boundary layer. The importance of cell Reynolds number, grid density, and thermal boundary conditions to the accurate prediction of the flowfield are examined through numerical emamples.

Huebner, Lawrence D.↗

Numerical simulation of the viscous flow around a simplified F/A-18 at high angles of attack

A numerical method developed for solving the viscous flow around three-dimensional complex configurations is presently used to simulate the flow around a simplified F/A-18 configuration encompassing forebody, wing, leading-edge extension, faired-over inlet, and deflected wing leading-edge flaps, at Mach 0.243 and 30.3 deg angle of attack. The computational results show the details of the flowfield structure, including primary, secondary, and tertiary separation lines, the development of forebody and leading-edge extension vortex, and the burst of this vortex. A grid-refinement study is conducted to assess the effect of grid characteristics on solution accuracy. Substantial agreement is obtained between these results and flight data.

Rizk, Yehia M.↗

Navier-Stokes analysis of three-dimensional unsteady flows inside turbine stages

This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

Navier-Stokes analysis of three-dimensional flow and heat transfer inside turbine blade rows

A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

Convergence of Defect-Correction and Multigrid Iterations for Inviscid Flows

Convergence of multigrid and defect-correction iterations is comprehensively studied within different incompressible and compressible inviscid regimes on high-density grids. Good smoothing properties of the defect-correction relaxation have been shown using both a modified Fourier analysis and a more general idealized-coarse-grid analysis. Single-grid defect correction alone has some slowly converging iterations on grids of medium density. The convergence is especially slow for near-sonic flows and for very low compressible Mach numbers. Additionally, the fast asymptotic convergence seen on medium density grids deteriorates on high-density grids. Certain downstream-boundary modes are very slowly damped on high-density grids. Multigrid scheme accelerates convergence of the slow defect-correction iterations to the extent determined by the coarse-grid correction. The two-level asymptotic convergence rates are stable and significantly below one in most of the regions but slow convergence is noted for near-sonic and very low-Mach compressible flows. Multigrid solver has been applied to the NACA 0012 airfoil and to different flow regimes, such as near-tangency and stagnation. Certain convergence difficulties have been encountered within stagnation regions. Nonetheless, for the airfoil flow, with a sharp trailing-edge, residuals were fast converging for a subcritical flow on a sequence of grids. For supercritical flow, residuals converged slower on some intermediate grids than on the finest grid or the two coarsest grids.

Diskin, Boris↗

A direct-inverse transonic wing-design method in curvilinear coordinates including viscous-interaction

Progress in the direct-inverse wing design method in curvilinear coordinates has been made. A spanwise oscillation problem and proposed remedies are discussed. Test cases are presented which reveal the approximate limits on the wing's aspect ratio and leading edge wing sweep angle for a successful design, and which show the significance of spanwise grid skewness, grid refinement, viscous interaction, the initial airfoil section and Mach number-pressure distribution compatibility on the final design. Furthermore, preliminary results are shown which indicate that it is feasible to successfully design a region of the wing which begins aft of the leading edge and terminates prior to the trailing edge.

Ratcliff, Robert R.↗

Unstructured multigrid through agglomeration

In this work the compressible Euler equations are solved using finite volume techniques on unstructured grids. The spatial discretization employs a central difference approximation augmented by dissipative terms. Temporal discretization is done using a multistage Runge-Kutta scheme. A multigrid technique is used to accelerate convergence to steady state. The coarse grids are derived directly from the given fine grid through agglomeration of the control volumes. This agglomeration is accomplished by using a greedy-type algorithm and is done in such a way that the load, which is proportional to the number of edges, goes down by nearly a factor of 4 when moving from a fine to a coarse grid. The agglomeration algorithm has been implemented and the grids have been tested in a multigrid code. An area-weighted restriction is applied when moving from fine to coarse grids while a trivial injection is used for prolongation. Across a range of geometries and flows, it is shown that the agglomeration multigrid scheme compares very favorably with an unstructured multigrid algorithm that makes use of independent coarse meshes, both in terms of convergence and elapsed times.

Venkatakrishnan, V.↗

Evaluation of Wall-Modeled Les for Flow Over A Multi-Element Airfoil

We perform wall-modeled large-eddy simulations (WMLES) of turbulent flow over a 30P30N multielement airfoil at various angles of attack to assess its ability to accurately predict lift using an unstructured node-based finite-volume methodology. The flow conditions are based on the experiments reported in Klausmeyer & Lin [NASA /TM-112858]. The oncoming flow Mach number is 0.2, and the stowed chord-based Reynolds number (𝑅𝑒𝑐) is 9 million. The angle of attack (𝛼) is varied between 8◦ and 23◦ with the expected stall angle of around 21◦. The simulated geometry is periodic in the span, thus any three-dimensional effects present in the experiment are ignored. We first analyze the predictions with baseline isotropic and anisotropic grids containing over 16 grid points per boundary-layer thickness. The near-wall streamwise/nominal spacing ratio is 0.4 at the leading/trailing edges of the slat, main element and flap, and 1 elsewhere, and the first wall-normal/nominal spacing ratio is 0.375 over most of the airfoil for the anisotropic grid, while these are unity for the isotropic grid. It is found that the predicted lift and wall pressures, near-wall velocity profiles, and the location of onset of resolved turbulence are sensitive to the near-wall grid anisotropy. Further sensitivities were explored at lower angles of attack by independently varying the streamwise and wall-normal grid anisotropies. The effects of both off-wall and near-wall grid refinement were also quantified. For the grid resolutions used here, our results suggest that the first wall-normal spacing, which coincides with the WMLES exchange location, has a dominant effect on the predictions for this flow configuration.

Computational Fluid Dynamics↗

Reynolds-Averaged Navier-Stokes Simulations of Two Partial-Span Flap Wing Experiments

Structured Reynolds Averaged Navier-Stokes simulations of two partial-span flap wing experiments were performed. The high-lift aerodynamic and aeroacoustic wind-tunnel experiments were conducted at both the NASA Ames 7-by 10-Foot Wind Tunnel and at the NASA Langley Quiet Flow Facility. The purpose of these tests was to accurately document the acoustic and aerodynamic characteristics associated with the principle airframe noise sources, including flap side-edge noise. Specific measurements were taken that can be used to validate analytic and computational models of the noise sources and associated aerodynamic for configurations and conditions approximating flight for transport aircraft. The numerical results are used to both calibrate a widely used CFD code, CFL3D, and to obtain details of flap side-edge flow features not discernible from experimental observations. Both experimental set-ups were numerically modeled by using multiple block structured grids. Various turbulence models, grid block-interface interaction methods and grid topologies were implemented. Numerical results of both simulations are in excellent agreement with experimental measurements and flow visualization observations. The flow field in the flap-edge region was adequately resolved to discern some crucial information about the flow physics and to substantiate the merger of the two vortical structures. As a result of these investigations, airframe noise modelers have proposed various simplified models which use the results obtained from the steady-state computations as input.

Takalluk, M. A.↗

Navier-Stokes calculations and turbulence modeling in the trailing edge region of a circulation control airfoil

The accurate prediction of turbulent flows over curved surfaces in general and over the trailing edge region of circulation control airfoils in particular requires the coupled efforts of turbulence modelers, numerical analysts and experimentalists. The purpose of the research program in this area is described. Then, the influence on turbulence modeling of the flow characteristics over a typical circulation control wing is discussed. Next, the scope of this effort to study turbulence in the trailing edge region of a circulation control airfoil is presented. This is followed by a brief overview of the computation scheme, including the grid, governing equations, numerical method, boundary conditions and turbulence models applied to date. Then, examples of applications of two algebraic eddy viscosity models to the trailing edge region of a circulation control airfoil is presented. The results from the calculations is summarized, and conclusions drawn based on examples. Finally, the future directions of the program is outlined.

Viegas, John R.↗

Revolutionizing thermal Management in Next-Generation AI data centers: Challenges and breakthrough innovations

Data centers (DCs) serve as critical infrastructure for powering the growth and evolution of AI. Next-generation AI DCs present unique challenges in thermal management driven by unprecedented computational demands. This paper provides a comprehensive summary of key stakeholder perspectives on technology gaps, infrastructure requirements, test bed needs, emerging opportunities, and preliminary solutions related to thermal management for AI DCs. It establishes six strategic pillars of thermal management for next generation AI DC: reliability, deployability, efficiency, resilience, measurability, and valorization. The discussion spans a range of critical topics, including advanced cooling technologies, thermal strategies for emerging modular and edge DCs, system-level optimization and control frameworks, infrastructure planning and grid integration designs, benchmarking approaches, and pathways for waste heat recovery and reuse. The proposed research, development, and demonstration efforts are aimed at accelerating the deployment of AI DCs while ensuring energy efficiency, reliability, safety, and regulatory compliance.

Wang, Pengtao [ORNL] (ORCID:0000000214713429)↗

Physics-informed heterogeneous graph neural networks for DC blocker placement

The threat of geomagnetic disturbances (GMDs) to the reliable operation of the bulk energy system has spurred the development of effective strategies for mitigating their impacts. One such approach involves placing transformer neutral blocking devices, which interrupt the path of geomagnetically induced currents (GICs) to limit their impact. The high cost of these devices and the sparsity of transformers that experience high GICs during GMD events, however, calls for a sparse placement strategy that involves high computational cost. To address this challenge, we developed a physics-informed heterogeneous graph neural network (PIHGNN) for solving the graph-based dc-blocker placement problem. Our approach combines a heterogeneous graph neural network (HGNN) with a physics-informed neural network (PINN) to capture the diverse types of nodes and edges in ac/dc networks and incorporates the physical laws of the power grid. We train the PIHGNN model using a surrogate power flow model and validate it using case studies. Results demonstrate that PIHGNN can effectively and efficiently support the deployment of GIC dc-current blockers, ensuring the continued supply of electricity to meet societal demands. Furthermore, our approach has the potential to contribute to the development of more reliable and resilient power grids capable of withstanding the growing threat that GMDs pose.

24 POWER TRANSMISSION AND DISTRIBUTION↗