In fluid dynamics, Beltrami flows are flows in which the vorticity vector 
  
    
      
        
          ω
        
      
    
    {\displaystyle \mathbf {\omega } }
  
 and the velocity vector 
  
    
      
        
          v
        
      
    
    {\displaystyle \mathbf {v} }
  
 are parallel to each other. In other words, Beltrami flow is a flow where Lamb vector is zero. It is named after the Italian mathematician Eugenio Beltrami due to his derivation of the Beltrami vector field, while initial developments in fluid dynamics were done by the Russian scientist Ippolit S. Gromeka in 1881.
 Since the vorticity vector 
  
    
      
        
          ω
        
      
    
    {\displaystyle {\boldsymbol {\omega }}}
  
 and the velocity vector 
  
    
      
        
          v
        
      
    
    {\displaystyle \mathbf {v} }
  
 are parallel to each other, we can write
 where 
  
    
      
        α
        (
        
          x
        
        ,
        t
        )
      
    
    {\displaystyle \alpha (\mathbf {x} ,t)}
  
 is some scalar function. One immediate consequence of Beltrami flow is that it can never be a planar or axisymmetric  flow because in those flows, vorticity is always perpendicular to the velocity field. The other important consequence will be realized by looking at the incompressible vorticity equation
 where 
  
    
      
        
          f
        
      
    
    {\displaystyle \mathbf {f} }
  
 is an external body forces such as gravitational field, electric field etc., and 
  
    
      
        ν
      
    
    {\displaystyle \nu }
  
 is the kinematic viscosity. Since 
  
    
      
        
          ω
        
      
    
    {\displaystyle {\boldsymbol {\omega }}}
  
 and 
  
    
      
        
          v
        
      
    
    {\displaystyle \mathbf {v} }
  
 are parallel, the non-linear terms in the above equation are identically zero 
  
    
      
        (
        
          v
        
        ⋅
        ∇
        )
        
          ω
        
        =
        (
        
          ω
        
        ⋅
        ∇
        )
        
          v
        
        =
        0
      
    
    {\displaystyle (\mathbf {v} \cdot \nabla ){\boldsymbol {\omega }}=({\boldsymbol {\omega }}\cdot \nabla )\mathbf {v} =0}
  
. Thus Beltrami flows satisfies the linear equation
 When 
  
    
      
        
          f
        
        =
        0
      
    
    {\displaystyle \mathbf {f} =0}
  
, the components of vorticity satisfies a simple heat equation.
 Viktor Trkal considered the Beltrami flows without any external forces in 1919 for the scalar function 
  
    
      
        α
        (
        
          x
        
        ,
        t
        )
        =
        c
        =
        
          constant
        
      
    
    {\displaystyle \alpha (\mathbf {x} ,t)=c={\text{constant}}}
  
, i.e.,
 Introduce the following separation of variables
 then the equation satisfied by 
  
    
      
        
          g
        
        (
        
          x
        
        )
      
    
    {\displaystyle \mathbf {g} (\mathbf {x} )}
  
 becomes
 The Chandrasekhar–Kendall functions satisfy this equation.
 Ratip Berker obtained the solution in Cartesian coordinates for 
  
    
      
        
          g
        
        (
        
          x
        
        )
      
    
    {\displaystyle \mathbf {g} (\mathbf {x} )}
  
 in 1963,[dubious  – discuss]
 The generalized Beltrami flow satisfies the condition
 which is less restrictive than the Beltrami condition 
  
    
      
        
          v
        
        ×
        
          ω
        
        =
        0
      
    
    {\displaystyle \mathbf {v} \times {\boldsymbol {\omega }}=0}
  
. Unlike the normal Beltrami flows, the generalized Beltrami flow can be studied for planar and axisymmetric flows.
 For steady generalized Beltrami flow, we have 
  
    
      
        
          ∇
          
            2
          
        
        
          ω
        
        =
        0
        ,
         
        ∇
        ×
        (
        
          v
        
        ×
        
          ω
        
        )
        =
        0
      
    
    {\displaystyle \nabla ^{2}{\boldsymbol {\omega }}=0,\ \nabla \times (\mathbf {v} \times {\boldsymbol {\omega }})=0}
  
 and since it is also planar we have 
  
    
      
        
          v
        
        =
        (
        u
        ,
        v
        ,
        0
        )
        ,
         
        
          ω
        
        =
        (
        0
        ,
        0
        ,
        ζ
        )
      
    
    {\displaystyle \mathbf {v} =(u,v,0),\ {\boldsymbol {\omega }}=(0,0,\zeta )}
  
