Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cross flow turbine”

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.

155 records · Page 9

Optical Sensing of Combustion Instabilities in Gas Turbines

In a continuing program of research and development, a system has been demonstrated that makes high-speed measurements of thermal infrared radiance from gas-turbine engine exhaust streams. When a gas-turbine engine is operated under conditions that minimize the emission of pollutants, there is a risk of crossing the boundary from stable to unstable combustion. Combustion instability can lead to engine damage and even catastrophic failure. Sensor systems of the type under development could provide valuable data during the development testing of gas-turbine engines or of engine components. A system of the type under development makes high-speed measurements of thermal infrared radiance from the engine exhaust stream. The sensors of this system can be mounted outside the engine, which eliminates the need for engine case penetrations typical with other engine dynamics monitors. This is an important advantage in that turbine-engine manufacturers consider such penetrations to be very undesirable. A prototype infrared sensor system has been built and demonstrated on a turbine engine. This system includes rugged and inexpensive near-infrared sensors and filters that select wavelengths of infrared radiation for high sensitivity. In experiments, low-frequency signatures were consistently observed in the detector outputs. Under some conditions, the signatures also included frequency components having one or two radiance cycles per engine revolution. Although it has yet to be verified, it is thought that the low-frequency signatures may be associated with bulk-mode combustion instabilities or flow instabilities in the compressor section of the engine, while the engine- revolution-related signatures may be indicative of mechanical problems in the engine. The system also demonstrated the ability to detect transient high-radiance events. These events indicate hot spots in the exhaust stream and were found to increase in frequency during engine acceleration.

Markham, James R.↗

Advanced Turbulence Modeling Concepts

The ZCET program developed at NASA Glenn Research Center is to study hydrogen/air injection concepts for aircraft gas turbine engines that meet conventional gas turbine performance levels and provide low levels of harmful NOx emissions. A CFD study for ZCET program has been successfully carried out. It uses the most recently enhanced National combustion code (NCC) to perform CFD simulations for two configurations of hydrogen fuel injectors (GRC- and Sandia-injector). The results can be used to assist experimental studies to provide quick mixing, low emission and high performance fuel injector designs. The work started with the configuration of the single-hole injector. The computational models were taken from the experimental designs. For example, the GRC single-hole injector consists of one air tube (0.78 inches long and 0.265 inches in diameter) and two hydrogen tubes (0.3 inches long and 0.0226 inches in diameter opposed at 180 degree). The hydrogen tubes are located 0.3 inches upstream from the exit of the air element (the inlet location for the combustor). To do the simulation, the single-hole injector is connected to a combustor model (8.16 inches long and 0.5 inches in diameter). The inlet conditions for air and hydrogen elements are defined according to actual experimental designs. Two crossing jets of hydrogen/air are simulated in detail in the injector. The cold flow, reacting flow, flame temperature, combustor pressure and possible flashback phenomena are studied. Two grid resolutions of the numerical model have been adopted. The first computational grid contains 0.52 million elements, the second one contains over 1.3 million elements. The CFD results have shown only about 5% difference between the two grid resolutions. Therefore, the CFD result obtained from the model of 1.3-million grid resolution can be considered as a grid independent numerical solution. Turbulence models built in NCC are consolidated and well tested. They can handle both coarse and fine grids near the wall. They can model the effect of anisotropy of turbulent stresses and the effect of swirling. The chemical reactions of Magnusson model and ILDM method were both used in this study.

Shih, Tsan-Hsing↗

Dual Active Surface Heat Flux Gage Probe

A unique plug-type heat flux gage probe was tested in the NASA Ames Research Center 2x9 turbulent flow duct facility. The probe was fabricated by welding a miniature dual active surface heat flux gage body to the end of a hollow metal cylindrical bolt containing a metal inner tube. Cooling air flows through the inner tube, impinges onto the back of the gage body and then flows out through the annulus formed between the inner tube and the hollow bolt wall. Heat flux was generated in the duct facility with a Huels arc heater. The duct had a rectangular cross section and one wall was fabricated from 2.54 centimeter thick thermal insulation rigid surface material mounted onto an aluminum plate. To measure heat flux, the probe was inserted through the plate and insulating materials with the from of the gage located flush with the hot gas-side insulation surface. Absorbed heat fluxes measured with the probe were compared with absorbed heat fluxes measured with six water-cooled reference calorimeters. These calorimeters were located in a water-cooled metal duct wall which was located across from the probe position. Correspondence of transient and steady heat fluxes measured with the reference calorimeters and heat flux gage probe was generally within a satisfactory plus or minus 10 percent. This good correspondence was achieved even though the much cooler probe caused a large surface temperature disruption of 1000K between the metal gage and the insulation. However, this temperature disruption did not seriously effect the accuracy of the heat flux measurement. A current application for dual active surface heat flux gages is for transient and steady absorbed heat flux, surface temperature and heat transfer coefficient measurements on the surface of an oxidizer turbine inlet deflector operating in a space shuttle test bed engine.

