A method, device, electronic device and storage medium for determining contact angle

Abstract

The present application discloses a contact angle determination method, device, electronic device and storage medium, the method comprising: constructing a flow field; the flow field comprises multiple components and a flow field grid constructed based on the flow field, the multiple components comprising a first fluid component, a second fluid component and a solid component; initializing a flow field model based on the flow field; determining a first parameter and a second parameter; determining a third parameter and a fourth parameter according to the first parameter and the second parameter; determining a first force of each grid point in the flow field grid according to the flow field model, the first parameter and the second parameter, and determining a second force of each grid point according to the flow field model, the third parameter and the fourth parameter; determining a current contact angle according to the first force and the second force; when the flow field model is at the current contact angle, the first fluid component is stationary on the fixed component, and the current contact angle is determined to be the target contact angle. In this way, the contact angle calculation process can be simplified, and the contact angle simulation calculation result is stable and accurate.

Description

Technical Field

[0001] The present application relates to the field of computational fluid dynamics, and in particular to a contact angle determination method, device, electronic device and storage medium. Background Art

[0002] Surface wettability is an important characteristic in multiphase flow and is used to characterize the ductility of liquid on the solid surface. Surface wettability is measured by the size of the contact angle. The contact angle is the characteristic angle (angle) formed by the contact position between gas, liquid and solid, specifically the angle formed by the tangent of the interface formed by the gas phase and the liquid phase and the tangent of the interface formed by the solid phase and the liquid phase. Measuring the contact angle can obtain a lot of information on the interaction between the gas-liquid, solid-gas, and solid-liquid interfaces on the material surface, determine the surface wettability of the material, and has extremely high application value in industrial fields such as microfluidic chips, medical health, oil extraction, and new materials.

[0003] At present, the methods for measuring contact angle are mainly based on experiments or numerical simulations of droplet wetting phenomena, among which numerical simulations are widely welcomed due to their advantages of speed, efficiency and low cost. However, for multi-component models including contact angle, since multi-component models involve more interactions between different components and different phases, they are more complicated than single-component models. It is often difficult to select appropriate parameters for the strength of fluid-solid interaction in simulations, thus affecting the stability of the system and the accuracy of the simulation results. Summary of the invention

[0004] The embodiments of the present application provide a contact angle determination method, device, electronic device and storage medium, which can simplify the contact angle calculation process, and the contact angle simulation calculation results are stable and accurate.

[0005] On the one hand, an embodiment of the present application provides a contact angle method, the method comprising:

[0006] Constructing a flow field; the flow field includes a plurality of components and a flow field grid constructed based on the flow field, the plurality of components including a first fluid component, a second fluid component and a solid component;

[0007] Initialize the flow field model based on the flow field;

[0008] Determine a first parameter and a second parameter; the first parameter is the strength of the first fluid component on the second fluid component; the second parameter is the strength of the second fluid component on the first fluid component;

[0009] Determine a third parameter and a fourth parameter according to the first parameter and the second parameter; the third parameter is the strength of the solid component on the first fluid component; the fourth parameter is the strength of the solid component on the second fluid component;

[0010] The first force of each grid point in the flow field grid is determined according to the flow field model, the first parameter and the second parameter, and the second force of each grid point is determined according to the flow field model, the third parameter and the fourth parameter; the first force is the force of the second fluid component and the third fluid component on the first fluid component or the force of the first fluid component and the fourth fluid component on the second fluid component, and the second force is the force of the solid component on the first fluid component or the second fluid component; the third fluid component and the first fluid component are the same fluid component, and the fourth fluid component and the second fluid component are the same fluid component;

[0011] Determine a current contact angle according to the first force and the second force; the current contact angle is the contact angle between the first fluid component and the solid component;

[0012] When the flow field model is located at the current contact angle, the first fluid component is stationary on the fixed component, and the current contact angle is determined to be the target contact angle.

[0013] Furthermore, the flow field model is initialized based on the flow field, including:

[0014] Setting initial parameters of the flow field grid; the initial parameters include any one or more of the solution range of the flow field grid, the component type of each grid point, the initial density of the first fluid component, the initial velocity of the first fluid component, the initial density of the second fluid component and the initial velocity of the second fluid component;

[0015] determining a distribution function of the first fluid component according to an initial density of the first fluid component and an initial velocity of the first fluid component, and determining a distribution function of the second fluid component according to an initial density of the second fluid component and an initial velocity of the second fluid component;

[0016] Determine the migration step and collision step of each grid point according to the distribution function;

[0017] A distribution function of the boundary of the first fluid component is determined based on the migration step and the collision step of each grid point.

[0018] Furthermore, after determining the current contact angle according to the first force and the second force, the method further includes:

[0019] If the first fluid component moves on the fixed component when the flow field model is at the current contact angle, the grid point density and the grid point velocity of each grid point are determined according to the flow field model, the first force and the second force;

[0020] Update the flow field model according to the grid point density and grid point velocity of each grid point; repeat the steps of: determining the first force of each grid point of the flow field according to the updated flow field model, the first parameter and the second parameter, and determining the second force of each grid point of the flow field according to the updated flow field model, the third parameter and the fourth parameter; determine the current contact angle according to the first force and the second force;

[0021] Until the first fluid component is stationary on the fixed component when the flow field model is located at the current contact angle, the current contact angle is determined to be the target contact angle.

[0022] Furthermore, the flow field model is updated according to the grid density and grid velocity of each grid point, including:

[0023] Update the distribution function of the first fluid component and the distribution function of the second fluid component according to the grid point density and the grid point velocity of each grid point;

[0024] Update the migration step and collision step of each grid point according to the updated distribution function of the first fluid component and the distribution function of the second fluid component;

[0025] updating a distribution function of a boundary of a first fluid component according to a migration step and a collision step at each grid point;

[0026] The flow field model is updated according to the updated distribution function of the first fluid component, the distribution function of the second fluid component, and the distribution function of the boundary of the first fluid component.

[0027] Further, the distribution function of the boundary of the first fluid component is determined according to the migration step and the collision step of each grid point, including:

[0028] Setting boundary conditions for the first fluid component;

[0029] Based on the boundary conditions, a distribution function of the boundary of the first fluid component is determined according to the migration step and the collision step of each grid point.

[0030] Further, setting the boundary conditions of the first fluid component includes:

[0031] Set the left and right boundaries of the first fluid component as periodic boundaries;

[0032] Set the lower boundary of the first fluid component to a non-slip boundary using the standard rebound format;

