Optimization Design Method, Medium and Equipment for Capillary Groove Cross-section Based on LBM

Abstract

The present invention relates to the technical field of optimized design of capillary groove cross-sections, and relates to an optimized design method for capillary groove cross-sections based on LBM, including: S1: constructing a numerical model; obtaining the initial values of capillary groove cross-section parameters and the initial cross-section shape; S2: calculating the reduced area of the capillary groove; S3: if the reduced area is smaller than that of the previous generation, adjust any one of the capillary groove cross-section parameters, and return to S2; otherwise, enter S4; S4: draw the cross-section image of the capillary groove according to the capillary groove cross-section parameters of this generation, and use image recognition technology to obtain a 0-1 matrix reflecting the boundary information of the capillary groove; S5: import the matrix into the numerical model to obtain the time taken for water to fill the entire capillary groove; S6: calculate the water absorption efficiency of the capillary groove. If the water absorption efficiency of this generation is smaller than that of the previous generation or smaller than Q, adjust the width of the water absorption slot or the width of the diversion groove, and return to S2; if the water absorption efficiency of this generation is greater than that of the previous generation and greater than Q, an optimized capillary groove cross-section with better drainage effect is obtained.

Description

Technical Field

[0001] The present invention relates to the technical field of optimized design of capillary groove cross-sections, and particularly relates to an optimized design method for capillary groove cross-sections based on LBM. Background Art

[0002] Currently, drainage components that utilize capillary action for active drainage are various drainage belts, drainage pipes, and water nail pipes produced based on patents with patent numbers CN971 15515.1 and CN99126199.2. Their technical characteristics are that, by means of the tiny capillary slots densely distributed on the outer surface of the components, water in the soil is sucked into the diversion groove through the water absorption slots under the action of capillary force, and then, through the synergistic action of gravity and siphonage, the water is discharged from the soil through the internal drainage pipe or directly discharged from the soil.

[0003] The patent with patent number CN202020131255.2 has optimized the cross-sectional shape of the capillary groove part of this component, designed the diversion groove into a water droplet shape, and as a result, enhanced the capillary action of the component and increased the drainage volume. However, the prior art has not optimized parameters such as the cross-sectional shape and size of the capillary groove, the width and depth of the slot opening, etc. Compared with the existing capillary grooves, it fails to achieve better water absorption, drainage, or water permeability effects, and there are certain limitations. Summary of the Invention

[0004] The purpose of the present invention is to provide an optimized design method for capillary groove cross-sections based on LBM that can design capillary groove cross-sections with better drainage effects. The specific technical solutions are as follows:

[0005] The present invention provides an optimized design method for capillary groove cross-sections based on LBM, including the following steps:

[0006] S1: Obtain the reduced area A of the existing capillary groove 0 and the water absorption efficiency η of the capillary groove 0 ; construct a numerical model based on LBM;

[0007] Obtain the initial values and the initial cross-sectional shape of the capillary groove cross-sectional parameters, where: the capillary groove cross-sectional parameters include the water absorption slot width, the water absorption slot length, the diversion groove width, and the diversion groove length;

[0008] S2: Calculate the reduced area A of the capillary groove according to the capillary groove cross-sectional parameters k , k is the number of iterations, which takes a natural number greater than or equal to 1;

[0009] S3: If A k < A k-1 , then let k = k + 1, and adjust any one of the water absorption slot width, the water absorption slot length, the diversion groove width, and the diversion groove length in the capillary groove cross-sectional parameters using the first adjustment formula, and return to S2; if Ak ≥A k-1 , then proceed to S4;

[0010] The first adjustment formula includes: b k = b k-1 + 2n; d k = d k-1 - 2n; h k = h k-1 + 2n; l k = l k-1 + 2n; where: b k is the width of the water absorption slit in the k-th iteration; d k is the length of the water absorption slit in the k-th iteration, h k is the width of the flow guiding groove in the k-th iteration; l k is the length of the flow guiding groove in the k-th iteration; n is the grid length in the numerical model;

[0011] S4: Draw the cross-sectional image of the capillary groove according to the cross-sectional parameters of the capillary groove in the k-th iteration, and use image recognition technology to binarize the cross-sectional image of the capillary groove to obtain a 0-1 matrix reflecting the boundary information of the capillary groove;

