Particle motion simulation method, device and equipment and storage medium

Abstract

The invention provides a particle motion simulation method, device and equipment and a storage medium, and relates to the technical field of hydromechanics, and the method comprises the steps: constructing a geometric model of an inhaler, and carrying out the meshing of the geometric model; determining initial parameter information of particulate matters in the geometric model; determining fluid parameter information of the fluid phase; and according to the initial parameter information and the fluid parameter information, utilizing a pre-constructed motion simulation model to solve the motion trail of the particulate matter in the inhaler. According to the invention, the geometric model of the inhaler is constructed by using OpenFOAM software, and gas and solid are captured at the same time according to the initial parameter information and fluid parameter information of the particulate matters, so that the movement of the particulate matters and the movement track of the landing position after stopping are simulated.

Description

Technical Field

[0001] The present invention relates to the field of fluid mechanics technology, and in particular to a particle motion simulation method, device, equipment and storage medium. Background Art

[0002] The Dose Per Inhalation (DPI) is a medical device used to inhale solid drug powder directly into the lungs. The performance evaluation of the powder inhaler has always been a key link and needs to rely on actual experimental tests for evaluation. These experiments are designed to directly observe the dispersibility and deposition efficiency of drug particles, as well as the aerodynamic behavior of patients when using the inhaler. Although these experiments can provide data that is highly consistent with the actual usage scenario, thereby accurately evaluating the performance of the inhaler, the evaluation is often costly and time-consuming. With the development of CFD technology, some preliminary numerical simulations have begun to be applied to the performance prediction of powder inhalers. The simulation process usually requires the use of commercial software, such as ANSYS Fluent, COMSOL Multiphysics and other commercial software, which provide powerful solvers and user interfaces, but the cost is relatively high, and the flexibility may be limited by the customization options of the software vendor. In the development process of traditional powder inhalers, a large number of experimental tests are required to evaluate their performance, which makes the experimental tests not only costly but also time-consuming. Summary of the invention

[0003] The present invention provides a particle motion simulation method, device, equipment and storage medium, aiming to solve at least one of the above-mentioned technical problems.

[0004] The present invention provides a particle motion simulation method, comprising:

[0005] constructing a geometric model of the inhaler and meshing the geometric model;

[0006] Determining initial parameter information of particles in the geometric model;

[0007] determining fluid parameter information of the fluid phase;

[0008] According to the initial parameter information and the fluid parameter information, a pre-constructed motion simulation model is used to solve the motion trajectory of the particles in the inhaler.

[0009] According to a particle motion simulation method provided by the present invention, the step of constructing a geometric model of an inhaler and meshing the geometric model comprises:

[0010] According to the model size of the inhaler, a plurality of coordinate positions are determined;

[0011] Determining a physical space of the geometric model according to the multiple coordinate positions;

[0012] The number of grids and the grid expansion rate are determined in the physical space to grid the geometric model.

[0013] According to a particle motion simulation method provided by the present invention, solving the motion trajectory of particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information includes:

[0014] Construct the NS equations, turbulence models, and equations of motion for particles;

[0015] According to the initial parameter information and the fluid parameter information, physical parameters are calculated using the NS equations and the turbulence model;

[0016] The motion trajectory of the particles in the inhaler is obtained by solving the motion equation according to the physical parameters, the initial parameter information and the fluid parameter information.

[0017] According to a particle motion simulation method provided by the present invention, after constructing the geometric model of the inhaler, the method further comprises:

[0018] Determine geometric boundary conditions corresponding to the geometric model, wherein the geometric boundary conditions include pressure values ​​corresponding to an inlet and an outlet in the geometric model, and the pressure values ​​corresponding to the inlet and the outlet are used as parameters for solving the motion trajectory of the particles in the inhaler.

[0019] According to a particle motion simulation method provided by the present invention, the initial parameter information includes at least one of the starting emission position of the particle, the initial velocity field, the initial pressure field, the volume fraction boundary conditions of air and particles, the magnitude and direction of gravity, fluid properties, control parameters, parallel parameters and discrete formats; wherein the volume fraction of the particle refers to the content of the particle;

