Pulse extraction column process simulation method, device, verification method and system

Abstract

The present invention discloses a pulse extraction column process simulation method, device, verification method and system. The simulation method comprises: obtaining equipment attribute information, initialization information and extraction information of the pulse extraction column, wherein the equipment attribute information comprises the number of plates and equipment parameters. According to the number of plates of the pulse extraction column, a plurality of single plate models are established. The plurality of single plate models are combined to obtain an extraction calculation model of the pulse extraction column. According to the equipment parameters, initialization information and extraction information, the extraction calculation model is solved to obtain operation simulation data of the pulse extraction column. The simulation method can simulate the process of the pulse extraction column to predict the actual production status of the pulse extraction column and reduce the prediction cost.

Description

Technical Field

[0001] The invention belongs to the technical field of nuclear industry, and specifically relates to a pulse extraction column process simulation method and device, and a pulse extraction column design verification method and system. Background Art

[0002] The pulse extraction column model is a common extraction equipment in the post-processing production process of nuclear power plants. In the post-processing production process, production personnel need to understand and grasp the impact of changes in different production parameters on the extraction equipment and upstream and downstream process parameters of production, so as to select a more suitable and efficient extraction solution. Therefore, before the actual extraction, it is necessary to build a pulse extraction column model dedicated to the post-processing process to intuitively display and predict the dynamic changes of the pulse extraction column at different times in the production process, so as to provide reliable design verification and technical support for the production design of the post-processing plant.

[0003] However, the materials of each batch of post-processing are different, and the structure of the corresponding pulse extraction column is also different. If a physical pulse extraction column model is built before each extraction to verify the extraction effect, the cost will be greatly increased. Therefore, it is urgent to propose a method that can simulate and analyze the process of pulse extraction column in post-processing. Summary of the invention

[0004] The technical problem to be solved by the present invention is to provide a pulse extraction column process simulation method, device, verification method and system in view of the above-mentioned deficiencies in the prior art. The simulation method can simulate the process of the pulse extraction column to predict the actual production status of the pulse extraction column and reduce the prediction cost.

[0005] According to an embodiment of the first aspect of the present invention, a pulse extraction column process simulation method is provided, comprising:

[0006] S1: Acquire device attribute information, initialization information and extraction information of the pulse extraction column, wherein the device attribute information includes the number of plates and device parameters.

[0007] S2: Establishing a plurality of single plate models according to the number of plates of the pulse extraction column.

[0008] S3: combining the plurality of single plate models to obtain an extraction calculation model of the pulse extraction column.

[0009] S4: solving the extraction calculation model according to the equipment parameters, initialization information and the extraction information to obtain operation simulation data of the pulse extraction column.

[0010] Preferably, the equipment parameters include: plate segment diameter, plate segment height, thickness, height from the ground, volume, upper and lower limits of liquid level measurement, diameter of the upper expansion section, height of the upper expansion section, diameter of the lower expansion section, height of the lower expansion section, inner diameter of the baffle disk, inner diameter of the ring; the initialization information includes: initial pressure, initial liquid level, initial component information; the extraction information includes: equilibrium constant, mass transfer coefficient, dispersed phase selection information, dispersed phase slip velocity.

[0011] Preferably, step S2 specifically includes: according to the number of plates of the pulse extraction column, splitting the pulse extraction column into multiple single plate areas, wherein the single plate area is the area between two adjacent plates; respectively establishing the mass conservation equation, energy conservation equation, mass transfer equation, phase equilibrium equation, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation for each single plate area; and combining the above equations to obtain a single plate model for each single plate area.

[0012] Preferably, the mass transfer equation includes: a mass transfer rate equation and a total volume mass transfer coefficient correlation equation; wherein the expression of the mass transfer rate equation is:

[0013] N=K d α(yy * )Sh

[0014] N=K c α(x * -x)Sh

[0015] Where: N is the total mass transfer rate, mol / s; K d is the total mass transfer coefficient expressed by the dispersed phase solute concentration, m / s; K c is the total mass transfer coefficient expressed as the concentration of the continuous phase solute, m / s; (yy * ) is the mass transfer driving force expressed as the dispersed phase solute concentration, mol / m 3 ;(x * -x) is the mass transfer driving force expressed as the concentration of the continuous phase solute, mol / m 3 ; α is the contact surface area of ​​the two phases per unit volume, m 2 / m 3 ; S is the cross-sectional area of ​​the pulse extraction column, m 2 ; h is the total height of the liquid level in the single plate area of ​​the pulse extraction column, m;

[0016] The phase equilibrium equation adopts Richardson semi-empirical distribution ratio model;

[0017] The extraction influencing factor correlation formula includes: a dispersed phase droplet diameter correlation formula, a dispersed phase retention fraction correlation formula, and an axial mixing coefficient correlation formula.

[0018] Preferably, the step S4 specifically includes the following steps: initializing the extraction calculation model according to the initialization information to obtain a static device model; substituting the device parameters and the extraction information into the static device model, performing iterative calculations on the static device model, and obtaining operation simulation data of the pulse extraction column.

[0019] Preferably, the operation simulation data includes: liquid level height, boundary data of liquid surface boundary, temperature data, pressure data, upper liquid volume, lower liquid volume, component ratio data of logistics, and flow data of logistics.

[0020] According to an embodiment of the second aspect of the present invention, a design verification method for a pulse extraction column is provided, comprising: calculating operation simulation data of the pulse extraction column according to the above-mentioned pulse extraction column process simulation method; and judging whether the pulse extraction column meets the design requirements according to the operation simulation data.

[0021] Preferably, the operation simulation data includes component ratio data of the logistics in the pulse extraction column; and judging whether the pulse extraction column meets the design requirements based on the operation simulation data specifically includes: when the component ratio data in the pulse extraction column is within a preset component ratio range, judging that the pulse extraction column meets the design requirements; otherwise, judging that the pulse extraction column does not meet the design requirements.

[0022] According to an embodiment of the third aspect of the present invention, there is provided a pulse extraction column process simulation device, comprising: an acquisition module, a processing module, a control module and a calculation module; the acquisition module is used to acquire equipment attribute information, initialization information and extraction information of the pulse extraction column, the equipment attribute information including the number of plates and equipment parameters; the processing module is electrically connected to the acquisition module, and is used to establish a plurality of single-plate models according to the number of plates of the pulse extraction column; the control module is electrically connected to the processing module, and is used to combine a plurality of the single-plate models to obtain an extraction calculation model of the pulse extraction column; the calculation module is electrically connected to the acquisition module and the control module, and is used to solve the extraction calculation model according to the equipment parameters, initialization information and the extraction information to obtain operation simulation data of the pulse extraction column.

