Complex reinforcement structure concrete electrochemical dechlorination computer simulation prediction method

Abstract

A computer simulation prediction method for electrochemical dechlorination of complex steel bar structure concrete comprises the following steps: step 1) geometric modeling: determining a geometric shape: determining a simulated geometric model according to an actual electrochemical dechlorination device; defining dimensions and boundaries: defining dimension parameters of the geometric model, setting boundary conditions, for electrochemical dechlorination, defining boundaries of an electrode surface and an interface of an electrolyte and an electrode, defining the electrode surface as a current input boundary, and defining the interface of the electrolyte and the electrode as a charge transfer boundary; 2) setting a physical field; step 3) chemical reaction setting; step 4) performing grid division; step 5) solving setting: selecting a solver: selecting the solver according to the established physical model and mathematical equation; setting solving parameters; simulating a dynamic process of electrochemical dechlorination, and setting a time step length. According to the method, the research and development period is remarkably shortened, the research and development cost is reduced, and high flexibility and customizability are achieved.

Description

Technical Field

[0001] The invention belongs to the field of concrete testing and relates to a method for predicting electrochemical dechlorination of concrete with a complex reinforced structure. Background Art

[0002] In order to prevent and control steel corrosion caused by chloride salt erosion and improve the durability of concrete, researchers have proposed many repair methods, such as physical repair technology, surface coating technology and electrochemical dechlorination technology. Physical repair mainly involves removing the cracked concrete protective layer and re-casting new concrete. It is not only expensive and causes great damage to the reinforced concrete structure, but also electrochemical incompatibility between the new concrete and the old concrete. Alternatively, concrete surface coating technology is used, but most coatings have problems such as short service life, easy wear and poor economic efficiency. Physical repair technology and surface coating technology can play a certain role, but they cannot fundamentally eliminate the problem of steel corrosion. Considering cost factors and adaptability in actual engineering, these repair methods are not the best choice. In comparison, electrochemical dechlorination technology has a shorter treatment cycle and does not cause direct damage to the concrete structure, which is more advantageous. The main principle of electrochemical dechlorination technology is to apply a direct current of a certain current density between the steel bar and the external anode immersed in the electrolyte solution to remove the chloride ions that have penetrated into the concrete, increase the concrete pH near the steel bar, and regenerate a passivation film on the surface of the rusted steel bar.

[0003] Since electrochemical dechlorination is a relatively long process and it is difficult to consider all kinds of influencing factors in specimen preparation, the use of computer software modeling has become a better research method. At present, the main research on electrochemical dechlorination by computer software is the current density of hydrogen evolution reaction, oxygen reduction reaction, iron oxidation reaction and the influence of oxygen content on the above reactions. The research on electrochemical dechlorination of concrete is mainly aimed at the internal ion distribution and electric field distribution of concrete. The more complex the structure, the more complex the configuration of its steel bars, and the longer the calculation time required. If the concrete grid differentiation can be replaced by other methods without affecting the calculation accuracy, or the appropriate steel bar shape can be selected, the operation efficiency can be greatly improved. Summary of the invention

[0004] In order to overcome the shortcomings of the existing technology, the present invention provides a computer simulation prediction method for electrochemical dechlorination of complex reinforced concrete structures, which can significantly shorten the R&D cycle, reduce R&D costs, and has high flexibility and customizability. The method can simultaneously and comprehensively calculate the influence of the electric field, ion concentration distribution and various boundary conditions (such as potential, exchange current density and concentration) on the concrete during the period of power-on, so as to determine whether to continue dechlorination.

[0005] The technical solution adopted by the present invention to solve its technical problem is:

[0006] A computer simulation prediction method for electrochemical dechlorination of complex reinforced concrete structure, comprising the following steps:

[0007] Step 1), geometric modeling, the process is as follows:

[0008] Determine the geometric shape: Determine the simulation geometric model based on the actual electrochemical dechlorination device.

