Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogenization theory”

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 667 records · Page 37

An atmospheric analysis of the carbon-rich white dwarf G35 - 26

The first detailed atmospheric analysis of the DQ white dwarf G35 - 26 is presented. The analysis is based on published spectra and on helium-rich, homogeneous, blanketed LTE models containing carbon, hydrogen, and small amounts of N, O, and Ca, and on the assumption that G35 - 26 is a single star. The observed spectra are fitted against the model spectra, and T(eff), log(g), and abundances of H, He, C are determined along with upper limits on N, O, and Ca. The mass and zero temperature Hamada-Salpeter (ZTHS) radius, distance, parallax, and tangential velocity are derived. The mass, 1.2-1.33 solar, is the highest yet observed for a single white dwarf if the star is on the ZTHS relation. O + Ne + Mg white dwarf in this mass range are predicted by theory and G35 - 26 may be the first such object discovered. Given an ZTHS radius, the distance is between 33 and 69 pc.

Thejll, Peter↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

Supported Single‐Atom Manganese Catalysts for the Trimerization of Ethylene

Selective ethylene oligomerization via oxidative cyclization, forming metallacyclic intermediates, is typically catalyzed by molecular titanium and chromium complexes to produce butenes, hexenes, or octenes, depending on the supporting ligand framework. However, this mechanism requires significant electron density at the metal active site and is not known to be generalizable to other first-row transition metals. In this work, we computationally investigate the electronic modulation of five transition metals (Mn, Fe, Co, Ni, and Cu) supported on titania (TiO₂) through reductive lithium intercalation to promote selective oligomerization via oxidative cyclization, using density functional theory (DFT). Our findings predict that Mn/LiTiO₂ exhibits high catalytic activity due to the exergonic nature of oxidative cyclization with two ethylene molecules. Additionally, lithium titanate (LiTiO₂) supports enhance catalytic performance compared to TiO₂. Experimental validation confirms that Mn/LiTiO₂ achieves higher conversion rates and improved selectivity toward hexene (C₄:C₆ = 1:2.6). The enhanced activity is attributed to lithiation, which alters the electronic environment around Mn active sites. Mechanistic studies reveal that the formation of a seven-membered ring, a key intermediate for hexene formation, is more favorable on LiTiO₂ than TiO₂. This work provides the first evidence of Mn catalyzing selective ethylene oligomerization via oxidative cyclization in either homogeneous or heterogeneous catalysis.