[0023] Preferably, the processing module includes a first processing unit, a second processing unit and a third processing unit; the first processing unit is electrically connected to the acquisition module, and is used to split the pulse extraction column into multiple single-plate areas according to the number of plates of the pulse extraction column, and the single-plate area is the area between two adjacent upper and lower plates; the second processing unit is electrically connected to the first processing unit, and is used to establish the mass conservation equation, energy conservation equation, mass transfer equation, phase equilibrium equation, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation for each single-plate area respectively; the third processing unit is electrically connected to the second processing unit, and combines the above equations to obtain a single-plate model for each single-plate area.

[0024] Preferably, the calculation module includes a first calculation unit and a second calculation unit; the first calculation unit is electrically connected to the acquisition module and the control module, and is used to initialize the extraction calculation model according to the initialization information to obtain a static device model; the second calculation unit is electrically connected to the first calculation unit and the acquisition module, and is used to substitute the device parameters and the extraction information into the static device model, perform iterative calculations on the static device model, and obtain operation simulation data of the pulse extraction column.

[0025] According to an embodiment of the fourth aspect of the present invention, a design verification system for a pulse extraction column is provided, comprising a judgment module and the above-mentioned pulse extraction column process simulation device, wherein the pulse extraction column process simulation device is used to calculate operation simulation data of the pulse extraction column; the judgment module is electrically connected to the pulse extraction column process simulation device, and is used to judge whether the pulse extraction column meets the design requirements based on the operation simulation data.

[0026] The pulse extraction column process simulation method of the present invention establishes a single plate model between adjacent plates of the pulse extraction column, and combines multiple single plate models to obtain an extraction calculation model. The extraction calculation model is solved through the equipment attribute information, initialization information and extraction information of the pulse extraction column to obtain the operation simulation data of the pulse extraction column, and then the actual production status of the pulse extraction column can be predicted based on the operation simulation data. In addition, this method predicts by mathematical model calculation, without the need to establish a physical model of the pulse extraction column, which can reduce the prediction cost.

[0027] In addition, this method can be applied to the dynamic simulation process monitoring and prediction of the post-processing plant process, thereby realizing the design verification of the extraction production plan of a large post-processing plant. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Figure 1 is a schematic diagram of the module construction of the pulse extraction column in some embodiments of the present invention;

[0029] Figure 2 is a schematic diagram of an input interface for basic parameters of a pulse extraction column in some embodiments of the present invention;

[0030] Figure 3 is a schematic diagram of an input interface of equipment parameters of a pulse extraction column in some embodiments of the present invention;

[0031] Figure 4 is a schematic diagram of an initial property setting interface of a pulse extraction column in some embodiments of the present invention;

[0032] Figure 5 is a schematic diagram of a parameter configuration interface of a pulse extraction column in some embodiments of the present invention;

[0033] Figure 6 is a schematic diagram of a reaction configuration interface of a pulse extraction column in some embodiments of the present invention;

[0034] Figure 7 It is a schematic diagram of the simulation result output interface of the pulse extraction column in some embodiments of the present invention. DETAILED DESCRIPTION

[0035] The following will be combined with the accompanying drawings in the present invention to clearly and completely describe the technical solutions in the invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of the present invention.

[0036] In the description of the present invention, the terms “first” and “second” are used for descriptive purposes only and should not be understood as indicating or implying relative importance.

[0037] In the description of the present invention, each unit and module involved may correspond to only one physical structure, or may be composed of multiple physical structures, or multiple units and modules may be integrated into one physical structure; the units and modules involved may be implemented by software or by hardware, for example, the units and modules may be located in a processor.

[0038] In the description of the present invention, the functions and steps marked in the flowcharts and block diagrams of the present invention may occur in an order different from that marked in the drawings if there is no conflict.

[0039] Example 1

[0040] See also Figure 1 The present invention discloses a pulse extraction column process simulation method, comprising the following steps:

[0041] S1: Obtaining device attribute information, initialization information and extraction information of the pulse extraction column, where the device attribute information includes the number of plates and device parameters;

[0042] S2: Establish multiple single-plate models according to the number of plates of the pulse extraction column;

[0043] S3: combining multiple single-plate models to obtain an extraction calculation model of a pulse extraction column;

[0044] S4: Solve the extraction calculation model according to the equipment parameters, initialization information and extraction information to obtain the operation simulation data of the pulse extraction column.

[0045] Before explaining the pulse extraction column process simulation method in this embodiment, it is necessary to explain the structure and working principle of the pulse extraction column: the pulse extraction column is a column-type extraction equipment with energy input. By inputting pulse energy, the fluid in the column can be pulsated rapidly back and forth, which can not only crush the droplets and increase the retained fraction of the dispersed phase, thereby greatly increasing the contact area between the two phases, but also increase the turbulence of the phase interface and promote the mass transfer process.

[0046] like Figure 1 As shown, Figure 1 The simplified model of the pulse extraction column is shown. It can be seen that the column body of the pulse extraction column is provided with multiple plates, which are arranged at intervals along the axial direction of the main body. The upper part of the pulse extraction column is provided with a heavy phase inlet, and the lower part is provided with a light phase inlet. The simplified model includes information such as the number of plates, the position of the heavy phase inlet, the position of the light phase inlet, the position of the heavy phase outlet, and the position of the light phase outlet of the pulse extraction column. By simplifying the structure of the pulse extraction column, simulation calculations can be performed more conveniently.

[0047] During the operation of the pulse extraction column, liquid continuously enters the column (the heavy phase enters through the heavy phase inlet, and the light phase enters through the light phase inlet). In the column, in addition to the countercurrent contact of the two phases due to the density difference, the dispersed phase is better broken by inputting pulse energy, which increases the mass transfer specific surface area of ​​the two phases, thereby improving the mass transfer efficiency of the two phases. During the operation, the pulse frequency and amplitude are adjusted by the pulse device and the pulse speed controller to improve the extraction efficiency of the pulse extraction column. Since the pulse extraction column has the advantages of high extraction efficiency and small footprint, it has been widely used in the field of spent fuel reprocessing.

[0048] In this embodiment, the pulse extraction column process simulation method is a simulation algorithm of a pulse extraction column device based on the Purex process. The Purex process is a spent fuel post-processing process, which mainly uses tributyl phosphate as an extractant to separate and recover uranium and plutonium from spent fuel. This simulation method can realize the simulation calculation of the pulse extraction column in the dynamic mechanism simulation software of Zhejiang Zhongkong Intelligent Simulation Platform (VxOTS).