[0033] The upper boundary of the first fluid component is set using a non-equilibrium extrapolation method.

[0034] Furthermore, before determining the third parameter and the fourth parameter according to the first parameter and the second parameter, the method further includes:

[0035] Determine a fifth parameter and a sixth parameter; the fifth parameter is the strength of the action of the third fluid component on the first fluid component; the sixth parameter is the strength of the action of the fourth fluid component on the second fluid component;

[0036] determining a first pseudopotential function of the first fluid component based on the fifth parameter and the equation of state of the first fluid component;

[0037] determining a second pseudopotential function of the second fluid component based on the sixth parameter and the equation of state of the second fluid component;

[0038] The first pseudo potential function, the second pseudo potential function, the fifth parameter and the sixth parameter are used to determine the numerical relationship between the first parameter, the second parameter and the third parameter and the fourth parameter.

[0039] On the other hand, an embodiment of the present application provides a contact angle determination device, the device comprising:

[0040] An initialization module is used to construct a flow field; a flow field model is initialized based on the flow field;

[0041] A parameter determination module, used to determine a first parameter and a second parameter; the first parameter is the strength of the first fluid component on the second fluid component; the second parameter is the strength of the second fluid component on the first fluid component; and is also used to determine a third parameter and a fourth parameter based on the first parameter and the second parameter; the third parameter is the strength of the solid component on the first fluid component; the fourth parameter is the strength of the solid component on the second fluid component;

[0042] A force determination module is used to determine the first force of each grid point in the flow field grid according to the flow field model, the first parameter and the second parameter, and to determine the second force of each grid point according to the flow field model, the third parameter and the fourth parameter; the first force is the force of the second fluid component and the third fluid component on the first fluid component or the force of the first fluid component and the fourth fluid component on the second fluid component, and the second force is the force of the solid component on the first fluid component or the second fluid component; the third fluid component and the first fluid component are the same fluid component, and the fourth fluid component and the second fluid component are the same fluid component;

[0043] A current contact angle determination module, used to determine a current contact angle according to the first force and the second force; the current contact angle is a contact angle between the first fluid component and the solid component;

[0044] The target contact angle determination module is used to determine the current contact angle as the target contact angle when the flow field model is located at the current contact angle and the first fluid component is stationary on the fixed component.

[0045] In an optional embodiment, the initialization module 501 is used to set the initial parameters of the flow field grid; the initial parameters include the solution range of the flow field grid, the component type of each grid point, the initial density of the first fluid component, the initial velocity of the first fluid component, the initial density of the second fluid component and the initial velocity of the second fluid component; the distribution function of the first fluid component is determined according to the initial density of the first fluid component and the initial velocity of the first fluid component, and the distribution function of the second fluid component is determined according to the initial density of the second fluid component and the initial velocity of the second fluid component; the migration step and collision step of each grid point are determined according to the distribution function; the distribution function of the boundary of the first fluid component is determined according to the migration step and collision step of each grid point.

[0046] Furthermore, the device also includes:

[0047] A model updating module is used to determine the grid density and grid speed of each grid point according to the flow field model, the first force and the second force when the first fluid component moves on the fixed component when the flow field model is at the current contact angle; update the flow field model according to the grid density and grid speed of each grid point; repeat the steps of: determining the first force of each grid point of the flow field according to the updated flow field model, the first parameter and the second parameter, determining the second force of each grid point of the flow field according to the updated flow field model, the third parameter and the fourth parameter; and determining the current contact angle according to the first force and the second force.

[0048] Furthermore, the model updating module is also used to update the distribution function of the first fluid component and the distribution function of the second fluid component according to the grid density and grid velocity of each grid point; update the migration step and collision step of each grid point according to the updated distribution function of the first fluid component and the distribution function of the second fluid component; update the distribution function of the boundary of the first fluid component according to the migration step and collision step of each grid point; update the flow field model according to the updated distribution function of the first fluid component, the distribution function of the second fluid component and the distribution function of the boundary of the first fluid component.

[0049] Furthermore, the initialization module is also used to set the boundary conditions of the first fluid component; based on the boundary conditions, the distribution function of the boundary of the first fluid component is determined according to the migration step and collision step of each grid point.

[0050] Furthermore, the initialization module is used to set the left and right boundaries of the first fluid component as periodic boundaries; set the lower boundary of the first fluid component as a non-slip boundary using a standard rebound format; and set the upper boundary of the first fluid component using a non-equilibrium extrapolation format.

[0051] Furthermore, a parameter determination module is used to determine a fifth parameter and a sixth parameter; the fifth parameter is the intensity of the action of the third fluid component on the first fluid component; the sixth parameter is the intensity of the action of the fourth fluid component on the second fluid component; the first pseudopotential function of the first fluid component is determined based on the fifth parameter and the state equation of the first fluid component; the second pseudopotential function of the second fluid component is determined based on the sixth parameter and the state equation of the second fluid component; the first pseudopotential function, the second pseudopotential function, the fifth parameter and the sixth parameter are used to determine the numerical relationship between the first parameter, the second parameter and the third parameter, the fourth parameter.

[0052] On the other hand, an embodiment of the present application provides an electronic device, which includes a processor and a memory, wherein the memory stores at least one instruction or at least one program, and the at least one instruction or at least one program is loaded by the processor and executes any of the contact angle determination methods described above.

[0053] On the other hand, an embodiment of the present application provides a computer-readable storage medium, in which at least one instruction or at least one program is stored, and the at least one instruction or at least one program is loaded and executed by a processor to implement any of the contact angle determination methods described above.

[0054] The contact angle determination and allocation method, device, electronic device, and storage medium provided in the embodiments of the present application have the following technical effects:

[0055] Construct a flow field; the flow field includes multiple components and a flow field grid constructed based on the flow field, and the multiple components include a first fluid component, a second fluid component and a solid component; initialize a flow field model based on the flow field; determine a first parameter and a second parameter; determine a third parameter and a fourth parameter based on the first parameter and the second parameter; determine a first force for each grid point in the flow field grid based on the flow field model, the first parameter and the second parameter, and determine a second force for each grid point based on the flow field model, the third parameter and the fourth parameter; determine a current contact angle based on the first force and the second force; when the flow field model is located at the current contact angle, the first fluid component is stationary on the fixed component, and the current contact angle is determined to be the target contact angle. In this way, the third parameter and the fourth parameter are determined by the first parameter and the second parameter, which simplifies the determination process of the third parameter and the fourth parameter, and the first force and the second force can be directly obtained based on the third parameter and the fourth parameter, and the current contact angle can be further determined. The current contact angle can be determined through three steps, which simplifies the calculation process of the contact angle; and the appropriate third parameter and the fourth parameter can be determined based on the first parameter and the second parameter, which improves the stability of the simulation, thereby improving the accuracy of the contact angle simulation result. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present application or the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0057] Figure 1 It is a schematic diagram of an application scenario provided by an embodiment of the present application;

