A dynamic transfer method of finite element magnetic field data and finite volume flow field data

Abstract

The application discloses a dynamic transmission method of finite element magnetic field data and finite volume flow field data, and is used for the driving process of a magnetic fluid pump, and comprises the following steps: constructing a geometric model of the magnetic fluid pump; performing finite element calculation on the magnetic fluid micropump by using Maxwell software; transmitting the magnetic field data in the Maxwell software and the flow field data in Fluent software; adding a magnetic field source term force to Fluent flow control equations; performing coordinate system conversion to describe a circular motion traveling wave magnetic field; defining a space-time function to realize dynamic update of the magnetic field position; performing magnetic field data update and interpolation to process the position change of the dynamic magnetic field; and storing the fixed magnetic field data and the moving magnetic field data in different UDMs respectively to realize independent control of the magnetic field behavior. The application couples the magnetic field result of the finite element method with the flow field result of the finite volume method, realizes data transmission of the finite element magnetic field data and the finite volume flow field data, and realizes dynamic magnetic field control in the flow field of the finite volume method.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of flow characteristics of magnetic fluid, and particularly relates to a dynamic transfer method of finite element magnetic field data and finite volume flow field data taking a magnetic fluid pump driving process as an object. BACKGROUND

[0002] As a functional material, the magnetic fluid is widely used in many frontier fields such as magnetic fluid driving, magnetic fluid sealing, magnetic fluid sensor, magnetic fluid heat transfer and the like. However, due to the unique physical properties (combining fluid flowability, magnetic performance and opaque appearance), it is difficult to observe the flow characteristics of the magnetic fluid such as velocity and pressure distribution law through experimental methods, which limits the understanding of the flow characteristics of the magnetic fluid and hinders the further expansion and application of the magnetic fluid products.

[0003] With the development of computer technology and image processing technology, numerical calculation has become an important means to study the flow mechanism of fluid in complex structure. Some researchers try to capture the flow characteristics of ferromagnetic fluid in small-scale flow space such as sealing and driving through numerical calculation. However, the multi-field coupling (magnetic field-flow field data transfer technology) has always been a big problem in this field. Some researchers simplify the magnetic characteristics of the magnetic fluid and use numerical calculation method to predict the performance of the magnetic fluid product. Although the simulation results provide reference value for the performance of the product, the technical drawbacks are also very obvious-the motion state of the magnetic fluid under certain conditions cannot be obtained.

[0004] More importantly, some magnetic fluid products have dynamic magnetic fields such as magnetic fluid pump driving process, which involves periodic motion of permanent magnets or regular arrangement and connection of electromagnets. Under the action of the dynamic magnetic field, the magnetic fluid produces a flow state and follows the motion of the dynamic magnetic field. This scenario cannot be realized by numerical simulation method so far. The main technical problems are the data transfer of magnetic field-flow field and the control method of dynamic magnetic field.

[0005] At present, Comsol software can be used for numerical calculation of multi-field coupling in the market. However, the Comsol software is based on the finite element method to complete the model process, and for fluid, the finite volume method has obvious advantages in flow process simulation.

[0006] Based on the magnetic field results of Maxwell software and the flow field results of Fluent software, the present application proposes a numerical method of magnetic field-flow field data transfer and dynamic magnetic field control taking the magnetic fluid driving as an object, solves the technical bottleneck problem of magnetic field-flow field data transfer, and provides technical support for further research and development of magnetic fluid related products. SUMMARY

[0007] The application aims to provide a dynamic transmission method of finite element magnetic field data and finite volume flow field data.

[0008] To achieve the above-mentioned application purposes, the application adopts the following technical solutions.

[0009] A dynamic transmission method of finite element magnetic field data and finite volume flow field data comprises the following steps.

[0010] (A) Constructing a geometric model of the magnetic fluid pump, including a micro-scale flow channel of the magnetic fluid pump, a rotating permanent magnet and a fixed permanent magnet;

[0011] (B) Performing finite element calculation on the magnetic fluid micro-pump through Maxwell software, setting parameters of the fixed permanent magnet and the rotating permanent magnet, and obtaining magnetic field results;

[0012] (C) Transmitting the magnetic field data in Maxwell software and the flow field data in Fluent software;

[0013] (D) Adding a magnetic field source term force to the Fluent flow control equation, representing the influence of the magnetic field on the magnetic fluid in the form of the source term force;

[0014] (E) Coordinate system conversion, converting the magnetic field data from the Cartesian coordinate system to the polar coordinate system to describe the circular motion traveling wave magnetic field;

[0015] (F) Defining a space-time function, associating the initial magnetic field coordinates with the current time to realize dynamic update of the magnetic field position;

[0016] (G) Magnetic field data update and interpolation, assigning the magnetic field value to the fluid grid through matching or inverse distance interpolation algorithm to process the position change of the dynamic magnetic field;

[0017] (H) Magnetic field data storage allocation, storing the fixed magnetic field and the moving magnetic field data in different user-defined memories (UDM) to realize independent control of the magnetic field behavior.

[0018] Further, the specific process of step (C) is as follows:

[0019] (1) Selecting the fluid domain model of the magnetic fluid pump in the Fluent software, drawing the fluid domain grid, recording the center node position coordinates of the fluid domain grid, writing them out, and importing them into the Maxwell software;

[0020] (2) Then, the magnetic field data of the Maxwell software is exported as a txt file according to the center node coordinates of the fluid domain grid;

[0021] (3) Finally, the data in the txt file, i.e. the magnetic field results, are matched with the coordinates of the center nodes of the fluid grid, and are stored in the user-defined memory (UDM) of the center of the fluid grid, thus completing the data transfer process from the magnetic field to the flow field.

[0022] Further, in step (3), the position coordinates in the txt document are matched with the position coordinates in the Fluent grid, and the following algorithm is used:

[0023] An error value Error_1 is defined, which is one-tenth of the minimum side length of the fluid domain grid; the distance between the position coordinates in the txt document and the position coordinates in the Fluent grid is calculated, and if the distance value is less than Error_1, it is considered that the position coordinates in the txt document are the original position coordinates in the Fluent grid, and the magnetic field value is assigned to the UDM of the grid center, otherwise the distance calculation is continued for the next Fluent grid coordinates until the assignment is successful.

[0024] Further, the specific process of step (E) is as follows:

[0025] For the processing of the rotating magnetic field, the Cartesian coordinate system (x, y) in the txt document is first converted into the polar coordinate system (ρ r , θ) to describe the circular motion traveling wave magnetic field, and the corresponding conversion relationship is shown in equations (1) and (2).

[0026]

[0027] In the equations, ρ r represents the polar radius, and θ represents the polar angle.

[0028] Further, the specific process of step (F) is as follows:

[0029] The spatial coordinates are represented as a function of time, with the coordinate point (x0, y0) being the UDM grid center coordinate at the initial time when the magnetic field data is stored, and the coordinate point (x t , y t ) being the UDM grid center coordinate at the current time when the magnetic field data is stored. According to the mathematical relationship of polar coordinates, combined with the trigonometric function and the difference integral formula, the initial coordinate point (x0, y0) is expressed as the function of the current coordinate point (x t , y t ) and time t using equations (3) and (4), which is called the space-time function.

[0030]

[0031] In the equations, t is the running time of the magnetic field, and ω is the running angular velocity.

[0032] Further, the specific process of step (G) is as follows:

[0033] This step is to contain the initial magnetic field data in the flow field grid, and the magnetic field data storage location changes over time; after each time iteration process, a new set of magnetic field position coordinates will be obtained, and the new magnetic field position coordinates will again traverse all fluid domain grids to find the fluid grid that is infinitely close in space position, when the spatial position coordinates calculated by the space-time function do not meet the distance discrimination standard with all fluid grid centers, a search radius r e The inverse distance interpolation algorithm is used to weight average within the search radius, and the weighted average magnetic field value is assigned to the surrounding fluid grid, and the inverse distance weight coefficient w i The calculation method is shown in formula (5), and the weighted average processing method is shown in formula (6);

[0034]

[0035] In the formula, (x f , y f ) represents the fluid grid center coordinates, (x i , y i ) represents the magnetic field coordinates, Nc represents the number of grids around the target grid center, p c is the calculation coefficient, and H i is the magnetic field strength of the surrounding grid; according to formula (6), the weighted average magnetic field strength H is obtained.

[0036] Further, the specific process of step (H) is as follows:

[0037] The fixed magnetic field data is directly stored in the fluid grid center UDM0, the magnetic field derivative of the fixed magnetic field data in the x direction is stored in UDM1, the magnetic field derivative in the y direction is stored in UDM2, and the magnetic field derivative in the z direction is stored in UDM3;

[0038] The moving magnetic field data is recorded in the fluid grid center UDM4 in a dynamic magnetic field-flow field coupling manner, the magnetic field derivative of the moving magnetic field data in the x direction is stored in UDM5, the magnetic field derivative in the y direction is stored in UDM6, and the magnetic field derivative in the z direction is stored in UDM7;

[0039] At this time, the program for controlling the movement of the magnetic field only acts on UDM5, UDM6 and UDM7, and the magnetic field of UDM1, UDM2 and UDM3 remains unchanged, and similar magnetic field behavior to the actual micro-pump driving process is realized in numerical calculation.

[0040] The beneficial effects of the present application are:

[0041] The application is based on the magnetic field-flow field results of the finite element method Maxwell software and the finite volume method Fluent software, and the data sharing process of the magnetic field-flow field is completed through data processing technology, not only the data transmission of the finite element magnetic field data-flow field data of the finite volume is realized, but also the dynamic magnetic field control in the flow field of the finite volume method is realized.

[0042] The application couples the magnetic field results of the finite element method with the flow field results of the finite volume method, guarantees the reliability of the magnetic field calculation, guarantees the high precision of the fluid calculation, and is a novel high-precision magnetic field-flow field coupling simulation method. BRIEF DESCRIPTION OF DRAWINGS

[0043] Figure 1 It is a schematic diagram of the numerical model of the magnetic fluid pump and the actual driving process;

[0044] Figure 2 It is a schematic diagram of the magnetic field results calculated by the Maxwell software, wherein (a) is the magnetic field results of the fixed permanent magnet, and (b) is the magnetic field results of the rotating permanent magnet;

[0045] Figure 3 It is a schematic diagram of the comparison between the magnetic field results calculated by the Maxwell software and the magnetic field results imported into the Fluent software;

[0046] Figure 4 It is a schematic diagram of the dynamic magnetic field operation process;

[0047] Figure 5 It is a schematic diagram of the comparison between the experimental results and the numerical results of the macroscopic flow phenomenon, wherein (a) is a schematic diagram of the experimental results of the actual magnetic fluid pump, and (b) is a schematic diagram of the numerical results of the method of the application;

[0048] Figure 6 It is a schematic diagram of the comparison between the experimental results and the numerical results of the magnetic fluid pump at a rotation speed of 1 rpm.

[0049] The drawings are only used for illustrative description, and cannot be understood as a limitation on the application; in order to better illustrate the embodiment, some components in the drawings may be omitted, enlarged or reduced, and do not represent the size of the actual product; for those skilled in the art, it is understandable that some known structures and their descriptions in the drawings may be omitted. DETAILED DESCRIPTION

[0050] The technical solutions in the embodiments of the application will be clearly and completely described below with reference to the drawings in the embodiments of the application.

[0051] As shown in the drawings, Figure 1 The embodiment discloses a dynamic transmission method of finite element magnetic field data and finite volume flow field data, which comprises the following steps:

[0052] (A) First, the geometric model of the magnetic fluid pump is constructed, including the magnetic fluid pump micro-scale flow channel, rotating permanent magnet and fixed permanent magnet, as shown in Figure 1 .

[0053] (B) The magnetic fluid micro-pump is calculated by Maxwell software to obtain the magnetic field results, the fixed permanent magnet and the rotating permanent magnet are set to residual magnetism 400 Gs and coercive force 837000 A / m, the fixed permanent magnet and the rotating permanent magnet are arranged according to the position shown in the figure, and the movement process of the rotating permanent magnet is not considered.

[0054] Figure 2 The obtained fixed permanent magnet magnetic field and rotating permanent magnet magnetic field results, wherein (a) is the magnetic field result of the fixed permanent magnet magnetic field, and (b) is the magnetic field result of the rotating permanent magnet.

[0055] (C) The transfer of the magnetic field data in Maxwell software and the flow field data in Fluent software is performed.

[0056] The transfer process is divided into three steps: (1) selecting the magnetic fluid pump fluid domain model in the Fluent software, drawing the fluid domain grid, recording the center node position coordinates of the fluid domain grid, writing out, and importing into the Maxwell software; (2) then exporting the magnetic field data of the Maxwell software into a txt file according to the center node coordinates of the fluid domain grid; (3) finally, the data (magnetic field results) in the txt file are matched according to the center node coordinates of the fluid grid, and stored in the user-defined memory (User-Defined Memory, UDM) of the fluid grid center, completing the data transfer process from the magnetic field to the flow field.