[0049] According to the Pruex process, the pulse extraction tower equipment principle is as follows Figure 2 As shown, Figure 2 The extraction column shown in the figure is a baffle pulse extraction column. Under the action of the pulse, the dispersed phase liquid flow is dispersed into small droplets when passing through the baffle holes, and two-phase turbulent motion is generated at the same time, which improves the mass transfer efficiency. When using dynamic simulation development, the dispersed phase longitudinal diffusion and axial backmixing of the pulse extraction column can be ignored. The single plate model is used, that is, the space between the two plates is used as a calculation unit, and the distribution coefficient and mass transfer coefficient of different components on each tower plate are calculated to establish the liquid-liquid two-phase mass transfer equation.

[0050] The simulation method first obtains the device attribute information, initialization information and extraction information of the pulse extraction column to be simulated. Then, according to the number of plates of the pulse extraction column, the pulse extraction column is divided into multiple single plate areas, and a single plate model is established for each single plate area. After combining multiple single plate models, the extraction calculation model of the pulse extraction column is obtained. Then, the device parameters, initialization information and extraction information are substituted into the extraction calculation model for iterative calculation. After a period of time, the running simulation data is output.

[0051] It should be noted that the advantage of using the single plate model to establish equations for each single plate region of the pulse extraction column and then solving them jointly is that it can calculate various types of pulse extraction columns. The types of pulse extraction columns mainly include annular baffle pulse extraction columns, sieve plate pulse extraction columns and conical flow pulse extraction columns. Of course, it is understandable that the expression of the established single plate model will be different for each pulse extraction column with different structures.

[0052] The following will further explain the various steps of this simulation method:

[0053] like Figure 2 and Figure 3As shown, when simulating the pulse extraction column, the device attribute information, initialization information and extraction information of the pulse extraction column are first obtained. Among them, the device attribute information includes: the number of tower plates and device parameters. The equipment parameters include: plate segment diameter, plate segment height, thickness, height from the ground, volume, upper and lower limits of liquid level measurement, diameter of the upper expansion section, height of the upper expansion section, diameter of the lower expansion section, height of the lower expansion section, inner diameter of the baffle disc, inner diameter of the ring. Initialization information includes: initial pressure, initial liquid level, initial component information. Extraction information includes: equilibrium constant, mass transfer coefficient, dispersed phase selection information, dispersed phase slip velocity. Specifically, the device attribute information, initialization information and extraction information of the pulse extraction column can be manually input into the simulation device by the staff after manual measurement.

[0054] In this embodiment, step S2: establish multiple single plate models according to the number of plates of the pulse extraction column. Specifically, it includes: according to the number of plates of the pulse extraction column, split the pulse extraction column into multiple single plate areas, and the single plate area is the area between two adjacent plates. Establish the mass conservation equation, energy conservation equation, mass transfer equation, phase equilibrium equation, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation for each single plate area respectively. The above equations are combined to obtain the single plate model of each single plate area.

[0055] Specifically, the mass transfer equations in the single plate model include: mass transfer rate equation, total volume mass transfer coefficient correlation. The phase equilibrium equation adopts Richardson's semi-empirical distribution ratio model. The extraction influencing factor correlations include: dispersed phase droplet diameter correlation, dispersed phase retention fraction correlation, and axial mixing coefficient correlation.

[0056] It should be noted that each single plate model is a set of equations established based on mass balance, energy balance, pressure flow, component balance, and phase balance, which are used to solve and calculate the single plate model. The principle is that when performing dynamic simulation, the longitudinal diffusion of the dispersed phase and the shaft system back mixing of the pulse extraction column can be ignored, and then the single plate model is adopted, that is, the space between the two plates is used as a calculation unit, and the distribution coefficient and mass transfer coefficient of different components on each plate are calculated to establish the liquid-liquid two-phase mass transfer equation.

[0057] In summary, the single-plate model includes the following equations: mass conservation equation, energy conservation equation, mass transfer rate equation, total volume mass transfer coefficient correlation equation (including baffle pulse extraction column, sieve plate pulse extraction column, tapered flow pulse extraction column), Richardson distribution ratio model (material phase equilibrium in the n-plate region), dispersed phase droplet diameter correlation equation (including baffle pulse extraction column, sieve plate pulse extraction column, tapered flow pulse extraction column), dispersed phase retention fraction correlation equation, axial mixing coefficient correlation equation, normalization equation and reaction kinetics equation.

[0058] Among them, the general mass conservation equation is used to calculate the material balance of the single plate, the general energy conservation equation is used to calculate the energy balance of the single plate, the general mass transfer equation is used to calculate the mass transfer within the single plate, the Richardson distribution ratio model is used to calculate the liquid-liquid phase balance, the extraction influencing factor correlation calculation, the reaction kinetics calculation (strip extraction) and the normalization calculation. The mass balance calculation, energy balance and phase balance calculation are the balance calculations of energy and components on the single plate of the pulse extraction column, that is, the temperature, concentration, phase fraction, density, enthalpy value and other information of the inlet and outlet logistics are calculated according to the principles of energy conservation, component conservation and phase balance.

[0059] It should also be noted that the three main types of pulse extraction columns, namely, annular baffle pulse extraction column, sieve plate pulse extraction column, and conical flow pulse extraction column, have different single-plate models, which are mainly reflected in the different expressions of the total volume mass transfer coefficient correlation equation and the extraction influencing factor correlation equation (i.e., dispersed phase droplet diameter correlation equation, dispersed phase retention fraction correlation equation, and axial mixing coefficient correlation equation). Therefore, when establishing a single-plate model, it is necessary to consider the type of pulse extraction column to be simulated and calculated to select different total volume mass transfer coefficient correlation equations and extraction influencing factor correlation equations.

[0060] The specific expressions of each of the above equations are presented below and explained:

[0061] (1) Mass conservation calculation: The mass conservation equation is used to calculate the mass conservation in the n-plate area (where subscript n represents the current single plate area, n-1 represents the previous single plate area, and n+1 represents the next single plate area). The material on each plate in the extraction tower is conserved, that is, the mass of the material entering and leaving the single plate area is conserved. The mass conservation equation expresses the conservation of solute transport in the two phases.

[0062]

[0063]

[0064] V n ·z n =V c,n ·x n +V d,n ·y n

[0065] V n =V c,n +V d,n

[0066] Where: u c ——apparent velocity of continuous phase, m / s;

[0067] u d——apparent velocity of dispersed phase, m / s;

[0068] x——continuous phase molar concentration mol / m 3 ;

[0069] y——dispersed phase molar concentration mol / m 3 ;