Liebert, Curt H.↗

Navier-Stokes analysis and experimental data comparison of compressible flow within ducts

Many aircraft employ ducts with centerline curvature or changing cross-sectional shape to join the engine with inlet and exhaust components. S-ducts convey air to the engine compressor from the intake and often decelerate the flow to achieve an acceptable Mach number at the engine compressor by increasing the cross-sectional area downstream. Circular-to-rectangular transition ducts are used on aircraft with rectangular exhaust nozzles to connect the engine and nozzle. To achieve maximum engine performance, the ducts should minimize flow total pressure loss and total pressure distortion at the duct exit. Changes in the curvature of the duct centerline or the duct cross-sectional shape give rise to streamline curvature which causes cross stream pressure gradients. Secondary flows can be caused by deflection of the transverse vorticity component of the boundary layer. This vortex tilting results in counter-rotating vortices. Additionally, the adverse streamwise pressure gradient caused by increasing cross-sectional area can lead to flow separation. Vortex pairs have been observed in the exit planes of both duct types. These vortices are due to secondary flows induced by pressure gradients resulting from streamline curvature. Regions of low total pressure are produced when the vortices convect boundary layer fluid into the main flow. The purpose of the present study is to predict the measured flow field in a diffusing S-duct and a circular-to-rectangular transition duct with a full Navier-Stokes computer program, PARC3D, and to compare the numerical predictions with new detailed experimental measurements. The work was undertaken to extend previous studies and to provide additional CFD validation data needed to help model flows with strong secondary flow and boundary layer separation. The S-duct computation extends the study of Smith et al, and Harloff et al, which concluded that the computation might be improved by using a finer grid and more advanced turbulence models. The present study compares results for both the Baldwin-Lomas and k-epsilon turbulence models and is conducted with a refined grid. For the transition duct, two inlet conditions were considered, the first with straight flow and the second with swirling flow. The first case permits examination of the effects of the geometric transition on the flow field, while the second case includes the rotational flow effect characteristic of a gas turbine engine.

Harloff, G. J.↗

Glenn Heat Transfer Simulation and Solver Graphical User Interface: Development and Testing

In the Tui ine Branch of the Turbomachinery and Propulsion Systems Division, researching and developing efficient turbine aerothermodynamics technologies is the main objective. Creating effective turbines for jet engines is a process which, if based purely on physical experimental testing, would be extremely expensive. It is for this reason, and also for the reasons of speed and ease, that the Turbine Branch spends a large amount of effort working with simulations of turbines. Specifically, they focus their work on two main fields: Computational Field Dynamics (CFD), and Experimental data analysis. The experimental field involves comparing experimental results to simulated results, whereas the CFD field involves running these simulations. The simulations are applied to aerodynamics and heat transfer cases, for both steady and unsteady flow conditions. By and large this work is applied to the domain of flow and heat transfer in axial turbines. The main application used to run these heat flow simulations is GlennHT. This program, recently rewritten in FORTRAN 90, allows the user to input a job file which specifies all the necessary parameters needed to simulate flow through a user-defined grid. There are several other executables used as well, ranging in application from converting grid files to and from particular formats, to merging blocks in a connectivity file, to converting connectivity files to a GlennHT compatible format. All of these executables are run from the command line in a terminal; some of them have interactive prompts where the user must specify the files to be manipulated after the program starts, while others take all of their parameters from the command line. With this amount of variation comes a good deal of commands and formats to memorize, which can cause slower and less efficient work, as users may forget how to execute a certain program, or not remember the pathnames of the files they wish to use. Two years ago, steps were made to expedite this process with a graphical user interface (GUI) that combines the functionality of all the executables along with adding some new functionality, such as residuals graphing and boundary conditions creation. Upon my beginning here at Glenn, many parts of the GUI, which was developed in Java, were nonfunctional. There were also issues with cross-platforming, as systems in the branch were transitioning from Silicon Graphics (SGI) machines to Linux machines. My goals this summer are to finish the parts of the GUI that are not yet completed, fix parts that did not work correctly, expand the functionality to include other useful features, such as grid surface highlighting, and make the system compatible with both Linux and SGI. I will also be heavily testing the system and providing sufficient documentation on how to use the GUI, as no such documentation existed previously.