[0009] Define size and boundaries: clarify the size parameters of the geometric model and set boundary conditions. For electrochemical dechlorination, define the boundaries of the electrode surface and the interface between the electrolyte and the electrode. Define the electrode surface as the current input boundary and the interface between the electrolyte and the electrode as the charge transfer boundary.

[0010] Step 2) Physical field setting, the process is as follows:

[0011] Select the physical field module: Electrochemical dechlorination involves the physical processes of current conduction and material diffusion. Select the physical field module, including the "Transport of Diluted Species" and "Current" interfaces in the "Electrochemistry" module;

[0012] Set physical field parameters: according to the electrochemical system and the dechlorination process, set the electrolyte conductivity, diffusion coefficient, initial chlorine concentration parameters, as well as the material properties and electrode potential of the electrode;

[0013] Step 3) Chemical reaction setup, the process is as follows:

[0014] Define the electrochemical reaction. In the "Electrochemistry" module, define the electrochemical reaction through the "Chemical Reaction" node. The electrochemical dechlorination reaction is the oxidation reaction of chloride ions at the anode. Enter the reaction equation and reaction rate constant information.

[0015] Coupling physical fields and chemical reactions: coupling the defined electrochemical reactions with the current and material transfer physics fields to achieve the interaction between the physical process and the chemical reaction;

[0016] Step 4) Grid division, the process is as follows:

[0017] Select the mesh type: Select the mesh type according to the characteristics of the geometric model and the calculation accuracy requirements;

[0018] Generate mesh: mesh the geometric model and generate discrete mesh units.

[0019] Step 5) Solve the settings. The process is as follows:

[0020] Select solver: Select the solver based on the established physical model and mathematical equations;

[0021] Set solution parameters: Simulate the dynamic process of electrochemical dechlorination and set the time step to ensure the stability and accuracy of the calculation.

[0022] Further, in step 5), different electrode kinetic expression types are used for the cathode and anode of the electrochemical dechlorination device, and the concentration-dependent kinetic expression is selected for the cathode to couple the electric field with the concentration. The calculation formula for the local current density is:

[0023]

[0024] Where i0 is the cathode exchange current density, unit A / m 3 , is the exchange current density when iron undergoes oxidation reaction, C Rexp is the local reduction product expression, which is c / c0, c is the instantaneous ion concentration in mol / L, c0 is the initial ion concentration in mol / L, C 0exp is the oxide expression, 1, α a is the anode transfer coefficient, α C is the cathode transfer coefficient, α a +α C =1, F is Faraday's constant, which is 9.648×10 4 C / mol, R is the ideal gas constant, which is 8.134 J / (mol·K), T is the absolute temperature, which is 298 K, where η is the potential difference from anode to cathode, in V;

[0025] For the anode-selective anode Tafel equation kinetic expression, the local current density is calculated as:

[0026]

[0027] In the formula, A a is the anode Tafel slope, is the Tafel slope of the copper oxidation reaction, i0 is the anode exchange current density, unit A / m 3 , is the exchange current density when copper undergoes oxidation reaction;

[0028] The ion concentration distribution is calculated using the dilute species transfer model. Based on the mass conservation equation, the transmission of different ion types is expressed as:

[0029]

[0030] In the formula, c i The concentration of ions in concrete in mol / L, t is time in s, J i It is the type of ions in concrete, unit is mol / (m 2 ·s), N is the total number of ionic species in concrete;

[0031] The removal of chloride ions from chloride-contaminated concrete follows the principle of ion mass transfer under the influence of an externally applied electric field. The establishment of the mass conservation equation is based on the unit volume of concrete, so the moles of ions passing through the unit area of ​​concrete per unit time are defined as the ion flux, which is described by the Nernst-Planck equation:

[0032] J i =-D i ▽c i -z i u m，， Fc i ▽V