[0070] z——total molar concentration of two phases, mol / m 3 ;

[0071] S——Cylinder cross-sectional area m 2 ;

[0072] V——total volume of two phases, m 3 ;

[0073] V c ——volume of continuous phase, m 3 ;

[0074] V d ——dispersed phase volume, m 3 ;

[0075] N——total mass transfer rate, mol / s;

[0076] r——reaction rate, mol / (m 3 s);

[0077] (2) Energy balance calculation: Use the energy conservation equation to calculate the energy balance of the current single plate area. The energy on each plate in the extraction tower is conserved, that is, the energy of the materials entering and leaving the single plate area is conserved.

[0078]

[0079] Where: T——Liquid temperature in the equipment, K;

[0080] C P ——Specific heat capacity at constant pressure, J·Kg -1 ·K -1 ;

[0081] F——mass flow rate, Kg / s;

[0082] m——mass, Kg;

[0083] Q r ——reaction heat, J / s;

[0084] Q 散 ——heat dissipation of equipment, J / s;

[0085] The subscript in represents the logistics entering the current board area, for example, F inis the mass flow rate of the logistics entering the current single board area; out represents the logistics flowing out of the current single board area, for example, F out It is the mass flow rate of the material flowing out of the current board area. The meanings of other parameters with in and out subscripts are similar and will not be repeated here.

[0086] (3) Mass transfer calculation: The mass transfer rate equation and the total volume mass transfer coefficient correlation formula are used to perform mass transfer calculations.

[0087] 1) Mass transfer rate equation:

[0088] N=K d α(yy * )Sh

[0089] N=K c α(x * -x)Sh

[0090] Where: N——total mass transfer rate, mol / s;

[0091] K d ——Total mass transfer coefficient expressed as dispersed phase solute concentration, m / s;

[0092] K c ——Total mass transfer coefficient expressed as the concentration of the continuous phase solute, m / s;

[0093] (yy * )——Mass transfer driving force expressed as dispersed phase solute concentration, mol / m 3 ;

[0094] (x * -x)——Mass transfer driving force expressed as the concentration of the continuous phase solute, mol / m 3 ;

[0095] α——Specific surface area of ​​two-phase contact per unit volume, m 2 / m 3 ;

[0096] S——pulse extraction column cross-sectional area, m 2 ;

[0097] h——Total height of liquid level in the single plate area of ​​the pulse extraction column, m;

[0098] Among them, y * 、x * It is the solute concentration of the dispersed phase and the continuous phase after the two-phase equilibrium calculated by the phase equilibrium equation. The total mass transfer rate needs to be calculated by inputting the total mass transfer coefficient K of one phase and the contact specific surface area α of the two phases per unit volume.

[0099] 2) Total volume mass transfer coefficient correlation formula:

[0100] (Applicable to annular baffle pulse extraction column)

[0101] (Applicable to sieve plate pulse extraction column)

[0102] (Applicable to conical flow-through pulse extraction columns)

[0103]

[0104]

[0105] Where: K oc ——The total volumetric mass transfer coefficient based on the continuous phase and the dispersed terms;

[0106] ψ v ——Power dissipated per unit volume of fluid, W / m 3 ;

[0107] ——Average density of fluid in the tower, Kg / m 3 ;

[0108] k0, k1, k2, k3, k4——association parameters;

[0109] a——Specific surface area of ​​dispersed phase droplets in contact with the continuous phase;

[0110] Af——pulse extraction column pulse intensity;

[0111] (Af) t ——Critical value of pulse intensity between mixing and clarifying zone and dispersion zone, i.e. pulse intensity corresponding to minimum retention fraction, m / s;

[0112] μ c ——Continuous phase viscosity, Pa·s;

[0113] u c ——apparent velocity of continuous phase, m / s;

[0114] u d ——apparent velocity of dispersed phase, m / s;

[0115] σ——interfacial tension, N / m;

[0116] φ——dispersed phase retention fraction;

[0117] ρ c ——Density of continuous phase, Kg / m 3 ;

[0118] ρ d ——dispersed phase density;

[0119] α——Percentage of the opening area of ​​the tower plate;

[0120] ε——Cross-sectional integral ratio (D a / D) 2 ;

[0121] C0——discharge coefficient, m;

[0122] h——Total height of liquid level in the single plate area of ​​pulse extraction column, m.

[0123] (4) Phase equilibrium calculation:

[0124] The liquid-liquid two-phase equilibrium is calculated using Richardson's semi-empirical distribution ratio model, which is mainly used to calculate the material phase equilibrium in the nth single board area.

[0125] The extraction reaction of the ions to be extracted in the TBP-HNO3 system (TBP is tributyl phosphate) is generally m+ Will form M(NO3) m nTBP o The complex is then extracted into TBP, and the reaction expression is:

[0126]

[0127] The reaction thermodynamic equilibrium constant K of this reaction is:

[0128]

[0129]

[0130]

[0131]

[0132]

[0133]

[0134]

[0135]

[0136] x n =x H +2x U +4x Pu

[0137]

[0138]

[0139]

[0140]

[0141]

[0142]

[0143]

[0144] C o (TBP) is the concentration of free TBP in the organic phase under equilibrium conditions, and its expression is:

[0145]

[0146] Where: C——concentration, mol / L;

[0147] γ——activity coefficient;

[0148] Subscript o——organic phase

[0149] Subscript a——Inorganic phase

[0150] M——Metal M

[0151] C o (M(NO3) m nTBP) – the concentration of the complex in the organic phase;

[0152] C a (M m+ )——the concentration of M ions in the inorganic phase under equilibrium conditions;

[0153] ——the concentration of nitrate in the inorganic phase under equilibrium conditions;

[0154] C o n (TBP) – the concentration of free TBP in the organic phase under equilibrium conditions;

[0155] γ o ——activity coefficient in the organic phase;

[0156] γ Ma ——activity coefficient of M ions in the inorganic phase;

[0157] ——activity coefficient of nitrate ion;

[0158] ——activity coefficient of TBP;

[0159] U(VI)——U 6+ ;

[0160] PU(IV)——Pu 4+ ;

[0161] PU(III)——Pu 3+ ;

[0162] K——thermodynamic equilibrium constant;

[0163] ——apparent equilibrium constant;

[0164] ——Apparent equilibrium constant at 25℃ and 100% TBP;

[0165] τ——temperature influence factor;

[0166] F——volume fraction of TBP;

[0167] C T,TBP ——is the total molar concentration of TBP in the organic phase, mol / L, which can be obtained from C T,TBP =3.651F is obtained.

