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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Interaction

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, N. V.↗
NASA NTRS2013

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Migration

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, Nelson V.↗
NASA NTRS1991

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulence models. The difficulty lies in the fact that the reaction rate is in general an exponential function of the temperature, and the higher order correlations in the conventional moment closure models of the chemical source term can not be neglected, making the applications of such models impractical. The probability density function (pdf) method offers an attractive alternative: in a pdf model, the chemical source terms are closed and do not require additional models. A grid dependent Monte Carlo scheme was studied, since it is a logical alternative, wherein the number of computer operations increases only linearly with the increase of number of independent variables, as compared to the exponential increase in a conventional finite difference scheme. A new algorithm was devised that satisfies a restriction in the case of pure diffusion or uniform flow problems. Although for nonuniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗
NASA NTRS1991

Current collection and current closure in the Tethered Satellite System

Current collection and closure-path modeling are examined analytically with respect to the Tethered Satellite System (TSS). A particle-in cell code is compared with a one-dimensional unmagnetized fluid code to model the behavior of a positively charged satellite in the ionosphere. The morphology of the sheath and the sheath-region processes are thus examined, and the influence of ions leaving the sheath region is found to cause the attraction of an electron current that is 40 times greater than the steady state value. The enhancement is transient and enhances the acceleration of the electrons in the sheath. A set of modified MHD equations, including those for ion inertia, quasineutrality, and electron drift, is employed to model TSS current closure. Whistler modes are found to exist and can be excited as the TSS passes through the ionosphere. Important conclusions include a significant fluctuation level in the steady state sheath, an ion void which affects the electron population, and some long-lived electrons trapped in the settled sheath with respect to two directions.

Drobot, Adam↗
NASA NTRS2010

Toward Realistic Simulation of low-Level Clouds Using a Multiscale Modeling Framework With a Third-Order Turbulence Closure in its Cloud-Resolving Model Component

This study presents preliminary results from a multiscale modeling framework (MMF) with an advanced third-order turbulence closure in its cloud-resolving model (CRM) component. In the original MMF, the Community Atmosphere Model (CAM3.5) is used as the host general circulation model (GCM), and the System for Atmospheric Modeling with a first-order turbulence closure is used as the CRM for representing cloud processes in each grid box of the GCM. The results of annual and seasonal means and diurnal variability are compared between the modified and original MMFs and the CAM3.5. The global distributions of low-level cloud amounts and precipitation and the amounts of low-level clouds in the subtropics and middle-level clouds in mid-latitude storm track regions in the modified MMF show substantial improvement relative to the original MMF when both are compared to observations. Some improvements can also be seen in the diurnal variability of precipitation.

Xu, Kuan-Man↗
NASA NTRS1994

Comparison of PDF and Moment Closure Methods in the Modeling of Turbulent Reacting Flows

In modeling turbulent reactive flows, Probability Density Function (PDF) methods have an advantage over the more traditional moment closure schemes in that the PDF formulation treats the chemical reaction source terms exactly, while moment closure methods are required to model the mean reaction rate. The common model used is the laminar chemistry approximation, where the effects of turbulence on the reaction are assumed negligible. For flows with low turbulence levels and fast chemistry, the difference between the two methods can be expected to be small. However for flows with finite rate chemistry and high turbulence levels, significant errors can be expected in the moment closure method. In this paper, the ability of the PDF method and the moment closure scheme to accurately model a turbulent reacting flow is tested. To accomplish this, both schemes were used to model a CO/H2/N2- air piloted diffusion flame near extinction. Identical thermochemistry, turbulence models, initial conditions and boundary conditions are employed to ensure a consistent comparison can be made. The results of the two methods are compared to experimental data as well as to each other. The comparison reveals that the PDF method provides good agreement with the experimental data, while the moment closure scheme incorrectly shows a broad, laminar-like flame structure.

Norris, Andrew T.↗
NASA NTRS1992

A mapping closure for turbulent scalar mixing using a time-evolving reference field

A general mapping-closure approach for modeling scalar mixing in homogeneous turbulence is developed. This approach is different from the previous methods in that the reference field also evolves according to the same equations as the physical scalar field. The use of a time-evolving Gaussian reference field results in a model that is similar to the mapping closure model of Pope (1991), which is based on the methodology of Chen et al. (1989). Both models yield identical relationships between the scalar variance and higher-order moments, which are in good agreement with heat conduction simulation data and can be consistent with any type of epsilon(phi) evolution. The present methodology can be extended to any reference field whose behavior is known. The possibility of a beta-pdf reference field is explored. The shortcomings of the mapping closure methods are discussed, and the limit at which the mapping becomes invalid is identified.

