Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid integration study”

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 91 records · Page 5

Use of finite volume schemes for transition simulation

The use of finite-volume methods in the study of spatially and temporally evolving transitional flows over a flat plate is investigated. Schemes are developed with both central and upwind differencing. The compressible Navier-Stokes equations are solved with a Runge-Kutta time stepping scheme. Disturbances are determined using linear theory and superimposed at the inflow boundary. Time accurate integration is then used to allow temporal and spatial disturbance evolution. Characteristic-based boundary conditions are employed. The requirements of using finite-volume algorithms are studied in detail. Special emphasis is placed on difference schemes, grid resolution, and disturbance amplitudes. Moreover, comparisons are made with linear theory for small amplitude disturbances. Both subsonic and supersonic flows are considered, and it is shown that the locations of branch 1 and branch 2 of the neutral stability curve are well predicted, given sufficient resolution.

Fenno, Charles C., Jr.↗

Multigrid properties of upwind-biased data reconstructions

The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.

Warren, Gary P.↗

Validation of a Pressure-Based Combustion Simulation Tool Using a Single Element Injector Test Problem

The traditional design and analysis practice for advanced propulsion systems, particularly chemical rocket engines, relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment by non-CFD specialists. A computational tool, called Loci-STREAM is being developed for this purpose. It is a pressure-based, Reynolds-averaged Navier-Stokes (RANS) solver for generalized unstructured grids, which is designed to handle all-speed flows (incompressible to hypersonic) and is particularly suitable for solving multi-species flow in fixed-frame combustion devices. Loci-STREAM integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective of the ongoing work is to develop a robust simulation capability for combustion problems in rocket engines. As an initial step towards validating this capability, a model problem is investigated in the present study which involves a gaseous oxygen/gaseous hydrogen (GO2/GH2) shear coaxial single element injector, for which experimental data are available. The sensitivity of the computed solutions to grid density, grid distribution, different turbulence models, and different near-wall treatments is investigated. A refined grid, which is clustered in the vicinity of the solid walls as well as the flame, is used to obtain a steady state solution which may be considered as the best solution attainable with the steady-state RANS methodology. From a design point of view, quick turnaround times are desirable; with this in mind, coarser grids are also employed and the resulting solutions are evaluated with respect to the fine grid solution.

Thakur, Siddarth↗

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Newman, James C., III↗

Studies of Shock Wave Interactions with Homogeneous and Isotropic Turbulence

A nearly homogeneous nearly isotropic compressible turbulent flow interacting with a normal shock wave has been studied experimentally in a large shock tube facility. Spatial resolution of the order of 8 Kolmogorov viscous length scales was achieved in the measurements of turbulence. A variety of turbulence generating grids provide a wide range of turbulence scales. Integral length scales were found to substantially decrease through the interaction with the shock wave in all investigated cases with flow Mach numbers ranging from 0.3 to 0.7 and shock Mach numbers from 1.2 to 1.6. The outcome of the interaction depends strongly on the state of compressibility of the incoming turbulence. The length scales in the lateral direction are amplified at small Mach numbers and attenuated at large Mach numbers. Even at large Mach numbers amplification of lateral length scales has been observed in the case of fine grids. In addition to the interaction with the shock the present work has documented substantial compressibility effects in the incoming homogeneous and isotropic turbulent flow. The decay of Mach number fluctuations was found to follow a power law similar to that describing the decay of incompressible isotropic turbulence. It was found that the decay coefficient and the decay exponent decrease with increasing Mach number while the virtual origin increases with increasing Mach number. A mechanism possibly responsible for these effects appears to be the inherently low growth rate of compressible shear layers emanating from the cylindrical rods of the grid.

Briassulis, G.↗

A Central Asia Hydrologic Monitoring Dataset for Food and Water Security Applications in Afghanistan

From the Hindu Kush mountains to the Registan Desert, Afghanistan is a diverse landscape where droughts, floods, conflict, and economic market accessibility pose challenges for agricultural livelihoods and food security. The ability to remotely monitor environmental conditions is critical to support decision making for humanitarian assistance. The Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS) global and Central Asia data streams provide information on hydrologic states for routine integrated food security analysis. While developed for a specific project, these data are publicly available and useful for other applications that require hydrologic estimates of the water and energy balance. These two data streams are unique because of their suitability for routine monitoring, as well as for being a historical record for computing relative indicators of water availability. The global stream is available at ∼ 1-month latency, and monthly average outputs are on a 10 km grid from 1982–present. The second data stream, Central Asia (21–56°N, 30–100°E), at ∼ 1 d latency, provides daily average outputs on a 1 km grid from 2000–present. This paper describes the configuration of the two FLDAS data streams, background on the software modeling framework, selected meteorological inputs and parameters, and results from previous evaluation studies. We also provide additional analysis of precipitation and snow cover over Afghanistan. We conclude with an example of how these data are used in integrated food security analysis. For use in new and innovative studies that will improve understanding of this region, these data are hosted by U.S. Geological Survey data portals and the National Aeronautics and Space Administration (NASA). The Central Asia data described in this paper can be accessed via the NASA repository at https://doi.org/10.5067/VQ4CD3Y9YC0R (Jacob and Slinski, 2021), and the global data described in this paper can be accessed via the NASA repository at https://doi.org/10.5067/5NHC22T9375G (McNally, 2018).

