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 415 records · Page 23

The growth of small cracks in 4340 steel

A semiempirical crack-growth model incorporating crack-closure effects has been used to predict crack growth rates and total fatigue life in notched specimens of 4340 steel having an initial defect typical of those previously noted in this steel upon subjection to both constant-amplitude and spectrum loadings. Model predictions for fatigue life are found to be in reasonable agreement with experiments and, in particular, the model was able to accurately predict total fatigue life for specimens subjected to the 'Felix/28' spectrum loading sequence on the sole basis of constant-amplitude test test results.

Swain, M. H.↗

Modeling Fatigue Damage Onset and Progression in Composites Using an Element-Based Virtual Crack Closure Technique Combined With the Floating Node Method

A new methodology is proposed to model the onset and propagation of matrix cracks and delaminations in carbon-epoxy composites subject to fatigue loading. An extended interface element, based on the Floating Node Method, is developed to represent delaminations and matrix cracks explicitly in a mesh independent fashion. Crack propagation is determined using an element-based Virtual Crack Closure Technique approach to determine mixed-mode energy release rates, and the Paris-Law relationship to obtain crack growth rate. Crack onset is determined using a stressbased onset criterion coupled with a stress vs. cycle curve and Palmgren-Miner rule to account for fatigue damage accumulation. The approach is implemented in Abaqus/Standard® via the user subroutine functionality. Verification exercises are performed to assess the accuracy and correct implementation of the approach. Finally, it was demonstrated that this approach captured the differences in failure morphology in fatigue for two laminates of identical stiffness, but with layups containing θdeg plies that were either stacked in a single group, or distributed through the laminate thickness.

De Carvalho, Nelson V.↗

Analysis of Highly-Resolved Simulations of 2-D Humps Toward Improvement of Second-Moment Closures

Fully resolved simulation data of flow separation over 2-D humps has been used to analyze the modeling terms in second-moment closures of the Reynolds-averaged Navier- Stokes equations. Existing models for the pressure-strain and dissipation terms have been analyzed using a priori calculations. All pressure-strain models are incorrect in the high-strain region near separation, although a better match is observed downstream, well into the separated-flow region. Near-wall inhomogeneity causes pressure-strain models to predict incorrect signs for the normal components close to the wall. In a posteriori computations, full Reynolds stress and explicit algebraic Reynolds stress models predict the separation point with varying degrees of success. However, as with one- and two-equation models, the separation bubble size is invariably over-predicted.

Jeyapaul, Elbert↗

Numerical analysis and FORTRAN program for the computation of the turbulent wakes of turbomachinery rotor blades, isolated airfoils and cascade of airfoils

Turbulent wakes of turbomachinery rotor blades, isolated airfoils, and a cascade of airfoils were investigated both numerically and experimentally. Low subsonic and incompressible wake flows were examined. A finite difference procedure was employed in the numerical analysis utilizing the continuity, momentum, and turbulence closure equations in the rotating, curvilinear, and nonorthogonal coordinate system. A nonorthogonal curvilinear coordinate system was developed to improve the accuracy and efficiency of the numerical calculation. Three turbulence models were employed to obtain closure of the governing equations. The first model was comprised to transport equations for the turbulent kinetic energy and the rate of energy dissipation, and the second and third models were comprised of equations for the rate of turbulent kinetic energy dissipation and Reynolds stresses, respectively. The second model handles the convection and diffusion terms in the Reynolds stress transport equation collectively, while the third model handles them individually. The numerical results demonstrate that the second and third models provide accurate predictions, but the computer time and memory storage can be considerably saved with the second model.

Hah, C.↗

The Closure of the Ocean Mixed Layer Temperature Budget using Level-Coordinate Model Fields

Entrainment is an important element of the mixed layer mass, heat, and temperature budgets. Conventional procedures to estimate entrainment heat advection often do not permit the closure of heat and temperature budgets because of inaccuracies in its formulation. In this study a rigorous approach to evaluate the effect of entrainment using the output of a general circulation model (GCM) that does not have an explicit prognostic mixed layer model is described. The integral elements of the evaluation are 1) the rigorous estimates of the temperature difference between mixed layer water and entrained water at each horizontal grid point, 2) the formulation of the temperature difference such that the budget closes over a volume greater than one horizontal grid point, and 3) the apparent warming of the mixed layer during the mixed layer shoaling to account for the weak vertical temperature gradient within the mixed layer. This evaluation of entrainment heat advection is compared with the estimates by other commonly used ad hoc formulations by applying them in three regions: the north-central Pacific, the Kuroshio Extension, and the Nino-3 areas in the tropical Pacific. In all three areas the imbalance in the mixed layer temperature budget by the ad hoc estimates is significant, reaching a maximum of about 4 K yr(exp -1).

