METHOD FOR SIMULATING PERFORMANCE OF LNG AMBIENT AIR VAPORIZER UNDER FROSTING CONDITION

Abstract

The present disclosure discloses a method for simulating performance and optimizing design of an LNG ambient air vaporizer under a frosting condition, performing site operation test on the LNG ambient air vaporizer to obtain fitting relationship between the outer wall temperature of frosted finned tube and frost layer; transforming a sum of increased frost layer thermal resistance and thermal resistance into equivalent thermal contact resistance, and representing the equivalent thermal contact resistance as a function of the outer wall temperature; establishing a calculation model of the LNG ambient air vaporizer, performing simulation calculation by equivalent heat conduction coefficient, so as to obtain fluid-solid conjugate heat transfer characteristics and vaporization performance of the LNG ambient air vaporizer under the frosting condition for optimization design of the LNG ambient air vaporizer.

Description

TECHNICAL FIELD

The present disclosure belongs to methods for simulating performance and optimizing design of low-temperature heat exchange equipment, and particularly relates to a method for simulating performance and optimizing design of an LNG ambient air vaporizer under a frosting condition.

BACKGROUND

An LNG (liquefied natural gas) ambient air vaporizer is low-temperature heat exchange equipment that absorbs heat from air and transmits the heat to low-temperature LNG in tubes so as to vaporize the LNG into NG (natural gas) through phase change. Compared with vaporizers of other types, the LNG ambient air vaporizer is simple in structure and good in economical efficiency, meets the requirements of energy conservation and environmental protection, and is widely applied to LNG receiving stations and vaporizing stations. However, outer surfaces of finned tubes of the LNG ambient air vaporizer are easily affected by ambient air temperature and humidity, surrounding ambient temperature is reduced due to the absorption of surrounding heat during vaporization, and surfaces of the finned tubes frost, which affects the heat exchange effect and vaporization performance of the vaporizer; and in a severe case, the vaporizer may be non-uniformly stressed, and lateral tension is generated, resulting in pipeline breakage and LNG leakage, which leads to safety accidents.

The reason for frosting of the surfaces of the finned tubes of the LNG ambient air vaporizer is mainly that water vapor in air is formed into solid crystals through phase change, and geometric structures and stacking rules of the solid crystals change with the progress of frosting, which is a complex heat and mass transfer process; in the past performance simulation for the LNG ambient air vaporizer, the effect of frosting on the vaporizer was generally ignored, and ambient air was simplified as dry air ideal gas; and during actual operation of the vaporizer, frosting has a non-negligible effect on performance of the LNG ambient air vaporizer, which will directly reduce a heat transfer coefficient of an air side by up to 85%. Currently, research on the frosting performance of the LNG ambient air vaporizer falls into two main categories. One approach focuses on unstable-state frosting numerical simulation analysis for a local region of a single fin. Although some models can accurately predict the growth situation of a frost layer, the models cannot be applied to simulation of the whole vaporizer considering the complexity of these models, calculation duration and convergence of calculation results; and the other approach is to establish frosting heat transfer mathematical models to determine the influence of the frost layer on the total heat transfer coefficient of the ambient air vaporizer over different frosting times through numerical calculation. However, this method only calculates the overall performance change of the ambient air vaporizer under a frosting condition and neglects interactions between different finned tubes. Consequently, it cannot analyze the heat transfer performance of the different finned tubes or low-temperature media inside the tubes of the vaporizer. CN112580272A has provided an optimal design method of an LNG ambient air vaporizer based on numerical simulation. The method includes the steps: firstly, primarily designing the LNG ambient air vaporizer by adopting methods such as an empirical formula, and then simulating the primarily designed ambient air vaporizer by adopting Fluent software. Although the influence of the number of finned tubes on heat exchange efficiency of the vaporizer is considered in the method, the method cannot calculate the performance of the vaporizer under a frosting condition when treating ambient air as dry air. CN115114815A has provided a simulation method for predicting cold surface frosting by using surface properties of frost layers. The method includes the steps: establishing a calculation model for a condition of frosting over a flat plate through Eulerian multi-phase, performing repeated iterative calculations to solve a control equation for a calculation domain through user-defined functions (UDF) to simulate a growth of the frost layers. However, the method can only simulate frosting degrees of the frost layers in local regions of the flat plate, but cannot obtain the performance change of overall heat exchange equipment in a frosting state. In addition, frosting mainly occurs in a liquid phase section and a two-phase section, the influence of fluid-solid conjugate heat transfer needs to be considered during simulation calculation in the frosting state, however, there is no technical method to explain it.

SUMMARY

The present disclosure aims at providing a method for simulating performance and optimizing design of an LNG ambient air vaporizer under a frosting condition to overcome the defects in the prior art. Interactions between different finned tubes and the influence of frosting of the finned tubes are considered, heat transfer in gas-liquid phase change flow and fluid-solid conjugate heat transfer in the tubes of the vaporizer are coupled, as a result, models are more realistic, and simulation calculation results are more accurate; and a transient process of frosting of the finned tubes is simplified into a quasi-stable state process, and the method can save a large amount of calculation duration and is high in calculation efficiency and reliability.

A method for simulating performance and optimizing design of an LNG ambient air vaporizer under a frosting condition, includes the following steps, executed by a processor:

step 1, performing site operation test on the LNG ambient air vaporizer, measuring and obtaining site ambient temperature T₀, humidity H₀ and atmospheric pressure P₀, and measuring and obtaining, at a certain operating moment t_n, pressure P_in of liquefied natural gas at a vaporizer inlet and pressure P_out of natural gas at a vaporizer outlet, flow velocity V_in of the liquefied natural gas at the vaporizer inlet and temperature T_in of the liquefied natural gas at the vaporizer inlet, outer wall temperature T_s at different positions of each frosted finned tube of the vaporizer, frost layer temperature T_f, frost layer thickness d_f, and air flow velocity V_a outside a frost layer on a surface of the finned tube, the different positions referring to at least three equidistant point positions of an outer wall of a fin of each frosted finned tube from top to bottom; wherein, the site ambient temperature T₀, humidity H₀ and atmospheric pressure P₀ are measured by a temperature and humidity and atmospheric pressure instrument; pressure P_in and pressure P_out are measured by a pressure gauge disposed at the inlet and outlet of the air temperature gasifier, respectively; the flow velocity V_in is measured by a flow meter provided at the inlet of the air temperature gasifier; the temperature T_in and outer wall temperature T_s are measured by temperature sensors, which are arranged at an inlet of the gasifier and a different location of each finned tube of the gasifier; the frost layer temperature T_f is measured by temperature gauges; the frost layer thickness d_f is measured by vernier calipers; the air flow velocity V_a is measured by anemometers.

