Method and device for determining average diameter of bubbles in epoxy resin wastewater
By calculating the Froud number, Morton number, Reynolds number and Weber number, combined with the preset correction coefficient, the average diameter of the bubble group during the oxidation process of epoxy resin wastewater is solved, and the problem of difficulty in accurately calculating the bubble size in the prior art is improved, and the accuracy of reactor design and mass transfer efficiency are improved.
Patent Information
- Application Number
- CN202410017100.9
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-01-04
- Publication Date
- 2025-07-11
- Estimated Expiration
- 2044-01-04
AI Technical Summary
The prior art is difficult to accurately determine the size of bubbles during oxidation of epoxy resin wastewater, which makes it impossible to effectively guide the design and optimization of the reactor.
By obtaining the characteristic length, fluid viscosity, fluid density and surface tension of the fluid in the reactor, calculate the Froud number, Morton number, Reynolds number and Weber number, use the preset correction coefficient to determine the average diameter of the bubble group according to the formula.
The accurate calculation of the average diameter of the bubble group during the oxidation process of epoxy resin wastewater is achieved, and the design and optimization of the reactor is guided, and the mass transfer efficiency is improved.
Smart Images

Figure CN120293065A_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the field of environmental protection technologies, and particularly to a method, device and storage medium for determining the average diameter of bubbles in epoxy resin wastewater. Background Art
[0002] Epoxy resin is an important type of thermosetting resin, which has excellent chemical stability and corrosion resistance. It has good adhesion performance with various materials, high mechanical strength, and flexible processing technology. Therefore, it is widely used in various fields such as chemical industry, light industry, machinery, electronics, water conservancy, transportation, automobiles, household appliances, and aerospace. However, a large amount of wastewater is generated during the production process of epoxy resin. Its complex components (mainly unreacted monomer epichlorohydrin, solvent toluene or methyl isobutyl ketone (MIBK), lye, and by-products such as NaCl, solid polymers, glycerol, and isopropanol generated during the reaction process), high salt content (more than 10%), high COD (more than 10000 ppm), and the halogenated organic compounds therein have stable structures, high inhibitory and toxic effects on organisms. It is very difficult for sewage treatment plants to effectively treat this type of wastewater, making the problem of epoxy resin wastewater treatment a bottleneck restricting the development of epoxy resin.
[0003] Based on the above problems, in the existing technology for epoxy resin wastewater treatment solutions, wet oxidation technology is usually adopted. In the presence of high temperature, high pressure, and a catalyst, a free radical chain reaction with hydroxyl radicals as the core oxidant is initiated. Hydroxyl radicals have a very high oxidation potential and can attack organic molecules, degrading them into water and carbon dioxide. By using reaction intensification technology, oxygen-rich is dispersed into micro-nano bubbles, changing the oxygen molecule transport interface scale from the centimeter-millimeter interface to the micron-scale molecular transport, achieving a several-fold or even hundreds-fold increase in the mass transfer rate, significantly enhancing the concentration of hydroxyl radicals and the degradation rate of wastewater molecules, and thus achieving the purpose of oxidizing and removing organic matter and meeting the application requirements of ion-exchange membrane caustic soda plants. However, the parameters of the wet oxidation reaction of liquid epoxy resin wastewater are complex and the action mechanism is non-linear, resulting in the inability to accurately determine the size of bubbles during the oxidation reaction. Summary of the Invention
[0004] The purpose of the embodiments of the present application is to provide a method and device for determining the average diameter of bubbles in epoxy resin wastewater, so as to solve the problem in the existing technology that it is difficult to accurately calculate the size of bubbles during the oxidation process of epoxy resin wastewater.
[0005] To achieve the above purpose, the first aspect of the present application provides a method for determining the average diameter of bubbles in epoxy resin wastewater, and the method includes:
[0006] Obtain the characteristic length in the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid;
[0007] Determine the Froude number inside the reactor based on the superficial gas velocity on the gas surface and the characteristic length;
[0008] Determine the Morton number inside the reactor based on the fluid viscosity, fluid density, and surface tension of the fluid;
[0009] Determine the Reynolds number inside the reactor based on the fluid density, superficial liquid gas velocity, characteristic length, and fluid viscosity;
[0010] Determine the Weber number inside the reactor based on the superficial gas velocity, characteristic length, fluid density, and surface tension of the fluid;
[0011] Based on a preset correction coefficient, determine the average diameter of the bubble swarm inside the reactor according to the Froude number, Morton number, Reynolds number, and Weber number.
[0012] In the embodiment of the present application, determining the average diameter of the bubble swarm inside the reactor based on a preset correction coefficient according to the Froude number, Morton number, Reynolds number, and Weber number includes determining the average diameter according to formula (1):
[0013] d 32 = λ3exp(8.0170 × We 0.02545 × Re -0.02485 )*exp(16.9690 × We -7.5106 × Re -1.4644 × Mo 0.5624 × Fr 8.4796 ); (1)
[0014] Wherein, d 32 is the average diameter of the bubble swarm, λ3 is the preset correction coefficient, Fr is the Froude number, Mo is the Morton number, Re is the Reynolds number, and We is the Weber number.
[0015] In the embodiment of the present application, the method further includes:
[0016] Obtain the diameter of the reactor, the inlet gas volume flow rate, and the inlet liquid volume flow rate;
[0017] Determine the superficial gas velocity inside the reactor according to the diameter and the inlet gas volume flow rate; wherein, the superficial gas velocity satisfies formula (2):
[0018]
[0019] Wherein, U g is the superficial gas velocity, Q G is the inlet gas volume flow rate of the reactor, and D is the diameter of the reactor.
[0020] In the embodiment of the present application, the superficial liquid gas velocity inside the reactor satisfies formula (3):
[0021]
[0022] Among them, U l is the gas superficial velocity at the liquid surface, Q L is the liquid inlet volume flow rate of the reactor, and D is the diameter.
[0023] In the embodiments of the present application, determining the Froude number inside the reactor according to the gas superficial velocity and the characteristic length includes determining the Froude number according to formula (4):
[0024]
[0025] Among them, Fr is the Froude number, and U g is the gas superficial velocity, and d p is the characteristic length.
[0026] In the embodiments of the present application, determining the Morton number inside the reactor according to the fluid viscosity, fluid density, and surface tension of the fluid includes determining the Morton number according to formula (5):
[0027]
[0028] Among them, Mo is the Morton number, and μ l is the liquid viscosity, ρ l is the fluid density, and σ l is the surface tension of the fluid.
[0029] In the embodiments of the present application, determining the Reynolds number inside the reactor according to the fluid density, gas superficial velocity at the liquid surface, characteristic length, and fluid viscosity includes determining the Reynolds number according to formula (6):
[0030]
[0031] Among them, Re is the Reynolds number, and μ l is the liquid viscosity, ρ l is the fluid density, and d p is the characteristic length.
[0032] In the embodiments of the present application, determining the Weber number inside the reactor according to the gas superficial velocity, characteristic length, fluid density, and surface tension of the fluid includes determining the Weber number according to formula (7):
[0033]
[0034] Among them, We is the Weber number, ρ l is the fluid density, d p is the characteristic length, σ l is the surface tension of the fluid, and U g is the gas superficial velocity.
[0035] The second aspect of the present application provides a device for determining the average diameter of bubbles in epoxy resin wastewater. The device includes:
[0036] A memory configured to store instructions; and
[0037] A processor configured to call instructions from the memory and capable of implementing a method for determining the average diameter of bubbles in epoxy resin wastewater when executing the instructions.
[0038] The third aspect of the present application provides a storage medium, on which instructions are stored, and the instructions are used to cause a machine to execute a method for determining the average diameter of bubbles in epoxy resin wastewater.
[0039] Through the above technical solutions, first, the characteristic length in the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid are obtained. The Froude number inside the reactor is determined according to the superficial gas velocity on the gas surface and the characteristic length in the reactor. The Morton number inside the reactor is determined according to the fluid viscosity, the fluid density, and the surface tension of the fluid. The Reynolds number inside the reactor is determined according to the fluid density, the superficial liquid velocity in the reactor, the characteristic length, and the fluid viscosity. The Weber number inside the reactor is determined according to the superficial gas velocity on the gas surface in the reactor, the characteristic length, the fluid density, and the surface tension of the fluid. Based on a preset correction coefficient, the average diameter of the bubble swarm in the reactor is determined according to the Froude number, the Morton number, the Reynolds number, and the Weber number. Through this technical solution, the average diameter of the bubble swarm during the oxidation of epoxy resin wastewater can be accurately calculated, and subsequently, the design and optimization scheme of the reactor can be guided according to the average diameter of the bubble swarm.
[0040] Other features and advantages of the embodiments of the present application will be described in detail in the subsequent specific implementation section. Description of the Drawings
[0041] The drawings are used to provide a further understanding of the embodiments of the present application, and constitute a part of the specification. Together with the following specific implementation, they are used to explain the embodiments of the present application, but do not constitute a limitation to the embodiments of the present application. In the drawings:
[0042] Figure 1 Schematically shows a flowchart of a method for determining the average diameter of bubbles in epoxy resin wastewater according to an embodiment of the present application;
[0043] Figure 2 Schematically shows a schematic diagram of the fitting result of the standard deviation under the liquid epoxy resin wastewater wet oxidation system according to an embodiment of the present application;
[0044] Figure 3 Schematically shows a schematic diagram of the fitting result of the natural logarithm average value under the liquid epoxy resin wastewater wet oxidation system according to an embodiment of the present application;
[0045] Figure 4 A schematic diagram showing the fitting result of the average bubble diameter in the wet air oxidation system for liquid epoxy resin wastewater according to an embodiment of the present application;
[0046] Figure 5 A structural diagram schematically showing a device for determining the average bubble diameter of epoxy resin wastewater according to an embodiment of the present application.
[0047] Description of reference numerals
[0048] 510 Memory 520 Processor Detailed implementation manners
[0049] To make the objectives, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. It should be understood that the specific implementation manners described herein are only used to illustrate and explain the embodiments of the present application, and are not used to limit the embodiments of the present application. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present application without creative efforts shall fall within the scope of protection of the present application.
[0050] It should be noted that if there are directional indications (such as up, down, left, right, front, back,...) involved in the embodiments of the present application, the directional indications are only used to explain the relative positional relationship and movement conditions between components in a specific posture (as shown in the drawings). If the specific posture changes, the directional indications will also change accordingly.
[0051] In addition, if descriptions such as "first" and "second" are involved in the embodiments of the present application, the descriptions of "first", "second", etc. are only for descriptive purposes, and cannot be understood as indicating or implying their relative importance or implicitly indicating the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include at least one such feature. In addition, the technical solutions between various embodiments can be combined with each other, but it must be based on the fact that those of ordinary skill in the art can implement them. When the combination of technical solutions results in contradictions or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present application.
[0052] Figure 1 A flowchart schematically showing a method for determining the average bubble diameter of epoxy resin wastewater according to an embodiment of the present application. As Figure 1 shown, the embodiments of the present application provide a method for determining the average bubble diameter of epoxy resin wastewater, and the method may include the following steps.
[0053] Step 101: Obtain the characteristic length inside the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid;
[0054] Step 102: Determine the Froude number inside the reactor based on the superficial gas velocity at the gas-liquid interface and the characteristic length;
[0055] Step 103: Determine the Morton number inside the reactor based on the fluid viscosity, the fluid density, and the surface tension of the fluid;
[0056] Step 104: Determine the Reynolds number inside the reactor based on the fluid density, the superficial gas velocity at the liquid surface, the characteristic length, and the fluid viscosity;
[0057] Step 105: Determine the Weber number inside the reactor based on the superficial gas velocity at the gas-liquid interface, the characteristic length, the fluid density, and the surface tension of the fluid;
[0058] Step 106: Based on a preset correction factor, determine the average diameter of the bubble swarm inside the reactor according to the Froude number, the Morton number, the Reynolds number, and the Weber number.
[0059] In one embodiment, first obtain the characteristic length inside the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid, which are used to provide the necessary basic data for calculating the Froude number, the Morton number, the Reynolds number, and the Weber number later. In one embodiment, the characteristic length is related to the geometric shape and size of the pipeline and the reactor. For a circular pipe, the characteristic length is usually selected as the pipe diameter (diameter). This is because in a circular pipe, the diameter is an obvious geometric feature that can be used to describe the fluid flow and heat and mass transfer processes. In one embodiment, the fluid viscosity of the fluid is a property that describes the internal resistance of the fluid, that is, the resistance of the fluid to flow. The higher the viscosity, the more difficult it is for the fluid to flow. In the wet air oxidation reaction of epoxy resin wastewater, the fluid viscosity is related to the mixing, mass transfer, and heat transfer processes of the reactants. A fluid with a high viscosity may result in poor mixing performance and affect the mass transfer rate. In one embodiment, the fluid density refers to the mass of the fluid per unit volume and is a property that describes the compressibility and inertia of the fluid. The greater the density, the heavier the fluid. In the wet air oxidation reaction of epoxy resin wastewater, the fluid density is involved in the movement, floating, or sinking of the fluid in the reactor and has an impact on the flow properties of the fluid. In one embodiment, the fluid surface tension is a property that describes the intermolecular force on the fluid surface. The greater the surface tension, the more difficult it is for the surface to be broken. In the wet air oxidation reaction of epoxy resin wastewater, the surface tension may affect the formation, floating, and dispersion of gas bubbles, and thus affect the gas-liquid interaction inside the reactor. A reactor is a device used to control chemical reactions. It is the core component of the chemical reaction process and is used to convert raw materials (reactants) into the desired products. In one embodiment, the epoxy resin wastewater undergoes a wet air oxidation reaction in the reactor.
[0060] In one embodiment, the Froude number inside the reactor is determined based on the gas superficial velocity and the characteristic length. The Froude number refers to a dimensionless parameter in fluid mechanics that characterizes the relative magnitudes of the inertial force and the gravitational force of a fluid, and it represents the ratio of the magnitudes of the inertial force and the gravitational force. The gas superficial velocity refers to the gas velocity at the gas surface inside the reactor. In one embodiment, the Morton number inside the reactor can be determined based on the fluid viscosity, the fluid density, and the surface tension of the fluid. The Morton number is a dimensionless number in fluid mechanics and, together with the Eötvös number, describes the shape of a bubble or a water droplet when it moves in a fluid or a continuous phase. The Reynolds number inside the reactor can be determined based on the fluid density, the liquid superficial velocity, the characteristic length, and the fluid viscosity. The Reynolds number in fluid mechanics is the ratio of the inertial force of a fluid to the viscous force, and it is a dimensionless quantity. The liquid superficial velocity refers to the gas velocity at the liquid surface inside the reactor. The Weber number inside the reactor can be determined based on the gas superficial velocity, the characteristic length, the fluid density, and the surface tension of the fluid. The Weber number represents the ratio of the inertial force and the surface tension effect. The smaller the Weber number, the more important the surface tension. Further, based on a preset correction coefficient, the average diameter of the bubble swarm inside the reactor can be determined according to the Froude number, the Morton number, the Reynolds number, and the Weber number. Among them, the preset correction coefficient is obtained by fitting experimental data.
[0061] In one embodiment, a calculation model for the average diameter of the bubble swarm in the wet oxidation system of liquid epoxy resin wastewater is established by fitting experimental data, and the average diameter of the bubble swarm during the oxidation process of epoxy resin wastewater can be accurately calculated.
[0062] In one embodiment, determining the average diameter of the bubble swarm inside the reactor based on the preset correction coefficient, the Froude number, the Morton number, the Reynolds number, and the Weber number may further include determining the average diameter according to formula (1):
[0063] d 32 =λ3exp(8.0170×We 0.02545 ×Re -0.02485 )*exp(16.9690×We -7.5106 ×Re -1.4644 ×Mo 0.5624 ×Fr 8.4796 );(1)
[0064] Wherein, d 32 is the average diameter of the bubble swarm, λ3 is the preset correction coefficient, Fr is the Froude number, Mo is the Morton number, Re is the Reynolds number, and We is the Weber number.
[0065] In one embodiment, the superscripts of the four numbers, namely the Froude number, Morton number, Reynolds number, and Weber number, and the preset correction coefficient λ3 are all obtained by fitting experimental data. Preferably, the preset correction coefficient λ3 can take values between 1 and 1.5. In one embodiment, the average diameter calculation model of the bubble swarm in the wet oxidation system of epoxy resin wastewater can accurately calculate the average diameter of the bubble swarm during the oxidation process of epoxy resin wastewater through formula (1).
[0066] In one embodiment, the method may further include: obtaining the diameter of the reactor, the inlet gas volume flow rate, and the inlet liquid volume flow rate; determining the gas superficial velocity in the reactor according to the diameter and the inlet gas volume flow rate; wherein, the gas superficial velocity satisfies formula (2):
[0067]
[0068] wherein, U g is the gas superficial velocity, Q G is the inlet gas volume flow rate of the reactor, and D is the diameter of the reactor.
[0069] The liquid superficial velocity in the reactor satisfies formula (3):
[0070]
[0071] wherein, U l is the liquid superficial velocity, Q L is the inlet liquid volume flow rate of the reactor, and D is the diameter of the reactor.
[0072] In one embodiment, it is necessary to calculate the superficial gas velocity and the superficial liquid velocity in the reactor based on the inlet gas volume flow rate of the reactor, the diameter of the reactor, and the inlet liquid volume flow rate of the reactor, which are also important parameters affecting the superficial gas velocity and the superficial liquid velocity in the reactor. Changes in these parameters can affect the generation, floating, and coalescence of bubbles, thereby affecting the average diameter of the bubble swarm. In one embodiment, considering the inlet gas volume flow rate of the reactor, the diameter of the reactor, and the inlet liquid volume flow rate of the reactor, the average diameter of the bubble swarm can be calculated more accurately. In one embodiment, the superficial gas velocity refers to the volume flow rate of gas passing through a unit cross-sectional area and is used to characterize the flow velocity in a pipeline or a reactor. In the wet air oxidation reaction of epoxy resin wastewater, an appropriate superficial gas velocity can affect the size and distribution of bubbles, and thus affect the gas-liquid interaction and the mass transfer performance of the gas. Too small a superficial gas velocity may result in insufficient gas-liquid mixing, while too large a superficial gas velocity may result in overly fine bubbles. In one embodiment, the inlet liquid volume flow rate of the reactor represents the volume of gas entering the reactor per unit time and determines the rate at which gas enters the reactor. In the wet air oxidation reaction of epoxy resin wastewater, the inlet gas volume flow rate directly determines the rate at which gas enters the reactor, thereby affecting the supply rate of reactants. Under aerobic conditions, sufficient oxygen supply is the key to generating oxidation reactions in the wet air oxidation reaction. In one embodiment, the diameter of the reactor refers to the diameter of the cross-section of the reactor, and the size of the reactor directly affects the volume of the reactor and the distribution of gas in the liquid. In the wet air oxidation reaction of epoxy resin wastewater, the diameter of the reactor is crucial for providing sufficient reaction volume to ensure sufficient mixing and gas mass transfer. An appropriate reactor diameter helps to avoid overly large or overly small bubbles, thereby affecting the gas mass transfer performance and the reaction effect. In one embodiment, the Froude number inside the reactor can also be determined based on the superficial gas velocity and the characteristic length, including determining the Froude number according to formula (4):
[0073]
[0074] where Fr is the Froude number, U g is the superficial gas velocity, and d p is the characteristic length.
[0075] In one embodiment, the Morton number inside the reactor can also be determined based on the fluid viscosity, the fluid density, and the surface tension of the fluid, including determining the Morton number according to formula (5):
[0076]
[0077] where Mo is the Morton number, μ l is the liquid viscosity, ρ l is the fluid density, and σ lis the surface tension of the fluid.
[0078] In one embodiment, the Reynolds number inside the reactor can also be determined according to the fluid density, the superficial gas velocity of the liquid surface, the characteristic length, and the fluid viscosity, including determining the Reynolds number according to formula (6):
[0079]
[0080] where Re is the Reynolds number, and μ l is the liquid viscosity, ρ l is the fluid density, d p is the characteristic length.
[0081] In one embodiment, the Weber number inside the reactor can also be determined according to the superficial gas velocity of the gas, the characteristic length, the fluid density, and the surface tension of the fluid, including determining the Weber number according to formula (7):
[0082]
[0083] where We is the Weber number, ρ l is the fluid density, d p is the characteristic length, σ l is the surface tension of the fluid, and U g is the superficial gas velocity of the gas.
[0084] In one embodiment, the magnitude of the Froude number is used to describe the relationship between the inertial force and the gravitational force in the fluid. Specifically, there are three cases of the Froude number: small Froude number, large Froude number, and medium Froude number. When the Froude number is a small Froude number, the gravitational force is relatively large compared to the inertial force. This may result in a smaller average diameter of the bubble swarm because the larger inertial force can counteract the smaller gravitational force, making it easier for the bubbles to float upward. When the Froude number is a large Froude number, the inertial force is relatively large compared to the gravitational force. This may result in a larger average diameter of the bubble swarm because the inertial force is large enough that it is more difficult for the bubbles to float upward, which may lead to a longer residence time of the bubbles in the liquid and affect the processes of bubble coalescence and dispersion. When the Froude number is a medium Froude number, the inertial force and the gravitational force are relatively balanced. This may result in the average diameter of the bubble swarm being affected by both of them and showing a certain equilibrium state. In one embodiment, by adjusting the superficial gas velocity of the gas and the characteristic length, the Froude number can be controlled, thereby affecting the average diameter of the bubble swarm.
[0085] In one embodiment, the magnitude of the Morton number is used to describe the relationship between the inertial force and the surface tension. Specifically, there are three cases of the Morton number: small Morton number, large Morton number, and medium Morton number. When the Morton number is a small Morton number, it indicates that the surface tension is relatively large compared to the inertial force. At this time, the hydrodynamics is mainly dominated by the surface tension, the average diameter of the bubble group may be small, showing a strong surface tension effect, resulting in bubbles being more likely to coalesce to form larger bubbles. When the Morton number is a large Morton number, it indicates that the inertial force is relatively large compared to the surface tension. At this time, the inertial effect dominates the hydrodynamics, the average diameter of the bubble group may be large, showing a strong inertial effect, resulting in bubbles being more difficult to coalesce and possibly forming smaller bubbles. When the Morton number is a medium Morton number, it indicates that the surface tension and the inertial force are relatively balanced. At this time, the influence of both is roughly equivalent, and the average diameter of the bubble group is affected by both the surface tension and the inertial force, and may show a certain equilibrium state. The liquid viscosity, fluid density, and surface tension of the fluid can affect the formation, coalescence, and dispersion processes of the bubble group, thereby affecting the average diameter of the bubbles.
[0086] In one embodiment, the Reynolds number is used to describe the relationship between the inertial force and the viscous force of a fluid. Specifically, there are three cases of the Reynolds number: small Reynolds number, large Reynolds number, and medium Reynolds number. When the Reynolds number is a small Reynolds number, the viscous force is relatively large compared to the inertial force. This may result in a smaller average diameter of the bubble group because the viscous force can inhibit the coalescence and dispersion of the bubbles. When the Reynolds number is a large Reynolds number, the inertial force is relatively large compared to the viscous force. This may result in a larger average diameter of the bubble group because the inertial force makes it easier for the bubbles to disperse and coalesce. When the Reynolds number is a medium Reynolds number, the inertial force and the viscous force are relatively balanced. This may result in the average diameter of the bubble group being affected by both of them, showing a certain equilibrium state.
[0087] In one embodiment, the Weber number is used to describe the relationship between the surface tension, inertia, and gravity in a liquid fluid. Specifically, the Weber number is divided into three cases: small Weber number, large Weber number, and medium Weber number. When the Weber number is a small Weber number, the surface tension is relatively large compared to the inertial force. This may result in a smaller average diameter of the bubble group because the surface tension promotes the easier formation and dispersion of the bubbles. When the Weber number is a large Weber number, the inertial force is relatively large compared to the surface tension. This may result in a larger average diameter of the bubble group because the inertial force makes it more difficult for the bubbles to form and be maintained. When the Weber number is a medium Weber number, the surface tension and the inertial force are relatively balanced. This may result in the average diameter of the bubble group being affected by both of them, showing a certain equilibrium state. By adjusting the gas superficial velocity, characteristic length, fluid density, and surface tension, the Weber number can be controlled, thereby affecting the average diameter of the bubble group.
[0088] In one embodiment, the effects of parameters such as the gas superficial velocity at the surface, characteristic length, fluid viscosity, fluid density, surface tension of the fluid, fluid density, and liquid superficial gas velocity on the average diameter of the bubble group are considered, thereby making the calculation of the average diameter of the bubble group more accurate.
[0089] The following takes a specific embodiment to illustrate the results of the above method for determining the average diameter of epoxy resin wastewater bubbles.
[0090] First, it is necessary to record the physical properties and operating parameters of the reaction process. Record the physical properties of liquid epoxy resin wastewater at 298.15K and 1atm. Among them, in a specific embodiment, the physical properties of liquid epoxy resin wastewater at 298.15K and 1atm are shown in Table 1:
[0091] Table 1
[0092] <![CDATA[σ L (N / m)]]> <![CDATA[ρ L (kg / m 3 )]]> <![CDATA[g (m / s 2 )]]> <![CDATA[μ L (Pa·s)]]> 0.07436 999.7 9.8 0.0018077 0.07362 999.1 9.8 0.0009404 0.07288 998.2 9.8 0.0010051 0.07214 997.0 9.8 0.0008937 0.07140 995.6 9.8 0.0008007
[0093] Record the reactor parameters and operating parameters of the wet oxidation reaction of this epoxy resin wastewater. Among them, the reactor parameters and operating parameters are shown in Table 2:
[0094] Table 2
[0095] D(m) <![CDATA[Q G (m 3 / s)]]> <![CDATA[Q L (m 3 / s)]]> 1.8 3.924 0.127 1.8 3.924 0.254 1.8 3.924 0.509 1.8 3.924 1.02 1.8 3.924 1.27
[0096] Substitute the data obtained from the above records into the following formula to conduct a fitting data experiment.
[0097] (A) Establish an equation for the generalized average bubble diameter:
[0098] For a group of bubble groups, let d, σ, and m be the diameter, standard deviation, and natural logarithm average value of the bubbles respectively. Then the probability density function of the bubble diameter d satisfies formula (8):
[0099]
[0100] The mathematical expectation (i.e., the mathematical average value) and variance of the bubble diameter respectively satisfy formula (9) and formula (10):
[0101]
[0102]
[0103] Kiyomi Akita deduced that the bubble diameter satisfies the above natural logarithm distribution formula, and the calculation formula for the average bubble diameter d 32 satisfies formula (11):
[0104] d 32 = exp(m + 2.5σ 2)); (11)
[0105] By combining Equation (9) and Equation (11), the equation for the generalized bubble mean diameter, Equation (12), can be obtained:
[0106]
[0107] (B) Establish the correlation equation for the standard deviation to satisfy Equation (13):
[0108] σ = λ1We -3.6053 × Re -0.7322 × Mo 0.2812 × Fr 4.2398 ; (13)
[0109] where σ is the standard deviation, λ1 is the correction coefficient, and its value range is between 2 and 3.
[0110] The Froude number, Morton number, Reynolds number, and Weber number satisfy Equations (4), (5), (6), and (7):
[0111]
[0112]
[0113]
[0114]
[0115] where U g is the gas superficial velocity at the gas-liquid interface, with the unit of m / s; U l is the liquid superficial velocity at the gas-liquid interface, with the unit of m / s; d p is the characteristic length, with the unit of m; g is the acceleration due to gravity, with the unit of m / s 2 ; μ L is the viscosity of the fluid, with the unit of Pa·s; ρ L is the density of the fluid, with the unit of kg / m 3 ; σ L is the surface tension of the fluid, with the unit of N / m.
[0116] The gas superficial velocity satisfies Equation (2):
[0117]
[0118] The liquid superficial velocity satisfies Equation (3):
[0119]
[0120] where the reactor diameter is D, with the unit of m; the inlet gas volume flow rate is Q G , m 3 / s; The liquid inflow volume flow rate is Q L , in the unit of m 3 / s.
[0121] (C) Establish the correlation equation for the natural logarithm average value to satisfy formula (14):
[0122] m = λ2We 0.02545 ×Re -0.02485 ; (14)
[0123] Where λ2 is a correction coefficient, and its value range is between 7 and 9.
[0124] (D) Establish the calculation model for the average bubble diameter in the wet oxidation system of liquid epoxy resin wastewater:
[0125] According to formulas (11), (12), (13) and (14), the calculation model for the average bubble diameter satisfies formula (1):
[0126] d 32 = λ3exp(8.0170×We 0.02545 ×Re -0.02485 )*exp(16.9690×
[0127] We - 7.5106×Re - 1.4644×Mo0.5624×Fr8.4796; (1)
[0128] Where λ3 is a correction coefficient, and its value is between 1 and 1.5.
[0129] (E) Establish the gas-liquid phase interface area model in the wet oxidation system of liquid epoxy resin wastewater to satisfy formula (15):
[0130] a = 6ε G / d 32 (15)
[0131] Where a is the gas-liquid phase interface area, in the unit of m 2 / m 3 ; ε G is the gas holdup obtained from the bubble size distribution, in the unit of m 3 / m 3 ; d 32 is the Sauter mean diameter of the bubble swarm, in the unit of m.
[0132] Conduct fitting tests according to the above formulas to obtain the following test results.
[0133] Record the Froude number, Morton number, Reynolds number and Weber number obtained after the experiment ends.
[0134] The dimensionless criteria obtained from the experiment are shown in Table 3:
[0135] Table 3
[0136] Fr Mo Re We 0.21259 2.55E-10 1023.97 0.040332 1.02041 1.92E-11 1639.31 0.135710 4.53515 2.59E-11 1379.30 0.493073 18.14059 1.67E-11 1549.22 1.990130 28.34467 1.11E-11 1726.73 3.137394
[0137] Record the prediction results of the interface of microbubbles with different diameters and the bubble average diameter model.
[0138] The prediction results of the bubble average diameter model are shown in Table 4:
[0139] Table 4
[0140] m σ <![CDATA[d 32 (Calculation) (μm)]]> <![CDATA[d 32 (Experiment) (μm)]]> Error (%) 6.2 0.02 726.06 660.05 9.1 6.3 0.04 768.25 731.66 4.7 6.5 0.06 907.12 898.14 1 6.6 0.05 1078.96 989.87 8.2 6.8 0.08 1196.46 1220.88 2
[0141] The prediction results of the mass transfer flux are shown in Table 5:
[0142] Table 5
[0143] <![CDATA[k GA (m 2 / s)]]> <![CDATA[k LA (m 2 / s)]]> <![CDATA[N A (mol / m 2 / s)]]> 0.01124 0.00169 0.016838 0.01063 0.00170 0.016919 0.00903 0.00159 0.015822 0.00762 0.00146 0.014552 0.00690 0.00140 0.013906
[0144] Through the above technical solution, first establish the generalized average diameter equation of the bubble swarm to obtain the relationship between the bubble average diameter and the natural logarithm average value and the standard deviation, then establish the correlation equation of the natural logarithm average value and the correlation equation of the standard deviation respectively, and finally obtain the calculation model of the bubble average diameter under the wet oxidation system of liquid epoxy resin wastewater, so that the average diameter of the bubble swarm under the wet oxidation system of liquid epoxy resin wastewater can be calculated more accurately. Furthermore, by adjusting the reactor diameter according to the calculation results of the gas holdup and the phase boundary area, it can also guide the design and optimization of the reactor for the wet oxidation system of liquid epoxy resin wastewater.
[0145] Figure 2 Schematically shows the fitting result of the standard deviation under the wet oxidation system of liquid epoxy resin wastewater according to an embodiment of the present application. As Figure 2 shown, in one embodiment, according to the recorded data, a fitting schematic diagram of the standard deviation under the wet oxidation system of liquid epoxy resin wastewater is obtained. The abscissa is the experimental value, and the ordinate is the calculated value of the fitting formula. The experimental value fluctuates within 20% above and below the fitting formula, indicating that the fitting formula of the standard deviation is relatively accurate. In the figure, the upper line represents the boundary line floating 20% above the theoretical value, the middle line represents the theoretical value of the standard deviation under the wet oxidation system of liquid epoxy resin wastewater, and the lower line represents the boundary line floating 20% below the theoretical value. Many dots in the figure represent the error values between the experimental results and the theoretical values.
[0146] Figure 3 Schematically shows the fitting result of the natural logarithm average value under the wet oxidation system of liquid epoxy resin wastewater according to an embodiment of the present application. As Figure 3As shown, in one embodiment, a fitting schematic diagram of the natural logarithm average value of the liquid epoxy resin wastewater wet oxidation system is obtained based on the recorded data. The abscissa is the experimental value, and the ordinate is the calculated value of the fitting formula. The experimental value fluctuates within 5% above and below the fitting formula, indicating that the fitting formula of the natural logarithm average value is relatively accurate. In the figure, the upper line represents the boundary line with a 5% upward float from the theoretical value, the middle line represents the theoretical value of the natural logarithm average value in the liquid epoxy resin wastewater wet oxidation system, and the lower line represents the boundary line with a 5% downward float from the theoretical value. Many dots in the figure represent the error values between the experimental results and the theoretical values.
[0147] Figure 4 Schematically shows the fitting result of the average bubble diameter in the liquid epoxy resin wastewater wet oxidation system according to an embodiment of the present application. As Figure 4 shown, in one embodiment, the abscissa is the experimental value (unit: μm), and the ordinate is the calculated value of the fitting formula (unit: μm). Further research and regression are carried out on the average bubble diameter model. As the working conditions change, the error calculated by the model is always within 20%, that is, the result shows that the accuracy of the prediction model has been improved to a certain extent compared with the prior art. At the same time, according to the calculation of the gas holdup and the phase boundary area, it further guides the design and optimization of the liquid epoxy resin wastewater wet oxidation reactor, and solves the pain points in the reactor design. In the figure, the upper line represents the boundary line with a 20% upward float from the theoretical value, the middle line represents the theoretical value of the average bubble diameter in the liquid epoxy resin wastewater wet oxidation system, and the lower line represents the boundary line with a 20% downward float from the theoretical value. Many dots in the figure represent the error values between the experimental results and the theoretical values.
[0148] Figure 5 Schematically shows the structural diagram of a device for determining the average bubble diameter of epoxy resin wastewater according to an embodiment of the present application. As Figure 5 shown, an embodiment of the present application provides a device for determining the average bubble diameter of epoxy resin wastewater, which may include:
[0149] A memory 510, configured to store instructions; and
[0150] A processor 520, configured to call instructions from the memory 510 and be able to implement the above method for determining the average bubble diameter of epoxy resin wastewater when executing the instructions.
[0151] Specifically, in one embodiment, the processor 520 may be configured to: obtain the characteristic length in the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid; determine the Froude number inside the reactor according to the superficial gas velocity on the gas surface and the characteristic length; determine the Morton number inside the reactor according to the fluid viscosity, the fluid density, and the surface tension of the fluid; determine the Reynolds number inside the reactor according to the fluid density, the superficial liquid velocity, the characteristic length, and the fluid viscosity; determine the Weber number inside the reactor according to the superficial gas velocity, the characteristic length, the fluid density, and the surface tension of the fluid; and determine the average diameter of the bubble swarm in the reactor based on a preset correction coefficient according to the Froude number, the Morton number, the Reynolds number, and the Weber number.
[0152] Further, the processor 520 may also be configured to:
[0153] Determining the average diameter of the bubble swarm in the reactor based on a preset correction coefficient according to the Froude number, the Morton number, the Reynolds number, and the Weber number includes determining the average diameter according to formula (1):
[0154] d 32 = λ3 exp(8.0170 × We 0.02545 × Re -0.02485 ) *
[0155] exp(16.9690 × We -7.5106 × Re -1.4644 × Mo 0.5624 × Fr 8.4796 ); (1)
[0156] wherein, d 32 is the average diameter of the bubble swarm, λ3 is the preset correction coefficient, Fr is the Froude number, Mo is the Morton number, Re is the Reynolds number, and We is the Weber number.
[0157] Further, the processor 520 may also be configured to: obtain the diameter of the reactor, the inlet gas volume flow rate, and the inlet liquid volume flow rate; determine the superficial gas velocity in the reactor according to the diameter and the inlet gas volume flow rate; wherein, the superficial gas velocity satisfies formula (2):
[0158]
[0159] wherein, U g is the superficial gas velocity, Q G is the inlet gas volume flow rate of the reactor, and D is the diameter of the reactor.
[0160] Further, the processor 520 may also be configured to:
[0161] The superficial liquid velocity in the reactor satisfies formula (3):
[0162]
[0163] Among them, U l is the gas superficial velocity at the liquid surface, Q L is the liquid inlet volume flow rate of the reactor, and D is the diameter.
[0164] Furthermore, the processor 520 can also be configured to:
[0165] Determining the Froude number inside the reactor according to the gas superficial velocity and the characteristic length includes determining the Froude number according to formula (4):
[0166]
[0167] Among them, Fr is the Froude number, U g is the gas superficial velocity, d p is the characteristic length.
[0168] Furthermore, the processor 520 can also be configured to:
[0169] Determining the Morton number inside the reactor according to the fluid viscosity, fluid density, and surface tension of the fluid includes determining the Morton number according to formula (5):
[0170]
[0171] Among them, Mo is the Morton number, μ l is the liquid viscosity, ρ l is the fluid density, σ l is the surface tension of the fluid.
[0172] Furthermore, the processor 520 can also be configured to:
[0173] Determining the Reynolds number inside the reactor according to the fluid density, gas superficial velocity at the liquid surface, characteristic length, and fluid viscosity includes determining the Reynolds number according to formula (6):
[0174]
[0175] Among them, Re is the Reynolds number, μ l is the liquid viscosity, ρ l is the fluid density, d p is the characteristic length.
[0176] Furthermore, the processor 520 can also be configured to:
[0177] Determining the Weber number inside the reactor according to the gas superficial velocity, characteristic length, fluid density, and surface tension of the fluid includes determining the Weber number according to formula (7):
[0178]
[0179] Among them, We is the Weber number, ρ l is the fluid density, d p is the characteristic length, σ l is the surface tension of the fluid, U g is the gas superficial velocity at the gas surface.
[0180] Through the above technical solution, first, the processor 520 obtains the characteristic length in the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid. Secondly, the processor 520 determines the Froude number inside the reactor according to the gas superficial velocity and the characteristic length in the reactor. Then, the processor 520 determines the Morton number inside the reactor according to the fluid viscosity, the fluid density, and the surface tension of the fluid. Furthermore, the processor 520 determines the Reynolds number inside the reactor according to the fluid density, the liquid superficial velocity in the reactor, the characteristic length, and the fluid viscosity. After that, the processor 520 determines the Weber number inside the reactor according to the gas superficial velocity, the characteristic length, the fluid density, and the surface tension in the reactor. Finally, through this technical solution, the processor 520 can accurately calculate the average diameter of the bubble swarm during the oxidation process of epoxy resin wastewater, and subsequently, the design and optimization scheme of the reactor can be guided according to the average diameter of the bubble swarm.
[0181] The embodiment of the present application also provides a machine-readable storage medium, and instructions are stored on the machine-readable storage medium, and the instructions are used to cause the machine to execute the above method for determining the average diameter of epoxy resin wastewater bubbles.
[0182] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0183] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, and the combination of processes and / or blocks in the flowchart and / or block diagram can be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate for realizing in the process Figure 1 one process or multiple processes and / or blocksFigure 1 means for the functions specified in one or more boxes.
[0184] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction means that implements the functions specified in one Figure 1 process or more processes and / or boxes Figure 1 means for the functions specified in one or more boxes.
[0185] These computer program instructions may also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are performed on the computer or other programmable device to produce a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one Figure 1 process or more processes and / or boxes Figure 1 means for the functions specified in one or more boxes.
[0186] In a typical configuration, a computing device includes one or more processors (CPUs), an input / output interface, a network interface, and memory.
[0187] The memory may include non-permanent memory in the form of computer-readable media, random access memory (RAM) and / or non-volatile memory such as read-only memory (ROM) or flash memory (flash RAM). The memory is an example of computer-readable media.
[0188] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology for information storage. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette tapes, magnetic disk storage or other magnetic storage devices, or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory computer-readable media such as modulated data signals and carrier waves.
[0189] It should also be noted that the term "comprise", "include" or any other variant thereof is intended to cover non-exclusive inclusion, such that a process, method, commodity or device comprising a series of elements not only includes those elements but also includes other elements not expressly listed, or further includes elements inherent to such process, method, commodity or device. Without further limitation, an element defined by the statement "comprising an..." does not exclude the existence of additional identical elements in the process, method, commodity or device comprising the element.
[0190] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, various modifications and changes can be made to the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included within the scope of the claims of the present application.
Claims
1. A method for determining the average diameter of bubbles in epoxy resin wastewater, characterized in that The method includes: obtaining the characteristic length in the reactor, the fluid viscosity of the fluid, the fluid density, and the surface tension of the fluid; determining the Froude number inside the reactor according to the superficial gas velocity at the gas-liquid interface and the characteristic length; determining the Morton number inside the reactor according to the fluid viscosity, the fluid density, and the surface tension of the fluid; determining the Reynolds number inside the reactor according to the fluid density, the superficial liquid velocity, the characteristic length, and the fluid viscosity; determining the Weber number inside the reactor according to the superficial gas velocity at the gas-liquid interface, the characteristic length, the fluid density, and the surface tension of the fluid; determining the average diameter of the bubble swarm in the reactor based on a preset correction coefficient according to the Froude number, the Morton number, the Reynolds number, and the Weber number.
2. The method according to claim 1, wherein The step of determining the average diameter of the bubble swarm in the reactor based on a preset correction coefficient according to the Froude number, the Morton number, the Reynolds number, and the Weber number includes determining the average diameter according to formula (1): d 32 = λ3exp(8.0170 × We 0.02545 × Re -0.02485 ) * exp(16.9690 × We -7.5106 × Re -1.4644 × Mo 0.5624 × Fr 8.4796 ); (1) where d 32 is the average diameter of the bubble swarm, λ3 is the preset correction coefficient, Fr is the Froude number, Mo is the Morton number, Re is the Reynolds number, and We is the Weber number.
3. The method according to claim 1, wherein The method further includes: obtaining the diameter of the reactor, the inlet gas volume flow rate, and the inlet liquid volume flow rate; determining the superficial gas velocity at the gas-liquid interface in the reactor according to the diameter and the inlet gas volume flow rate; wherein, the superficial gas velocity satisfies formula (2): Among them, U g is the superficial gas velocity on the gas surface, Q G is the inlet gas volume flow rate of the reactor, and D is the diameter of the reactor.
4. The method according to claim 3, characterized in that, The superficial liquid velocity at the gas-liquid interface in the reactor satisfies formula (3): Among them, U l is the superficial gas velocity of the liquid surface, Q L is the liquid inlet volume flow rate of the reactor, and D is the diameter.
5. The method according to claim 1, wherein The step of determining the Froude number inside the reactor according to the superficial gas velocity at the gas-liquid interface and the characteristic length includes determining the Froude number according to formula (4): where Fr is the Froude number, U g is the superficial gas velocity, d p is the characteristic length.
6. The method according to claim 1, characterized in that, The step of determining the Morton number inside the reactor according to the fluid viscosity, the fluid density, and the surface tension of the fluid includes determining the Morton number according to formula (5): where Mo is the Morton number, μ l is the liquid viscosity, ρ l is the fluid density, σ l is the surface tension of the fluid.
7. The method according to claim 1, characterized in that, The step of determining the Reynolds number inside the reactor according to the fluid density, the superficial liquid velocity, the characteristic length, and the fluid viscosity includes determining the Reynolds number according to formula (6): where Re is the Reynolds number, μ l is the liquid viscosity, ρ l is the fluid density, d p is the characteristic length.
8. The method according to claim 1, characterized in that, The step of determining the Weber number inside the reactor according to the superficial gas velocity at the gas-liquid interface, the characteristic length, the fluid density, and the surface tension of the fluid includes determining the Weber number according to formula (7): where, We is the Weber number, ρ l is the fluid density, d p is the characteristic length, σ l is the surface tension of the fluid, U g is the superficial gas velocity at the gas surface.
9. An apparatus for determining the average diameter of bubbles in epoxy resin wastewater, characterized in that, including: a memory configured to store instructions; and a processor configured to call the instructions from the memory and, when executing the instructions, be capable of implementing the method for determining the average diameter of epoxy resin wastewater bubbles according to any one of claims 1 to 8.
10. A machine-readable storage medium, characterized in that, Instructions are stored on the machine-readable storage medium, and the instructions are used to cause the machine to execute the method for determining the average diameter of epoxy resin wastewater bubbles according to any one of claims 1 to 8.
Citation Information
Patent Citations
Method for collecting underwater bubbles and method for monitoring underwater bubble flux under hydrodynamic conditions
CN102608273A
Gas and liquid two-phase flow bubble diameter detection method
CN103591901A
Micro-interface enhanced reactor bubble scale structure-activity regulation and control model building method
CN107563051A
Method and device for determining bubble floating speed in gas-liquid separator
CN112149238A
Gas jet optimization method of gas-liquid bubbling fluidized bed
CN114950284A