temperature gradient↗

Reheating the intergalactic medium under extremal conditions

Sherman (1982) commenced the investigation of a hypothetical closure density intergalactic medium (IGM) in the recombination epoch and found present acceptable temperatures in the range from 160,000 to 4,200,000 K for IGM heating. The present investigation has the objective to complete these calculations by finding extreme bounds on the reheating epoch which preceded recombination. The final results are presented in a table. Three models have been generated. The presented models suggest that a closure density IGM is constrained to be heated by a source in the range from 10 to the -31st to 10 to the -27th ergs per s per cu cm. The heating had to originate no earlier than z approximately 28.

Sherman, R. D.↗

A review of Reynolds stress models for turbulent shear flows

A detailed review of recent developments in Reynolds stress modeling for incompressible turbulent shear flows is provided. The mathematical foundations of both two-equation models and full second-order closures are explored in depth. It is shown how these models can be systematically derived for two-dimensional mean turbulent flows that are close to equilibrium. A variety of examples are provided to demonstrate how well properly calibrated versions of these models perform for such flows. However, substantial problems remain for the description of more complex turbulent flows where there are large departures from equilibrium. Recent efforts to extend Reynolds stress models to nonequilibrium turbulent flows are discussed briefly along with the major modeling issues relevant to practical naval hydrodynamics applications.

Speziale, Charles G.↗

Non-Equilibrium Turbulence and Two-Equation Modeling

Two-equation turbulence models are analyzed from the perspective of spectral closure theories. Kolmogorov theory provides useful information for models, but it is limited to equilibrium conditions in which the energy spectrum has relaxed to a steady state consistent with the forcing at large scales; it does not describe transient evolution between such states. Transient evolution is necessarily through nonequilibrium states, which can only be found from a theory of turbulence evolution, such as one provided by a spectral closure. When the departure from equilibrium is small, perturbation theory can be used to approximate the evolution by a two-equation model. The perturbation theory also gives explicit conditions under which this model can be valid, and when it will fail. Implications of the non-equilibrium corrections for the classic Tennekes-Lumley balance in the dissipation rate equation are drawn: it is possible to establish both the cancellation of the leading order Re1/2 divergent contributions to vortex stretching and enstrophy destruction, and the existence of a nonzero difference which is finite in the limit of infinite Reynolds number.

Rubinstein, Robert↗

Combustion research in the Internal Fluid Mechanics Division

The goal of this research is to bring computational fluid dynamics to a state of practical application for the aircraft engine industry. The approach is to have a strongly integrated computational and experimental program for all the disciplines associated with the gas turbine and other aeropropulsion systems by advancing the understanding of flow physics, heat transfer, and combustion processes. The computational and experimental research is integrated in the following way: the experiments that are performed provide an empirical data set so that physical models can be formulated to describe the processes that are occurring - for example, turbulence or chemical reaction. These experiments also form a data base for those who are doing code development by providing experimental data against which the codes can be verified and assesed. Models are generated as closure to some of the numerical codes, and they also provide physical insight for experiments. At the same time, codes which solve the complete Navier-Stokes equations can be used as a kind of numerical experiment from which far more extensive data can be obtained than ever could be obtained experimentally. This could provide physical insight into the complex processes that are taking place. These codes are also exercised against experimental data to assess the accuracy and applicability of models.

Mularz, Edward J.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Multiple-Time-Scale Concepts in Turbulent Transport Modeling

Progress in developing a turbulence closure employing two or more independently calculated time scales with which to characterize the rates of progress of different turbulent interactions is reported. The energy containing part of the spectrum is divided into two regions which respond at different rates and in different ways to changes in the environment. The scheme may be regarded as providing an intermediate level of approximation between the relatively simple, but fallible, single-point closures and the vastly more elaborate two-point closures which have so far been applied only to simulating homogeneous flows. The proposed approach requires only slightly more computational effort than single-scale schemes. 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.↗

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗

Space ecosynthesis: An approach to the design of closed ecosystems for use in space

The use of closed ecological systems for the regeneration of wastes, air, and water is discussed. It is concluded that such systems, if they are to be used for the support of humans in space, will require extensive mechanical and physico-chemical support. The reason for this is that the buffering capacity available in small systems is inadequate, and that natural biological and physical regulatory mechanisms rapidly become inoperative. It is proposed that mathematical models of the dynamics of a closed ecological system may provide the best means of studying the initial problems of ecosystem closure. A conceptual and mathematical model of a closed ecosystem is described which treats the biological components as a farm, calculates the rates of flow of elements through the system by mass-balance techniques and control theory postulates, and can evaluate the requirements for mechanical buffering activities. It is suggested that study of the closure of ecosystems can significantly aid in the establishment of general principles of ecological systems.

Macelroy, R. D.↗

Alfven Wave Reflection Model of Field-Aligned Currents at Mercury

An Alfven Wave Reflection (AWR) model is proposed that provides closure for strong field-aligned currents (FACs) driven by the magnetopause reconnection in the magnetospheres of planets having no significant ionospheric and surface electrical conductance. The model is based on properties of the Alfven waves, generated at high altitudes and reflected from the low-conductivity surface of the planet. When magnetospheric convection is very slow, the incident and reflected Alfven waves propagate along approximately the same path. In this case, the net field-aligned currents will be small. However, as the convection speed increases. the reflected wave is displaced relatively to the incident wave so that the incident and reflected waves no longer compensate each other. In this case, the net field-aligned current may be large despite the lack of significant ionospheric and surface conductivity. Our estimate shows that for typical solar wind conditions at Mercury, the magnitude of Region 1-type FACs in Mercury's magnetosphere may reach hundreds of kilo-Amperes. This AWR model of field-aligned currents may provide a solution to the long-standing problem of the closure of FACs in the Mercury's magnetosphere. c2009 Elsevier Inc. All rights reserved.

Lyatsky, Wladislaw↗

Computation of separated transonic turbulent flows

The two-dimensional Reynolds-averaged compressible Navier-Stokes equations are solved using MacCormack's second-order-accurate explicit finite difference method to simulate the separated transonic turbulent flow field over an airfoil. Four different algebraic eddy viscosity models are tested for viability to achieve turbulence closure for the class of flows considered. These models range from an unmodified boundary-layer mixing-length model to a relaxation model incorporating special considerations for the separation bubble region. Results of this study indicate the necessity for special attention to the separated flow region and suggest limits of applicability of algebraic turbulence models to these separated flow fields.

Deiwert, G. S.↗

Numerical experiments in homogeneous turbulence

The direct simulation methods developed by Orszag and Patternson (1972) for isotropic turbulence were extended to homogeneous turbulence in an incompressible fluid subjected to uniform deformation or rotation. The results of simulations for irrotational strain (plane and axisymmetric), shear, rotation, and relaxation toward isotropy following axisymmetric strain are compared with linear theory and experimental data. Emphasis is placed on the shear flow because of its importance and because of the availability of accurate and detailed experimental data. The computed results are used to assess the accuracy of two popular models used in the closure of the Reynolds-stress equations. Data from a variety of the computed fields and the details of the numerical methods used in the simulation are also presented.

Rogallo, R. S.↗

Direct numerical simulation of reacting flows

The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Riley, J. J.↗

Intraplate deformation and closure of the Australia-Antarctica-Africa plate circuit

Plate motion data along the Southeast, Southwest, and Central Indian ridges have been reduced to 67 spreading rates, 38 transform fault azimuths, and 135 earthquake slip vectors in order to study the current motion between the Australian, Antarctic, and African plates and to investigate whether this plate circuit obeys closure. Magnetic profiles are modeled to determine rates consistently over a 3-m.y. time-averaging interval, and the new rates are shown to differ from published rates by as much as 5 mm/yr. The results indicate that Indian Ocean plate circuit nonclosure and the deformation that it suggests are much smaller than previously supposed, and support a model in which the significant deformation occurs in a diffuse plate boundary along the equatorial Indian Ocean between the Central Indian Ridge and the Sumatra Trench.

Demets, Charles↗