Search NASA⌕ Search

SEARCH · Search NASA

Results for “normalizing flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

Analysis and Modeling of a Two-Phase Jet Pump of a Flow Boiling Test Facility for Aerospace Applications

Jet pumps are devices capable of pumping fluids to a higher pressure employing a nozzle/diffuser/mixing chamber combination. A primary fluid is usually allowed to pass through a converging-diverging nozzle where it can accelerate to supersonic speeds at the nozzle exit. The relatively high kinetic energy that the primary fluid possesses at the nozzle exit is accompanied by a low pressure region in order to satisfy Bernoulli's equation. The low pressure region downstream of the nozzle exit permits a secondary fluid to be entrained into and mixed with the primary fluid in a mixing chamber located downstream of the nozzle. Several combinations may exist in terms of the nature of the primary and secondary fluids in so far as whether they are single or two-phase fluids. Depending on this, the jet pump may be classified as gas/gas, gas/liquid, liquid/liquid, two-phase/liquid, or similar combinations. The mixing chamber serves to create a homogeneous single-phase or two-phase mixture which enters a diffuser where the high kinetic energy of the fluid is converted into pressure energy. If the fluid mixture entering the diffuser is in the supersonic flow regime, a normal shock wave usually develops inside the diffuser. If the fluid mixture is one that can easily change phase, a condensation shock would normally develop. Because of the overall rise in pressure in the diffuser as well as the additional rise in pressure across the shock layer, condensation becomes more likely. Associated with the pressure rise across the shock is a velocity reduction from the supersonic to the subsonic range. If the two-phase flow entering the diffuser is predominantly gaseous with liquid droplets suspended in it, it will transform into a predominantly liquid flow containing gaseous bubbles (bubbly flow) somewhere in the diffuser. While past researchers have been able to model the two-phase flow jet pump using the one-dimensional assumption with no shock waves and no phase change, there is no research known to the authors apart from that of Anand (1992) which accounted for condensation shocks. One of the objectives of this research effort is to develop a comprehensive model in which the effects of phase slip and inter-phase heat transfer as well as the wall friction and shock waves are accounted for. While this modeling effort is predominantly analytical in nature and is primarily intended to provide a parametric understanding of the jet pump performance under different operating scenarios, another parallel effort employing a commercial CFD code is also implemented. The latter effort is primarily intended to model an axisymmetric counterpart of the problem in question. The viability of using the CFD code to model a two-phase flow jet pump will be assessed by attempting to recreate some of the existing performance data of similar jet pumps. The code will eventually be used to generate the jet pump performance characteristics of several scenarios involving jet pump geometries as well as flow regimes in order to be able to determine an optimum design which would be suitable for a two-phase flow boiling test facility at NASA-Marshall. Because of the extensive nature of the analytical model developed, the following section will only provide very brief highlights of it, while leaving the details to a more complete report submitted to the NASA colleague. This report will also contain some of the simulation results obtained using the CFD code.

Sherif, S. A.↗

Particle acceleration in step function shear flows - A microscopic analysis

The transport of energetic particles in a moving, scattering fluid, which has a large shear in its velocity over a distance small compared with the scattering mean free path is discussed. The analysis is complementary to an earlier paper by Earl, Jokipii, and Morfill (1988), which considered effects of more-gradual shear in the diffusion approximation. The case in which the scattering fluid undergoes a step function change in velocity, in the direction normal to the flow is considered. An analytical, approximate calculation and a Monte Carlo analysis of particle motion are presented. It is found that particles gain energy at a rate proportional to the square of the magnitude of the velocity change.

Jokipii, J. R.↗

Surface Tension Driven Convection Experiment-2 (STDCE-2)

