Bounded domain modeling with specified boundary conditions and mass balancing

Computational fluid dynamics modeling of a bounded domain is provided which includes solving iteratively a computational fluid dynamics model of the bounded domain using mass flow boundary conditions that are specified. The solving includes automatically adjusting the specified mass flow boundary conditions for at least one iteration of the solving. The automatically adjusting includes mass balancing flows at a boundary of the bounded domain to provide corrected mass flow boundary conditions for the solving. The mass balancing includes applying different corrections to mass in-flow across the domain boundary, compared with mass out-flow across the boundary. The mass balancing may include multiplying mass in-flows across the boundary by (1−CF), and mass out-flows across the boundary by (1+CF), where CF is a determined correction factor, and net mass flow is positive for mass out-flow exiting the domain, and negative for mass in-flow entering the domain.