Cellular foams have found a predominant application in automotive industry for efficient energy absorption to meet stringent and continuously improving vehicle crashworthiness and occupant protection criteria. The recent inclusion of pedestrian protection regulations mandate the use of foams of different densities for impact energy absorption at identified impact locations; this has paved way for significant advancements in molding techniques such as dual density and tri-density molding. With increased emphasis on light-weighting for improved fuel-efficiency, solutions involving the use of polymeric or metallic foams as fillers in hollow structures and foam encapsulated metal structures are being explored. Another major automotive application of foams is in the seat comfort area which again involves intricate shapes and sizes. Also, a few recently developed foams are anisotropic, adding on to the existing complexities. In addition, complexities associated with controlled/ uncontrolled spatial variation in density and the geometry of molded parts and use of foams in sandwich composites offer several challenges for the CAE community to model the foam parts. As a first step to capture these complexities, we need utilize the available LS-DYNA modeling features to enable effective finite element analysis (FEA) of foam components optimally. This paper aims to investigate the various underlying parameters such as element formulations and size, contact stiffness and hourglass control, governing stability and accuracy of foam models, with the aim of identifying optimal settings of these parameters. The optimal settings for the identified parameters are elucidated in the context of LS-DYNA *MAT_FU-CHANG_FOAM; however, these learning apply equally to any LS-DYNA material model for cellular materials with zero Poisson’s ratio.