Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Stress”

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 73 records · Page 4

Numerical methods for large-scale, time-dependent partial differential equations

A survey of numerical methods for time dependent partial differential equations is presented. The emphasis is on practical applications to large scale problems. A discussion of new developments in high order methods and moving grids is given. The importance of boundary conditions is stressed for both internal and external flows. A description of implicit methods is presented including generalizations to multidimensions. Shocks, aerodynamics, meteorology, plasma physics and combustion applications are also briefly described.

Turkel, E.↗

The DE-AC accelerator system for high beam current density

A radically different ion accelerator system design has been developed to increase the extracted beam current density of ion thrusters. By independently controlling plasma ion extraction and acceleration, the Decoupled Extraction and Amplified Current (DE-AC) accelerator system gives higher beam currents at lower electric field stress levels. Results are presented for DE-AC grid system operation with the Test Bed Ion Engine. Using argon propellant, feasibility of a 60 mA/sq cm beam current density is shown. A discussion of power supply requirements, specific impulse operation and grid system geometrical design is presented also.

Aston, G.↗

Calculation of 3D turbulent jets in crossflow with a multigrid method and a second-moment closure model

A multigrid method is presented for calculating turbulent jets in crossflow. Fairly rapid convergence is obtained with the k-epsilon turbulence model, but computations with a full Reynolds stress turbulence model (RSM) are not yet very efficient. Grid dependency tests show that there are slight differences between results obtained on the two finest grid levels. Computations using the RSM are significantly different from those with k-epsilon model and compare better to experimental data. Some work is still required to improve the efficiency of the computations with the RSM.

Demuren, A. O.↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

Implementation of the Baldwin-Barth turbulence model into the ZETA code and its diagnosis

The Baldwin-Barth turbulence model was implemented into Zeta, a time-accurate, zonal, integro-differential code for incompressible laminar and turbulent flows. The implementation procedure is patterned after the model subroutine in ARC2D. The results of ZETA with the Baldwin-Barth turbulence model were compared with experimental data, with ZETA using Baldwin-Lomax model, and with ARC2D using the Baldwin-Barth model. The Baldwin-Barth model subroutine was tested by inputting an ARC2D velocity solution of an NACA-0012 airfoil at R(sub e) = 3.9 x 10(exp 6) and alpha = 5 deg. The resultant turbulent viscosity and Reynolds stresses compared favorably with the original data. For the same grid having grid points inside the laminar sublayer, which is necessary due to the one-equation nature of the model, ZETA however predicts early separation. It was found that the current ZETA has problem with such a fine grid. Further work is in progress to solve this problem.

Low, Scott L.↗

Shear Band Formation in Plane Strain Experiments of Sand

A series of biaxial (plane strain) experiments were conducted on three sands under low (15 kPa) and high (100 kPa) confining pressure conditions to investigate the effects of specimen density, confining pressure, and sand grain size and shape on the constitutive and stability behavior of granular materials. The three sands used in the experiments were fine-, medium-, and coarse-grained uniform silica sands with rounded, subangular, and angular grains, respectively. Specimen deformation was readily monitored and analyzed with the help of a grid pattern imprinted on the latex membrane. The overall stress-strain behavior is strongly dependent on the specimen density, confining pressure, sand grain texture, and the resulting failure mode(s). That became evident in different degrees of softening responses at various axial strains. The relationship between the constitutive behavior and the specimens' modes of instability is presented. The failure in all specimens was characterized by two distinct and opposite shear bands. It was found that the measured dilatancy angles increase as the sand grains' angularities and sizes increase. The measured shear band inclination angles are also presented and compared with classical Coulomb and Roscoe solutions.

Alshibli, Khalid A.↗

Assessment of Constitutive and Stability Behavior of Sands Under Plane Strain Condition