[0058] Figure 2 It is a schematic flow chart of a contact angle determination method provided in an embodiment of the present application;

[0059] Figure 3 is a flow field schematic diagram of a contact angle determination method provided in an embodiment of the present application;

[0060] Figure 4 It is a schematic diagram of a flow chart of a contact angle determination method provided in an embodiment of the present application;

[0061] Figure 5 is a structural schematic diagram of a contact angle determination device provided in an embodiment of the present application;

[0062] Figure 6 This is a hardware structure block diagram of a server that provides a contact angle determination method according to an embodiment of the present application. DETAILED DESCRIPTION

[0063] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions 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 making creative work are within the scope of protection of this application.

[0064] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or server that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0065] See also Figure 1 , Figure 1 1 is a schematic diagram of an application scenario provided by an embodiment of the present application, and the schematic diagram includes a first fluid component 101, a second fluid component 102, and a solid component 103. In this application scenario, a flow field model is established, and the first fluid component 101 and the second fluid component 102 can be miscible fluid components or immiscible fluid components. The first fluid component 101 and the second fluid component 102 are arranged on the solid component 103, and the contact angle between the first fluid component 101 and the solid component 103 is simulated and determined by the flow field model, thereby determining the surface wettability of the solid component 103.

[0066] Specifically, the contact angle determination method provided in the embodiment of the present application constructs a flow field and initializes a flow field model based on the flow field; wherein the flow field includes a first fluid component 101, a second fluid component 102 and a solid component 103, and a flow field grid constructed based on the flow field; determines a first parameter and a second parameter; wherein the first parameter is the intensity of the action of the first fluid component 101 on the second fluid component 102, and the second parameter is the intensity of the action of the second fluid component 102 on the first fluid component 101; determines a third parameter and a fourth parameter based on the first parameter and the second parameter; wherein the third parameter is the intensity of the action of the solid component 103 on the first fluid component 101, and the fourth parameter is the intensity of the action of the solid component 103 on the second fluid component 102; determines each parameter in the flow field grid based on the flow field model, the first parameter and the second parameter. The first force of the grid point, and the second force of each grid point are determined according to the flow field model, the third parameter and the fourth parameter; wherein the first force is the force of the second fluid component and the third fluid component on the first fluid component or the force of the first fluid component and the fourth fluid component on the second fluid component, and the second force is the force of the solid component on the first fluid component or the second fluid component; the third fluid component and the first fluid component are the same fluid components, and the fourth fluid component and the second fluid component are the same fluid components; the current contact angle is determined according to the first force and the second force, and the current contact angle is the contact angle between the first fluid component 101 and the solid component 103; when the flow field model is at the current contact angle, the first fluid component 101 is stationary on the solid component 103, and the current contact angle is determined to be the target contact angle.

[0067] Among them, in the embodiments of the present application, Figure 1The first fluid component 101 shown can be water. In this case, the flow field model can be used to determine whether the solid component 103 is hydrophilic or hydrophobic. When the contact angle is less than 90°, it means that the solid component 103 is hydrophilic. When the contact angle is greater than 90°, it means that the solid component 103 is hydrophobic. The first fluid component 101 can also be other gases or liquids. The flow field model is used to determine the surface wetting degree of different gases or liquids on the solid component 103, and the surface characteristics of the solid component 103 under different gas or liquid environments are determined to assist production researchers in determining the gas or liquid applied to the surface of the solid component 103 or the surrounding environment. The first fluid component 101 can also be a fixed gas or liquid. The solid component 103 can be different types of materials. The flow field model is used to determine the surface wetting degree of the first fluid component 101 on different types of materials, and assist production researchers in determining the application of the first fluid component 101 on different materials.

[0068] The following describes a specific embodiment of a contact angle determination method of the present application. Figure 2 It is a flow chart of a contact angle determination method provided in an embodiment of the present application. This specification provides method operation steps such as the embodiment or flow chart, but may include more or fewer operation steps based on conventional or non-creative labor. The order of steps listed in the embodiment is only one way of executing the steps among many orders, and does not represent the only order of execution. When the actual system or server product is executed, it can be executed in the order shown in the embodiment or the drawings or in parallel (for example, in a parallel processor or multi-threaded processing environment). Specifically, Figure 2 As shown, the method may include:

[0069] S201: Construct a flow field.

[0070] In the embodiment of the present application, the flow field includes multiple components and a flow field grid constructed based on the flow field, and the multiple components include a first fluid component, a second fluid component and a solid component. Determine that the first fluid component is σ and the second fluid component is The solid component is ads.

[0071] Optionally, the flow field may be two-dimensional, and the corresponding flow field grid is a two-dimensional plane. Specifically, Figure 3 A flow field schematic diagram of a contact angle determination method provided in an embodiment of the present application is shown, such as Figure 3As shown in the figure, the flow field grid is a hypothetical plane composed of multiple square grid points, that is, the grid step lengths in the x-direction and the y-direction are equal, and multiple components are distributed on the grid points of the plane. In a two-dimensional flow field, the flow field grid represents the planar effect of the droplet wetting phenomenon, which is equivalent to observing the droplet wetting phenomenon from the side. However, unlike the macroscopic droplet wetting phenomenon, the flow field grid represents the droplet wetting phenomenon from the perspective of microscopic particles. Through the movement of particles on the grid points in the flow field grid, a flow field model is established to simulate the droplet wetting phenomenon and determine the size of the contact angle.

[0072] Optionally, the flow field can be three-dimensional, and the corresponding flow field grid is a three-dimensional space. In this case, the flow field grid is a space assumed to be composed of multiple cubic grid points, that is, the grid steps in the x-direction, y-direction and z-direction are equal. In the three-dimensional flow field, the flow field grid characterizes the spatial effect of the droplet wetting phenomenon.

[0073] S203: Initialize the flow field model based on the flow field.