step 2, performing, by data processing software, fitting analysis on the collected frost layer temperature T_f and frost layer thickness d_f at the different positions of all the frosted finned tubes at the certain operating moment t_n and the air flow velocity V_a outside the frost layer on the surface of the finned tube, and data of the outer wall temperature T_s of the finned tube respectively, and obtaining, by a least square method, fitting relational expressions d_f=f_d(T_s)=A₁T_s²+B₁T_s+C₁, T_f=f_T(T_s)=A₂T_s+B₂, and V_a=f_V(T_s)=A₃T_s²+B₃T_s+C₃ between the outer wall temperature T_s of all the frosted finned tubes of the whole vaporizer and the frost layer thickness d_f, the frost layer temperature T_f and the air flow velocity V_a outside the frost layer on the surface of the finned tube at the certain operating moment t_n, A₁, B₁, C₁, A₂, B₂, A₃, B₃, and C₃ in the formula being respectively fitted constants;

step 3, establishing a calculation model for an equivalent thermal conductivity coefficient of the LNG ambient air vaporizer at the certain operating moment t_n under frost conditions: transforming a sum of increased frost layer thermal resistance R_f of each frosted finned tube in unit length at the certain operating moment t_n and thermal resistance R_o of the finned tube body in unit length into equivalent thermal contact resistance R_e of the finned tube in unit length under a non-frosting condition, and representing an equivalent thermal conductivity coefficient λ_e of the equivalent thermal contact resistance R_e as a function of the outer wall temperature T_s of the frosted finned tube in unit length, the unit length being a length of a minimum mesh during geometric meshing of the vaporizer; wherein

${\lambda }_e=F\left(T_s\right)=\frac{1}{\frac{1}{\lambda }+\frac{1}{d_{in}\ln \frac{d_{out}}{d_{in}}}·\frac{Z\left(T_s\right)}{\frac{A_2^{\prime }\beta }{A_2}+\frac{A_2^{″}\beta }{A_2}·\frac{\mathrm{th}\left(A_m\sqrt{f_V\left(T_s\right)}\right)}{A_m\sqrt{f_V\left(T_s\right)}}}}$

wherein, parameters λ, d_in, d_out, A₂, A₂′, A₂″, A_m, and β in the formula are all constant values; λ is a thermal conductivity coefficient of a vaporizer material, namely aluminum alloy; d_in is an internal diameter of the finned tube, and d_out is an external diameter of the finned tube; A₂ is an external surface area of the finned tube in unit length; A₂′ is a surface area of a non-fin part outside the finned tube in unit length; A₂″ is a surface area of a fin part outside the finned tube in unit length; A_m=l*[36/(δ*80 )]^1/2, where l is a fin height, and δ is a thickness of the fin; β is a finning coefficient of the finned tube, β=A₀/A₂, A₀ is an internal surface area of the finned tube, and A₂ is an external surface area of the finned tube; f_V(T_s) is a function expression relational expression of the air flow velocity V_a outside the frost layer, V_a=f_V(T_s), and T_s is the outer wall temperature at different positions of each frosted finned tube of the vaporizer; and Z(T_s) is a function expression relational expression of the frost layer thermal resistance R_f, R_f=Z(T_s)=d_f/λ_f=f_d(T_s)/g(T_s), where, g(T_s) is a function expression relational expression of the frost layer thermal conductivity coefficient λ_f, λ_f=g(T_s)=0.001202×(650e^{0.277[(fT(Ts)−273.15)]})^0.963;

step 4, establishing an overall geometric model of the LNG ambient air vaporizer in simulation software, performing meshing and dividing of computational domain, selecting physical models and equations, setting material attributes and boundary conditions of the computational domain, adopting the equivalent thermal conductivity coefficient λ_e of the finned tube of the vaporizer as the thermal conductivity coefficient during frosting of the vaporizer material, performing solving and initialized setting, and then performing simulation calculation, specifically as follows:

S1: establishing, by three-dimensional geometric modeling software, the overall geometric model of the LNG ambient air vaporizer in a ratio of 1:1, and then performing, by finite element meshing software, meshing and dividing of the computational domain of the overall geometric model; dividing the computational domain into an LNG fluid domain, a vaporizer solid domain and an air fluid domain; the LNG fluid domain is a flow region of LNG in an internal channel of the vaporizer; the vaporizer solid domain is a vaporizer body; and the air fluid domain is an air flow region outside the vaporizer body;

S2: importing the meshed overall geometric model of the LNG ambient air vaporizer into fluid analysis software, and adopting the LNG fluid domain, the vaporizer solid domain and the air fluid domain as the computational domain; setting a contact surface between the LNG fluid domain and the vaporizer solid domain and a contact surface between the vaporizer solid domain and the air fluid domain as Interface surfaces, and checking a Couple option in Interface setting, so that the corresponding contact surface can complete heat transfer;

S3: enabling a gravity model, a multi-phase model, a turbulence model, a boiling phase change model, a continuity equation, a momentum equation, an energy equation and a component transport equation in the fluid analysis software, and adopting a standard wall function method for near-wall processing; adopting a Mixture model as the multi-phase model; adopting a Realizable k-ϵ turbulence model as the turbulence model, and adopting an evaporation-condensation Lee model as the boiling phase change model;

S4: setting the material attributes of the computational domain:

respectively introducing LNG and NG fluid materials in the fluid analysis software, then setting the LNG fluid material as a first term in the multi-phase model, setting phase change from LNG to NG, and selecting the evaporation-condensation Lee model for a reaction mechanism; introducing an aluminum alloy solid material in the fluid analysis software, adopting physical data in a software material library as parameters of the aluminum alloy solid material, then modifying the thermal conductivity coefficient of the aluminum alloy solid material from the constant value λ to a value represented by a piecewise polynomial temperature function method, setting the thermal conductivity coefficient of the vaporizer material within a frosting temperature range as the equivalent thermal conductivity coefficient λ_e=F(T_s), and setting the thermal conductivity coefficient of the vaporizer material within a non-frosting temperature range as the constant value λ; introducing a wet air mixed material in the fluid analysis software, the wet air mixed material including air and water vapor, and physical property data in the software material library are adopted as material attributes of the wet air mixed material;

