I'd like to know if anyone has any books/links /article on "solving boussinesq equations using numerical methods". Particularly some stuff on "MOL( Method of Lines)".
Very generally speaking, the Method of Lines approach is to first spatially discretize the PDEs (in this case, the Boussinesq equations) and then use a time-integrator to advance the resulting system of ODEs. The choice of discretization method depends on the geometry and the Prandtl number of the Boussinesq equations. For example, people use finite differences and even Chebyshev pseudospectral methods for this purpose. However, the equations may have complicated boundary conditions, in which case people use things like the Influence Matrix method.
For solving Boussinesq equations with an infinite Prandtl number in a spherical shell, see Yoshida, M., and A. Kageyama (2004), Application of the Yin‐Yang grid to a thermal convection of a Boussinesq fluid
with infinite Prandtl number in a three-dimensional spherical shell. See also this paper by Wright, Flyer and Yuen which combines RBF-PS and Chebyshev-PS to accomplish this (http://scholarworks.boisestate.edu/cgi/viewcontent.cgi?article=1027&context=math_facpubs)
Very generally speaking, the Method of Lines approach is to first spatially discretize the PDEs (in this case, the Boussinesq equations) and then use a time-integrator to advance the resulting system of ODEs. The choice of discretization method depends on the geometry and the Prandtl number of the Boussinesq equations. For example, people use finite differences and even Chebyshev pseudospectral methods for this purpose. However, the equations may have complicated boundary conditions, in which case people use things like the Influence Matrix method.
For solving Boussinesq equations with an infinite Prandtl number in a spherical shell, see Yoshida, M., and A. Kageyama (2004), Application of the Yin‐Yang grid to a thermal convection of a Boussinesq fluid
with infinite Prandtl number in a three-dimensional spherical shell. See also this paper by Wright, Flyer and Yuen which combines RBF-PS and Chebyshev-PS to accomplish this (http://scholarworks.boisestate.edu/cgi/viewcontent.cgi?article=1027&context=math_facpubs)
The Boussines equations are implemented in the OpenFoam software packeage, which is open source. With some experience in CFD you will get very fast nice results.
An efficient numerical technique is used to obtain solutions of the Boussinesq equation for problems of groundwater recession and groundwater flow in response to changes in stream stage. Comparison of the numerical results with the analytical solution of Boussinesq [1904] attests to the accuracy of the former. The solution of a recession problem indicates that the form of the Werner and Sundquist [1951] model of groundwater discharge recession is appropriate. The solution of a problem involving changing boundary conditions on an aquifer because of a flood wave provides data relative to groundwater outflow. A comparison with the results of Cooper and Rorabaugh [1963] evaluates the applicability of a linear model of unconfined flow.