Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure correlations”

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 127 records · Page 7

Troctolite 76535 - A study in the preservation of early isotopic records

The lunar rock considered in the present investigation is a coarse-grained troctolite granulite containing about 58(vol)% plagioclase, 37% olivine, 4% pyroxene, and less than 1% accessory phases with a texture which indicates formation as a cumulate at depths between 10 and 30 km followed by an extended period of slow cooling. A description is presented of noble gas studies of separated minerals from 76535. The quantity of fission xenon from the in situ decay of Pu-244 provides further evidence for different, mineral-specific, isotopic closure times. The presented data shows that 76535 loses its surface-correlated xenon component upon disaggregation. No other xenon component is lost. The presence of solar gases in 76535 would seem to argue in favor of the external acquisition of the parentless extinct isotope effects and consequently favor 'thermal diffusion' and 'adsorption' over local redistribution models.

Caffee, M.↗

Effect of Stability on Mixing in Open Canopies

In open canopies, the within-canopy flux from the ground surface and understory can account for a significant fraction of the total flux above the canopy. This study incorporates the important influence of within-canopy stability on turbulent mixing and subcanopy fluxes into a first-order closure scheme. Toward this goal, we analyze within-canopy eddy-correlation data from the old aspen site in the Boreal Ecosystem - Atmosphere Study (BOREAS) and a mature ponderosa pine site in Central Oregon, USA. A formulation of within-canopy transport is framed in terms of a stability- dependent mixing length, which approaches Monin-Obukhov similarity theory above the canopy roughness sublayer. The new simple formulation is an improvement upon the usual neglect of the influence of within-canopy stability in simple models. However, frequent well-defined cold air drainage within the pine subcanopy inversion reduces the utility of simple models for nocturnal transport. Other shortcomings of the formulation are discussed.

Lee, Young-Hee↗

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.↗

Decay of homogeneous turbulence from a specified state

The homogeneous turbulence problem is formulated by first specifying the multipoint velocity correlations or their spectral equivalents at an initial time. Those quantities, together with the correlation or spectral equations, are then used to calculate initial time derivatives of correlations or spectra. The derivatives in turn are used in time series to calculate the evolution of turbulence quantities with time. When the problem is treated in this way, the correlation equations are closed by the initial specification of the turbulence and no closure assumption is necessary. An exponential series which is an iterative solution of the Navier stokes equations gave much better results than a Taylor power series when used with the limited available initial data. In general, the agreement between theory and experiment was good.

Deissler, R. G.↗

Finite element analysis of aeroacoustic jet-flap flows

A computational analysis was performed on the steady, turbulent aerodynamic flowfields associated with a jet-blown flap. For regions devoid of flow separation, a parabolic approximation to the governing time-averaged Navier-Stokes equations was applied. Numerical results are presented for the symmetry plane flow of a slot-nozzle planar jet flap geometry, including prediction of flowfield evolution within the secondary mixing region immediately downstream of the trailing edge. Using a two equation turbulence kinetic energy closure model, rapid generation and decay of large spatial gradients in mean and correlated fluctuating velocity components within the immediate wake region were predicted. Modifications to the turbulent flow structure, as induced by porous surface treatment of the flap, were evaluated. The recirculating flow within a representative discrete slot in the surface was evaluated, using the two dimensional, time-averaged Navier-Stokes equations.

Baker, A. J.↗

Effect of Liquid Droplets on Turbulence Structure in a Round Gaseous Jet

A second-order model which predicts the modulation of turbulence in jets laden with uniform size solid particles or liquid droplets is discussed. The approach followed is to start from the separate momentum and continuity equations of each phase and derive two new conservation equations. The first is for the carrier fluid's kinetic energy of turbulence and the second for the dissipation rate of that energy. Closure of the set of transport equations is achieved by modeling the turbulence correlations up to a third order. The coefficients (or constants) appearing in the modeled equations are then evaluated by comparing the predictions with LDA-measurements obtained recently in a turbulent jet laden with 200 microns solid particles. This set of constants is then used to predict the same jet flow but laden with 50 microns solid particles. The agreement with the measurement in this case is very good.

Elghobashi, S. E.↗

A two-equation turbulence model for two-phase flows

A two-equation turbulence model has been developed for predicting two-phase flows. The two equations describe the conservation of turbulence kinetic energy and dissipation rate of that energy for the carrier fluid in a two-phase flow. They have been derived rigorously from the momentum equations of the carrier fluid. Closure of the time-mean equations is achieved by modeling the turbulent correlations up to third order. The new model eliminates the need to simulate in an ad hoc manner the effects of the dispersed phase on turbulence structure. Preliminary testing indicates that the model is successful in predicting the main features of a round gaseous jet laden with uniform-size solid particles.

Elghobashi, S. E.↗

Advancements in engineering turbulence modeling

Some new developments in two-equation models and second order closure models are presented. Two-equation models (k-epsilon models) have been widely used in computational fluid dynamics (CFD) for engineering problems. Most of low-Reynolds number two-equation models contain some wall-distance damping functions to account for the effect of wall on turbulence. However, this often causes the confusion and difficulties in computing flows with complex geometry and also needs an ad hoc treatment near the separation and reattachment points. A set of modified two-equation models is proposed to remove the aforementioned shortcomings. The calculations using various two-equation models are compared with direct numerical simulations of channel flow and flat boundary layers. Development of a second order closure model is also discussed with emphasis on the modeling of pressure related correlation terms and dissipation rates in the second moment equations. All the existing models poorly predict the normal stresses near the wall and fail to predict the 3-D effect of mean flow on the turbulence (e.g. decrease in the shear stress caused by the cross flow in the boundary layer). The newly developed second order near-wall turbulence model is described and is capable of capturing the near-wall behavior of turbulence as well as the effect of 3-D mean flow on the turbulence.

Shih, T.-H.↗

Gear crack propagation investigations

Analytical and experimental studies were performed to investigate the effect of gear rim thickness on crack propagation life. The FRANC (FRacture ANalysis Code) computer program was used to simulate crack propagation. The FRANC program used principles of linear elastic fracture mechanics, finite element modeling, and a unique re-meshing scheme to determine crack tip stress distributions, estimate stress intensity factors, and model crack propagation. Various fatigue crack growth models were used to estimate crack propagation life based on the calculated stress intensity factors. Experimental tests were performed in a gear fatigue rig to validate predicted crack propagation results. Test gears were installed with special crack propagation gages in the tooth fillet region to measure bending fatigue crack growth. Good correlation between predicted and measured crack growth was achieved when the fatigue crack closure concept was introduced into the analysis. As the gear rim thickness decreased, the compressive cyclic stress in the gear tooth fillet region increased. This retarded crack growth and increased the number of crack propagation cycles to failure.

Lewicki, David G.↗

An Evaluation of the Plasticity-Induced Crack-Closure Concept and Measurement Methods

An assessment of the plasticity-induced crack-closure concept is made, in light of some of the questions that have been raised on the validity of the concept, and the assumptions that have been made concerning crack-dp damage below the crack-opening stress. The impact of using other crack-tip parameters, such as the cyclic crack-tip displacement, to model crack-growth rate behavior was studied. Crack-growth simulations, using a crack-closure model, showed a close relation between traditional Delta K eff, and the cyclic crack-tip displacement (Delta eff) for an aluminum alloy and a steel. Evaluations of the cyclic hysteresis energy demonstrated that the cyclic plastic damage below the crack-opening stress was negligible in the Paris crack-growth regime. Some of the standard and newly proposed remote measurement methods to determine the 'effective' crack-tip driving parameter were evaluated on middle-crack tension specimens. A potential source of the Kmax effect on crack-growth rates was studied on an aluminum alloy. Results showed that the ratio of Kmax to Kc had a strong effect on crack-growth rates at high stress ratios and at low stress ratios for very high stress levels. The crack-closure concept and the traditional crack-growth rate equations were able to correlate and predict crack-growth rates under these extreme conditions.

Newman, James C., Jr.↗

Dynamical System Analysis of Reynolds Stress Closure Equations

In this paper, we establish the causality between the model coefficients in the standard pressure-strain correlation model and the predicted equilibrium states for homogeneous turbulence. We accomplish this by performing a comprehensive fixed point analysis of the modeled Reynolds stress and dissipation rate equations. The results from this analysis will be very useful for developing improved pressure-strain correlation models to yield observed equilibrium behavior.

Girimaji, Sharath S.↗

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗

Charge Radius of the Short-Lived Ni 68 and Correlation with the Dipole Polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus 68 Ni at the N = 40 subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α$ D has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus 48 Ca. Three-particle–three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of 68 Ni.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hf-W chronology of ordinary chondrites

Fifteen H, L, and LL ordinary chondrites of petrologic types 4–6 have been analyzed for Hf-W isotope systematics to constrain the chronology, internal structure, and thermal history of their parent bodies. For most samples coarse-grained metals plot below the isochrons defined by silicate-dominated fractions which consist of variable mixtures of silicate minerals with tiny metal inclusions. This offset results from an earlier Hf-W closure in the large metal grains and provides a new means for simultaneously determining cooling rates and Hf-W closure ages for individual samples. For most type 5 and 6 samples, cooling rates and Hf-W ages are inversely correlated, indicating that these samples derive from concentrically zoned bodies in which more strongly metamorphosed samples derive from greater depth. These data, therefore, provide strong evidence for a common ‘onion shell’ structure for the H, L, and LL chondrite parent bodies. The cooling rates and Hf-W ages of some type 5 and 6 chondrites overlap, indicating that the Hf-W systematics provide a more robust measure of the thermal history and burial depth of a given sample than the simple petrographic distinction between types 5 and 6. Two type 6 samples deviate from the correlation between cooling rates and Hf-W ages and cooled much faster than expected for their Hf-W age. These samples likely were excavated by impacts that occurred during high-temperature metamorphism and prior to complete closure of the Hf-W system at ~10 Ma after CAI formation. As these impacts would have disturbed the asteroid’s cooling history, these samples likely derive from different bodies than samples with undisturbed cooling histories, implying that ordinary chondrites derive from more than just three parent bodies. The Hf-W data reveal that metal-silicate fractionation among the H, L, and LL groups occurred between ~2 and ~2.7 Ma after CAI formation and, hence, was about coeval to chondrule formation. As both metal-silicate fractionation and chondrule formation occurred prior to chondrite parent body accretion, there should be no ordinary chondrite chondrules that are younger than ~2.7 Ma. In conclusion, ordinary chondrite precursors had lower Hf/W ratios than carbonaceous chondrites, suggesting that inner and outer solar system materials, respectively, were chemically distinct even for refractory elements.

