Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “zonal flows”

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 433 records · Page 24

A simulation of rotor-stator interaction using the Euler equations and patched grids

An unsteady Euler code to study rotor-stator interaction problem was developed. The code uses patched grids that move relative to each other to simulate the motion of the rotor airfoils with respect to the stator airfoils. The Osher integration scheme is used in conjunction with an implicit relaxation approach. The scheme is second order accurate in space and time, and is also TVD in each spatial direction. The numerical results were found to be periodic in time, thus demonstrating the capability of the integration and zonal schemes in simulating periodic time dependent flow. The pressure contours obtained are almost oscillation free because of the TVD nature of the scheme. A new procedure was developed to simulate flows about bodies that move relative to each other. This capability should prove to be very useful in the areas of rotor-stator interaction, propeller-nacelle interaction, and helicopter rotor-fuselage interaction.

Rai, M. M.↗

A composite structured/unstructured-mesh Euler method for complex airfoil shapes

A general two-dimensional Euler zonal method has been developed for computing flows about complex airfoil geometries such as multielement and iced airfoils. The method utilizes a composite structured and unstructured grid generated using conformal mapping and Delaunay triangulation, respectively. The finite-volume Euler method is then modified to couple solutions in the zones with structured and unstructured grids. Solutions about an iced airfoil and a multielement airfoil are given as examples of applications of the scheme.

Hefazi, H.↗

Angle-of-attack validation of a new zonal CFD method for airfoil simulations

The angle-of-attack validation of a new concept suggested by Summa (1990) for coupling potential and viscous flow methods has been investigated for two-dimensional airfoil simulations. The fully coupled potential/Navier-Stokes code, ZAP2D (Zonal Aerodynamics Program 2D), has been used to compute the flow field around an NACA 0012 airfoil for a range of angles of attack up to stall at a Mach number of 0.3 and a Reynolds number of 3 million. ZAP2D calculation for various domain sizes from 25 to 0.12 chord lengths are compared with the ARC2D large domain solution as well as with experimental data.

Yoo, Sungyul↗

Computational aspects of zonal algorithms for solving the compressible Navier-Stokes equations in three dimensions

Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The flow field immediately adjacent to the wing surface is resolved with fine grid zones and solved using a Navier-Stokes algorithm. Flow field regions removed from the wing are resolved with less finely clustered grid zones and are solved with an Euler algorithm. Computational issues associated with this zonal approach, including data base management aspects, are discussed. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects. Additional cases with significant shock induced separation on the upper wing surface are also presented.

Holst, T. L.↗

MPD thruster performance - Propellant distribution and species effects

Previous research has shown that MPD thruster performance, as measured by thrust, specific impulse, and efficiency, is limited by a so-called 'onset' current level. Experimental results are presented which show that the onset current level varies with the radial position of mass injection. Specifically, mass injection into either of two zones results in drastically different performance for both argon and xenon propellants. This difference in performance is evidenced by the onset current levels, the onset V-J signatures, and the overall V-J characteristics. The behavior of the thruster with the total mass flow split between the two zones follows a transition between the zonal behavioral limits. The behavior of propellant mixtures also follows a transition between the behavior of the individual propellant species. The flow region between the two zones is characterized by a particle number density discontinuity, a large axial current density, and a large electron Hall parameter, suggesting a possible correlation between observed onset phenomena and a shear-flow discontinuity located in the region between the two zones.

Kelly, A. J.↗

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗

Transmission Capacity Allocation in Zonal Electricity Markets

In this work, we propose a novel framework for modeling zonal electricity markets, based on projecting the constraints of the nodal network onto the space of the zonal aggregation of the network. The framework avoids circular definitions and discretionary parameters, which are recurrent in the implementation and study of zonal markets. Using this framework, we model and analyze two zonal market designs currently present in Europe: flow-based market coupling (FBMC) and available-transfer-capacity market coupling (ATCMC). We develop cutting-plane algorithms for simulating FBMC and ATCMC while accounting for the robustness of imports/exports to single element failures, and we conduct numerical simulations of FBMC and ATCMC on a realistic instance of the Central Western European system under the equivalent of 100 years of operating conditions. We find that FBMC and ATCMC are unable to anticipate the congestion of branches interconnecting zones and branches within zones and that both zonal designs achieve similar overall cost efficiencies (0.01% difference), whereas a nodal market design largely outperforms both of them (5.09% better than FBMC). These findings raise the question of whether it is worth it for more European countries to switch from ATCMC to FBMC, instead of advancing directly toward a nodal market design.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A solution-adaptive hybrid-grid method for the unsteady analysis of turbomachinery

