Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution grid”

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 757 records · Page 42

Heat transfer predictions for two turbine nozzle geometries at high Reynolds and Mach numbers

Predictions of turbine vane and endwall heat transfer and pressure distributions are compared with experimental measurements for two vane geometries. The differences in geometries were due to differences in the hub profile, and both geometries were derived from the design of a high rim speed turbine (HRST). The experiments were conducted in the Isentropic Light Piston Facility (ILPF) at Pyestock at a Reynolds number of 5.3 x 10(exp 6), a Mach number of 1.2, and a wall-to-gas temperature ratio of 0.66. Predictions are given for two different steady-state three-dimensional Navier-Stokes computational analyses. C-type meshes were used, and algebraic models were employed to calculate the turbulent eddy viscosity. The effects of different turbulence modeling assumptions on the predicted results are examined. Comparisons are also given between predicted and measured total pressure distributions behind the vane. The combination of realistic engine geometries and flow conditions proved to be quite demanding in terms of the convergence of the CFD solutions. An appropriate method of grid generation, which resulted in consistently converged CFD solutions, was identified.

Boyle, R. J.↗

A Study of Heavy Precipitation Events in Taiwan During 10-13 August, 1994: Mesoscale Model Simulations - Part 2

Several heavy precipitation episodes occurred over Taiwan from August 10 to 13, 1994. Precipitation patterns and characteristics are quite different between the precipitation events that occurred from August 10 and I I and from August 12 and 13. In Part I (Chen et al. 2001), the environmental situation and precipitation characteristics are analyzed using the EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. In this study (Part II), the Penn State/NCAR Mesoscale Model (MM5) is used to study the precipitation characteristics of these heavy precipitation events. Various physical processes (schemes) developed at NASA Goddard Space Flight Center (i.e., cloud microphysics scheme, radiative transfer model, and land-soil-vegetation surface model) have recently implemented into the MM5. These physical packages are described in the paper, Two way interactive nested grids are used with horizontal resolutions of 45, 15 and 5 km. The model results indicated that Cloud physics, land surface and radiation processes generally do not change the location (horizontal distribution) of heavy precipitation. The Goddard 3-class ice scheme produced more rainfall than the 2-class scheme. The Goddard multi-broad-band radiative transfer model reduced precipitation compared to a one-broad band (emissivity) radiation model. The Goddard land-soil-vegetation surface model also reduce the rainfall compared to a simple surface model in which the surface temperature is computed from a Surface energy budget following the "force-re store" method. However, model runs including all Goddard physical processes enhanced precipitation significantly for both cases. The results from these runs are in better agreement with observations. Despite improved simulations using different physical schemes, there are still some deficiencies in the model simulations. Some potential problems are discussed. Sensitivity tests (removing either terrain or radiative processes) are performed to identify the physical processes that determine the precipitation patterns and characteristics for heavy rainfall events. These sensitivity tests indicated that terrain can play a major role in determining the exact location for both precipitation events. The terrain can also play a major role in determining the intensity of precipitation for both events. However, it has a large impact on one event but a smaller one on the other. The radiative processes are also important for determining, the precipitation patterns for one case but. not the other. The radiative processes can also effect the total rainfall for both cases to different extents.

Tao, Wei Kuo↗

Certification of CFD heat transfer software for turbine blade analysis

Accurate modeling of heat transfer effects is a critical component of the Turbine Branch of the Turbomachinery and Propulsion Systems Division. Being able to adequately predict and model heat flux, coolant flows, and peak temperatures are necessary for the analysis of high pressure turbine blades. To that end, the primary goal of my internship this summer will be to certify the reliability of the CFD program GlennHT for the purpose of turbine blade heat transfer analysis. GlennHT is currently in use by the engineers in the Turbine Branch who use the FORTRAN 77 version of the code for analysis. The program, however, has been updated to a FORTRAN 90 version which is more robust than the older code. In order for the new code to be distributed for use, its reliability must first be certified. Over the course of my internship I will create and run test cases using the FORTRAN 90 version of GlennHT and compare the results to older cases which are known to be accurate, If the results of the new code match those of the sample cases then the newer version will be one step closer to certification for distribution. In order to complete these it will first be necessary to become familiar with operating a number of other programs. Among them are GridPro, which is used to create a grid mesh around a blade geometry, and FieldView, whose purpose is to graphically display the results from the GlennHT program. Once enough familiarity is established with these programs to render them useful, then the work of creating and running test scenarios will begin. The work is additionally complicated by a transition in computer hardware. Most of the working computers in the Turbine Branch are Silicon Graphics machines, which will soon be replaced by LINUX PC's. My project is one of the first to make use the new PC's. The change in system architecture however, has created several software related issues which have greatly increased the time and effort investments required by the project.Although complications with the project continue to arise, it is expected that the goal of my internship can still be achieved within the remaining time period. Critical steps have been achieved and test scenarios can now be designed and run. At the completion of my internship, the FORTRAN 90 version of GlennHT should be well on its way to certification.

