Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large grids”

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 775 records · Page 43

Unstructured LES of Reacting Multiphase Flows in Realistic Gas Turbine Combustors

As part of the Accelerated Strategic Computing Initiative (ASCI) program, an accurate and robust simulation tool is being developed to perform high-fidelity LES studies of multiphase, multiscale turbulent reacting flows in aircraft gas turbine combustor configurations using hybrid unstructured grids. In the combustor, pressurized gas from the upstream compressor is reacted with atomized liquid fuel to produce the combustion products that drive the downstream turbine. The Large Eddy Simulation (LES) approach is used to simulate the combustor because of its demonstrated superiority over RANS in predicting turbulent mixing, which is central to combustion. This paper summarizes the accomplishments of the combustor group over the past year, concentrating mainly on the two major milestones achieved this year: 1) Large scale simulation: A major rewrite and redesign of the flagship unstructured LES code has allowed the group to perform large eddy simulations of the complete combustor geometry (all 18 injectors) with over 100 million control volumes; 2) Multi-physics simulation in complex geometry: The first multi-physics simulations including fuel spray breakup, coalescence, evaporation, and combustion are now being performed in a single periodic sector (1/18th) of an actual Pratt & Whitney combustor geometry.

Ham, Frank↗

Simulating Transonic Buffet Aerodynamics for the Boeing Transonic Truss-Braced Wing Aircraft

The flow past the Boeing Transonic Truss-Braced Wing (TTBW) aircraft is simulated with a Hybrid Reynolds-Averaged Navier Stokes Large-Eddy Simulations (HRLES) turbulence modeling approach to study transonic buffet onset and evaluate the predictive capability of the numerical approach. All simulations used structured overset curvilinear grids for the free-air configuration, with all simulations performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Aerodynamic loads and surface pressure obtained from the HRLES solutions are compared to results from unsteady RANS (URANS) and experiments conducted at NASA Ames 11-by 11-Foot Transonic Wind Tunnel. Overall good agreement is obtained in the predicted loads and surface pressure. Unsteady pressure data from HRLES and URANS simulations is utilized to compute power spectral density (PSD) to predict transonic buffet onset and compared to the Kulite data obtained from the experiment. Comparisons of the PSD spectrum also show reasonable agreement with the wind tunnel data, with HRLES and URANS predicting buffet onset at a slightly earlier angle-of-attack than the experiment.

AATT↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

A Wall-Modeled LES Perspective for the High Lift Common Research Model Using LAVA

A new immersed boundary Wall-Modelled Large Eddy Simulation (WMLES) formulation is developed to study high-lift aerodynamics on the NASA High-Lift Common Research Model (HL-CRM). A sequence of Cartesian Octree grids with sizes ranging from 100 Million through 2.02 Billion grid points is utilized to systematically assess grid-sensitivity and convergence for the in-tunnel (QinetiQ) configuration of the model, and remarkable agreement between the immersed boundary and the curvilinear body-aligned WMLES formulations is reported on grids with comparable resolutions. In the free-air configuration of the model, consistent predictions between the Curvilinear Overset and the Cartesian Octree formulations are reported for angles of attack up to C(L,max) at a=19.57. However, some differences in the onset of stall are seen between the two methods for a>20°: while the curvilinear WMLES experiences wing-root separation with increasing angle of attack (Topology A), the Cartesian Octree formulation shows a different flow topology characterized by boundary layer weakness on the main element, emanating from the pylon-wing attachment (Topology B). In order to obtain further insight into the two-distinct topologies, carefully designed numerical experiments to isolate effects of the model standoff and the tunnel wall-boundary layers are conducted using the immersed boundary WMLES formulation. The increased incidence angle-of-attack on the inboard portion of the wing due to the standoff is shown to be sufficient for triggering a switch from Topology B to Topology A in Cartesian WMLES. The role of the floor boundary layer is further examined in detail by identification of additional corner-flow vorticity generated by the viscous juncture flow interactions between the floor boundary layer and the standoff leading to formation of a strong coherent and persistent vortex on the belly-side of the fuselage. The intensity of this vortex is shown to increase with the thickness of the floor boundary layer. A further increase in the incidence angle of attack near the leading-edge strake caused by the presence of this belly-side vortex is quantified for two-distinct floor boundary layers. Both of the floor boundary layers considered result in the onset of large scale wing-root separation at a=21.47in non-confined (free-air) configurations.

TTT↗

Numerical Investigation of C L,max Prediction on the NASA High-Lift Common Research Model

An assessment of a Hybrid RANS/LES (HRLES) approach for C L,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both free air and wind tunnel configurations of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES is shown to address some of the known shortcomings with RANS and prevent inboard and outboard flow separation particularly in the region of C L,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly underperformed RANS. For the in-tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.

TTT↗

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗

Using Multi-scale Modeling System to Study the Interactions between Clouds, Precipitation, Aerosols, Radiation and Land Surface

