Elliptical tanks were used as an alternative to circular tanks in order to improve mixing efficiency and reduce mixing time in unbaffled stirred tanks(USTs). Five different aspect ratios of elliptical vessels were des...Elliptical tanks were used as an alternative to circular tanks in order to improve mixing efficiency and reduce mixing time in unbaffled stirred tanks(USTs). Five different aspect ratios of elliptical vessels were designed to compare their mixing time and flow field. Computational fluid dynamics(CFD) simulations were performed using the k–ε model to calculate the mixing time and simulate turbulent flow field features, such as streamline shape, velocity distribution, vortex core region distribution, and turbulent kinetic energy(TKE) transfer. Visualization was also carried out to track the tinctorial evolution of the liquid phase. Results reveal that elliptical stirred tanks can significantly improve mixing performance in USTs. Specifically, the mixing time at an aspect ratio of 2.00 is only 45.3% of the one of a circular stirred tank. Furthermore, the secondary flow is strengthened and the vortex core region increases with the increase of aspect ratio. The axial velocity is more sensitive to the aspect ratio than the circumferential and radial velocity. Additionally, the TKE transfer in elliptical vessels is altered. These findings suggest that elliptical vessels offer a promising alternative to circular vessels for enhancing mixing performance in USTs.展开更多
The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computatio...The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computational fluid dynamics(CFD) software package Fluent 6.2, the mixing characteristics in a tank agitated by dual six-blade-Rushton-turbines(6-DT) are predicted using the detached eddy simulation(DES) method. A sliding mesh(SM) approach is adopted to solve the rotation of the impeller. The simulated flow patterns and liquid velocities in the agitated tank are verified by experimental data in the literature. The simulation results indicate that the DES method can obtain more flow details than Reynolds-averaged Navier-Stokes(RANS) model. Local and global mixing time in the agitated tank is predicted by solving a tracer concentration scalar transport equation. The simulated results show that feeding points have great influence on mixing process and mixing time. Mixing efficiency is the highest for the feeding point at location of midway of the two impellers. Two methods are used to determine global mixing time and get close result. Dimensionless global mixing time remains unchanged with increasing of impeller speed. Parallel, merging and diverging flow pattern form in the agitated tank, respectively, by changing the impeller spacing and clearance of lower impeller from the bottom of the tank. The global mixing time is the shortest for the merging flow, followed by diverging flow, and the longest for parallel flow. The research presents helpful references for design, optimization and scale-up of agitated tanks with multi-impeller.展开更多
The macroscopic mixing in a stirred tank with different tracer injection locations, impeller speeds and impeller positions is simulated numerically by solving the transport equation of the tracer based on the whole fl...The macroscopic mixing in a stirred tank with different tracer injection locations, impeller speeds and impeller positions is simulated numerically by solving the transport equation of the tracer based on the whole flow field in the baffled tank with a Rushton disk turbine numerically resolved using the improved inner-outer iterative procedure. Predicted mixing time is compared well with the literature correlations. The predicted residence time distribution of the stirred tank is very close to the present experimental results. The effect of the installation of a draft tube on the mixing time and residence time distributions is addressed.展开更多
The effect of mixing on the precipitation of barium sulfate in a continuous stirred tank is simulated nu- merically with different feeding location, feed concentration, impeller speed and residence time through solvin...The effect of mixing on the precipitation of barium sulfate in a continuous stirred tank is simulated nu- merically with different feeding location, feed concentration, impeller speed and residence time through solving the standard momentum and mass transport equations in combination with the moment equations for crystal population balance. The numerical method was validated with the literature data. The simulation results including the distribu- tion of the local supersaturation ratio distribution in the precipitator, mean crystal size and coefficient of variation under different operating conditions compared well with experimental data in the literature. The effect of the pres- ence of a draft tube on precipitation were also investigated, and it is suggested that the installation of a draft tube increased the mean crystal size, in general agreement with experimental work in the literature.展开更多
Laminar mixing in the stirred tank is widely encountered in chemical and biological industries.Isolated mixing regions(IMRs)usually exist when the fluid medium has high viscosity,which are not conducive to mixing.In t...Laminar mixing in the stirred tank is widely encountered in chemical and biological industries.Isolated mixing regions(IMRs)usually exist when the fluid medium has high viscosity,which are not conducive to mixing.In this work,the researches on IMRs,enhancement of laminar mixing and the phenomenon of particle clustering within IMRs are reviewed.For most studies,the aim is to destroy IMRs and improve the chaotic mixing.To this end,the mechanism of chaotic mixing and the structure of IMRs were well investigated.The methods developed to destroy IMRs include off-centered agitation,dynamic mixing protocol,special designs of impellers,baffles,etc.In addition,the methods to characterize the shape and size of IMRs as well as mixing effect by experiments and simulations are summarized.However,IMRs are not always nuisance,and it may be necessary in some situations.Finally,the present engineering applications are summarized,and the prospect of the future application is predicted.For example,particle clustering will form in the co-existing system of chaotic mixing and IMRs,which can be used for solid–liquid separation and recovery of particles from high viscosity fluid.展开更多
A novel control strategy for a continuous stirred tank reactor(CSTR)system,which has the typical characteristic of strongly pronounced nonlinearity,multiple operating points,and a wide operating range,is initiated fro...A novel control strategy for a continuous stirred tank reactor(CSTR)system,which has the typical characteristic of strongly pronounced nonlinearity,multiple operating points,and a wide operating range,is initiated from the point of hybrid systems.The proposed scheme makes full use of the modeling power of mixed logical dy- namical(MLD)systems to describe the highly nonlinear dynamics and multiple operating points in a unified framework as a hybrid system,and takes advantage of the good control quality of model predictive control(MPC) to design a controller.Thus,this approach avoids oscillation during switching between sub-systems,helps to relieve shaking in transition,and augments the stability robustness of the whole system,and finally achieves optimal(i.e. fast and smooth)transition between operating points.The simulation results demonstrate that the presented ap- proach has a satisfactory performance.展开更多
Mechanical agitation in baffled vessels with turbines plays a vital role in achieving homogeneous fluid mixing and promoting various transfer operations.Therefore,designing vessels with optimal energy efficiency and f...Mechanical agitation in baffled vessels with turbines plays a vital role in achieving homogeneous fluid mixing and promoting various transfer operations.Therefore,designing vessels with optimal energy efficiency and flow dynamics is essential to enhance operational performance and eliminate flow perturbations.Hence,the present research focuses on a numerical investigation of the impact of inclined slots with different angles installed at the side-wall of a cylindrical vessel equipped with a Rushton turbine.This study explores power consumption and vortex size while considering various rotation directions of the impeller with different rotation speeds.The numerical simulations are conducted for Reynolds numbers ranging from 104 to 105,using the RANS k-εturbulence model to govern the flow inside the stirred vessel,accounting for mass and momentum balances.The results have shown that the installation of slots reduces power consumption and vortex size compared to conventional vessel configu-rations.Moreover,increasing the slot angle from 0 to 32.5°further reduces energy consumption and vortex size,especially with negative rotation speeds.On the other hand,increasing the Reynolds numbers leads to a decrease in power consumption and an increase in vortex size.The present research therefore proposes a design for con-structing Rushton-turbine stirred vessels offering optimal operation,characterized by reduced energy consumption and minimized vortex size.展开更多
The discrete element method(DEM)was used in this study to numerically simulate the mixing process and motion law of particles in brown rice germination device.And the reliability of simulation experiments was verified...The discrete element method(DEM)was used in this study to numerically simulate the mixing process and motion law of particles in brown rice germination device.And the reliability of simulation experiments was verified through physical experiments.In the discrete element simulation experiment,there were three mixing stages in the mixing process of the particles.The particle motion conditions at different rotational speeds were rolling,cascading,cataracting and centrifuging.The lower the filling degree,the higher the particle mixing efficiency.The radial trajectory of the particles was approximated as an elliptical helix that continuously shrank towards the axis.The research results indicated that under the same speed and filling conditions,the motion of brown rice particles in both the simulated and physical test environments is rolling and the drop height is the same.展开更多
Macro-instability(MI)is an important natural phenomenon affecting the mixing performance in stirred tank significantly,which results from movement and evolution of large eddies.Factors associated with MI frequency in ...Macro-instability(MI)is an important natural phenomenon affecting the mixing performance in stirred tank significantly,which results from movement and evolution of large eddies.Factors associated with MI frequency in mixing operation and energy dissipation related to MI were reviewed.Flow pattern in turbulent regime contains many coherent structures,which have much energy vortices and may result in flow field MI.Adjustment of coherent structure and controlling MI may contribute to energy saving and good mixing performance.Controlling methods for MI was prospected.展开更多
The multiphase flow in the solid-liquid tank stirred with a new structure of Intermig impeller was analyzed by computational fluid dynamics(CFD).The Eulerian multiphase model and standard k-ε turbulence model were ad...The multiphase flow in the solid-liquid tank stirred with a new structure of Intermig impeller was analyzed by computational fluid dynamics(CFD).The Eulerian multiphase model and standard k-ε turbulence model were adopted to simulate the fluid flow,turbulent kinetic energy distribution,mixing performance and power consumption in a stirred tank.The simulation results were also verified by the water model experiments,and good agreement was achieved.The solid-liquid mixing performances of Intermig impeller with different blade structures were compared in detail.The results show that the improved Intermig impeller not only enhances the solid mixing and suspension,but also saves more than 20% power compared with the standard one.The inner blades have relatively little influence on power and the best angle of inner blades is 45°,while the outer blades affect greatly the power consumption and the optimized value is 45°.展开更多
The aim of this paper is to study the effect of agitator’s types on the turbulent flows in stirred tanks without and with baffles.The hydrodynamics behavior induced by four different agitator’s types:a Rushton turbi...The aim of this paper is to study the effect of agitator’s types on the turbulent flows in stirred tanks without and with baffles.The hydrodynamics behavior induced by four different agitator’s types:a Rushton turbine(RT),a circular blade turbine(CBT),a diverging triangular blade turbine(DTBT)and converging triangular blade turbine(CTBT)are numerically predicted by solving the Navier-Stokes equations and RNGκ–εturbulent model.The simulations are carried out using the Multi Reference Frame(MRF)approach.The numerical results showed good agreement with experiment.We find that the agitator CTBT gives an important profit on the power consumption per report/ratio the others and DTBT give a good reduction of the vortex size of the impeller angles.展开更多
For the mixing-sensitive reactions, both chemical kinetics and mixing conditions of the reactants determine the distributions of products. The direct quadrature method of moments combining with the interaction by exch...For the mixing-sensitive reactions, both chemical kinetics and mixing conditions of the reactants determine the distributions of products. The direct quadrature method of moments combining with the interaction by exchange with the mean micro-mixing model(DQMOM-IEM) has been validated for the chemical reacting flows in microreactors. Quite encouraging simulation results offer great promise, but the applicability of this method is needed to be explored furthermore, such as in stirred reactors. In this work, the two-environment DQMOM-IEM model was created with C language and used to customize Fluent through the user-defined functions. The mixing effects on the course of parallel competing chemical reactions carried out in a semi-batch single-phase stirred reactor were predicted. The simulation results show that the rising feed velocity enlarges the volume of reaction zone and maximize the yield of the by-product, which also indicates that the feed stream is more difficultly dispersed into the main stream and the zone surrounding feedpipe exit with high turbulent kinetic dissipation rate cannot be efficiently used.展开更多
Liquefied petroleum gas (LPG) is an important urban gas source in China. Before supplied to customers by pipeline supply systems, LPG is stored in tanks in LPG vaporization stations. Designers usually decide the numbe...Liquefied petroleum gas (LPG) is an important urban gas source in China. Before supplied to customers by pipeline supply systems, LPG is stored in tanks in LPG vaporization stations. Designers usually decide the number and the size of storage tanks by their experience during constructions of vaporization stations. These decisions are usually not best and most economical. To solve the problem, a compact mixed integer nonlinear programming model has been developed in this paper. The objective is to minimize annual storage cost of the vaporization station. The model has been transformed into a general nonlinear programming model by transforming integer variables and 0-1 variables into continuous variables. One LPG vaporization station was taken as an example to illustrate the usage of the model. The results show that the optimal storage scheme can be determined accurately and quickly by the model and about 15% of storage cost can be saved every year after optimization.展开更多
基金supported by the National Key Research and Development Project(2022YFB3504305,2019YFC1905802)National Natural Science Foundation of China(22078030)+2 种基金Joint Funds of the National Natural Science Foundation of China(U1802255)Key Project of Independent Research Project of State Key Laboratory of Coal Mine Disaster Dynamics and Control(2011DA105287-zd201902)Three Gorges Laboratory Open Fund of Hubei Province(SK211009,SK215001).
文摘Elliptical tanks were used as an alternative to circular tanks in order to improve mixing efficiency and reduce mixing time in unbaffled stirred tanks(USTs). Five different aspect ratios of elliptical vessels were designed to compare their mixing time and flow field. Computational fluid dynamics(CFD) simulations were performed using the k–ε model to calculate the mixing time and simulate turbulent flow field features, such as streamline shape, velocity distribution, vortex core region distribution, and turbulent kinetic energy(TKE) transfer. Visualization was also carried out to track the tinctorial evolution of the liquid phase. Results reveal that elliptical stirred tanks can significantly improve mixing performance in USTs. Specifically, the mixing time at an aspect ratio of 2.00 is only 45.3% of the one of a circular stirred tank. Furthermore, the secondary flow is strengthened and the vortex core region increases with the increase of aspect ratio. The axial velocity is more sensitive to the aspect ratio than the circumferential and radial velocity. Additionally, the TKE transfer in elliptical vessels is altered. These findings suggest that elliptical vessels offer a promising alternative to circular vessels for enhancing mixing performance in USTs.
基金Supported by Key Scientific Research Project of Sichuan Provincial Education Department(Grant No.15ZA0107)Doctor Foundation of Southwest University of Science and Technology(Grant No.11zx7162)
文摘The Issue of mixing efficiency in agitated tanks has drawn serious concern in many industrial processes. The turbulence model is very critical to predicting mixing process in agitated tanks. On the basis of computational fluid dynamics(CFD) software package Fluent 6.2, the mixing characteristics in a tank agitated by dual six-blade-Rushton-turbines(6-DT) are predicted using the detached eddy simulation(DES) method. A sliding mesh(SM) approach is adopted to solve the rotation of the impeller. The simulated flow patterns and liquid velocities in the agitated tank are verified by experimental data in the literature. The simulation results indicate that the DES method can obtain more flow details than Reynolds-averaged Navier-Stokes(RANS) model. Local and global mixing time in the agitated tank is predicted by solving a tracer concentration scalar transport equation. The simulated results show that feeding points have great influence on mixing process and mixing time. Mixing efficiency is the highest for the feeding point at location of midway of the two impellers. Two methods are used to determine global mixing time and get close result. Dimensionless global mixing time remains unchanged with increasing of impeller speed. Parallel, merging and diverging flow pattern form in the agitated tank, respectively, by changing the impeller spacing and clearance of lower impeller from the bottom of the tank. The global mixing time is the shortest for the merging flow, followed by diverging flow, and the longest for parallel flow. The research presents helpful references for design, optimization and scale-up of agitated tanks with multi-impeller.
基金Supported by the National Natural Science Foundation of China (No.20236050 50134020) The National Basic Research Priorities Program (973) (No.2004CB217604)
文摘The macroscopic mixing in a stirred tank with different tracer injection locations, impeller speeds and impeller positions is simulated numerically by solving the transport equation of the tracer based on the whole flow field in the baffled tank with a Rushton disk turbine numerically resolved using the improved inner-outer iterative procedure. Predicted mixing time is compared well with the literature correlations. The predicted residence time distribution of the stirred tank is very close to the present experimental results. The effect of the installation of a draft tube on the mixing time and residence time distributions is addressed.
基金National Natural Science Foundation of China (Nos.20236050, 50134020) the Special Funds for MajorState Basic Research Program of China (973 Program, 2004CB217604)
文摘The effect of mixing on the precipitation of barium sulfate in a continuous stirred tank is simulated nu- merically with different feeding location, feed concentration, impeller speed and residence time through solving the standard momentum and mass transport equations in combination with the moment equations for crystal population balance. The numerical method was validated with the literature data. The simulation results including the distribu- tion of the local supersaturation ratio distribution in the precipitator, mean crystal size and coefficient of variation under different operating conditions compared well with experimental data in the literature. The effect of the pres- ence of a draft tube on precipitation were also investigated, and it is suggested that the installation of a draft tube increased the mean crystal size, in general agreement with experimental work in the literature.
基金supports from National Key Research and Develop-ment Program(2020YFA0906804)the National Natural Science Foundation of China(21776282,21978296 and 22078229)+4 种基金the NSFC Key Program(21938009)major project(91934301)the National Key R&D Program of China(2019YFC1905805)Chemistry and Chemical Engineering Guangdong Laboratory Shantou(1922006)Innovation Academy for Green Manufacture,Chinese Academy of Sciences(IAGM2020C06)。
文摘Laminar mixing in the stirred tank is widely encountered in chemical and biological industries.Isolated mixing regions(IMRs)usually exist when the fluid medium has high viscosity,which are not conducive to mixing.In this work,the researches on IMRs,enhancement of laminar mixing and the phenomenon of particle clustering within IMRs are reviewed.For most studies,the aim is to destroy IMRs and improve the chaotic mixing.To this end,the mechanism of chaotic mixing and the structure of IMRs were well investigated.The methods developed to destroy IMRs include off-centered agitation,dynamic mixing protocol,special designs of impellers,baffles,etc.In addition,the methods to characterize the shape and size of IMRs as well as mixing effect by experiments and simulations are summarized.However,IMRs are not always nuisance,and it may be necessary in some situations.Finally,the present engineering applications are summarized,and the prospect of the future application is predicted.For example,particle clustering will form in the co-existing system of chaotic mixing and IMRs,which can be used for solid–liquid separation and recovery of particles from high viscosity fluid.
基金Supported by the National Natural Science Foundation of China (No.60404018) and the State Key Development Program for Basic Research of China (No.2002CB312200).
文摘A novel control strategy for a continuous stirred tank reactor(CSTR)system,which has the typical characteristic of strongly pronounced nonlinearity,multiple operating points,and a wide operating range,is initiated from the point of hybrid systems.The proposed scheme makes full use of the modeling power of mixed logical dy- namical(MLD)systems to describe the highly nonlinear dynamics and multiple operating points in a unified framework as a hybrid system,and takes advantage of the good control quality of model predictive control(MPC) to design a controller.Thus,this approach avoids oscillation during switching between sub-systems,helps to relieve shaking in transition,and augments the stability robustness of the whole system,and finally achieves optimal(i.e. fast and smooth)transition between operating points.The simulation results demonstrate that the presented ap- proach has a satisfactory performance.
文摘Mechanical agitation in baffled vessels with turbines plays a vital role in achieving homogeneous fluid mixing and promoting various transfer operations.Therefore,designing vessels with optimal energy efficiency and flow dynamics is essential to enhance operational performance and eliminate flow perturbations.Hence,the present research focuses on a numerical investigation of the impact of inclined slots with different angles installed at the side-wall of a cylindrical vessel equipped with a Rushton turbine.This study explores power consumption and vortex size while considering various rotation directions of the impeller with different rotation speeds.The numerical simulations are conducted for Reynolds numbers ranging from 104 to 105,using the RANS k-εturbulence model to govern the flow inside the stirred vessel,accounting for mass and momentum balances.The results have shown that the installation of slots reduces power consumption and vortex size compared to conventional vessel configu-rations.Moreover,increasing the slot angle from 0 to 32.5°further reduces energy consumption and vortex size,especially with negative rotation speeds.On the other hand,increasing the Reynolds numbers leads to a decrease in power consumption and an increase in vortex size.The present research therefore proposes a design for con-structing Rushton-turbine stirred vessels offering optimal operation,characterized by reduced energy consumption and minimized vortex size.
基金the National Natural Science Foundation of China(Grant No.32001423)Natural Science Foundation of Hubei Province(Grant No.2020CFB471)+2 种基金Huazhong Agricultural University College Students Science and Technology Innovation Fund Project(Grant No.2022255)Fundamental Research Funds for the Central Universities(Grant No.2662020GXPY017)First Division Alar City Science and Technology Plan Project(Grant No.2023ZB01)for financial support and all of the persons who assisted in this writing.
文摘The discrete element method(DEM)was used in this study to numerically simulate the mixing process and motion law of particles in brown rice germination device.And the reliability of simulation experiments was verified through physical experiments.In the discrete element simulation experiment,there were three mixing stages in the mixing process of the particles.The particle motion conditions at different rotational speeds were rolling,cascading,cataracting and centrifuging.The lower the filling degree,the higher the particle mixing efficiency.The radial trajectory of the particles was approximated as an elliptical helix that continuously shrank towards the axis.The research results indicated that under the same speed and filling conditions,the motion of brown rice particles in both the simulated and physical test environments is rolling and the drop height is the same.
基金financially supported by Natural Science Foundation of China(20806095)ScientificResearch Foundation for the Returned Overseas Chinese Scholars,State Education Ministry(教外司留20091341-2)+1 种基金China Postdoctoral Science Foundation(20080430747)National High Technology Researchand Development Program of China(2008AA031202)
文摘Macro-instability(MI)is an important natural phenomenon affecting the mixing performance in stirred tank significantly,which results from movement and evolution of large eddies.Factors associated with MI frequency in mixing operation and energy dissipation related to MI were reviewed.Flow pattern in turbulent regime contains many coherent structures,which have much energy vortices and may result in flow field MI.Adjustment of coherent structure and controlling MI may contribute to energy saving and good mixing performance.Controlling methods for MI was prospected.
基金supported by the National Natural Science Foundation of China (U1760120, U1508217)National Key R&D Program of China (2017YFC0210403, 2017YFC0210404)Shenyang Science & Technology Project (17-500-8-01)。
文摘The multiphase flow in the solid-liquid tank stirred with a new structure of Intermig impeller was analyzed by computational fluid dynamics(CFD).The Eulerian multiphase model and standard k-ε turbulence model were adopted to simulate the fluid flow,turbulent kinetic energy distribution,mixing performance and power consumption in a stirred tank.The simulation results were also verified by the water model experiments,and good agreement was achieved.The solid-liquid mixing performances of Intermig impeller with different blade structures were compared in detail.The results show that the improved Intermig impeller not only enhances the solid mixing and suspension,but also saves more than 20% power compared with the standard one.The inner blades have relatively little influence on power and the best angle of inner blades is 45°,while the outer blades affect greatly the power consumption and the optimized value is 45°.
文摘The aim of this paper is to study the effect of agitator’s types on the turbulent flows in stirred tanks without and with baffles.The hydrodynamics behavior induced by four different agitator’s types:a Rushton turbine(RT),a circular blade turbine(CBT),a diverging triangular blade turbine(DTBT)and converging triangular blade turbine(CTBT)are numerically predicted by solving the Navier-Stokes equations and RNGκ–εturbulent model.The simulations are carried out using the Multi Reference Frame(MRF)approach.The numerical results showed good agreement with experiment.We find that the agitator CTBT gives an important profit on the power consumption per report/ratio the others and DTBT give a good reduction of the vortex size of the impeller angles.
基金Supported by the National Basic Research Program of China(2012CB224806)the National Natural Science Foundation of China(21306197,91434126)+1 种基金the Major National Scientific Instrument Development Project(21427814)Jiangsu National Synergetic Innovation Center for Advanced Materials
文摘For the mixing-sensitive reactions, both chemical kinetics and mixing conditions of the reactants determine the distributions of products. The direct quadrature method of moments combining with the interaction by exchange with the mean micro-mixing model(DQMOM-IEM) has been validated for the chemical reacting flows in microreactors. Quite encouraging simulation results offer great promise, but the applicability of this method is needed to be explored furthermore, such as in stirred reactors. In this work, the two-environment DQMOM-IEM model was created with C language and used to customize Fluent through the user-defined functions. The mixing effects on the course of parallel competing chemical reactions carried out in a semi-batch single-phase stirred reactor were predicted. The simulation results show that the rising feed velocity enlarges the volume of reaction zone and maximize the yield of the by-product, which also indicates that the feed stream is more difficultly dispersed into the main stream and the zone surrounding feedpipe exit with high turbulent kinetic dissipation rate cannot be efficiently used.
文摘Liquefied petroleum gas (LPG) is an important urban gas source in China. Before supplied to customers by pipeline supply systems, LPG is stored in tanks in LPG vaporization stations. Designers usually decide the number and the size of storage tanks by their experience during constructions of vaporization stations. These decisions are usually not best and most economical. To solve the problem, a compact mixed integer nonlinear programming model has been developed in this paper. The objective is to minimize annual storage cost of the vaporization station. The model has been transformed into a general nonlinear programming model by transforming integer variables and 0-1 variables into continuous variables. One LPG vaporization station was taken as an example to illustrate the usage of the model. The results show that the optimal storage scheme can be determined accurately and quickly by the model and about 15% of storage cost can be saved every year after optimization.