Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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

Assimilation of scatterometer winds into surface pressure fields using a variational method

A variational formulation was used to assimilate Seasat-A scatterometer (SASS) surface wind measurements near and during a severe storm in the North Atlantic into conventional National Meteorological Center sea level pressure fields. An estimate of the relative vorticity at every point on a grid was calculated using each of these two data sets. A solution to a modified geostrophic stream function is found subject to the constraints that (1) the relative vorticities calculated from the data agree as closely as possible with the relative vorticities from the variational solution, and that (2) the average kinetic energy is a minimum. Results are obtained which support the idea that averaged satellite data can be treated as synoptic data. Direct substitution rather than a time-weighted insertion made from SASS winds generally resulted in more accurate pressure analyses. In addition, this relatively simple model provides surface pressure fields which agree extremely well with surface truth and the results of other investigators who required additional sources of input data into more complex models. It will be possible to obtain improved wind field maps from future scatterometer pressure fields in mid-latitudes.

Harlan, J., Jr.↗
NASA NTRS1989

Two-dimensional analysis of two-phase reacting flow in a firing direct-injection diesel engine

The flow field, spray penetration, and combustion in two-stroke diesel engines are described. Fuel injection begins at 345 degrees after top dead center (ATDC) and n-dodecane is used as the liquid fuel. Arrhenius kinetics is used to calculate the reaction rate term in the quasi-global combustion model. When the temperature, fuel, and oxygen mass fraction are within suitable flammability limits, combustion begins spontaneously. No spark is necessary to ignite a localized high temperature region. Compression is sufficient to increase the gaseous phase temperature to a point where spontaneous chemical reactions occur. Results are described for a swirl angle of 22.5 degrees.

Nguyen, H. Lee↗
NASA NTRS1990

Theoretical investigation of gas-surface interactions

A Dirac-Hartree-Fock code was developed for polyatomic molecules. The program uses integrals over symmetry-adapted real spherical harmonic Gaussian basis functions generated by a modification of the MOLECULE integrals program. A single Gaussian function is used for the nuclear charge distribution, to ensure proper boundary conditions at the nuclei. The Gaussian primitive functions are chosen to satisfy the kinetic balance condition. However, contracted functions which do not necessarily satisfy this condition may be used. The Fock matrix is constructed in the scalar basis and transformed to a jj-coupled 2-spinor basis before diagonalization. The program was tested against numerical results for atoms with a Gaussian nucleus and diatomic molecules with point nuclei. The energies converge on the numerical values as the basis set size is increased. Full use of molecular symmetry (restricted to D sub 2h and subgroups) is yet to be implemented.

Dyall, Kenneth G.↗
NASA NTRS1995

Nonaxisymmetric secular instabilities driven by star/disk coupling

We determine conditions for the onset of nonaxisymmetric secular instabilities in polytropes with a wide range of angular momentum distributions using Lagrangian techniques, and then calculate the growth rate of such instabilities when driven by the coupling of the perturbed star to a circumstellar disk. We use Langrangian displacement vectors with azimuthal dependence proportional to exp (im phi), where m is an integer and phi is the azimuthal coordinate. The onset of secular instability in terms of the quantity T/absolute value of W, the ratio of rotational kinetic energy to gravitational potential energy, is affected by both the compressibility and the angular momentum distribution of the polytrope. The largest effects occur when the angular momentum distribution is varied. For polytropic index n = 3/2, the onset of secular instability for the m = 2 mode (the bar mode), as determined by its neutral point, shifts from T/absolute value of W = 0.141 to 0.093, while the m = 5 mode neutral point shifts from T/absolute value of W = 0.088 to 0.031 over the range of angular momentum distributions we consider. The smallest critical T/absolute value of W-values occur for the angular momentum distributions which are the most peaked toward the equator. For the angular momentum distribution of a Maclaurin spheroid, as the polytropic index n is increased from 3/2 to 5/2, the neutral point for m = 2 shifts from T/absolute value of W = 0.141 to 0.144 and the netural point for m = 5 shifts from T/absolute value of W = 0.069 to 0.078. The netural points for m = 2 and 5 for the Maclaurin sequence (n = 0) are 0.137 and 0.0629, respectively. As the angular momentum distribution becomes more peaked toward the equatorial radius of the polytropes, the critical T/absolute value of W-values generally become less sensitive to the compressibility of the polytrope. Star/disk coupling can drive the secular instability in systems where the star is surrounded by a massive disk and, if the instability can grow to moderate amplitude, then the coupling can transport significant amounts of angular momentum from the star into the circumstellar disk. We find, for the particular case of rotating protostars during the accretion phase, that the instability growth time can be shorter than the accretion time. Further, if the instability can grow to amplitudes on the order of several percent, the star/disk coupling can remove angular momentum from the forming star faster than it is added by accretion.

