Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid codes”

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 865 records · Page 48

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interpretable, Extensible Linear Regression Model for Electron Density Prediction

This code allows us to generate a model to predict ground state electron densities directly from atomic structure thereby bypassing the need to perform direct DFT calculations. It is based on many body correlation descriptors defined in https://arxiv.org/abs/2004.14442 (IM release numbers: LLNL-JRNL- 808791 and LLNL-JRNL-840698). The code reads in electron densities, in CHGCAR format used by the DFT code VASP, calculates many-body correlation descriptors based on atoms surrounding a grid point and then uses linear regression to obtain the fitting parameters.

Kumar, Shashikant↗

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High sensitivity of simulated fog properties to parameterized aerosol activation in case studies from ParisFog

Aerosols influence fog properties such as visibility and lifetime by affecting fog droplet number concentrations (N d ). Numerical weather prediction (NWP) models often represent aerosol–fog interactions using highly simplified approaches. Incorporating prognostic size-resolved aerosol microphysics from climate models could allow them to simulate N d and aerosol–fog interactions without incurring excessive computational expense. However, microphysics code designed for coarse spatial resolution may struggle with sub-kilometer-scale grid spacings. Here, we test the ability of the UK Met Office Unified Model to simulate aerosol and fog properties during case studies from the ParisFog field campaign in 2011. We examine the sensitivity of fog properties to variations in N d caused by modifications to simulated aerosol activation. Our model, with a 500 m horizontal resolution and interactive aerosol and cloud microphysics, significantly underpredicts N d , although it only slightly underestimates the cloud condensation nuclei concentration. With an updated version of the Abdul-Razzak and Ghan (2000) activation scheme, we produce N d that are more consistent with those predicted by a cloud parcel model under fog-like conditions. We activate droplets only by adiabatic cooling. We incorporate more realistic hygroscopicities for sulfate and organic aerosols and explore the sensitivity of simulated N d to unresolved updrafts. We find that both N d and simulated fog liquid water content are very sensitive to the updated activation scheme but remain less affected by the update to hygroscopicities. Our improvements offer insights into the physical processes regulating N d in stable conditions, potentially laying foundations for improved operational fog forecasts that incorporate interactive aerosol simulations or aerosol climatologies.

