The linear operator plays an important role in the computational process of Homotopy Analysis Method (HAM). In HAM frame any kind of linear operator can be chosen to develop a solution. Hence, it is easy to introduce the modified/physical parameter dependent linear operators. The effective use of these operators has been judged through solving fluid flow problems. Modification in linear operators affects the solution and improves the computational efficiency of HAM for larger values of parameters. The convergence rate of the solution is rapid and several times higher resulting in lesser computational time.

We couple different flow models, i.e. a finite element solver for the Navier-Stokes equations and a Lattice Boltzmann automaton, using the framework Peano as a common base. The new coupling strategy between the meso- and macroscopic solver is presented and validated in a 2D channel flow scenario. The results are in good agreement with theory and results obtained in similar works by Latt et al. In addition, the test scenarios show an improved stability of the coupled method compared to pure Lattice Boltzmann simulations.