[0057] In the third step, the position coordinates in the txt document need to be matched with the position coordinates in the Fluent grid, so a specific algorithm needs to be used. In the present application, an error value Error_1 is defined, which is one tenth of the minimum side length of the fluid domain grid. The distance between the position coordinates in the txt document and the position coordinates in the Fluent grid is calculated, if the distance value is less than Error_1, it is considered that the position coordinates in the txt document are the original position coordinates in the Fluent grid, and the magnetic field value is assigned to the UDM of the grid center, otherwise the distance calculation of the next Fluent grid coordinates is continued until the assignment is successful.

[0058] Because all the position coordinates in the txt document are provided by the Fluent grid model, the matching position coordinates are equal, that is, the error value is close to 0, and defining Error_1 to be one-tenth of the minimum grid length can ensure that each Fluent grid is assigned only once, and the remaining txt document position coordinates do not meet the condition.

[0059] Figure 3 The comparison of the magnetic field results calculated by Maxwell software and imported into Fluent software shows that the magnetic field results of the two are almost consistent.

[0060] (D) Adding a magnetic field source term force to the Fluent flow control equation to represent the influence of the magnetic field on the magnetic fluid in the form of a source term force.

[0061] (E) Coordinate system conversion, converting the magnetic field data from the Cartesian coordinate system to the polar coordinate system to describe the circular motion of the traveling wave magnetic field. For the processing of rotating magnetic fields (dynamic magnetic fields), first convert the Cartesian coordinate system (x, y) in the txt document to the polar coordinate system (ρ r , θ), to describe the circular motion of the traveling wave magnetic field. The corresponding conversion relationship is shown in equations (1) and (2), where ρ r represents the polar radius, and θ represents the polar angle.

[0062]

[0063] (F) Define a space-time function to associate the initial magnetic field coordinates with the current time to achieve dynamic updating of the magnetic field position.

[0064] Express the spatial coordinates in terms of the time function. Let the coordinate point (x0, y0) be the UDM grid center coordinate of the stored magnetic field data at the initial time, and the coordinate point (x t , y t ) be the UDM grid center coordinate of the stored magnetic field data at the current time. The magnetic field running time is t, and the running angular velocity is ω. According to the polar coordinate mathematical relationship, combined with the trigonometric function and the difference integral formula, the initial coordinate point (x0, y0) can be expressed as the current coordinate point (x t , y t ) and the function of time t, which is called a space-time function.

[0065]

[0066] (G) Magnetic field data updating and interpolation, assigning the magnetic field value to the fluid grid through matching or inverse distance interpolation algorithm to process the position change of dynamic magnetic field.

[0067] ​At this time, the initial magnetic field data is contained in the flow field grid, and the magnetic field data storage location changes over time. After each time iteration, a new set of magnetic field position coordinates is obtained, which is again traversed through all fluid domain grids to find the fluid grid that is infinitely close to its spatial position (infinitely close is less than Error_1). When the spatial position coordinates calculated by the space-time function do not meet the criteria for the distance from all fluid grid centers, a search radius r e is set, and the weighted average magnetic field value is assigned to the surrounding fluid grid using the inverse distance interpolation algorithm within the search radius. The inverse distance weight coefficient w i is calculated as shown in equation (5), and the weighted average processing method is shown in equation (6).

[0068]

[0069] where (x f , y f ) represents the fluid grid center coordinates, (x i , y i ) represents the magnetic field coordinates, Nc represents the number of grids around the target grid center, p c is the calculation coefficient, which is 2, H i is the surrounding grid magnetic field strength, and the weighted average magnetic field strength H can be obtained according to equation (6). The specific process is shown in Figure 4 .

[0070] (H) Magnetic field data storage allocation, fixed magnetic field and moving magnetic field data are stored in different user-defined memories (UDM), which realizes independent control of the magnetic field behavior.

[0071] The fixed magnetic field data is directly stored in the fluid grid center UDM0 (the magnetic field derivative in the x direction of the fixed magnetic field data is stored in UDM1, the magnetic field derivative in the y direction is stored in UDM2, and the magnetic field derivative in the z direction is stored in UDM3), and the moving magnetic field data is recorded in the fluid grid center UDM4 in a dynamic magnetic field-flow field coupling manner (the magnetic field derivative in the x direction of the moving magnetic field data is stored in UDM5, the magnetic field derivative in the y direction is stored in UDM6, and the magnetic field derivative in the z direction is stored in UDM7). At this time, the program that controls the magnetic field motion only acts on UDM5, UDM6, and UDM7, while the magnetic field of UDM1, UDM2, and UDM3 remains unchanged, which can realize similar magnetic field behavior to the actual micro-pump driving process in numerical calculation.

[0072] The present application takes the magnetic fluid pump (micro-pump) as an example to experimentally verify the magnetic field-flow field coupling method and the dynamic magnetic field control method.

[0073] Firstly, the macroscopic flow phenomenon driven by magnetic fluid, such as Figure 5 Figure 5 Fig. 1 shows the experimental results of the physical magnetic fluid pump (a) and the numerical results of the method of the present application (b), and it can be seen that the numerical results are almost consistent with the experimental results and highly similar in microscopic flow details, verifying the feasibility of the method of the present application.

[0074] In addition, the present application performs numerical calculation on the transport speed of the magnetic fluid pump and compares it with the experimental data of the laser Doppler velocimeter.

[0075] As shown in Fig. 2, due to the long time consumption of numerical calculation, the results only involve 0-0.7 periods, but the calculation time domain involves all main processes of the micropump. Figure 6 The numerical results of the method of the present application and the experimental results of the physical magnetic fluid pump both show strong speed pulsation in the ferromagnetic fluid fusion stage, and due to the addition of the rectangular glass tube as the observation flow channel in the experiment, the long-distance transport of liquid has a buffering effect on the liquid pulsation, so the strong pulsation near 0.3 period is not captured in the experiment.In the remaining time domain, the speed value of the micropump outlet is small-amplitude pulsation around a certain average value, and the numerical results and the experimental results are consistent in the pulsation law of the speed-time curve, and the speed amplitude is close, verifying the correctness of the numerical model. Therefore, although the time domain of the numerical results is not completely coincided with the time domain of the experimental results, it is still scientific, further verifying the feasibility of the method of the present application.

[0076] The above examples are only used to illustrate but not to limit the technical solutions of the present application, and although the present application has been described in detail with reference to the above examples, those skilled in the art should understand that the present application can still be modified or equivalently replaced without departing from the spirit and scope of the present application, and any modification or partial replacement should be covered in the scope of the claims of the present application.

[0077] If the words "first", "second", etc. are used to limit parts in this document, those skilled in the art should know that the use of "first", "second" is only for the convenience of describing the present application and simplifying the description, and the above words have no special meaning unless otherwise stated.

[0078] In the description of the application, it needs to be understood that the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential" and the like indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, and are only for the purpose of facilitating the description of the application and simplifying the description, and do not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation on the application.

[0079] In the description of the application, it needs to be understood that unless otherwise explicitly specified and limited, the terms "mounting", "connecting", "connecting" should be understood broadly, for example, it can be fixedly connected, or it can be detachably connected, or integrally connected; it can be mechanically connected, or it can be electrically connected; it can be directly connected, or it can be indirectly connected through an intermediate medium; it can be the communication inside two elements. For those skilled in the art, the specific meaning of the above terms in the application can be understood according to the specific circumstances.

Claims