Kim, Yu Lim [Argonne National Laboratory (ANL), Ar↗

Non-equilibrium one-dimensional two-phase flow in variable area channels

A one-dimensional nonequilibrium flow analysis has been formulated for a one component two phase flow. The flow is considered homogeneous and essentially isothermal. Phase change is assumed to occur at heterogeneous nucleation sites and the growth of the vapor bubbles is governed by heat conduction from the liquid to the bubble. The analysis adjusted for friction is applied to liquid nitrogen flow in a venturi and comparison is made with the NASA experimental results of Simoneau. Good agreement with the experiments is obtained when one assumes the effective activation energy for nucleus formation to be small but nonzero. The computed pressure distributions deviate from the experimental results in the throat region of the venturi in a manner consistent with centrifugal effects not accounted for in the one-dimensional theory. The results are shown to depend not only on cavitation number but on additional dimensionless parameters governing the nonequilibrium production and subsequent growth of nuclei.

Rohatgi, U. S.↗

Analytic modeling of soot nucleation under fuel rich conditions

The objective of the present research is to construct a soot nucleation model according to a proposed chemical kinetic scheme to delineate quantitatively the nucleation mechanism in the soot formation process. Instead of following the traditional views which generally associate sooting with the homogeneous nucleation process in phase transformation or polymerization, we choose a chemical kinetic approach. In our proposed scheme the number of carbon atoms in the intermediate species between the fuel molecule and soot nuclei is continuously increased by radical additions. The number of hydrogen atoms in the intermediate species on the other hand is steadily decreased by radical dehydrogenation. When the number of carbon atoms in each of the intermediate molecules has exceeded a certain limit and the number of hydrogen atoms has fallen below a certain level, they may coagulate with one and another to form a larger molecule which is regarded as the initial soot nuclei in the present theory. Further coagulation and surface growth of the nuclei will lead to observable soot particles.

Yang, C. H.↗

Nonlinear decay of magnetic helicity in magnetohydrodynamic turbulence with a mean magnetic field

We show that the magnetic helicity associated with fluctuations in homogeneous incompressible magnetohydrodynamic (MHD) turbulence with a mean magnetic field decays in time because of nonlinear processes. Evidence is obtained numerically, by use of both dissipative and nondissipative spectral method simulations. The described effect stands in contrast to expectations based on studies of MHD turbulence without an applied mean field, in which magnetic helicity is transferred nonlinearly to long wavelengths and is preserved in time because of selective decay when dissipation is present. The process of nonlinear decay is described in terms of a generalized ideal helicity invariant and is characterized by the transient production of a mean induced electric field aligned with the applied magnetic field, an effect reminiscent of the alpha effect in dynamo theory. A simple phenomenological model for the decay process is proposed.

Stribling, Troy↗

Anisotropic Developments for Homogeneous Shear Flows

The general decomposition of the spectral correlation tensor R(sub ij)(k) by Cambon et al. (J. Fluid Mech., 202, 295; J. Fluid Mech., 337, 303) into directional and polarization components is applied to the representation of R(sub ij)(k) by spherically averaged quantities. The decomposition splits the deviatoric part H(sub ij)(k) of the spherical average of R(sub ij)(k) into directional and polarization components H(sub ij)(sup e)(k) and H(sub ij)(sup z)(k). A self-consistent representation of the spectral tensor in the limit of weak anisotropy is constructed in terms of these spherically averaged quantities. The directional polarization components must be treated independently: models that attempt the same representation of the spectral tensor using the spherical average H(sub ij)(k) alone prove to be inconsistent with Navier-Stokes dynamics. In particular, a spectral tensor consistent with a prescribed Reynolds stress is not unique. The degree of anisotropy permitted by this theory is restricted by realizability requirements. Since these requirements will be less severe in a more accurate theory, a preliminary account is given of how to generalize the formalism of spherical averages to higher expansion of the spectral tensor. Directionality is described by a conventional expansion in spherical harmonics, but polarization requires an expansion in tensorial spherical harmonics generated by irreducible representations of the spatial rotation group SO(exp 3). These expansions are considered in more detail in the special case of axial symmetry.

Cambon, Claude↗

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗

Forced oscillations in quadratically damped systems

Bayliss (1975) has studied the question whether in the case of linear differential equations the relationship between the stability of the homogeneous equations and the existence of almost periodic solutions to the inhomogeneous equation is preserved by finite difference approximations. In the current investigation analogous properties are considered for the case in which the damping is quadratic rather than linear. The properties of the considered equation for arbitrary forcing terms are examined and the validity is proved of a theorem concerning the characteristics of the unique solution. By using the Lipschitz continuity of the mapping and the contracting mapping principle, almost periodic solutions can be found for perturbations of the considered equation. Attention is also given to the Lipschitz continuity of the solution operator and the results of numerical tests which have been conducted to test the discussed theory.

Bayliss, A.↗

A study of the application of power-spectral methods of generalized harmonic analysis to gust loads on airplanes

The applicability of some results from the theory of generalized harmonic analysis (or power-spectral analysis) to the analysis of gust loads on airplanes in continuous rough air is examined. The general relations for linear systems between power spectrums of a random input disturbance and an output response are used to relate the spectrum of airplane load in rough air to the spectrum of atmospheric gust velocity. The power spectrum of loads is shown to provide a measure of the load intensity in terms of the standard deviation (root mean square) of the load distribution for an airplane in flight through continuous rough air. For the case of a load output having a normal distribution, which appears from experimental evidence to apply to homogeneous rough air, the standard deviation is shown to describe the probability distribution of loads or the proportion of total time that the load has given values. Thus, for airplane in flight through homogeneous rough air, the probability distribution of loads may be determined from a power-spectral analysis. In order to illustrate the application of power-spectral analysis to gust-load analysis and to obtain an insight into the relations between loads and airplane gust-response characteristics, two selected series of calculations are presented. The results indicate that both methods of analysis yield results that are consistent to a first approximation.

Press, Harry↗

Analysis of the two-point velocity correlations in turbulent boundary layer flows

The general objective of the present work is to explore the use of Rapid Distortion Theory (RDT) in analysis of the two-point statistics of the log-layer. RDT is applicable only to unsteady flows where the non-linear turbulence-turbulence interaction can be neglected in comparison to linear turbulence-mean interactions. Here we propose to use RDT to examine the structure of the large energy-containing scales and their interaction with the mean flow in the log-region. The contents of the work are twofold: First, two-point analysis methods will be used to derive the law-of-the-wall for the special case of zero mean pressure gradient. The basic assumptions needed are one-dimensionality in the mean flow and homogeneity of the fluctuations. It will be shown that a formal solution of the two-point correlation equation can be obtained as a power series in the von Karman constant, known to be on the order of 0.4. In the second part, a detailed analysis of the two-point correlation function in the log-layer will be given. The fundamental set of equations and a functional relation for the two-point correlation function will be derived. An asymptotic expansion procedure will be used in the log-layer to match Kolmogorov's universal range and the one-point correlations to the inviscid outer region valid for large correlation distances.

Oberlack, M.↗

Polarized View of Supercooled Liquid Water Clouds

Supercooled liquid water (SLW) clouds, where liquid droplets exist at temperatures below 0 C present a well known aviation hazard through aircraft icing, in which SLW accretes on the airframe. SLW clouds are common over the Southern Ocean, and climate-induced changes in their occurrence is thought to constitute a strong cloud feedback on global climate. The two recent NASA field campaigns POlarimeter Definition EXperiment (PODEX, based in Palmdale, California, January-February 2013) and Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC4RS, based in Houston, Texas in August- September 2013) provided a unique opportunity to observe SLW clouds from the high-altitude airborne platform of NASA's ER-2 aircraft. We present an analysis of measurements made by the Research Scanning Polarimeter (RSP) during these experiments accompanied by correlative retrievals from other sensors. The RSP measures both polarized and total reflectance in 9 spectral channels with wavelengths ranging from 410 to 2250 nm. It is a scanning sensor taking samples at 0.8deg intervals within 60deg from nadir in both forward and backward directions. This unique angular resolution allows for characterization of liquid water droplet size using the rainbow structure observed in the polarized reflectances in the scattering angle range between 135deg and 165deg. Simple parametric fitting algorithms applied to the polarized reflectance provide retrievals of the droplet effective radius and variance assuming a prescribed size distribution shape (gamma distribution). In addition to this, we use a non-parametric method, Rainbow Fourier Transform (RFT),which allows retrieval of the droplet size distribution without assuming a size distribution shape. We present an overview of the RSP campaign datasets available from the NASA GISS website, as well as two detailed examples of the retrievals. In these case studies we focus on cloud fields with spatial features varying between glaciated and liquid phases at altitudes as high as 10 km, which correspond to temperatures close to the homogeneous freezing temperature of pure water drops (about -35 C or colder). The multimodal droplet size distributions retrieved from RSP data in these cases are consistent with the multi-layer cloud structure observed by correlative Cloud Physics Lidar (CPL) measurements.