Girimaji, Sharath S.↗
NASA NTRS1988

Reynolds stress closure in jet flows using wave models

Ways of implementing the turbulence closure scheme based on modeling the large scale coherent structures as instability waves were sought. The computational tools necessary to apply this scheme to jets of arbitrary geometry were developed. The model, developed earlier, was extended to the shock structure of supersonic jets of arbitrary geometry and multiple jets. It was found that though the qualititate features of the unsteady flow field could be predicted there were always difficulties with some of the quantitative features. This led to the new formation of the closure scheme. The schemes for computations tools which were developed are efficient and represent the application of the very powerful mathematical tools to the problems of practical significance.

Morris, P. J.↗
NASA NTRS1985

A comparison of three algebraic stress closures for combustor flow calculations

A comparison is made of the performance of two locally nonequilibrium and one equilibrium algebraic stress closures in calculating combustor flows. Effects of four different pressure-strain models on these closure models are also analyzed. The results show that the pressure-strain models have a much greater influence on the calculated mean velocity and turbulence field than the algebraic stress closures, and that the best mean strain model for the pressure-strain terms is that proposed by Launder, Reece and Rodi (1975). However, the equilibrium algebraic stress closure with the Rotta return-to-isotropy model (1951) for the pressure-strain terms gives as good a correlation with measurements as when the Launder et al. mean strain model is included in the pressure-strain model. Finally, comparison of the calculations with the standard k-epsilon closure results show that the algebraic stress closures are better suited for simple turbulent flow calculations.

Nikjooy, M.↗
NASA NTRS2001

Turbulence Modeling Effects on the Prediction of Equilibrium States of Buoyant Shear Flows