S5: setting the boundary conditions of the computational domain:

setting an outlet of the LNG fluid domain as a pressure outlet boundary, and adopting the pressure P_out at the vaporizer outlet tested on site as pressure; setting an inlet of the LNG fluid domain as a velocity inlet boundary, and adopting the flow velocity V_in and temperature T_in at the vaporizer inlet tested on site as velocity and temperature; setting a top surface of the air fluid domain above the vaporizer and a side surface of the air fluid domain around the vaporizer as pressure inlet boundaries, adopting the atmospheric pressure P₀ and ambient temperature T₀ tested on site as pressure and temperature, and setting air humidity of the air fluid domain for simulating the air humidity according to the ambient humidity tested on site; setting a bottom surface of the air fluid domain at a bottom of the vaporizer as a pressure outlet boundary;

defining the air fluid domain in this step as a hexahedron capable of surrounding the vaporizer, space, except for the vaporizer, inside the hexahedron representing air outside the vaporizer, and the top surface, the side surface and the bottom surface of the air fluid domain referring to a top surface, a side surface and a bottom surface of the hexahedron outside the whole vaporizer;

S6: performing initialization setting by adopting a SIMPLE algorithm in the fluid analysis software as a solving method for the geometric model meshes divided in step S1, and then performing calculation simulation on the geometric model established in step S1; stopping calculation and outputting result data from numerical simulation in heat transfer of the vaporizer if a residual variance curve converges and the monitored NG outlet temperature and flow velocity do not change any more; or else, continuing to operate; and

step 5, importing the result data from numerical simulation into post-processing software for analysis, displaying a temperature cloud diagram on a surface of the LNG ambient air vaporizer, and a temperature cloud diagram, a velocity cloud diagram and a component cloud diagram of the LNG fluid domain in the tube of the vaporizer in the post-processing software, selecting an outlet section of the LNG fluid domain, to obtain LNG outlet temperature and outlet flow velocity, and selecting the component cloud diagram of the LNG fluid domain, to obtain proportions of a liquid phase section, a two-phase section and a gas-phase section in the fluid domain; and checking temperature, heat flux and other thermal parameters at different position points or sections of the surface of the vaporizer, so as to visually acquire fluid-solid conjugate heat transfer characteristics and vaporization performance of the LNG ambient air vaporizer under the frosting condition.

The present disclosure has the beneficial effects:

1. According to the method for simulating and optimizing the performance of the LNG ambient air vaporizer under the frosting condition provided by the present disclosure, the transient process of frosting of the finned tubes is simplified into the quasi-stable-state process, that is, it is considered that the process is stable in a certain operating time step length; the sum of the increased frost layer thermal resistance at different operating moments under the frosting condition and the thermal resistance of the finned tube body is transformed into the equivalent thermal contact resistance of the materials of the finned tubes of the vaporizer under the non-frosting condition, and the thermal conductivity coefficient of the equivalent thermal contact resistance of the finned tubes is represented as the function of the outer wall temperature of the finned tubes; and by means of the method, simulation and optimization for overall performance of the vaporizer in the frosting state is achieved, the interactions between the different finned tubes are considered, the simulation and optimization results are more accurate, the simulation and optimization calculation duration can be greatly shortened, and the simulation and optimization efficiency is improved.

2. According to the method for simulating and optimizing the performance of the LNG ambient air vaporizer under the frosting condition provided by the present disclosure, the influences of heat transfer in gas-liquid phase change flow and fluid-solid conjugate heat transfer in the tubes of the vaporizer on the performance of the vaporizer under the frosting condition are considered, the models are more realistic, the simulation and optimization calculation results are more accurate, and the method can be used for simulation and optimization research on the heat transfer performance and the vaporization performance of the vaporizer under the frosting condition and design optimization for the vaporizer, and has good practicality.

BRIEF DESCRIPTION OF THE DRAWINGS

Drawings of the specification constituting a part of the present disclosure are described for further understanding the present disclosure. Schematic embodiments of the present disclosure and descriptions thereof are schematic of the present disclosure, and are not construed as an improper limitation to the present disclosure.

FIG. 1 is a flowchart of a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure;

FIG. 2 is a quarter cross-section diagram of a finned tube of an LNG ambient air vaporizer used for performance simulation under a frosting condition according to the present disclosure;

FIG. 3 is a principle diagram of transforming a sum of frost layer thermal resistance R_f and thermal resistance R_o of a finned tube body into an equivalent thermal contact resistance R_e of the finned tube adopting a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure;

FIG. 4 is a schematic diagram of a relationship between frost layer thermal resistance R_f at different positions of a finned tube and outer wall temperature T_s of the finned tube adopting a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure;

FIG. 5 is a diagram of a fitting curve of frost layer thermal resistance and outer wall temperature of a finned tube at different operating times adopting a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure;

FIG. 6 is a schematic diagram of a geometric model of a vaporizer adopting a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure;

FIG. 7 is a top view of the geometric model of the vaporizer shown in FIG. 6;

FIG. 8 is a front view of the geometric model of the vaporizer shown in FIG. 6;

FIG. 9 is a side view of the geometric model of the vaporizer shown in FIG. 6;

FIG. 10 is a schematic diagram of meshing of a geometric model of a vaporizer adopting a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure;

FIG. 11 is a proportion comparison diagram of liquid phase sections, two-phase sections and gas liquid sections in branch tubes of different vaporizer finned tubes in a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure; and

FIG. 12 is a comparison diagram of simulation results versus actual measurement results of vaporizer outlet temperature at different operating times in a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure.

In the drawings, 1: frost layer outside fin; 2: finned tube of vaporizer; and 3: LNG.

DETAILED DESCRIPTION OF THE PRESENT DISCLOSURE

It should be noted that the following detailed descriptions are exemplary, which are intended to further explain the present application. Unless otherwise indicated, all technical and scientific terms used herein have the same meaning as commonly understood by those ordinarily skilled in the prior art to which the present application pertains.

It should be noted that the terms used here are not intended to limit the exemplary implementations according to the present application, but are merely descriptive of the specific implementations. Unless otherwise directed by the context, singular forms of terms used here are intended to include plural forms. Besides, it should be also appreciated that, when the terms “comprise” and/or “include” are used in the specification, it is indicated that characteristics, steps, operations, devices, assemblies, and/or combinations thereof exist.