A series of biaxial (plane strain) experiments were conducted on three sands under low (15 kPa) and high (100 kPa) confining pressure conditions to investigate the effects of specimen density, confining pressure, and sand grains size and shape on the constitutive and stability behavior of granular materials. The three sands used in the experiments were fine, medium, and coarse-grained uniform silica sands with rounded, sub-angular, and angular grains, respectively. Specimen deformation was readily monitored and analyzed with the help of a grid pattern imprinted on the latex membrane. The overall stress-strain behavior is strongly dependent on the specimen density, confining pressure, sand grain texture, and the resulting failure mode(s). That became evident in different degrees of softening responses at various axial strains. The relationship between the constitutive behavior and the specimens' modes of instability is presented. The failure in all specimens was characterized by two distinct and opposite shear bands. It was found that the measured dilatancy angles increase as the sand grains' angularity and size increase. The measured shear band inclination angles are also presented and compared with classical Coulomb and Roscoe solutions.

Alshibli, Khalid A.↗

Numerical Investigation of Flow Over Periodic Hill

Hybrid turbulence modeling strategies that blend between large-eddy simulation and Reynolds-averaged Navier-Stokes promise to make optimal use of the available grid resolution. A delayed detached eddy simulation hybrid turbulence modeling strategy, based on the shear-stress transport turbulence model, was implemented into a finite-volume compressible Navier-Stokes code. Then the flow through a two-dimensional periodic hill channel was simulated. This geometry is a widely accepted benchmark test case for validating hybrid turbulence modeling strategies. The computational grid for the simulation was downloaded from a hybrid turbulence model validation website. The present hybrid turbulence model simulation predicts the size of the separated flow region as well as velocity and Reynolds stress profiles in good agreement with large-eddy simulation reference data.

Computational Fluid Dynamics↗

NASTRAN modifications for recovering strains and curvatures

Modifications to the NASTRAN structural analysis computer program are described. The modifications allow the recovery of strain and curvature data for the general two-dimensional elements, in addition to the usual stress data. Option features allow the transformation of the strain/curvature (or stress) data to a common coordinate system and representation at the grid points of the structural model rather than at the conventional element center locations. Usage information is provided which will allow present users of NASTRAN to easily utilize the new capability.

Hennrich, C. W.↗

Overflow Juncture Flow Computations Compared with Experimental Data

