Method and device for calculating porosity and tortuosity and electronic equipment
By constructing a three-dimensional porous dielectric material structure and reconstructing, the problems of unstable and cost of calculating porosity and tortuosity in the prior art are solved, and efficient and accurate calculation results are achieved.
Patent Information
- Application Number
- CN202411927058.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-25
- Publication Date
- 2025-04-15
AI Technical Summary
When calculating porosity and tortuosity in the prior art, there is an unclear empirical constant that leads to unstable results. The geometric analytical method lacks a global microscopic mechanism, and the three-dimensional CT scanning method is costly and difficult to widely use.
By obtaining solid particle structure information, spherical spherical particles and their stacking aggregates are constructed, a three-dimensional porous dielectric material structure is obtained, the porosity is determined, and a smooth interface structure is generated through reconstruction, and the tortuosity is calculated.
Improve the accuracy and efficiency of calculations, reduce costs, and quickly obtain stable porosity and tortuosity calculation results without high experimental equipment.
Smart Images
Figure CN119833045A_ABST
Abstract
Description
Technical Field
[0001] The present disclosure relates to the field of computer technology, and in particular to a method, device and electronic device for calculating porosity and tortuosity. Background Art
[0002] Porosity is the percentage of the pore volume inside a material to its total volume, and tortuosity is the ratio of the actual length of the seepage channel to the apparent length through the seepage medium. In order to calculate porosity and tortuosity, the model method used in the relevant technology has many empirical constants when calculating. The physical meaning of these constants is unclear, which can easily lead to unstable results and large errors from the true values; the geometric analytical method model is mostly derived from macro formulas, and rarely considers the global microscopic mechanism inside the porous medium. Therefore, it is too far from the actual situation, and the tortuosity calculation expression is not easy to generalize; and the time cost and economic cost of the three-dimensional CT scanning method are high, making it difficult to be widely used. Summary of the invention
[0003] In order to overcome the problems existing in the related art, the present disclosure provides a method, device and electronic device for calculating porosity and tortuosity.
[0004] According to a first aspect of an embodiment of the present disclosure, a method for calculating porosity and tortuosity is provided, comprising:
[0005] Obtain solid particle structure information;
[0006] According to the solid particle structure information, spherical ball particles and their stacking aggregates are constructed to obtain a three-dimensional porous medium material structure;
[0007] Determining the porosity of the porous medium material according to the three-dimensional porous medium material structure;
[0008] According to the three-dimensional porous medium material structure, the three-dimensional porous medium material structure is reconstructed to generate a smooth interface structure;
[0009] The tortuosity of the three-dimensional porous medium material structure is determined according to the porosity of the porous medium material and the smooth interface structure.
[0010] In one embodiment, the solid particle structure information includes the type of spherical particles, the diameter of the spherical particles, and the mass fraction of the spherical particles.
[0011] In one embodiment, determining the tortuosity of the three-dimensional porous medium material structure includes:
[0012] Determine the intrinsic conductivity of the liquid phase based on the electric potential at a position where the current and material length are L;
[0013] Determine the liquid phase equivalent conductivity according to a first potential at a position where the material length is 0, a second potential at a position where the material length is L, a distance between the position where the material length is 0 and the position where the material length is L, and a current;
[0014] The tortuosity is determined according to the liquid phase equivalent conductivity, the liquid phase intrinsic conductivity and the porosity.
[0015] In one embodiment, the tortuosity of the three-dimensional porous medium material structure is determined using the following calculation formula:
[0016]
[0017] Among them, σ bulk is the intrinsic conductivity of the liquid phase, is the potential at the position where the material length is 0, I is the current, σ eff is the liquid equivalent conductivity, is the electric potential at a position where the material length is L.
[0018] In one embodiment, determining the tortuosity of the three-dimensional porous medium material structure includes:
[0019] Determine the liquid phase diffusion coefficient based on the convection coefficient, liquid phase concentration, ambient concentration and direction vector;
[0020] Determine an equivalent diffusion coefficient according to a liquid concentration at a position where the material length is 0, a liquid concentration at a position where the material length is L, a distance between the position where the material length is 0 and the position where the material length is L, and a flux;
[0021] The tortuosity is determined based on the equivalent diffusion coefficient, the liquid phase diffusion coefficient and the porosity.
[0022] In one embodiment, the tortuosity of the three-dimensional porous medium material structure is determined using the following calculation formula:
[0023]
[0024] c| x＝0 = c0;
[0025]
[0026] Among them, D eff is the equivalent diffusion coefficient, c out is the liquid concentration at the position where the material length is L, J is the flux, c is the liquid concentration, t is the time, D bulk is the liquid diffusion coefficient, c0 is the liquid concentration at the position where the material length is 0, is the direction vector, k is the convection coefficient, c envis the ambient concentration; ε is the porosity, and τ is the tortuosity.
[0027] In one embodiment, reconstructing the three-dimensional porous medium material structure includes:
[0028] The boundary is smoothed by the distance function Ψ, which is calculated using the following formula:
[0029]
[0030] Among them, d is the distance between the centers of the two spheres, R is the radius of the sphere, and A is a fixed coefficient.
[0031] In one embodiment, the three-dimensional porous medium material structure is reconstructed using a boundary smoothing method, including:
[0032] The boundary is smoothed by the phase field function, and the calculation formula of the phase field function is as follows:
[0033]
[0034] Among them, Ψ is the distance function, t is the time, ω is the potential well function, and κ is the gradient energy coefficient.
[0035] In a second aspect, the present application also proposes a porosity and tortuosity calculation device, comprising:
[0036] An acquisition module, used to obtain solid particle structure information;
[0037] A three-dimensional construction module is used to construct spherical small ball particles and their stacking aggregates according to the solid particle structure information to obtain a three-dimensional porous medium material structure;
[0038] A first determination module is used to determine the porosity of the porous medium material according to the three-dimensional porous medium material structure;
[0039] A reconstruction module, used to reconstruct the three-dimensional porous medium material structure according to the three-dimensional porous medium material structure to generate a smooth interface structure;
[0040] The second determination module is used to determine the tortuosity of the three-dimensional porous medium material structure according to the porosity of the porous medium material and the smooth interface structure.
[0041] The second determination module is also used to determine the tortuosity of the three-dimensional porous medium material structure using the following calculation formula:
[0042]
[0043]
[0044] Among them, σ bulk is the intrinsic conductivity of the liquid phase, is the potential at the position where the material length is 0, I is the current, σ eff is the liquid equivalent conductivity, is the electric potential at a position where the material length is L.
[0045] The second determination module is further used to determine the liquid phase diffusion coefficient according to the convection coefficient, the liquid phase concentration, the ambient concentration and the direction vector;
[0046] Determine the equivalent diffusion coefficient according to the liquid concentration at the position where the material length is 0 in the material, the liquid concentration at the position where the material length is L, the distance L and the flux;
[0047] The tortuosity is determined based on the above-mentioned equivalent diffusion coefficient, liquid phase diffusion coefficient and porosity.
[0048] The second determination module is also used to determine the tortuosity of the three-dimensional porous medium material structure using the following calculation formula:
[0049]
[0050] c| x＝0 = c0;
[0051]
[0052] Among them, D eff is the equivalent diffusion coefficient, c out is the liquid concentration at the position where the material length is L, J is the flux, c is the liquid concentration, t is the time, D bulk is the liquid diffusion coefficient, c0 is the liquid concentration at the position where the material length is 0, is the direction vector, k is the convection coefficient, c env is the ambient concentration; ε is the porosity, and τ is the tortuosity.
[0053] The reconstruction module is also used to perform boundary smoothing using a distance function Ψ, which is calculated using the following formula:
[0054]
[0055] Among them, d is the distance between the centers of the two spheres, R is the radius of the sphere, and A is a fixed coefficient.
[0056] The reconstruction module is also used to perform boundary smoothing through the phase field function. The calculation formula of the above phase field function is as follows:
[0057]
[0058] Among them, Ψ is the distance function, t is the time, ω is the potential well function, and κ is the gradient energy coefficient.
[0059] In a third aspect, the present application also proposes an electronic device comprising at least one processor and at least one memory; the memory is used to store one or more program instructions; the processor is used to run one or more program instructions to execute the porosity and tortuosity calculation method described in any one of the above items.
[0060] According to a fourth aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium contains one or more program instructions, and the one or more program instructions are used to execute any of the above-mentioned methods for calculating porosity and tortuosity.
[0061] The technical solution disclosed herein may have the following beneficial effects:
[0062] The present application constructs spherical ball particles and their stacked aggregates based on the solid particle structure information to obtain a three-dimensional porous medium material structure; determines the porosity of the porous medium material based on the three-dimensional porous medium material structure; and then calculates the tortuosity. Compared with related technologies, this method saves costs and improves accuracy.
[0063] It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present disclosure. BRIEF DESCRIPTION OF THE DRAWINGS
[0064] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and, together with the description, serve to explain the principles of the present disclosure.
[0065] Figure 1 is a flow chart showing a method for calculating porosity and tortuosity according to an exemplary embodiment;
[0066] Figure 2 is a three-dimensional porous medium material structure diagram according to an exemplary embodiment;
[0067] Figure 3 is a diagram showing a smooth interface structure according to an exemplary embodiment;
[0068] Figure 4 is a block diagram of a porosity and tortuosity calculation device according to an exemplary embodiment;
[0069] Figure 5 It is a block diagram of an electronic device according to an exemplary embodiment. DETAILED DESCRIPTION
[0070] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.
[0071] The terms "first" or "second" and the like used in this specification to refer to numbers or ordinals are used for descriptive purposes only and are not to be understood as explicitly or implicitly indicating relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of this specification, "plurality" means at least two, such as two, three or more, etc., unless otherwise clearly and specifically defined.
[0072] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.
[0073] This application proposes a method for calculating porosity and tortuosity, see Appendix Figure 1 , the method may include the following steps:
[0074] In step S102, solid particle structure information is obtained.
[0075] In this embodiment, the solid particle structure information includes the type of spherical particles, the diameter of the spherical particles, and the mass fraction of the spherical particles.
[0076] Construct the structure.csv file. The structure.csv file contains the above-mentioned solid particle structure information. Specifically, the structure.csv file contains 3 columns of solid particle structure data. The first column of solid particle structure is the type of small ball particles, and the name of the small ball particle type can be set arbitrarily. The second column of solid particle structure data is the diameter of the small ball particles, which is input in microns (μm); the third column of solid particle structure data is the mass fraction of small ball particles. The input example of structure data of the structure.csv file is shown in Table 1.
[0077] Particle Type Diameter(μm) Quality score A D G B E H C F I
[0078] In step S104, spherical beads and their stacked aggregates are constructed according to the above solid particle structure information to obtain a three-dimensional porous medium material structure.
[0079] In this embodiment, see the attached Figure 2 According to the solid particle structure information in the structure.csv file above, the discrete element method is used to construct spherical beads and their stacking aggregates. The spherical beads are randomly and irregularly arranged to form a three-dimensional porous medium material structure.
[0080] In step S106, the porosity of the porous medium material is determined according to the three-dimensional porous medium material structure.
[0081] In this embodiment, according to the three-dimensional porous medium material structure, the porosity of the porous medium material is determined by a statistical calculation method. The porosity calculation formula is shown in Formula 1:
[0082]
[0083] Where ε is the porosity, V s is the solid phase volume, V Ω is the total volume of the constructed three-dimensional porous medium material structure.
[0084] In step S108, the three-dimensional porous medium material structure is reconstructed according to the three-dimensional porous medium material structure to generate a smooth interface structure.
[0085] In this embodiment, see the attached Figure 3 According to the above three-dimensional porous medium material structure, the boundary smoothing method is used to reconstruct the structure and generate a smooth interface structure. There are two schemes for the boundary smoothing method, namely the distance function method and the phase field function method.
[0086] In step S110, the tortuosity of the three-dimensional porous medium material structure is determined according to the porosity of the porous medium material and the smooth interface structure.
[0087] In this embodiment, the structural tortuosity of the three-dimensional porous medium material is calculated based on the porosity of the material structure and the smooth interface structure.
[0088] The above-mentioned technical solution of the present invention obtains the three-dimensional porous medium material structure according to the solid particle structure information. The porosity of the porous medium material is determined according to the above-mentioned three-dimensional porous medium material structure. The tortuosity of the above-mentioned three-dimensional porous medium material structure is further determined according to the porosity of the above-mentioned porous medium material. Compared with the related art, the above-mentioned technical solution of the present application is not only conducive to improving the calculation efficiency and the accuracy of the calculation, but also quickly obtaining the three-dimensional structure, porosity and tortuosity of the material, which is convenient for characterizing the macroscopic transport performance of the material and guiding production development. It is also conducive to reducing the calculation cost, and stable porosity and tortuosity calculation results can be quickly obtained without expensive experimental equipment.
[0089] In some embodiments, in step S110, the determining the tortuosity of the three-dimensional porous medium material structure may further include the following steps:
[0090] The intrinsic conductivity of the liquid phase is determined by the potential at a position where the length L of the material is determined based on the current.
[0091] In this embodiment, any two positions in the three-dimensional porous medium material structure can be determined, and the tortuosity can be determined based on the potential of the two positions and the current between the two positions. Specifically, the liquid phase intrinsic conductivity is obtained by dividing the current by the change in the potential at the position where the material length is L.
[0092] The liquid phase equivalent conductivity is determined according to a first potential at a position in the material where the material length is 0, a second potential at a position where the material length is L, the distance L and the current.
[0093] In this embodiment, the potential difference between the first potential at the position where the material length is 0 and the second potential at the position where the material length is L can be calculated; the distance L between the position where the material length is 0 and the position where the material length is L is calculated; the above distance L is divided by the above potential difference, and multiplied by the current I to obtain the liquid phase equivalent conductivity.
[0094] The tortuosity is determined according to the liquid phase equivalent conductivity, the liquid phase intrinsic conductivity and the porosity.
[0095] In this embodiment, the tortuosity can be obtained by multiplying the liquid phase intrinsic conductivity by the porosity and dividing by the liquid phase equivalent conductivity.
[0096] In some embodiments, in step S110, the tortuosity of the three-dimensional porous medium material structure is determined using the following calculation formula:
[0097]
[0098]
[0099] Among them, σ bulk is the intrinsic conductivity of the liquid phase, is the potential at the position where the material length is 0, I is the current, σ eff is the liquid equivalent conductivity, is the electric potential at a position where the material length is L.
[0100] In some embodiments, in step S110, the determining the tortuosity of the three-dimensional porous medium material structure may further include the following steps:
[0101] The liquid phase diffusion coefficient is determined from the convection coefficient, liquid phase concentration, ambient concentration, and direction vector.
[0102] The equivalent diffusion coefficient is determined according to the liquid concentration at a position where the material length is 0, the liquid concentration at a position where the material length is L, the distance L between the position where the material length is 0 and the position where the material length is L, and the flux.
[0103] The tortuosity is determined based on the equivalent diffusion coefficient, the liquid phase diffusion coefficient and the porosity.
[0104] In some embodiments, the tortuosity of the three-dimensional porous medium material structure is determined using the following calculation formula:
[0105]
[0106] c| x＝0 = c0;
[0107]
[0108] Among them, D eff is the equivalent diffusion coefficient, c out is the liquid concentration at the position where the material length is L, J is the flux, c is the liquid concentration, t is the time, D bulk is the liquid diffusion coefficient, c0 is the liquid concentration at the position where the material length is 0, is the direction vector, k is the convection coefficient, c env is the ambient concentration; ε is the porosity, and τ is the tortuosity.
[0109] In some embodiments, reconstructing the three-dimensional porous medium material structure may include the following steps:
[0110] The boundary is smoothed by the distance function Ψ, which is calculated using the following formula:
[0111]
[0112] Among them, d is the distance between the centers of the two spheres, R is the radius of the sphere, and A is a fixed coefficient, which can be set to 4.8.
[0113] In some embodiments, reconstructing the three-dimensional porous medium material structure using a boundary smoothing method may include the following steps:
[0114] The boundary is smoothed by the phase field function, and the calculation formula of the phase field function is as follows:
[0115]
[0116] Among them, Ψ is the distance function, t is the time, ω is the potential well function, and κ is the gradient energy coefficient.
[0117] The present application has the following advantages: 1) Wide scope of application: For any material type, the present invention only needs to specify the diameter and mass fraction of the small spherical particles that constitute the three-dimensional structure to generate a three-dimensional porous medium structure and calculate the porosity and tortuosity. 2) Efficient calculation and stable results: For the three-dimensional porous structure that has been generated, multiple simulations and the use of two different numerical algorithms can quickly obtain porosity and tortuosity results with basically consistent values. 3) Low research cost: Stable porosity and tortuosity calculation results can be quickly obtained without expensive experimental equipment. 4) Easy for researchers to understand and master: The present invention calculates porosity and tortuosity with the characteristics of standardization and process, and does not require researchers to have too much professional knowledge to obtain the required porosity and tortuosity results.
[0118] Based on the same inventive concept, the present application also proposes a porosity and tortuosity calculation device 2, see the attached Figure 4 A schematic diagram of a porosity and tortuosity calculation device is shown; the device comprises:
[0119] An acquisition module 21 is used to acquire solid particle structure information;
[0120] A three-dimensional construction module 22 is used to construct spherical small ball particles and their stacking aggregates according to the above solid particle structure information to obtain a three-dimensional porous medium material structure;
[0121] A first determination module 23 is used to determine the porosity of the porous medium material according to the three-dimensional porous medium material structure;
[0122] A reconstruction module 24 is used to reconstruct the three-dimensional porous medium material structure according to the three-dimensional porous medium material structure to generate a smooth interface structure;
[0123] The second determination module 25 is used to determine the tortuosity of the three-dimensional porous medium material structure according to the porosity of the porous medium material and the smooth interface structure.
[0124] The second determination module 25 is further used to determine the liquid phase intrinsic conductivity according to the current and the potential at the position where the material length is L; determine the liquid phase equivalent conductivity according to the first potential at the position where the material length is 0, the second potential at the position where the material length is L, the distance L between the position where the material length is 0 and the position where the material length is L, and the current;
[0125] The tortuosity is determined based on the above-mentioned liquid phase equivalent conductivity, liquid phase intrinsic conductivity and porosity.
[0126] The second determination module 25 is further used to determine the tortuosity of the three-dimensional porous medium material structure using the following calculation formula:
[0127]
[0128] Among them, σ bulk is the intrinsic conductivity of the liquid phase, is the potential at the position where the material length is 0, I is the current, σ eff is the liquid equivalent conductivity, is the electric potential at a position where the material length is L.
[0129] The second determination module 25 is further used to determine the liquid phase diffusion coefficient according to the convection coefficient, the liquid phase concentration, the ambient concentration and the direction vector;
[0130] Determine the equivalent diffusion coefficient according to the liquid concentration at the position where the material length is 0, the liquid concentration at the position where the material length is L, the distance between the position where the material length is 0 and the position where the material length is L, and the flux;
[0131] The tortuosity is determined based on the above-mentioned equivalent diffusion coefficient, liquid phase diffusion coefficient and porosity.
[0132] The second determination module 25 is further used to determine the tortuosity of the three-dimensional porous medium material structure using the following calculation formula:
[0133]
[0134] c| x＝0 = c0;
[0135]
[0136] Among them, D eff is the equivalent diffusion coefficient, c out is the liquid concentration at the position where the material length is L, J is the flux, c is the liquid concentration, t is the time, D bulk is the liquid diffusion coefficient, c0 is the liquid concentration at the position where the material length is 0, is the direction vector, k is the convection coefficient, c env is the ambient concentration; ε is the porosity, and τ is the tortuosity.
[0137] The reconstruction module 24 is also used to perform boundary smoothing using a distance function Ψ. The distance function Ψ is calculated using the following formula:
[0138]
[0139] Among them, d is the distance between the centers of the two spheres, R is the radius of the sphere, and A is a fixed coefficient.
[0140] The reconstruction module 24 is also used to perform boundary smoothing by using a phase field function. The calculation formula of the phase field function is as follows:
[0141]
[0142] Among them, Ψ is the distance function, t is the time, ω is the potential well function, and κ is the gradient energy coefficient.
[0143] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.
[0144] According to the third aspect of this application, see the attached Figure 5 , an electronic device is provided, comprising at least one processor 31 and at least one memory 32. The memory 32 is used to store one or more program instructions; the processor 31 is used to run one or more program instructions to execute any one of the above-mentioned methods for calculating porosity and tortuosity.
[0145] For example Figure 1 Alternatively, when the processor 30 executes the computer program 32, the functions of the modules in the above-mentioned device embodiments are realized, for example, Figure 4 The functions of the acquisition module 21 to the second determination module 25 are shown.
[0146] Exemplarily, the computer program 32 may be divided into one or more modules / units, which are stored in the memory 31 and executed by the processor 30 to complete the present application. The one or more modules / units may be a series of computer program instruction segments capable of completing specific functions, which are used to describe the execution process of the computer program 32 in the electronic device 3.
[0147] The electronic device 3 may be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The electronic device may include, but is not limited to, a processor 30 and a memory 31. Those skilled in the art will appreciate that Figure 5 It is only an example of the electronic device 3 and does not constitute a limitation of the electronic device 3. It may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the above-mentioned electronic device may also include input and output devices, network access devices, buses, etc.
[0148] The processor 30 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.
[0149] The memory 31 may be an internal storage unit of the electronic device 3, such as a hard disk or memory of the electronic device 3. The memory 31 may also be an external storage device of the electronic device 3, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the electronic device 3. Further, the memory 31 may include both an internal storage unit and an external storage device of the electronic device 3. The memory 31 is used to store the computer program and other programs and data required by the electronic device. The memory 31 may also be used to temporarily store data that has been output or is to be output.
[0150] In a fourth aspect, the present application further proposes a computer-readable storage medium, which contains one or more program instructions, and the one or more program instructions are used to execute the porosity and tortuosity calculation method described in any one of the above items.
[0151] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims
1. A method for calculating porosity and tortuosity, characterized in that: include: Obtain solid particle structure information; According to the solid particle structure information, spherical ball particles and their stacking aggregates are constructed to obtain a three-dimensional porous medium material structure; Determining the porosity of the porous medium material according to the three-dimensional porous medium material structure; According to the three-dimensional porous medium material structure, the three-dimensional porous medium material structure is reconstructed to generate a smooth interface structure; The tortuosity of the three-dimensional porous medium material structure is determined according to the porosity of the porous medium material and the smooth interface structure.
2. The method for calculating porosity and tortuosity according to claim 1, characterized in that: The solid particle structure information includes the type of spherical particles, the diameter of the spherical particles and the mass fraction of the spherical particles.
3. The method for calculating porosity and tortuosity according to claim 1, characterized in that: Determining the tortuosity of the three-dimensional porous medium material structure includes: Determine the intrinsic conductivity of the liquid phase based on the electric potential at a position where the current and material length are L; Determine the liquid phase equivalent conductivity according to the potential at a position where the material length is 0, the potential at a position where the material length is L, the distance L and the current; The tortuosity is determined according to the liquid phase equivalent conductivity, the liquid phase intrinsic conductivity and the porosity.
4. The method for calculating porosity and tortuosity according to claim 3, characterized in that: The tortuosity of the three-dimensional porous medium material structure is determined using the following calculation formula: Among them, σ bulk is the intrinsic conductivity of the liquid phase, is the potential at the position where the material length is 0, I is the current, σ eff is the liquid equivalent conductivity, is the electric potential at a position where the material length is L.
5. The method for calculating porosity and tortuosity according to claim 1, characterized in that: Determining the tortuosity of the three-dimensional porous medium material structure includes: Determine the liquid phase diffusion coefficient based on the convection coefficient, liquid phase concentration, ambient concentration and direction vector; Determine the equivalent diffusion coefficient according to the liquid concentration at the position where the material length is 0 in the material, the liquid concentration at the position where the material length is L, the distance L and the flux; The tortuosity is determined based on the equivalent diffusion coefficient, the liquid phase diffusion coefficient and the porosity.
6. The method for calculating porosity and tortuosity according to claim 5, characterized in that: The tortuosity of the three-dimensional porous medium material structure is determined using the following calculation formula: c| x＝0 =c0; Among them, D eff is the equivalent diffusion coefficient, c out is the liquid concentration at the position where the material length is L, J is the flux, c is the liquid concentration, t is the time, D bulk is the liquid diffusion coefficient, c0 is the liquid concentration at the position where the material length is 0, is the direction vector, k is the convection coefficient, c env is the ambient concentration; ε is the porosity, and τ is the tortuosity.
7. The method for calculating porosity and tortuosity according to claim 1, characterized in that: Reconstructing the three-dimensional porous medium material structure includes: The boundary is smoothed by the distance function Ψ, which is calculated using the following formula: Among them, d is the distance between the centers of the two spheres, R is the radius of the sphere, and A is a fixed coefficient.
8. The method for calculating porosity and tortuosity according to claim 1, characterized in that: The three-dimensional porous medium material structure is reconstructed using a boundary smoothing method, including: The boundary is smoothed by the phase field function, and the calculation formula of the phase field function is as follows: Among them, Ψ is the distance function, t is the time, ω is the potential well function, and k is the gradient energy coefficient.
9. A device for calculating porosity and tortuosity, characterized in that: include: An acquisition module, used to obtain solid particle structure information; A three-dimensional construction module is used to construct spherical small ball particles and their stacking aggregates according to the solid particle structure information to obtain a three-dimensional porous medium material structure; A first determination module is used to determine the porosity of the porous medium material according to the three-dimensional porous medium material structure; A reconstruction module, used to reconstruct the three-dimensional porous medium material structure according to the three-dimensional porous medium material structure to generate a smooth interface structure; The second determination module is used to determine the tortuosity of the three-dimensional porous medium material structure according to the porosity of the porous medium material and the smooth interface structure.
10. An electronic device, characterized in that: The method comprises at least one processor and at least one memory; the memory is used to store one or more program instructions; the processor is used to run one or more program instructions to execute the method according to any one of claims 1 to 8.