[0074] In the embodiment of the present application, by setting the initial parameters of the flow field grid, including the initial density and initial velocity of the first fluid component σ, the second fluid component The initial density and initial velocity of the first fluid component σ are determined; the distribution function of the first fluid component σ is determined according to the initial density and initial velocity of the first fluid component σ, and the distribution function of the second fluid component σ is determined according to the initial density and initial velocity of the first fluid component σ. The initial density and initial velocity of the second fluid component determine distribution function; determine the migration step and collision step of each grid point according to the above distribution function; determine the distribution function of the boundary of the first fluid component σ according to the migration step and collision step of each grid point, thereby completing the initialization of the flow field model.

[0075] In the embodiment of the present application, the first fluid component σ and the second fluid component All satisfy the lattice Boltzmann evolution equation (LBM evolution equation), which can be expressed as formula (1):

[0076] f α (x+e α δt,t+δt)=f α (x, t)-Ω α +δtF α (x, t)……(1)

[0077] Among them, f α (x, t) is the first fluid component σ or the second fluid component The distribution function of , δt is the time step, Ω α is the collision term of the multi-relaxation time model of the first fluid component σ, F α (x, t) is the external force term in velocity space, e αis the discrete velocity vector.

[0078] Optionally, for two-dimensional flow, the discrete velocity vector e α It satisfies the D2Q9 lattice model, that is, there are 9 discrete velocity directions in the two-dimensional plane, which can be specifically expressed as formula (2):

[0079]

[0080] Among them, c = δx / δt, c is the grid constant, and δx is the grid step size.

[0081] Optionally, for three-dimensional flow, the discrete velocity vector e α It satisfies the D3Q19 lattice model, that is, there are 19 discrete velocity directions in three-dimensional space.

[0082] In the embodiment of the present application, the initial parameters of the flow field grid include the solution range of the flow field grid, the component type of each grid point, the initial density of the first fluid component σ, the initial velocity of the first fluid component σ, the second fluid component The initial density and the second fluid composition The initial velocity.

[0083] The solution range of the flow field grid is the size of the flow field grid. For example, setting the solution range to Lx=200, Ly=200 indicates that the flow field grid is 200 grid units long in the x direction and 200 grid units long in the y direction.

[0084] Optionally, the component types of the grid points include the first fluid component node, the second fluid component node, the solid component node, the first fluid component and the second fluid component boundary point, the first fluid component and the solid component boundary point, and the second fluid component and the solid component boundary point. Among them, the first fluid component node is a grid point in the first fluid component σ that does not contact other components, and the second fluid component node is a grid point in the second fluid component. The solid component node is the grid point that does not contact other components in the solid component ads. The boundary point between the first fluid component σ and the second fluid component The contact lattice point of the first fluid component and the solid component is the lattice point where the first fluid component σ contacts the solid component ads, and the contact lattice point of the second fluid component and the solid component is the lattice point where the second fluid component The lattice points in contact with the solid component ads. Figure 3 As shown, the grid points of the component types, namely, the intersection points of the first fluid component and the second fluid component and the intersection points of the first fluid component and the solid component, constitute the boundary of the first fluid component, representing the droplets in the macroscopic phenomenon. The more grid points the first fluid component σ contains, the larger the flow field grid range covered, and the more accurate the simulation of the boundary of the first fluid component σ.

[0085] In the embodiment of the present application, the initial density of the first fluid component σ is determined as Determine the initial velocity of the first fluid component σ as The initial density of the first fluid component σ Initial density and initial velocity Discretize to each discrete velocity direction and obtain the distribution function of the first fluid component σ Determine the second fluid composition The initial density is Determine the second fluid composition The initial velocity is The second fluid component The initial density and initial velocity Discrete into each discrete velocity direction to obtain the second fluid component The distribution function of

[0086] Optionally, the initial density of the first fluid component σ and a second fluid component The initial density It can be set according to the actual simulation requirements. The initial velocity of the first fluid component σ and a second fluid component Initial velocity Can be set to zero.

[0087] Furthermore, particles of each component in the flow field model migrate from one grid point to another, changing the distribution function of the particles, generating a migration step, and collide under the premise of conservation of mass, momentum and energy, changing the distribution function of the particles, generating a collision step. Specifically, determining the migration step and collision step of each grid point according to the distribution function includes: Determine the migration and collision steps of the grid points of the first fluid component nodes, and the second fluid component The distribution function of Determine the migration and collision steps for the grid points of the nodes of the second fluid component.

[0088] According to the distribution function of the first fluid component σ The migration step of the grid point for determining the first fluid component node can be expressed as formula (3):

[0089]

[0090] in, is the distribution function of the particles of the first fluid component σ migrating to the grid point at the previous time step, is the distribution function of the particles of the first fluid component σ at the current grid point.

[0091] According to the distribution function of the first fluid component σ The collision step of the grid point of the first fluid component node can be expressed as formula (4):

[0092]

[0093] According to the second fluid component The distribution function of The migration step of the grid points for determining the second fluid component nodes can be expressed as formula (5):

[0094]

[0095] in, The second fluid component The distribution function of the particle migrated to the grid point at the previous time step, The second fluid component The distribution function of the particles at the current grid point.

[0096] According to the second fluid component The distribution function of The collision step for determining the grid points of the second fluid component nodes can be expressed as formula (6):

[0097]

[0098] Furthermore, a boundary condition of the first fluid component σ is set, and based on the boundary condition, a distribution function of the boundary of the first fluid component σ is determined according to the migration step and the collision step of each grid point.

[0099] In the embodiment of the present application, the left and right boundaries of the first fluid component σ are set as periodic boundaries, the lower boundary is a non-slip boundary, a standard rebound format is adopted, and the upper boundary adopts a non-equilibrium extrapolation format. Optionally, the lower boundary can also adopt a half-step rebound format.

[0100] S205: Determine the first parameter G 12 and the second parameter G 21 .

[0101] In the embodiment of the present application, the first parameter G 12 is the first fluid component σ to the second fluid component The intensity of the action, the second parameter G 21 The second fluid component The intensity of the action on the first fluid component σ.

[0102] In the embodiment of the present application, the first parameter G 12 and the second parameter G21 It can be any value between 0.001 and 0.1. The first parameter G 12 With the second parameter G 21 When the first parameter G 12 With the second parameter G 21 When it is a negative number, it means the first fluid component σ and the second fluid component are miscible fluid components, when the first parameter G 12 With the second parameter G 21 When it is a positive number, it means the first fluid component σ and the second fluid component Are immiscible fluid components.

[0103] S207: According to the first parameter G 12 and the second parameter G 21 Determine the third parameter G ads，1 and the fourth parameter G ads，2 .

[0104] In the embodiment of the present application, the third parameter G ads，1 is the strength of the solid component ads on the first fluid component σ, and the fourth parameter G ads，2 The solid component ads is the second fluid component intensity of action.

[0105] In the embodiment of the present application, according to the first parameter G 12 and the second parameter G 21 Determine the third parameter G ads，1 and the fourth parameter G ads，2 Before, you can also include:

[0106] Determine the fifth parameter G 11 and the sixth parameter G 22 , according to the fifth parameter G 11 and the state equation of the first fluid component σ Determine the first pseudopotential function ψ of the first fluid component σ σ (x), according to the sixth parameter G 22 and a second fluid component The state equation Determine the second fluid composition The second pseudopotential function The first pseudopotential function ψ σ (x), the second pseudopotential function The fifth parameter G 11 and the sixth parameter G 22 To determine the first parameter G 12 , the second parameter G 21 and the third parameter G ads，1 , the fourth parameter G ads，2The numerical relationship between .

[0107] The third fluid component and the first fluid component σ are the same fluid component, and the fourth fluid component and the second fluid component The fifth parameter G 11 is the strength of the third fluid component on the first fluid component σ, and the fifth parameter G 11 The sixth parameter G represents the intensity of the effect of particles in the first fluid component σ on other particles in the first fluid component σ. 22 The fourth fluid component is the second fluid component The sixth parameter G 22 Represents the second fluid component Internal particles to the second fluid component The intensity of the interaction with other particles inside.

[0108] Optionally, according to the fifth parameter G 11 and the state equation of the first fluid component σ Determine the first pseudopotential function ψ of the first fluid component σ σ (x) can be expressed as formula (7):

[0109]

[0110] According to the sixth parameter G 22 and a second fluid component The state equation Determine the second fluid composition The second pseudopotential function It can be expressed as formula (8):

[0111]

[0112] in, is the lattice sound velocity, the state equation of the first fluid component σ and a second fluid component The state equation There are many options. In the embodiment of the present application, the ideal gas state equation is selected to describe the ideal fluid, and the Carnahan-Starling (CS) state equation is selected to describe the non-ideal fluid.

[0113] In the embodiment of the present application, the first parameter G 12 , the second parameter G 21 , the third parameter G ads，1 and the fourth parameter G ads，2 The numerical relationship between them can be expressed as formula (9). Based on formula (9), according to the first parameter G 12 , the second parameter G 21 Determine the third parameter Gads，1 and the fourth parameter G ads，2 :

[0114]

[0115] Among them, a is an indefinite constant.

[0116] The first parameter G provided in the embodiment of the present application 12 , the second parameter G 21 , the third parameter G ads，1 and the fourth parameter G ads，2 The numerical relationship between them can simplify the third parameter G ads，1 and the fourth parameter G ads，2 The calculation of the contact angle is convenient, and the simulation results of the flow field model are more stable.

[0117] S209: According to the flow field model, the first parameter G 12 and the second parameter G 21 Determine the first force F1(x, t) at each grid point in the flow field grid, according to the flow field model and the third parameter G ads，1 and the fourth parameter G ads，2 Determine the second force F2(x, t) at each grid point.

[0118] Among them, the first force F1(x, t) is the second fluid component and the force of the third fluid component on the first fluid component σ Or the first fluid component σ and the fourth fluid component to the second fluid component The force The first force F1(x, t) represents the interaction force between the fluid components. The second force F2(x, t) is the force of the solid component ads on the first fluid component σ or solid component ads to the second fluid component The force The second force F2(x, t) represents the force of the solid component on the fluid component.

[0119] In the embodiment of the present application, according to the flow field model, the first parameter G 12 and the second parameter G 21 Determine the second fluid composition and the force of the third fluid component on the first fluid component σ It can be expressed as formula (10):

[0120] According to the flow field model, the first parameter G 12 and the second parameter G 21Determine the effect of the first fluid component σ and the fourth fluid component on the second fluid component The force It can be expressed as formula (11):

[0121]

[0122] Among them, w(|e α | 2 ) are the weights of particles in different discrete speed directions, w(1) = 1 / 3, w(2) = 1 / 12.

[0123] According to the flow field model, the third parameter G ads，1 and the fourth parameter G ads，2 Determine the force exerted by the solid component ads on the first fluid component σ It can be expressed as formula (12):

[0124]

[0125] According to the flow field model, the third parameter G ads，1 and the fourth parameter G ads，2 Determine the solid component ads to the second fluid component The force It can be expressed as formula (13):

[0126]

[0127] Among them, s(x+e α δt) is the index function. For the solid component node, the boundary between the first fluid component and the solid component, and the boundary between the second fluid component and the solid component, the index function s(x+e α δt) is 1, for the first fluid component node and the second fluid component node, the index function s(x+e α δt) is 0.

[0128] S211: Determine the current contact angle according to the first force F1 (x, t) and the second force F2 (x, t)

[0129] In the present application embodiment, the current contact angle is the contact angle between the first fluid component σ and the solid component ads.

[0130] Specifically, on the surface of the solid component completely wetted by the fluid component, in the y direction, the adhesion force between the solid component and the fluid component is equal to the cohesive force of the fluid component, that is, in the y direction, the first force F1(x, t) and the second force F2(x, t) are equal in magnitude and opposite in direction, which can be expressed by formula (14). According to formula (14), the current contact angle can be determined as Size:

[0131] n y F1(x, t)+n y F2(x, t) = 0 ……(14)

[0132] S213: If the flow field model is at the current contact angle In the case of , the first fluid component σ is stationary on the solid component ads, and the current contact angle is determined to be the target contact angle θ.

[0133] In the embodiment of the present application, the first fluid component σ is stationary on the solid component ads, indicating that the flow field model has stabilized and the simulation of the flow field model has been completed. At this time, the target contact angle This is the target contact angle θ.

[0134] Figure 4 is a schematic diagram of a contact angle determination method provided in an embodiment of the present application, such as Figure 4 As shown, the current contact angle is determined according to the first force F1 (x, t) and the second force F2 (x, t). After that, it also includes:

[0135] S214: If the flow field model is at the current contact angle In the case of , the first fluid component σ moves on the solid component ads, and the grid density ρ and the grid velocity v of each grid point are determined according to the flow field model, the first force F1 (x, t) and the second force F2 (x, t).

[0136] In the embodiment of the present application, the grid point density ρ of each grid point is the distribution function of the grid point in each discrete velocity direction. or Specifically, the grid density ρ includes the first fluid component density ρ σ and the density of the second fluid component