Imamura, James, N.↗
DOE OSTI2024

Direct observation of a magnetic-field-induced Wigner crystal

Wigner predicted that when the Coulomb interactions between electrons become much stronger than their kinetic energy, electrons crystallize into a closely packed lattice. A variety of two-dimensional systems have shown evidence for Wigner crystals (WCs). However, a spontaneously formed classical or quantum WC has never been directly visualized. Neither the identification of the WC symmetry nor direct investigation of its melting has been accomplished. Here we use high-resolution scanning tunnelling microscopy measurements to directly image a magnetic-field-induced electron WC in Bernal-stacked bilayer graphene and examine its structural properties as a function of electron density, magnetic field and temperature. At high fields and the lowest temperature, we observe a triangular lattice electron WC in the lowest Landau level. The WC possesses the expected lattice constant and is robust between filling factor ν ≈ 0.13 and ν ≈ 0.38 except near fillings where it competes with fractional quantum Hall states. Increasing the density or temperature results in the melting of the WC into a liquid phase that is isotropic but has a modulated structure characterized by the Bragg wavevector of the WC. At low magnetic fields, the WC unexpectedly transitions into an anisotropic stripe phase, which has been commonly anticipated to form in higher Landau levels. Furthermore, analysis of individual lattice sites shows signatures that may be related to the quantum zero-point motion of electrons in the WC lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Pynta-An Automated Workflow for Calculation of Surface and Gas-Surface Kinetics

Many important industrial processes rely on heterogeneous catalytic systems. However, given all possible catalysts and conditions of interest, it is impractical to optimize most systems experimentally. Automatically generated microkinetic models can be used to efficiently consider many catalysts and conditions. However, these microkinetic models require accurate estimation of many thermochemical and kinetic parameters. Manually calculating these parameters is tedious and error prone, involving many interconnected computations. Here, we present Pynta, a workflow software for automating the calculation of surface and gas–surface reactions. Pynta takes the reactants, products, and atom maps for the reactions of interest, generates sets of initial guesses for all species and saddle points, runs all optimizations, frequency, and IRC calculations, and computes the associated thermochemistry and rate coefficients. It is able to consider all unique adsorption configurations for both adsorbates and saddle points, allowing it to handle high index surfaces and bidentate species. Pynta implements a new saddle point guess generation method called harmonically forced saddle point searching (HFSP). HFSP defines harmonic potentials based on the optimized adsorbate geometries and which bonds are breaking and forming that allow initial placements to be optimized using the GFN1-xTB semiempirical method to create reliable saddle point guesses. This method is reaction class agnostic and fast, allowing Pynta to consider all possible adsorbate site placements efficiently. We demonstrate Pynta on 11 diverse reactions involving monodenate, bidentate, and gas-phase species, many distinct reaction classes, and both a low and a high index facet of Cu. Our results suggest that it is very important to consider reactions between adsorbates adsorbed in all unique configurations for interadsorbate group transfers and reactions on high index surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
NASA NTRS1983

NASA redox storage system development project

The operating temperature was raised from 25 C to 65 C, which enhanced the kinetics of the chromium electrode charging reactions. The design of the auxiliary electrochemical cell, which is used to keep both reactants at the same state of charge, was modified, leading to better and more stable performance. Preliminary testing has shown that the four tank mode of operation improves energy efficiency as much as 5 percentage points over the conventional two tank mode. Another variation in operating mode, the use of mixed reactants, potentially offers several very attractive advantages. Preliminary reactant cost studies lend further weight to the feasibility of the mixed reactant concept. Electrode studies show that reproducibility of performance is very dependent on the pyrolysis temperature at which the carbon/graphite felt substrate is formed. Membrane development work continued to concentrate on cost reduction and the enhancement of resistivity and selectivity.

Source record↗
NASA NTRS2000

Ground Vibration Test Planning and Pre-Test Analysis for the X-33 Vehicle