[0020] The fluid parameter information includes at least one of kinematic viscosity, fluid density, surface tension, and whether the fluid is a Newtonian fluid.

[0021] According to a particle motion simulation method provided by the present invention, the starting emission position of the particle is determined according to the exit position of the capsule compartment in the geometric model.

[0022] The present invention also provides a particle motion simulation device, comprising:

[0023] A construction module, used for constructing a geometric model of the inhaler and meshing the geometric model;

[0024] A first determination module, used to determine initial parameter information of particles in the geometric model;

[0025] A second determination module, used to determine fluid parameter information of the fluid phase;

[0026] The solution module is used to solve the motion trajectory of the particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information.

[0027] According to a particle motion simulation device provided by the present invention, the solution module is also used for:

[0028] Construct the NS equations, turbulence models, and equations of motion for particles;

[0029] According to the initial parameter information and the fluid parameter information, physical parameters are calculated using the NS equations and the turbulence model;

[0030] The motion trajectory of the particles in the inhaler is obtained by solving the motion equation according to the physical parameters, the initial parameter information and the fluid parameter information.

[0031] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, any of the particle motion simulation methods described above is implemented.

[0032] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the particle motion simulation method as described above is implemented.

[0033] The present invention also provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the computer program implements any of the above-mentioned particle motion simulation methods.

[0034] The particle motion simulation method, device, equipment and storage medium provided by the present invention include: constructing a geometric model of an inhaler, and meshing the geometric model; determining initial parameter information of particles in the geometric model; determining fluid parameter information of a fluid phase; and solving the motion trajectory of particles in the inhaler using a pre-constructed motion simulation model according to the initial parameter information and the fluid parameter information. The present invention constructs a geometric model of the inhaler using OpenFOAM software, combines the initial parameter information of the particles and the fluid parameter information to achieve simultaneous capture of gas and solid, simulates the motion of particles and the motion trajectory of the landing position after stopping, and then determines the performance of the inhaler according to the motion trajectory of the particles, thereby improving the efficiency of the inhaler test, and does not need to use commercial software such as ANSYS Fluent and COMSOL Multiphysics, and has low cost. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced one by one below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0036] Figure 1 It is a schematic diagram of the process of the particle motion simulation method provided by the present invention;

[0037] Figure 2 It is a structural schematic diagram of the particle motion simulation device provided by the present invention;

[0038] Figure 3 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0039] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0040] The terms used in one or more embodiments of the present invention are only for the purpose of describing specific embodiments, and are not intended to limit one or more embodiments of the present invention. The singular forms of "a", "said" and "the" used in one or more embodiments of the present invention are also intended to include plural forms, unless the context clearly indicates other meanings. It should also be understood that the term "and / or" used in one or more embodiments of the present invention refers to and includes any or all possible combinations of one or more associated listed items.

[0041] It should be understood that, although the terms first, second, etc. may be used to describe various information in one or more embodiments of the present invention, these information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of one or more embodiments of the present invention, the first may also be referred to as the second, and similarly, the second may also be referred to as the first. Depending on the context, the word "if" as used herein may be interpreted as "when..." or "when...".

[0042] The Dose Per Inhalation (DPI) is a medical device used to inhale solid drug powder directly into the lungs. The performance evaluation of the powder inhaler has always been a key link and needs to rely on actual experimental tests for evaluation. These experiments are designed to directly observe the dispersibility and deposition efficiency of drug particles, as well as the aerodynamic behavior of patients when using the inhaler. Although these experiments can provide data that is highly consistent with the actual usage scenario, thereby accurately evaluating the performance of the inhaler, the evaluation is often costly and time-consuming. With the development of CFD technology, some preliminary numerical simulations have begun to be applied to the performance prediction of powder inhalers. The simulation process usually requires the use of commercial software, such as ANSYS Fluent, COMSOL Multiphysics and other commercial software, which provide powerful solvers and user interfaces, but the cost is relatively high, and the flexibility may be limited by the customization options of the software vendor. In the development process of traditional powder inhalers, a large number of experimental tests are required to evaluate their performance, which makes the experimental tests not only costly but also time-consuming.