NASA’s Transformational Tools and Technologies Project (T(cu.) is supporting a substantial effort to investigate the formation and origin of separation bubbles found on wing-body juncture zones: the Juncture Flow experiment. The first phase of the Juncture Flow experiment, performed in NASA Langley’s 14- by 22-Foot Subsonic Tunnel, has been completed. This paper documents the CFD analysis done in conjunction with the experiment. Comparisons between CFD simulations and wind tunnel experimental results will be shown. Oil flow results, surface pressure cuts, velocity profiles, and Reynolds stress profiles will be compared. Preliminary results analyzing the effect of grid resolution, wind tunnel walls, and turbulence models on the above data will be presented. The results are not meant to be a validation study, but more as an evaluation of the current status of Reynolds averaged Navier Stokes (RANS) CFD simulations.

Lee, Henry C.↗

Shape optimization including finite element grid adaptation

The prediction of optimal shape design for structures depends on having a sufficient level of precision in the computation of structural response. These requirements become critical in situations where the region to be designed includes stress concentrations or unilateral contact surfaces, for example. In the approach to shape optimization discussed here, a means to obtain grid adaptation is incorporated into the finite element procedures. This facility makes it possible to maintain a level of quality in the computational estimate of response that is surely adequate for the shape design problem.

Kikuchi, N.↗

On modelling the pressure-strain correlations in wall bounded flows

Turbulence models for the pressure-strain term of the Reynolds-stress equations in the vicinity of a moving wall are evaluated for a high Reynolds number flow using decaying grid turbulence as a model problem. The data of Thomas and Hancock are used as a base for evaluating the different turbulence models. In particular, the Rotta model for return-to-isotropy is evaluated both in its inclusion into the Reynolds-stress equation model and in comparison to a nonlinear model advanced by Sarkar and Speziale. Further, models for the wall correction to the transfer term advanced by Launder et al., Shir, and Shih and Lumley are compared. Initial data using the decaying grid turbulence experiment as a base suggests that the coefficients proposed for these models are high perhaps by as much as an order of magnitude. The Shih and Lumley model which satisfies realizability constraints, in particular, seems to hold promise in adequately modeling the Reynolds stress components of this flow. Extensions of this work are to include testing the homogeneous transfer model by Shih and Lumley and the testing of the wall transfer models using their proposed coefficients and the coefficients chosen from this work in a flow with mean shear component.

Peltier, L. J.↗

An Advanced User Interface Approach for Complex Parameter Study Process Specification in the Information Power Grid

The creation of parameter study suites has recently become a more challenging problem as the parameter studies have now become multi-tiered and the computational environment has become a supercomputer grid. The parameter spaces are vast, the individual problem sizes are getting larger, and researchers are now seeking to combine several successive stages of parameterization and computation. Simultaneously, grid-based computing offers great resource opportunity but at the expense of great difficulty of use. We present an approach to this problem which stresses intuitive visual design tools for parameter study creation and complex process specification, and also offers programming-free access to grid-based supercomputer resources and process automation.

Yarrow, Maurice↗

Upgrades of Two Computer Codes for Analysis of Turbomachinery

Major upgrades have been made in two of the programs reported in "ive Computer Codes for Analysis of Turbomachinery". The affected programs are: Swift -- a code for three-dimensional (3D) multiblock analysis; and TCGRID, which generates a 3D grid used with Swift. Originally utilizing only a central-differencing scheme for numerical solution, Swift was augmented by addition of two upwind schemes that give greater accuracy but take more computing time. Other improvements in Swift include addition of a shear-stress-transport turbulence model for better prediction of adverse pressure gradients, addition of an H-grid capability for flexibility in modeling flows in pumps and ducts, and modification to enable simultaneous modeling of hub and tip clearances. Improvements in TCGRID include modifications to enable generation of grids for more complicated flow paths and addition of an option to generate grids compatible with the ADPAC code used at NASA and in industry. For both codes, new test cases were developed and documentation was updated. Both codes were converted to Fortran 90, with dynamic memory allocation. Both codes were also modified for ease of use in both UNIX and Windows operating systems.

Chima, Rodrick V.↗

Compensation for use of monthly-averaged winds in numerical modeling

Ratios R of the monthly averaged wind speeds to the magnitudes of the monthly averaged wind vectors are presented over a 41 x 41 grid covering the southern Ocean and the Antarctic continent. The ratio is found to vary from 1 to over 1000, with an average value of 1.86. These ratios R are relevant for converting from sensible and latent heats calculated with mean monthly data to those calculated with 12 hourly data. The corresponding ratios alpha for wind stress, along with the angle deviations involved, are also presented over the same 41 x 41 grid. The values of alpha generally exceed those for R and average 2.66. Regions in zones of variable wind directions have larger R and alpha ratios, over the ice-covered portions of the southern Ocean averaging 2.74 and 4.35 for R and alpha respectively. Thus adjustments to compensate for the use of mean monthly wind velocities should be stronger for wind stress than for turbulent heats and stronger over ice covered regions than over regions with more persistent wind directions, e.g., those in the belt of mid-latitude westerlies.

Parkinson, C. L.↗

Computation of turbulent boundary layers employing the defect wall-function method

In order to decrease overall computational time requirements of spatially-marching parabolized Navier-Stokes finite-difference computer code when applied to turbulent fluid flow, a wall-function methodology, originally proposed by R. Barnwell, was implemented. This numerical effort increases computational speed and calculates reasonably accurate wall shear stress spatial distributions and boundary-layer profiles. Since the wall shear stress is analytically determined from the wall-function model, the computational grid near the wall is not required to spatially resolve the laminar-viscous sublayer. Consequently, a substantially increased computational integration step size is achieved resulting in a considerable decrease in net computational time. This wall-function technique is demonstrated for adiabatic flat plate test cases from Mach 2 to Mach 8. These test cases are analytically verified employing: (1) Eckert reference method solutions, (2) experimental turbulent boundary-layer data of Mabey, and (3) finite-difference computational code solutions with fully resolved laminar-viscous sublayers. Additionally, results have been obtained for two pressure-gradient cases: (1) an adiabatic expansion corner and (2) an adiabatic compression corner.

Brown, Douglas L.↗