This paper describes the results of the modal test planning and the pre-test analysis for the X-33 vehicle. The pre-test analysis included the selection of the target modes, selection of the sensor and shaker locations and the development of an accurate Test Analysis Model (TAM). For target mode selection, four techniques were considered, one based on the Modal Cost technique, one based on Balanced Singular Value technique, a technique known as the Root Sum Squared (RSS) method, and a Modal Kinetic Energy (MKE) approach. For selecting sensor locations, four techniques were also considered; one based on the Weighted Average Kinetic Energy (WAKE), one based on Guyan Reduction (GR), one emphasizing engineering judgment, and one based on an optimum sensor selection technique using Genetic Algorithm (GA) search technique combined with a criteria based on Hankel Singular Values (HSV's). For selecting shaker locations, four techniques were also considered; one based on the Weighted Average Driving Point Residue (WADPR), one based on engineering judgment and accessibility considerations, a frequency response method, and an optimum shaker location selection based on a GA search technique combined with a criteria based on HSV's. To evaluate the effectiveness of the proposed sensor and shaker locations for exciting the target modes, extensive numerical simulations were performed. Multivariate Mode Indicator Function (MMIF) was used to evaluate the effectiveness of each sensor & shaker set with respect to modal parameter identification. Several TAM reduction techniques were considered including, Guyan, IRS, Modal, and Hybrid. Based on a pre-test cross-orthogonality checks using various reduction techniques, a Hybrid TAM reduction technique was selected and was used for all three vehicle fuel level configurations.

Bedrossian, Herand↗
NASA NTRS2011

Differential Processing of Low and High LET Radiation Induced DNA Damage: Investigation of Switch from ATM to ATR Signaling

The members of the phosphatidylinositol kinase-like kinase family of proteins namely ataxia-telangiectasia mutated (ATM) and ATM- and Rad3-related (ATR) are directly responsible for the maintenance of genomic integrity by mounting DDR through signaling and facilitating the recruitment of repair factors at the sites of DNA damage along with coordinating the deployment of cell cycle checkpoints to permit repair by phosphorylating Checkpoint kinase Chk1, Chk2 and p53. High LET radiation from GCR (Galactic Cosmic Rays) consisting mainly of protons and high energy and charged (HZE) particles from SPE (Solar Particle Event) pose a major health risk for astronauts on their space flight missions. The determination of these risks and the design of potential safeguards require sound knowledge of the biological consequences of lesion induction and the capability of the cells to counter them. We here strive to determine the coordination of ATM and ATR kinases at the break sites directly affecting checkpoint signaling and DNA repair and whether differential processing of breaks induced by low and high LET radiation leads to possible augmentation of swap of these damage sensors at the sites of DNA damage. Exposure of cells to IR triggers rapid autophosphorylation of serine-1981 that causes dimer dissociation and initiates monomer formation of ATM. ATM kinase activity depends on the disruption of the dimer, which allows access and phosphorylation of downstream ATM substrates like Chk2. Evidence suggests that ATM is activated by the alterations in higher-order chromatin structure although direct binding of ATM to DSB ends may be a crucial step in its activation. On the other hand, in case of ATR, RPA (replication protein A)-coated ssDNA (single-stranded DNA) generated as a result of stalled DNA replication or during processing of chromosomal lesions is crucial for the localization of ATR to sites of DNA damage in association with ATR-interacting protein (ATRIP). Although the majority of RPA-coated ssDNA is generally present only during DNA replication, ATR activation in G1 and G2-phase might still require formation of RPA-coated ssDNA, probably initiated by the MRN-CtIP complex and then extended by the Exo1- or BLM-dependent mechanisms at the sites of DSBs. Evidence accumulates that activation of ATM and ATR are oppositely regulated by the length of single stranded overhangs generated at the break sites by processes mentioned above and these stretches of single stranded overhangs hold the clue for ATM to ATR switch at broken DNA ends. We irradiated 82-6hTERT human fibroblast cells with low LET gamma-rays and high LET Fe and Si particles. Preliminary results with cells exposed to 1Gy gamma-rays show that the kinetics of pChk2-pT68 foci formation is comparable to that of gamma-H2AX although they appear to recede quicker. The number and intensity of observed foci reaches a maximum at 30 min and 60 min post IR for Chk2-pT68 and gamma-H2AX foci respectively and all Chk2-pT68 foci colocalize with gamma-H2AX foci. The kinetics of Chk1-pS345 and ATRIP are being determined. Results of Chk2-pT68 foci kinetics was also corroborated by western blot experiments, although phosphorylation was detected as early as 10 min and started receding 30 min post IR with 2Gy of gamma-rays. On the other hand, level of ATR-pS428 reached its maximum between 60 and 120 min and was maintained until the last measured time point of 4 hours post IR as determined by western blotting. Experiments performed with high LET Fe and Si particles will be reported.

Saha, Janapriya↗
DOE OSTI2020

Reaction mechanisms of a cyclic ether intermediate: Ethyloxirane

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β -hydroperoxyalkyl radicals ( β -˙QOOH) or from reactions of HOÓ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n -butane oxidation. While rate coefficients for β -˙QOOH → ethyloxirane + ÓH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. Here we provide fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ˙R + O 2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ˙QOOH → cyclic ether + ÓH and ˙R + O 2 → conjugate alkene + HOÓ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ˙QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ˙QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1971

Gain kinetics of CO2 gasdynamic laser mixtures at high pressure.

Description of a combined analytical and experimental investigation of rapidly expanded CO2 laser mixtures. Study of the vibrational kinetics indicates that population inversions with high vibrational energy density can be produced at high pressures by utilizing low (starved) concentrations of CO2. Scaling laws for predicting the gain of starved systems are developed. A description is given of gain measurements carried out using a CO2 laser to probe the flow in a nozzle at the point where the Mach number is approximately four. Population inversions at static pressures as high as about 1 atm have been observed and found to be in reasonable agreement with predictions.

Christiansen, W. H.↗
NASA NTRS1977

Localized induced electric field within the magnetotail

The conversion of magnetic energy, stored within the magnetotail during the growth phase, into particle kinetic energy is studied by taking induced electric fields, due to explicit time dependence of magnetic fields, into consideration. The polarization of a plasma in the presence of such a field is discussed, and the topological difference between a polarized and an induced electric field, namely that one is irrotational and the other is not, is pointed out. A localized perturbation in neutral sheet current is also discussed.

Heikkila, W. J.↗
NASA NTRS1932

The Gaseous Explosive Reaction : the Effect of Pressure on the Rate of Propagation of the Reaction Zone and upon the Rate of Molecular Transformation

This study of gaseous explosive reaction has brought out a number of important fundamental characteristics of the explosive reaction indicating that the basal processes of the transformation are much simpler and corresponds more closely to the general laws and principles of ordinary transformations than is usually supposed. The report calls attention to the point that the rate of molecular transformation within the zone was found in all cases to be proportional to pressure, that the transformation within the zone is the result of binary impacts. This result is of unusual interest in the case of the reaction of heavy hydrocarbon fuels and the reaction mechanism proposed by the recent kinetic theory of chain reactions.

Stevens, F W↗
NASA NTRS2008

The New NASA Orbital Debris Mitigation Procedural Requirements and Standards

NASA has issued major updates to its principal orbital debris mitigation policy directive and standards. The new NASA Procedural Requirements for Limiting Orbital Debris (NPR 8715.6), with its supporting NASA Standard 8719.14, both refine earlier orbital debris mitigation documents and in some areas expand their applicability. Organizational and individual responsibilities along with general directives are set forth in NPR 8715.6. New requirements include routine conjunction assessments for all maneuverable NASA spacecraft in LEO and GEO, prompt notifications of intended or unintended debris generation, preparation and maintenance of formal end-of-mission plans, and disposal of vehicles in operation around the Moon and Mars and at the Earth-Sun Lagrangian points. NASA Standard 8719.14 replaces the 1995 NASA Safety Standard 1740.14 with no major new requirements but with several refinements and additions, some of which had already been adopted. Compliance with human casualty risk limitations from reentering debris will be calculated explicitly and not be expressed in terms of average debris casualty area. Moreover, the minimum kinetic energy threshold for potentially injurious reentering debris is set at 15 Joules. The overarching requirement for the disposal of GEO spacecraft and launch vehicle orbital stages is to ensure that the vehicles do not come within GEO + 200 km for at least 100 km after end of mission, rather than setting specific requirements for the disposal orbit. Spacecraft operating in or routinely transiting LEO must remain in the region for no more than 25 years after end of mission or 30 years after launch, whichever occurs sooner. A comprehensive new NASA handbook on orbital debris has also been prepared to provide background on the orbital debris environment and the related NASA mitigation requirements and standards.

Johnson, Nicholas L.↗
DOE OSTI2023

Characterizing and Engineering Trehalose Biosynthesis Pathways to Improve Acetate Tolerance in Pseudomonas putida KT2440

Alkaline pretreatment can depolymerize lignin into a diverse mixture of monomers and smaller compounds, enabling bioconversion of these compounds to value-added chemicals in engineered bacterial hosts such as Pseudomonas putida KT2440. Acetate and salts comprise approximately 8% and 32% of the lignin stream, respectively. High concentrations of these components can cause cell death in P. putida, so lignin-derived streams must be fed in lower amounts to decrease stress. Investigation of the mechanisms of bacterial tolerance to chemicals in lignin-rich substrates is therefore required to overcome this obstacle in lignin bioconversion. Previous work has shown that trehalose, an endogenously produced disaccharide, is involved in tolerance to heat, cold, and osmotic stress, as well as implicated in acid stress in other bacteria. As such, trehalose biosynthesis presents a useful target to increase acetate and osmotic tolerance in P. putida. This work investigates the native TreSA, TreSB, and TreY/TreZ trehalose biosynthesis pathways in this bacterial strain as well as the OtsAB pathway of Sphingobium sp. SYK-6. Overexpression of the OtsAB pathway in P. putida increased tolerance to acetate relative to wild-type. Additionally, functional knockouts of treSA, treSB, treY, and treZ in P. putida showed that interruption of any trehalose biosynthesis pathway, and, in the case of TreY/TreZ, interruption of the pathway at any point, significantly hindered growth in acetate as well as on glucose and LB broth. Each trehalose biosynthesis pathway was then engineered for overexpression in P. putida to determine whether overexpression of these genes conferred additional tolerance to acetate. The TreSA, TreSB, OtsA, and OtsB enzymes were also expressed and purified to determine the reaction rates and kinetics of these pathways. The hypothetical glycoside hydrolase family 15 protein, encoded in the gene located in between otsB and otsA in the Sphingobium operon, was purified as well to investigate its potential involvement in the OtsAB pathway. The findings of this work illuminated the mechanisms of trehalose biosynthesis in P. putida and identified genetic targets to improve strain tolerance to common components of alkaline-pretreated lignin streams. These findings could ultimately improve yields of desired products from lignin in engineered strains of P. putida.

alkaline pretreatment↗
NASA NTRS1985

Undercooling and solidification behavior in the InSb-Sb system

Use of the droplet emulsion technique has been successful in studying the undercooling and crystallization behavior of Sb, InSb, and an InSb-Sb eutectic alloy. Both droplet size and imposed cooling rate were influential in controlling the extent of liquid undercooling. The droplet surface coating was of significant importance in determining the resultant solidification product structure through its effect on nucleation kinetics. The maximum undercooling for pure Sb was extended from 0.08 to 0.23 T sub m. While simple crushing techniques provided a dramatic increase in droplet undercooling over the bulk material, emulsification treatments both enhanced this undercooling and allowed successful formation of a metastable simple cubic Sb phase. This phase was stable to temperatures approaching the melting point. The simple cubic phase was detected in droplet samples processed using DTA, air and water quenching, and drop tube processing under a helium gas atmosphere. A deviation in the InSb parent ingot composition limited interpretation of the line compound results, however, emulsification techniques extended the undercooling of this material to 0.17 T sub L and provided a stable, protective surface coating for the droplets. Emulsification of the eutectic alloy was effective at producing various levels of undercooling from 0.1 to 0.2 T sub E. Microstructural examination revealed a normal-type eutectic structure in the undercooled droplets indicating that solidification occurred within the coupled zone and that this zone is somewhat symmetric about the eutectic composition.

Graves, J. A.↗
NASA NTRS2001

Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation

Subgrid-scale (SGS) models for large. eddy simulation (LES) have generally been evaluated by their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. Recent application- of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem in which time, correlations determine the frequency distribution of acoustic radiation, suggest that subgrid models should also be evaluated by their ability to predict time correlations in turbulent flows. This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay more slowly than the corresponding time. correlations in the DNS fields. This observation is confirmed by theoretical estimates of time correlations using the Taylor expansion technique. Tile reason for the slower decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy viscosity associated with energy transfer is a leading order approximation.

He, Guowei↗