Thermocapillary flows are known to become oscillatory (time-periodic), but how and when they become oscillatory in containers of unit-order aspect ratio are not yet fully understood. The present work is a part of our continuous effort to obtain a better understanding of the phenomenon. Thermocapillary flow experiments in normal gravity are limited to a narrow parametric range in order to minimize gravity and buoyancy effects, which is an important reason for our lack of full understanding of the oscillation phenomenon. One important unanswered question is what role, if any, free surface deformation plays in the oscillation mechanism. For that reason we performed thermocapillary flow experiments, called the Surface Tension Driven Convection Experiment-2 (STDCE-2), aboard the USML-2 Spacelab in 1995. The main objectives of the experiments were to investigate oscillatory thermocapillary flows in microgravity and to clarify the importance of free surface deformation in such flows. Steady and oscillatory thermocapillary flows were generated in cylindrical containers by employing two heating modes. A CO2 laser with adjustable power and beam diameter was used in the Constant Flux (CF) configuration to heat the free surface. The other configuration investigated in STDCE-2 was the Constant Temperature (CT) configuration in which a submerged cylindrical cartridge heater placed at the symmetry (axial) axis of the test container heated the fluid. Both heating modes cause non-uniform temperature distributions on the free surface, which generates thermocapillary flow. The flow field was investigated by flow visualization, and the temperature field was measured by thermistors and an infrared imager. The free surface shape and motion were measured by a Ronchi system. The hardware performed well and we were able to conduct more tests than originally planned. From the successful experiments a large amount of data was acquired. The analysis of the data is now nearly complete. Some important results are presented and discussed herein.

Masud, J.↗

Flow Visualization by Elastic Light Scattering in the Boundary Layer of a Supersonic Flow

We demonstrate instantaneous flow visualization of the boundary layer region of a Mach 2.5 supersonic flow over a flat plate that is interacting with an impinging shock wave. Tests were performed in the Unitary Plan Wind Tunnel (UPWT) at NASA Langley Research Center. The technique is elastic light scattering using 10-nsec laser pulses at 532 nm. We emphasize that no seed material of any kind, including water (H2O), is purposely added to the flow. The scattered light comes from a residual impurity that normally exists in the flow medium after the air drying process. Thus, the technique described here differs from the traditional vapor-screen method, which is typically accomplished by the addition of extra H2O vapor to the airflow. The flow is visualized with a series of thin two-dimensional light sheets (oriented perpendicular to the streamwise direction) that are located at several positions downstream of the leading edge of the model. This geometry allows the direct observation of the unsteady flow structure in the spanwise dimension of the model and also allows the indirect observation of the boundary layer growth in the streamwise dimension.

Herring, G. C.↗

Mechanical and Thermal Forcing for Upslope Flows and Cumulus Convection over the Sierras de Córdoba

Abstract The upslope flow processes affecting the vertical extent of orographic cumulus convection are examined using observations from the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Specifically, clear air returns from the U.S. Department of Energy (DOE) second-generation C-band scanning Atmospheric Radiation Measurement (ARM) precipitation radar (CSAPR2) are used to characterize the structure and variability of the ridge-normal (i.e., up/downslope) flow components, which transport mass to the crest of Argentina’s Sierras de Córdoba and contribute to convective initiation. Data are compiled for the entire CACTI period (October–April), including days with clear skies, shallow cumuli, cumulus congestus, and deep convection. To examine shared variability among >70 000 radar scans, we use (i) a principal component analysis (PCA) to isolate modes of variability in the upslope flow and (ii) composite analysis based on convective outcomes, determined from GOES-16 satellite observations. These data are contextualized with observed surface sensible heat fluxes, thermodynamic profiles, and synoptic-scale analysis. Results indicate distinct thermally and mechanically forced upslope flow modes, modulated by diurnal heating and synoptic-scale variations, respectively. In some instances, there is a superposition of thermal and mechanical forcing, yielding either deeper or shallower upslope flow. The composite analyses based on satellite data show that successively deeper convective outcomes are associated with successively deeper upslope flow layers that more readily transport mass to the ridge crest in conjunction with lower lifting condensation levels, facilitating convective initiation. These results help to isolate the forcing mechanisms for orographic convection and thus provide a foundation for parameterizing orographic convective processes in coarse resolution models.

Meteorology & Atmospheric Sciences↗

Numerical study of twin-jet impingement upwash flow

Two horizontally spaced jets impinging normally on a flat surface create a fountain upwash flow due to the collision of the radially flowing wall jets. This fountain flow is of importance to the dynamics and propulsion of STOVL aircraft. The fountain flow influences the lift forces on the aircraft and the ingestion of hot gases and debris by the engine inlet. In this paper, a multigrid based finite-difference numerical procedure has been applied to solve the equations governing this three-dimensional flow. The standard k-epsilon turbulence model has been used. Comparisons with experimental data reveal that while the mean velocities are predicted with reasonable accuracy, the turbulent kinetic energies are seriously in error. The reasons for this discrepancy could be the intense unsteadiness and large-scale structures of the flow in the near-wall region, which cannot be captured well by any Reynolds-averaged turbulence model.