Numerical cloud models, which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Because cloud-scale dynamics are treated explicitly, uncertainties stemming from convection that have to be parameterized in (hydrostatic) large-scale models are obviated, or at least mitigated, in cloud models. Global models will use the non-hydrostatic framework when their horizontal resolution becomes about 10 kilometers, the theoretical limit for the hydrostatic approximation. This juncture will be reached one to two decades from now. Over the past generation, voluminous datasets on atmospheric convection have been accumulated from radar, instrumented aircraft, satellites, and rawinsonde measurements in field campaigns, enabling the detailed evaluation of models. Improved numerical methods have resulted in more accurate and efficient dynamical cores in models. Improvements have been made in the parameterizations of microphysical processes, radiation, boundary-layer effects, and turbulence; however, microphysical parameterizations remain a major source of uncertainty in all classes of atmospheric models. In recent years, exponentially increasing computer power has extended cloud-resolving-model integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-kilometer scales are resolved in horizontal domains as large as 10,000 kilometers in two-dimensions, and 1,000 x 1,000 square kilometers in three-dimensions. Cloud models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale model can be run in grid size similar to cloud resolving model through nesting technique. A review of developments, improvements and applications of cloud models (GCE and WRF) at Goddard will be presented in this talk. In particular, a new approach to using multi-scale modeling system to study the interactions between clouds, precipitation, aerosols and land will be presented.

Tao, Wei-Kuo↗

Using Multi-scale Modeling System to Study the Interactions between Clouds, Precipitation, Aerosols, Radiation and Land Surface

Numerical cloud models, which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Because cloud-scale dynamics are treated explicitly, uncertainties stemming from convection that have to be parameterized in (hydrostatic) large-scale models are obviated, or at least mitigated, in cloud models. Global models will use the non-hydrostatic framework when their horizontal resolution becomes about 10 kilometers, the theoretical limit for the hydrostatic approximation. This juncture will be reached one to two decades from now. Over the past generation, voluminous datasets on atmospheric convection have been accumulated from radar, instrumented aircraft, satellites, and rawinsonde measurements in field campaigns, enabling the detailed evaluation of models. Improved numerical methods have resulted in more accurate and efficient dynamical cores in models. Improvements have been made in the parameterizations of microphysical processes, radiation, boundary-layer effects, and turbulence; however, microphysical parameterizations remain a major source of uncertainty in all classes of atmospheric models. In recent years, exponentially increasing computer power has extended cloud-resolving-model integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial-lkm scales are resolved in horizontal domains as large as 10,000 kilometers in two-dimensions, and 1,000 x 1,000 square kilometers in three-dimensions. Cloud models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale model can be run in grid size similar to cloud resolving model through nesting technique. A review of developments, improvements and applications of cloud models (GCE and WRF) at Goddard wlll be is presented in this talk. In particular, a new approach to using multi-scale modeling system to study the interactions between clouds, precipitation, aerosols and land will be presented.

Tao, Wei-Kuo↗

Using Multi-scale Modeling System to Study the Interactions between Clouds, Precipitation, Aerosols, Radiation and Land Surface

Numerical cloud models, which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Because cloud-scale dynamics are treated explicitly, uncertainties stemming from convection that have to be parameterized in (hydrostatic) large-scale models are obviated, or at least mitigated, in cloud models. Global models will use the non-hydrostatic framework when their horizontal resolution becomes about 10 km, the theoretical limit for the hydrostatic approximation. This juncture will be reached one to two decades from now. Over the past generation, voluminous datasets on atmospheric convection have been accumulated from radar, instrumented aircraft, satellites, and rawinsonde measurements in field campaigns, enabling the detailed evaluation of models. Improved numerical methods have resulted in more accurate and efficient dynamical cores in models. Improvements have been made in the parameterizations of microphysical processes, radiation, boundary-layer effects, and turbulence; however, microphysical parameterizations remain a major source of uncertainty in all classes of atmospheric models. In recent years, exponentially increasing computer power has extended cloud-resolving-model integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-km scales are resolved in horizontal domains as large as 10,000 km in two-dimensions, and 1,000 x 1,000 sq km in three-dimensions. Cloud models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that NWP and mesoscale model can be run in grid size similar to cloud resolving model through nesting technique. A review of developments, improvements and applications of cloud models (GCE and WRF) at Goddard will be presented in this talk. In particular, a new approach to using multi-scale modeling system to study the interactions between clouds, precipitation, aerosols and land will be presented.

Tao, Wei-Kuo↗

Results of a utility survey of the status of large wind turbine development

Wind energy conversion systems were surveyed from a utility viewpoint to establish the state of the art with regard to: (1) availability of the type of machines; (2) quality of power generation; (3) suitability for electrical grid; (4) reliability; and (5) economics. Of the 23 designs discussed, 7 have vertical axis wind turbines, 9 have upwind horizontal axis turbines, and 7 have downwind horizontal axis turbines.

