Numerical Method for Multi-Dimensional Freezing Problems in Arbitrary Domains | Synapse