Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure model”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Sharply curved turn around duct flow predictions using spectral partitioning of the turbulent kinetic energy and a pressure modified wall law

Computational predictions of turbulent flow in sharply curved 180 degree turn around ducts are presented. The CNS2D computer code is used to solve the equations of motion for two-dimensional incompressible flows transformed to a nonorthogonal body-fitted coordinate system. This procedure incorporates the pressure velocity correction algorithm SIMPLE-C to iteratively solve a discretized form of the transformed equations. A multiple scale turbulence model based on simplified spectral partitioning is employed to obtain closure. Flow field predictions utilizing the multiple scale model are compared to features predicted by the traditional single scale k-epsilon model. Tuning parameter sensitivities of the multiple scale model applied to turn around duct flows are also determined. In addition, a wall function approach based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients is tested. Turn around duct flow characteristics utilizing this modified wall law are presented and compared to results based on a standard wall treatment.

Santi, L. Michael↗
NASA NTRS2000

Tar Production from Biomass Pyrolysis in a Fluidized Bed Reactor: A Novel Turbulent Multiphase Flow Formulation

A novel multiphase flow model is presented for describing the pyrolysis of biomass in a 'bubbling' fluidized bed reactor. The mixture of biomass and sand in a gaseous flow is conceptualized as a particulate phase composed of two classes interacting with the carrier gaseous flow. The solid biomass is composed of three initial species: cellulose, hemicellulose and lignin. From each of these initial species, two new solid species originate during pyrolysis: an 'active' species and a char, thus totaling seven solid-biomass species. The gas phase is composed of the original carrier gas (steam), tar and gas; the last two species originate from the volumetric pyrolysis reaction. The conservation equations are derived from the Boltzmann equations through ensemble averaging. Stresses in the gaseous phase are the sum of the Newtonian and Reynolds (turbulent) contributions. The particulate phase stresses are the sum of collisional and Reynolds contributions. Heat transfer between phases, and heat transfer between classes in the particulate phase is modeled, the last resulting from collisions between sand and biomass. Closure of the equations must be performed by modeling the Reynolds stresses for both phases. The results of a simplified version (first step) of the model are presented.

Bellan, J.↗
NASA NTRS1990

Application of a Reynolds stress turbulence model to the compressible shear layer

Theoretically based turbulence models have had success in predicting many features of incompressible, free shear layers. However, attempts to extend these models to the high-speed, compressible shear layer have been less effective. In the present work, the compressible shear layer was studied with a second-order turbulence closure, which initially used only variable density extensions of incompressible models for the Reynolds stress transport equation and the dissipation rate transport equation. The quasi-incompressible closure was unsuccessful; the predicted effect of the convective Mach number on the shear layer growth rate was significantly smaller than that observed in experiments. Having thus confirmed that compressibility effects have to be explicitly considered, a new model for the compressible dissipation was introduced into the closure. This model is based on a low Mach number, asymptotic analysis of the Navier-Stokes equations, and on direct numerical simulation of compressible, isotropic turbulence. The use of the new model for the compressible dissipation led to good agreement of the computed growth rates with the experimental data. Both the computations and the experiments indicate a dramatic reduction in the growth rate when the convective Mach number is increased. Experimental data on the normalized maximum turbulence intensities and shear stress also show a reduction with increasing Mach number.

Sarkar, S.↗
NASA NTRS1990

Application of a Reynolds stress turbulence model to the compressible shear layer

Theoretically based turbulence models have had success in predicting many features of incompressible, free shear layers. However, attempts to extend these models to the high-speed, compressible shear layer have been less effective. In the present work, the compressible shear layer was studied with a second-order turbulence closure, which initially used only variable density extensions of incompressible models for the Reynolds stress transport equation and the dissipation rate transport equation. The quasi-incompressible closure was unsuccessful; the predicted effect of the convective Mach number on the shear layer growth rate was significantly smaller than that observed in experiments. Having thus confirmed that compressibility effects have to be explicitly considered, a new model for the compressible dissipation was introduced into the closure. This model is based on a low Mach number, asymptotic analysis of the Navier-Stokes equations, and on direct numerical simulation of compressible, isotropic turbulence. The use of the new model for the compressible dissipation led to good agreement of the computed growth rates with the experimental data. Both the computations and the experiments indicate a dramatic reduction in the growth rate when the convective Mach number is increased. Experimental data on the normalized maximum turbulence intensities and shear stress also show a reduction with increasing Mach number.