The effects of turbulence modeling on the prediction of equilibrium states of turbulent buoyant shear flows were investigated. The velocity field models used include a two-equation closure, a Reynolds-stress closure assuming two different pressure-strain models and three different dissipation rate tensor models. As for the thermal field closure models, two different pressure-scrambling models and nine different temperature variance dissipation rate, Epsilon(0) equations were considered. The emphasis of this paper is focused on the effects of the Epsilon(0)-equation, of the dissipation rate models, of the pressure-strain models and of the pressure-scrambling models on the prediction of the approach to equilibrium turbulence. Equilibrium turbulence is defined by the time rate (if change of the scaled Reynolds stress anisotropic tensor and heat flux vector becoming zero. These conditions lead to the equilibrium state parameters. Calculations show that the Epsilon(0)-equation has a significant effect on the prediction of the approach to equilibrium turbulence. For a particular Epsilon(0)-equation, all velocity closure models considered give an equilibrium state if anisotropic dissipation is accounted for in one form or another in the dissipation rate tensor or in the Epsilon(0)-equation. It is further found that the models considered for the pressure-strain tensor and the pressure-scrambling vector have little or no effect on the prediction of the approach to equilibrium turbulence.

Zhao, C. Y.↗
NASA NTRS1997

Direct Aerosol Radiative Forcing: Calculations and Measurements from the Tropospheric

Radiative forcing is defined as the change in the net (downwelling minus upwelling) radiative flux at a given level in the atmosphere. This net flux is the radiative power density available to drive climatic processes in the earth-atmosphere system below that level. Recent research shows that radiative forcing by aerosol particles is a major source of uncertainty in climate predictions. To reduce those uncertainties, TARFOX was designed to determine direct (cloud-free) radiative forcing by the aerosols in one of the world's major industrial pollution plumes--that flowing from the east coast of the US over the Atlantic Ocean. TARFOX measured a variety of aerosol radiative effects (including direct forcing) while simultaneously measuring the chemical, physical, and optical properties of the aerosol particles causing those effects. The resulting data sets permit a wide variety of tests of the consistency, or closure, among the measurements and the models that link them. Because climate predictions use the same or similar model components, closure tests help to assess and reduce prediction uncertainties. In this work we use the TARFOX-determined aerosol, gas, and surface properties to compute radiative forcing for a variety of aerosol episodes, with inadvisable optical depths ranging from 0.07 to 0.6. We calculate forcing by several techniques with varying degrees of sophistication, in part to test the range of applicability of simplified techniques--which are often the only ones feasible in climate predictions by general circulation models (GCMs). We then compare computed forcing to that determined from: (1) Upwelling and downwelling fluxes (0.3-0.7 mm and 0.7-3.0 mm) measured by radiometers on the UK MRF C-130. and (2) Daily average cloud-free absorbed solar and emitted thermal radiative flux at the top of the atmosphere derived from the AVHRR radiometer on the NOAA- 14 satellite. The calculations and measurements all yield aerosol direct radiative forcing in the range -50 to -190 W sq m per unit inadvisable optical depth. The magnitudes are about 15 to 100 times larger than the global-average direct forcing expected for the global-average sulfate aerosol optical depth of 0.04. The reasons for the larger forcing in TARFOX include the relatively large optical depths and the focus on cloud-free, daytime conditions over the dark ocean surface. These are the conditions that produce the actual major radiative forcing events that contribute to any global-average climate effect. Detailed comparisons of calculated and measured forcings for specific events are used for more refined tests of closure.

Russell, P. B.↗
NASA NTRS1991

Experiments and modeling

The usual approach in establishing the correctness and accuracy of turbulence models is to numerically solve the modeled differential equations and then compare the results with the experiment. However, in the case of a discrepancy, this procedure does not pinpoint where in the model the drawback lies. It is also possible that the model overcompensates one physical phenomenon and undercompensates the other so that the net result is a good agreement between the two. Therefore, a more desirable approach is to directly compare the individual terms in the equations with their models. To achieve this objective, primary physical experiments were used to carry out the second moment budgets. These can then be used to analyze and assess various models and closure assumptions and seek improvements/modifications where models prove deficient.

Shabbir, Aamir↗
NASA NTRS2015

How Difficult is it to Reduce Low-Level Cloud Biases With the Higher-Order Turbulence Closure Approach in Climate Models?

Low-level clouds cover nearly half of the Earth and play a critical role in regulating the energy and hydrological cycle. Despite the fact that a great effort has been put to advance the modeling and observational capability in recent years, low-level clouds remains one of the largest uncertainties in the projection of future climate change. Low-level cloud feedbacks dominate the uncertainty in the total cloud feedback in climate sensitivity and projection studies. These clouds are notoriously difficult to simulate in climate models due to its complicated interactions with aerosols, cloud microphysics, boundary-layer turbulence and cloud dynamics. The biases in both low cloud coverage/water content and cloud radiative effects (CREs) remain large. A simultaneous reduction in both cloud and CRE biases remains elusive. This presentation first reviews the effort of implementing the higher-order turbulence closure (HOC) approach to representing subgrid-scale turbulence and low-level cloud processes in climate models. There are two HOCs that have been implemented in climate models. They differ in how many three-order moments are used. The CLUBB are implemented in both CAM5 and GDFL models, which are compared with IPHOC that is implemented in CAM5 by our group. IPHOC uses three third-order moments while CLUBB only uses one third-order moment while both use a joint double-Gaussian distribution to represent the subgrid-scale variability. Despite that HOC is more physically consistent and produces more realistic low-cloud geographic distributions and transitions between cumulus and stratocumulus regimes, GCMs with traditional cloud parameterizations outperform in CREs because tuning of this type of models is more extensively performed than those with HOCs. We perform several tuning experiments with CAM5 implemented with IPHOC in an attempt to produce the nearly balanced global radiative budgets without deteriorating the low-cloud simulation. One of the issues in CAM5-IPHOC is that cloud water content is much higher than in CAM5, which is combined with higher low-cloud coverage to produce larger shortwave CREs in some low-cloud prevailing regions. Thus, the cloud-radiative feedbacks are exaggerated there. The turning exercise is focused on microphysical parameters, which are also commonly used for tuning in climate models. The results will be discussed in this presentation.

Xu, Kuan-Man↗
NASA NTRS1978

Closure of a kinetic model of plasma in strong turbulence by relaxation

A Fokker-Planck kinetic equation for a turbulent plasma is derived by a repeated cascade decomposition. Calculation of the propagator and the kinetic equation determine the transport coefficients (diffusivity and turbulent viscosity) by means of a closure based on a relaxation procedure governing the approach to equilibrium. The k to the minus third power spectral law is obtained, which governs the coupling between the velocity and the electrostatic field fluctuations.

Tchen, C. M.↗
NASA NTRS1980

Progress in turbulence modeling for complex flow fields including effects of compressibility

Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.

Wilcox, D. C.↗
NASA NTRS1977

Turbulent transport model of a thunderstorm gust front

A second-order closure turbulent model is adopted to predict complex fluid flows; few a priori assumptions about the internal turbulent and mean flow structures are required. The flow simulations are applied to the analysis of local gust fronts generated by rain-cooled outflows from severe thunderstorms. Results of the simulations provide a means for discussing the behavior of an aircraft passing through a thunderstorm-related downdraft.

Teske, M. E.↗
NASA NTRS1992

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