Pegues, W. J.↗

Design, Construction and Testing of Annular Diffusers for High Speed Civil Transportation Combustor Applications

A theoretical and an experimental design study of subsonic flow through curved-wall annular diffusers has been initiated under this award in order to establish the most pertinent design parameters and hence performance characteristics for such devices, an the implications of their application in the design of engine components in the aerospace industries. The diffusers under this study are expected to contain appreciable regions of stall and the effects of swirl on their performance are being studied. The experimental work involves the application of Computer Aided Design software tool to the development of a suitable annular diffuse geometry and the subsequent downloading of such data to a CNC machine at Central State University (CSU). Two experimental run segments have been completed so far during FY-95 involving flow visualization and diffuser performance evaluation based on Kinetic Energy Dissipation. The method of calculation of the performance of diffusers based on pressure recovery coefficient has been shown to have some shortcomings and so the kinetic energy dissipation approach has been introduced in the run segment two with some success. The application of the discretized, full Navier Stokes and Continuity equations to the numerical study of the problem described above for the time-mean flow is expected to follow. Various models of turbulence are being evaluated for adoption throughout the study and comparisons would be made with experimental data where they exist. Assessment of diffuser performance based on the dissipated mechanical energy would also be made. The result of the investigations are expected to indicate that more cost effective component design of such devices as diffusers which normally contain complex flows can still be achieved.

Okhio, Cyril B.↗

Initial Developments of a Low Drag Variable Depth Acoustic Liner

The purpose of this study is to assess the acoustic performance and relative drag of a new type of variable depth liner containing pairs of resonators coupled together by shared inlet volumes just below the facesheet. This type of liner has the potential to achieve a targeted impedance with fewer openings in the facesheet, and therefore less drag, than previous designs. To better understand the limitations of the concept, three sets of samples were designed and tested in the NASA Langley normal incidence and grazing flow impedance tubes. Each set of samples consisted of a baseline variable depth liner with straight chambers, and a prototype liner with shared inlet volumes. Measurements conducted in both the normal incidence and grazing flow impedance tubes confirm that the prototype liner can achieve the same impedance as the baseline variable depth liner at discrete frequencies with 50% less open area. This results in a 75% reduction in the drag penalty relative to the baseline design. The proposed concept is not, however, able to achieve a flat, broadband impedance as effectively as the baseline liner. So, while the proposed concept is well suited for multitonal design metrics, tradeoffs must be made between liner drag and acoustic performance when broadband attenuation is required.

Schiller, Noah H.↗

Aerodynamic Characteristics of a Circular Cylinder at Mach Number of 6.86 and Angles of Attack up to 90 Degrees

Pressure-distribution and force tests of a circular cylinder have been made in the Langley 11-inch hypersonic tunnel at a Mach number of 6.86, a Reynolds number of 129,000 based on diameter, and angles of attack up to 90 degrees. The results are compared with the hypersonic approximation of Grimminger, Williams, and Young and with a simple modification of the Newtonian flow theory. The comparison of experimental results shows that either theory gives adequate general aerodynamic characteristics but that the modified Newtonian theory gives a more accurate prediction of the pressure distribution. The calculated crossflow drag coefficients plotted as a function of crossflow Mach number were found to be in reasonable agreement with similar results obtained from other investigations at lower supersonic Mach numbers. Comparison of the results of this investigation with data obtained at a lower Mach number indicates that the drag coefficient of a cylinder normal to the flow is relatively constant for Mach numbers above about 4.

Penland, Jim A↗

Low Stretch PMMA Burning in Microgravity: Status of the Ground-Based Program and New ISS Glovebox Experiment SALSA