[0043] In view of the above problems, the present invention provides the following embodiments. Specifically, Figure 1 Schematic diagram of the particle motion simulation method provided by the present invention. Figure 1 As shown, the particle motion simulation method includes:

[0044] Step S11, constructing a geometric model of the inhaler and meshing the geometric model;

[0045] It should be noted that the root inhaler model is surrounded by a wall, an inlet and an outlet, and the inhaler contains a capsule compartment.

[0046] Specifically, according to the model size of the inhaler, multiple coordinate positions are determined; for example, 8 coordinate positions are extracted, wherein the coordinate positions of the 8 points are based on the model size of the inhaler. In a specific example, the xyz coordinates of the eight points are as follows: vertices((-0.0210474, -0.036108, -0.0098439), (0.0210474, -0.036108, -0.0098439), (-0.0210474, 0.0310308, -0.00984 39), (0.0210474, 0.0310308, -0.0098439), (-0.0210474, -0.036108, 0.0098439), (0.0210474, -0.036108, 0.0098439), (-0.0210474, 0.0310308, 0.0098439), (0.0210474, 0.0310308, 0.0098439)), the coordinate positions of the 8 points can completely contain the inhaler model. Then, according to the multiple coordinate positions, the physical space of the geometric model is determined. Furthermore, in order to perform numerical simulation, the computational domain needs to be meshed. For example, the number of meshes and the mesh expansion rate are set in the z-axis direction (vertical direction). In a specific example, hex(0, 1, 3, 2, 4, 5, 7, 6) defines the index of the eight vertices constituting the hexahedral mesh block; (80, 180, 80) defines the number of division units of the mesh block in the x, y, and z directions, that is, the hexahedral mesh block is divided into 80x180x80 small hexahedral units. simpleGrading(1 1 1) is used to control the non-uniform distribution of the mesh to achieve meshing of the geometric model.

[0047] In addition, in other embodiments, it is necessary to determine the geometric boundary conditions corresponding to the geometric model, wherein the geometric boundary conditions include the pressure values ​​corresponding to the inlet and outlet in the geometric model, optionally, the inlet and outlet pressure difference is constant to 4kpa, and the pressure values ​​corresponding to the inlet and outlet are used as the parameters for solving the motion trajectory of particles in the inhaler. Understandably, it is necessary to set the parameters of the inlet and outlet in the computational domain. For example, the simulation operating condition is that the pressure difference between the inlet and the outlet needs to be constant at 4kpa, therefore, when setting the set boundary conditions, the inlet needs to be set to a pressure, and the outlet sets a pressure value, thereby realizing the purpose of the pressure difference 4kpa.

[0048] In addition, the multiphase fluid domain of the capsule chamber is set. For example, the inlet of the inhaler is set as the air inlet area. Air enters from the inlet and drives the powder in the capsule chamber to move after passing through the capsule chamber.

[0049] Step S12, determining initial parameter information of particles in the geometric model;

[0050] Step S13, determining fluid parameter information of the fluid phase;

[0051] It should be noted that the initial parameter information includes the starting launch position of the particles, the initial velocity field, the initial pressure field, the volume fraction boundary conditions of air and particles, the magnitude and direction of gravity, fluid properties, control parameters, parallel parameters and discrete formats, among which the starting launch position of the particles is determined according to the exit position of the capsule compartment in the geometric model. The initial velocity field refers to the velocity distribution of the fluid and particles at various positions at the beginning of the simulation. The initial pressure field refers to the pressure distribution of the fluid at various positions at the beginning of the simulation, and the pressure difference in the initial pressure field is set to 4kpa; the volume fraction of particles refers to the content of particles; the volume fraction boundary condition of air is the content of air in the calculation domain, and optionally, the volume fraction of particles is 80%. The magnitude of gravity is 9.81 dimensions [0 1 -2 0 0 0 0]: This is the dimensional representation of gravitational acceleration. In OpenFOAM, dimensions are usually represented by a 7-dimensional vector, corresponding to mass (M), length (L), time (T), temperature (Θ), amount of substance (N), current (I) and light intensity (J). [0 1 -2 0 0 0 0] represents the massless dimension (0), length dimension (1), time dimension (-2), and the remaining dimensions are all 0, i.e. m / s2. Fluid properties refer to the physical and chemical properties of the fluid, including density, viscosity, thermal conductivity, etc. Control parameters refer to various parameters used to control the simulation process, including time step, number of iterations, convergence criterion, etc. Control parameters: The denseParticleFoam solver is a solver in OpenFOAM that is specifically used to simulate dense particle flow. The denseParticleFoam solver can handle the movement of particles in fluids, including collisions between particles and fluid-particle interactions. Parallel parameters refer to various parameters used for parallel computing, including the number of processors, task allocation strategy, etc. Discrete format refers to the method of discretizing continuous physical equations, including spatial discretization and temporal discretization.