A solution-adaptive method for the time-accurate analysis of two-dimensional flows in turbomachinery is described. The method employs a hybrid structured-unstructured zonal grid topology in conjunction with appropriate modeling equations and solution techniques in each zone. The viscous flow region in the immediate vicinity of the airfoils is resolved on structured O-type grids while the rest of the domain is discretized using an unstructured mesh of triangular cells. Implicit, third-order accurate, upwind solutions of the Navier-Stokes equations are obtained in the inner regions. In the outer regions, the Euler equations are solved using an explicit upwind scheme that incorporates a second-order reconstruction procedure. An efficient and robust grid adaptation strategy, including both grid refinement and coarsening capabilities, is developed for the unstructured grid regions. Grid adaptation is also employed to facilitate information transfer at the interfaces between unstructured grids in relative motion. Results for grid adaptation to various features pertinent to turbomachinery flows are presented. Good comparisons between the present results and experimental measurements and earlier structured-grid results are obtained.

Mathur, Sanjay R.↗

Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint

Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.

atmospheric boundary layer↗

Navier-Stokes solutions for flows related to store separation

The objective is developing CFD capabilities to obtain solutions for viscous flows about generic configurations of internally and externally carried stores. The emphasis is placed on the supersonic flow regime with extensions being made to the transonic regime. The project is broken into four steps: (1) Cavity flows for internal carriage configurations; (2) High angle of attack flows, which may be experienced during the separation of the stores: (3) Flows about a body near a flat plate for external carriage configurations; and (4) Flows about a body inside or in the proximity of a cavity. Three-dimensional unsteady cavity flow solutions are obtained by an explicit, MacCormack algorithm, EMCAV3, for open, close, and transitional cavities. High angle of attack flows past cylinders are solved by an implicit, upwind algorithm. All the results compare favorably with the experimental data. For flows about multiple body configurations, the Chimera embedding scheme is modified for finite-volume and multigrid algorithms, MaGGiE. Then a finite volume, implicit, upwind, multigrid Navier-Stokes solver which uses on overlapped/embedded and zonal grids, VUMXZ3, is developed from the CFL3D code. Supersonic flows past a cylinder near a flat plate are computed using this code. The results are compared with the experimental data. Currently the VUMXZ3 code is being modified to accomplish step 4 of this project. Wind tunnel experiments are also being conducted for validation purposes.

Baysal, Oktay↗

The role of forced planetary waves in the annual cycle of the zonal mean circulation of the middle atmosphere

A severely truncated semispectral numerical model is used to simulate the annual cycle of the zonally averaged circulation in the middle atmosphere (16-96 km). The model includes only a single zonal harmonic wave component which interacts with the mean flow; the circulation is driven by diabatic heating and by a specified perturbation in the topography of the lower boundary, which is taken to be the 100 mb surface. A comparison of the annual cycle simulated by this model with the results of an analogous two-dimensional model indicates that planetary waves have relatively little influence on the zonal mean temperature profiles and on the solstice mean zonal winds at high latitudes. The primary effects of the forced waves are in decelerating the mean winds at low latitudes in the winter hemisphere to produce a region of weak westerlies, and in generating final warmings at the spring equinoxes.

Holton, J. R.↗

Tip vortices of wings in subsonic and transonic flow: A numerical simulation

Thin layer Navier-Stokes and Euler equations are numerically solved using a multi-block zonal approach to simulate the formation and roll up of tip vortices of wings in subsonic and transonic flows. Several wing planforms were considered to examine the influence of tip-cap shape, planform geometry and free stream Mach number on the formation process. A good definition of the formation and qualitative roll up of tip vortices was achieved.

Srinivasan, G. R.↗

VHF radar measurements in the summer polar mesosphere

Measurements in the mesosphere over Andoya/Norway (69 N, 16 E) were carried out using the mobile SOUSY-VHF radar with an extended beam configuration during the MAC/SINE campaign in summer 1987. First results of a 48 h and a 3 h observational period for heights between about 83 and 91 km are presented. Zonal mean winds are characterized by a strong westward flow of up to 50/ms, whereas the equatorward directed meridional component is weaker. The dominating semidiurnal tide has amplitudes up to 30/ms and a vertical wavelength of about 55 km. The diurnal tide is less pronounced. The total upward flux of horizontal momentum takes values of -2 sq m/sq s near 84 km and increases with increasing height, reaching a maximum value of 22 sq m/sqs for both the zonal and meridional components. However, measurements of the horizontal isotropy of the wave field suggest significant anisotropy. The major contribution to the momentum flux is from the 10 min to 1 h period range below about 87 km, and from the 1 to 6 h period range above this height.