Jordan, William A.↗

Particle Morphology and Size Results from the Smoke Aerosol Measurement Experiment-2

Results are presented from the Reflight of the Smoke Aerosol Measurement Experiment (SAME-2) which was conducted during Expedition 24 (July-September 2010). The reflight experiment built upon the results of the original flight during Expedition 15 by adding diagnostic measurements and expanding the test matrix. Five different materials representative of those found in spacecraft (Teflon, Kapton, cotton, silicone rubber and Pyrell) were heated to temperatures below the ignition point with conditions controlled to provide repeatable sample surface temperatures and air flow. The air flow past the sample during the heating period ranged from quiescent to 8 cm/s. The smoke was initially collected in an aging chamber to simulate the transport time from the smoke source to the detector. This effective transport time was varied by holding the smoke in the aging chamber for times ranging from 11 to 1800 s. Smoke particle samples were collected on Transmission Electron Microscope (TEM) grids for post-flight analysis. The TEM grids were analyzed to observe the particle morphology and size parameters. The diagnostics included a prototype two-moment smoke detector and three different measures of moments of the particle size distribution. These moment diagnostics were used to determine the particle number concentration (zeroth moment), the diameter concentration (first moment), and the mass concentration (third moment). These statistics were combined to determine the diameter of average mass and the count mean diameter and, by assuming a log-normal distribution, the geometric mean diameter and the geometric standard deviations can also be calculated. Overall the majority of the average smoke particle sizes were found to be in the 200 nm to 400 nm range with the quiescent cases producing some cases with substantially larger particles.

Urban, David L.↗

Cloud Regimes as a Tool for Systematic Study of Various Aerosol-Cloud-Precipitation Interactions

Systematic changes of clouds and precipitation are notoriously difficult to ascribe to aerosols. This presentation will showcase yet one more attempt to at least credibly detect the signal of aerosol-cloud-precipitation interactions. We surmise that the concept of cloud regimes (CRs) is appropriate to conduct such an investigation. Previous studies focused on what we call here dynamical CRs, and while we continue to adopt those too for our analysis, we have found that a different way of organizing cloud systems, namely via microphysical regimes is also promising. Our analysis relies on MODIS Collection 6 Level-3 data for clouds and aerosols, and TRMM-TMPA data for precipitation. The regimes are derived by applying clustering analysis on MODIS joint histograms, and once each grid cell is assigned a regime, aerosol and precipitation data can be spatiotemporally matched and composited by regime. The composites of various cloud and precipitation variables for high (upper quartile of distribution) and low (lower quartile) aerosol loadings can then be contrasted. We seek evidence of aerosol effects both in regimes with large fractions of deep ice-rich clouds, as well as regimes where low liquid phase clouds dominate. Signals can be seen, especially when the analysis is broken by land-ocean and when additional filters are applied, but there are of course caveats which will be discussed.

cloud regimes↗

Simulation of transonic viscous flow over a fighter-like configuration including inlet