[0052] It should be noted that the fluid parameter information includes information such as kinematic viscosity, fluid density, surface tension, and whether it is a Newtonian fluid. Among them, kinematic viscosity is the ratio of the dynamic viscosity of the fluid to its density, indicating the ease with which the fluid flows under the action of gravity. Surface tension is the tension generated on the surface of a liquid due to the mutual attraction between molecules. A Newtonian fluid refers to a fluid whose shear stress is linearly related to the shear rate.

[0053] Alternatively, refer to the following specific example for numerical simulation in OpenFOAM, involving the movement and distribution of particles in a fluid. These parameters define the characteristics of the normal distribution and the physical properties of the fluid.

[0054] Defines the parameters of a normal distribution, including expectation (mean), variance, minimum and maximum values. These parameters are often used to generate random numbers that conform to a normal distribution in numerical simulation, statistical analysis, or random process simulation. normalDistribution / / declares a normal distribution type or object

[0055] {

[0056] expectation 100e-6; / / defines the expectation (mean) of the normal distribution, set to 100 microns (100e-6 meters)

[0057] variance 25e-6; / / defines the variance of the normal distribution, set to 25 microns squared (where the variance unit is the square of the original unit)

[0058] minValue 20e-6; / / defines the minimum value when generating random numbers, set to 20 microns

[0059] maxValue 180e-6; / / defines the maximum value when generating random numbers, set to 180 microns}

[0060] It can be understood that the expectation (mean) is the central value of the normal distribution, which represents the average value of the particle size. Variance describes the dispersion of the particle size. The larger the variance, the wider the distribution of the particle size. Minimum and maximum values ​​limit the range of the generated random numbers to ensure that the generated particle size is within a reasonable range.

[0061] The process of setting fluid property parameters is as follows:

[0062] continuousPhaseName air; / / Set the name of the continuous phase to "air", that is, the fluid in the simulation is air.

[0063] rho.air 1.2; / / Set the density (rho) of air (air) to 1.2kg / m3.

[0064] transportModel Newtonian; / / Set the transport model of the fluid to "Newtonian", that is, Newtonian fluid. Newtonian fluid refers to a fluid that satisfies Newton's viscosity law, and its shear stress is proportional to the shear strain rate.

[0065] nu 1.568e-05; / / Set the dynamic viscosity (nu, also known as viscosity coefficient) of the fluid to 1.568e-05m2 / s. Dynamic viscosity is a physical quantity that describes the viscosity of a fluid, which determines the flow behavior of the fluid under shear force.

[0066] Step S14, solving the motion trajectory of the particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information.

[0067] Specifically, the Euler method is used to solve the equations of motion of the air (such as the NS equations) and calculate the velocity field, pressure field, etc. of the air. Particle tracking: Use the Lagrangian method to track the movement of powder particles. The equations of motion of particles usually include the effects of aerodynamics, gravity, buoyancy and other forces. In addition, the movement of particles needs to be coupled with the flow of air to ensure that the movement of particles affects the flow of air, and the flow of air also affects the movement of particles. Optionally, both modeling and calculation work are completed using the open source software Openfoam, and the source code program can be secondary developed and modified, reducing the acquisition cost of commercial software.