Sarkar, S.↗
NASA NTRS1988

A numerical method for computing unsteady 2-D boundary layer flows

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.

Krainer, Andreas↗
NASA NTRSDate not supplied

Hybrid NEP-Chemical Vehicle and Propulsion Technology Study for Crewed Mars Missions

Nuclear Electric Propulsion (NEP) is a suite of technologies with the potential to enable crewed opposition-class missions to Mars. Each subsystem comprising an NEP system has multiple technology options that present tradeoffs in vehicle sizing and technology development risk. Multiple vehicle concepts have been studied in the past and are currently being studied at NASA. Three models have been developed to illustrate the impact of key technology parameters on overall sizing for a crewed Mars transportation vehicle. Vehicle mass required to close the mission as a function of NEP system alpha (kWe/kg) and specific impulse (Isp) is estimated with a detailed trajectory model. A power system mass model estimates alpha as a function of power and radiator mass assumption. A power conversion system thermodynamic model predicts the radiator area required to close a Brayton cycle power conversion system. Combined with the architecture mass required for mission closure, the radiator area and alpha model provide insight into how technology development may impact the mission. Our focus is on parametric sweeps of the whole design space rather than any particular point design. Results from the mission model comparing electric propulsion technology (Hall, Ion, Magnetoplasmadynamic) and both 2039 and 2042 opposition-class mission launch windows are included. The results of the study illustrate the technology parameters that can result in mission closure and illustrate where the architecture is sensitive to variations in technology performance or requirements.

Matthew Duchek↗
NASA NTRS2020

Juncture Flow Computations using kL-Based Turbulence Models

The development and implementation of kL-based Reynolds-Average Navier-Stokes (RANS) two-equation turbulence models are reported in this paper. The kL model is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta showed that a reliable transport equation can be formed from the turbulence length scale L, and the turbulence kinetic energy k. Rotta's kL equation is well suited for a term-by-term modeling and displays useful features compared to other scale formulations. One of the important differences is the inclusion of higher-order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The k-kL scheme has been modified to include the Algebraic Reynolds Stress Model (ARSM) and the Quadratic Constitutive Relation (QCR) as nonlinear models. These models have gone through extensive validations using two-dimensional benchmark problems. The present study documents the application of the k-kL schemes to simulate flow around the Juncture Flow Model (JFM). The k-kL prediction results show generally good comparisons with measurements. The results from this formulation are similar to, or better than results using the SA-RC-QCR2000 two-equation turbulence model. The ARSM and the QCR formulations of k-kL show promise with a similar level of computational resources as basic two-equation turbulence models.

K S Abdol-Hamid↗
NASA NTRS2004

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗
NASA NTRS1978

Dual-loop model of the human controller

A dual-loop model of the human controller in single-axis compensatory tracking tasks is introduced. This model possesses an inner-loop closure that involves feeding back that portion of controlled element output rate that is due to control activity. A novel feature of the model is the explicit appearance of the human's internal representation of the manipulator-controlled element dynamics in the inner loop. The sensor inputs to the human controller are assumed to be system error and control force. The former can be sensed via visual, aural, or tactile displays, whereas the latter is assumed to be sensed in kinesthetic fashion. A set of general adaptive characteristics for the model is hypothesized, including a method for selecting simplified internal models of the manipulator-controlled element dynamics. It is demonstrated that the model can produce controller describing functions that closely approximate those measured in four laboratory tracking tasks in which the controlled element dynamics vary considerably in terms of ease of control. An empirically derived expression for the normalized injected error remnant spectrum is introduced.

Hess, R. A.↗
NASA NTRSDate not supplied

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗
NASA NTRS1979

Multiple-time-scale concepts in turbulent transport modelling