[0137] For the grid points of the first fluid component σ, that is, the grid points whose component type is the first fluid component node, there is a first density ρ σ , according to the distribution function of the corresponding first fluid component σ Determine the first density ρ σ It can be expressed as formula (15):

[0138]

[0139] For the second fluid component The grid points, that is, the grid points whose component type is the second fluid component node, have the second density According to the corresponding second fluid component The distribution function of Angular second density It can be expressed as formula (16):

[0140]

[0141] In the embodiment of the present application, the grid point velocity v of each grid point is determined according to the momentum p and the grid point density ρ of the grid point, and the distribution function of the grid point is determined according to the distribution function of the grid point. or The first force F1 (x, t) and the second force F2 (x, t) determine the momentum p of the grid point.

[0142] The grid point velocity v includes the first velocity v σ and second speed Specifically, for the grid points of the first fluid component σ, there is a first velocity v σ , according to the corresponding first density ρ σ , the distribution function of the first fluid component σ First force and the second force Determine the first speed v σ , which can be expressed as formula (17):

[0143]

[0144] For the second fluid component The grid points have the second velocity According to the corresponding second density Second fluid component The distribution function of First force and the second force Determine the second speed It can be expressed as formula (18):

[0145]

[0146] S216: Update the flow field model according to the grid point density ρ and the grid point velocity v of each grid point; repeat the steps: update the flow field model according to the updated flow field model, the first parameter G 12 and the second parameter G 21 Determine the first force F1(x, t) of each grid point in the flow field, and calculate the third parameter G according to the updated flow field model. ads，1 and the fourth parameter G ads，2 Determine the second force F2(x, t) of each grid point of the flow field; determine the current contact angle according to the first force F1(x, t) and the second force F2(x, t)

[0147] In the embodiment of the present application, for a two-dimensional flow field, based on the D2Q9 model, the grid density ρ and grid velocity v of each grid point are discretized into each discrete velocity direction to obtain an updated distribution function: or Optionally, for a three-dimensional flow field, based on the D3Q19 model, the grid density ρ and grid velocity v of each grid point are discretized into each discrete velocity direction to obtain the updated distribution function or

[0148] In the embodiment of the present application, for the lattice points of the first fluid component σ, the lattice point density ρ σ and the grid point velocity v σ Discretize to each discrete velocity direction and get the updated distribution function For the second fluid component The grid points of and grid velocity Discretize to each discrete velocity direction and get the updated distribution function

[0149] According to the updated distribution function of the first fluid component σ Update the migration step and collision step of the corresponding grid points, as shown in the above formula (3) and formula (4). According to the updated second fluid composition The distribution function of Update the migration step and collision step of the corresponding grid point, as shown in the above formula (5) and formula (6). Update the distribution function of the boundary of the first fluid component σ according to the updated migration step and collision step of each grid point, thereby completing the update of the flow field model.

[0150] Repeat the steps: According to the updated flow field model, the first parameter G 12 and the second parameter G 21 Determine the first force F1(x, t) of each grid point in the flow field, and calculate the third parameter G according to the updated flow field model. ads，1 and the fourth parameter G ads，2 Determine the second force F2(x, t) of each grid point of the flow field; determine the current contact angle according to the first force F1(x, t) and the second force F2(x, t) If the flow field model is at the current contact angle In the case of the first fluid component σ being stationary on the solid component ads, determine the current contact angle is the target contact angle θ.

[0151] If the flow field model is at the current contact angle In the case of , the first fluid component σ moves on the solid component ads, and the above steps S214-S216 are repeated until the flow field model is at the current contact angle In the case of the first fluid component σ being stationary on the solid component ads, determine the current contact angle is the target contact angle θ.

[0152] Optionally, the preset number of cycles can be set to n. When the number of cycles exceeds the preset number of cycles n, the above steps are stopped to determine the current contact angle. is the target contact angle θ.

[0153] Optionally, you can also set the maximum number of loops to n. max , determine the convergence conditions, such as the current contact angle The fluctuation range is within 0.1°. When the number of cycles reaches the maximum number of cycles n max When , stop executing the above steps and determine whether the flow field model meets the convergence conditions at this time. If the convergence conditions are met, determine the current contact angle is the target contact angle θ. If the convergence condition is not met, it means that the first parameter G is selected. 12 and the second parameter G 21 It is not the optimal value, resulting in too long a number of cycles. Re-determine the first parameter G. 12 and the second parameter G 21 To redefine the target contact angle θ.

[0154] The present application also provides a contact angle determination device. Figure 5 is a schematic diagram of the structure of a contact angle determination device provided in an embodiment of the present application, such as Figure 5 As shown, the device comprises:

[0155] The initialization module 501 is used to construct a flow field and initialize a flow field model based on the flow field.

[0156] The parameter determination module 502 is used to determine a first parameter and a second parameter; the first parameter is the strength of the first fluid component on the second fluid component; the second parameter is the strength of the second fluid component on the first fluid component; and is used to determine a third parameter and a fourth parameter based on the first parameter and the second parameter; the third parameter is the strength of the solid component on the first fluid component; and the fourth parameter is the strength of the solid component on the second fluid component.

[0157] The force determination module 503 is used to determine the first force of each grid point in the flow field grid according to the flow field model, the first parameter and the second parameter, and to determine the second force of each grid point according to the flow field model, the third parameter and the fourth parameter; the first force is the force of the second fluid component and the third fluid component on the first fluid component or the force of the first fluid component and the fourth fluid component on the second fluid component, and the second force is the force of the solid component on the first fluid component or the second fluid component; the third fluid component and the first fluid component are the same fluid components, and the fourth fluid component and the second fluid component are the same fluid components.

[0158] The current contact angle determination module 504 is used to determine the current contact angle according to the first force and the second force; the current contact angle is the contact angle between the first fluid component and the solid component.

[0159] The target contact angle determination module 505 is used to determine the current contact angle as the target contact angle if the flow field model is located at the current contact angle and the first fluid component is stationary on the fixed component.

[0160] In an optional embodiment, the initialization module 501 is used to set the initial parameters of the flow field grid; the initial parameters include the solution range of the flow field grid, the component type of each grid point, the initial density of the first fluid component, the initial velocity of the first fluid component, the initial density of the second fluid component and the initial velocity of the second fluid component; the distribution function of the first fluid component is determined according to the initial density of the first fluid component and the initial velocity of the first fluid component, and the distribution function of the second fluid component is determined according to the initial density of the second fluid component and the initial velocity of the second fluid component; the migration step and collision step of each grid point are determined according to the distribution function; the distribution function of the boundary of the first fluid component is determined according to the migration step and collision step of each grid point.

