Abstract: A mathematic model for cooling water discharge to forked bays is established by using the boundary fitted orthogonal coordinate system to simulate the complex boundaries, the full implicit coupled solution model to discrete the base equations, and the matrix chase-after method to solve the algebraic equations. Examples show that the envelope method to evaluate scales of thermal pollution can be applied to assess the engineering design and environment.