The objective of this ground-based program is to study low stretch diffusion flames burning PMMA as the solid fuel to determine the relationship between buoyant low stretch burning in normal gravity and forced flow low stretch burning in microgravity. The low stretch is generated in normal gravity by using the buoyant convection induced by burning the bottom of a large radius of curvature sample. Low stretch is also generated using the Combustion Tunnel drop tower rig (2.2 and 5.2 second facilities), which provides a forced convective low velocity flow past smaller radius of curvature samples. Lastly, an ISS glovebox investigation is being developed to study low stretch burning of PMMA spheres to obtain long duration testing needed to accurately assess the flammability and burning characteristics of the material in microgravity. A comparison of microgravity experiment results with normal gravity test results allows us to establish a direct link between a material's burning characteristics in normal gravity (easily measured) with its burning characteristics in extraterrestrial environments, including microgravity forced convective environments. Theoretical predictions and recent experimental results indicate that it should be possible to understand a material's burning characteristics in the low stretch environment of spacecraft (non-buoyant air movement induced by fans and crew disturbances) by understanding its burning characteristics in an equivalent Earth-based low stretch environment (induced by normal gravity buoyancy). Similarly, Earth-based stretch environments can be made equivalent to those in Lunar- and Martian-surface stretch environments (which would induce partial-gravity buoyancy).

Olson, S. L.↗

Pulmonary blood flow distribution in sheep: effects of anesthesia, mechanical ventilation, and change in posture

BACKGROUND: Recent studies providing high-resolution images of pulmonary perfusion have questioned the classical zone model of pulmonary perfusion. Hence the present work was undertaken to provide detailed maps of regional pulmonary perfusion to examine the influence of anesthesia, mechanical ventilation, and posture. METHODS: Pulmonary perfusion was analyzed with intravenous fluorescent microspheres (15 microm) in six sheep studied in four conditions: prone and awake, prone with pentobarbital-anesthesia and breathing spontaneously, prone with anesthesia and mechanical ventilation, and supine with anesthesia and mechanical ventilation. Lungs were air dried at total lung capacity and sectioned into approximately 1,100 pieces (about 2 cm3) per animal. The pieces were weighed and assigned spatial coordinates. Fluorescence was read on a spectrophotometer, and signals were corrected for piece weight and normalized to mean flow. Pulmonary blood flow heterogeneity was assessed using the coefficient of variation of flow data. RESULTS: Pentobarbital anesthesia and mechanical ventilation did not influence perfusion heterogeneity, but heterogeneity increased when the animals were in the supine posture (P < 0.01). Gravitational flow gradients were absent in the prone position but present in the supine (P < 0.001 compared with zero). Pulmonary perfusion was distributed with a hilar-to-peripheral gradient in animals breathing spontaneously (P < 0.05). CONCLUSIONS: The influence of pentobarbital anesthesia and mechanical ventilation on pulmonary perfusion heterogeneity is small compared with the effect of changes in posture. Analysis of flow gradients indicate that gravity plays a small role in determining pulmonary blood flow distribution.

NASA Discipline Cardiopulmonary↗

Hot-jet simulation in cryogenic wind tunnels

In order to evaluate hot jet simulation capability in cryogenic wind tunnel testing, simple theoretical calculations were performed. The similarity parameters, isentropic flow properties, and normal shock relations were calculated for a variety of jet simulation techniques. The results were compared with those estimated for a full scale flight condition. It was shown that the cryogenic wind tunnel testing provides an opportunity for the most accurate hot jet simulation technique. By using a compressed nitrogen gas at ambient or moderately elevated temperatures as a jet gas, most all of the relevant similarity parameters including the jet temperature and velocity ratios and the Reynolds numbers, can be set to the full scale flight values. The only exception is the ratio of specific heats for jet flow. In an attempt to match the ratio of specific heats for the turbojet flow, gases other than pure nitrogen were considered. It was found that a nitrogen/methane mixture at moderately elevated temperature behaves like the real combustion gas. Using this mixture as a jet gas, complete simulation of the full scale turbojet exhaust becomes possible in cryogenic wind tunnels.

Asai, Keisuke↗

Design, construction and testing of annular diffusers for high speed civil transportation combustor applications