The paper reports progress in developing a closure employing two or more independently calculated time scales with which to characterize the rates of progress of different turbulent interactions. The approach contrasts with that used by earlier single-point models which adopt just a single time scale, proportional to the turbulence energy turnover time. The present treatment divides the energy containing part of the spectrum into two regions which respond at different rates and in different ways to changes in the environment. Computational results are reported for several thin shear flows which show striking improvement in the level of agreement with experiment over that obtained with models employing only one time scale.

Hanjalic, K.↗
NASA NTRS1996

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗
NASA NTRSDate not supplied

Recent Advances in Discrete Crack Modeling Applied to Laminated Composites with Emphasis on: Floating Node Method, VCCT and Cohesive Zone Modeling

The present talk will provide an overview of the work performed during the Advanced Composites Project (ACP) on the development, and verification and validation of the Floating Node Method (FNM) as well as the Virtual Crack Closure Technique (VCCT) and cohesive zone modeling (CZM). The FNM is a finite element-based technique to represent crack networks. The complex nature of matrix crack-delamination interactions observed in unidirectional (UD) tape laminates suggests that a methodology such as the FNM method may be required to be able to accurately simulate damage progression in these laminates. Simulating crack onset and growth within the context of the FNM relies on techniques such as VCCT and CZM. The talk is organized to provide, via select examples, an overview of the breadth of the Verification & Validation (V&V) exercises performed during the ACP, and how these challenged the state-of-the art and guided further developments in discrete crack modeling, while helping to establish confidence in the progress made and map the challenges ahead. The performance of the VCCT and CZM individually, and in combination with the FNM can be assessed through verification exercises. These exercises typically consist of a comparison of simulation results to known numerical or analytical solutions. Verification is key to identify implementation issues and limitations that, otherwise, may remain undetected and cloud any subsequent validation efforts. Indeed, a subset of these numerical exercises led to further developments of the VCCT and the FNM method as will be illustrated. Before embarking on the subsequent validation of the framework, it is critical to have adequate characterization data. However, the testing campaign conducted revealed material responses that challenged the state-of-the-art and required further developments. The developments in CZM technology associated with the modeling of the responses of hybrid interfaces (fabric/UD) will be given as an example. Finally, the talk will conclude with a summary of the validation exercises performed under quasi-static and fatigue loadings, highlighting some of the key achievements, outstanding challenges and lessons learned.

finite elements↗
NASA NTRS1999

Mesoscale Simulations of a Florida Sea Breeze Using the PLACE Land Surface Model Coupled to a 1.5-Order Turbulence Parameterization

A sophisticated land-surface model, PLACE, the Parameterization for Land Atmospheric Convective Exchange, has been coupled to a 1.5-order turbulent kinetic energy (TKE) turbulence sub-model. Both have been incorporated into the Penn State/National Center for Atmospheric Research (PSU/NCAR) mesoscale model MM5. Such model improvements should have their greatest effect in conditions where surface contrasts dominate over dynamic processes, such as the simulation of warm-season, convective events. A validation study used the newly coupled model, MM5 TKE-PLACE, to simulate the evolution of Florida sea-breeze moist convection during the Convection and Precipitation Electrification Experiment (CaPE). Overall, eight simulations tested the sensitivity of the MM5 model to combinations of the new and default model physics, and initialization of soil moisture and temperature. The TKE-PLACE model produced more realistic surface sensible heat flux, lower biases for surface variables, more realistic rainfall, and cloud cover than the default model. Of the 8 simulations with different factors (i.e., model physics or initialization), TKE-PLACE compared very well when each simulation was ranked in terms of biases of the surface variables and rainfall, and percent and root mean square of cloud cover. A factor separation analysis showed that a successful simulation required the inclusion of a multi-layered, land surface soil vegetation model, realistic initial soil moisture, and higher order closure of the planetary boundary layer (PBL). These were needed to realistically model the effect of individual, joint, and synergistic contributions from the land surface and PBL on the CAPE sea-breeze, Lake Okeechobee lake breeze, and moist convection.

Lynn, Barry H.↗
NASA NTRS1992

Exploiting similarity in turbulent shear flows for turbulence modeling