[0033] Where D i is the diffusion coefficient of different ions in concrete, in m 2 ·s,c i is the concentration of the ion at that point in mol / L, z i is the ionic charge, cations are positive, anions are negative, and F is the Faraday constant, which is 9.648×10 4 C / mol, V is the potential difference, unit V, u m，i is the mobility of concrete when electrified, which is described by the Nernst-Planck equation:

[0034]

[0035] Where D C is the ion diffusion coefficient during the power-on process, in m 2 ·s, R is the ideal gas constant, which is 8.134 J / (mol·K), and F is the Faraday constant, which is 9.648×10 4 C / mol;

[0036] The ion concentration of the concrete can be coupled with the electroneutrality condition and current protection:

[0037]

[0038] Where I is the current density applied to the pore solution in concrete, which is converted into an external voltage to obtain:

[0039]

[0040] The above formula represents the equilibrium relationship between the internal ion concentration and potential at any time during the ion transport process in concrete. The ion concentration and potential gradient are calculated based on the initial conditions, boundary conditions and the concentration of bound ions. F is the Faraday constant, which is 9.648×10 4C / mol, R is the ideal gas constant, which is 8.134 J / (mol·K), T is the absolute temperature, which is 298 K, I is the current density applied to the pore solution in the concrete, z i is the ionic charge, D i is the diffusion coefficient of different ions in concrete, c i is the concentration of ions at that point, Represents the potential difference across a point.

[0041] Furthermore, the method further comprises the following steps:

[0042] Step 6) Result analysis and post-processing, the process is as follows:

[0043] Visualization results: After the solution is completed, the calculation results can be visualized and analyzed. The chloride ion concentration distribution, current density distribution and potential distribution cloud map can be drawn to intuitively understand the changes of various physical quantities in the electrochemical dechlorination process;

[0044] Data extraction and analysis: Extract physical quantity data at key locations, such as the current density on the electrode surface and the change in chloride ion concentration in the electrolyte over time, and conduct further analysis and processing, including calculating performance indicators such as chlorine removal efficiency and energy consumption.

[0045] The technical concept of the present invention is: the electric field calculation during the concrete electrification adopts a tertiary current distribution module, because the tertiary current distribution not only covers the potential factor in the primary current distribution and the mass transfer factor in the secondary current distribution when considering the current distribution, but also focuses on the deep polarization effect. When the current density is extremely high, the concentration of the reactants near the electrode will be significantly reduced, resulting in the influence of the mass transfer process on the current distribution becoming extremely significant. This deep polarization effect makes the current distribution more complicated, and comprehensively considers the material transfer process in the electrolyte, including diffusion, electromigration and convection. The flux of ions is described by the Nernst-Planck equation, so as to accurately calculate the current density. Compared with the assumption that the reactants are uniformly distributed and can be supplied to the reaction site infinitely fast in the primary current distribution, and the secondary current distribution only focuses on a single mass transfer process such as diffusion, it is closer to the actual situation and can more accurately reflect the physical phenomena in the electrochemical system. In the past, the tertiary current distribution mainly considered the primary cell reaction, and the present invention eliminated the influence of the convection factor and is suitable for concrete structures.

[0046] The present invention can keep the final calculation result close when the steel bars are removed, the steel bar boundary is used as the cathode boundary, the influence of steel bar corrosion on chlorine removal is not considered, and the concrete working surface is used as the anode boundary, thereby fully balancing the calculation efficiency and calculation accuracy. Secondly, the present invention can simultaneously calculate the concrete electric field distribution and ion concentration distribution, taking into account the bidirectional coupling of electric field migration and mass transfer of porous media, and has a relatively sufficient theoretical basis.

[0047] The beneficial effects of the present invention are mainly manifested in:

[0048] 1. With the powerful multi-physics coupling simulation capability of COMSOL software, the present invention can deeply explore the microscopic physical mechanism in the electrochemical dechlorination process, so as to more accurately optimize the dechlorination process parameters. Compared with the traditional trial and error method, it can significantly shorten the R&D cycle and reduce R&D costs. For example, a large number of parameter combinations can be screened out in a short time through simulation, while traditional experiments may take weeks or even months to test one by one.