[0012] S5: Import the 0-1 matrix reflecting the boundary information of the capillary groove into the numerical model to simulate the rising process of water in the capillary groove, and obtain the time t E ;

[0013] S6: Calculate the water absorption efficiency η E of the capillary groove according to t k . If η k < η k-1 or η k < Q, then let k = k + 1, and adjust the width of the water absorption slit or the width of the flow guiding groove using the second adjustment formula, and return to S2; If η k ≥ η k-1 and η k ≥ Q, then obtain the optimized cross-section of the capillary groove; where: Q is the set threshold;

[0014] The second adjustment formula includes: b k = b k-1 - n; h k = h k-1 - n.

[0015] Optionally, in S3, the specific calculation formula for the reduced area A k of the capillary groove is as follows:

[0016]

[0017] Wherein: A is the cross-sectional area of the capillary groove, h k is the width of the diversion groove, and s is the spacing of the diversion grooves;

[0018] Among any four adjacent adjustments, the cross-sectional parameters of the capillary grooves selected are all different.

[0019] Optionally, in step S4, using image recognition technology to binarize the cross-sectional image of the capillary groove to obtain a 0-1 matrix reflecting the boundary information of the capillary groove specifically includes:

[0020] Obtaining the image threshold using the graythresh function;

[0021] According to the image threshold, using the im2bw function to convert the cross-sectional image into a binary image, obtaining a 0-1 matrix reflecting the boundary information of the capillary groove, wherein: the solid region is set to 1 and the fluid region is set to 0.

[0022] Optionally, step S5 includes:

[0023] S5.1 Importing the 0-1 matrix reflecting the boundary information of the capillary groove into the numerical model in the form of a two-dimensional array;

[0024] S5.2 Assigning initial values of macroscopic quantities to the fluid region, and the macroscopic quantities include macroscopic density and macroscopic velocity, wherein: the initial density value is p 1 , and the initial velocity value is u 1 ;

[0025] S5.3 Calculating the equilibrium state function of the multi-relaxation model moment space in the j-th cycle, and substituting the initial density value and the initial velocity value into the equilibrium state function to calculate the initial velocity distribution function. The equilibrium state function m j eq The specific calculation formula is as follows:

[0026] m j eq =ρ j (1, -2 + 3u j 2 , 1 - 3u j 2 , u x , -u x , u y , -u y , u x 2 , -u y 2 , u x , u y ) T ;

[0027] Wherein: j is the number of cycles, taking a natural number greater than or equal to 1; ρ jis the macroscopic density at the j-th cycle; u j is the macroscopic velocity at the j-th cycle, u j i.e., u x and u y is the resultant velocity; u x is the horizontal component velocity of u j ; u y is the vertical component velocity of u j ;

[0028] S5.4. Calculate the interaction potential ψ j and calculate the pseudo force based on the interaction potential using the common Shan-Chen pseudo force calculation formula. The interaction potential ψ j The specific calculation formula is as follows:

[0029]

[0030] where: P EOS is the gas equation of state; c s is the lattice sound speed in the numerical model, G is the interaction strength between fluid nodes; c is the lattice width in the numerical model;

[0031] S5.5. Calculate the distribution function m j * in the moment space after collision at the j-th cycle based on the pseudo force and the equilibrium state function in the multi-relaxation model moment space. The specific calculation formula is:

[0032] m j * = m - Λ[m - m j eq + δ t (I - 0.5Λ)S;

[0033] where: m is the distribution function before collision in the moment space; Λ is the relaxation matrix of the model; δ t (I - 0.5Λ)S is the expression of the force term in the moment space; δ t is the time step;

[0034] S5.6. Convert the distribution function after collision in the moment space to the velocity space to obtain the distribution function f j * in the velocity space after collision at the j-th cycle. The specific calculation formula is:

[0035] f j * = M -1 ·m j * ;

[0036] where: M-1 is the inverse matrix of the transformation matrix M;

[0037] S5.7. Calculate the migrated velocity distribution function based on the distribution function of the velocity space after the collision. The specific calculation formula is as follows:

[0038] f i (r + e i δ t , t + δ t ) j = f * i (r, t) j ;

[0039] where: f i (r + e i δ t , t + δ t ) j is the migrated velocity distribution function in the direction of the i-th velocity component of a certain grid in the j-th cycle, r is the particle position vector, and t is the time variable;

[0040] S5.8. If ∑ i f * i (r, t) j - f * i (r, t) j-1 < Z, or j ≥ N, then jump out of the loop, and obtain the time t taken for the water to fill the entire capillary groove by analyzing the changes in the macroscopic quantities obtained in each cycle E ; Otherwise, let j = j + 1, and update the macroscopic density and macroscopic velocity respectively, and return to S5.3, where: Z and N are set thresholds;

[0041] The update formulas for the macroscopic density and macroscopic velocity are as follows:

[0042] ρ j = ∑ i f i (r + e i δ t , t + δ t ) j-1 ;

[0043]

[0044] where: e i is the velocity configuration of the D2Q9 model in the direction of the i-th velocity component of a certain grid.

[0045] Optionally, in the step S5.2, the initial density values include the initial liquid phase density value and the initial gas phase density value. The initial liquid phase density value is taken as 0.45455, and the initial gas phase density value is taken as 0.00060511. The initial velocity values include the initial liquid phase velocity value and the initial gas phase velocity value, and both the initial liquid phase velocity value and the initial gas phase velocity value are taken as 0.

[0046] In the step S5.7, boundary conditions need to be considered when calculating the migration step. Periodic boundaries are adopted on the left and right of the numerical model, non-equilibrium extrapolation boundaries are adopted at the gas-liquid inlet and outlet, and half-step bounce-back boundaries are adopted between the gas-liquid and solid phases.

[0047] In the step S5.8, Z is 10 -6 ; N is 300000.

[0048] Optionally, in the step S6, the capillary groove cross-section parameters selected in any two adjacent adjustments are not the same.

[0049] Water absorption efficiency η k The specific calculation formula is:

[0050] η k = 1 / t E .

[0051] Optionally, in the step S1, the contour line of the initial cross-sectional shape is a convex curve with a smooth transition.

[0052] Optionally, in the step S1, the initial value range of the water absorption slit width is 0.25 mm - 0.35 mm; the value range of the water absorption slit length is 0.7 mm - 1 mm; the value range of the diversion groove width is 0.8 mm - 1.2 mm, and the initial value of the diversion groove length is l 1 , 1 mm ≤ l 1 ≤ 1.3 mm; d 1 + l 1 = 2.

[0053] The present invention also provides a computer-readable storage medium, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the above-mentioned capillary groove cross-section optimization design method based on LBM is realized.

[0054] The present invention also provides an electronic device, including: at least one processor, at least one memory, and computer program instructions stored in the memory. When the computer program instructions are executed by the processor, the above-mentioned capillary groove cross-section optimization design method based on LBM is realized.

[0055] Applying the technical solution of the present invention, the rising process of water in the capillary groove is simulated through a numerical model constructed based on LBM to study the influence of the cross-sectional shape of the capillary groove on its water absorption efficiency, and by adjusting the parameters of the cross-section of the capillary groove, an optimized cross-section of the capillary groove with better drainage effect is obtained.

[0056] In addition to the purposes, features and advantages described above, the present invention has other purposes, features and advantages. The present invention will be further described in detail below with reference to the drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0057] The drawings forming a part of this application are used to provide a further understanding of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation to the present invention. In the drawings:

[0058] Figure 1 is a schematic structural diagram of the capillary groove in an embodiment of the present invention;

[0059] Figure 2 is a schematic diagram of the numerical model calculation domain and boundary conditions in an embodiment of the present invention;

[0060] Figure 3 is a schematic structural diagram of capillary grooves with different cross-sections in an embodiment of the present invention;

[0061] Figure 4 is a comparison diagram of the water absorption efficiency of capillary grooves with different cross-sections in an embodiment of the present invention;

[0062] Figure 5 is a curve graph showing the change of the reduced water absorption amount of capillary grooves with different cross-sections over time steps in an embodiment of the present invention.

[0063] Among them, 1 is the diversion groove, and 2 is the water absorption slit. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0064] The following embodiments of the present invention will be described in detail with reference to the drawings, but the present invention can be implemented in many different ways defined and covered by the claims.

