In the vehicle suspension system field, smart materials are contributing their top-notch characteristics in building up magneto-rheological (MR) dampers. Though MR fluid has good characteristics to absorb shock, sedimentation problem is a significant issue while performing an MR damper. MR grease eliminates the sedimentation problem and can be used a MR damper in vehicle suspensions system, in vibration control of large structures, biomedical, etc. This paper aims to conduct a numerical study based on coupled FEA and CFD and evaluate the damping force induced in the MR grease gap region for a different input current ranging from 0.1A to 0.4A at frequencies 1 Hz and 2 Hz. ANSYS CFX solver has been used with SST k-ω turbulence model and ‘Herschel-Bulkley’ model for the numerical simulation. At first the simulation approach/model was validated with the experimental results. Initially, the simulation was assumed to be a laminar flow, but as the piston of the damper moves for the stroke, the flow becomes turbulent. Without Shear Stress Transport (SST) turbulence model the simulation was unable to capture laminar to turbulent transition and suddenly stopped when the transition occurred without writing any result. Also, a comparison analysis between ‘Herschel-Bulkley’ fluid model and ‘Bingham’ model was done to prove that only the ‘Herschel-Bulkley’ fluid model can capture the rectangular curve of the damping force for MR Grease as it is stated by the authors who validated the experimental results for MR Grease. No numerical simulation has been found for characterizing MR Grease damper. The FEA has been used to determine the fluid properties due to various currents and frequencies, which were input for the CFD analysis to determine the damping forces and Quasi-static analysis. This simulation model could capture the hysteresis behavior correctly also. [ABSTRACT FROM AUTHOR]