[0049] 2. The present invention is highly flexible and customizable. Since the software model can easily modify the geometric structure, material properties and physical field parameters, it can quickly adapt to different application scenarios and demand changes. For example, when treating wastewater with different chloride ion concentrations, it is only necessary to adjust the electrolyte concentration parameters in the model and re-simulate to obtain the corresponding optimization solution without large-scale transformation of the entire reaction device.

[0050] 3. Through the feedback correction mechanism of real-time monitoring data and COMSOL model, the stability and reliability of the electrochemical dechlorination system in the actual operation process can be effectively improved, ensuring that the dechlorination effect always remains at a high level, and reducing the problems of decreased dechlorination efficiency and increased energy consumption caused by factors such as operating condition fluctuations. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a reinforced concrete structure modeling diagram.

[0052] Figure 2 It is the vector diagram of the electrolyte potential and current density in the plane concrete. DETAILED DESCRIPTION

[0053] The present invention will be further described below in conjunction with the accompanying drawings.

[0054] Reference Figure 1 and Figure 2 , a computer simulation prediction method for electrochemical dechlorination of complex reinforced concrete structure, comprising the following steps:

[0055] Step 1), geometric modeling, the process is as follows:

[0056] Determine the geometric shape: Determine the simulation geometric model based on the actual electrochemical dechlorination device.

[0057] Define size and boundaries: clarify the size parameters of the geometric model and set boundary conditions. For electrochemical dechlorination, define the boundaries of the electrode surface and the interface between the electrolyte and the electrode. For example, define the electrode surface as the current input boundary and the interface between the electrolyte and the electrode as the charge transfer boundary.

[0058] Step 2) Physical field setting, the process is as follows:

[0059] Select the physical field module: In COMSOL software, electrochemical dechlorination involves physical processes such as current conduction and material diffusion, so it is necessary to select the corresponding physical field module, such as the "Transport of Diluted Species" and "Current" interfaces in the "Electrochemistry" module;

[0060] Setting physical field parameters: According to the electrochemical system and the dechlorination process, set the relevant physical field parameters, for example, set the conductivity, diffusion coefficient, initial chlorine concentration and other parameters of the electrolyte, as well as the material properties of the electrode, electrode potential, etc.;

[0061] Step 3) Chemical reaction setup, the process is as follows:

[0062] Define the electrochemical reaction. In the "Electrochemistry" module, define the electrochemical reaction through the "Chemical Reaction" node. For electrochemical dechlorination, the main reaction is the oxidation reaction of chloride ions at the anode. Accurately enter the reaction equation, reaction rate constant and other related information;

[0063] Coupling physical fields and chemical reactions: coupling defined electrochemical reactions with physical fields such as current and material transfer to achieve interaction between physical processes and chemical reactions;

[0064] Step 4) Grid division, the process is as follows:

[0065] Select the mesh type: Select the mesh type, such as triangular mesh, quadrilateral mesh, etc., according to the characteristics of the geometric model and the requirements for calculation accuracy. For areas with complex geometric shapes or drastic changes in physical fields, the local mesh encryption method can be used to improve the calculation accuracy.

[0066] Generate mesh: mesh the geometric model and generate discrete mesh units. When meshing, it is necessary to control the size and quality of the mesh to ensure the accuracy and stability of the calculation results.

[0067] Step 5) Solve the settings. The process is as follows:

[0068] Select solver: Select a solver based on the established physical model and mathematical equations. Computer software provides a variety of solvers, such as direct solvers, iterative solvers, etc., which can be selected according to the scale and nature of the problem;

[0069] Set solution parameters: including solution time, time step, convergence criteria and other parameters. For the dynamic process simulation of electrochemical dechlorination, the time step needs to be set reasonably to ensure the stability and accuracy of the calculation;