[0161] In an optional embodiment, the device further comprises:

[0162] The model updating module 506 is used to determine the grid density and grid velocity of each grid point according to the flow field model, the first force and the second force when the first fluid component moves on the fixed component when the flow field model is at the current contact angle; update the flow field model according to the grid density and grid velocity of each grid point; repeat the steps of: determining the first force of each grid point of the flow field according to the updated flow field model, the first parameter and the second parameter, determining the second force of each grid point of the flow field according to the updated flow field model, the third parameter and the fourth parameter; and determining the current contact angle according to the first force and the second force.

[0163] In an optional embodiment, the model updating module 506 is used to update the distribution function of the first fluid component and the distribution function of the second fluid component according to the grid density and grid velocity of each grid point; update the migration step and collision step of each grid point according to the updated distribution function of the first fluid component and the distribution function of the second fluid component; update the distribution function of the boundary of the first fluid component according to the migration step and collision step of each grid point; update the flow field model according to the updated distribution function of the first fluid component, the distribution function of the second fluid component and the distribution function of the boundary of the first fluid component.

[0164] In an optional implementation, the initialization module 501 is used to set the boundary conditions of the first fluid component; based on the boundary conditions, the distribution function of the boundary of the first fluid component is determined according to the migration step and collision step of each grid point.

[0165] In an optional embodiment, the initialization module 501 is used to set the left and right boundaries of the first fluid component as periodic boundaries; set the lower boundary of the first fluid component as a non-slip boundary using a standard rebound format; and set the upper boundary of the first fluid component using a non-equilibrium extrapolation format.

[0166] In an optional embodiment, the parameter determination module 502 is used to determine the fifth parameter and the sixth parameter; the fifth parameter is the intensity of the effect of the third fluid component on the first fluid component; the sixth parameter is the intensity of the effect of the fourth fluid component on the second fluid component; the first pseudopotential function of the first fluid component is determined according to the fifth parameter and the state equation of the first fluid component; the second pseudopotential function of the second fluid component is determined according to the sixth parameter and the state equation of the second fluid component; the first pseudopotential function, the second pseudopotential function, the fifth parameter and the sixth parameter are used to determine the numerical relationship between the first parameter, the second parameter and the third parameter, the fourth parameter.

[0167] The device and method embodiments in the embodiments of this application are based on the same application concept.

[0168] The method embodiments provided in the embodiments of the present application can be executed in a computer terminal, a server or a similar computing device. Taking running on a server as an example, Figure 6 : is a hardware structure block diagram of a server of a contact angle determination method provided in an embodiment of the present application. Figure 6As shown, the server 600 may have relatively large differences due to different configurations or performances, and may include one or more central processing units (CPU) 610 (the processor 610 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA), a memory 630 for storing data, and one or more storage media 620 (such as one or more mass storage devices) for storing application programs 623 or data 622. Among them, the memory 630 and the storage medium 620 can be short-term storage or permanent storage. The program stored in the storage medium 620 may include one or more modules, each of which may include a series of instruction operations on the server. Furthermore, the central processing unit 610 can be configured to communicate with the storage medium 620 and execute a series of instruction operations in the storage medium 620 on the server 600. The server 600 may also include one or more power supplies 660, one or more wired or wireless network interfaces 650, one or more input and output interfaces 640, and / or one or more operating systems 621, such as Windows Server™, Mac OS X™, Unix™, Linux™, FreeBSD™, etc.

[0169] The input / output interface 640 may be used to receive or send data via a network. The specific example of the network may include a wireless network provided by a communication provider of the server 600. In one example, the input / output interface 640 includes a network adapter (Network Interface Controller, NIC), which may be connected to other network devices via a base station so as to communicate with the Internet. In one example, the input / output interface 640 may be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0170] It can be understood by those skilled in the art that Figure 6 The structure shown is for illustration only and does not limit the structure of the above electronic device. Figure 6 More or fewer components as shown, or with Figure 6 Different configurations shown.

[0171] An embodiment of the present application also provides a computer storage medium, which can be set in a server to store at least one instruction, at least one program, code set or instruction set related to a contact angle determination method in a method embodiment, and the at least one instruction, the at least one program, the code set or instruction set is loaded and executed by the processor to implement the above-mentioned contact angle determination method.

[0172] Optionally, in this embodiment, the storage medium may be located in at least one of the multiple network servers of the computer network. Optionally, in this embodiment, the storage medium may include, but is not limited to, various media that can store program codes, such as a USB flash drive, a read-only memory (ROM), a random access memory (RAM), a mobile hard disk, a magnetic disk, or an optical disk.

[0173] It can be seen from the embodiments of the contact angle method, device or storage medium provided by the present application that a flow field is constructed in the present application; the flow field includes multiple components and a flow field grid constructed based on the flow field, and the multiple components include a first fluid component, a second fluid component and a solid component; the flow field model is initialized based on the flow field; the first parameter and the second parameter are determined; the third parameter and the fourth parameter are determined according to the first parameter and the second parameter; the first force of each grid point in the flow field grid is determined according to the flow field model, the first parameter and the second parameter, and the second force of each grid point is determined according to the flow field model, the third parameter and the fourth parameter; the current contact angle is determined according to the first force and the second force; when the flow field model is located at the current contact angle, the first fluid component is stationary on the fixed component, and the current contact angle is determined to be the target contact angle. In this way, the calculation process of the contact angle can be simplified, and the contact angle simulation calculation result is stable and accurate.

[0174] It should be noted that the above-mentioned sequence of the embodiments of the present application is for description only and does not represent the advantages and disadvantages of the embodiments. The above-mentioned specific embodiments of this specification are described. Other embodiments are within the scope of the attached claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0175] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0176] A person skilled in the art will understand that all or part of the steps to implement the above embodiments may be accomplished by hardware or by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a disk or an optical disk, etc.

[0177] The above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.

Claims

