Acoustic streaming is a steady flow in a fluid driven by the absorption of high amplitude acoustic oscillations.  This phenomenon can be observed near sound emitters, or in the standing waves within a Kundt's tube. Acoustic streaming was explained first by Lord Rayleigh in 1884.
It is the less-known opposite of sound generation by a flow.
 There are two situations where sound is absorbed in its medium of propagation:
 Consider a plane standing sound wave that corresponds to the velocity field 
  
    
      
        U
        (
        x
        ,
        t
        )
        =
        
          v
          
            0
          
        
        cos
        ⁡
        k
        x
        cos
        ⁡
        ω
        t
        =
        ε
        cos
        ⁡
        k
        x
        ℜ
        (
        
          e
          
            −
            i
            ω
            t
          
        
        )
      
    
    {\displaystyle U(x,t)=v_{0}\cos kx\cos \omega t=\varepsilon \cos kx\Re (e^{-i\omega t})}
  
 where 
  
    
      
        k
        =
        2
        π
        
          /
        
        λ
        =
        ω
        
          /
        
        c
      
    
    {\displaystyle k=2\pi /\lambda =\omega /c}
  
. Let the characteristic (transverse) dimension of the problem be 
  
    
      
        l
      
    
    {\displaystyle l}
  
. The flow field just described corresponds to inviscid flow. However viscous effects will be important close to a solid wall; there then exists a boundary layer of thickness or, penetration depth 
  
    
      
        δ
        =
        (
        2
        ν
        
          /
        
        ω
        
          )
          
            1
            
              /
            
            2
          
        
      
    
    {\displaystyle \delta =(2\nu /\omega )^{1/2}}
  
. Rayleigh streaming is best visualized in the approximation 
  
    
      
        λ
        ≫
        l
        ≫
        δ
        .
      
    
    {\displaystyle \lambda \gg l\gg \delta .}
  
 As in 
  
    
      
        U
        (
        x
        ,
        t
        )
      
    
    {\displaystyle U(x,t)}
  
, the velocity components 
  
    
      
        (
        u
        ,
        v
        )
      
    
    {\displaystyle (u,v)}
  
 are much less than 
  
    
      
        c
      
    
    {\displaystyle c}
  
. In addition, the characteristic time scale within the boundary layer is very large (because of the smallness of 
  
    
      
        δ
      
    
    {\displaystyle \delta }
  
) in comparison with the acoustic time scale 
  
    
      
        l
        
          /
        
        c
      
    
    {\displaystyle l/c}
  
. These observations imply that the flow in the boundary layer may be regarded as incompressible. 
 The unsteady, incompressible boundary-layer equation is
 where the right-hand side terms correspond to the pressure gradient imposed on the boundary layer. The problem can be solved using the stream function 
  
    
      
        ψ
      
    
    {\displaystyle \psi }
  
 that satisfies 
  
    
      
        u
        =
        ∂
        ψ
        
          /
        
        ∂
        y
      
    
    {\displaystyle u=\partial \psi /\partial y}
  
 and 
  
    
      
        v
        =
        −
        ∂
        ψ
        
          /
        
        ∂
        x
        .
      
    
    {\displaystyle v=-\partial \psi /\partial x.}
  
 Since by definition, velocity field 
  
    
      
        U
      
    
    {\displaystyle U}
  
 in the sound wave is very small, we can formally obtain the solution for the boundary layer equation by introducing the asymptotic series for 
  
    
      
        ε
        →
        0
      
    
    {\displaystyle \varepsilon \rightarrow 0}
  
 as 
  
    
      
        u
        =
        ε
        
          u
          
            1
          
        
        +
        
          ε
          
            2
          
        
        
          u
          
            2
          
        
        +
        ⋯
      
    
    {\displaystyle u=\varepsilon u_{1}+\varepsilon ^{2}u_{2}+\cdots }
  
, 
  
    
      
        ψ
        =
        ε
        
          ψ
          
            1
          
        
        +
        
          ε
          
            2
          
        
        
          ψ
          
            2
          
        
        ⋯
      
    
    {\displaystyle \psi =\varepsilon \psi _{1}+\varepsilon ^{2}\psi _{2}\cdots }
  
 etc.
 In the first approximation, one obtains
 The solution that satisfies the no-slip condition at the wall 
  
    
      
        y
        
          /
        
        δ
        =
        0
      
    
    {\displaystyle y/\delta =0}
  
 and approaches 
  
    
      
        U
      
    
    {\displaystyle U}
  
 as 
  
    
      
        y
        
          /
        
        δ
        →
        ∞
      
    
    {\displaystyle y/\delta \rightarrow \infty }
  
 is given by
 where 
  
    
      
        κ
        =
        (
        1
        −
        i
        )
        
          /
        
        δ
      
    
    {\displaystyle \kappa =(1-i)/\delta }
  
 and 
  
    
      
        
          ζ
          
            1
          
        
        =
        y
        +
        (
        
          e
          
            −
            κ
            y
          
        
        −
        1
        )
        
          /
        
        κ
        .
      
    
    {\displaystyle \zeta _{1}=y+(e^{-\kappa y}-1)/\kappa .}
  
 
 The equation at the next order is
 Since each term on the right-hand side is quadratic, it will result in terms with frequencies 
  
    
      
        ω
        +
        ω
        =
        2
        ω
      
    
    {\displaystyle \omega +\omega =2\omega }
  
 and 
  
    
      
        ω
        −
        ω
        =
        0.
      
    
    {\displaystyle \omega -\omega =0.}
  
 The 
  
    
      
        ω
        =
        0
      
    
    {\displaystyle \omega =0}
  
 terms correspond to time independent forcing for 
  
    
      
        
          u
          
            2
          
        
      
    
    {\displaystyle u_{2}}
  
. Let us find solution that corresponds only to this time-independent part. This leads to 
  
    
      
        
          ψ
          
            2
          
        
        =
        sin
        ⁡
        2
        k
        x
        
        
          ζ
          
            2
          
        
        (
        y
        )
        
          /
        
        c
      
    
    {\displaystyle \psi _{2}=\sin 2kx\,\zeta _{2}(y)/c}
  
 where 
  
    
      
        
          ζ
          
            2
          
        
      
    
    {\displaystyle \zeta _{2}}
  
 satisfies the equation
 where prime denotes differentiation with respect to 
  
    
      
        y
        .
      
    
    {\displaystyle y.}
  
 The boundary condition at the wall implies that 
  
    
      
        ζ
        (
        0
        )
        =
        
          ζ
          ′
        
        (
        0
        )
        =
        0.
      
    
    {\displaystyle \zeta (0)=\zeta '(0)=0.}
  
 As 
  
    
      
        y
        
          /
        
        δ
        →
        ∞
      
    
    {\displaystyle y/\delta \rightarrow \infty }
  
, 
  
    
      
        
          ζ
          
            2
          
        
      
    
    {\displaystyle \zeta _{2}}
  
 must be finite. Integrating the above equation twice gives
 As 
  
    
      
        y
        
          /
        
        δ
        →
        ∞
      
    
    {\displaystyle y/\delta \rightarrow \infty }
  
, 
  
    
      
        
          ζ
          ′
        
        (
        ∞
        )
        =
        3
        
          /
        
        8
      
    
    {\displaystyle \zeta '(\infty )=3/8}
  
 leading to the result that 
  
    
      
        
          v
          
            2
          
        
        (
        x
        ,
        ∞
        ,
        t
        )
        =
        (
        3
        
          /
        
        8
        c
        )
        sin
        ⁡
        2
        k
        x
        .
      
    
    {\displaystyle v_{2}(x,\infty ,t)=(3/8c)\sin 2kx.}
  
 Thus, at the edge of the boundary, there is a steady fluid motion superposed on the oscillating motion. This velocity forcing will drive a steady streaming motion outside the boundary layer. The interesting result is that since 
  
    
      
        
          v
          
            2
          
        
        (
        ∞
        )
      
    
    {\displaystyle v_{2}(\infty )}
  
 is independent of 
  
    
      
        ν
      
    
    {\displaystyle \nu }
  
, the steady streaming motion happening outside the boundary layer is also independent of viscosity, although its origin of existence due to the viscous boundary layer. 
 The outer steady streaming incompressible motion will depend on the geometry of the problem. If there are two walls one at 
  
    
      
        y
        =
        0
      
    
    {\displaystyle y=0}
  
 and 
  
    
      
        y
        =
        2
        h
      
    
    {\displaystyle y=2h}
  
, then the solution is
 which corresponds a periodic array of counter-rotating vortices, as shown in the figure.  
 Acoustic streaming is a non-linear effect.  

We can decompose the velocity field in a vibration part and a steady part  
  
    
      
        
          u
        
        =
        v
        +
        
          
            u
            ¯
          
        
      
    
    {\displaystyle {u}=v+{\overline {u}}}
  
.
The vibration part 
  
    
      
        v
      
    
    {\displaystyle v}
  
 is due to sound, while the steady part is the acoustic streaming velocity (average velocity).
The Navier–Stokes equations implies for the acoustic streaming velocity:
 The steady streaming originates from a steady body force 
  
    
      
        
          f
          
            i
          
        
        =
        −
        
          ∂
        
        (
        
          
            
              ρ
              
                v
                
                  i
                
              
              
                v
                
                  j
                
              
            
            ¯
          
        
        )
        
          /
        
        
          ∂
          
            x
            
              j
            
          
        
      
    
    {\displaystyle f_{i}=-{\partial }({\overline {\rho v_{i}v_{j}}})/{\partial x_{j}}}
  
 that appears on the right hand side. This force is a function of what is known as the Reynolds stresses in turbulence 
  
    
      
        −
        
          
            
              ρ
              
                v
                
                  i
                
              
              
                v
                
                  j
                
              
            
            ¯
          
        
      
    
    {\displaystyle -{\overline {\rho v_{i}v_{j}}}}
  
. The Reynolds stress depends on the amplitude of sound vibrations, and the body force reflects diminutions in this sound amplitude.
 We see that this stress is non-linear (quadratic)  in the velocity amplitude. It is non-vanishing only where the velocity amplitude varies. 
If  the velocity of the fluid oscillates because of sound as 
  
    
      
        ϵ
        cos
        ⁡
        (
        ω
        t
        )
      
    
    {\displaystyle \epsilon \cos(\omega t)}
  
,  the quadratic non-linearity generates a  steady force proportional to 

  
    
      
        
          
            
              
                
                  ϵ
                  
                    2
                  
                
                
                  cos
                  
                    2
                  
                
                ⁡
                (
                ω
                t
                )
              
              ¯
            
          
          =
          
            ϵ
            
              2
            
          
          
            /
          
          2
        
      
    
    {\displaystyle \scriptstyle {\overline {\epsilon ^{2}\cos ^{2}(\omega t)}}=\epsilon ^{2}/2}
  
.
 Even if viscosity is responsible for acoustic streaming, the value of viscosity disappears from the resulting streaming velocities in the case of near-boundary acoustic steaming.
 The order of magnitude of streaming velocities are:
 with 
  
    
      
        
          v
          
            0
          
        
      
    
    {\displaystyle v_{0}}
  
 the sound vibration velocity and 
  
    
      
        x
      
    
    {\displaystyle x}
  
 along the wall boundary. The flow is directed towards decreasing sound vibrations (vibration nodes).