SEARCH · Engineering Papers

Results for “finite scale equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗

NASA NTRS1991

Numerical analysis of ignition transients in solid rocket motors

A model to analyze the unstready multidimensional turbulent flow in the head-end, star slot section of a solid rocket motor during the ignition transient is developed. The present paper examines the fluid dynamic aspects of the starting transient. A finite-difference solution of the unsteady, compressible, full Navier-Stokes equations, together with a two-equation k-epsilon turbulence model, is obtained utilizing the MacCormack explicit, predictor-corrector technique. Computed results for the flow in the slot are compared with experimental data obtained in cold flow tests of a scaled model of the Space Shuttle SRM head-end, star grain section. The agreement between the calculated flowfield and existing experimental data is very good.

Ciucci, A.↗

NASA NTRS1997

Large-Eddy Simulation of a High Reynolds Number Flow Around a Cylinder Including Aeroacoustic Predictions

The dynamic subgrid-scale model is employed in large-eddy simulations of flow over a cylinder at a Reynolds number, based on the diameter of the cylinder, of 90,000. The Centric SPECTRUM(trademark) finite element solver is used for the analysis. The far field sound pressure is calculated from Lighthill-Curle's equation using the computed fluctuating pressure at the surface of the cylinder. The sound pressure level at a location 35 diameters away from the cylinder and at an angle of 90 deg with respect to the wake's downstream axis was found to have a peak value of approximately 110 db. Slightly smaller peak values were predicted at the 60 deg and 120 deg locations. A grid refinement study suggests that the dynamic model demands mesh refinement beyond that used here.

Spyropoulos, Evangelos T.↗

NASA NTRS2008

Direct-View Multi-Point Two-Component Interferometric Rayleigh Scattering Velocimeter

This paper describes an instantaneous velocity measurement system based on the Doppler shift of elastically scattered laser light from gas molecules (Rayleigh scattering) relative to an incident laser. The system uses a pulsed laser as the light source, direct-viewing optics to collect the scattered light, an interferometer to analyze spectrally the scattered light mixed with the incident laser light, and a CCD camera to capture the resulting interferogram. The system is capable of simultaneous, spatially (approximately 0.2 mm(exp 3)) and temporally (approximately 40 ns) resolved, multiple point measurements of two orthogonal components of flow velocity in the presence of background scattered light, acoustic noise and vibrations, and flow particulates. Measurements in a large-scale axi-symmetric Mach 1.6 H2-air combustion-heated jet running at a flow sensible enthalpy specific to Mach 5.5 hypersonic flight are performed to demonstrate the technique. The measurements are compared with CFD calculations using a finite-volume discretization of the Favre-averaged Navier-Stokes equations (VULCAN code).

Bivolaru, Daniel↗

NASA NTRS1987

Experiments on the applicability of MAE techniques for predicting sound diffraction by irregular terrains

The sound field diffracted by a single smooth hill of finite impedance is studied both analytically, within the context of the theory of Matched Asymptotic Expansions (MAE), and experimentally, under laboratory scale modeling conditions. Special attention is given to the sound field on the diffracting surface and throughout the transition region between the illuminated and the shadow zones. The MAE theory yields integral equations that are amenable to numerical computations. Experimental results are obtained with a spark source producing a pulse of 42 microsec duration and about 130 Pa at 1 m. The insertion loss of the hill is inferred from measurements of the acoustic signals at two locations in the field, with subsequent Fourier analysis on an IBM PC/AT. In general, experimental results support the predictions of the MAE theory, and provide a basis for the analysis of more complicated geometries.

Berthelot, Yves H.↗

NASA NTRSDate not supplied

Understanding the Role of CFD, Including RANS, DES, LES and DNS, in Hypersonic Propulsion Design and Analysis

Computational fluid dynamics (CFD) has matured to the point where it is becoming ubiquitous in the design and analysis of aerospace vehicles and their propulsion systems. The current state-of-the-art of CFD in hypersonic propulsion design and analysis as well as a historical perspective of its development over the last 50 years is presented. The current state-of-the-art for the design of finite rate chemically reacting turbulent hypersonic flow is shown to currently be Reynolds averaged simulations (RAS) where the Reynolds averaged Navier-Stokes (RAS) equations are solved to obtain steady state solutions. An overview of the physical modeling techniques employed by RAS for the simulation of hypersonic chemically reacting flow are presented. However, scale resolving simulation techniques, such as detached eddy simulation (DES), large eddy simulations (LES), hybrid RAS/LES and direct numerical simulations (DNS), are becoming becoming powerful tools for point design analysis. One example of such a scale resolving technique, hybrid RAS/LES, will be briefly described as well as how it can be used to help address model form uncertainties in RAS.

Hypersonic↗

NASA NTRS1980

An analysis of turbulent diffusion flame in axisymmetric jet

