A two-dimensional nonlinear PBL numerical model using an energy closure (E-ε)method has been employed to study the sea breeze circulation and TIBL in coastal areas. The main characteristics of sea breeze obtained from numerical experiments agree with those from general observation facts. The depth of sea breeze ranges from 300 to 900 meters, maximum velocity from 1,5 m/s to 4 m/s, and its height from 100 m to 300 m. The agreement between comparisons reveals that the performance of the model is good, and the selected experiment conditions are reasonable. This paper refits the function of TIBL profiles using the Weisman’s formula and the exponent value is considered to change with the different states of sea breeze. Numerical experiment results indicate that the exponent of the TIBL profile, ranging from 0.4 to 1. 1, is related to the strength and depth of the sea breeze. The exponent of 0.5 is suitable only when sea breeze is fully developed. This paper also gives various exponents under different sea breezes. A two-dimensional nonlinear PBL numerical model using an energy closure (E-ε) method has been employed to study the sea breeze circulation and TIBL in coastal areas. The main characteristics of sea breeze obtained from numerical experiments agree with those from general observation facts The depth of the sea breeze ranges from 300 to 900 meters, maximum velocity from 1,5 m / s to 4 m / s, and its height from 100 m to 300 m. The agreement between comparisons reveals that the performance of the model is This paper refits the function of TIBL profiles using the Weisman’s formula and the exponent value is considered to change with the different states of sea breeze. Numerical experiment results indicate that the exponent of the TIBL profile, ranging from 0.4 to 1. 1, is related to the strength and depth of the sea breeze. The exponent of 0.5 is suitable only when sea breeze is fully developed. This paper also gives various exponents und er different sea breezes.