. Introduce the stream function
 Integration of 
  
    
      
        ∇
        ×
        (
        
          v
        
        ×
        
          ω
        
        )
        =
        0
      
    
    {\displaystyle \nabla \times (\mathbf {v} \times {\boldsymbol {\omega }})=0}
  
 gives 
  
    
      
        ζ
        =
        −
        f
        (
        ψ
        )
      
    
    {\displaystyle \zeta =-f(\psi )}
  
. So, complete solution is possible if it satisfies all the following three equations
 A special case is considered when the flow field has uniform vorticity 
  
    
      
        f
        (
        ψ
        )
        =
        C
        =
        
          constant
        
      
    
    {\displaystyle f(\psi )=C={\text{constant}}}
  
. Wang (1991) gave the generalized solution as
 assuming a linear function for 
  
    
      
        A
        (
        x
        ,
        y
        )
      
    
    {\displaystyle A(x,y)}
  
. Substituting this into the vorticity equation and introducing the separation of variables 
  
    
      
        ψ
        (
        x
        ,
        y
        )
        =
        X
        (
        x
        )
        Y
        (
        y
        )
      
    
    {\displaystyle \psi (x,y)=X(x)Y(y)}
  
 with the separating constant 
  
    
      
        
          λ
          
            2
          
        
      
    
    {\displaystyle \lambda ^{2}}
  
 results in
 The solution obtained for different choices of 
  
    
      
        a
        ,
         
        b
        ,
         
        λ
      
    
    {\displaystyle a,\ b,\ \lambda }
  
 can be interpreted differently, for example, 
  
    
      
        a
        =
        0
        ,
         
        b
        =
        −
        U
        ,
        
          λ
          
            2
          
        
        >
        0
      
    
    {\displaystyle a=0,\ b=-U,\lambda ^{2}>0}
  
 represents a flow downstream a uniform grid, 
  
    
      
        a
        =
        −
        U
        ,
         
        b
        =
        0
        ,
        
          λ
          
            2
          
        
        =
        0
      
    
    {\displaystyle a=-U,\ b=0,\lambda ^{2}=0}
  
 represents a flow created by a stretching plate, 
  
    
      
        a
        =
        −
        U
        ,
         
        b
        =
        U
        ,
        
          λ
          
            2
          
        
        =
        0
      
    
    {\displaystyle a=-U,\ b=U,\lambda ^{2}=0}
  
 represents a flow into a corner, 
  
    
      
        a
        =
        −
        V
        ,
         
        b
        =
        −
        U
        ,
        
          λ
          
            2
          
        
        =
        0
      
    
    {\displaystyle a=-V,\ b=-U,\lambda ^{2}=0}
  
 represents an Asymptotic suction profile etc.
 Here,
 G. I. Taylor gave the solution for a special case where 
  
    
      
        ζ
        =
        K
        ψ
      
    
    {\displaystyle \zeta =K\psi }
  
, where 
  
    
      
        K
      
    
    {\displaystyle K}
  
 is a constant in 1923.  He showed that the separation 
  
    
      
        ψ
        =
        
          e
          
            −
            ν
            λ
            t
          
        
        Ψ
        (
        x
        ,
        y
        )
      
    
    {\displaystyle \psi =e^{-\nu \lambda t}\Psi (x,y)}
  
 satisfies the equation and also
 Taylor also considered an example, a decaying system of eddies rotating alternatively in opposite directions and arranged in a rectangular array
 which satisfies the above equation with 
  
    
      
        λ
        =
        2
        
          π
          
            2
          
        
        
          /
        
        
          d
          
            2
          
        
      
    
    {\displaystyle \lambda =2\pi ^{2}/d^{2}}
  
, where 
  
    
      
        d
      
    
    {\displaystyle d}
  
 is the length of the square formed by an eddy. Therefore, this system of eddies decays as
 O. Walsh generalized Taylor's eddy solution in 1992. Walsh's solution is of the form 
  
    
      
        
          v
        
        =
        
          e
          
            −
            ν
            λ
            t
          
        
        
          u
        
        (
        x
        ,
        y
        )
      
    
    {\displaystyle \mathbf {v} =e^{-\nu \lambda t}\mathbf {u} (x,y)}
  