Additionally, any directional indication (such as upper, lower, left, right, front, back, or the like) involved in the embodiments of the present disclosure is only used for explaining relative position relations, movement conditions and the like of components in a certain specific posture (as shown in figures). If the specific posture changes, the directional indications may change accordingly.

As shown in FIG. 1, a method for simulating performance of an LNG ambient air vaporizer under a frosting condition according to the present disclosure, includes the following steps:

step 1, site operation test is performed on the LNG ambient air vaporizer, site ambient temperature T₀, humidity H₀ and atmospheric pressure P₀ are measured, and pressure P_in of liquefied natural gas at a vaporizer inlet and pressure P_out of natural gas at a vaporizer outlet, flow velocity V_in of the liquefied natural gas at the vaporizer inlet and temperature T_in of the liquefied natural gas at the vaporizer inlet, outer wall temperature T_s at different positions of each frosted finned tube of the vaporizer, frost layer temperature T_f, frost layer thickness d_f, and air flow velocity V_a outside a frost layer on a surface of the finned tube at a certain operating moment t_n (such as 1 h, 2 h, 4 h or 8 h) are measured, the different positions referring to at least three equidistant point positions of an outer wall of a fin of each frosted finned tube from top to bottom;

step 2, fitting analysis is performed on the collected frost layer temperature T_f and frost layer thickness d_f at the different positions of all the frosted finned tubes at the certain operating moment t_n and the air flow velocity V_a on the outer side of the frost layer on the surface of the finned tube, and data of the outer wall temperature T_s of the finned tube respectively is performed by data processing software (such as SPSS or Origin), and fitting relational expressions d_f=f_d(T_s)=A₁T_s²+B₁T_s+C₁, T_f=f_T(T_s)=A₂T_s+B₂, and V_a=f_V(T_s)=A₃T_s²+B₃T_s+C₃ between the outer wall temperature T_s of all the frosted finned tubes of the whole vaporizer and the frost layer thickness d_f, the frost layer temperature T_f and the air flow velocity V_a outside the frost layer on the surface of the finned tube at the certain operating moment t_n are obtained by a least square method, A₁, B₁, C₁, A₂, B₂, A₃, B₃, and C₃ in the formulas being respectively fitted constants.

In this step, a relationship between the outer wall temperature of each frosted finned tube and the frost layer thickness are integrated to fit a relational expression, that is, data of all the finned tubes at the same operating time t_n is collected, the data of all the frosted finned tubes of the vaporizer is fitted, and each relational expression (all the finned tubes are fitted in this relational expression) of d_f=f_d(T_s); T_f=f_T(T_s); V_a=f_V(T_s) is fitted; and different relational expressions may correspond to different operating times. Since temperature distribution of each finned tube is different, a temperature distribution range can be widened by collecting the data of all the finned tubes, which is conducive to improving fitting accuracy. FIG. 2 is a quarter cross-section diagram of a finned tube of an LNG ambient air vaporizer used for performance simulation under a frosting condition according to the present disclosure, which shows the frost layer thickness d_f, an external diameter d_out of the finned tube, an internal diameter d_in of the finned tube, and a thickness δ of a fin.

Step 3, a calculation model for an equivalent thermal conductivity coefficient, of the LNG ambient air vaporizer during frosting at the certain operating moment t_n is established: as shown in FIG. 3, a sum of increased frost layer thermal resistance R_f of each frosted finned tube in unit length (the unit length is a length of a minimum mesh during geometric meshing of the vaporizer) at the certain operating moment t_n and thermal resistance R_o of the finned tube body in unit length is transformed into equivalent thermal contact resistance R_e of the finned tube in unit length under a non-frosting condition, and an equivalent thermal conductivity coefficient λ_e of the equivalent thermal contact resistance R_e is represented as a function of the outer wall temperature T_s of the frosted finned tube in unit length:

${\lambda }_e=F\left(T_s\right)=\frac{1}{\frac{1}{\lambda }+\frac{1}{d_{in}\ln \frac{d_{out}}{d_{in}}}·\frac{Z\left(T_s\right)}{\frac{A_2^{\prime }\beta }{A_2}+\frac{A_2^{″}\beta }{A_2}·\frac{\mathrm{th}\left(A_m\sqrt{f_V\left(T_s\right)}\right)}{A_m\sqrt{f_V\left(T_s\right)}}}}.$

wherein, parameters λ, d_in, d_out, A₂, A₂′, A₂″, A_m, and β in the formula are all constant values; λ is a thermal conductivity coefficient of a vaporizer material, namely aluminum alloy, W/(m·K), which can be found in professional books (such as Practical Handbook of Nonferrous Metal Materials, Guangdong Science and Technology Press, 2006); d_in is the internal diameter of the finned tube, and d_out is the external diameter of the finned tube, m, which can be obtained through a design drawing; A₂ is an external surface area of the finned tube in unit length, m²; A₂′ is a surface area of a non-fin part outside the finned tube in unit length, m²; A₂″ is a surface area of a fin part outside the finned tube in unit length, m²; A₂, A₂′, and A₂″ can be calculated from the design drawing; A_m=l*[36/(λ*δ)]^1/2, where l is a fin height, δ is the thickness of the fin, and m, l, and δ can be obtained from the design drawing; β is a finning coefficient of the finned tube, β=A₀/A₂, A₀ is an internal surface area of the finned tube, and A₂ is an external surface area of the finned tube, m², which can be calculated from the design drawing. f_V(T_s) is a function expression relational expression of the air flow velocity Va outside the frost layer, and Ts is the outer wall temperature K at different positions of each frosted finned tube of the vaporizer; and Z(T_s) is a function expression relational expression of the frost layer thermal resistance R_f, R_f=Z(T_s)=d_f/λ_f=f_d(T_s)/g(T_s), where, g(T_s) is a function expression relational expression of the frost layer thermal conductivity coefficient λ_f, λ_f=g(T_s)=0.001202×(650 e^{0.277[(fT(Ts)−273.15)]})^0.963.

Since the outer wall temperature T_s and the frost layer thickness d_f on each section of the finned tube are different, the frost layer thermal resistance R_e of the finned tube in unit length is also different (as shown in FIG. 4), the relational expression is established between the frost layer thermal resistance R_f and T_s, and the frost layer thermal resistance R_f in each section of the finned tube in unit length may be represented by the outer wall temperature T_s of the finned tube; therefore, the increased frost layer thermal resistance R_e of the whole finned tube is not a constant value for each section any more, but is a linear value changing with the outer wall temperature T_s of each section of the finned tube; finally, the thermal conductivity coefficient λ_e of the equivalent thermal contact resistance R_e of the frosted finned tube is represented as the function of the outer wall temperature T_s of the finned tube; and the increased frost layer thermal resistance R_f in each section of each finned tube in unit length and the thermal conductivity coefficient λ_e are different, however, the thermal conductivity coefficients of all the sections are calculated through the function relational expression λ_e=F(T_s).