Clouds↗

A Novel NASICON‐Type Na 4 MnCr(PO 4 ) 3 Demonstrating the Energy Density Record of Phosphate Cathodes for Sodium‐Ion Batteries

Abstract Sodium‐ion batteries (SIBs) have attracted incremental attention as a promising candidate for grid‐scale energy‐storage applications. To meet practical requirements, searching for new cathode materials with high energy density is of great importance. Herein, a novel Na superionic conductor (NASICON)‐type Na 4 MnCr(PO 4 ) 3 is developed as a high‐energy cathode for SIBs. The Na 4 MnCr(PO 4 ) 3 nanoparticles homogeneously embedded in a carbon matrix can present an extraordinary reversible capacity of 160.5 mA h g −1 with three‐electron reaction at ≈3.53 V during the Na + extraction/insertion process, realizing an unprecedentedly high energy density of 566.5 Wh kg −1 in the phosphate cathodes for SIBs. It is intriguing to reveal the underlying mechanism of the unique Mn 2+ /Mn 3+ , Mn 3+ /Mn 4+ , and Cr 3+ /Cr 4+ redox couples via X‐ray absorption near‐edge structure spectroscopy. The whole electrochemical reaction undergoes highly reversible single‐phase and biphasic transitions with a moderate volume change of 7.7% through in situ X‐ray diffraction and ex situ high‐energy synchrotron X‐ray diffraction. Combining density functional theory (DFT) calculations with the galvanostatic intermittent titration technique, the superior performance is ascribed to the low ionic‐migration energy barrier and desirable Na‐ion diffusion kinetics. The present work can offer a new insight into the design of multielectron‐reaction cathode materials for SIBs.

Zhang, Jian↗

Scale-free, programmable design of morphable chain loops of kilobots and colloidal motors

Micron-scale robots require systems that can morph into arbitrary target configurations controlled by external agents such as heat, light, electricity, and chemical environment. Achieving this behavior using conventional approaches is challenging because the available materials at these scales are not programmable like their macroscopic counterparts. To overcome this challenge, we propose a design strategy to make a robotic machine that is both programmable and compatible with colloidal-scale physics. Our strategy uses motors in the form of active colloidal particles that constantly propel forward. We sequence these motors end-to-end in a closed chain forming a two-dimensional loop that folds under its mechanical constraints. We encode the target loop shape and its motion by regulating six design parameters, each scale-invariant and achievable at the colloidal scale. We demonstrate the plausibility of our design strategy using centimeter-scale robots called kilobots . We use Brownian dynamics simulation to explore the large design space beyond that possible with kilobots, and present an analytical theory to aid the design process. Multiple loops can also be fused together to achieve several complex shapes and robotic behaviors, demonstrated by folding a letter shape “M,” a dynamic gripper , and a dynamic pacman . The material-agnostic, scale-free, and programmable nature of our design enables building a variety of reconfigurable and autonomous robots at both colloidal scales and macroscales.

36 MATERIALS SCIENCE↗

Radiative transfer calculated from a Markov chain formalism

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Esposito, L. W.↗

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗

Photoinduced electron transfer from zinc meso -tetraphenylporphyrin to a one-dimensional perylenediimide aggregate: Probing anion delocalization effects

Organic photovoltaics incorporating non-fullerene acceptors based on perylenediimide (PDI) now rival fullerene acceptor-based devices in performance, although the mechanisms of charge generation in PDI-based devices are not yet fully understood. Fullerene-based systems are proposed to undergo electron transfer directly from the photoexcited donor into a band of delocalized acceptor states, thus increasing charge generation efficiency. Similarly, anion delocalization has been shown to enhance the rate of electron transfer from a photoexcited donor to two electronically coupled PDI acceptors. Here we investigate how additional electron acceptors may further increase the rate of electron transfer from the donor zinc meso-tetraphenylporphyrin (ZnTPP) to an aggregate of PDI acceptors (PDI[Formula: see text]. Femtosecond transient visible and mid-infrared absorption spectroscopies show that the rate of electron transfer from 1 *ZnTPP to the PDI assembly ZnTPP 2 -PDI 3 is statistically identical to that of the previously examined ZnTPP-PDI 2 . A Marcus theory analysis indicates that the parameters governing electron transfer are nearly identical for the two molecules, suggesting that the maximum electron transfer rate enhancement has been achieved in a cofacial PDI dimer because the ZnTPP directly couples to the first two PDI acceptors whereas the coupling to the third PDI is too weak.

bio-inspired↗

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