Ruester, R.↗

Momentum flux measurements: Techniques and needs, part 4.5A

The vertical flux of horizontal momentum by internal gravity waves is now recognized to play a significant role in the large-scale circulation and thermal structure of the middle atmosphere. This is because a divergence of momentum flux due to wave dissipation results in an acceleration of the local mean flow towards the phase speed of the gravity wave. Such mean flow acceleration are required to offset the large zonal accelerations driven by Coriolis torques acting on the diabatic meridional circulation. Techniques and observations regarding the momentum flux distribution in the middle atmosphere are discussed.

Fritts, D. C.↗

Some aspects of stratospheric dynamics

Observations of the stratosphere are analyzed. The behavior of the wind and temperature in the stratosphere is compared to that of the troposphere. Techniques for observing the circulation and temperatures of the stratosphere and the basic circulation and temperature characteristics are described. The propagation of tropospheric planetary waves to the stratosphere and their interaction with the mean flow are studied. The stability of the zonal-mean state in the stratosphere is investigated by calculating the meridional gradient of zonal-mean quasi-geostrophic potential vorticity.

Hartmann, D. L.↗

Changes in Jupiter's Zonal Wind Profile Preceding and During the Juno Mission

We present five epochs of WFC3 HST Jupiter observations taken between 2009-2016 and extract global zonal wind profiles for each epoch. Jupiter's zonal wind field is globally stable throughout these years, but significant variations in certain latitude regions persist. We find that the largest uncertainties in the wind field are due to vortices or hot-spots, and show residual maps which identify the strongest vortex flows. The strongest year-to-year variation in the zonal wind profiles is the 24 deg N jet peak. Numerous plume outbreaks have been observed in the Northern Temperate Belt and are associated with decreases in the zonal velocity and brightness. We show that the 24 deg N jet peak velocity and brightness decreased in 2012 and again in late 2016, following outbreaks during these years. Our February 2016 zonal wind profile was the last highly spatially resolved measurement prior to Juno s first science observations. The final 2016 data were taken in conjunction with Juno's perijove 3 pass on 11 December 2016, and show the zonal wind profile following the plume outbreak at 24 deg N in October 2016.

Jupiter↗

Characterization of Unsteady Flow Structures Near Landing-Edge Slat: 2D Computations - Part 2

In our previous computational studies of a generic high-lift configuration, quasi-laminar (as opposed to fully turbulent) treatment of the slat cove region proved to be an effective approach for capturing the unsteady dynamics of the cove flow field. Combined with acoustic propagation via Ffowes Williams and Hawkings formulation, the quasi-laminar simulations captured some important features of the slat cove noise measured with microphone array techniques. However. a direct assessment of the computed cove flow field was not feasible due to the unavailability of off-surface flow measurements. To remedy this shortcoming, we have undertaken a combined experiment and computational study aimed at characterizing the flow structures and fluid mechanical processes within the slat cove region. Part I of this paper outlines the experimental aspects of this investigation focused on the 30P30N high-lift configuration; the present paper describes the accompanying computational results including a comparison between computation and experiment at various angles of attack. Even through predictions of the time-averaged flow field agree well with the measured data, the study indicates the need for further refinement of the zonal turbulence approach in order to capture the full dynamics of the cove's fluctuating flow field.

Khorrami, Mehdi↗

Transonic Navier-Stokes wing solutions using a zonal approach. Part 2: High angle-of-attack simulation

A computer code is under development whereby the thin-layer Reynolds-averaged Navier-Stokes equations are to be applied to realistic fighter aircraft configurations. This transonic Navier-Stokes code (TNS) utilizes a zonal approach in order to treat complex geometries and satisfy in-core computer memory constraints. The zonal approach was applied to isolated wing geometries in order to facilitate code development. The TNS finite difference algorithm, zonal methodology, and code validation with experimental data is addressed. Also addressed are some numerical issues such as code robustness, efficiency, and accuracy at high angles of attack. Special free-stream-preserving metrics proved an effective way to treat H-mesh singularities over a large range of severe flow conditions, including strong leading edge flow gradients, massive shock induced separation, and stall. Furthermore, lift and drag coefficients were computed for a wing up through CLmax. Numerical oil flow patterns and particle trajectories are presented both for subcritical and transonic flow. These flow simulations are rich with complex separated flow physics and demonstrate the efficiency and robustness of the zonal approach.

Chaderjian, Neal M.↗