The kinetic theory of turbulent flow was employed to study the mixing limited combustion of hydrogen in axisymmetric jets. The integro-differential equations in two spatial and three velocity coordinates describing the combustion were reduced to a set of hyperbolic partial differential equations in the two spatial coordinates by a binodal approximation. The MacCormick's finite difference method was then employed for solution. The flame length was longer than that predicted by the flame-sheet analysis, and was found to be in general agreement with a recent experimental result. Increase of the turbulence energy and scale resulted in an enhancement of the combustion rate and, hence, in a shorter flame length. Details of the numerical method as well as of the physical findings are discussed.

Chung, P. M.↗

NASA NTRS1989

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, Dimitri J.↗

NASA NTRS1990

Multigrid solution of the Navier-Stokes equations on triangular meshes

Mavriplis, D. J.↗

NASA NTRS1982

Entrainment and detrainment in a simple cumulus cloud model

A cumulus cloud model, analogous to the mixed-layer models of the planetary boundary layer and the upper ocean, is developed using a single, unitary entrainment process in which the motion of the cloud boundary relative to the mean flow is permitted, produced, and controlled by turbulent processes. An alternate theory to the mixing-length theory of Asai and Kasahara (1967) is proposed which completely removes the strong scale-dependence of the Asai-Kasahara model. The model reintroduces scale-dependence by introducing including the pe5turbation pressure term of the equation of vertical motion. It is shown that the model predicts deeper clouds than the Asai-Kasahara model for a given sounding, due to the entrainment assumption and the effects of the perturbation pressure. Lateral entrainment dominates cloud-top entrainment, although finite-difference errors increase the cloud-top entrainment rate from zero to a positive value in actual situations. The fractional entrainment rate for updrafts is determined to vary only slightly with height and to decrease only slowly as the cloud radius increases, while the fractional detrainment rate for updrafts increases with height.

Randall, D. A.↗

NASA NTRS1987

Computation of turbulent flows using an extended k-epsilon turbulence closure model

An extended kappa-epsilon turbulence model is proposed and tested with successful results. An improved transport equation for the rate of dissipation of the turbulent kinetic energy, epsilon, is proposed. The proposed model gives more effective response to the energy production rate than does the standard kappa-epsilon turbulence model. An extra time scale of the production range is included in the dissipation rate equation. This enables the present model to perform equally well for several turbulent flows with different characteristics, e.g., plane and axisymmetric jets, turbulent boundary layer flow, turbulent flow over a backward-facing step, and a confined turbulent swirling flow. A second-order accurate finite difference boundary layer code and a nearly second-order accurate finite difference elliptic flow solver are used for the present numerical computations.

Chen, Y.-S.↗

NASA NTRS1988

Dynamic structural correlation via nonlinear programming techniques

A solution to the correlation between structural dynamic test results and finite element analyses of the same components is presented in this paper. Basically, the method can be categorized as a Levenberg-Marquardt type Gauss-Newton method which requires only the differences between FE modal analyses and test results and their first derivatives with respect to preassigned design variables. With proper variable normalization and equation scaling, the method has been made numerically better-conditioned and the inclusion of the Levenberg-Marquardt technique overcomes any remaining difficulty encountered in inverting singular or near-singular matrices. An important feature is that each iteration requires only one function evaluation along with the associated design sensitivity analysis and so the procedure is computationally efficient.

Ting, T.↗

NASA NTRS1994

Unstructured-grid large-eddy simulation of flow over an airfoil

Historically, large-eddy simulations (LES) have been restricted to simple geometries where spectral or finite difference methods have dominated due to their efficient use of structured grids. Structured grids, however, not only difficulty representing complex domains and adapting to complicated flow features, but also are rather inefficient for simulating flows at high Reynolds numbers. The lack of efficiency stems from the need to resolve the viscous sublayer, which requires very fine resolution in all three directions near the wall. Structured grids make use of a stretching to reduce the normal grid spacing but must carry the fine resolution in the streamwise and spanwise directions throughout the domain. The unnecessarily fine grid for much of the domain leads to disturbingly high grid estimates. Chapman (1979), and later Moin & Jimenez (1993), pointed out that, in order to advance the technology to airfoils at flight Reynolds numbers, structured grids must be abandoned in lieu of what are known as nested or unstructured grids. The finite element method can efficiently solve the Navier-Stokes equations on unstructured grids. Although the CPU cost per time step per element is somewhat higher than structured grid methods, this effect is more than offset by the reduction in the number of elements. The use of unstructured grids, coupled with the advances in dynamic subgrid-scale modeling such as those made by Germano et al. (1991) and Ghosal et al. (1994), make LES of an airfoil tractable. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation since this flow has not been successfully simulated with the Reynolds-averaged Navier-Stokes equations.

Jansen, Kenneth↗

NASA NTRS1988

A finite element computation of turbulent boundary layer flows with an algebraic stress turbulence model