A theoretical and an experimental design study of subsonic flow through curved-wall annular diffusers is being carried out in order to establish the most pertinent design parameters for such devices and the implications of their application in the design of engine components in the aerospace industries. This investigation consists of solving numerically the full Navier Stokes and Continuity equations for the time-mean flow. Various models of turbulence are being evaluated for adoption throughout the study and comparisons would be made with experimental data where they exist. Assessment of diffuser performance based on the dissipated mechanical energy would also be made. The experimental work involves the application of Computer Aided Design software tool to the development of a suitable annular diffuser geometry and the subsequent downloading of such data to a CNC machine at Central State University. The results of the investigations are expected to indicate that more cost effective component design of such devices as effective component design of such devices as diffusers which normally contain complex flows can still be achieved. In this regard a review paper was accepted and presented at the First International Conference on High Speed Civil Transportation Research held at North Carolina A&T in December of 1994.

Okhio, Cyril B.↗

Modeling scramjet combustor flowfields with a grid adaptation scheme

The accurate description of flow features associated with the normal injection of fuel into supersonic primary flows is essential in the design of efficient engines for hypervelocity aerospace vehicles. The flow features in such injections are complex with multiple interactions between shocks and between shocks boundary layers. Numerical studies of perpendicular sonic N2 injection and mixing in a Mach 3.8 scramjet combustor environment are discussed. A dynamic grid adaptation procedure based on the equilibration of spring-mass system is employed to enhanced the description of the complicated flow features. Numerical results are compared with experimental measurements and indicate that the adaptation procedure enhances the capability of the modeling procedure to describe the flow features associated with scramjet combustor components.

Ramakrishnan, R.↗

Effect of Slow External Flow on Flame Spreading over Solid Material: Opposed Spreading over Polyethylene Wire Insulation

The effect of slow external flow on solid combustion is very important from the view of fire safety in space because the solid material in spacecraft is generally exposed to the low air flow for ventilation. Further, the effect of low external flow on fuel combustion is generally fundamental information for industrial combustion system, such as gas turbine, boiler incinerator and so on. However, it is difficult to study the effect of low external flow on solid combustion in normal gravity, because the buoyancy-induced flow strongly disturbs the flow field, especially for low flow velocity. In this research therefore, the effect of slow external flow on opposed flame spreading over polyethylene (PE) wire insulation have been investigated in microgravity. The microgravity environment was provided by Japan Microgravity Center (JAMIC) in Japan and KC-135 at NASA GRC. The tested flow velocity range is 0-30cm/s with different oxygen concentration and inert gas component.

Fujita, O.↗

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