The simulation of transonic viscous flow over a modified F-16A including inlet is presented. A zonal approach is utilized which allows appropriate clustering suitable for viscous calculations on all solid surfaces. Computational efficiency is enhanced by solving the thin-layer Navier-Stokes equations in viscous zones adjacent to the aircraft and the Euler equations on those coarse zones away from the aircraft. The flow conditions for this transonic case are freestream Mach number of 0.9, angle of attack of 4.12 deg, and a Reynolds number (based on root chord) of 4.5 million. A total of 19 zones are utilized yielding a total of 350,000 grid points. This case required about 3000 iterations to reduce the residual by three orders, which takes about 10 hr of CPU time on the Cray X-MP/48 computer. Pressure distributions on the wing and on cross sections through the inlet region compare favorably with the experimental data for this transonic case.

Flores, Jolen↗

Steady state cooling flow models for normal elliptical galaxies

An extensive grid of steady state cooling flow models for the hot gas in elliptical galaxies is calculated. The effects of allowing inflow of surrounding gas and varying the supernova heating rate, the mass distribution of the galaxy, the rate of stellar mass loss, and the outer boundary pressure are studied. The X-ray luminosities, spectra, and surface brightness profiles of the galaxies are computed. In all of the models, the gas density varies approximately as r exp -3/2. In the models with heavy halos, the X-ray and optical surface brightnesses are nearly proportional to one another within individual galaxies, and the temperature profiles are roughly isothermal.

Sarazin, Craig L.↗

Head-on parallel blade-vortex interaction

An experimental and computational study was carried out to investigate the parallel head-on blade-vortex interaction (BVI) and its noise generation mechanism. A shock tube, with an enlarged test section, was used to generate a compressible starting vortex which interacted with a target airfoil. The dual-pulsed holographic interferometry (DPHI) technique and airfoil surface pressure measurements were employed to obtain quantitative flow data during the BVI. A thin-layer Navier-Stokes code (BV12D), with a high-order upwind-biased scheme and a multizonal grid, was also used to simulate numerically the phenomena occurring in the head-on BVI. The detailed structure of a convecting vortex was studied through independent measurements of density and pressure distributions across the vortex center. Results indicate that, in a strong head-on BVI, the opposite pressure peaks are generated on both sides of the leading edge as the vortex approaches. Then, as soon as the vortex passes by the leading edge, the high-pressure peak suddenly moves toward the low-peak-reducing in magnitude as it moves--simultaneously giving rise to the initial sound wave. In both experiment and computation, it is shown that the viscous effect plays a significant role in head-on BVIs.

Lee, Soogab↗

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗

FIDDLE: A Computer Code for Finite Difference Development of Linear Elasticity in Generalized Curvilinear Coordinates

A three-dimensional numerical solver based on finite-difference solution of three-dimensional elastodynamic equations in generalized curvilinear coordinates has been developed and used to generate data such as radial and tangential stresses over various gear component geometries under rotation. The geometries considered are an annulus, a thin annular disk, and a thin solid disk. The solution is based on first principles and does not involve lumped parameter or distributed parameter systems approach. The elastodynamic equations in the velocity-stress formulation that are considered here have been used in the solution of problems of geophysics where non-rotating Cartesian grids are considered. For arbitrary geometries, these equations along with the appropriate boundary conditions have been cast in generalized curvilinear coordinates in the present study.

Kaul, Upender K.↗

Wall-Resolved Large Eddy Simulations of Transonic Shock-Induced Flow Separation

This paper reports the wall-resolved large eddy simulations of shock-induced boundary layer separation over an axisymmetric bump for a flow Mach number of 0.875 and a chord-based Reynolds number of 2.763 million. The incoming boundary layer has a momentum-thickness Reynolds number of 6600 at one and a half chord lengths upstream of the leading edge. The calculations simulate the experiment by Bachalo and Johnson (AIAA Journal, Vol. 24, No. 3, 1986), except that the tunnel walls are ignored and the simulations are performed assuming free air with as many as 24 billion grid points. The effects of domain span, grid resolution and time step on the predictions are examined. The results are found to show some sensitivity to the studied parameters. Owing to the outer boundary conditions, the predicted surface pressure distribution as well as the flow separation and reattachment locations tend to agree better with the experimental results from the larger (6 × 6 ft) tunnel than those from the smaller (2 × 2 ft) tunnel. The predicted Reynolds shear stress profiles in the separated region differ by as much as 31%from the experimental results that were only obtained in the smaller tunnel. The most accurate surface pressure distribution obtained in this study lies within the scatter of the measurements taken in the two facilities.