An algebraic stress turbulence model and a computational procedure for turbulent boundary layer flows which is based on the semidiscrete Galerkin FEM are discussed. In the algebraic stress turbulence model, the eddy viscosity expression is obtained from the Reynolds stress turbulence model, and the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale. Good agreement with experimental data is found for the examples of a fully developed channel flow, a fully developed pipe flow, a flat plate boundary layer flow, a plane jet exhausting into a moving stream, a circular jet exhausting into a moving stream, and a wall jet flow.

Kim, Sang-Wook↗

NASA NTRS1988

The meridional scale of baroclinic waves with latent heat release

The control on the meridional scale of a class of the baroclinic waves exercised by latent heat release is analyzed. A meridional-scale equation is derived, in which the dry model and the moist model without the meridional variation of the baroclinic waves are revealed. It is shown that, in the dry model the stability analysis cannot determine the meridional scale of the baroclinic waves. When latent heat release is included, the meridional variation of the waves either vanishes or is finite. When the waves have the meridional variation with latent heat release, the growth rate increases as the heating increases for a given Froude number, and there are two modes - the first mode has a small ascending region and a large descending region, while the second mode has a small descending region and a large ascending region.

Tang, Chung-Muh↗

NASA NTRS1979

Sound production due to large-scale coherent structures

The acoustic pressure fluctuations due to large-scale finite amplitude disturbances in a free turbulent shear flow are calculated. The flow is decomposed into three component scales; the mean motion, the large-scale wave-like disturbance, and the small-scale random turbulence. The effect of the large-scale structure on the flow is isolated by applying both a spatial and phase average on the governing differential equations and by initially taking the small-scale turbulence to be in energetic equilibrium with the mean flow. The subsequent temporal evolution of the flow is computed from global energetic rate equations for the different component scales. Lighthill's theory is then applied to the region with the flowfield as the source and an observer located outside the flowfield in a region of uniform velocity. Since the time history of all flow variables is known, a minimum of simplifying assumptions for the Lighthill stress tensor is required, including no far-field approximations. A phase average is used to isolate the pressure fluctuations due to the large-scale structure, and also to isolate the dynamic process responsible. Variation of mean square pressure with distance from the source is computed to determine the acoustic far-field location and decay rate, and, in addition, spectra at various acoustic field locations are computed and analyzed. Also included are the effects of varying the growth and decay of the large-scale disturbance on the sound produced.

Gatski, T. B.↗

NASA NTRS1987

Re-analysis procedure based on the mixed formulation

An analysis procedure is presented for large-scale structural systems, with large numbers of degrees of freedom and design variables. The procedure uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. Other elements of the procedure include: (1) lumping the design variables into a single tracing parameter; (2) operator splitting or additive decomposition of different arrays in the finite element equations into the corresponding arrays of the original structure plus correction terms; and (3) application of a reduction method through the use of the finite element method and the classical Bubnov-Galerkin technique. The re-analysis procedure is applied to the linear static and free vibration problems of plate and shell structures.

Noor, Ahmed K.↗

NASA NTRS2013

Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates

Progressive damage and failure analysis (PDFA) tools are needed to predict the nonlinear response of advanced fiber-reinforced composite structures. Predictive tools should incorporate the underlying physics of the damage and failure mechanisms observed in the composite, and should utilize as few input parameters as possible. The purpose of the Enhanced Schapery Theory (EST) was to create a PDFA tool that operates in conjunction with a commercially available finite element (FE) code (Abaqus). The tool captures the physics of the damage and failure mechanisms that result in the nonlinear behavior of the material, and the failure methodology employed yields numerical results that are relatively insensitive to changes in the FE mesh. The EST code is written in Fortran and compiled into a static library that is linked to Abaqus. A Fortran Abaqus UMAT material subroutine is used to facilitate the communication between Abaqus and EST. A clear distinction between damage and failure is imposed. Damage mechanisms result in pre-peak nonlinearity in the stress strain curve. Four internal state variables (ISVs) are utilized to control the damage and failure degradation. All damage is said to result from matrix microdamage, and a single ISV marks the micro-damage evolution as it is used to degrade the transverse and shear moduli of the lamina using a set of experimentally obtainable matrix microdamage functions. Three separate failure ISVs are used to incorporate failure due to fiber breakage, mode I matrix cracking, and mode II matrix cracking. Failure initiation is determined using a failure criterion, and the evolution of these ISVs is controlled by a set of traction-separation laws. The traction separation laws are postulated such that the area under the curves is equal to the fracture toughness of the material associated with the corresponding failure mechanism. A characteristic finite element length is used to transform the traction-separation laws into stress-strain laws. The ISV evolution equations are derived in a thermodynamically consistent manner by invoking the stationary principle on the total work of the system with respect to each ISV. A novel feature is the inclusion of both pre-peak damage and appropriately scaled, post-peak strain softening failure. Also, the characteristic elements used in the failure degradation scheme are calculated using the element nodal coordinates, rather than simply the square root of the area of the element.

Pineda, Evan J.↗