[0070] Different electrode kinetic expression types are used for the cathode and anode of the electrochemical dechlorination device. For the cathode, the concentration-dependent kinetic expression is selected to couple the electric field with the concentration. The calculation formula for the local current density is:

[0071]

[0072] Where i0 is the cathode exchange current density, unit A / m 3 , is the exchange current density when iron undergoes oxidation reaction, C Rexp is the local reduction product expression, which is c / c0, c is the instantaneous ion concentration in mol / L, c0 is the initial ion concentration in mol / L, C 0exp is the oxide expression, 1, α a is the anode transfer coefficient, α C is the cathode transfer coefficient, α a +α C =1, F is Faraday's constant, which is 9.648×10 4 C / mol, R is the ideal gas constant, which is 8.134 J / (mol·K), T is the absolute temperature, which is 298 K, where η is the potential difference from anode to cathode, in V;

[0073] For the anode-selective anode Tafel equation kinetic expression, the local current density is calculated as:

[0074]

[0075] In the formula, A a is the anode Tafel slope, is the Tafel slope of the copper oxidation reaction, i0 is the anode exchange current density, unit A / m 3 , is the exchange current density when copper undergoes oxidation reaction;

[0076] The ion concentration distribution is calculated using the dilute species transfer model. Based on the mass conservation equation, the transmission of different ion types is expressed as:

[0077]

[0078] In the formula, ci The concentration of ions in concrete in mol / L, t is time in s, J i It is the type of ions in concrete, unit is mol / (m 2 ·s), N is the total number of ionic species in concrete;

[0079] The removal of chloride ions from chloride-contaminated concrete follows the principle of ion mass transfer under the influence of an externally applied electric field. The establishment of the mass conservation equation is based on the unit volume of concrete, so the moles of ions passing through the unit area of ​​concrete per unit time are defined as the ion flux, which is described by the Nernst-Planck equation:

[0080] J i =-D i ▽c i -z i u m，， Fc i ▽V

[0081] Where D i is the diffusion coefficient of different ions in concrete, in m 2 ·s,c i is the concentration of the ion at that point in mol / L, z i is the ionic charge, cations are positive, anions are negative, and F is the Faraday constant, which is 9.648×10 4 C / mol, V is the potential difference, unit V, u m，i is the mobility of concrete when electrified, which is described by the Nernst-Planck equation:

[0082]

[0083] Where D C is the ion diffusion coefficient during the power-on process, in m 2 ·s, R is the ideal gas constant, which is 8.134 J / (mol·K), and F is the Faraday constant, which is 9.648×10 4 C / mol;

[0084] The ion concentration of the concrete can be coupled with the electroneutrality condition and current protection:

[0085]

[0086] Where I is the current density applied to the pore solution in concrete, which is converted into an external voltage to obtain:

[0087]

[0088] The above formula represents the equilibrium relationship between the internal ion concentration and potential at any time during the ion transport process in concrete. The ion concentration and potential gradient are calculated based on the initial conditions, boundary conditions and the concentration of bound ions. F is the Faraday constant, which is 9.648×10 4 C / mol, R is the ideal gas constant, which is 8.134 J / (mol·K), T is the absolute temperature, which is 298 K, I is the current density applied to the pore solution in the concrete, z i is the ionic charge, D i is the diffusion coefficient of different ions in concrete, c i is the concentration of ions at that point, Represents the potential difference between the two sides of a point;

[0089] During the electrochemical dechlorination process, the electrode reaction occurs on the surface of the cathode steel bar and the anode copper sheet. The oxidation and hydrogen evolution reaction on the cathode steel bar surface is expressed as:

[0090] 2H2O+O2+4e - =4OH -

[0091] 2H2O+2e - =2OH - +H2

[0092] Only hydroxyl ions are generated on the cathode steel rod surface, while the magnetic flux of other species is zero. Therefore, the magnetic flux of hydroxyl ions generated is equal to the externally applied current density. The volume of the electrolyte is large enough. Therefore, although chloride ions are transported into the electrolyte over time, the chloride ions in the electrolyte are still negligible compared to the chloride content in the corrosion layer. In addition, in actual ECE treatment, the anodic electrolyte must be replaced regularly. Therefore, during the ECE process, the ion concentration at the anode is assumed to be constant. Due to the lack of current density and ion transport, the boundary conditions without anode or cathode are assumed to be zero flow.

[0093] Step 6) Result analysis and post-processing, the process is as follows:

[0094] Visualization of results: After the solution is completed, the calculation results can be visualized and analyzed using the post-processing function of the software. Cloud maps such as chloride ion concentration distribution, current density distribution, and potential distribution can be drawn to intuitively understand the changes in various physical quantities during the electrochemical dechlorination process.

[0095] Data extraction and analysis: Extract physical quantity data at key locations, such as the current density on the electrode surface, the change in chloride ion concentration in the electrolyte over time, etc., and conduct further analysis and processing, such as calculating performance indicators such as chlorine removal efficiency and energy consumption.

[0096] Through the above steps, this embodiment can use computer software to perform a more comprehensive simulation and analysis of the electrochemical dechlorination process, and provide theoretical support for the optimization and design of the electrochemical dechlorination technology. It should be noted that in the actual operation process, it is necessary to appropriately simplify and adjust the model according to specific problems and research purposes, and verify and correct the model in combination with experimental data to improve the accuracy and reliability of the model.

[0097] The contents described in the embodiments of this specification are merely enumerations of implementation forms of the inventive concept and are for illustrative purposes only. The protection scope of the present invention should not be considered to be limited to the specific forms described in this embodiment, and the protection scope of the present invention also extends to equivalent technical means that can be thought of by ordinary technicians in this field based on the inventive concept.

Claims

1. A computer simulation prediction method for electrochemical dechlorination of complex reinforced concrete, characterized in that: The method comprises the following steps: Step 1), geometric modeling, the process is as follows: Determine the geometric shape: Determine the simulation geometric model based on the actual electrochemical dechlorination device; Define size and boundaries: clarify the size parameters of the geometric model and set boundary conditions. For electrochemical dechlorination, define the boundaries of the electrode surface and the interface between the electrolyte and the electrode. Define the electrode surface as the current input boundary and the interface between the electrolyte and the electrode as the charge transfer boundary. Step 2) Physical field setting, the process is as follows: Select the physical field module: Electrochemical dechlorination involves the physical processes of current conduction and material diffusion. Select the physical field module, including the "Transport of Diluted Species" and "Current" interfaces in the "Electrochemistry" module; Set physical field parameters: according to the electrochemical system and the dechlorination process, set the electrolyte conductivity, diffusion coefficient, initial chlorine concentration parameters, as well as the material properties and electrode potential of the electrode; Step 3) Chemical reaction setup, the process is as follows: Define the electrochemical reaction. In the "Electrochemistry" module, define the electrochemical reaction through the "Chemical Reaction" node. The electrochemical dechlorination reaction is the oxidation reaction of chloride ions at the anode. Enter the reaction equation and reaction rate constant information. Coupling physical fields and chemical reactions: coupling the defined electrochemical reactions with the current and material transfer physics fields to achieve the interaction between the physical process and the chemical reaction; Step 4) Grid division, the process is as follows: Select the mesh type: Select the mesh type according to the characteristics of the geometric model and the calculation accuracy requirements; Generate mesh: mesh the geometric model and generate discrete mesh units. Step 5) Solve the settings. The process is as follows: Select solver: Select the solver based on the established physical model and mathematical equations; Set solution parameters: Simulate the dynamic process of electrochemical dechlorination and set the time step to ensure the stability and accuracy of the calculation.