Uzun, Ali↗

A Conservative, Scalable, Space-Time Blade Element Rotor Model for Multi-Rotor Vehicles

The development of a parallel blade-element rotor model and its implementation into an adaptive Cartesian method is described. The unsteady version of the rotor model applies a body force to all cells contained in the swept space-time volume at each timestep and special care is taken to maintain axisymmetry on the Cartesian grid. Mesh convergence of rotor thrust and torque is obtained with around 10000 cells in the disk for the steady model. Parallelization is accomplished using OpenMP and the rotor force computation is distributed across all available nodes. Simulations of an isolated XV-15 rotor in hover show good correlation with experimental data and predictions of multi-rotor thrust variation closely match previous high fidelity simulations. The final paper will also include results from the unsteady rotor model and parallel scaling tests.

Conservative↗

Parameterization-Induced Uncertainties and Impacts of Crop Management Harmonization in a Global Gridded Crop Model Ensemble

Global gridded crop models (GGCMs) combine agronomic or plant growth models with gridded spatial input data to estimate spatially explicit crop yields and agricultural externalities at the global scale. Differences in GGCM outputs arise from the use of different biophysical models, setups, and input data. GGCM ensembles are frequently employed to bracket uncertainties in impact studies without investigating the causes of divergence in outputs. This study explores differences in maize yield estimates from five GGCMs based on the public domain field-scale model Environmental Policy Integrated Climate (EPIC) that participate in the AgMIP Global Gridded Crop Model Intercomparison initiative. Albeit using the same crop model, the GGCMs differ in model version, input data, management assumptions, parameterization, and selection of subroutines affecting crop yield estimates via cultivar distributions, soil attributes, and hydrology among others. The analyses reveal inter-annual yield variability and absolute yield levels in the EPIC-based GGCMs to be highly sensitive to soil parameterization and crop management. All GGCMs show an intermediate performance in reproducing reported yields with a higher skill if a static soil profile is assumed or sufficient plant nutrients are supplied. An in-depth comparison of setup domains for two EPIC-based GGCMs shows that GGCM performance and plant stress responses depend substantially on soil parameters and soil process parameterization, i.e. hydrology and nutrient turnover, indicating that these often neglected domains deserve more scrutiny. For agricultural impact assessments, employing a GGCM ensemble with its widely varying assumptions in setups appears the best solution for coping with uncertainties from lack of comprehensive global data on crop management, cultivar distributions and coefficients for agro-environmental processes. However, the underlying assumptions require systematic specifications to cover representative agricultural systems and environmental conditions. Furthermore, the interlinkage of parameter sensitivity from various domains such as soil parameters, nutrient turnover coefficients, and cultivar specifications highlights that global sensitivity analyses and calibration need to be performed in an integrated manner to avoid bias resulting from disregarded core model domains. Finally, relating evaluations of the EPIC-based GGCMs to a wider ensemble based on individual core models shows that structural differences outweigh in general differences in configurations of GGCMs based on the same model, and that the ensemble mean gains higher skill from the inclusion of structurally different GGCMs. Although the members of the wider ensemble herein do not consider crop-soil-management interactions, their sensitivity to nutrient supply indicates that findings for the EPIC-based sub-ensemble will likely become relevant for other GGCMs with the progressing inclusion of such processes.

Folberth, Christian↗

Winter Arctic Sea Ice Thickness from ICESat-2 Freeboards

NASA’s ICESat-2 mission was launched in September 2018 with the primary goal of monitoring our rapidly changing polar regions. The sole instrument onboard, the Advanced Topographic Laser Altimeter System (ATLAS), is now providing routine, very high-resolution, surface elevation data across the globe, including the Arctic and Southern Oceans. In this study we demonstrate our new processing chain for converting the along-track ICESat-2 sea ice freeboard product (ATL10) into sea ice thickness, focusing our initial efforts on the Arctic Ocean. For this conversion we primarily make use of snow depth and density data from the NASA Eulerian Snow on Sea Ice Model (NESOSIM). The coarse resolution (~100 km) snow data are redistributed onto the high resolution (~30 - 100 m) ATL10 freeboards using relationships obtained from snow depth and freeboard data collected by NASA’s Operation IceBridge mission. We present regional sea ice thickness distributions and highlight their seasonal evolution through our first winter-season of data collection. We include ice thickness uncertainty estimates, while also acknowledging the limitations of these estimates. We generate a gridded monthly thickness product and compare this with various monthly sea ice thickness estimates obtained from ESA’s CryoSat-2 satellite mission, with ICE 33 Sat-2 showing consistently lower thicknesses. Finally, we compare our February/March 2019 thickness estimates to ICESat February/March (19th February – 21st March) 2008 ice thickness estimates using the same input assumptions, which show a ~0.40 m or ~21% thinning across an inner Arctic Ocean domain in this 11-year time period.