Watts, A.↗

Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques

Refinements to the modeling procedure developed to compute the magnetic fields at satellite orbit due to current distributions in the ionosphere and magnetosphere are described. The modeling technique utilizes a linear current element representation of the large scale space current system. A model polar current system is presented and magnetic field perturbations resulting from this system are computed along two hypothetical satellite orbits.

Klumpar, D. M.↗

Measured effects of wind turbine generation at the Block Island Power Company

Data measurements made on the NASA MOD-OA 200-kw wind-turbine generator (WTG) installed on a utility grid form the basis for an overall performance analysis. Fuel displacement/-savings, dynamic interactions, and WTG excitation (reactive-power) control effects are studied. Continuous recording of a large number of electrical and mechanical variables on FM magnetic tape permit evaluation and correlation of phenomena over a bandwidth of at least 20 Hz. Because the wind-power penetration reached peaks of 60 percent, the impact of wind fluctuation and wind-turbine/diesel-utility interaction is evaluated in a worst-case scenario. The speed-governor dynamics of the diesel units exhibited an underdamped response, and the utility operation procedures were not altered to optimize overall WTG/utility performance. Primary findings over the data collection period are: a calculated 6.7-percent reduction in fuel consumption while generating 11 percent of the total electrical energy; acceptable system voltage and frequency fluctuations with WTG connected; and applicability of WTG excitation schemes using voltage, power, or VARS as the controlled variable.

Wilreker, V. F.↗

Numerical simulation of high-incidence flow over the F-18 fuselage forebody

As part of the NASA High Alpha Technology Program, fine-grid Navier-Stokes solutions have been obtained for flow over the fuselage forebody and wing leading-edge extension of the F/A-18 High Alpha Research Vehicle at large incidence. The resulting flows are complex and exhibit cross-flow separation from the sides of the forebody and from the leading-edge extension. A well-defined vortex pattern is observed in the leeward-side flow. Results obtained for laminar flow show good agreement with flow visualizations obtained in ground-based experiments. Further, turbulent flows computed at high-Reynolds-number flight-test conditions show good agreement with surface and off-surface visualizations obtained in flight.

Schiff, Lewis B.↗

Asymmetric supersonic flow around cones with noncircular sections

The unsteady, compressible, single and double thin-layer, Navier-Stokes equations are used to solve for steady and unsteady, asymmetric, supersonic flow around pointed bodies with noncircular sections at high incidence and zero side slip. The equations are solved by using an implicit, upwind, flux-difference splitting finite-volume scheme. Since the flow is locally conical, the solutions are presented on a cross-flow plane at the axial station of unit. The grid is generated by using a modified Joukowski transformation. The computational applications cover noncircular sections with elliptic and diamond shapes. Unsteady asymmetric vortex shedding has been captured at large angles of attack. It is shown that for the same flow conditions and same cross-section fineness ratio, the diamond-section cones develop less flow asymmetry than the elliptic-section cones. Passive control of flow asymmetry was also demonstrated computationally.

Kandil, Osama A.↗

Numerical simulation of rarefied flow through a slit. I - Direct simulation Monte Carlo results

The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.

Wadsworth, D. C.↗

Runtime volume visualization for parallel CFD

This paper discusses some aspects of design of a data distributed, massively parallel volume rendering library for runtime visualization of parallel computational fluid dynamics simulations in a message-passing environment. Unlike the traditional scheme in which visualization is a postprocessing step, the rendering is done in place on each node processor. Computational scientists who run large-scale simulations on a massively parallel computer can thus perform interactive monitoring of their simulations. The current library provides an interface to handle volume data on rectilinear grids. The same design principles can be generalized to handle other types of grids. For demonstration, we run a parallel Navier-Stokes solver making use of this rendering library on the Intel Paragon XP/S. The interactive visual response achieved is found to be very useful. Performance studies show that the parallel rendering process is scalable with the size of the simulation as well as with the parallel computer.

Ma, Kwan-Liu↗

Discovery of a New Cassini-Size Basin on Mars from MOLA Topographic Data

MOLA profile data collected during the Science Phasing Operations have revealed a previously unknown, 450 km wide, 2 km deep basin on Mars centered at 30N, 312W near the Phison Rupes. This structure, as large as and somewhat deeper than the very obvious Cassini impact basin located 1000 km to the SW, is not apparent in the existing good quality Viking imagery. Gridded MOLA data show the basin as a closed depression, but elevation contours show only weak correlation with what little structure exists in the area and with mapped geologic units. From analysis of slope breaks readily visible in two MOLA profiles we suggest this is at least a three-ring basin. Portions of rings from concentric fits to slope breaks align with some of the linear ridge-like structures of the Phison Rupes, and outline a region of lower crater density and smoother inter-crater plains. The discovery of such a pronounced topographic depression which lacks obvious visible structure suggests more such previously unknown features may exist on Mars, and that MOLA data may be useful in finding them.

Frey, H.↗