CN113533138A - Quantification method for reducing ion transmission rate of cement-based material pore channel - Google Patents
Quantification method for reducing ion transmission rate of cement-based material pore channel Download PDFInfo
- Publication number
- CN113533138A CN113533138A CN202110648065.7A CN202110648065A CN113533138A CN 113533138 A CN113533138 A CN 113533138A CN 202110648065 A CN202110648065 A CN 202110648065A CN 113533138 A CN113533138 A CN 113533138A
- Authority
- CN
- China
- Prior art keywords
- cement
- potential energy
- molecules
- based material
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011148 porous material Substances 0.000 title claims abstract description 41
- 239000004568 cement Substances 0.000 title claims abstract description 37
- 230000005540 biological transmission Effects 0.000 title claims abstract description 30
- 239000000463 material Substances 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000011002 quantification Methods 0.000 title claims abstract description 9
- 150000002500 ions Chemical class 0.000 claims abstract description 22
- 239000000243 solution Substances 0.000 claims abstract description 19
- 239000012266 salt solution Substances 0.000 claims abstract description 17
- ISYWECDDZWTKFF-UHFFFAOYSA-N nonadecanoic acid Chemical compound CCCCCCCCCCCCCCCCCCC(O)=O ISYWECDDZWTKFF-UHFFFAOYSA-N 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000000329 molecular dynamics simulation Methods 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 8
- 238000005381 potential energy Methods 0.000 claims description 32
- 230000035515 penetration Effects 0.000 claims description 14
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 12
- 230000033001 locomotion Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- CSNNHWWHGAXBCP-UHFFFAOYSA-L Magnesium sulfate Chemical compound [Mg+2].[O-][S+2]([O-])([O-])[O-] CSNNHWWHGAXBCP-UHFFFAOYSA-L 0.000 claims description 6
- 238000005452 bending Methods 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 235000002639 sodium chloride Nutrition 0.000 claims description 5
- 239000011780 sodium chloride Substances 0.000 claims description 5
- PMZURENOXWZQFD-UHFFFAOYSA-L Sodium Sulfate Chemical compound [Na+].[Na+].[O-]S([O-])(=O)=O PMZURENOXWZQFD-UHFFFAOYSA-L 0.000 claims description 4
- 150000001450 anions Chemical class 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 150000001768 cations Chemical class 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 4
- 229910052938 sodium sulfate Inorganic materials 0.000 claims description 4
- 235000011152 sodium sulphate Nutrition 0.000 claims description 4
- UXVMQQNJUSDDNG-UHFFFAOYSA-L Calcium chloride Chemical compound [Cl-].[Cl-].[Ca+2] UXVMQQNJUSDDNG-UHFFFAOYSA-L 0.000 claims description 3
- AIYUHDOJVYHVIT-UHFFFAOYSA-M caesium chloride Chemical compound [Cl-].[Cs+] AIYUHDOJVYHVIT-UHFFFAOYSA-M 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000013078 crystal Substances 0.000 claims description 3
- 238000005421 electrostatic potential Methods 0.000 claims description 3
- 229910052943 magnesium sulfate Inorganic materials 0.000 claims description 3
- 235000019341 magnesium sulphate Nutrition 0.000 claims description 3
- 230000002040 relaxant effect Effects 0.000 claims description 2
- 230000037427 ion transport Effects 0.000 claims 2
- 239000000126 substance Substances 0.000 claims 1
- 150000001735 carboxylic acids Chemical class 0.000 abstract description 19
- 238000010276 construction Methods 0.000 abstract description 3
- 230000036571 hydration Effects 0.000 description 6
- 238000006703 hydration reaction Methods 0.000 description 6
- 230000003628 erosive effect Effects 0.000 description 5
- 230000009471 action Effects 0.000 description 4
- 239000007863 gel particle Substances 0.000 description 4
- 239000000758 substrate Substances 0.000 description 4
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 239000011398 Portland cement Substances 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical group [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000004593 Epoxy Substances 0.000 description 1
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 description 1
- 238000005411 Van der Waals force Methods 0.000 description 1
- OBNDGIHQAIXEAO-UHFFFAOYSA-N [O].[Si] Chemical compound [O].[Si] OBNDGIHQAIXEAO-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000002019 doping agent Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000002209 hydrophobic effect Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- MKTRXTLKNXLULX-UHFFFAOYSA-P pentacalcium;dioxido(oxo)silane;hydron;tetrahydrate Chemical compound [H+].[H+].O.O.O.O.[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-][Si]([O-])=O.[O-][Si]([O-])=O.[O-][Si]([O-])=O.[O-][Si]([O-])=O.[O-][Si]([O-])=O.[O-][Si]([O-])=O MKTRXTLKNXLULX-UHFFFAOYSA-P 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- -1 silylene, octadecanecarboxylic acid Chemical compound 0.000 description 1
- 229910001415 sodium ion Inorganic materials 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010189 synthetic method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/04—Investigating osmotic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N2013/003—Diffusion; diffusivity between liquids
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a quantification method for reducing the transmission rate of cement-based material pore ions, which comprises the steps of molecular dynamics model construction, calculation simulation and transmission rate quantification. The molecular dynamics model consists of three parts: C-S-H gel pore canal, octadecane carboxylic acid and salt solution. The method can simulate the transmission behavior of ions in the C-S-H pore channel under different solution environments, quantify the influence of octadecanoic carboxylic acid on the transmission rate of the C-S-H pore channel, and calculate the permeation rate of the ions in the C-S-H pore channel under the established environment, thereby providing a basis for predicting the durability life of the cement-based material.
Description
Technical Field
The invention relates to calculation of ion transmission rate, in particular to a quantification method for reducing ion transmission rate of a cement-based material pore channel.
Background
Cement concrete is the main material for constructing highways, bridges, houses, tunnels, etc. Cement concrete structures exposed to severe environments such as sea salt areas, frozen areas, saline soil and the like are often degraded prematurely due to invasion of erosion media, and finally the structures are damaged or even fail. Thus, structural durability is closely linked to the microscopic resistance of the cement-based structure to ionic attack. The C-S-H gel is a main product of Portland cement hydration, is a main body of a cement hydration gel phase, exists in a form of gel particles of 1-100 nm, pores among the gel particles provide channels for transmission of an erosion medium, and the hydrophilic property of the surface of the C-S-H gel provides power for capillary transmission, so that the transmission resistance of the traditional cement concrete is further reduced.
For this reason, more and more foreign and domestic scholars start to reduce the transmission of cement concrete by means of nano-doping and surface coating. Mainly comprises the following steps: 1) before the cement-based material is formed, graphene oxide, epoxy, silylene, octadecanecarboxylic acid and the like are doped to participate in the cement hydration process and interact with cement hydration products; 2) after the cement-based material is formed, the nano-dopant is made into a solution to be coated on the surface of the cement-based material. The method can improve the anti-transmission capability of the cement-based material to a certain extent. Particularly, when the octadecanoic carboxylic acid is coated on the cement-based material, the octadecanoic carboxylic acid molecules can be attached to the surfaces of the cement-based material pore canals, the hydrophilicity of the surfaces of the pore canals is reduced to a certain degree, the capillary transmission rate is reduced, and therefore the transmission performance of the cement-based material is greatly reduced, and the research significance is high.
However, the principle of the octadecanoic carboxylic acid for reducing the transmission inhibition of the erosion medium in the cement-based material is still unclear at present, the existing experimental means cannot fully illustrate the action mechanism, and the existing experimental means cannot quantitatively characterize the erosion property of the octadecanoic carboxylic acid for inhibiting the cement-based material, so that the related engineering application of the octadecanoic carboxylic acid for reducing the transmission of the cement-based material is few.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a quantification method for reducing the ion transmission rate of a cement-based material pore passage, and solves the problem that the corrosion of a cement-based material inhibited by octadecanoic carboxylic acid cannot be quantified and characterized.
The technical scheme is as follows: the invention aims to provide a quantification method for reducing the ion transmission rate of a cement-based material pore channel, which comprises the following steps:
(1) establishing a molecular dynamics model of C-S-H gel molecules, octadecane carboxylic acid molecules and salt solution, constructing the C-S-H molecules, the octadecane carboxylic acid molecules and the salt solution molecules, and relaxing under an NVT (noise, vibration and harshness) ensemble to enable all parts of the model to be close to each other to form a stable molecular structure;
(2) applying a force to each atom in the molecular model created above, the force field being of the form shown below:
Etotal=Ebond+Eangle+Etorsion+Evdw+Ecoul
in the formula, EtotalRepresents the total potential energy, EbondIndicating key expansion potential energy, EangleIndicating the bending potential energy of the key angle, EtorsionRepresenting dihedral torsional potential energy, EvdwDenotes van der Waals' force, EcoulRepresenting the coulomb electrostatic potential;
(3) setting preset environment temperature, a calculation time step and the outer edge of an integral molecular model as periodic boundary conditions, ensuring the system stability in the simulation process by using a Nose-Hoover constant temperature calculation method, and obtaining an atom motion trajectory by using a Verlet algorithm;
(4) and (4) according to the atom motion coordinate track obtained in the step (3), obtaining the positions and the penetration depths of the ions in the pore channels at different times, drawing a correlation curve of the penetration depths and the time, and fitting the obtained correlation curve of the penetration depths and the time to obtain the diffusion coefficient.
When the molecular dynamics model is established in the step (1), octadecanecarboxylic acid molecules are placed on the surface of the inner wall of the C-S-H molecule, and the salt solution molecules are fully distributed at the spatial position below the C-S-H molecule.
The octadecane carboxylic acid model in the step (1) is the same as the actual molecular structure.
In the step (1), the salt solution is at least one of a sodium chloride solution, a calcium chloride solution, a cesium chloride solution, a sodium sulfate solution and a magnesium sulfate solution, and the construction of the salt solution model is to mix a certain amount of anions and cations with water molecules so as to obtain the salt solution with a certain concentration.
The vertical formula of the step (2) is developed as follows:
in the formula, the first, second, third and fourth terms on the right respectively represent key expansion potential energy, key angle bending potential energy, dihedral angle twisting potential energy and non-key energy, specifically, k1i、k2i、k3iRepresenting the elastic constant of the key expansion term; Δ l represents the bond length; Δ θ represents a bond angle; ofijAnd σijRepresents a potential energy parameter (epsilon)ijPotential energy difference between reaction atoms, sigmaijReflecting equilibrium distance between atoms); r isijRepresents an interatomic distance; q. q.si(j)Represents the charge of the ith (j) ion in the molecule; cos (nw-r) represents the dihedral angle parameter; epsilon0Representing a potential energy parameter.
Has the advantages that: the invention can reveal the interaction mechanism of octadecanoic carboxylic acid and C-S-H gel pore canal on the molecular scale, quantify the influence of octadecanoic carboxylic acid on the transmission rate of the C-S-H gel pore canal, make up the defect of the molecular scale representation by the experimental means, provide the basis for the durability life of the cement-based material by the obtained data, and have guiding significance for the experimental research and the engineering application of the octadecanoic carboxylic acid for improving the anti-transmission performance of the cement-based material, particularly for the research on the influence of the durability of the cement-based material.
Drawings
FIG. 1 is a general schematic diagram of a molecular dynamics model;
FIG. 2 is a schematic diagram of a C-S-H molecule;
FIG. 3 shows (a) a model of the solution transport through the C-S-H channels, (b) a depth of penetration as a function of time, the percentages representing different octadecanecarboxylic acid loadings, and the C-S-H blank.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention discloses a quantification method for reducing the ion transmission rate of a cement-based material pore passage, which comprises the following steps:
(1) establishing a molecular dynamics model of C-S-H gel molecules, octadecane carboxylic acid molecules and salt solution, wherein the model comprises three parts as shown in figure 1: C-S-H gel molecules, octadecanecarboxylic acid molecules and salt solution (water molecules and ions), wherein the octadecanecarboxylic acid molecules are placed on the inner wall surface of the C-S-H molecules, the salt solution molecules are fully distributed at the spatial position below the C-S-H molecules, and all parts of the model are relaxed by 100-200 ps under an NVT (noise, vibration and harshness) ensemble to be close to each other, so that a stable molecular structure is formed;
the C-S-H gel is a main product of Portland cement hydration, is a main body of a cement hydration gel phase and exists in the form of gel particles of 1-100 nm, pores among the gel particles provide channels for transmission of an erosion medium, and the hydrophilicity of the surface of the C-S-H gel provides power for capillary transmission. The C-S-H gel pore channel model is a pore channel model formed by parallelly placing C-S-H gel molecule models serving as substrates and reserving a transmission pore channel in the middle, and is called as a C-S-H gel pore channel model, and the diameter and the length of the pore channel can be adjusted by translating and extending the C-S-H substrates;
the saline solution can be a variety of aggressive media including, but not limited to: sodium chloride solution, calcium chloride solution, cesium chloride solution, sodium sulfate solution and magnesium sulfate solution;
the construction of the saline solution model is to mix a certain amount of anions and cations with water molecules so as to obtain a saline solution with a certain concentration;
the octadecanoic carboxylic acid can be dispersed in ethanol in a nanoscale, one end of the molecular structure of the octadecanoic carboxylic acid is hydrophilic, the other end of the molecular structure of the octadecanoic carboxylic acid is hydrophobic, and the hydrophilic end of the molecular structure of the octadecanoic carboxylic acid can interact with a C-S-H gel pore channel and is adsorbed on the pore channel. The octadecanoic carboxylic acid model has the same structure as the actual molecule, the constructed model C-S-H gel pore canal is perpendicular to the salt solution, and a layer of octadecanoic carboxylic acid molecules is attached to the C-S-H gel pore canal.
(2) Applying force to each atom in the established molecular model, wherein the force field describes potential energy for driving the atom to move in the system, converting the potential energy into acceleration of the movement of each atom through Newton's second law, and under the condition of giving an initial position and speed, the coordinate of each atom at each time point can be calculated, and the force field adopts an empirical force field, and the fundamental form of the empirical force field is shown as formula 1 and comprises bond energy among atoms, van der Waals force action and long coulomb force action;
Etotal=Ebond+Eangle+Etorsion+Evdw+Ecoulequation 1
The expansion formula is as follows:
in the formula 1, EtotalRepresents the total potential energy, EbondIndicating key expansion potential energy, EangleIndicating the bending potential energy of the key angle, EtorsionRepresenting dihedral torsional potential energy, EvdwDenotes van der Waals' force, EcoulRepresenting the coulomb electrostatic potential.
In the expansion formula, the first term, the second term, the third term and the fourth term on the right respectively represent key expansion potential energy, key angle bending potential energy, dihedral angle twisting potential energy and non-key energy. Specifically, k1i、k2i、k3iRepresenting the elastic constant of the key expansion term; Δ l represents the bond length; Δ θ represents a bond angle; ofijAnd σijRepresents a potential energy parameter (epsilon)ijPotential energy difference between reaction atoms, sigmaijReflecting equilibrium distance between atoms); r isijRepresents an interatomic distance; q. q.si(j)Represents the charge of the ith (j) ion in the molecule; cos (nw-r) represents the dihedral angle parameter; epsilon0Representing a potential energy parameter.
(3) The method comprises the steps of giving preset environment temperature, calculation time step and periodic boundary conditions, ensuring system stability in a simulation process by using a Nose-Hoover constant temperature calculation method, and obtaining an atomic motion trajectory by using a Verlet algorithm;
(4) and (4) according to the atom motion coordinate track obtained in the step (3), obtaining the positions and the penetration depths of the ions in the pore channels at different times, drawing a correlation curve of the penetration depths and the time, and fitting the obtained correlation curve of the penetration depths and the time to obtain the diffusion coefficient.
When the method is adopted to carry out actual simulation calculation, a molecular dynamics model is established firstly, wherein C-S-H gel is usedThe tobermorite crystal substitution is shown in fig. 2, and the crystallographic parameters of the cells are as follows: α ═ β ═ γ ═ 90 °, including 43 calcium atoms, 24 silicon-oxygen tetrahedrons, and 16 water molecules. The C-S-H gel pore channel is formed by two C-S-H substrates which are arranged in parallel, and the width of the constructed pore channel is 3.4nm and the length is 7 nm; saline solution model: mixing 50 sulfate ions and 50 sodium ions with water molecules to obtain a sodium sulfate solution of 0.6 mol/L; octadecanecarboxylic acid model: the octadecanecarboxylic acid model has the same structure as the actual molecule, and the molecular structural formula is [ CH ]2]18COOH. Octadecanecarboxylic acid molecules are attached to the C-S-H pore canals, and the attachment rate is 25%.
Establishing a molecular dynamics model, applying a force field on each atom, and performing simulation calculation on LAMMPS, wherein the system temperature is 300K, the pressure is one atmosphere, and the time step is 1 fs. The periodic boundaries are set at X, Y, Z in three directions. The initial velocity of each atom is randomly generated based on the initial temperature. The position of the next moment of the atom is calculated using the Verlet algorithm with a long-range force cutoff radius ofThe simulation process is divided into three steps. Firstly, energy optimization is carried out on the structure through energy minimization; then, a barrier is arranged between the solution and the C-S-H substrate, so that the solution can not enter the pore channel temporarily, and meanwhile, the whole system moves for 1000ps under a normal synthetic method (NVT) to achieve the balance of each part of structure; finally, under the regular system, the barrier is deleted and dissolvedThe liquid enters the pore channel under the action of capillary force and moves 1000 ps. In the simulation process, thermodynamic parameters such as temperature, pressure and the like and the size of the box are output every 1ps, coordinates of all atoms are output once, and the number of atomic coordinate frames is output for 1000 frames in total. The process of solution transmission in the C-S-H pore channel is shown in fig. 3(a), and the ion penetration depth at different times is calculated according to the atomic coordinates, and the obtained graph of the relation between the penetration depth and the time is obtained by an exponential equation (y ═ x)D+ b) fitting the curve of penetration depth and time to obtain the diffusion coefficient D, wherein b: a regulatory factor; y: the diffusion depth; x: time. And the diffusion coefficient D is used as a characterizing parameter for the transmission rate.
In order to further compare the effect of different amounts of octadecanoic carboxylic acid on the transmission efficiency, the relationship between the penetration depth and time was determined by the method of the present invention for 50%, 75%, and 100% of octadecanoic carboxylic acid, and compared with the blank case, the result is shown in fig. 3(b), from which it can be seen that the ion transmission inhibition effect is larger as the amount of octadecanoic carboxylic acid is increased.
Claims (6)
1. A quantification method for reducing the ion transmission rate of a cement-based material pore passage is characterized by comprising the following steps:
(1) establishing a molecular dynamics model of C-S-H gel molecules, octadecane carboxylic acid molecules and salt solution, constructing the C-S-H molecules, the octadecane carboxylic acid molecules and the salt solution molecules, and relaxing under an NVT (noise, vibration and harshness) ensemble to enable all parts of the model to be close to each other to form a stable molecular structure;
(2) applying a force to each atom in the molecular model created above, the force field being of the form shown below:
Etotal=Ebond+Eangle+Etorsion+Evdw+Ecoul
in the formula, EtotalRepresents the total potential energy, EbondIndicating key expansion potential energy, EangleIndicating the bending potential energy of the key angle, EtorsionRepresenting dihedral torsional potential energy, EvdwDenotes van der Waals' force, EcoulRepresenting electrostatic potential of coulombs;
(3) Setting preset environment temperature, calculation time step length and the outer edge of the whole molecular model as periodic boundary conditions, ensuring system stability in the simulation process by using a Nose-Hoover constant temperature calculation method, and obtaining an atom motion coordinate trajectory by using a Verlet algorithm;
(4) and (4) according to the atom motion coordinate track obtained in the step (3), obtaining the positions and the penetration depths of the ions in the pore channels at different times, drawing a correlation curve of the penetration depths and the time, and fitting the obtained correlation curve of the penetration depths and the time to obtain the diffusion coefficient.
2. The method for quantifying the reduction in the ion transfer rate of the cement-based material pore canals according to claim 1, wherein when the molecular dynamics model is established in the step (1), octadecanecarboxylic acid molecules are placed on the surface of the inner wall of the C-S-H molecules, and the salt solution molecules are distributed at the spatial positions below the C-S-H molecules.
4. The method for quantifying the reduction in the ion transport rate of the cement-based material pore channels according to claim 1, wherein the octadecanecarboxylic acid model in the step (1) is the same as the actual molecular structure.
5. The method for quantifying the reduction in the ion transmission rate of the cement-based material pore channels according to claim 1, wherein the salt solution in the step (1) is at least one of a sodium chloride solution, a calcium chloride solution, a cesium chloride solution, a sodium sulfate solution and a magnesium sulfate solution, and the salt solution model is constructed by mixing anions and cations with water molecules, wherein the ratio of the anions to the cations satisfies a chemical molecular formula, so that the salt solution with a certain concentration is obtained.
6. The method for quantifying and reducing the pore ion transport rate of cement-based materials according to claim 1, wherein the field formula of the step (2) is developed as follows:
in the formula, the first, second, third and fourth terms on the right respectively represent key expansion potential energy, key angle bending potential energy, dihedral angle twisting potential energy and non-key energy, specifically, k1i、k2i、k3iRepresenting the elastic constant of the key expansion term; Δ l represents the bond length; Δ θ represents a bond angle; ofijAnd σijRepresents a potential energy parameter (epsilon)ijPotential energy difference between reaction atoms, sigmaijReflecting equilibrium distance between atoms); r isijRepresents an interatomic distance; q. q.si(j)Represents the charge of the ith (j) ion in the molecule; cos (nw-r) represents the dihedral angle parameter; epsilon0Representing a potential energy parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110648065.7A CN113533138B (en) | 2021-06-10 | 2021-06-10 | Quantification method for reducing ion transmission rate of cement-based material pore channel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110648065.7A CN113533138B (en) | 2021-06-10 | 2021-06-10 | Quantification method for reducing ion transmission rate of cement-based material pore channel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113533138A true CN113533138A (en) | 2021-10-22 |
CN113533138B CN113533138B (en) | 2023-03-21 |
Family
ID=78124776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110648065.7A Active CN113533138B (en) | 2021-06-10 | 2021-06-10 | Quantification method for reducing ion transmission rate of cement-based material pore channel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113533138B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5405509A (en) * | 1989-05-08 | 1995-04-11 | Ionex | Remediation of a bulk source by electropotential ion transport using a host receptor matrix |
CA2349153A1 (en) * | 2001-05-24 | 2002-11-24 | Jacques Marchand | Method for modeling the transport of ions in hydrated cement systems |
CN107561100A (en) * | 2017-09-06 | 2018-01-09 | 东南大学 | A kind of method of continuous monitoring cement based porous materials Water Transport |
CN108129051A (en) * | 2017-12-29 | 2018-06-08 | 江苏苏博特新材料股份有限公司 | A kind of long-term efficiently concrete admixture of curing of chloride ion and application |
CN108333082A (en) * | 2018-01-23 | 2018-07-27 | 广西大学 | The construction method of unsaturation Chloride Diffusion Coefficient in Concrete multi-scale prediction model |
CN109872778A (en) * | 2019-01-30 | 2019-06-11 | 东南大学 | A kind of evaluation method of regenerative agent diffusion in waste asphalt mixture |
CN111310360A (en) * | 2020-03-26 | 2020-06-19 | 东南大学 | Concrete damage transmission assessment method under sulfate erosion in splash zone |
CN111899806A (en) * | 2020-06-11 | 2020-11-06 | 东南大学 | Intelligent design method of graphene oxide/C-S-H gel composite material based on molecular dynamics theory |
-
2021
- 2021-06-10 CN CN202110648065.7A patent/CN113533138B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5405509A (en) * | 1989-05-08 | 1995-04-11 | Ionex | Remediation of a bulk source by electropotential ion transport using a host receptor matrix |
CA2349153A1 (en) * | 2001-05-24 | 2002-11-24 | Jacques Marchand | Method for modeling the transport of ions in hydrated cement systems |
CN107561100A (en) * | 2017-09-06 | 2018-01-09 | 东南大学 | A kind of method of continuous monitoring cement based porous materials Water Transport |
CN108129051A (en) * | 2017-12-29 | 2018-06-08 | 江苏苏博特新材料股份有限公司 | A kind of long-term efficiently concrete admixture of curing of chloride ion and application |
CN108333082A (en) * | 2018-01-23 | 2018-07-27 | 广西大学 | The construction method of unsaturation Chloride Diffusion Coefficient in Concrete multi-scale prediction model |
CN109872778A (en) * | 2019-01-30 | 2019-06-11 | 东南大学 | A kind of evaluation method of regenerative agent diffusion in waste asphalt mixture |
CN111310360A (en) * | 2020-03-26 | 2020-06-19 | 东南大学 | Concrete damage transmission assessment method under sulfate erosion in splash zone |
CN111899806A (en) * | 2020-06-11 | 2020-11-06 | 东南大学 | Intelligent design method of graphene oxide/C-S-H gel composite material based on molecular dynamics theory |
Non-Patent Citations (5)
Title |
---|
乔宏霞等: "玄武岩石粉对C50混凝土氯离子渗透性的影响", 《非金属矿》 * |
王彩辉等: "动载-环境耦合作用下氯离子在混凝土中的扩散性能研究", 《工业建筑》 * |
王彩辉等: "水泥基复合材料在多尺度方面的研究进展", 《硅酸盐学报》 * |
蒋金洋等: "多尺度分析方法在土木工程材料教学中的探索与实践", 《高等建筑教育》 * |
钱如胜等: "水泥-粉煤灰体系早龄期液相离子浓度与电导率的关系", 《材料导报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113533138B (en) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hou et al. | Modified Lucas-Washburn function of capillary transport in the calcium silicate hydrate gel pore: A coarse-grained molecular dynamics study | |
Hou et al. | Molecular dynamics study on the structure and dynamics of NaCl solution transport in the nanometer channel of CASH gel | |
Soliman et al. | Using Particle Packing and Statistical Approach to Optimize Eco-Efficient Ultra-High-Performance Concrete. | |
Espinosa-Marzal et al. | Advances in understanding damage by salt crystallization | |
CN111310360B (en) | Concrete damage transmission assessment method under sulfate erosion in splash zone | |
Basista et al. | Micromechanical modeling of sulphate corrosion in concrete: influence of ettringite forming reaction | |
Zacchei et al. | Chloride diffusion assessment in RC structures considering the stress-strain state effects and crack width influences | |
Zhou et al. | The inhibiting effect and mechanisms of smart polymers on the transport of fluids throughout nano-channels | |
Wang et al. | Effect of temperature on the capillary transport of sodium sulfate solution in calcium silicate hydrate nanopore: A molecular dynamics study | |
Tu et al. | Molecular dynamics simulation study of the transport of pairwise coupled ions confined in CSH gel nanopores | |
Zehtab et al. | Diffusion study for chloride ions and water molecules in CSH gel in nano-scale using molecular dynamics: Case study of tobermorite | |
Zhou et al. | The design and evaluation of a smart polymer-based fluids transport inhibitor | |
Yaphary et al. | Mechanical properties of colloidal calcium-silicate-hydrate gel with different gel-pore ionic solutions: A mesoscale study | |
Chen et al. | A multi-phase mesoscopic simulation model for the long-term chloride ingress and electrochemical chloride extraction | |
Liang et al. | Freezing mechanism of NaCl solution ultra-confined on surface of calcium-silicate-hydrate: A molecular dynamics study | |
CN113533138B (en) | Quantification method for reducing ion transmission rate of cement-based material pore channel | |
Hou et al. | Molecular-scale insights on structure-efficiency relationship of silane-based waterproofing agents | |
Li et al. | Surface modification strategy for controlling wettability and ionic diffusion behaviors of calcium silicate hydrate | |
Luo et al. | Molecular simulation of calcium-silicate-hydrate and its applications: a comprehensive review | |
CN117198424B (en) | Concrete diffusion prediction method under sulfate-chloride coupling erosion | |
Hayat et al. | Atomic-level investigation into the transport of NaCl solution in porous cement paste: The effects of pore size and temperature | |
Chen et al. | Molecular dynamics study on the transport of water molecules and chloride ions in graphene oxide-modified cement composites | |
Liu et al. | Water transport in calcium silicate hydrate nanoscale pores by molecular dynamics: Unsaturated steady-state diffusion and suction | |
Wan et al. | Molecular dynamics simulations of the interface friction behavior between fiber-reinforced polymer pile and sand | |
Ferjaoui | Nanoscale modelling of ionic transport in the porous CSH network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |