Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid code”

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 1,621 records · Page 90

ninterp: N-dimensional interpolation Rust crate [SWR-25-25]

The ninterp crate provides multivariate interpolation over a regular, sorted, nonrepeating grid of any dimensionality. A variety of interpolation strategies are implemented, however more are likely to be added. Extrapolation beyond the range of the supplied coordinates is supported for 1-D linear interpolators, using the slope of the nearby points. There are hard-coded interpolators for lower dimensionalities (up to N = 3) for better runtime performance. All interpolation is handled through instances of the Interpolator enum, with the selected tuple variant containing relevant data. Interpolation is executed by calling Interpolator::interpolate.

Carow, Kyle↗

PATCHWORK: A Multipatch Infrastructure for Multiphysics/Multiscale/Multiframe Fluid Simulations

Here, we present a "multipatch" infrastructure for the numerical simulation of fluid problems in which subregions require different grid scales, different grid geometries, different physical equations, or different reference frames. Its key element is a sophisticated client–router–server framework for efficiently linking processors supporting different regions ("patches") that must exchange boundary data. This infrastructure may be used with a wide variety of fluid dynamics codes; the only requirement is that their primary dependent variables be the same in all patches, e.g., fluid mass density, internal energy density, and velocity. Its structure can accommodate either Newtonian or relativistic dynamics. The overhead imposed by this system is both problem and computer cluster architecture dependent. Compared to a conventional simulation using the same number of cells and processors employed on a problem not requiring multipatch methods, the cell update per processor rate decreases by an amount that can range from negligible to a factor of a few; however, even in these problems, the infrastructure can permit substantial decreases in the total number of cell updates required.

97 MATHEMATICS AND COMPUTING↗

Numerical simulation of a conducting disk partially covered by an insulator interacting with a plasma

A numerical simulation has been made of the 'pinhole effect' which produces the characteristic 'S-shaped' current-voltage curve. A disk-shaped conducting probe immersed in a plasma is modeled using a particle-in-cell (PIC) code. A probe partially covered by a very thin insulating layer is considered, as well as a probe mounted on an insulating disk. The simulation uses a cylindrical particle mover and allows for a variable number of particles in the system. The simulation space grid uses three different mesh sizes, the coarsest being away from the probe and the finest near the probe, in order to accurately calculate the trajectories of the simulation particles contributing the current to the probe and to the surface charge density on the dielectric. The calculation of the electrostatic potential is done self-consistently using successive over-relaxation (SOR). Backscattering and secondary electron emission are included for the case of positive probe voltage.

Chaky, R. C.↗

Nonreflecting boundary conditions for the complete unsteady transonic small-disturbance equation

Nonreflecting far-field boundary conditions that are consistent with the complete transonic small-disturbance (TSD) equations are derived. They are implemented in a new code for solving the complete TSD equation and are tested for a harmonically oscillating NACA 64A010 airfoil in transonic flow and for a flat-plate airfoil with a pulse in the angle of attack. Using the new boundary conditions on a relatively small grid, solutions for the airfoil that are obtained that agree with large-grid calculations, resulting in a 44 percent savings in computer time. Frequency responses for the flat plate show that most of the disturbances incident on the computational boundaries are absorbed by the boundary conditions.

Whitlow, W., Jr.↗

CFD Predictions of Sonic-Boom Characteristics for Unmodified and Modified SR-71 Configurations

Shaped sonic-boom signatures refer to signatures that look something other than the typical N-waves. Shaped sonic-boom signatures such as "flat-top," "ramp-type," or "hybrid-type" waveforms have been shown to reduce the subjective loudness without requiring reductions in overpressure peaks. The shaping of sonic-boom signatures requires increasing the shock rise time and changes in frequency spectra. So far, a flat-top waveform was shown to be achievable in wind tunnels; however, the influence of long propagation distance and real atmosphere on shaped signatures should be addressed using flight tests. Two different approaches have been proposed for sonic-boom minimization flight tests. The first approach, proposed by Eagle Aerospace, is for a flight test using a modified BQM-34 "FIREBEE" remotely piloted vehicle. The 30-foot long FIREBEE has a steady state flight condition at the Mach number and altitude of interest, and it can be recovered by helicopter from the water. As an alternative approach, a modified SR-71 vehicle has been proposed by the McDonnell Douglas Corporation. Benefits of the SR-71 include its variable geometry supersonic inlets, small cockpit bulge, higher Mach number capabilities, slender design, and longer length (105 foot). The present investigation addresses the sonic-boom analysis for the second vehicle.The objective of the current investigation is to assess the feasibility of a modified SR-71 configuration, with McDonnell Douglas-designed fuselage modifications, intended to produce shaped sonic-boom signatures on the ground. The present study describes the use of a higher-order computational fluid dynamics (CFD) method to predict the sonic-boom characteristics for both unmodified and modified SR-71 configurations. An Euler unstructured grid methodology is used to predict the near-field, three-dimensional pressure patterns generated by both SR-71 models. The computed near-field pressure signatures are extrapolated to specified distances below the aircraft down to impingement on the ground using the code MDBOOM. Comparisons of the near-field pressure signatures with available flight-test data are presented in the current paper.

Fouladi, Kamran↗

The NASA/Ames Mars General Circulation Model: Model Improvements and Comparison with Observations

For many years, the NASA/Ames Mars General Circulation Model (GCM) has been built around the UCLA B-grid dynamical core. An attached tracer transport scheme based on the aerosol microphysical model of Toon et al. (1988) provided a tool for studying dust storm transport and feedbacks (Murphy et al., 1995). While we still use a B-grid version of the model, the Ames group is now transitioning to the ARIES/GEOS Goddard C-grid dynamical core (Suarez and Takacs, 1995). The C-grid produces smoother fields when the model top is raised above 50 km, and has a built in transport scheme for an arbitrary number of tracers. All of our transport simulations are now carried out with the C-grid. We have also been updating our physics package. Several years ago we replaced our bulk boundary layer scheme with a level 2 type diffusive scheme, and added a multi-level soil model (Haberle et al., 2000). More recently we replaced our radiation code with a more generalized two-stream code that accounts for aerosol multiple scattering and gaseous absorption. This code gives us much more flexibility in choosing aerosol optical properties and radiatively active gases.

Haberle, R. M.↗

Fluid Structure Interaction in a Turbine Blade

An unsteady, three dimensional Navier-Stokes solution in rotating frame formulation for turbomachinery applications is presented. Casting the governing equations in a rotating frame enabled the freezing of grid motion and resulted in substantial savings in computer time. The turbine blade was computationally simulated and probabilistically evaluated in view of several uncertainties in the aerodynamic, structural, material and thermal variables that govern the turbine blade. The interconnection between the computational fluid dynamics code and finite element structural analysis code was necessary to couple the thermal profiles with the structural design. The stresses and their variations were evaluated at critical points on the Turbine blade. Cumulative distribution functions and sensitivity factors were computed for stress responses due to aerodynamic, geometric, mechanical and thermal random variables.

Gorla, Rama S. R.↗

AEROELASTIC SIMULATION TOOL FOR INFLATABLE BALLUTE AEROCAPTURE

A multidisciplinary analysis tool is under development for predicting the impact of aeroelastic effects on the functionality of inflatable ballute aeroassist vehicles in both the continuum and rarefied flow regimes. High-fidelity modules for continuum and rarefied aerodynamics, structural dynamics, heat transfer, and computational grid deformation are coupled in an integrated multi-physics, multi-disciplinary computing environment. This flexible and extensible approach allows the integration of state-of-the-art, stand-alone NASA and industry leading continuum and rarefied flow solvers and structural analysis codes into a computing environment in which the modules can run concurrently with synchronized data transfer. Coupled fluid-structure continuum flow demonstrations were conducted on a clamped ballute configuration. The feasibility of implementing a DSMC flow solver in the simulation framework was demonstrated, and loosely coupled rarefied flow aeroelastic demonstrations were performed. A NASA and industry technology survey identified CFD, DSMC and structural analysis codes capable of modeling non-linear shape and material response of thin-film inflated aeroshells. The simulation technology will find direct and immediate applications with NASA and industry in ongoing aerocapture technology development programs.

Liever, P. A.↗

Analysis of a Cylindrical Specimen Heated by an Impinging Hot Hydrogen Jet

A computational conjugate heat transfer methodology was developed, as a first step towards an efficient and accurate multiphysics, thermo-fluid computational methodology to predict environments for hypothetical solid-core, nuclear thermal engine thrust chamber and components. A solid conduction heat transfer procedure was implemented onto a pressure-based, multidimensional, finite-volume, turbulent, chemically reacting, thermally radiating, and unstructured grid computational fluid dynamics formulation. The conjugate heat transfer of a cylindrical material specimen heated by an impinging hot hydrogen jet inside an enclosed test fixture was simulated and analyzed. The solid conduction heat transfer procedure was anchored with a standard solid heat transfer code. Transient analyses were then performed with ,variable thermal conductivities representing three composites of a material utilized as flow element in a legacy engine test. It was found that material thermal conductivity strongly influences the transient heat conduction characteristics. In addition, it was observed that high thermal gradient occur inside the cylindrical specimen during an impulsive or a 10 s ramp start sequence, but not during steady-state operations.

Wang, Ten-See↗

Assessing dual continuum method for multicomponent reactive transport

The subsurface flow and reactive transport code PFLOTRAN can accommodate fracture–matrix interaction for multicomponent reactive transport using a Dual Continuum formulation. In this work we present two examples to assess this approach. Through these examples, we also illustrate how this Dual Continuum model can be parameterised and what kind of grid spacing is needed to accurately capture reactive transport in the secondary continuum. The results of the calculations show that when geochemical reactions in the matrix (or the secondary continuum) are included, a very fine grid spacing close to the fracture–matrix interface is required to properly capture geochemical gradients at the interface. We also show that one can reduce the computational cost by using a variable grid spacing, that limits the total number in the secondary continuum, and yet obtain accurate results. Furthermore, this feature makes the Dual Continuum approach suitable for large-scale reactive transport modelling in fractured crystalline rock.

58 GEOSCIENCES↗

High Performance Approximate Computing

This code repository contains the implementation of the "High-Performance Approximate Computing" (HPAC) toolkit. The toolkit allows you to approximate your own C/C++. The developer uses "pragma's" to annotate code regions as approximate. The compiler extensions lower these pragmas to either compiletime approximate techniques or runtime approximation techniques. At execution time, the implemented runtime system decides which annotated regions it should approximate. HPAC also provides a set of script utilities. The utilities perform a grid search within approximation parameters and performance. The user can analyze the raw data to identify optimal approximation techniques for the application.

Parasyris, Konstantinos↗

Impact of LWR Assembly Structural Features on Cladding Burst Behavior under LOCA Conditions

This report builds on the previous work by assessing the impact of assembly structural features on cladding balloon and burst behavior. The BISON fuel performance code was used to generate two dimensional RZ meshes containing structural features commonly used in nuclear fuel assemblies. The results were then compared to the cladding burst temperature and balloon strain results from the previous work. A mesh sensitivity study was performed to ensure spacer grid effects were appropriately considered, resulting in a more refined mesh than the previous study. The balloon deformation and burst times of the cladding were compared. Consideration was also given to the effect of rod initial pressure. Overall, the BISON simulations indicated the spacer grids can locally influence cladding performance, however, the impact on the cladding performance was minor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Monthly Quality-filtered Aggregation of NOAA Climate Data Record (CDR) of AVHRR Leaf Area Index (LAI) and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), Version 5

This dataset contains gridded monthly Leaf Area Index (LAI) derived from the daily NOAA Climate Data Record (CDR) of AVHRR Leaf Area Index (LAI) and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), Version 5. This data record spans from 1981 to 2018 using data from eight NOAA polar orbiting satellites: NOAA-7, -9, -11, -14, -16, -17, -18 and -19. The data are projected on a 0.05 degree x 0.05 degree global grid, as in the original CDR. The original CDR is one of the Land Surface CDR Version 5 products produced by the NASA Goddard Space Flight Center (GSFC) and the University of Maryland (UMD), which is accompanied by algorithm documentation, data flow diagram and source code for the NOAA CDR Program. This dataset is in the netCDF-4 file format following ACDD and CF Conventions. This dataset has applied quality assurance information to only include "OK" data from the original CDR in the monthly aggregation.

Vermote, Eric [NASA Goddard Space Flight Center (G↗

Documentation of the Fourth Order Band Model

A general circulation model is presented which uses quadratically conservative, fourth order horizontal space differences on an unstaggered grid and second order vertical space differences with a forward-backward or a smooth leap frog time scheme to solve the primitive equations of motion. The dynamic equations for motion, finite difference equations, a discussion of the structure and flow chart of the program code, a program listing, and three relevent papers are given.

Kalnay-Rivas, E.↗

A perspective of theoretical and applied computational fluid dynamics

The present effectiveness and future promise of computational fluid dynamics (CFD) are discussed. The CFD ideal of computer-designed flight vehicles and the progress that has been made toward its achievement, expressed in terms of primary and secondary pacing items, are also addressed. The primary pacing items listed by Chapman - grid generation, turbulence modeling, solution methodology development, and mainframe computer methods and architecture - are modified and reviewed. Important secondary pacing items in CFD development are also examined: algorithm development, complex geometry definition, and predata and postata processing. In addition, user requirements placed on CFD computer codes are defined and some problems associated with the dissemination of CFD technology are identified.

Kutler, P.↗

The CMC/3DPNS computer program for prediction of three-dimension, subsonic, turbulent aerodynamic juncture region flow. Volume 2: Users' manual

The CMC fluid mechanics program system was developed to transmit the theoretical solution of finite element numerical solution methodology, applied to nonlinear field problems into a versatile computer code for comprehensive flow field analysis. Data procedures for the CMC 3 dimensional Parabolic Navier-Stokes (PNS) algorithm are presented. General data procedures a juncture corner flow standard test case data deck is described. A listing of the data deck and an explanation of grid generation methodology are presented. Tabulations of all commands and variables available to the user are described. These are in alphabetical order with cross reference numbers which refer to storage addresses.

Manhardt, P. D.↗

Aerodynamic analysis for aircraft with nacelles, pylons, and winglets at transonic speeds

A computational method has been developed to provide an analysis for complex realistic aircraft configurations at transonic speeds. Wing-fuselage configurations with various combinations of pods, pylons, nacelles, and winglets can be analyzed along with simpler shapes such as airfoils, isolated wings, and isolated bodies. The flexibility required for the treatment of such diverse geometries is obtained by using a multiple nested grid approach in the finite-difference relaxation scheme. Aircraft components (and their grid systems) can be added or removed as required. As a result, the computational method can be used in the same manner as a wind tunnel to study high-speed aerodynamic interference effects. The multiple grid approach also provides high boundary point density/cost ratio. High resolution pressure distributions can be obtained. Computed results are correlated with wind tunnel and flight data using four different transport configurations. Experimental/computational component interference effects are included for cases where data are available. The computer code used for these comparisons is described in the appendices.

Boppe, Charles W.↗

Wall functions for the kappa-epsilon turbulence model in generalized nonorthogonal curvilinear coordinates

A k-epsilon turbulence model suitable for compressible flow, including the new wall function formulation, has been incorporated into an existing compressible Reynolds-averaged Navier-Stokes code, F3D. The low Reynolds number k-epsilon model of Chien (1982) was added for comparison with the present method. A number of features were added to the F3D code including improved far-field boundary conditions and viscous terms in the streamwise direction. A series of computations of increasing complexity was run to test the effectiveness of the new formulation. Flow over a flat plate was computed by using both orthogonal and nonorthogonal grids, and the friction coefficients and velocity profiles compared with a semi-empirical equation. Flow over a body of revolution at zero angle of attack was then computed to test the method's ability to handle flow over a curved surface. Friction coefficients and velocity profiles were compared to test data. All models gave good results on a relatively fine grid, but only the wall function formulation was effective with coarser grids. Finally, in order to demonstrate the method's ability to handle complex flow fields, separated flow over a prolate spheroid at angle of attack was computed, and results were compared to test data. The results were also compared to a k-epsilon model by Kim and Patel (1991), in which one equation model patched in at the wall was employed. Both models gave reasonable solutions, but improvement is required for accurate prediction of friction coefficients in the separated regions.

Sondak, D. L.↗