This paper presents a new approach to verify the accuracy of computational simulations. We develop mathematical theorems which can serve as robust a posteriori error estimation techniques to identify numerical pollution, check the performance of adaptive meshes, and verify numerical solutions. We demonstrate performance of this methodology on problems from flow thorough porous media. However, one can extend it to other models. We construct mathematical properties such that the solutions to Darcy and Darcy-Brinkman equations satisfy them. The mathematical properties include the total minimum mechanical power, minimum dissipation theorem, reciprocal relation, and maximum principle for the vorticity. All the developed theorems have firm mechanical bases and are independent of numerical methods. So, these can be utilized for solution verification of finite element, finite volume, finite difference, lattice Boltzmann methods and so forth. In particular, we show that, for a given set of boundary conditions, Darcy velocity has the minimum total mechanical power of all the kinematically admissible vector fields. We also show that a similar result holds for Darcy-Brinkman velocity. We then show for a conservative body force, the Darcy and Darcy-Brinkman velocities have the minimum total dissipation among their respective kinematically admissible vector fields. Using numerical examples, we show that the minimum dissipation and total mechanical power theorems can be utilized to identify pollution errors in numerical solutions. The solutions to Darcy and Darcy-Brinkman equations are shown to satisfy a reciprocal relation, which has the potential to identify errors in the numerical implementation of boundary conditions. It is also shown that the vorticity under both steady and transient Darcy-Brinkman equations satisfy maximum principles if the body force is conservative and the permeability is homogeneous and isotropic. A discussion on the nature of vorticity under steady and transient Darcy equations is also presented. Using several numerical examples, we will demonstrate the predictive capabilities of the proposed a posteriori techniques in assessing the accuracy of numerical solutions for a general class of problems, which could involve complex domains and general computational grids.

In recent decades, problems related to the squeeze of fluid films in the presence of aporous medium draw attention of researchers and are the subject of many applied studiesfor industry and biomechanics. Our concerns in this paper are the numerical simulation ofthe viscous shear stresses effects on the fluid film characteristics between two discswith one porous. This study is based on the coupling, at the fluid film-porous discinterface, of the Darcy-Brinkman equations in the porous medium and the modified Reynoldsequation describing the flow in the fluid film. The system of equations obtained isdiscretized by the means of finite differences method and solved numerically using thetechnique of Successive Over-Relaxation (SOR). The results show that the viscous sheareffects increase the radial and the axial fluid film velocities as well as the squeezefilm velocity but decrease the response time. Moreover, these effects are enlarged forsmaller viscous shear parameter and for smaller fluid film thickness.