It is well known that current k-epsilon models cannot predict the flow over a flat plate and its wake. In an effort to address this issue and other issues associated with turbulence closure, a new approach for turbulence modeling is proposed which exploits similarities in the flow field. Thus, if we consider the flow over a flat plate and its wake, then in addition to taking advantage of the log-law region, we can exploit the fact that the flow becomes self-similar in the far wake. This latter behavior makes it possible to cast the governing equations as a set of total differential equations. Solutions of this set and comparison with measured shear stress and velocity profiles yields the desired set of model constants. Such a set is, in general, different from other sets of model constants. The rational for such an approach is that if we can correctly model the flow over a flat plate and its far wake, then we can have a better chance of predicting the behavior in between. It is to be noted that the approach does not appeal, in any way, to the decay of homogeneous turbulence. This is because the asymptotic behavior of the flow under consideration is not representative of the decay of homogeneous turbulence.

Robinson, David F.↗
NASA NTRS1999

The Impact of Buoyancy and Flame Structure on Soot, Radiation and NOx Emissions from a Turbulent Diffusion Flame

It is hypothesized that the spatial structure of a turbulent diffusion flame plays an important role in determining the emissions of radiative energy, soot and NO, from a combustor. This structure, manifested in the two point statistics, is influenced by buoyancy. Radiation, soot and NOx emissions are the cumulative result of processes that occur throughout a flame. For example, radiation fluxes along a line of sight can be found from summing up the contributions from sources in individual pockets of hot soot that emit, and from sinks in cold soot that absorb. Soot and NOx are both the results of slow chemistry and are not equilibrium products. The time that is available for production and burnout is crucial in determining the eventual emissions of these pollutants. Turbulence models generally rely on a single point closure of the appropriate time averaged equations. Hence, spatial information is lost and needs to be modeled using solution variables such as turbulence kinetic energy and dissipation rate, often with the assumption of isotropy. However, buoyancy can affect the physical structure of turbulent flames and can change the spatial extent of soot bearing regions. Theoretical comparisons with models are best done in the limit of infinite Froude number because the inclusion of buoyancy in flow models introduces significant uncertainties. Hence, LII measurements of soot, measurements of radiation fluxes from soot, Particle Imaging Velocimetry (PIV) of the flow field and measurements of post flame NOX will be carried out on the NASA Lewis 2.2 sec drop tower and eventually on the parabolic flight aircraft. The drop rig will be a modified version of a unit that has been successfully used at Lewis in the past.

Kennedy, I. M.↗
NASA NTRS1992

Studies in turbulence

Various papers on turbulence are presented. Individual topics addressed include: modeling the dissipation rate in rotating turbulent flows, mapping closures for turbulent mixing and reaction, understanding turbulence in vortex dynamics, models for the structure and dynamics of near-wall turbulence, complexity of turbulence near a wall, proper orthogonal decomposition, propagating structures in wall-bounded turbulence flows. Also discussed are: constitutive relation in compressible turbulence, compressible turbulence and shock waves, direct simulation of compressible turbulence in a shear flow, structural genesis in wall-bounded turbulence flows, vortex lattice structure of turbulent shear slows, etiology of shear layer vortices, trilinear coordinates in fluid mechanics.

Gatski, Thomas B.↗
NASA NTRS1987

Low Reynolds number modeling of turbulent flows with and without wall transpiration

A full Reynolds-stress closure that is capable of describing the flow all the way to the wall is formulated. The closure is based on the conventional high Reynolds number form of the redistribution model, the inclusion of molecular diffusion, and a modified dissipation model to account for viscous effects near a wall. Two dissipation models are investigated along with two gradient diffusion and two redistribution models. Their respective effects on the calculated flow properties are assessed by comparing them with the data of fully developed turbulent flows and a developing pipe flow with wall transpiration. The near-wall behavior is very well predicted; however, the wall correction to the redistribution modeling is found to have little effect on the calculated results. The overall behavior of the fully developed turbulent flows is best described by a nonisotropic gradient diffusion model, a return-to-isotropy redistribution model, and a dissipation model that accounts for viscous behavior near a wall. This same closure also gives the best prediction of the axial pressure drop behavior along a pipe with a uniform wall suction. Furthermore, the near-wall behavior of such a flow is very well predicted by this closure.

So, Ronald M. C.↗