[0168] (5) Calculation of correlation formula of extraction influencing factors:

[0169] In consideration of other influencing factors of extraction, the dispersed phase droplet diameter correlation, dispersed phase retention fraction correlation and axial mixing coefficient correlation were used for calculation.

[0170] 1) Correlation formula for dispersed phase droplet diameter

[0171] (Applicable to baffled pulse extraction columns)

[0172] (Applicable to sieve plate pulse extraction column)

[0173] (Applicable to conical flow-through pulse extraction columns)

[0174] Where: d 32 ——Average droplet diameter, m;

[0175] σ——interfacial tension, N / m;

[0176] h——board spacing, m;

[0177] D——tower inner diameter, m;

[0178] D c ——Inner diameter of the disk, m;

[0179] D a ——Inner diameter of the ring, m;

[0180] ε——Cross-sectional integral ratio (D a / D) 2 ;

[0181] μ——viscosity, Pa·s;

[0182] ρ——density, Kg / m 3 ;

[0183] Δρ——Difference in density between two phases, Kg / m 3 ;

[0184] 2) Dispersed phase retention fraction correlation

[0185] (Applicable to annular baffle pulse extraction column)

[0186] (Applicable to sieve plate pulse extraction column)

[0187] (Applicable to conical flow-through pulse extraction columns)

[0188]

[0189]

[0190] Where: V m ——molar volume, m 3 / mol;

[0191] φ——dispersed phase retention fraction;

[0192] (Af) t ——Critical value of pulse intensity between mixing and clarifying zone and dispersion zone, i.e. pulse intensity corresponding to minimum retention fraction, m / s;

[0193] ρ——density, Kg / m 3 ;

[0194] Δρ——Difference in density between two phases, Kg / m 3 ;

[0195] ρ c ——Density of continuous phase, Kg / m 3 ;

[0196] ρ d ——dispersed phase density;

[0197] μ d ——Viscosity of dispersed phase, Pa·s;

[0198] μc ——Continuous phase viscosity, Pa·s;

[0199] u c ——apparent velocity of continuous phase, m / s;

[0200] u d ——apparent velocity of dispersed phase, m / s;

[0201] 3) Axial mixing coefficient correlation

[0202] (Applicable to annular baffle pulse extraction column)

[0203] (Applicable to sieve plate pulse extraction column)

[0204] (Applicable to conical flow-through pulse extraction columns)

[0205] In the formula, E c ——axial mixing coefficient;

[0206] Δρ——Difference in density between two phases, Kg / m 3 ;

[0207] μ c ——Continuous phase viscosity, Pa·s;

[0208] u c ——apparent velocity of continuous phase, m / s;

[0209] u d ——apparent velocity of dispersed phase, m / s;

[0210] φ——dispersed phase retention fraction;

[0211] σ——interfacial tension, N / m;

[0212] Af——pulse extraction column pulse intensity;

[0213] (6) Normalized equation calculation:

[0214] The normalized equation is as follows:

[0215]

[0216] Where: x i ——molar concentration of component i in the liquid phase;

[0217] x 总 ——Total molar concentration of liquid components;

[0218] y i ——molar concentration of component i in the gas phase;

[0219] y 总 ——Total molar concentration of gas phase components;

[0220] (7) Reaction kinetic equation:

[0221] For the redox reaction of the material contained in the extraction column during stripping, add the reaction kinetic equation:

[0222] Reaction kinetics equation:

[0223]

[0224]

[0225] Where: C——concentration, mol / L;

[0226] Subscript a,b,c,d: substance a,b,c,d.

[0227] k – reaction rate constant;

[0228] α, β, χ, δ——reaction order;

[0229] A——pre-exponential factor, also known as Arrhenius constant, with the same unit as k;

[0230] E——activation energy, J / mol;

[0231] R——molar gas constant, J / (mol·K);

[0232] It should be noted that the superscripts, subscripts and suffixes in the above mass conservation equation, energy conservation equation, mass transfer calculation, distribution ratio model, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation are as follows:

[0233] c – continuous phase;

[0234] d——dispersed phase;

[0235] o——organic phase;

[0236] a——water phase.

[0237] So far, the above equations constitute a complete data model of a single plate area of ​​the extraction tower. The physical quantities of the single plate model are related to the states of the upper and lower plates. To calculate the distribution of various physical quantities in the pulse extraction tower, it is necessary to jointly calculate all single plate equations to obtain the final extraction calculation model.

[0238] In this embodiment, step S4: solving the extraction calculation model according to the equipment parameters, initialization information and extraction information to obtain the operation simulation data of the pulse extraction column. Specifically, it includes the following steps:

[0239] Initialize the extraction calculation model according to the initialization information to obtain a static device model;

[0240] The extraction information is substituted into the static equipment model, and the static equipment model is iteratively calculated to obtain the operation simulation data of the pulse extraction column.

[0241] like Figure 4 As shown in FIG. 1 , a specific example of the initial properties of a device is provided. The initialization information includes the number of moles of the initial component, the initial pressure, the initial liquid level, etc. The initialization information is substituted into the extraction calculation model to obtain a static device model.

[0242] like Figure 5 The above provides a specific example of setting extraction parameters. The extraction parameters include: equilibrium constant, mass transfer coefficient, dispersed phase selection, and dispersed phase slip velocity. Among them, the equilibrium constant refers to: the concentration of the substance in the upper solvent / the concentration in the lower solvent (molar ratio). For different components, users can select the corresponding dynamic distribution ratio empirical correlation formula. During the development and design, the user-defined equipment is opened to customize the equilibrium constant equation, making the constant more operational. The equilibrium constant is K in the phase equilibrium calculation. The mass transfer coefficient refers to the coefficient of the mass transfer rate equation, and the specific surface area is the ratio of the contact area between the dispersed phase and the continuous phase to the total volume.

[0243] Furthermore, the static equipment model of the pulse extraction column can be iteratively calculated by the solving software of the computer equipment. The operation simulation data obtained by solving include: liquid level height, boundary position data of liquid surface boundary, temperature data, pressure data, upper layer liquid volume, lower layer liquid volume, component ratio data of logistics, and flow data of logistics.

[0244] In summary, the pulse extraction column process simulation method can simulate the process of the pulse extraction column and predict the actual production status of the pulse extraction column. It does not need to build a physical pulse extraction column model, which can reduce the prediction cost.

[0245] Example 2

[0246] The present invention also discloses a design verification method for a pulse extraction column, comprising:

[0247] According to the pulse extraction column process simulation method in Example 1, the operation simulation data of the pulse extraction column is calculated.