1. A method for determining a contact angle, It is characterized in that The method comprises: Constructing a flow field; the flow field includes a plurality of components and a flow field grid constructed based on the flow field, the plurality of components including a first fluid component, a second fluid component and a solid component; Initializing a flow field model based on the flow field; Determine a first parameter and a second parameter; the first parameter is the intensity of the effect of the first fluid component on the second fluid component; the second parameter is the intensity of the effect of the second fluid component on the first fluid component; Determine a third parameter and a fourth parameter according to the first parameter and the second parameter; the third parameter is the strength of the action of the solid component on the first fluid component; the fourth parameter is the strength of the action of the solid component on the second fluid component; The first force of each grid point in the flow field grid is determined according to the flow field model, the first parameter and the second parameter, and the second force of each grid point is determined according to the flow field model, the third parameter and the fourth parameter; the first force is the force of the second fluid component and the third fluid component on the first fluid component or the force of the first fluid component and the fourth fluid component on the second fluid component, and the second force is the force of the solid component on the first fluid component or the second fluid component; the third fluid component and the first fluid component are the same fluid component, and the fourth fluid component and the second fluid component are the same fluid component; Determine a current contact angle according to the first force and the second force; the current contact angle is the contact angle between the first fluid component and the solid component; When the flow field model is located at the current contact angle, the first fluid component is stationary on the solid component, and the current contact angle is determined to be a target contact angle.

2. The contact angle determination method according to claim 1, It is characterized in that The initializing the flow field model based on the flow field comprises: Setting initial parameters of the flow field grid; the initial parameters include any one or more of the solution range of the flow field grid, the component type of each grid point, the initial density of the first fluid component, the initial velocity of the first fluid component, the initial density of the second fluid component and the initial velocity of the second fluid component; determining a distribution function of the first fluid component according to an initial density of the first fluid component and an initial velocity of the first fluid component, and determining a distribution function of the second fluid component according to an initial density of the second fluid component and an initial velocity of the second fluid component; Determining the migration step and the collision step of each grid point according to the distribution function of the first fluid component and the distribution function of the second fluid component; A distribution function of the boundary of the first fluid component is determined according to the migration step and the collision step of each grid point.

3. The contact angle determination method according to claim 1, It is characterized in that After determining the current contact angle according to the first force and the second force, the method further includes: If the flow field model is at the current contact angle, the first fluid component moves on the solid component, and a grid point density and a grid point velocity of each grid point are determined according to the flow field model, the first force, and the second force; Update the flow field model according to the grid point density and the grid point velocity of each grid point; repeat the steps of: determining the first force of each grid point of the flow field according to the updated flow field model, the first parameter and the second parameter, determining the second force of each grid point of the flow field according to the updated flow field model, the third parameter and the fourth parameter; determine the current contact angle according to the first force and the second force; Until the first fluid component is stationary on the solid component when the flow field model is located at the current contact angle, the current contact angle is determined to be the target contact angle.

4. The contact angle determination method according to claim 3, It is characterized in that The updating of the flow field model according to the grid point density and the grid point velocity of each grid point comprises: updating the distribution function of the first fluid component and the distribution function of the second fluid component according to the grid point density and the grid point velocity of each grid point; updating the migration step and the collision step of each grid point according to the updated distribution function of the first fluid component and the distribution function of the second fluid component; updating a distribution function of a boundary of the first fluid component according to a migration step and a collision step of each grid point; The flow field model is updated according to the updated distribution function of the first fluid component, the distribution function of the second fluid component, and the distribution function of the boundary of the first fluid component.

5. The contact angle determination method according to claim 2, It is characterized in that Determining a distribution function of a boundary of the first fluid component according to the migration and collision of each grid point includes: Setting boundary conditions for the first fluid component; Based on the boundary conditions, a distribution function of the boundary of the first fluid component is determined according to the migration step and the collision step of each grid point.

6. The contact angle determination method according to claim 5, It is characterized in that The step of setting the boundary conditions of the first fluid component comprises: Setting the left and right boundaries of the first fluid component as periodic boundaries; Setting the lower boundary of the first fluid component as a non-slip boundary, using a standard rebound format; The upper boundary of the first fluid component is set using a non-equilibrium extrapolation method.

7. The contact angle determination method according to claim 1, It is characterized in that Before determining the third parameter and the fourth parameter according to the first parameter and the second parameter, the method further includes: Determine a fifth parameter and a sixth parameter; the fifth parameter is the strength of the action of the third fluid component on the first fluid component; the sixth parameter is the strength of the action of the fourth fluid component on the second fluid component; determining a first pseudopotential function of the first fluid component according to the fifth parameter and the state equation of the first fluid component; determining a second pseudopotential function of the second fluid component based on the sixth parameter and the state equation of the second fluid component; The first pseudo potential function, the second pseudo potential function, the fifth parameter and the sixth parameter are used to determine the numerical relationship between the first parameter, the second parameter and the third parameter and the fourth parameter.

8. A contact angle determination device, It is characterized in that The device comprises: An initialization module is used to construct a flow field; the flow field includes a plurality of components and a flow field grid constructed based on the flow field, the plurality of components include a first fluid component, a second fluid component and a solid component; a flow field model is initialized based on the flow field; A parameter determination module, used to determine a first parameter and a second parameter; the first parameter is the intensity of the effect of the first fluid component on the second fluid component; the second parameter is the intensity of the effect of the second fluid component on the first fluid component; and is also used to determine a third parameter and a fourth parameter according to the first parameter and the second parameter; the third parameter is the intensity of the effect of the solid component on the first fluid component; the fourth parameter is the intensity of the effect of the solid component on the second fluid component; A force determination module, used to determine the first force of each grid point in the flow field grid according to the flow field model, the first parameter and the second parameter, and determine the second force of each grid point according to the flow field model, the third parameter and the fourth parameter; the first force is the force of the second fluid component and the third fluid component on the first fluid component or the force of the first fluid component and the fourth fluid component on the second fluid component, and the second force is the force of the solid component on the first fluid component or the second fluid component; the third fluid component and the first fluid component are the same fluid component, and the fourth fluid component and the second fluid component are the same fluid component; a current contact angle determination module, used to determine a current contact angle according to the first force and the second force; the current contact angle is a contact angle between the first fluid component and the solid component; The target contact angle determination module is used to determine the current contact angle as the target contact angle when the flow field model is located at the current contact angle and the first fluid component is stationary on the solid component.

9. An electronic device, It is characterized in that The electronic device comprises a processor and a memory, wherein the memory stores at least one instruction or at least one program, and the at least one instruction or at least one program is loaded by the processor and executes the contact angle determination method as described in any one of claims 1-7.

10. A computer storage medium, It is characterized in that The computer storage medium stores at least one instruction or at least one program, and the at least one instruction or at least one program is loaded and executed by the processor to implement the contact angle determination method according to any one of claims 1 to 7.

Images