2. A computer simulation prediction method for electrochemical dechlorination of complex reinforced concrete structure as claimed in claim 1, characterized in that: In the step 5), different electrode kinetic expression types are used for the cathode and anode of the electrochemical dechlorination device. For the cathode, a concentration-dependent kinetic expression is selected to couple the electric field with the concentration. The calculation formula for the local current density is: Where i0 is the cathode exchange current density, unit A / m 3 , is the exchange current density when iron undergoes oxidation reaction, C Rexp is the local reduction product expression, which is c / c0, c is the instantaneous ion concentration in mol / L, c0 is the initial ion concentration in mol / L, C 0exp is the oxide expression, 1, α a is the anode transfer coefficient, α C is the cathode transfer coefficient, α a +α C =1, F is Faraday's constant, which is 9.648×10 4 C / mol, R is the ideal gas constant, which is 8.134 J / (mol·K), T is the absolute temperature, which is 298 K, where η is the potential difference from anode to cathode, in V; For the anode-selective anode Tafel equation kinetic expression, the local current density is calculated as: In the formula, A a is the anode Tafel slope, is the Tafel slope of the copper oxidation reaction, i0 is the anode exchange current density, unit A / m 3 , is the exchange current density when copper undergoes oxidation reaction; The ion concentration distribution is calculated using the dilute species transfer model. Based on the mass conservation equation, the transmission of different ion types is expressed as: In the formula, c i The concentration of ions in concrete in mol / L, t is time in s, J i It is the type of ions in concrete, unit is mol / (m 2 ·s), N is the total number of ionic species in concrete; The removal of chloride ions from chloride-contaminated concrete follows the principle of ion mass transfer under the influence of an externally applied electric field. The establishment of the mass conservation equation is based on the unit volume of concrete, so the moles of ions passing through the unit area of ​​concrete per unit time are defined as the ion flux, which is described by the Nernst-Planck equation: Where D i is the diffusion coefficient of different ions in concrete, in m 2 ·s,c i is the concentration of the ion at that point in mol / L, z i is the ionic charge, cations are positive, anions are negative, and F is the Faraday constant, which is 9.648×10 4 C / mol, V is the potential difference, unit V, u m，i is the mobility of concrete when electrified, which is described by the Nernst-Planck equation: Where D C is the ion diffusion coefficient during the power-on process, in m 2 ·s, R is the ideal gas constant, which is 8.134 J / (mol·K), and F is the Faraday constant, which is 9.648×10 4 C / mol; The ion concentration of the concrete can be coupled with the electroneutrality condition and current protection: Where I is the current density applied to the pore solution in concrete, which is converted into an external voltage to obtain: The above formula represents the equilibrium relationship between the internal ion concentration and potential at any time during the ion transport process in concrete. The ion concentration and potential gradient are calculated based on the initial conditions, boundary conditions and the concentration of bound ions. F is the Faraday constant, which is 9.648×10 4 C / mol, R is the ideal gas constant, which is 8.134 J / (mol·K), T is the absolute temperature, which is 298 K, I is the current density applied to the pore solution in the concrete, z i is the ionic charge, D i is the diffusion coefficient of different ions in concrete, c i is the concentration of ions at that point, Represents the potential difference across a point.

3. A computer simulation prediction method for electrochemical dechlorination of complex reinforced concrete structure as claimed in claim 1 or 2, characterized in that: The method further comprises the following steps: Step 6) Result analysis and post-processing, the process is as follows: Visualization results: After the solution is completed, the calculation results can be visualized and analyzed. The chloride ion concentration distribution, current density distribution and potential distribution cloud map can be drawn to intuitively understand the changes of various physical quantities in the electrochemical dechlorination process; Data extraction and analysis: Extract physical quantity data at key locations, such as the current density on the electrode surface and the change in chloride ion concentration in the electrolyte over time, and conduct further analysis and processing, including calculating performance indicators such as chlorine removal efficiency and energy consumption.

Images