36 MATERIALS SCIENCE↗

Cyclic Permeability Evolution During Repose Then Reactivation of Fractures and Faults

Abstract Cyclic growth and decay of permeability in fractures is explored during repeated reactivation and repose of saw‐cut fractures of Green River shale. These slide‐hold‐slide experiments are supported by measurements of fracture normal deformation and optical surface profilometry. Overall, we observe continuous permeability decay during repose (holds) and significant permeability enhancement during slow reactivation (slide). The permeability decay is accompanied by fault compaction. Both hydraulic aperture change (Δ b h ) and measured compaction (Δ b s ) are consistent with time‐dependent power law closure with a power exponent of ~0.2–0.4. These dual compaction magnitudes are positively correlated but Δ b h > Δ b s in late stage holds. Permeability enhancement during reactivation is typically also accompanied by fault dilation. However, we also observe some cases where hydraulic aperture change decouples from the measured deformation, conceivably driven by mobilization of wear products and influenced by the development of flow bottlenecks. Pretest and posttest surface profiles show that the surface topography of the fractures is planed down by shear removal. The shear removal is significant with initial laboratory prepared surface (~10 μm of aperture height) but less significant following consecutive reactivations (~2 μm). The flattened surfaces retain small‐scale, ~10–20 μm wavelength, roughness. Flow simulations, constrained by the surface topography and measured deformation, indicate that small‐scale roughness may control permeability at flow bottlenecks within a dominant flow channel. These results suggest cycles of permeability creation and destruction are an intrinsic component of the natural hydraulic system present in faults and fractures and provide an improved mechanistic understanding of the evolution of permeability during fault repose and reactivation.

Im, Kyungjae↗

Risk-Informed MultiPphysics BestEestimate Plus Uncertainties (BEPU) Application Development of RELAP5-3D Perturbation Model

The United States nuclear industry is facing a strong challenge to ensure maximum safety while enhancing economic benefit. Safety is a key parameter to all aspects related to light water reactor (LWR) nuclear power plants (NPPs), especially cost savings. Since the goal is to extend the lifetimes of these NPPs, the traditional deterministic safety concept may not guarantee a current economic asset. The Light Water Reactor Sustainability (LWRS) Program has been promoting a wide range of research and development (R&D) in this field to maximize the safety, economics, and performance of these NPPs through improved scientific understanding. One of the best practices to achieve this goal is to identify and optimize safety margins, which can lead to cost reduction. To do this, under the LWRS framework, the Risk-Informed Systems Analysis (RISA) Pathway will focus on the optimization of safety margin and minimization of uncertainties to ensure both safety and economics at the highest level. The RISA Pathway will provide enhanced capabilities for analyzing and characterizing LWR systems performance by developing and demonstrating methods, tools, and data to enable risk-informed margins management (RIMM). The goals of the RISA Pathway are twofold: (1) deploy the risk-informed tools and methods that enable better representation of safety margins and factors that contribute to cost and safety; and (2) conduct advanced risk assessment applications with industry to support margin management strategies that enable more cost-effective plant operation. The tools and methods provided by the RISA Pathway will support effective margin management for both active and passive safety systems, structures, and components (SSC) of an NPP. The tools and methods used in the RISA Pathway should have high confidence and highest technical maturity for and implementation to industry at its current setting. They should also have a capability to support risk-informed decision making for both probabilistic and deterministic elements of safety. The RISA Pathway will, therefore, perform a comprehensive assessment of verification and validation (V&V) status of RISA Toolkit to enhance credibility RISA Toolkit which be used by industry. This report summarizes RELAP5-3D development activities for best estimated and uncertainty (BEPU) capability to support ongoing RISA Pathway pilot projects for risk informed uncertainty quantification applications. Work scope includes (1) development of initial perturbation model for selected closure laws and (2) testing of suitable probability distribution functions (PDF) for the selected correlations. Developed BEPU modules are tested with selected separate effect cases and compared with conventional thermal-hydraulics model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Solid motor aft closure insulation erosion

The erosion rate of aft closure insulation in a number of large solid propellant motors was empirically analyzed by correlating the average ablation rate with a number of variables that had previously been demonstrated to affect heat flux. The main correlating parameter was a heat flux based on the simplified Bartz heat transfer coefficient corrected for two-dimensional effects. A multiplying group contained terms related to port-to-throat ratio, local wall angle, grain geometry and nozzle cant angle. The resulting equation gave a good correlation and is a useful design tool.

Stampfl, E.↗

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