[0068] The embodiment of the present invention, through the above scheme, includes: constructing a geometric model of the inhaler, and meshing the geometric model; determining the initial parameter information of the particles in the geometric model; determining the fluid parameter information of the fluid phase; and solving the motion trajectory of the particles in the inhaler using a pre-constructed motion simulation model according to the initial parameter information and the fluid parameter information. The embodiment of the present invention constructs the geometric model of the inhaler by using OpenFOAM software, combines the initial parameter information of the particles and the fluid parameter information to achieve simultaneous capture of gas and solid, simulates the motion of the particles and the motion trajectory of the landing position after stopping; and then determines the performance of the inhaler according to the motion trajectory of the particles, thereby improving the efficiency of the inhaler test, and does not need to use commercial software such as ANSYS Fluent and COMSOL Multiphysics, and the cost is low.

[0069] In one embodiment of the present invention, solving the motion trajectory of particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information includes:

[0070] Constructing a group of NS equations, a turbulence model and a motion equation of particles; calculating physical parameters using the group of NS equations and the turbulence model according to the initial parameter information and the fluid parameter information; solving the motion equation according to the physical parameters, the initial parameter information and the fluid parameter information to obtain the motion trajectory of the particles in the inhaler.

[0071] It should be noted that the turbulence model can be set as a k-ε model, where the k-ε model consists of two main equations, the k equation and the ε equation. The formula form of the k equation (turbulent kinetic energy equation) is as follows:

[0072]

[0073] Where ρ is the fluid density, k is the turbulent kinetic energy, which represents the average turbulent kinetic energy per unit mass of fluid, t is the time, and u i is the velocity component, x i is the spatial coordinate, μ is the molecular viscosity, μ t is the turbulent viscosity, σ k is the turbulent Prandtl number of the k-equation, G k is the turbulent kinetic energy generated by the mean velocity gradient, and ε is the turbulent kinetic energy dissipation rate.

[0074] The formula form of the ε equation (turbulent kinetic energy dissipation rate equation) is as follows:

[0075]

[0076] Among them, σ ∈ is the turbulent Prandtl number of the ε equation, and C1 and C2 are empirical constants.

[0077] The formula of the NS equations is as follows:

[0078]

[0079] Where ρ is the density of the fluid, u is the velocity vector, p is the pressure, τ is the shear stress tensor, and f is the external force vector.

[0080] The equation of motion is given by:

[0081]

[0082] Where m is the mass of the particle, r is the position vector of the particle, and F is the force vector acting on the particle.

[0083] Specifically, according to the initial parameter information and the fluid parameter information, the NS equations are solved using a preset numerical method (e.g., finite volume method or finite element method, etc.) to obtain physical parameters such as velocity and pressure of the fluid phase; according to the selected turbulence model, the physical parameters of turbulence (such as turbulent kinetic energy k and turbulent kinetic energy dissipation rate ∈) are solved. Then, according to the fluid parameter information, the initial parameter information, the velocity and pressure of the fluid phase, and the physical parameters such as the physical parameters of turbulence, the motion equations of the particles in the inhaler are solved using a preset numerical integration method (e.g., Euler method, Runge-Kutta method, etc.), and the physical quantities such as displacement and velocity of the particles are obtained, thereby determining the motion trajectory of the particles. In the solution process of this embodiment, the interaction force between the particles and the fluid needs to be considered to improve the accuracy of the motion trajectory of the simulated particles.

[0084] In addition, it should be noted that the interaction between the continuous phase and the particle phase is achieved through phase coupling. This is usually achieved by adding source terms to the equations solved for the fluid phase and the particle phase.

[0085] The modeling and calculation work in this embodiment are all completed using the open source software Openfoam, and the source code program is secondary developed and modified to reduce the purchase cost of commercial software.

[0086] The embodiment of the present invention, through the above scheme, includes: constructing NS equations, turbulence model and motion equations of particles; according to the initial parameter information and the fluid parameter information, using the NS equations and the turbulence model to calculate physical parameters; according to the physical parameters, the initial parameter information and the fluid parameter information, using the motion equation to solve the motion trajectory of the particles in the inhaler. It is realized that the NS equations are established to solve the physical parameters of the fluid phase (continuous phase), and then the motion trajectory of each particle in the flow field is tracked and solved by Newton's second law to reflect the entire particle phase discreteness.

[0087] In one embodiment of the present invention, a geometric model of an inhaler is constructed and meshed, including:

[0088] According to the model size of the inhaler, a plurality of coordinate positions are determined; according to the plurality of coordinate positions, the physical space of the geometric model is determined; and the number of grids and the grid expansion rate are determined in the physical space to grid the geometric model.

[0089] Specifically, according to the design drawings or three-dimensional model of the inhaler, its key size parameters, such as length, width, height, etc., can be obtained. Select multiple key points on the inhaler model, which can be the outlet, inlet, vertex, edge, center point, etc. of geometric features. The coordinates of these points will be used to describe the geometry of the model. Select a suitable coordinate system (such as Cartesian coordinate system), and determine the coordinate position of each key point, thereby through multiple coordinate positions, determine the scope of the geometric model in the physical space. Further, determine the number of grids and grid expansion rate in the physical space, optionally, according to the complexity of the model and the computing resources, determine the number of grids in the entire physical space, usually the more the number of grids, the higher the calculation accuracy, but the computing cost is also higher; the grid expansion rate refers to the ratio of the grid gradually becoming larger or smaller in a certain direction, thereby realizing that the physical space of the geometric model is divided into multiple grid units, so that numerical simulation is carried out.

[0090] The embodiment of the present invention, through the above scheme, includes: determining multiple coordinate positions according to the model size of the inhaler; determining the physical space of the geometric model according to the multiple coordinate positions; determining the number of grids and the grid expansion rate in the physical space to grid the geometric model. By constructing a geometric model of the inhaler and gridding the geometric model, the motion trajectory of the particles in the inhaler can be simulated according to the geometric model.

[0091] The particle motion simulation device provided by the present invention is described below. The particle motion simulation device described below and the particle motion simulation method described above can be referenced to each other.

[0092] Figure 2 Schematic diagram of the structure of the particle motion simulation device provided by the present invention. Figure 2 As shown, a particle motion simulation device according to an embodiment of the present invention comprises:

[0093] A construction module 21 is used to construct a geometric model of the inhaler and mesh the geometric model;

[0094] A first determination module 22, used to determine initial parameter information of particles in the geometric model;

[0095] A second determination module 23, used to determine fluid parameter information of the fluid phase;

[0096] The solving module 24 is used to solve the motion trajectory of the particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information.

[0097] The building block 21 is also used to:

[0098] According to the model size of the inhaler, a plurality of coordinate positions are determined;

[0099] Determining a physical space of the geometric model according to the multiple coordinate positions;

[0100] The number of grids and the grid expansion rate are determined in the physical space to grid the geometric model.

[0101] The solution module 24 is also used for:

[0102] Construct the NS equations, turbulence models, and equations of motion for particles;

[0103] According to the initial parameter information and the fluid parameter information, physical parameters are calculated using the NS equations and the turbulence model;

[0104] The motion trajectory of the particles in the inhaler is obtained by solving the motion equation according to the physical parameters, the initial parameter information and the fluid parameter information.

[0105] The particle motion simulation device also includes:

[0106] Determine geometric boundary conditions corresponding to the geometric model, wherein the geometric boundary conditions include pressure values ​​corresponding to an inlet and an outlet in the geometric model, and the pressure values ​​corresponding to the inlet and the outlet are used as parameters for solving the motion trajectory of the particles in the inhaler.

[0107] The particle motion simulation device also includes:

[0108] The initial parameter information includes at least one of the starting emission position of the particles, the initial velocity field, the initial pressure field, the volume fraction boundary conditions of air and particles, the magnitude and direction of gravity, fluid properties, control parameters, parallel parameters and discrete formats; wherein the volume fraction of particles refers to the content of particles;

[0109] The fluid parameter information includes at least one of kinematic viscosity, fluid density, surface tension, and whether the fluid is a Newtonian fluid.

[0110] The particle motion simulation device also includes:

[0111] The starting emission position of the particles is determined according to the exit position of the capsule compartment in the geometric model.

[0112] It should be noted here that the above-mentioned device provided in the embodiment of the present invention can implement all the method steps implemented in the above-mentioned method embodiment, and can achieve the same technical effect. The parts and beneficial effects that are the same as the method embodiment in this embodiment will not be described in detail here.

