A new hybrid method, Monte-Carlo-Heat-Flux (MCHF) method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo meth...A new hybrid method, Monte-Carlo-Heat-Flux (MCHF) method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo method with the heat flux method. Its accuracy and reliability was proved by comparing the computational results with exact results from classical "Zone Method".展开更多
This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip veloci...This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.展开更多
During the second course of USA - PRC joint air sea interaction experiment in 1986, the temperature structure parameters CT2 were measured by sodar over the Western Pacific Ocean. Based on similarity theory, a method ...During the second course of USA - PRC joint air sea interaction experiment in 1986, the temperature structure parameters CT2 were measured by sodar over the Western Pacific Ocean. Based on similarity theory, a method is discussed to calculate the sensible heat flux over the ocean in unstable stratification. Becausehumidity is great over the ocean, so we have to consider the influence of water vapor structure parameter Ce2and the correlation coefficient betweene and T on the calculation of sensible heat flux using CT2 profiles measured by sodar. A new formula is suggested in terms of parameterization. The sensible heat flux calculated by sodar measurements is compared with that by bulk transfer method, and the results agree well.展开更多
Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolat...Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to moving heat source.展开更多
In this paper, the analytical solution of a viscous and incompressible fluid towards an exponentially stretching porous sheet with surface heat flux in porous medium, for the boundary layer and heat transfer flow, is ...In this paper, the analytical solution of a viscous and incompressible fluid towards an exponentially stretching porous sheet with surface heat flux in porous medium, for the boundary layer and heat transfer flow, is presented. The equations of continuity, momentum and the energy are transformed into non-linear ordinary differential by using similarity transformation. The solutions of these highly non-linear ordinary differential equations are found analytically by means of Homotopy Analysis Method (HAM). The result obtained by HAM is compared with numerical results presented in the literature. The accuracy of the HAM is indicated by close agreement of the two sets of results. By this method, an expression is obtained which is admissible for all values of effective parameters. This method has the ability to control the convergence of the solution.展开更多
A variational technique (VT) is applied to estimate surface sensible and latent heat fluxes based on observations of air temperature, wind speed, and humidity, respectively, at three heights (1 m, 4 m, and 10 m), and ...A variational technique (VT) is applied to estimate surface sensible and latent heat fluxes based on observations of air temperature, wind speed, and humidity, respectively, at three heights (1 m, 4 m, and 10 m), and the surface energy and radiation budgets by the surface energy and radiation system (SERBS). The method fully uses all information provided by the measurements of air temperature, wind, and humidity profiles, the surface energy budget, and the similarity profile formulae as well. Data collected at Feixi experiment station installed by the China Heavy Rain Experiment and Study (HeRES) Program are used to test the method. Results show that the proposed technique can overcome the well-known unstablility problem that occurs when the Bowen method becomes singular; in comparison with the profile method, it reduces both the sensitivities of latent heat fluxes to observational errors in humidity and those of sensible heat fluxes to observational errors in temperature, while the estimated heat fluxes approximately satisfy the surface energy budget. Therefore, the variational technique is more reliable and stable than the two conventional methods in estimating surface sensible and latent heat fluxes.展开更多
Three approaches, i.e., the harmonic analysis (HA) technique, the thermal diffusion equation and correction (TDEC) method, and the calorimetric method used to estimate ground heat flux, are evaluated by using obse...Three approaches, i.e., the harmonic analysis (HA) technique, the thermal diffusion equation and correction (TDEC) method, and the calorimetric method used to estimate ground heat flux, are evaluated by using observations from the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) in July, 2008. The calorimetric method, which involves soil heat flux measurement with an HFP01SC self-calibrating heat flux plate buried at a depth of 5 cm and heat storage in the soil between the plate and the surface, is here called the ITHP approach. The results show good linear relationships between the soil heat fluxes measured with the HFP01SC heat flux plate and those calculated with the HA technique and the TDEC method, respectively, at a depth of 5 cm. The soil heat fluxes calculated with the latter two methods well follow the phase measured with the HFP01SC heat flux plate. The magnitudes of the soil heat flux calculated with the HA technique and the TDEC method are close to each other, and they are about 2 percent and 6 percent larger than the measured soil heat flux, respectively, which mainly occur during the nighttime. Moreover, the ground heat fluxes calculated with the TDEC method and the HA technique are highly correlated with each other (R2= 0.97), and their difference is only about 1 percent. The TDEC-calculated ground heat flux also has a good linear relationship with the ITttP-calculated ground heat flux (R2 = 0.99), but their difference is larger (about 9 percent). Furthermore, compared to the HFP01SC direct measurements at a depth of 5 cm, the ground heat flux calculated with the HA technique, the TDEC method, and the ITHP approach can improve the surface energy budget closure by about 6 percent, 7 percent, and 6 percent at SACOL site, respectively. Therefore, the contribution of ground heat flux to the surface energy budget is very important for the semi-arid grassland over the Loess Plateau in China. Using turbulent heat fluxes with common corrections, soil heat storage between the surface and the heat flux plate can improve the surface energy budget closure by about 6 to 7 percent, resulting in a closure of 82 to 83 percent at the SACOL site.展开更多
Oceanic heat flux(Fw) is the vertical heat flux that is transmitted to the base of sea ice. It is the main source of sea ice bottom melting. The residual method was adopted to study oceanic heat flux under sea ice. Th...Oceanic heat flux(Fw) is the vertical heat flux that is transmitted to the base of sea ice. It is the main source of sea ice bottom melting. The residual method was adopted to study oceanic heat flux under sea ice. The data acquired by 28 ice mass balance buoys(IMBs) deployed over the period of 2004 to 2013 in the Arctic Ocean were used. Fw values presented striking seasonal and spatial variations. The average summer Fw values for the Canada Basin, Transpolar Drift, and Multiyear Ice area were 16.8, 7.7, and 5.9 W m^-2, respectively. The mean summer F-w for the whole Arctic was 10.1 W m^-2, which was equivalent to a bottom melt of 0.4 m. Fw showed an autumn peak in November in the presence of the near-surface temperature maximum(NSTM). The average Fw for October to December was 3.7 W m^-2. And the average Fw for January to March was 1.0 W m^-2, which was approximately one third of the average Fw in the presence of NSTM. The summer Fw was almost wholly attributed to the incident solar radiation that enters the upper ocean through leads and the open water. Fw calculated through the residual method using IMB data was compared with that calculated through the parameterization method using Autonomous Ocean Flux Buoy data. The results revealed that the Fw provided by the two methods were consistent when the sea ice concentration exceeded 70% and mixing layer temperature departure from freezing point was less than 0.15℃. Otherwise, the Fw yielded by the residual method was approximately one third smaller than that provided by the parameterization method.展开更多
Research on heat and mass flux yielded by modern seafloor hydrothermal activity is very important, because it is involved not only in the base of ocean environment research, but also in the historical evolution of sea...Research on heat and mass flux yielded by modern seafloor hydrothermal activity is very important, because it is involved not only in the base of ocean environment research, but also in the historical evolution of seawater properties. Currently, estimating heat flux is based on the observation data of hydrothermal smokers, low-temperature diffusive flow and mid-ocean ridge mainly. But there are some faults, for example, there is lack of a concurrent conductive item in estimating the heat flux by smokers and the error between the half-space cooling model and the observation data is too large. So, three kinds of methods are applied to re-estimating the heat flux of hydrothermal activity resepectively, corresponding estimation is 97. 359 GW by hydrothermal smoker and diffusive flow, 84.895 GW by hydrothermal plume, and 4. 11 TW by exponential attenuation method put forward by this paper. Research on mass flux estimation is relatively rare, the main reason for this is insufficient field observation data. Mass fluxes of different elements are calculated using hydrothermal vent fluid data from the TAG hydrothermal area on the Mid-Atlantic Ridge for the first time. Difference of estimations by different methods reflects the researching extent of hydrothermal activity, and systematically in - situ observation will help to estimate the contribution of hydrothermal activity to ocean chemical environment, ocean circulation and global climate precisely.展开更多
The connecting rod is one of the most important moving components in an internal combustion engine. The present work determined the possibility of using aluminium alloy 7075 material to design and manufacture a connec...The connecting rod is one of the most important moving components in an internal combustion engine. The present work determined the possibility of using aluminium alloy 7075 material to design and manufacture a connecting rod for weight optimisation without losing the strength of the connecting rod. It considered modal and thermal analyses to investigate the suitability of the material for connecting rod design. The parameters that were considered under the modal analysis were: total deformation, and natural frequency, while the thermal analysis looked at the temperature distribution, total heat flux and directional heat flux of the four connecting rods made with titanium alloy, grey cast iron, structural steel and aluminium 7075 alloy respectively. The connecting rod was modelled using Autodesk inventor2017 software using the calculated parameters. The steady-state thermal analysis was used to determine the induced heat flux and directional heat flux. The study found that Aluminium 7075 alloy deformed more than the remaining three other materials but has superior qualities in terms of vibrational natural frequency, total heat flux and lightweight compared to structural steel, grey cast iron and titanium alloy.展开更多
The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the bea...The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.展开更多
The calculative method of heat transfer coefficient between ice cover and water is analyzed considering the heat balance at ice cover bottom firstly. The heat transfer coefficient is calculated with the meteorological...The calculative method of heat transfer coefficient between ice cover and water is analyzed considering the heat balance at ice cover bottom firstly. The heat transfer coefficient is calculated with the meteorological, oceanographic data and sea ice conditions measured on the JZ20-2 Oil/Gas Platform in the Bohai Sea during the winter of 1997/1998. From the results, it is shown that the heat transfer coefficient is smaller in the freezing and melting periods, which is about 0.16×10-3 and 0.04×10-3 respectively. In the middle of ice season, the heat transfer coefficient has a larger value, which is about 0.5×10-3. Lastly, the influences of ice thickness and ice type on the heat transfer coefficient are discussed. With the heat transfer coefficient determined above, the oceanic heat flux in the winter of 1997~1998 is calculated, and its trend in the winter is analyzed. This study can be referenced in the sea ice numerical simulation and prediction in the Bohai Sea.展开更多
基金financially supported by the National Natural Science Foundation of China (No.50464004)
文摘A new hybrid method, Monte-Carlo-Heat-Flux (MCHF) method, was presented to analyze the radiative heat transfer of participating medium in a three-dimensional rectangular enclosure using combined the Monte-Carlo method with the heat flux method. Its accuracy and reliability was proved by comparing the computational results with exact results from classical "Zone Method".
文摘This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.
文摘During the second course of USA - PRC joint air sea interaction experiment in 1986, the temperature structure parameters CT2 were measured by sodar over the Western Pacific Ocean. Based on similarity theory, a method is discussed to calculate the sensible heat flux over the ocean in unstable stratification. Becausehumidity is great over the ocean, so we have to consider the influence of water vapor structure parameter Ce2and the correlation coefficient betweene and T on the calculation of sensible heat flux using CT2 profiles measured by sodar. A new formula is suggested in terms of parameterization. The sensible heat flux calculated by sodar measurements is compared with that by bulk transfer method, and the results agree well.
文摘Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to moving heat source.
文摘In this paper, the analytical solution of a viscous and incompressible fluid towards an exponentially stretching porous sheet with surface heat flux in porous medium, for the boundary layer and heat transfer flow, is presented. The equations of continuity, momentum and the energy are transformed into non-linear ordinary differential by using similarity transformation. The solutions of these highly non-linear ordinary differential equations are found analytically by means of Homotopy Analysis Method (HAM). The result obtained by HAM is compared with numerical results presented in the literature. The accuracy of the HAM is indicated by close agreement of the two sets of results. By this method, an expression is obtained which is admissible for all values of effective parameters. This method has the ability to control the convergence of the solution.
基金This work was supported by the National Natural Science Foundation of China under Grant No.E-D0119-90202014the National Key Programme for Developing Basic Sciences of China under Great No.G1998040902.
文摘A variational technique (VT) is applied to estimate surface sensible and latent heat fluxes based on observations of air temperature, wind speed, and humidity, respectively, at three heights (1 m, 4 m, and 10 m), and the surface energy and radiation budgets by the surface energy and radiation system (SERBS). The method fully uses all information provided by the measurements of air temperature, wind, and humidity profiles, the surface energy budget, and the similarity profile formulae as well. Data collected at Feixi experiment station installed by the China Heavy Rain Experiment and Study (HeRES) Program are used to test the method. Results show that the proposed technique can overcome the well-known unstablility problem that occurs when the Bowen method becomes singular; in comparison with the profile method, it reduces both the sensitivities of latent heat fluxes to observational errors in humidity and those of sensible heat fluxes to observational errors in temperature, while the estimated heat fluxes approximately satisfy the surface energy budget. Therefore, the variational technique is more reliable and stable than the two conventional methods in estimating surface sensible and latent heat fluxes.
基金supported by the National Natural Science Foundation of China (GrantNo. 40725015)
文摘Three approaches, i.e., the harmonic analysis (HA) technique, the thermal diffusion equation and correction (TDEC) method, and the calorimetric method used to estimate ground heat flux, are evaluated by using observations from the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) in July, 2008. The calorimetric method, which involves soil heat flux measurement with an HFP01SC self-calibrating heat flux plate buried at a depth of 5 cm and heat storage in the soil between the plate and the surface, is here called the ITHP approach. The results show good linear relationships between the soil heat fluxes measured with the HFP01SC heat flux plate and those calculated with the HA technique and the TDEC method, respectively, at a depth of 5 cm. The soil heat fluxes calculated with the latter two methods well follow the phase measured with the HFP01SC heat flux plate. The magnitudes of the soil heat flux calculated with the HA technique and the TDEC method are close to each other, and they are about 2 percent and 6 percent larger than the measured soil heat flux, respectively, which mainly occur during the nighttime. Moreover, the ground heat fluxes calculated with the TDEC method and the HA technique are highly correlated with each other (R2= 0.97), and their difference is only about 1 percent. The TDEC-calculated ground heat flux also has a good linear relationship with the ITttP-calculated ground heat flux (R2 = 0.99), but their difference is larger (about 9 percent). Furthermore, compared to the HFP01SC direct measurements at a depth of 5 cm, the ground heat flux calculated with the HA technique, the TDEC method, and the ITHP approach can improve the surface energy budget closure by about 6 percent, 7 percent, and 6 percent at SACOL site, respectively. Therefore, the contribution of ground heat flux to the surface energy budget is very important for the semi-arid grassland over the Loess Plateau in China. Using turbulent heat fluxes with common corrections, soil heat storage between the surface and the heat flux plate can improve the surface energy budget closure by about 6 to 7 percent, resulting in a closure of 82 to 83 percent at the SACOL site.
基金supported by the National Basic Research Program of China (No. 2015CB953900)the National Natural Science Foundation of China (No. 413309 60)the National Key Research and Development Plan (No. 2016YFC1402705)
文摘Oceanic heat flux(Fw) is the vertical heat flux that is transmitted to the base of sea ice. It is the main source of sea ice bottom melting. The residual method was adopted to study oceanic heat flux under sea ice. The data acquired by 28 ice mass balance buoys(IMBs) deployed over the period of 2004 to 2013 in the Arctic Ocean were used. Fw values presented striking seasonal and spatial variations. The average summer Fw values for the Canada Basin, Transpolar Drift, and Multiyear Ice area were 16.8, 7.7, and 5.9 W m^-2, respectively. The mean summer F-w for the whole Arctic was 10.1 W m^-2, which was equivalent to a bottom melt of 0.4 m. Fw showed an autumn peak in November in the presence of the near-surface temperature maximum(NSTM). The average Fw for October to December was 3.7 W m^-2. And the average Fw for January to March was 1.0 W m^-2, which was approximately one third of the average Fw in the presence of NSTM. The summer Fw was almost wholly attributed to the incident solar radiation that enters the upper ocean through leads and the open water. Fw calculated through the residual method using IMB data was compared with that calculated through the parameterization method using Autonomous Ocean Flux Buoy data. The results revealed that the Fw provided by the two methods were consistent when the sea ice concentration exceeded 70% and mixing layer temperature departure from freezing point was less than 0.15℃. Otherwise, the Fw yielded by the residual method was approximately one third smaller than that provided by the parameterization method.
基金This study was supported by the Major State Basic Research Program of China under contract No.G2000078503the National Natural Science Foundation of China under contract No.40246024.
文摘Research on heat and mass flux yielded by modern seafloor hydrothermal activity is very important, because it is involved not only in the base of ocean environment research, but also in the historical evolution of seawater properties. Currently, estimating heat flux is based on the observation data of hydrothermal smokers, low-temperature diffusive flow and mid-ocean ridge mainly. But there are some faults, for example, there is lack of a concurrent conductive item in estimating the heat flux by smokers and the error between the half-space cooling model and the observation data is too large. So, three kinds of methods are applied to re-estimating the heat flux of hydrothermal activity resepectively, corresponding estimation is 97. 359 GW by hydrothermal smoker and diffusive flow, 84.895 GW by hydrothermal plume, and 4. 11 TW by exponential attenuation method put forward by this paper. Research on mass flux estimation is relatively rare, the main reason for this is insufficient field observation data. Mass fluxes of different elements are calculated using hydrothermal vent fluid data from the TAG hydrothermal area on the Mid-Atlantic Ridge for the first time. Difference of estimations by different methods reflects the researching extent of hydrothermal activity, and systematically in - situ observation will help to estimate the contribution of hydrothermal activity to ocean chemical environment, ocean circulation and global climate precisely.
文摘The connecting rod is one of the most important moving components in an internal combustion engine. The present work determined the possibility of using aluminium alloy 7075 material to design and manufacture a connecting rod for weight optimisation without losing the strength of the connecting rod. It considered modal and thermal analyses to investigate the suitability of the material for connecting rod design. The parameters that were considered under the modal analysis were: total deformation, and natural frequency, while the thermal analysis looked at the temperature distribution, total heat flux and directional heat flux of the four connecting rods made with titanium alloy, grey cast iron, structural steel and aluminium 7075 alloy respectively. The connecting rod was modelled using Autodesk inventor2017 software using the calculated parameters. The steady-state thermal analysis was used to determine the induced heat flux and directional heat flux. The study found that Aluminium 7075 alloy deformed more than the remaining three other materials but has superior qualities in terms of vibrational natural frequency, total heat flux and lightweight compared to structural steel, grey cast iron and titanium alloy.
基金supported by the National Natural Science Foundation of China(No.11605234)the Foundation of ASIPP(No.DSJJ-15-GC02)
文摘The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.
基金co-sponsored by the National Foundation of Natural Sciencein part by the Bohai Offshore Oil Corporation(Grant No.59739170).
文摘The calculative method of heat transfer coefficient between ice cover and water is analyzed considering the heat balance at ice cover bottom firstly. The heat transfer coefficient is calculated with the meteorological, oceanographic data and sea ice conditions measured on the JZ20-2 Oil/Gas Platform in the Bohai Sea during the winter of 1997/1998. From the results, it is shown that the heat transfer coefficient is smaller in the freezing and melting periods, which is about 0.16×10-3 and 0.04×10-3 respectively. In the middle of ice season, the heat transfer coefficient has a larger value, which is about 0.5×10-3. Lastly, the influences of ice thickness and ice type on the heat transfer coefficient are discussed. With the heat transfer coefficient determined above, the oceanic heat flux in the winter of 1997~1998 is calculated, and its trend in the winter is analyzed. This study can be referenced in the sea ice numerical simulation and prediction in the Bohai Sea.