The numerical simulation of many aerodynamic non-periodic flows of practical interest involves discretized computational domains that often must be artificially truncated. Appropriate boundary conditions are required at these truncated domain boundaries, and ideally, these boundary conditions should be perfectly "absorbing" or "nonreflecting" so that they do not contaminate the flow field in the interior of the domain. The proper specification of these boundaries is critical to the stability, accuracy, convergence, and quality of the numerical solution, and has been the topic of considerable research. The need for accurate boundary specification has been underscored in recent years with efforts to apply higher-fidelity methods (DNS, LES) in conjunction with high-order low-dissipation numerical schemes to realistic flow configurations. One of the most popular choices for specifying these boundaries is the characteristics-based boundary condition where the linearized flow field at the boundaries are decomposed into characteristic waves using either one-dimensional Riemann or other multi-dimensional Riemann approximations. The values of incoming characteristics are then suitably modified. The incoming characteristics are specified at the in flow boundaries, and at the out flow boundaries the variation of the incoming characteristic is zeroed out to ensure no reflection. This, however, makes the problem ill-posed requiring the use of an ad-hoc parameter to allow small reflections that make the solution stable. Generally speaking, such boundary conditions work reasonably well when the characteristic flow direction is normal to the boundary, but reflects spurious energy otherwise. An alternative to the characteristic-based boundary condition is to add additional "buffer" regions to the main computational domain near the artificial boundaries, and solve a different set of equations in the buffer region in order to minimize acoustic reflections. One approach that has been used involves modeling the pressure fluctuations as acoustic waves propagating in the far-field relative to a single noise-source inside the buffer region. This approach treats vorticity-induced pressure fluctuations the same as acoustic waves. Another popular approach, often referred to as the "sponge layer," attempts to dampen the flow perturbations by introducing artificial dissipation in the buffer region. Although the artificial dissipation removes all perturbations inside the sponge layer, incoming waves are still reflected from the interface boundary between the computational domain and the sponge layer. The effect of these refkections can be somewhat mitigated by appropriately selecting the artificial dissipation strength and the extent of the sponge layer. One of the most promising variants on the buffer region approach is the Perfectly Matched Layer (PML) technique. The PML technique mitigates spurious reflections from boundaries and interfaces by dampening the perturbation modes inside the buffer region such that their eigenfunctions remain unchanged. The technique was first developed by Berenger for application to problems involving electromagnetic wave propagation. It was later extended to the linearized Euler, Euler and Navier-Stokes equations by Hu and his coauthors. The PML technique ensures the no-reflection property for all waves, irrespective of incidence angle, wavelength, and propagation direction. Although the technique requires the solution of a set of auxiliary equations, the computational overhead is easily justified since it allows smaller domain sizes and can provide better accuracy, stability, and convergence of the numerical solution. In this paper, the PML technique is developed in the context of a high-order spectral-element Discontinuous Galerkin (DG) method. The technique is compared to other approaches to treating the in flow and out flow boundary, such as those based on using characteristic boundary conditions and sponge layers. The superiority of the current PML technique over other approaches is demonstrated for a range of test cases, viz., acoustic pulse propagation, convective vortex, shear layer flow, and low-pressure turbine cascade flow. The paper is structured as follows. We first derive the PML equations from the non{linear Euler equations. A short description of the higher-order DG method used is then described. Preliminary results for the four test cases considered are then presented and discussed. Details regarding current work that will be included in the final paper are also provided.

Spectral↗

Large-Vortex Capture by a Wing at Very High Angles of Attack

In generating the lift on a wing, the static stall is a severe barrier. As the angle of attack, alpha, increases to the stall angle, alpha(sub stall) the flow separation point on the upper surface of the wing moves to the leading edge, so that on a two-dimensional airfoil or a large-aspect-ratio wing, the lift abruptly drops to a very low level. Therefore, the first generation of aeronautical flow type, i.e., the attached steady flow, has been limited to alpha less than alpha(sub stall). Owing to the obvious importance in applications, therefore, a great effort has been made in the past two decades to enlarge the range of usable angles of attack by various flow controls for a large-aspect-ratio wing. Basically, relevant works fall into two categories. The first category is usually refereed to as separation control, which concentrates on partially separated flow at alpha less than alpha(sub stall). Since the first experimental study of Collins and Zelenevitz, there has been ample literature showing that a partially separated flow can be turned to almost fully attached by flow controls, so that the lift is recovered and the stall is delayed (for a recent work see Seifert et al.). It has been well established that, in this category, unsteady controls are much more effective than steady ones and can be realized at a very low power-input level (Wu et al.; Seifert et al.). The second and more ambitious category of relevant efforts is the post-stall lift enhancement. Its possibility roots at the existence of a second lift peak at a very high angle of attack. In fact, As alpha further increases from alpha(sub stall), the completely separated flow develops and gradually becomes a bluff-body flow. This flow gives a normal force to the airfoil with a lift component, which reaches a peak at a maximum utilizable angle of attack, alpha(sub m) approx.= 40 deg. This second peak is of the same level as the first lift peak at alpha(sub stall). Meanwhile, the drag is also quickly increased (e.g., Fage and Johansen ; Critzos et al.). Figure 1 shows a typical experimental lift and drag coefficients of NACA-0012 airfoil in this whole range of angle of attack. Obviously, without overcoming the lift crisis at alpha(sub stall) the second lift peak is completely useless. Thus, the ultimate goal of post-stall lift enhancement is to fill the lift valley after stall by flow controls, so that a wing and/or flap can work at the whole range of 0 deg less than alpha less than alpha(sub m). Relevant early experimental studies have been extensively reviewed by Wu et al., who concluded that, first, similar to the leading-edge vortex on a slender wing, the lift enhancement on a large-aspect-ratio wing should be the result of capturing a vortex on the upper surface of the wing; and, second, using steady controls cannot reach the goal, and one must rely on unsteady controls with low-level power input as well. Wu et al. also conjectured that the underlying physics of post-stall lift enhancement by unsteady controls consists of a chain of mechanisms: vortex layer instability - receptivity resonance - nonlinear streaming.

Wu, J. M.↗