Kardamis, Joseph R.↗

Design and Modeling of an Open-Source Baseline Floating Marine Turbine: Preprint

Marine energy resources such as river, tidal, and ocean currents can provide abundant, predictable, and clean power to both densely populated and remote communities globally. However, marine turbine technologies are still in an early stage of development. Recently, features have been added to the National Renewable Energy Laboratory's open-source wind turbine modeling tool OpenFAST to support the simulation of axial-flow marine turbines. These new features enable marine energy developers to predict the performance and safety of their turbines and allow researchers and students to advance the field through the exploration of novel designs. Developing and demonstrating the new capabilities of OpenFAST for marine turbines requires a baseline design that can be used as an example and test case. There are limited open-source turbine models suitable for this purpose. One existing baseline design is the Reference Model 1 (RM1), a fixed-bottom, dual-rotor marine turbine that was published in 2014 as part of the Department of Energy's Water Power Reference Model Project. To test and demonstrate OpenFAST's ability to model floating marine turbines, several modifications were made to the RM1. A floating platform and mooring system were designed to support a single RM1 rotor in a wide range of possible current and wave conditions. The system and model are intended to be a stable reference case and starting point for research efforts. The platform features all cylindrical members and familiar design elements from the offshore wind industry; the mooring system uses simple catenary chains. While not optimized for the lowest platform cost, the design is understandable and robust and can serve as a baseline for floating marine turbine technologies. This is a useful starting point for studies to explore new technologies or design methodologies. Additionally, the floating RM1 design enables cross-collaboration among industry, academia, and national lab researchers and serves as an educational platform for newcomers to marine energy. This short paper presents the specifications of the turbine, including its floating support structure and mooring system. The rationale behind the selection of certain design parameters is described. To enable adoption of this baseline case for research and educational efforts, an OpenFAST model of the full system has been made publicly available at https://github.com/OpenFAST/r-test/tree/main/glue-codes/openfast/MHK_RM1_Floating.

floating↗

Abradable dual-density ceramic turbine seal system

A plasma sprayed dual density ceramic abradable seal system for direct application to the HPT seal shroud of small gas turbine engines. The system concept is based on the thermal barrier coating and depends upon an additional layer of modified density ceramic material adjacent to the gas flow path to provide the desired abradability. This is achieved by codeposition of inert fillers with yttria stabilized zirconia (YSZ) to interrupt the continuity of the zirconia struture. The investigation of a variety of candidate fillers, with hardness values as low as 2 on Moh's scale, led to the conclusion that solid filler materials in combination with a YSZ matrix, regardless of their hardness values, have a propensity for compacting rather than shearing as originally expected. The observed compaction is accompanied by high energy dissipation in the rub interaction, usually resulting in the adhesive transfer of blade material to the stationary seal member. Two YSZ based coating systems which incorported hollow alumino silicate spheres as density reducing agents were surveyed over the entire range of compositions from 100 percent filler to 100 percent YSZ. Abradability and erosion characteristics were determined, hardness and permeability characterized, and engine experience acquired with several system configurations.

Clingman, D. L.↗

Effects of Compression, Staging, and Braid Angle on Braided Rope Seal Performance

