The foaming of hydrogels presents a promising strategy for tailoring mechanical and radiological properties to replicate biological soft tissues for biomedical phantoms. Achieving uniform and predictable void fraction distributions in alginate hydrogel foams remains a challenge due to the complex interplay between bubble dynamics, matrix rheology, and interfacial forces during the pre-gelation aeration stage. This thesis develops a Computational Fluid Dynamics framework using the transient Euler