[0248] Based on the running simulation data, determine whether the pulse extraction column meets the design requirements.

[0249] Furthermore, the running simulation data includes component ratio data of the flow in the pulse extraction column. When the component ratio data in the pulse extraction column is within a preset component ratio range, it is determined that the pulse extraction column meets the design requirements; otherwise, it is determined that the pulse extraction column does not meet the design requirements.

[0250] The design verification method can be used to verify the extraction effect of the pulse extraction column, so as to optimize the pulse extraction column in a large post-processing plant.

[0251] Example 3

[0252] See also Figure 1 The present invention also discloses a pulse extraction column process simulation device, including: an acquisition module, a processing module, a control module and a calculation module.

[0253] The acquisition module is used to obtain the device attribute information, initialization information and extraction information of the pulse extraction column, and the device attribute information includes the number of plates and device parameters. The processing module is electrically connected to the acquisition module, and is used to establish multiple single-plate models according to the number of plates of the pulse extraction column. The control module is electrically connected to the acquisition module and the processing module, and is used to combine multiple single-plate models to obtain the extraction calculation model of the pulse extraction column. The calculation module is electrically connected to the acquisition module and the control module, and is used to solve the extraction calculation model according to the device parameters, initialization information and extraction information, and obtain the operation simulation data of the pulse extraction column.

[0254] It should be noted that the pulse extraction column process simulation device is a pulse extraction column simulation device of the Purex process developed and established using dynamic mechanism simulation software based on the simulation method in Example 1. It is necessary to customize the development of the core algorithm according to the device principle of the pulse extraction column, that is, to develop each device module using the mass balance relationship, energy balance relationship, mass transfer relationship, component balance relationship and phase balance relationship in the pulse extraction column to solve and calculate the device.

[0255] During the simulation process, the device attribute information, initialization information and extraction information of the pulse extraction column to be simulated are first obtained through the acquisition module. Then, the pulse extraction column is split into multiple single-plate areas according to the number of plates of the pulse extraction column through the processing module, and a single-plate model is established for each single-plate area. After the control module combines multiple single-plate models, the extraction calculation model of the pulse extraction column is obtained. Next, the calculation module is run. The calculation module calls the extraction calculation model, iteratively calculates the extraction calculation model according to the equipment parameters, initialization information and extraction information, and monitors the real-time dynamic operation data of the equipment to check the dynamic impact of process production conditions on extraction. After a period of time, the operation simulation data of the pulse extraction column is output.

[0256] The simulation device can calculate the operation simulation data of the pulse extraction column, so that the staff can predict the actual production status of the pulse extraction column according to the operation simulation data.

[0257] The following will further explain the various parts of this simulation device:

[0258] In this embodiment, the processing module includes a first processing unit, a second processing unit and a third processing unit. Among them, the first processing unit is electrically connected to the acquisition module, and is used to split the pulse extraction column into multiple single plate areas according to the number of plates of the pulse extraction column, and the single plate area is the area between two adjacent upper and lower plates. The second processing unit is electrically connected to the first processing unit, and is used to establish the mass conservation equation, energy conservation equation, mass transfer equation, phase equilibrium equation, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation for each single plate area respectively. The third processing unit is electrically connected to the second processing unit, and the above equations are combined to obtain a single plate model for each single plate area.

[0259] Among them, the first processing unit, the second processing unit and the third processing unit can all adopt the modules of the computer's dynamic simulation software, and the staff only needs to add the various solution equations of the single-plate model to the processing module. Specifically, when performing development calculations, the following equations are added to the processing model: mass conservation in the n-plate area, dispersed phase droplet diameter correlation (including baffle pulse extraction column, sieve plate pulse extraction column, cone flow pulse extraction column), dispersed phase retention fraction correlation, axial mixing coefficient correlation, total volume mass transfer coefficient correlation, energy conservation equation, distribution ratio model (material phase balance in the n-plate area), mass transfer rate equation, reaction kinetics equation. The processing module establishes the single-plate model of each single-plate area of ​​the pulse extraction column according to the above equations.

[0260] In this embodiment, the calculation module includes a first calculation unit and a second calculation unit. The first calculation unit is electrically connected to the acquisition module and the control module, and is used to initialize the extraction calculation model according to the initialization information to obtain a static device model. The second calculation unit is electrically connected to the first calculation unit and the acquisition module, and is used to substitute the device parameters and extraction information into the static device model, perform iterative calculations on the static device model, and obtain operation simulation data of the pulse extraction column.

[0261] The first calculation unit and the second calculation unit can both use the solver in the existing computer simulation software to complete the solution calculation. During the operation of the calculation module, the actual production status can also be predicted by changing the process production conditions and checking the parameter changes after the equipment has been running for a period of time.

[0262] In the present embodiment, the simulation device is specifically designed and developed in the dynamic mechanism simulation software of Zhejiang Zhongkong Intelligent Simulation Platform (VxOTS), and its development process mainly includes module algorithm development, interface development, and parameter configuration development. The equipment module is established according to the equipment principle of the pulse extraction column, and the pressure flow relationship, mass balance, energy balance, component balance, and mass transfer are used to solve and calculate the equipment. A static equipment model is established in the modeling software. According to the process production requirements, the production parameters are entered in the equipment module. According to the process requirements, the import and export logistics of the equipment are connected. Run the equipment module, monitor the dynamic operation data of the equipment, and predict the actual production status by checking the parameter changes after a certain period of operation.

[0263] Specifically, first, the device module is algorithmically developed based on the simulation method in Example 1 to obtain a processing module, a control module, and a calculation module. After the device module algorithm is developed, an interface needs to be developed to obtain a module interface, which is displayed graphically in the model software. The module interface is electrically connected to the acquisition module, and the module interface is opened to externally input various parameters.

[0264] After the equipment module is developed, it is displayed graphically in the model software. Figure 1 As shown in the figure, the module construction process of the pulse extraction column equipment is demonstrated. In the dynamic simulation software, according to the production process of the pulse extraction column, the inlet and outlet materials of the heavy phase are connected to the equipment module, and the main function of the valve is to adjust the flow.

[0265] The module interface includes the module input interface and the module output interface. When the module input interface is opened, various parameters can be externally input / selected to solve the equation. On the module output interface, the liquid level, interface, temperature, pressure, upper and lower liquid volumes, current inlet and outlet material flow components, flow rate and other parameters of the equipment output by the calculation unit after a period of calculation can be viewed in real time.

[0266] The following will further explain the module interface of this simulation device with a specific example:

[0267] like Figure 2 As shown, the basic parameter interface is displayed, where you can enter basic parameters:

[0268] According to the process flow and equipment properties of the pulse extraction column equipment, input the basic parameters of the equipment, including: equipment name, inlet logistics selection and inlet height setting (including light phase inlet and heavy phase inlet), overflow logistics selection and height setting (including gas phase outlet, heavy phase outlet, light phase outlet). At the same time, according to the production process requirements, add multiple streams of materials as light phase or heavy phase materials, that is, the selection of fluid packages (that is, the selection of thermodynamic equations).

[0269] It should be noted that the staff can directly select the fluid package that comes with the simulation software, or upload the required fluid package to the simulation software before the simulation.

[0270] like Figure 3 As shown, the device property interface is displayed, where you can enter the device parameters, including: plate segment diameter, plate segment height, thickness, height from the ground, volume, upper and lower limits of liquid level measurement, diameter of the upper expansion section, height of the upper expansion section, diameter of the lower expansion section, height of the lower expansion section, inner diameter of the baffle disc, and inner diameter of the ring.

[0271] In addition, the heat transfer parameters of the pulse extraction column can be set in the lower part of the device properties interface. It should be noted that heat transfer is mainly to set the heat dissipation of the device to the environment. The larger the environmental heat dissipation coefficient, the greater the heat dissipation to the environment. If you check whether to start heat dissipation, the model will start the heat dissipation function.

[0272] like Figure 4 As shown in the figure, the initial property interface is displayed, where the initial properties of the pulse extraction column can be set. The initial properties refer to the initial state of the device during operation, including the status information of the device. The initial property interface can set the initial pressure, initial liquid level, initial components and other information of the device.

[0273] For example, the initial properties are set to 100% N2 content, 25°C temperature, and 100kPa pressure. If [Run from initial state] is checked, it means that the pulse extraction column is filled with N2 at 25°C and 100kPa pressure in the initial state, and the device will operate in this state. This function can restore different production states of the device, and multiple production states can be flexibly switched.

[0274] like Figure 5 As shown, the parameter configuration interface is displayed, where the extraction parameters of the pulse extraction column can be set.

[0275] The parameter configuration interface is mainly used to set the extraction effect, including: equilibrium constant, mass transfer coefficient, dispersed phase selection, and dispersed phase slip velocity. The equilibrium constant refers to: the concentration of the substance in the upper solvent / the concentration in the lower solvent (molar ratio). For different components, users can select the corresponding dynamic distribution ratio empirical correlation formula. During the development and design, the user-defined equipment is open to define the equilibrium constant equation, making the constant more operable.

[0276] In addition, if Figure 6As shown, considering the stripping of the pulse extraction column, a reaction configuration interface is developed. The configuration of the reaction completes the stoichiometric ratio by selecting reactants, products and corresponding coefficients, forward orders, and reverse orders. At the same time, the reaction configuration is completed by configuring the pre-exponential factor A of the forward reaction and reverse reaction, the activation energy E, etc.

[0277] This simulation device takes into account the complexity of the flow in the extraction process. For the dispersed phase droplet size, retention fraction, axial diffusion coefficient, and total volume coefficient of mass transfer, the corresponding empirical correlation formula can be selected for configuration. The device parameters not reflected in the device property page can be realized by adding a custom device connected to the device. In the custom device, the custom variables added by yourself are associated with the selected model variables by writing a calculation script. Specifically, an external development interface can be set up, which is electrically connected to the processing module. The staff can directly input other calculation formulas into the processing module to further improve the calculation accuracy of the single-board model and the extraction calculation model.

[0278] The above is the parameter setting process of this simulation device, among which the feed / discharge flow rate, composition, and component concentration of heavy and light phase materials are mainly used for mass balance calculations. The feed / discharge temperature / enthalpy value of heavy and light phase materials, gas phase discharge temperature / enthalpy value, and equipment heat transfer (environmental heat dissipation) are mainly used for energy balance calculations. Dispersed phase selection, pulse amplitude, and pulse frequency are mainly used for mass transfer calculations. Component correlation selection, equilibrium constant, and dispersed phase selection are mainly used for phase equilibrium calculations. Extraction column type, various correlation coefficients (k), viscosity, and surface tension are mainly used for extraction influencing factor correlations. Reactants and products, coefficients, forward / reverse reaction orders, reaction phases, maximum / lowest reaction temperatures, and forward / reverse reaction coefficients are mainly used for reaction kinetics calculations.

[0279] In addition, equipment parameters including plate segment diameter, plate segment height, thickness, height from the ground, volume, upper and lower limits of liquid level measurement, diameter of the upper expansion section, height of the upper expansion section, diameter of the lower expansion section, height of the lower expansion section, inner diameter of the baffle disc, and inner diameter of the ring mainly affect the liquid level and interface of the pulse extraction column.

[0280] After completing the setting of the above parameters, start the calculation module. After the calculation module is running, a series of simulation operations will be performed to obtain the calculation results (i.e., the running simulation data), and the changes of various process parameters will be fed back in real time. Of course, during the calculation process, the module supports flexible selection of a certain process parameter to track and view its trend changes over a period of time. The software also has state operation control, such as: termination, freezing, thawing, acceleration, deceleration, etc., to achieve different equipment debugging application scenarios. In the module interface, you can view the equipment's liquid level height, interface position, temperature, pressure, upper and lower layer liquid volume, current inlet and outlet material logistics components, flow rate and other parameters in real time.

[0281] like Figure 7 As shown, the basic information of the calculation results is displayed, including the upper liquid level height, the lower liquid level height, pressure, temperature, etc.

[0282] In summary, the pulse extraction column process simulation device establishes an extraction calculation model of the pulse extraction column in the dynamic simulation software, and accurately simulates the various production status data of the equipment at different times in the post-processing production process. After changing various production input parameters in the simulation software, the equipment module and its upstream and downstream production devices will respond dynamically to reflect the production data change process of the actual production site; check the production status after the module has run for a certain period of time, and predict the changes in the actual production status. Since the extraction calculation model is a general model established based on the extraction principle and equipment principle of the pulse extraction column, the process simulation prediction of various pulse extraction columns can be achieved by changing the configurable items (equipment parameters, extraction and separation effects, inlet and outlet logistics).

[0283] Example 4

[0284] The present invention also discloses a design verification system for a pulse extraction column, comprising: a judgment module and the pulse extraction column process simulation device in embodiment 3.

[0285] The pulse extraction column process simulation device is used to calculate the operation simulation data of the pulse extraction column. The judgment module is electrically connected to the pulse extraction column process simulation device and is used to judge whether the pulse extraction column meets the design requirements according to the operation simulation data.