[0065] In one embodiment, as Figure 1 shown, the capillary groove on the capillary drainage member is the core component for water absorption and drainage, which is composed of a water absorption slit and a diversion groove, and the two are interconnected.

[0066] A method for optimizing the cross-section of a capillary groove based on LBM includes the following steps:

[0067] S1: Obtain the reduced area A 0 of the existing capillary groove 0 and the water absorption efficiency η

[0068] Obtain the initial values of the capillary groove cross-section parameters and the initial cross-section shape, where: the capillary groove cross-section parameters include the width of the water absorption slit 2, the length of the water absorption slit 2, the width of the diversion groove 1, and the length of the diversion groove 1; preferably, the contour line of the initial cross-section shape is a convex curve with a smooth transition. The value range of the initial width of the water absorption slit 2 is 0.25 mm - 0.35 mm; the value range of the length of the water absorption slit 2 is 0.7 mm - 1 mm; the value range of the width of the diversion groove 1 is 0.8 mm - 1.2 mm, and the initial value of the length of the diversion groove 1 is l 1 , 1 mm ≤ l 1 ≤ 1.3 mm; d 1 + l 1 = 2.

[0069] S2: Calculate the reduced area A of the capillary groove according to the capillary groove cross-section parameters k , k is the number of iterations, which takes a natural number greater than or equal to 1;

[0070] S3: If A k < A k-1 , then let k = k + 1, and adjust any one of the width of the water absorption slit 2, the length of the water absorption slit 2, the width of the diversion groove 1, and the length of the diversion groove 1 in the capillary groove cross-section parameters using the first adjustment formula, and return to S2; if A k ≥ A k-1 , then enter S4;

[0071] The first adjustment formula includes: b k = b k-1 + 2n; d k = d k-1 - 2n; h k = h k-1 + 2n; l k = l k-1 + 2n; where: b k is the width of the water absorption slit 2 in the k-th iteration; d k is the length of the water absorption slit 2 in the k-th iteration, h k is the width of the diversion groove 1 in the k-th iteration; l k is the length of the diversion groove 1 in the k-th iteration; n is the grid length in the numerical model;

[0072] In S3, the specific calculation formula for the reduced area A of the capillary groove is as follows: k

[0073]

[0074] where: A is the cross-sectional area of the capillary groove, h k is the width of the diversion groove 1, s is the distance between the diversion grooves;

[0075] ​Among any four adjacent adjustments, the capillary groove cross-section parameters selected are all different.

[0076] S4: Draw the cross-sectional image of the capillary groove according to the capillary groove cross-section parameters in the k-th iteration, and use image recognition technology to binarize the cross-sectional image of the capillary groove to obtain a 0-1 matrix reflecting the boundary information of the capillary groove;

[0077] Using image recognition technology to binarize the cross-sectional image of the capillary groove to obtain a 0-1 matrix reflecting the boundary information of the capillary groove specifically includes:

[0078] Obtain the picture threshold using the graythresh function;

[0079] According to the picture threshold, use the im2b w function to convert the cross-sectional image into a binary image, and obtain a 0-1 matrix reflecting the boundary information of the capillary groove, where: the solid region is set to 1 and the fluid region is set to 0.

[0080] S5: Import the 0-1 matrix reflecting the boundary information of the capillary groove into the numerical model to simulate the rising process of water in the capillary groove, and obtain the time t taken for the water to fill the entire capillary groove E ;

[0081] S5 includes:

[0082] S5.1: Import the 0-1 matrix reflecting the boundary information of the capillary groove into the numerical model in the form of a two-dimensional array;

[0083] S5.2: Assign initial density values and initial velocity values to the fluid region;

[0084] Among them, the initial density values include the initial liquid phase density value and the initial gas phase density value. The initial liquid phase density value is taken as 0.45455, and the initial gas phase density value is taken as 0.00060511; the initial velocity values include the initial liquid phase velocity value and the initial gas phase velocity value, and both the initial liquid phase velocity value and the initial gas phase velocity value are taken as 0;

[0085] S5.3: Calculate the equilibrium state function of the multi-relaxation model moment space, and substitute the initial density values and initial velocity values into the equilibrium state function to calculate the initial velocity distribution function. The equilibrium state function m of the multi-relaxation model moment space eq The specific calculation formula is as follows:

[0086] m j eq =ρ j (1, -2 + 3u j 2 , 1 - 3u j 2 , u x , -u x , uy , -u y , u x 2 , -u y 2 , u xu , y) T ;

[0087] Among them: j is the number of cycles, taking natural numbers greater than or equal to 1; p j is the macroscopic density at the j-th cycle; u j is the macroscopic velocity at the j-th cycle, u j i.e., u x , u y 's combined velocity; u x is the horizontal component velocity of u j ; u y is the vertical component velocity of u j ;

[0088] S5.4. Calculate the interaction potential ψ j , and calculate the pseudo force based on the interaction potential using the common Shan-Chen pseudo force calculation formula. The interaction potential ψ j The specific calculation formula is as follows:

[0089]

[0090] Among them: P EOS is the gas state equation; c s is the lattice sound speed in the numerical model, G is the interaction strength between fluid nodes; c is the lattice width in the numerical model;

[0091] S5.5. Calculate the distribution function m j * after collision in the moment space at the j-th cycle based on the pseudo force and the equilibrium state function in the multi-relaxation model moment space. The specific calculation formula is:

[0092] m j * = m - Λ[m - m j eq + δ t (I - 0.5Λ)S;

[0093] Among them: m is the distribution function before collision in the moment space; Λ is the relaxation matrix of the model; δ t (I - 0.5Λ)S is the expression of the force term in the moment space; δ t is the time step; preferably, δ t takes 1.0.

[0094] S5.6. Convert the distribution function after collision in the moment space to the velocity space to obtain the distribution function \(f\) in the velocity space after collision. j * , and the specific calculation formula is:

[0095] \(f\) j * = \(M^{-1}\cdot m\) j * ;

[0096] where: \(M\) -1 is the inverse matrix of the transformation matrix \(M\);

[0097] S5.7. Calculate the migrated velocity distribution function based on the distribution function in the velocity space after collision, and calculate the macroscopic quantities in the current computational domain. Among them, the specific calculation formula for the migrated velocity distribution function is as follows:

[0098] \(f\) i (\(\vec{r}+e\) i \(\delta\) t , \(t + \delta\) t ) j = \(f\) * i (\(\vec{r}, t\)) j ;

[0099] where: \(f\) i (\(\vec{r}+e\) i \(\delta\) t , \(t + \delta\) t ) j is the migrated velocity distribution function in the \(i\)-th velocity component direction of a certain grid in the \(j\)-th cycle, \(\vec{r}\) is the particle position vector, and \(t\) is the time variable;

[0100] As Figure 2 shown, boundary conditions need to be considered when calculating the migration step. Periodic boundaries are adopted on the left and right of the numerical model, non-equilibrium extrapolation boundaries are adopted at the gas-liquid inlet and outlet, and half-step bounce boundaries are adopted between the gas-liquid and solid phases;

[0101] S5.8. If \(\sum\) i \(f\) * i (\(\vec{r}, t\)) j - \(f\) * i (\(\vec{r}, t\)) j-1 < \(Z\), or \(j\geq N\), then jump out of the loop, and obtain the time \(t\) it takes for water to fill the entire capillary groove by analyzing the changes in the macroscopic quantities obtained in each cycle E ; Otherwise, let \(j = j + 1\), and update the macroscopic density and macroscopic velocity respectively, and return to S5.4, where: \(Z\) and \(N\) are set thresholds. Preferably, \(Z\) is \(10^{-6}\); \(N\) is 300000.

[0102] The update formulas for macroscopic density and macroscopic velocity are specifically as follows:

[0103] ρ j = ∑ i f i (r + e i δ t , t + δ t ) j-1 ;

[0104]

[0105] Where: e i is the D2Q9 model velocity configuration in the direction of the i-th velocity component of a certain lattice.

[0106] S6: Calculate the water absorption efficiency η of the capillary groove according to t E , if η k < η k or η k-1 < Q, then let k = k + 1, and adjust the width of the water absorption slit 2 or the width of the diversion groove 1 using the second adjustment formula, and return to S2; if η k ≥ η k and η k-1 ≥ Q, then obtain the optimized cross-section of the capillary groove; where: Q is the set threshold; according to the drainage requirements of the actual application environment, different water absorption efficiency thresholds can be set to meet the drainage requirements of different environments.

[0107] The second adjustment formula includes: b k = b k-1 - n; h k = h k-1 - n. k

[0108] In S6, the cross-section parameters of the capillary grooves selected for any two adjacent adjustments are different;

[0109] The specific calculation formula for the water absorption efficiency η k is:

[0110] η k = 1 / t E .

[0111] Figure 3 In the embodiment of the present invention, the optimized cross-section of the capillary groove is as shown in Figure 3 (c), the width of the water absorption slit 2 is 0.25 mm, the length is 0.7 mm, and the diversion groove 1 is designed as a capsule shape. Figure 3 (a) is the existing capillary groove with a circular cross-section, (b) The capillary groove with the existing water droplet-shaped cross-section is provided as a comparative example for analyzing the water absorption effect of the new generation of cross-sections. In the embodiments of the present invention, the reduced drainage area of the optimized capillary groove cross-section is increased by 21% and 6% respectively compared with that of the circular and water droplet-shaped cross-sections, meeting the requirements of the optimized drainage area design.

[0112] As Figure 4 and Figure 5 shown, the water absorption efficiency of the optimized capillary groove cross-section in the embodiments of the present invention is increased by 6% and 22% respectively compared with that of the circular cross-section and the water droplet-shaped cross-section, indicating that the optimized capillary groove cross-section has improved both in terms of the reduced drainage area and the water absorption efficiency compared with the previous two generations of designs, and is a better design for the capillary groove cross-section of the water nail.

[0113] It should be noted that the device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. In addition, in the attached drawings of the device embodiments provided by the present invention, the connection relationship between the modules indicates that they have a communication connection, which can be specifically implemented as one or more communication buses or signal lines. Those of ordinary skill in the art can understand and implement it without creative efforts.

[0114] This embodiment further includes an electronic device, including: at least one processor, at least one memory, and computer program instructions stored in the memory. When the computer program instructions are executed by the processor, the above-mentioned method for analyzing the stability of the mining-induced slope is performed.

[0115] Exemplarily, the computer program can be divided into one or more modules / units. The one or more modules / units are stored in the memory and executed by the processor to complete the present invention. The one or more modules / units can be a series of computer program instruction segments capable of performing specific functions, and these instruction segments are used to describe the execution process of the computer program in the electronic device.

[0116] The electronic device can be a computing device such as a mobile phone, a desktop computer, a notebook, a handheld computer, and a cloud server. The electronic device may include, but is not limited to, a processor and a memory. For example, the electronic device may further include input / output devices, network access devices, a bus, etc.

[0117] The so-called processor may be a Central Processing Unit (CPU), or it may also be other general-purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field-Programmable Gate Arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc. The processor is the control center of the electronic device and connects all parts of the entire electronic device through various interfaces and circuits.

[0118] The memory can be used to store the computer program and / or module. The processor realizes the computer program by running or executing the computer program and / or module stored in the memory, and by calling the data stored in the memory. The memory mainly includes a program storage area and a data storage area. Among them, the program storage area can store an operating system, application programs required for at least one function (such as a sound playback function, an image playback function, etc.), etc.; the data storage area can store data created according to the use of the mobile phone (such as audio data, phone book, etc.), etc. In addition, the memory may include high-speed random access memory, and may also include non-volatile memory, such as a hard disk, memory, plug-in hard disk, Smart Media Card (SMC), Secure Digital (SD) card, Flash Card, at least one magnetic disk storage device, flash memory device, or other volatile solid-state storage devices.

[0119] Among them, if the modules / units integrated in the electronic device are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, to implement all or part of the processes in the above-described embodiment methods of the present invention, it can also be completed by a computer program instructing relevant hardware. The computer program can be stored in a readable storage medium. When the computer program is executed by a processor, the steps of the above-described various method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file, or some intermediate form, etc. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disc, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0120] The foregoing are only the preferred embodiments of the present invention and are not intended to limit the present invention. For those skilled in the art, the present invention can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims

1. A capillary groove cross-section optimization design method based on LBM, characterized in that: The steps include: S1: Obtain the equivalent area A0 of the existing capillary grooves and the water absorption efficiency η0 of the capillary grooves; construct a numerical model based on LBM; Obtaining initial values ​​and initial cross-sectional shapes of the capillary groove cross-sectional parameters, wherein the capillary groove cross-sectional parameters include water absorption slit width, water absorption slit length, guide groove width, and guide groove length; S2: Calculate the equivalent area A of the capillary groove based on the capillary groove cross-sectional parameters k , k is the number of iterations, which is a natural number greater than or equal to 1; S3: If A k <A k-1 , then let k = k + 1, and use the first adjustment formula to adjust any one of the water absorption slit width, water absorption slit length, guide groove width and guide groove length in the capillary groove cross-section parameters, and return to S2; if A k ≥A k-1 , then enter S4; The first adjustment formula includes: b k =b k-1 +2n;d k =d k-1 -2n;h k =h k-1 +2n;l k = l k-1 +2n; where: b k is the width of the water absorption gap in the kth iteration; d k is the length of the water absorption gap in the kth iteration, h k is the width of the guide groove in the kth iteration; l k is the length of the diversion channel in the kth iteration; n is the grid length in the numerical model; S4: drawing a cross-sectional image of the capillary groove according to the cross-sectional parameters of the capillary groove in the kth iteration, and binarizing the cross-sectional image of the capillary groove using image recognition technology to obtain a 0-1 matrix reflecting the boundary information of the capillary groove; S5: Import the 0-1 matrix reflecting the boundary information of the capillary groove into the numerical model to simulate the rising process of water in the capillary groove, and obtain the time t taken for the water to fill the entire capillary groove E ; S6: According to t E Calculate the water absorption efficiency η of the capillary groove k , if η k <η k-1 or η k <Q, then let k=k+1, and use the second adjustment formula to adjust the width of the water absorption slit or the width of the guide groove, and return to S2; if η k ≥η k-1 And η k ≥Q, the optimized capillary groove cross section is obtained; where: Q is the set threshold; The second adjustment formula includes: b k =b k-1 -n;h k =h k-1 -n.