Future turbine engines and industrial systems will be operating at increased temperatures to achieve more demanding efficiency and performance goals. In the highest temperature sections of the engine new material systems such as ceramics and intermetallics are being considered to withstand the harsh thermal environment. Components constructed of these low expansion-rate materials experience thermal strains and a resulting reduction of life when rigidly attached to high expansion-rate, superalloy support structures. Seals are being designed to both seal and to serve as compliant mounts allowing for relative thermal growths between high temperature but brittle primary structures and the surrounding support structures. Previous seal research yielded several braided rope seal designs which demonstrated the ability to both seal and serve as a compliant mount. The hybrid seal was constructed of an all-ceramic (alumina-silica) core overbraided with a superalloy wire sheath (cobalt based superalloy). The all ceramic seal was constructed of an all-ceramic (alumina-silica) core overbraided with multiple ceramic (alumina-silica) sheath layers. Program goals for braided rope seals are to improve flow resistance and/or seal resilience. To that end, the current report studies the test results of: baseline and modified hybrid seals; two stage hybrid and two stage all-ceramic seal configurations; and single stage hybrid and single stage all-ceramic seal configurations for a range of seal crush conditions. Hybrid seal modifications include increasing the sheath braid angle and core coverage. For the same percent seal cross-sectional crush, results show that increasing the hybrid seal braid angle increased seal stiffness and seal unit load, resulting in flows approximately one third of the baseline hybrid seal flows. For both hybrid and all-ceramic seals, two stage seal configurations significantly outperformed single stage configurations. Two stage seal flows were at least 30% less than the single stage seal flows for the same seal crush. Furthermore, test results of single stage seals indicate that for both all-ceramic and hybrid seals, a specific seal crush condition exists at which minimum flows are achieved (i.e. increasing seal crush beyond a certain point does not result in better flow performance). Flow results are presented for a range of pressures and temperatures from ambient to 1300 F, before and after scrubbing. Compression tests results show that for both all-ceramic and hybrid seals, seal preload and stiffness increase with seal crush, but residual seal interference remains constant.

Steinetz, Bruce M.↗

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗

Experimental Studies of Low-Pressure Turbine Flows and Flow Control. Streamwise Pressure Profiles and Velocity Profiles

This report summarizes research performed in support of the NASA Glenn Research Center (GRC) Low-Pressure Turbine (LPT) Flow Physics Program. The work was performed experimentally at the U.S. Naval Academy faculties. The geometry corresponded to "Pak B" LPT airfoil. The test section simulated LPT flow in a passage. Three experimental studies were performed: (a) Boundary layer measurements for ten baseline cases under high and low freestream turbulence conditions at five Reynolds numbers of 25,000, 50,000, 100,000, 200,000, and 300,000, based on passage exit velocity and suction surface wetted length; (b) Passive flow control studies with three thicknesses of two-dimensional bars, and two heights of three-dimensional circular cylinders with different spanwise separations, at same flow conditions as the 10 baseline cases; (c) Active flow control with oscillating synthetic (zero net mass flow) vortex generator jets, for one case with low freestream turbulence and a low Reynolds number of 25,000. The Passive flow control was successful at controlling the separation problem at low Reynolds numbers, with varying degrees of success from case to case and varying levels of impact at higher Reynolds numbers. The active flow control successfully eliminated the large separation problem for the low Reynolds number case. Very detailed data was acquired using hot-wire anemometry, including single and two velocity components, integral boundary layer quantities, turbulence statistics and spectra, turbulent shear stresses and their spectra, and intermittency, documenting transition, separation and reattachment. Models were constructed to correlate the results. The report includes a summary of the work performed and reprints of the publications describing the various studies.This report summarizes research performed in support of the NASA Glenn Research Center (GRC) Low-Pressure Turbine (LPT) Flow Physics Program. The work was performed experimentally at the U.S. Naval Academy faculties. The geometry corresponded to "Pak B" LPT airfoil. The test section simulated LPT flow in a passage. Three experimental studies were performed: (a) Boundary layer measurements for ten baseline cases under high and low freestream turbulence conditions at five Reynolds numbers of 25,000, 50,000, 100,000, 200,000, and 300,000, based on passage exit velocity and suction surface wetted length; (b) Passive flow control studies with three thicknesses of two-dimensional bars, and two heights of three-dimensional circular cylinders with different spanwise separations, at same flow conditions as the 10 baseline cases; (c) Active flow control with oscillating synthetic (zero net mass flow) vortex generator jets, for one case with low freestream turbulence and a low Reynolds number of 25,000. The Passive flow control was successful at controlling the separation problem at low Reynolds numbers, with varying degrees of success from case to case and varying levels of impact at higher Reynolds numbers. The active flow control successfully eliminated the large separation problem for the low Reynolds number case. Very detailed data was acquired using hot-wire anemometry, including single and two velocity components, integral boundary layer quantities, turbulence statistics and spectra, turbulent shear stresses and their spectra, and intermittency, documenting transition, separation and reattachment. Models were constructed to correlate the results. The report includes a summary of the work performed and reprints of the publications describing the various studies. The folders in this supplement contain processed data in ASCII format. Streamwise pressure profiles and velocity profiles are included. The velocity profiles were acquired using single sensor and cross sensor hot-wire probes which were traversed from the wall to the freestream at various streamwise locations. In some of the flow control cases (3D Trips and Jets) profiles were acquired at multiple spanwise locations.

Volino, Ralph↗

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