[0286] The specific judgment standard is: when the component ratio in the pulse extraction column is within the preset component ratio range, the pulse extraction column is judged to meet the design requirements; otherwise, the pulse extraction column is judged to not meet the design requirements.

[0287] The design verification system of the pulse extraction column can verify the extraction effect of the pulse extraction column.

[0288] It is to be understood that the above embodiments are merely exemplary embodiments used to illustrate the principles of the present invention, but the present invention is not limited thereto. For those of ordinary skill in the art, various modifications and improvements can be made without departing from the spirit and essence of the present invention, and these modifications and improvements are also considered to be within the scope of protection of the present invention.

Claims

1. A pulse extraction column process simulation method, characterized in that: include: S1: Acquire device attribute information, initialization information and extraction information of the pulse extraction column, wherein the device attribute information includes the number of plates and device parameters; S2: establishing a plurality of single plate models according to the number of plates of the pulse extraction column; S3: combining a plurality of the single plate models to obtain an extraction calculation model of the pulse extraction column; S4: solving the extraction calculation model according to the equipment parameters, initialization information and the extraction information to obtain operation simulation data of the pulse extraction column; The step S2 specifically includes: According to the number of plates of the pulse extraction column, the pulse extraction column is divided into a plurality of single plate regions, wherein the single plate region is the region between two adjacent plates; The mass conservation equation, energy conservation equation, mass transfer equation, phase equilibrium equation, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation were established for each single plate area respectively; The above equations are combined to obtain the single board model of each single board area.

2. The method according to claim 1, characterized in that The equipment parameters include: plate segment diameter, plate segment height, thickness, height from the ground, volume, upper and lower limits of liquid level measurement, upper expansion segment diameter, upper expansion segment height, lower expansion segment diameter, lower expansion segment height, baffle disc inner diameter, and ring inner diameter; The initialization information includes: initial pressure, initial liquid level, and initial component information; The extraction information includes: equilibrium constant, mass transfer coefficient, dispersed phase selection information, and dispersed phase slip velocity.

3. The method according to claim 1, characterized in that The mass transfer equation includes: a mass transfer rate equation and a total volume mass transfer coefficient correlation equation; The expression of the mass transfer rate equation is: N=K d α(yy * )Sh N=K c α(x * -x)Sh Where: N is the total mass transfer rate, mol / s; K d is the total mass transfer coefficient expressed by the dispersed phase solute concentration, m / s; K c is the total mass transfer coefficient expressed as the concentration of the continuous phase solute, m / s; (yy * ) is the mass transfer driving force expressed as the dispersed phase solute concentration, mol / m 3 ;(x * -x) is the mass transfer driving force expressed as the concentration of the continuous phase solute, mol / m 3 ; α is the contact surface area of ​​the two phases per unit volume, m 2 / m 3 ; S is the cross-sectional area of ​​the pulse extraction column, m 2 ; h is the total height of the liquid level in the single plate area of ​​the pulse extraction column, m; The phase equilibrium equation adopts Richardson semi-empirical distribution ratio model; The extraction influencing factor correlation formula includes: a dispersed phase droplet diameter correlation formula, a dispersed phase retention fraction correlation formula, and an axial mixing coefficient correlation formula.

4. The method according to claim 1, characterized in that: The step S4 specifically includes the following steps: Initializing the extraction calculation model according to the initialization information to obtain a static device model; The equipment parameters and the extraction information are substituted into the static equipment model, and the static equipment model is iteratively calculated to obtain operation simulation data of the pulse extraction column.

5. The method according to claim 4, characterized in that The operation simulation data include: liquid level height, boundary data of liquid surface boundary, temperature data, pressure data, upper layer liquid volume, lower layer liquid volume, component ratio data of logistics, and flow data of logistics.

6. A design verification method for a pulse extraction column, characterized in that: include: According to the pulse extraction column process simulation method according to any one of claims 1 to 5, calculating the operation simulation data of the pulse extraction column; According to the operation simulation data, it is determined whether the pulse extraction column meets the design requirements.

7. The method according to claim 6, characterized in that The operation simulation data includes component ratio data of the flow in the pulse extraction column; The step of judging whether the pulse extraction column meets the design requirements according to the operation simulation data specifically includes: When the component ratio data in the pulse extraction column is within a preset component ratio range, it is determined that the pulse extraction column meets the design requirements; Otherwise, it is determined that the pulse extraction column does not meet the design requirements.

8. A pulse extraction column process simulation device, characterized in that: include: Acquisition module, processing module, control module and calculation module; The acquisition module is used to acquire device attribute information, initialization information and extraction information of the pulse extraction column, wherein the device attribute information includes the number of plates and device parameters; The processing module is electrically connected to the acquisition module and is used to establish multiple single plate models according to the number of plates of the pulse extraction column; The control module is electrically connected to the processing module and is used to combine the multiple single-board models to obtain an extraction calculation model of the pulse extraction column; The calculation module is electrically connected to the acquisition module and the control module, and is used to solve the extraction calculation model according to the device parameters, initialization information and the extraction information to obtain the operation simulation data of the pulse extraction column; The processing module includes a first processing unit, a second processing unit and a third processing unit; The first processing unit is electrically connected to the acquisition module, and is used to split the pulse extraction column into a plurality of single plate areas according to the number of plates of the pulse extraction column, wherein the single plate area is an area between two adjacent upper and lower plates; The second processing unit is electrically connected to the first processing unit and is used to establish the mass conservation equation, energy conservation equation, mass transfer equation, phase equilibrium equation, extraction influencing factor correlation equation, normalization equation and reaction kinetics equation for each of the single plate regions; The third processing unit is electrically connected to the second processing unit, and combines the above equations to obtain a single board model for each single board area.

9. The pulse extraction column process simulation device according to claim 8, characterized in that: The calculation module includes a first calculation unit and a second calculation unit; The first calculation unit is electrically connected to the acquisition module and the control module, and is used to initialize the extraction calculation model according to the initialization information to obtain a static device model; The second calculation unit is electrically connected to the first calculation unit and the acquisition module, and is used to substitute the equipment parameters and the extraction information into the static equipment model, perform iterative calculation on the static equipment model, and obtain operation simulation data of the pulse extraction column.

10. A design verification system for a pulse extraction column, characterized in that: It comprises a judgment module and a pulse extraction column process simulation device according to any one of claims 8 to 9, The pulse extraction column process simulation device is used to calculate the operation simulation data of the pulse extraction column; The judgment module is electrically connected to the pulse extraction column process simulation device, and is used to judge whether the pulse extraction column meets the design requirements based on the operation simulation data.

Images