2. The LBM-based capillary groove cross-section optimization design method according to claim 1, characterized in that: In S2, the equivalent area A of the capillary groove is k The specific calculation formula is as follows: Where: A is the cross-sectional area of ​​the capillary groove, h k is the width of the guide groove, s is the distance between the guide grooves; The cross-sectional parameters of the capillary grooves selected in any four adjacent adjustments are all different.

3. The LBM-based capillary groove cross-section optimization design method according to claim 2, characterized in that: In S4, using image recognition technology to binarize the cross-sectional image of the capillary groove to obtain a 0-1 matrix reflecting the boundary information of the capillary groove specifically includes: Use graythresh function to obtain the image threshold; The cross-sectional image is converted into a binary image using the im2bw function according to the image threshold, and a 0-1 matrix reflecting the boundary information of the capillary groove is obtained, in which: the solid area is set to 1 and the fluid area is set to 0.

4. The LBM-based capillary groove cross-section optimization design method according to claim 3, characterized in that: The S5 includes: S5.1, importing the 0-1 matrix reflecting the boundary information of the capillary groove into the numerical model in the form of a two-dimensional array; S5.

2. Assign initial values ​​of macro quantities to the fluid region. The macro quantities include macro density and macro velocity, where the initial value of density is ρ1 and the initial value of velocity is u1; S5.

3. Calculate the equilibrium state function of the multi-relaxation model moment space of the j-th cycle, and substitute the initial density and velocity values ​​into the equilibrium state function to calculate the initial velocity distribution function. The equilibrium state function m of the multi-relaxation model moment space is j eq The specific calculation formula is as follows: m j eq =ρ j (1,-2+3 and j 2 ,1-3u j 2 ,and x ,-and x ,and y ,-and y ,and x 2 -and y 2 ,and x and y ) T ; Where: j is the number of cycles, which is a natural number greater than or equal to 1; ρ j is the macroscopic density of the jth cycle; u j is the macro speed of the jth cycle, u j That is u x 、u y The total speed of x for u j The horizontal component of velocity; u y for u j The vertical component of velocity; S5.