Alek A Petty↗

Multi-Dimensional Spectral Difference Method for Unstructured Grids

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. It combines the best features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. Universal reconstructions are obtained by distributing unknowns in a geometrically similar manner for all unstructured cells. Placements of the unknown and flux points with various order of accuracy are given for the line, triangular and tetrahedral elements. The data structure of the new method permits an optimum use of cache memory, resulting in further computational efficiency on modern computers. A new pointer system is developed that reduces memory requirements and simplifies programming for any order of accuracy. Numerical solutions are presented and compared with the exact solutions for wave propagation problems in both two and three dimensions to demonstrate the capability of the method. Excellent agreement has been found. The method is simpler and more efficient than previous discontinuous Galerkin and spectral volume methods for unstructured grids.

Liu, Yen↗

Three-dimensional inviscid flow in mixers. II - Analysis of turbofan forced mixers

A small disturbance formulation for the three-dimensional potential analysis of the inviscid flow over a turbofan forced mixer configuration in which the governing equations are reduced by means of a flux volume formulation along a Cartesian grid is presently extended to include the effects of power addition within the potential formulation. Calculations are presented for practical turbofan mixer designs, and comparison calculations are also given with measured surface pressure distributions and measured axial velocity profiles.

Barber, T. J.↗

Effects of cumulus entrainment and multiple cloud types on a January global climate model simulation

An improved version of the GISS Model II cumulus parameterization designed for long-term climate integrations is used to study the effects of entrainment and multiple cloud types on the January climate simulation. Instead of prescribing convective mass as a fixed fraction of the cloud base grid-box mass, it is calculated based on the closure assumption that the cumulus convection restores the atmosphere to a neutral moist convective state at cloud base. This change alone significantly improves the distribution of precipitation, convective mass exchanges, and frequencies in the January climate. The vertical structure of the tropical atmosphere exhibits quasi-equilibrium behavior when this closure is used, even though there is no explicit constraint applied above cloud base.

Yao, Mao-Sung↗

Cloud atlas for the FIRE Cirrus Intensive Field Observation (IFO)

An Intensive Field Observation (IFO) of cirrus clouds was conducted over the mid-western U.S. during the period October 13 to November 2, 1986. This activity, part of the First ISCCP Regional Experiment (FIRE), included measurements made from specially deployed instruments on the ground, balloons, and aircraft as well as observations from existing operational and experimental satellites. One of the sets of satellite observations was the radiance measurements made with the 5-channel AVHRR radiometer on the NOAA 9 polar orbiting meteorological satellite. The ground resolution of the measurements at nadir is approx. 1 km. It is these measurements, made once each day at approximately 2:30 p.m. local time, that were used in determining the present cloud atlas. The area covered by the atlas is slightly larger than the area specified for the IFO, in order to be in alignment with the grid that will be used in a forthcoming atlas for the larger, ETO region. The atlas contains four pages of information for each satellite pass. The 1st page of each group shows the distribution of measured radiances in channel 1 (normalized to the incoming solar flux multiplied by the cosine of the solar zenith angle) and in channel 4 for the area as a whole and for each analysis box. The 2nd page shows the images in: channels 1 and 2, channel 3R; and channel 4. The 3rd page shows the retrieved parameters in graphical form for the region as a whole and for each analysis box, where cloud fraction appears as a contour plot with respect to optical thickness and cloudtop temperature. The 4th page provides a statistical summary of the retrieved parameters in numerical form for each analysis box.

Arking, Albert↗