Ghosh, Pratapaditya [Carnegie Mellon University, P↗

A chimera grid scheme

A mesh system composed of multiple overset body-conforming grids is described for adapting finite-difference procedures to complex aircraft configurations. In this so-called 'chimera mesh,' a major grid is generated about a main component of the configuration and overset minor grids are used to resolve all other features. Methods for connecting overset multiple grids and modifications of flow-simulation algorithms are discussed. Computational tests in two dimensions indicate that the use of multiple overset grids can simplify the task of grid generation without an adverse effect on flow-field algorithms and computer code complexity.

Steger, J. L.↗

Complete flowfields over low and wide angle AOTV conceptual configurations

Flowfields over two conceptual aeroassisted orbital transfer vehicles, a moderate lift-to-drag ratio biconic, and an axisymmetric zero lift aerobrake are investigated. The solution procedures employ a compressible Navier-Stokes time-asymptotic technique for blunt body flow and a Parabolized Navier-Stokes spatial marching technique for the primarily supersonic shock layer flow over low-angle bodies. Emphasis is placed on the laminar convective heating predictions and comparisons to experimental data where available. Code robustness and solution times are discussed and guidelines are given for the specification of various numerical parameters. The Parabolized Navier-Stokes technique has a high probability of solution success for the range of applications considered herein, and very good comparisons with experimental heat-transfer data have been obtained. The blunt body Navier-Stokes code is more sensitive to parameter specification for the aerobrake calculations. Still, coarse grid solutions can be obtained quickly to get a good first cut at flowfield definition. Fair to good comparisons have been obtained between experimental heating and pressure data and coarse and fine grid solutions.

Gnoffo, P. A.↗

Discrete element weld model, phase 2

A numerical method was developed for analyzing the tungsten inert gas (TIG) welding process. The phenomena being modeled include melting under the arc and the flow in the melt under the action of buoyancy, surface tension, and electromagnetic forces. The latter entails the calculation of the electric potential and the computation of electric current and magnetic field therefrom. Melting may occur at a single temperature or over a temperature range, and the electrical and thermal conductivities can be a function of temperature. Results of sample calculations are presented and discussed at length. A major research contribution has been the development of numerical methodology for the calculation of phase change problems in a fixed grid framework. The model has been implemented on CHAM's general purpose computer code PHOENICS. The inputs to the computer model include: geometric parameters, material properties, and weld process parameters.

Prakash, C.↗

A parametric study of three-dimensional separation at a wing/body junction for supersonic free-stream conditions

The problem of three-dimensional separation at a wing/body junction has been investigated numerically using a three-dimensional Navier-Stokes code which employs the MacCormack's time split finite volume technique. An algebraic grid generation technique is used for generating the grid at a wing/body junction. Specific computational results on velocity and pressure distribution in the separated flow region are compared with the experimental results. A parametric study of flow parameters such as Mach number and Reynolds number have been carried out to understand their effect in interaction flow field. The parametric study indicates a strong dependency of the number of vortices at the junction on Mach number and Reynolds number.

Lakshmanan, B.↗

An implicit three-dimensional Navier-Stokes solver for compressible flow

A three-dimensional numerical method based on the lower-upper symmetric-Gauss-Seidel implicit scheme in conjunction with the flux-limited dissipation model is developed for solving the compressible Navier-Stokes equations. A new computer code which is based on this method requires only 9 microsec per grid-point per iteration on a single processor of a Cray YMP computer and executes at the sustained rate of 170 MFLOPS. A reduction of 4 orders of magnitude in the residual for a high Reynolds number flow using 230 K grid points is obtained in 24 minutes. The computational results compare well with available experimental data.

Yoon, Seokkwan↗

Grid generation and compressible flow computations about a high-speed civil transport configuration

Techniques and software are discussed for generating grids about a high-speed civil transport configuration. The configuration is defined by a computer-aided design system in wing, fuselage, tail and engine-nacelle components. Grid topology and the surfaces outlining the blocks of the topology are computed with interactive software. The volume grid is computed using software based on transfinite interpolation and Lagrangian blending functions. Several volume grids for inviscid and viscous flow have been generated using this system of codes. Demonstration flowfields around this vehicle are described.

Abolhassani, J. S.↗

Computer Simulation Performed for Columbia Project Cooling System

This demo shows a high-fidelity simulation of the air flow in the main computer room housing the Columbia (10,024 intel titanium processors) system. The simulation asseses the performance of the cooling system and identified deficiencies, and recommended modifications to eliminate them. It used two in house software packages on NAS supercomputers: Chimera Grid tools to generate a geometric model of the computer room, OVERFLOW-2 code for fluid and thermal simulation. This state-of-the-art technology can be easily extended to provide a general capability for air flow analyses on any modern computer room. Columbia_CFD_black.tiff

Ahmad, Jasim↗

Validation of Multi-Dimensional Stirling Engine Design Codes: Measurements in the 90-Degree Turn Test Section

The work to be presented herein was motivated largely by a desire to improve the understanding of oscillatory fluid mechanics inside a Stirling engine. To this end, a CFD project was undertaken at Cleveland State University with the goal of accurately predicting the fluid dynamics within an engine or engine component. Along with the CFD efforts, a code validation project was undertaken at the University of Minnesota. The material covered herein consists of four main parts. In section 1, an experimental investigation of a small aspect ratio impinging jet is discussed. Included in this discussion is a description of the test facilities and instrumentation. A presentation of the collected data is given and comments are made. Next, in section 2, a parallel experimental investigation is presented in which the same geometry as that of section 1 is used, but the flow conditions are changed from steady unidirectional flow to sinusoidally oscillating flow. In section Two, collected data are presented and comments are made. In section 3, a comparison is made between the results of sections 1 and 2, namely, sinusoidally oscillating flow results are compared to steady, unidirectional flow results from the same geometry. Finally, in section 4, a comparison is made between experimentally collected data (the main subject of this work) and CFD generated results. Furthermore, in appendix A, an introductory description of the primary measurement tool used in the experimental process the hot wire anemometer is given for the unfamiliar. The anemometer calibration procedure is described in appendix B. A portfolio of data reduction and data processing codes is provided in appendix C and lastly, a DVD and a roadmap of its contents is provided in an appendix D. 1.0 Unidirectional Flow Investigations 1.1 Introduction This unidirectional experimental program was undertaken to complement an oscillatory flow investigation conducted at the University of Minnesota. The oscillatory investigation is discussed thoroughly in section 2. We defer the description of the motivation behind these experiments until the introduction of section 2. The work that is discussed in this thesis began (chronologically) with oscillatory flow visualization experiments. It was decided that it would be valuable and important to investigate the flow under unidirectional conditions in the same geometry as that of the oscillatory experiments. The thought was that the unidirectional case would be less complicated to model with a CFD program (a moving boundary would be replaced with a steady state boundary condition). Thus, a series of unidirectional experiments were carried out to capture the important features of the flow within the test section. The purpose of these experiments was to provide a data set for comparison to CFD generated velocity fields. Hot-wire anemometry data were taken and flow visualization was conducted as a standard for code validation. The flow geometry was simple, such that it could be easily gridded in a CFD program. However, the geometry provided separation and transition zones, shear layers and recirculation zones. These characteristics made the flow complex and challenging for CFD computation. We comment that the order of experiments that produced this report is as follows: experimental flow visualization under oscillatory flow conditions was carried out; this was followed by unidirectional flow visualization and hot wire anemometry; finally, oscillatory hot wire anemometry was conducted. We present the results out of chronological order for the following reason: the unidirectional results are easier

Simon, Terrence W.↗

Rapid Prediction of Installed Jet Noise from RANS

A new method of computing jet noise, called mSrc, was developed on the general principles of acoustic analogies. In the method, the problem of translating turbulent flow energy into acoustic energy at a far-field observer is broken into two parts, a calculation of acoustic source strengths and then their propagation. The acoustic sources are related to turbulent quantities in the jet plume in a robust manner. The propagation, which is more properly computed using a Green's function accounting for nonuniform speed of sound and solid surfaces, is instead modeled using commonly observed features of jet acoustic far field directivity, and by diffraction barrier theory for surfaces. The mSrc method does not require gradients of the predicted flow field, allowing it to make use of robust unstructured RANS CFD methods, including embedded boundary codes. Such codes do not require specification of surface meshes, and auto-refine their grid to resolve flow gradients, putting resolution where it is required without a priori user input. The ability of mSrc to use such radically unstructured flow input results in an efficient method of estimating noise from jet flows from complex nozzles installed on aircraft. Many validation cases are presented to demonstrate the accuracy and range of applicability of the mSrc method for representative jet noise applications.

Bridges, James↗

Rapid Prediction of Installed Jet Noise From RANS

A new method of computing jet noise, called mSrc, was developed on the general principles of acoustic analogies. In the method, the problem of translating turbulent flow energy into acoustic energy at a far-field observer is broken into two parts, a calculation of acoustic source strengths and then their propagation. The acoustic sources are related to turbulent quantities in the jet plume in a robust manner. The propagation, which is more properly computed using a Green's function accounting for nonuniform speed of sound and solid surfaces, is instead modeled using commonly observed features of jet acoustic far field directivity, and by diffraction barrier theory for surfaces. The mSrc method does not require gradients of the predicted flow field, allowing it to make use of robust unstructured RANS CFD methods, including embedded boundary codes. Such codes do not require specification of surface meshes, and auto-refine their grid to resolve flow gradients, putting resolution where it is required without a priori user input. The ability of mSrc to use such radically unstructured flow input results in an efficient method of estimating noise from jet flows from complex nozzles installed on aircraft. Many validation cases are presented to demonstrate the accuracy and range of applicability of the mSrc method for representative jet noise applications.

Bridges, James↗

Formation of the Lyman Continuum during Solar Flares

The Lyman continuum (LyC; <911.12 Å) forms at the top of the chromosphere in the quiet Sun, making LyC a powerful tool for probing the chromospheric plasma during solar flares. To understand the effects of nonthermal energy deposition in the chromosphere during flares, we analyzed LyC profiles from a grid of field-aligned radiative-hydrodynamic models generated using the RADYN code as part of the F-CHROMA project. The spectral response of LyC, the temporal evolution of the departure coefficient of hydrogen, b 1 , and the color temperature, T c , in response to a range of nonthermal electron distribution functions, were investigated. The LyC intensity was seen to increase by 4–5.5 orders of magnitude during solar flares, responding most strongly to the nonthermal electron flux of the beam. Generally, b 1 decreased from 10 2 –10 3 to closer to unity during solar flares, indicating a stronger coupling to local conditions, while T c increased from 8–9 to 10–16 kK. T c was found to be approximately equal to the electron temperature of the plasma when b1 was at a minimum. Both optically thick and optically thin components of LyC were found to be in agreement with the interpretation of recent observations. The optically thick layer forms deeper in the chromosphere during a flare compared to quiescent periods, whereas the optically thin layers form at higher altitudes due to chromospheric evaporation, in low-temperature, high-density regions propagating upward. We put these results in the context of current and future missions.

Solar flares↗

Capsule_Grid User Guide

Software utility CAPSULE_GRID and several related procedures that automate grid generation for computational fluid dynamics calculations about axisymmetric atmospheric entry vehicles are described as a supplement to the documentation within the source codes. This application program and its ancillary software utilities are all in the public domain at the web site https://software.nasa.gov/software/category/all/arc/1/cfdtools with the exception of the 2D hyperbolic volume gridding step for which scripts driving either GRIDGEN1 or HYPGEN2 are available. CAPSULE_GRID was developed at NASA Ames Research Center in support of the DPLR3 real gas flow solver. Both forebodies (only) and full bodies are treated, for axisymmetric 2D flow calculations or for 3D calculations at angle of attack. The expected application is to the smooth outer mold line of a hypersonic atmospheric entry vehicle, for which a single layer of point-matched structured grid blocks normally suffices. Control files for a variety of representative cases are included with illustrations of results, and some insights into the subtleties of curvature-based grid point distribution are provided. Gridding of a capsule on a sting is supported. Additional procedures for treating arbitrarily long 2D wakes, prompted by asteroid studies and relevant to spacecraft such as Stardust, are also documented. Most recently, treating an asymmetric forebody (including an off-center nose, with or without an aft body) has been incorporated as an extra step.

spacecraft↗

Validation of a Climate-Data Record of the "Clear-Kky" Surface Temperature of the Greenland Ice Sheet

Surface temperatures on the Greenland Ice Sheet have been studied on the ground, using automatic weather station (AWS) data from the Greenland-Climate Network (GC-Net), and from analysis of satellite sensor data. Using Advanced Very High Frequency Radiometer (AVHRR) weekly surface temperature maps, warming of the surface of the Greenland Ice Sheet has been documented since 1981. We extended and refined this record using higher-resolution Moderate-Resolution Imaging Spectroradiometer (MODIS) data from March 2000 to the present. We developed a daily and monthly climate-data record (CDR) of the "clear-sky" surface temperature of the Greenland Ice Sheet using an ice-surface temperature (1ST) algorithm developed for use with MODIS data. Validation of this CDR is ongoing. MODIS Terra swath data are projected onto a polar stereographic grid at 6.25-km resolution to develop binary, gridded daily and mean-monthly 1ST maps. Each monthly map also has a color-coded image map that is available to download. Also included with the monthly maps is an accompanying map showing number of days in the month that were used to calculate the mean-monthly 1ST. This is important because no 1ST decision is made by the algorithm for cells that are considered cloudy by the internal cloud mask, so a sufficient number of days must be available to produce a mean 1ST for each grid cell. Validation of the CDR consists of several facets: 1) comparisons between ISTs and in-situ measurements; 2) comparisons between ISTs and AWS data; and 3) comparisons of ISTs with surface temperatures derived from other satellite instruments such as the Thermal Emission and Reflection Radiometer (ASTER) and Enhanced Thematic Mapper Plus (ETM+). Previous work shows that Terra MODIS ISTs are about 3 C lower than in-situ temperatures measured at Summit Camp, during the winter of 2008-09 under clear skies. In this work we begin to compare surface temperatures derived from AWS data with ISTs from the MODIS CDR. The Greenland Ice Sheet 1ST CDR will be useful for monitoring surface-temperature trends and can be used as input or for validation of climate models. The CDR can be extended into the future using MODIS Terra, Aqua and NPOESS Preparatory Project Visible Infrared Imager Radiometer Suite (VII RS) data.

Hall, Dorothy K.↗