A deduction process of the function relational expression λ_e=F(T_s) of the thermal conductivity coefficient λ_e of the equivalent thermal contact resistance R_e of the frosted finned tube and the outer wall temperature T_s of the frosted finned tube is as follows:

Firstly, let a total heat transfer coefficient K_f (the total heat transfer coefficient is obtained according to a calculation formula for a heat transfer process through a ribbed wall in Heat Transfer 5th Edition (Higher Education Press, 2019) be equivalent to a total heat transfer coefficient K of the finned tube under a non-frosting condition, and then an expression of the equivalent thermal conductivity coefficient λ_e of the vaporizer material, namely aluminum alloy is obtained through transposition, as follows:

$\frac{1}{\frac{1}{h_{in}}+\frac{d_{in}}{\lambda }\ln \frac{d_{out}}{d_{in}}+\frac{1}{\eta ·\beta }\left(\frac{1}{h_{out}}+R_f\right)}=K_f=K=\frac{1}{\frac{1}{h_{in}}+\frac{d_{in}}{{\lambda }_e}\ln \frac{d_{out}}{d_{in}}+\frac{1}{\eta ·\beta ·h_{out}}}$ ${\lambda }_e=\frac{1}{\frac{1}{\lambda }+\frac{1}{d_{in}\ln \frac{d_{\mathrm{out}}}{d_{in}}}·\frac{h_f}{\eta ·\beta }}$

In the formula, K_f is the total heat transfer coefficient of the finned tube of the LNG ambient air vaporizer under the frosting condition, W/(m²·K); K is the total heat transfer coefficient of the finned tube of the LNG ambient air vaporizer under the non-frosting condition, W/(m²·K); h_in is a surface heat transfer coefficient in the finned tube, W/(m²·K); d_in is the internal diameter of the finned tube, m; λ is the thermal conductivity coefficient of the vaporizer material, namely aluminum alloy, W/(m·K); d_out is the external diameter of the finned tube, m; η is fin efficiency (the fin efficiency=actual heat dissipating capacity of the surface of the fin/heat dissipating capacity assuming that the outer wall temperature of the fin is equal to fin root temperature); β is a finning coefficient of the finned tube; h_out is a heat transfer coefficient on an air side outside the finned tube, W/(m²·K); R_f is the frost layer thermal resistance, (m²·K)/W; and λ_e is the equivalent thermal conductivity coefficient of the vaporizer material, namely the aluminum alloy, W/(m·K).

Parameter λ can be obtained in professional books (such as Practical Handbook of Nonferrous Metal Materials, Guangdong Science and Technology Press, 2006); d_in and d_out can be obtained from the design drawing; η, β, h_out, and R_f can be calculated through the following formula; and h_in is eliminated in calculation, which does not need to be solved.

The heat transfer coefficient h_out on the air side outside the finned tube is represented as an additive value of a convection heat exchange coefficient h_out,d on the air side outside the tube and a radiation heat exchange coefficient h_out,r on the air side outside the tube: h_out=h_out,d+h_out,r; and since the outer wall temperature T_s of the finned tube has little influence on the convection heat exchange coefficient h_out,d on the air side, the heat transfer coefficient can be calculated through the air flow velocity V_a on the outer side of the frost layer on the surface of the finned tube measured on site: h_out=18×V_a.

A calculation formula of the fin efficiency η is:

$\eta =\frac{A_2^{\prime }+A_2^{″}{\eta }_f}{A_2}$

In the formula, A₂ is the surface area outside the finned tube, m²; A₂ is a surface area of a non-fin part outside the finned tube in unit length, m²; A₂″ is a surface area of a fin part outside the finned tube in unit length, m²; A₂, A₂′, and A₂″ can be calculated from the design drawing; and η_f is fin surface efficiency, and a calculation formula is:

${\eta }_f=\frac{th\left(ml\right)}{ml}$ $m=\sqrt{\frac{2h_{out}}{\lambda \delta }}$

In the formula, m is the fin coefficient; l is the fin height, m; λ is the thermal conductivity coefficient of the vaporizer material, namely the aluminum alloy, W/(m·K); δ is the fin thickness, m; both l and δ can be obtained from the design drawing; h_out is the heat transfer coefficient on the air side outside the finned tube, W/(m²·K), when radiation heat exchange on the air side is neglected, h_out can be calculated according to the air flow velocity V_a outside the frost layer h_out=18×V_a, and h_out is expressed as the function h_out=18×V_a=18×f_V(T_s) with the outer wall temperature T_s of the finned tube according to the fitting relationship V_a=f_V(T_s) between the outer wall temperature T_s of the finned tube and the air flow velocity outside the frost layer, m=[2*18*f_V(T_s)/(λ*δ)]^1/2. Let A_m=l*[36/(λ*δ)]^1/2, then η_j=th(A_m*f_V(T_s)^1/2)/(A_m*f_V(T_s)^1/2).

A calculation formula of the finning coefficient β of the finned tube is:

$\beta =\frac{A_0}{A_2}$

In the formula, A₀ is the internal surface area of the finned tube, m², and A₂ is the external surface area of the finned tube, m², which can be calculated from the design drawing.

A calculation formula of the frost layer thermal resistance R_f is:

$R_f=\frac{d_f}{{\lambda }_f}$

In the formula, d_f is the frost layer thickness, mm; λ_f is the frost layer thermal conductivity coefficient, W/(m·K), which mainly depends on frost layer density ρ_f and can be calculated through Sanders relational expression (reference: Seker D, Karatas H, Egrican N. Frost formation on fin-and-tube heat exchangers. Part I—Modeling of frost formation on fin-and-tube heat exchangers [J].International Journal of Refrigeration, 2004, 27 (4): 367-374:

${\lambda }_f=0.001202{\rho }_f^{0.963}$

In the formula, ρ_f is the frost layer density, kg/m³; and a calculation formula is as follows:

${\rho }_f=650e^{0.277\left(T_f-273.15\right)}$

In the formula, T_f is the frost layer temperature, K.

The above two formulas are combined to obtain a formula of the frost layer thermal conductivity coefficient λ_f changing with the frost layer temperature T_f, as follows:

${\lambda }_f=0.001202\times {\left(650e^{0.277\left(T_f-273.15\right)}\right)}^{0.963}$

The function expression of the frost layer thermal conductivity coefficient λ_f and the outer wall temperature T_s of the finned tube is obtained according to the fitting relational expression T_f=f_T(T_s) between the outer wall temperature T_s of the finned tube and the frost layer temperature T_f, as follows:

${\lambda }_f=g\left(T_s\right)=0.001202\times {\left(650e^{0.277\left(f_T\left(T_s\right)-273.15\right)}\right)}^{0.963}$

The frost layer thermal resistance R_f is expressed as the function R_f=d_f/λ_f=f_d(T_s)/g(T_s)=Z(T_s) of the outer wall temperature T_s of the finned tube. FIG. 5 is a diagram of a fitting curve of frost layer thermal resistance R_f and the outer wall temperature T_s of the finned tube at different operating times.

R_f=Z(T_s), η=(A₂′+A₂″n_f)/A₂, and η_f=th(A_mf_V(T_s)^½)/(A_mf_V(T_s)^½) are substituted into the expression of the equivalent thermal conductivity coefficient λ_e of the vaporizer material, namely the aluminum alloy, to obtain:

${\lambda }_e=\frac{1}{\frac{1}{\lambda }+\frac{1}{d_{in}\ln \frac{d_{\mathrm{out}}}{d_{in}}}·\frac{h_f}{\eta ·\beta }}=\frac{1}{\frac{1}{\lambda }+\frac{1}{d_{in}\ln \frac{d_{\mathrm{out}}}{d_{in}}}·\frac{Z\left(T_s\right)}{\frac{A_2^{\prime }+A_2^{″}{n}_f}{A_2}·\beta }}=\frac{1}{\frac{1}{\lambda }+\frac{1}{d_{in}\ln \frac{d_{\mathrm{out}}}{d_{in}}}·\frac{Z\left(T_s\right)}{\frac{A_2^{\prime }\beta }{A_2}+\frac{A_2^{″}\beta }{A_2}·\frac{\mathrm{th}\left(A_m\sqrt{f_V\left(T_s\right)}\right)}{A_m\sqrt{f_V\left(T_s\right)}}}}$

In the formula, μ, d_in, d_out, A₂, A₂′, A₂″, A_m, and β are all constant values; and the equivalent thermal conductivity coefficient λ_e of the finned tube of the vaporizer is expressed as the function λ_e=F(T_s) of the outer wall temperature T_s of the finned tube.

Step 4, an overall geometric model of the LNG ambient air vaporizer is established in simulation software, meshing and dividing of computational domain are performed, physical models and equations are selected, material attributes and boundary conditions of the computational domain are set, the equivalent thermal conductivity coefficient λ_e of the finned tube of the vaporizer is adopted as the thermal conductivity coefficient during frosting of the vaporizer material, solving and initialization setting are performed, and then simulation and optimization calculation is performed. Details are as follows:

S1: as shown in FIG. 6 to FIG. 9, the overall geometric model of the LNG ambient air vaporizer in a ratio of 1:1 is established through three-dimensional geometric modeling software (such as SolidWorks or ANSYS DesignModeler), then the overall geometric model is subjected to meshing (as shown in FIG. 10, which is a schematic diagram of geometric model meshing through ANSYS Meshing) and dividing of the computational domain through finite element meshing software (such as ICEM CFD or ANSYS Meshing). The meshing needs to take the number, density and quality of meshes into account, so as to improve calculation efficiency and accuracy; the computational domain are divided into an LNG fluid domain, a vaporizer solid domain and an air fluid domain; the LNG fluid domain is a flow region of LNG in an internal channel of the vaporizer; the vaporizer solid domain is a vaporizer body; and the air fluid domain is an air flow region outside the vaporizer body.

S2: the meshed overall geometric model of the LNG ambient air vaporizer is imported into fluid analysis software (such as ANSYS Fluent), and the LNG fluid domain, the vaporizer solid domain and the air fluid domain are adopted as the computational domain; and a contact surface between the LNG fluid domain and the vaporizer solid domain and a contact surface between the vaporizer solid domain and the air fluid domain are set as Interface surfaces, and a Couple option checked set in Interface setting, so that the corresponding contact surface can complete heat transfer.

S3: a gravity model, a multi-phase model, a turbulence model, a boiling phase change model, a continuity equation, a momentum equation, an energy equation and a component transport equation are enabled in the fluid analysis software, and a standard wall function method is adopted for near-wall processing; a Mixture model is adopted as the multi-phase model; and a Realizable k-epsilon turbulence model is adopted as the turbulence model, and an evaporation-condensation Lee model is adopted as the boiling phase change model.

S4: the material attributes of the computational domain are set:

LNG and NG fluid materials are respectively introduced in the fluid analysis software, physical parameter data about LNG and NG in relevant books (Technical Handbook of Liquefied Natural Gas, China Machine Press, 2010) are adopted as material parameters, then the LNG fluid material is set as a first term in the multi-phase model, phase change from LNG to NG is set, and the evaporation-condensation Lee model is selected for a reaction mechanism; an aluminum alloy solid material is introduced in the fluid analysis software, physical data in a software material library is adopted as parameters of the aluminum alloy solid material, then the thermal conductivity coefficient of the aluminum alloy solid material is modified from the constant value λ to a value represented by a piecewise polynomial temperature function method, the thermal conductivity coefficient of the vaporizer material within a frosting temperature range is set as the equivalent thermal conductivity coefficient λ_e=F(T_s), and the thermal conductivity coefficient of the vaporizer material with a non-frosting temperature range is set as the constant value λ; and a wet air mixed material is introduced in the fluid analysis software, the wet air mixed material including air and water vapor, and physical property data in the software material library are adopted as material attributes of the wet air mixed material.

S5: the boundary conditions of the computational domain are set:

An outlet of the LNG fluid domain is set as a pressure outlet boundary, and the pressure P_out at the vaporizer outlet tested on site is adopted as pressure; an inlet of the LNG fluid domain is set as a velocity inlet boundary, and the flow velocity V_in and temperature T_in at the vaporizer inlet tested on site are adopted as velocity and temperature; a top surface of the air fluid domain above the vaporizer and a side surface of the air fluid domain around the vaporizer are set as pressure inlet boundaries, the atmospheric pressure P₀ and ambient temperature T₀ tested on site are set as pressure and temperature, and air humidity of the air fluid domain for simulating the air humidity is set according to the ambient humidity tested on site; and a bottom surface of the air fluid domain at a bottom of the vaporizer is set as a pressure outlet boundary.

In this step, the air fluid domain is defined as a hexahedron capable of surrounding the vaporizer, space, except for the vaporizer, inside the hexahedron represents air outside the vaporizer, and the top surface, the side surface and the bottom surface of the air fluid domain refer to a top surface, a side surface and a bottom surface of the hexahedron outside the whole vaporizer.

S6: initialization setting is performed by adopting a SIMPLE algorithm in the fluid analysis software as a solving method for the geometric model meshes divided in step S1, and then calculation simulation is performed on the geometric model established in step S1; calculation stops and result data from numerical simulation in heat transfer of the vaporizer is output if a residual variance curve converges and the monitored NG outlet temperature and flow velocity do not change any more; otherwise, it proceeds to operate.

Furthermore, in order to make the initialization result as close as possible to the actual physical result to ensure the stability of the calculation process and increase the convergence velocity, an LNG volume fraction of an inlet of the LNG fluid domain is set as 1 and temperature is set as 123 K through a Patch function in the fluid analysis software; temperature of the vaporizer solid domain is set as 260 K; a volume fraction of air in the air fluid domain outside the finned tube is set as 1, and temperature is determined through a self-defined formula: T=279+7 z, where, T is air temperature, K; and z is a height in an z-axis direction (z-axis is set to be perpendicular to the bottom of the vaporizer, which is a vertically upward direction), m.

Step 5, the result data from numerical simulation is imported into post-processing software (such as Tecplot or Ensight) for analysis, a temperature cloud diagram on a surface of the LNG ambient air vaporizer, and a temperature cloud diagram, a velocity cloud diagram and a component cloud diagram of the LNG fluid domain in the tube of the vaporizer are displayed in the post-processing software, an outlet section of the LNG fluid domain is selected, to obtain LNG outlet temperature and outlet flow velocity, and the component cloud diagram of the LNG fluid domain is selected, to obtain proportions of a liquid phase section, a two-phase section and a gas-phase section in the fluid domain; and temperature, heat flux and other thermal parameters at different position points or sections of the surface of the vaporizer are checked, so as to visually acquire fluid-solid conjugate heat transfer characteristics and vaporization performance of the LNG ambient air vaporizer under the frosting condition.

Analysis results are shown in FIG. 11 and FIG. 12, and FIG. 11 shows proportions of liquid phase sections, two-phase sections and gas liquid sections in branch tubes of different vaporizer finned tubes under a frosting condition; and FIG. 12 shows comparison of simulation results versus actual measurement results of vaporizer outlet temperature at different operating times, and it can be shown that errors between the simulation results of the method and the actual measurement results are small, thereby realizing good accuracy.

The above description is only the preferred embodiments of the present application, and is not intended to limit the present application, and for those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall fall within the protection scope of the present application.

The present disclosure can be applied to optimize design of the vaporizer: according to the outputted result of Step 5, LNG outlet temperature T_out, outlet flow velocity v_out and proportions of a liquid phase section, a two-phase section and a gas-phase section in the fluid domain, modify design parameters of the LNG ambient air vaporizer to satisfy LNG ambient air vaporizer under a frosting condition; wherein the design parameters including an internal diameter d_in of the finned tube, number n of the fin, height l and a thickness δ of a fin.

The details are as follows: according to the output LNG outlet temperature T_out, outlet flow velocity v_out at the certain operating moment t_n (such as 1 h, 2 h, 4 h, 8 h), obtaining vaporization amount Q_out, and the calculation formula is Q_out=3600*v_out*A; where A in the formula is the cross-sectional area of the outlet of the vaporizer; determining whether the vaporization amount Q_out at the certain operating moment t_n satisfy the designed requirement or not; if the requirement is not met, increase the design parameter of the inner diameter d_in of the finned tube until the parameter satisfy the designed requirement; and determining whether the LNG outlet temperature T_out meets the designed requirement or not according to the LNG outlet temperature T_out, outlet flow velocity v_out and proportions of a liquid phase section, a two-phase section and a gas-phase section in the fluid domain obtained in step 5; if the requirement is not met and the liquid phase segment in the fluid domain is relatively high, increase design parameters of number n of the fin and fin height l, and decrease the thickness δ of a fin; and if the gas phase segment in the fluid domain is relatively high, decrease design parameters of number n of the fin and fin height l.

The present disclosure can be further applied in the following scenario: according to the outputted result of Step 5, evaluation of the vaporization performance of an LNG ambient air vaporizer under a frosting condition is performed as follows. According to the output LNG outlet temperature Tout, outlet flow velocity v_out at the certain operating moment t_n (such as 1 h, 2 h, 4 h, 8 h), obtaining vaporization amount Q_out, and the calculation formula is Q_out=3600*v_out*A; where A in the formula is the cross-sectional area of the outlet of the vaporizer; determining whether the vaporization amount Q_out at the certain operating moment t_n satisfy the designed requirement or not; if the requirement is not met, increase the design parameter of the outlet flow velocity v_out and pressure P_in of liquefied natural gas at a vaporizer inlet; And determining the difference between the LNG outlet temperature T_out and site ambient temperature T₀, if T_out≥T₀−10, no additional operation is required; otherwise, performing fin defrosting operation, defrosting frequency f=24/t_n is calculated; meanwhile, a reheater connected to the vaporizer outlet is started, so as to increase the temperature to +5° C. or above, and it is ensured that the outlet temperature of the reheater is not less than 5° C.

The present disclosure can be implemented by a processing system. The processing system includes a processor coupled to a computer-readable medium/memory. The processor is responsible for general processing, including the execution of software stored on the computer-readable medium/memory. The software, when executed by the processor, causes the processing system to perform the various functions described supra for any particular apparatus. The computer-readable medium/memory may also be used for storing data that is manipulated by the processor when executing software. The processing system 1114 further includes at least one or more components that may be software components running in the processor, resident/stored in the computer readable medium/memory, one or more hardware components coupled to the processor, or some combination thereof.

On a basis of the same concept, an electronic device is also provided in one or more embodiments of this specification, which includes a memory, a processor and a computer program stored on the memory and operable on the processor, and the processor implements the method described in any one of the above when executing the program.

The electronic device provided in this embodiment, which may include a processor, a memory, an input/output interface, a communication interface and a bus. A communication connection with each other is realized among the processor, the memory, the input/output interface and the communication interface in the device through the bus.

The processor can be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application specific integrated circuit (ASIC), or one or more integrated circuits, etc., and is used for executing relevant programs to realize technical schemes provided in the embodiments of this specification.

The memory can be implemented in a form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, or the like. The memory can store the operating system and other application programs. When the technical schemes provided in the embodiments of this specification is implemented by software or firmware, relevant program codes are stored in the memory and called and executed by the processor.

The input/output interface is used to connect the input/output module to realize information input and output. The input/output module can be configured as a component in the device (not shown in the figure), or it can be externally connected to the device to provide corresponding functions. An input device can include a keyboard, a mouse, a touch screen, a microphone and various sensors, and an output device can include a display, a speaker, a vibrator and an indicator.

The communication interface is used to connect with a communication module (not shown in the figure) to realize communication interaction between the device and other devices. The communication module can communicate by wired means (such as USB, network cable, etc.) or by wireless means (such as mobile network, WIFI, Bluetooth, etc.).

The bus includes a path for transmitting information among various components of the device, such as the processor, the memory, the input/output interface and the communication interface.

The electronic devices in the above embodiments are used to implement the corresponding method in any of the above embodiments, and present beneficial effects of the corresponding method embodiments, which will not be described in detail here.

It should be understood by those of ordinary skill in the art that discussion of any of the above embodiments is only exemplary, and is not intended to imply that the scope of the disclosure (including the claims) is limited to these examples; under the idea of this disclosure, the technical features in the above embodiments or different embodiments can also be combined, and the steps can be realized in any order; and there are many other changes in different aspects of the embodiments of this disclosure as described above, which are not provided in details for brevity.

In addition, in order to simplify the description and discussion, and in order not to make the embodiments of the present disclosure difficult to understand, well-known power/ground connections with integrated circuit (IC) chips and other components may or may not be shown in the provided drawings. In addition, the devices can be shown in the form of block diagrams in order to avoid making the embodiments of the present disclosure difficult to understand, and this also takes into account the fact that the details about the implementation of these devices in block diagram are highly dependent on a platform on which the embodiments of the present disclosure will be implemented (i.e., these details should be completely within an understanding range of those skilled in the art). While specific details (e.g., circuits) have been set forth to describe exemplary embodiments of the present disclosure, it will be apparent to those skilled in the art that the embodiments of the present disclosure may be practiced without these specific details or with variations in these specific details. Therefore, these descriptions should be regarded as illustrative rather than restrictive.

Although the present disclosure has been described in conjunction with specific embodiments thereof, many alternatives, modifications and variations of these embodiments will be apparent to those of ordinary skill in the art based on the foregoing description.

The embodiments of the present disclosure are intended to cover all such alternatives, modifications and variations that fall within the broad scope of the append claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of this disclosure shall be encompassed within the protection scope of this disclosure.

Claims

1. A method for simulating performance and optimizing design of an LNG ambient air vaporizer under a frosting condition, comprising the following steps, executed by a processor: step 1, performing site operation test on the LNG ambient air vaporizer, measuring site ambient temperature T0, humidity H0 and atmospheric pressure P0, and measuring, at a certain operating moment tn, pressure Pin of liquefied natural gas at a vaporizer inlet and pressure Pout of natural gas at a vaporizer outlet, flow velocity Vin of the liquefied natural gas at the vaporizer inlet and temperature Tin of the liquefied natural gas at the vaporizer inlet, outer wall temperature Ts at different positions of each frosted finned tube of the vaporizer, frost layer temperature Tf, frost layer thickness df, and air flow velocity Va outside a frost layer on a surface of the finned tube, the different positions referring to at least three equidistant point positions of an outer wall of a fin of each frosted finned tube from top to bottom, wherein, the site ambient temperature T0, humidity H0 and atmospheric pressure P0 are measured by a temperature and humidity and atmospheric pressure instrument; pressure Pin and pressure Pout are measured by a pressure gauge disposed at the inlet and outlet of the air temperature gasifier, respectively; the flow velocity Vin is measured by a flow meter provided at an inlet of the air temperature gasifier; the temperature Tin and outer wall temperature Ts are measured by temperature sensors, which are arranged at the inlet of the gasifier and a different location of each finned tube of the gasifier; the frost layer temperature Tf is measured by temperature gauges; the frost layer thickness df is measured by vernier calipers; and the air flow velocity Va is measured by anemometers;step 2, performing by data processing software, fitting analysis on the collected frost layer temperature Tf and frost layer thickness df at the different positions of all the frosted finned tubes at the certain operating moment tn and the air flow velocity Va outside the frost layer on the surface of the finned tube, and data of the outer wall temperature Ts of the finned tube respectively, and obtaining, by a least square method, fitting relational expressions df=fd(Ts)=A1Ts2+B1Ts+C1, TffT(Ts)=A2Ts+B2, and Va=fV(Ts)=A3Ts2+B3Ts+C3 between the outer wall temperature Ts of all the frosted finned tubes of the whole vaporizer and the frost layer thickness df, the frost layer temperature Tf and the air flow velocity Va outside the frost layer on the surface of the finned tube at the certain operating moment tn, A1, B1, C1, A2, B2, A3, B3, and C3 in the formula being respectively fitted constants;step 3, establishing a calculation model for an equivalent thermal conductivity coefficient of the LNG ambient air vaporizer during frosting at the certain operating moment tn:transforming a sum of increased frost layer thermal resistance Rf of each frosted finned tube in unit length at the certain operating moment tn and thermal resistance Ro of the finned tube body in unit length into equivalent thermal contact resistance Re of the finned tube in unit length under a non-frosting condition, and representing an equivalent thermal conductivity coefficient λe of the equivalent thermal contact resistance Re as a function of the outer wall temperature Ts of the frosted finned tube in unit length, the unit length being a length of a minimum mesh during geometric meshing of the vaporizer; wherein

2. The method according to claim 1, wherein, in the step S5 of the step 4, an LNG volume fraction of an inlet of the LNG fluid domain is set as 1 and temperature is set as 123 K through a Patch function in the fluid analysis software; temperature of the vaporizer solid domain is set as 260 K; a volume fraction of air in the air fluid domain outside the finned tube is set as 1, and temperature is determined through a self-defined formula: T=279+7 z, where, T is air temperature, K; and z is a height in an z-axis direction, and z-axis is set to be perpendicular to the bottom of the vaporizer, which is a vertically upward direction.