1. A method of dynamic transfer of finite element magnetic field data and finite volume flow field data, characterized by It comprises the following steps: (A) constructing a geometric model of the magnetic fluid pump, including the micro-scale flow channel of the magnetic fluid pump, the rotating permanent magnet and the fixed permanent magnet; (B) performing finite element calculation on the magnetic fluid micro-pump by Maxwell software, setting parameters of the fixed permanent magnet and the rotating permanent magnet, and obtaining magnetic field results; (C) transferring magnetic field data in Maxwell software to flow field data in Fluent software; (D) adding a magnetic field source term to the flow control equation of Fluent, representing the influence of the magnetic field on the magnetic fluid in the form of a source term; (E) coordinate system conversion, converting the magnetic field data from Cartesian coordinates to polar coordinates to describe the circular motion of the traveling wave magnetic field; (F) defining a space-time function to associate the initial magnetic field coordinates with the current time and dynamically update the magnetic field position; (G) magnetic field data updating and interpolation, assigning magnetic field values to fluid grids by matching or inverse distance interpolation algorithm to handle the position changes of dynamic magnetic field; (H) magnetic field data storage allocation, storing fixed magnetic field and moving magnetic field data in different user-defined memories (UDMs) to realize independent control of magnetic field behavior; The specific process of step (C) is as follows: (1) selecting the fluid domain model of the magnetic fluid pump in Fluent software, drawing the fluid domain grid, recording the center node position coordinates of the fluid domain grid, writing them out, and importing them into Maxwell software; (2) then exporting the magnetic field data of Maxwell software as a txt file according to the center node coordinates of the fluid domain grid; (3) finally, matching the data in the txt file, i.e. the magnetic field results, with the center node coordinates of the fluid grid, and storing them in the user-defined memory (UDM) of the fluid grid center to complete the data transfer process from magnetic field to flow field; In step (3), the position coordinates in the txt document are matched with the position coordinates in the Fluent grid, and the following algorithm is used: Define an error value Error_1, which is one tenth of the minimum side length of the fluid domain grid; calculate the distance between the position coordinates in the txt document and the position coordinates in the Fluent grid, if the distance value is less than Error_1, then the position coordinates in the txt document are considered as the original position coordinates in the Fluent grid, and the magnetic field value is assigned to the grid center UDM, otherwise, continue to calculate the distance for the next Fluent grid coordinates until the assignment is successful; The specific process of step (E) is as follows: For the processing of rotating magnetic field, the Cartesian coordinate system (x, y) in the txt document is first converted into polar coordinate system (p r , θ) to describe the circular motion traveling wave magnetic field, and the corresponding conversion relationship is shown in equations (1) and (2). wherein p r denotes the polar radius, denotes the polar angle; The specific process of step (F) is as follows: The space coordinates are expressed as a function of time, the coordinate point (x0, y0) is the center coordinate of the UDM grid storing the magnetic field data at the initial time, the coordinate point (x t , y t ) is the center coordinate of the UDM grid storing the magnetic field data at the current time, according to the mathematical relationship of polar coordinates, combined with the trigonometric function and the difference integral formula, the initial coordinate point (x0, y0) is expressed as the current coordinate point (x t , y t ) and the function of time t by equations (3) and (4), which is called a space-time function; Where t is the magnetic field running time and ω is the running angular velocity.

2. The method of dynamic transfer of finite element magnetic field data and finite volume flow field data according to claim 1, wherein, The specific process of step (G) is as follows: This step contains the initial time magnetic field data in the flow field grid, and the magnetic field data storage location changes over time; after each time iteration process, a new set of magnetic field position coordinates will be obtained, and the new magnetic field position coordinates will again traverse all fluid domain grids to find the fluid grid that is infinitely close to its spatial position. When the spatial position coordinates calculated by the space-time function do not meet the distance criterion for all fluid grid centers, set a search radius r e Within the search radius, use the inverse distance interpolation algorithm for weighted average, and assign the weighted average magnetic field value to the surrounding fluid grid. The inverse distance weight coefficient w i The calculation method is shown in formula (5), and the weighted average processing method is shown in formula (6); where (x f , y f ) represents the fluid grid center coordinates, (x i , y i ) represents the magnetic field coordinates, Nc represents the number of surrounding grids of the target grid center, p c is a calculation coefficient, and H i is the surrounding grid magnetic field strength; the weighted average magnetic field strength H is obtained according to formula (6).

3. The method of dynamic transfer of finite element magnetic field data and finite volume flow field data according to claim 2, wherein, The specific process of step (H) is as follows: The fixed magnetic field data is directly stored in the fluid grid center UDM0, the magnetic field derivative of the fixed magnetic field data in the x direction is stored in UDM1, the magnetic field derivative in the y direction is stored in UDM2, and the magnetic field derivative in the z direction is stored in UDM3; The moving magnetic field data is recorded in the fluid grid center UDM4 in a dynamic magnetic field-flow field coupling mode, the magnetic field derivative in x direction is stored in UDM5, the magnetic field derivative in y direction is stored in UDM6, and the magnetic field derivative in z direction is stored in UDM7; At this time, the program for controlling the movement of the magnetic field only acts on UDM5, UDM6 and UDM7, while the magnetic field of UDM1, UDM2 and UDM3 remains unchanged, and similar magnetic field behavior to the actual micro-pump driving process is realized in numerical calculation.

Images