4. Calculation of the interaction potential ψ j , and the pseudo force is calculated based on the interaction potential using the commonly used Shan-Chen pseudo force calculation formula, the interaction potential ψ j The specific calculation formula is as follows: Where: P EOS is the gas state equation; c s is the lattice sound speed in the numerical model, G is the interaction strength between fluid nodes; c is the grid width in the numerical model; S5.

5. Calculate the distribution function m after the collision in the moment space of the jth cycle based on the equilibrium state function of the moment space of the pseudo-force and multi-relaxation model j * , the specific calculation formula is: m j * =m-Λ[m-m j eq ]+δ t (I-0.5Λ)S; Where: m is the distribution function before collision in moment space; Λ is the relaxation matrix of the model; δ t (I-0.5Λ)S is the expression of the force term in moment space; δ t is the time step; S5.

6. Convert the distribution function after collision in moment space to velocity space to obtain the distribution function f in velocity space after collision in the jth cycle j * , the specific calculation formula is: f j * =M ﹣1 ·m j * ; Where: M ﹣1 is the inverse matrix of the transformation matrix M; S5.

7. The velocity distribution function after migration is calculated based on the distribution function of the velocity space after the collision. The specific calculation formula is as follows: f i (r+e i δ t ,t+δ t ) j =f * i (r,t) j ; Where: f i (r+e i δ t ,t+δ t ) j is the velocity distribution function after migration in the direction of the i-th velocity component of a certain grid in the j-th cycle, r is the particle position vector, and t is the time variable; S5.8, if ∑ i f * i (r,t) j -f * i (r,t) j-1 <Z, or j≥N, then the cycle is jumped out, and the time t taken for water to fill the entire capillary groove is obtained by analyzing the changes in the macroscopic quantity obtained in each cycle. E ; Otherwise, let j = j + 1, and update the macro density and macro speed respectively, and return to S5.3, where: Z and N are set thresholds; The update formulas for macro density and macro speed are as follows: r j =∑ i f i (r+e i d t ,t+δ t ) j-1 ; Where: e i It is the velocity configuration of the D2Q9 model in the direction of the i-th velocity component of a certain grid.

5. The LBM-based capillary groove cross-section optimization design method according to claim 4, characterized in that: In S5.2, the initial density value includes the initial value of liquid phase density and the initial value of gas phase density, the initial value of liquid phase density is 0.45455, and the initial value of gas phase density is 0.00060511; the initial velocity value includes the initial value of liquid phase velocity and the initial value of gas phase velocity, and the initial value of liquid phase velocity and the initial value of gas phase velocity are both 0; In S5.7, boundary conditions need to be considered when calculating the migration step. The left and right sides of the numerical model use periodic boundaries, the gas-liquid inlet and outlet use non-equilibrium extrapolation boundaries, and the gas-liquid and solid phases use half-step rebound boundaries. In S5.8, Z is 10 -6 ; N is 300000.

6. The LBM-based capillary groove cross-section optimization design method according to any one of claims 1 to 5, characterized in that: In S6, the cross-sectional parameters of the capillary grooves selected in any two adjacent adjustments are different; Water absorption efficiency η k The specific calculation formula is: or k =1 / t E 。 7. The LBM-based capillary groove cross-section optimization design method according to any one of claims 1 to 5, characterized in that: In S1, the contour line of the initial cross-sectional shape is a convex curve with a smooth transition.

8. The LBM-based capillary groove cross-section optimization design method according to any one of claims 1 to 5, characterized in that: In S1, the initial value of the water absorption slit width ranges from 0.25mm to 0.35mm; the water absorption slit length ranges from 0.7mm to 1mm; the guide groove width ranges from 0.8mm to 1.2mm, and the initial value of the guide groove length is l1, 1mm≤l1≤1.3mm; d1+l1=2.

9. A computer-readable storage medium, characterized in that: Computer program instructions are stored thereon, and when the computer program instructions are executed by a processor, the LBM-based capillary groove cross-section optimization design method as described in any one of claims 1 to 8 is implemented.

10. An electronic device, characterized in that: include: At least one processor, at least one memory and computer program instructions stored in the memory, when the computer program instructions are executed by the processor, the LBM-based capillary groove cross-section optimization design method as described in any one of claims 1 to 8.

Images