[0113] Figure 3 is a schematic diagram of the structure of the electronic device provided by the present invention, such as Figure 3 As shown, the electronic device may include: a processor 310, a memory 320, a communications interface 330 and a communication bus 340, wherein the processor 310, the memory 320 and the communications interface 330 communicate with each other through the communication bus 340. The processor 310 may call the logic instructions in the memory 320 to execute the particle motion simulation method, which includes: constructing a geometric model of the inhaler and meshing the geometric model; determining the initial parameter information of the particles in the geometric model; determining the fluid parameter information of the fluid phase; and solving the motion trajectory of the particles in the inhaler using a pre-constructed motion simulation model according to the initial parameter information and the fluid parameter information.

[0114] In addition, the logic instructions in the above-mentioned memory 320 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0115] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which is implemented when the computer program is executed by a processor to execute the particle motion simulation method provided by the above methods.

[0116] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the particle motion simulation method provided by the above-mentioned methods.

[0117] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0118] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0119] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A particle motion simulation method, characterized in that: include: constructing a geometric model of the inhaler and meshing the geometric model; Determining initial parameter information of particles in the geometric model; determining fluid parameter information of the fluid phase; According to the initial parameter information and the fluid parameter information, a pre-constructed motion simulation model is used to solve the motion trajectory of the particles in the inhaler.

2. The particle motion simulation method according to claim 1, characterized in that: The step of constructing a geometric model of the inhaler and meshing the geometric model comprises: According to the model size of the inhaler, a plurality of coordinate positions are determined; Determining a physical space of the geometric model according to the multiple coordinate positions; The number of grids and the grid expansion rate are determined in the physical space to grid the geometric model.

3. The particle motion simulation method according to claim 1, characterized in that: The method of solving the motion trajectory of the particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information includes: Construct the NS equations, turbulence models, and equations of motion for particles; According to the initial parameter information and the fluid parameter information, physical parameters are calculated using the NS equations and the turbulence model; The motion trajectory of the particles in the inhaler is obtained by solving the motion equation according to the physical parameters, the initial parameter information and the fluid parameter information.

4. The particle motion simulation method according to claim 1, characterized in that: After constructing the geometric model of the inhaler, the method further comprises: Determine geometric boundary conditions corresponding to the geometric model, wherein the geometric boundary conditions include pressure values ​​corresponding to an inlet and an outlet in the geometric model, and the pressure values ​​corresponding to the inlet and the outlet are used as parameters for solving the motion trajectory of the particles in the inhaler.

5. The particle motion simulation method according to claim 1, characterized in that: The initial parameter information includes at least one of the starting emission position of the particles, the initial velocity field, the initial pressure field, the volume fraction boundary conditions of air and particles, the magnitude and direction of gravity, fluid properties, control parameters, parallel parameters and discrete formats; wherein the volume fraction of particles refers to the content of particles; The fluid parameter information includes at least one of kinematic viscosity, fluid density, surface tension, and whether the fluid is a Newtonian fluid.

6. The particle motion simulation method according to claim 5, characterized in that: The starting emission position of the particles is determined according to the exit position of the capsule compartment in the geometric model.

7. A particle motion simulation device, characterized in that: include: A construction module, used for constructing a geometric model of the inhaler and meshing the geometric model; A first determination module, used to determine initial parameter information of particles in the geometric model; A second determination module, used to determine fluid parameter information of the fluid phase; The solution module is used to solve the motion trajectory of the particles in the inhaler using a pre-built motion simulation model according to the initial parameter information and the fluid parameter information.

8. The particle motion simulation device according to claim 7, characterized in that: The solution module is also used for: Construct the NS equations, turbulence models, and equations of motion for particles; According to the initial parameter information and the fluid parameter information, physical parameters are calculated using the NS equations and the turbulence model; The motion trajectory of the particles in the inhaler is obtained by solving the motion equation according to the physical parameters, the initial parameter information and the fluid parameter information.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that: When the processor executes the program, the particle motion simulation method according to any one of claims 1 to 6 is implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the particle motion simulation method according to any one of claims 1 to 6 is implemented.

Images