Amy McNally↗

On spurious behavior of CFD simulations

Spurious behavior in underresolved grids and/or semi-implicit temporal discretizations for four computational fluid dynamics (CFD) simulations are studied. The numerical simulations consist of (a) a 1-D chemically relaxed nonequilibrium model, (b) the direct numerical simulation (DNS) of 2-D incompressible flow over a backward facing step, (c) a loosely-coupled approach for a 2-D fluid-structure interaction, and (d) a 3-D compressible unsteady flow simulation of vortex breakdown in delta wings. Using knowledge from dynamical systems theory, various types of spurious behaviors that are numerical artifacts were systematically identified. These studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by the available computing power. In large scale computations underresolved grids, semi-implicit procedures, loosely-coupled implicit procedures, and insufficiently long time integration in DNS are most often unavoidable. Consequently, care must be taken in both computation and in interpretation of the numerical data. The results presented confirm the important role that dynamical systems theory can play in the understanding of the nonlinear behavior of numerical algorithms and in aiding the identification of the sources of numerical uncertainties in CFD.

Numerical Solution↗

Quantifying Food Security and Mitigation Risks Consequential to Climate Change Impacts on Crop Yields

Climate change is expected to impact crop yields globally, with some regions benefiting from favorable conditions and CO2 fertilization, while others face adverse effects from altered precipitation and higher temperatures. Changes in crop yields can destabilize the global food system and pose challenges to food security. Moreover, crop production is crucial, as biofuels are becoming increasingly important contributors to climate change mitigation measures aimed at limiting global warming. This study uses the Integrated Model to Assess the Global Environment integrated assessment model framework to analyze different indicators related to food security and climate change mitigation under varying climate change impacts on crop yields. Twelve spatially explicit crop productivity projections were taken from the full archive of the Global Gridded Crop Model Intercomparison of 120 climate-crop model combinations, forced by CMIP6-based climate scenarios. The selection includes two average-performing climate-crop model combinations, two pessimistic combinations that perform one standard deviation below the mean, and two optimistic model combinations that perform one standard deviation above the mean. To single out the effect of climate change on productivity changes, we drew samples from two representative concentration pathways (RCP2.6 and RCP8.5). These productivity projections were applied within an otherwise uniform scenario (SSP2) and analyzed for their effect on total calorie demand, crop prices, and number of people at risk of undernourishment to quantify food security. Risks to climate change mitigation targets were explored by modeling the total bioenergy supply, emissions, and global mean temperature. The results revealed significant differences in the risk of food security and mitigation potential between different regions and climate change scenarios. Across scenarios, the crop area extent can vary up to 2 million km2 due to changing crop yields. The projected change in global hunger ranges from 60 to 160 million undernourished people, indicating uncertainty between climate and crop model combinations. Low-income regions are especially impacted because of their high sensitivity to changes in food prices. Global climate change mitigation ambitions can also deviate by the latter part of the 21st century, as changes in yields will impact biofuel production as well as agriculture, forestry and other land use emissions. The quantitative insights generated by this study highlight the need for global policy efforts to make the agricultural system more adaptive to climate change to handle potential negative impacts.

crop yields↗

Application of aerial photography to water-related programs in Michigan

Aerial photography and information system technology were used to generate information required for the effective operation of three water-related programs in Michigan. Potential mosquito breeding sites were identified from specially acquired low altitude 70 mm color photography for the city of Lansing; the inventory identified 35% more surface water areas than indicated on existing field maps. A comprehensive inventory of surface water sources and potential access sites was prepared to assist fire departments in Antrim County with fire truck water-recharge operations. Remotely-sensed land cover/use data for Windsor Township, Eaton County, were integrated with other resource data into a computer-based information system for regional water quality studies. Eleven thematic maps focusing on landscape features affecting non-point water pollution and waste disposal were generated from analyses of a four-hectare grid-based data file containing land cover/use, soils, topographic and geologic (well-log) data.

Enslin, W. R.↗

A verification of unsteady Navier-Stokes solutions around oscillating airfoils

A finite difference solution code for the two dimensional Navier-Stokes equations was combined with a moving-grid system. The thin layer Navier-Stokes equations with a turbulence model are solved in a time-accurate manner in order to study the unsteady aerodynamics around airfoils undergoing small amplitude pitching or heaving motions in the transonic regime. The accuracy of the solutions obtained by the use of the present moving-grid technqiue is investigated. The effects of the minimum grid size and the integrating time-step size on the solutions are also checked. Some of the solutions obtained by the present method are compared with experimental results. It is demonstrated that the unsteady aerodynamics around oscillating airfoils can be predicted fairly well by the present code for cases in which the dynamic angle of attack or displacement is small.

Nakamichi, J.↗

Characterizing the Severe Turbulence Environments Associated With Commercial Aviation Accidents: A Real-Time Turbulence Model (RTTM) Designed for the Operational Prediction of Hazardous Aviation Turbulence Environments

Real-time prediction of environments predisposed to producing moderate-severe aviation turbulence is studied. We describe the numerical model and its postprocessing system designed for said prediction of environments predisposed to severe aviation turbulence as well as presenting numerous examples of its utility. The numerical model is MASS version 5.13, which is integrated over three different grid matrices in real time on a university work station in support of NASA Langley Research Center s B-757 turbulence research flight missions. The postprocessing system includes several turbulence-related products, including four turbulence forecasting indices, winds, streamlines, turbulence kinetic energy, and Richardson numbers. Additionally, there are convective products including precipitation, cloud height, cloud mass fluxes, lifted index, and K-index. Furthermore, soundings, sounding parameters, and Froude number plots are also provided. The horizontal cross-section plot products are provided from 16 000 to 46 000 ft in 2000-ft intervals. Products are available every 3 hours at the 60- and 30-km grid interval and every 1.5 hours at the 15-km grid interval. The model is initialized from the NWS ETA analyses and integrated two times a day.

Kaplan, Michael L.↗

Application of aerial photography to water-related programs in Michigan

The paper describes the use of aerial photography and information system technology in the provision of information required for the effective operation of three water-related programs in Michigan. Potential mosquito breeding sites were identified from specially acquired low altitude 70 mm color photography for the City of Lansing Vector Control Area. A comprehensive inventory of surface water sources and potential access sites was prepared to assist fire departments in Antrim County with fire truck water-recharge operations. Remotely-sensed land cover/use data for Windsor Township, Eaton County were integrated with other resource data into a computer-based information system for regional water quality studies. Eleven thematic maps specifically focussed on landscape features affecting non-point water pollution and waste disposal were generated from analyses of a four-hectare grid-based data file containing land cover/use, soils, topographic and geologic (well-log) data.

Enslin, W. R.↗

GEWEX SRB Release 4: Public Release

The NASA/GEWEX Surface Radiation Budget (SRB) project produces 3-hourly shortwave and longwave surface and top of atmosphere radiative fluxes for the 1983-nearpresent time period. Here we announce the public availability of the new Release 4 Integrated data product (R4-IP). R4-IP uses the newly recalibrated and processed ISCCP HXS product as its primary input for cloud and radiance data, replacing ISCCP DX with a ninefold increase in pixel count (10km instead of 30km). This version retains a 1°x1° resolution but benefits from a much larger number of samples per grid box than the SRB Release 3. ISCCP also provides an atmospheric temperature and moisture dataset known as nnHIRS which we use here, along with Seaflux and Landflux surface and near-surface meteorological parameters constituting an “Integrated Product” for energy and water cycle studies. In addition to the input data improvements, several important algorithm improvements have been made since Release 3. These include recalculated shortwave (SW or solar wavelengths) atmospheric transmissivities and reflectivities yielding a somewhat less transmissive atmosphere. Both SW and LW (longwave or thermal infrared) also include variable aerosol composition and radiative properties, allowing for the use of a detailed aerosol history from the Max Planck Institute Aerosol Climatology (MAC). Both LW and SW algorithms now produce pristine sky fluxes, allowing the aerosol flux effects to be quantified. For SW, ocean albedo and snow/ice albedo are improved from Release 3. Total solar irradiance is now variable, and reduced to an average of 1361 Wm-2. The radiative treatment of ice cloud is improved. For the LW, a climatological monthly varying spectral surface emissivity is added.Here we analyze the new SW SRB datasets, comparing them to the previous Release 3, BSRN and PMEL surface measurements. Validation results show improved agreement compared to Release 3. Comparisons to ERBE and CERES TOA and surface satellite datasets before and after overlap periods show the long-term agreement with radiometric uncertainty. Monthly anomalies at global and zonal regions show very good correlation with the shorter ERBE and CERES flux anomalies linking those to the longer term. The uncertainties are discussed and periods of potential issues in inputs and/or algorithms are identified.

Earth energy budget↗

GEWEX SRB Release 4: Public Release

The NASA/GEWEX Surface Radiation Budget (SRB) project produces 3-hourly shortwave and longwave surface and top of atmosphere radiative fluxes for the 1983-near present time period. Here we announce the public availability of the new Release 4 Integrated data product (R4-IP). R4-IP uses the newly recalibrated and processed ISCCP HXS product as its primary input for cloud and radiance data, replacing ISCCP DX with a ninefold increase in pixel count (10km instead of 30km). This version retains a 1°x1° resolution but benefits from a much larger number of samples per grid box than the SRB Release 3. ISCCP also provides an atmospheric temperature and moisture dataset known as nnHIRS which we use here, along with Seaflux and Landflux surface and near-surface meteorological parameters constituting an “Integrated Product” for energy and water cycle studies. In addition to the input data improvements, several important algorithm improvements have been made since Release 3. These include recalculated shortwave (SW or solar wavelengths) atmospheric transmissivities and reflectivities yielding a somewhat less transmissive atmosphere. Both SW and LW (longwave or thermal infrared) also include variable aerosol composition and radiative properties, allowing for the use of a detailed aerosol history from the Max Planck Institute Aerosol Climatology (MAC). Both LW and SW algorithms now produce pristine sky fluxes, allowing the aerosol flux effects to be quantified. For SW, ocean albedo and snow/ice albedo are improved from Release 3. Total solar irradiance is now variable, and reduced to an average of 1361 Wm-2. The radiative treatment of ice cloud is improved. For the LW, a climatological monthly varying spectral surface emissivity is added. Here we analyze the new SW SRB datasets, comparing them to the previous Release 3, BSRN and PMEL surface measurements. Validation results show improved agreement compared to Release 3. Comparisons to ERBE and CERES TOA and surface satellite datasets before and after overlap periods show the long-term agreement with radiometric uncertainty. Monthly anomalies at global and zonal regions show very good correlation with the shorter ERBE and CERES flux anomalies linking those to the longer term. The uncertainties are discussed and periods of potential issues in inputs and/or algorithms are identified.

Earth Radiation Budget↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the 'multiple grids' or 'zonal grids' approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-stage Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, M.↗

A conservative approach for flow field calculations on multiple grids

In the computation of flow fields about complex configurations, it is very difficult to construct body-fitted coordinate systems. An alternative approach is to use several grids at once, each of which is generated independently. This procedure is called the multiple grids or zonal grids approach and its applications are investigated in this study. The method follows the conservative approach and provides conservation of fluxes at grid interfaces. The Euler equations are solved numerically on such grids for various configurations. The numerical scheme used is the finite-volume technique with a three-state Runge-Kutta time integration. The code is vectorized and programmed to run on the CDC VPS-32 computer. Some steady state solutions of the Euler equations are presented and discussed.

Kathong, Monchai↗

Numerical investigation of turbulent channel flow

Fully developed turbulent channel flow was simulated numerically at Reynolds number 13800, based on centerline velocity and channel halt width. The large-scale flow field was obtained by directly integrating the filtered, three dimensional, time dependent, Navier-Stokes equations. The small-scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC IV computer with up to 516,096 grid points. The computed flow field was used to study the statistical properties of the flow as well as its time dependent features. The agreement of the computed mean velocity profile, turbulence statistics, and detailed flow structures with experimental data is good. The resolvable portion of the statistical correlations appearing in the Reynolds stress equations are calculated. Particular attention is given to the examination of the flow structure in the vicinity of the wall.

Moin, P.↗

The GEOS-iODAS: Description and Evaluation

This report documents the GMAO's Goddard Earth Observing System sea ice and ocean data assimilation systems (GEOS iODAS) and their evolution from the first reanalysis test, through the implementation that was used to initialize the GMAO decadal forecasts, and to the current system that is used to initialize the GMAO seasonal forecasts. The iODAS assimilates a wide range of observations into the ocean and sea ice components: in-situ temperature and salinity profiles, sea level anomalies from satellite altimetry, analyzed SST, and sea-ice concentration. The climatological sea surface salinity is used to constrain the surface salinity prior to the Argo years. Climatological temperature and salinity gridded data sets from the 2009 version of the World Ocean Atlas (WOA09) are used to help constrain the analysis in data sparse areas. The latest analysis, GEOS ODAS5.2, is diagnosed through detailed studies of the statistics of the innovations and analysis departures, comparisons with independent data, and integrated values such as volume transport. Finally, the climatologies of temperature and salinity fields from the Argo era, 2002-2011, are presented and compared with the WOA09.

Climate↗