, where 
  
    
      
        
          ∇
          
            2
          
        
        
          u
        
        =
        −
        λ
        
          u
        
      
    
    {\displaystyle \nabla ^{2}\mathbf {u} =-\lambda \mathbf {u} }
  
 and 
  
    
      
        ∇
        ⋅
        
          u
        
        =
        0.
      
    
    {\displaystyle \nabla \cdot \mathbf {u} =0.}
  

 Here we have 
  
    
      
        
          v
        
        =
        
          (
          
            
              u
              
                r
              
            
            ,
            0
            ,
            
              u
              
                z
              
            
          
          )
        
        ,
         
        
          ω
        
        =
        (
        0
        ,
        ζ
        ,
        0
        )
      
    
    {\displaystyle \mathbf {v} =\left(u_{r},0,u_{z}\right),\ {\boldsymbol {\omega }}=(0,\zeta ,0)}
  
. Integration of 
  
    
      
        ∇
        ×
        (
        
          v
        
        ×
        
          ω
        
        )
        =
        0
      
    
    {\displaystyle \nabla \times (\mathbf {v} \times {\boldsymbol {\omega }})=0}
  
 gives 
  
    
      
        ζ
        =
        r
        f
        (
        ψ
        )
      
    
    {\displaystyle \zeta =rf(\psi )}
  
 and the three equations are
 The first equation is the Hicks equation. Marris and Aswani (1977) showed that the only possible solution is 
  
    
      
        f
        (
        ψ
        )
        =
        C
        =
        
          constant
        
      
    
    {\displaystyle f(\psi )=C={\text{constant}}}
  
 and the remaining equations reduce to
 A simple set of solutions to the above equation is
 
  
    
      
        
          c
          
            1
          
        
        ,
        
          c
          
            4
          
        
        ≠
        0
        ,
         
        
          c
          
            2
          
        
        ,
        
          c
          
            3
          
        
        ,
        
          c
          
            5
          
        
        =
        0
      
    
    {\displaystyle c_{1},c_{4}\neq 0,\ c_{2},c_{3},c_{5}=0}
  
 represents a flow due to two opposing rotational stream on a parabolic surface, 
  
    
      
        
          c
          
            2
          
        
        ≠
        0
        ,
        
          c
          
            1
          
        
        ,
        
          c
          
            3
          
        
        ,
        
          c
          
            4
          
        
        ,
        
          c
          
            5
          
        
        =
        0
      
    
    {\displaystyle c_{2}\neq 0,c_{1},c_{3},c_{4},c_{5}=0}
  
 represents rotational flow on a plane wall, 
  
    
      
        
          c
          
            1
          
        
        ,
        
          c
          
            2
          
        
        ,
        
          c
          
            3
          
        
        ≠
        0
        ,
         
        
          c
          
            4
          
        
        ,
        
          c
          
            5
          
        
        =
        0
      
    
    {\displaystyle c_{1},c_{2},c_{3}\neq 0,\ c_{4},c_{5}=0}
  
 represents a flow ellipsoidal vortex (special case – Hill's spherical vortex), 
  
    
      
        
          c
          
            1
          
        
        ,
        
          c
          
            3
          
        
        ,
        
          c
          
            5
          
        
        ≠
        0
        ,
         
        
          c
          
            2
          
        
        ,
        
          c
          
            4
          
        
        =
        0
      
    
    {\displaystyle c_{1},c_{3},c_{5}\neq 0,\ c_{2},c_{4}=0}
  
 represents a type of toroidal vortex etc.
 The homogeneous solution for 
  
    
      
        C
        =
        0
      
    
    {\displaystyle C=0}
  
 as shown by Berker
 where 
  
    
      
        
          J
          
            1
          
        
        (
        k
        r
        )
        ,
        
          Y
          
            1
          
        
        (
        k
        r
        )
      
    
    {\displaystyle J_{1}(kr),Y_{1}(kr)}
  
 are the Bessel function of the first kind and Bessel function of the second kind respectively. A special case of the above solution is Poiseuille flow for cylindrical geometry with transpiration velocities on the walls. Chia-Shun Yih found a solution in 1958 for Poiseuille flow into a sink when 
  
    
      
        C
        =
        2
        ,
        
        
          c
          
            1
          
        
        =
        −
        1
        
          /
        
        4
        ,
        
        
          c
          
            3
          
        
        =
        1
        
          /
        
        2
        ,
        
        
          c
          
            2
          
        
        =
        
          c
          
            4
          
        
        =
        
          c
          
            5
          
        
        =
        
          B
          
            k
          
        
        =
        
          C
          
            k
          
        
        =
        0
      
    
    {\displaystyle C=2,\,c_{1}=-1/4,\,c_{3}=1/2,\,c_{2}=c_{4}=c_{5}=B_{k}=C_{k}=0}
  
.
 Beltrami fields are a classical steady solution to the Euler equation. Beltrami fields play an important role in (ideal) fluid mechanics in equilibrium, as complexity is only expected for these fields.