CN113051842A - non-Newtonian fluid simulation method and device - Google Patents
non-Newtonian fluid simulation method and device Download PDFInfo
- Publication number
- CN113051842A CN113051842A CN202110246996.4A CN202110246996A CN113051842A CN 113051842 A CN113051842 A CN 113051842A CN 202110246996 A CN202110246996 A CN 202110246996A CN 113051842 A CN113051842 A CN 113051842A
- Authority
- CN
- China
- Prior art keywords
- phase
- denotes
- suspended particle
- newtonian fluid
- liquid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 167
- 238000004088 simulation Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000002245 particle Substances 0.000 claims abstract description 135
- 239000012071 phase Substances 0.000 claims abstract description 133
- 239000007791 liquid phase Substances 0.000 claims abstract description 67
- 239000000725 suspension Substances 0.000 claims abstract description 42
- 230000008859 change Effects 0.000 claims abstract description 35
- 230000008878 coupling Effects 0.000 claims abstract description 16
- 238000010168 coupling process Methods 0.000 claims abstract description 16
- 238000005859 coupling reaction Methods 0.000 claims abstract description 16
- 230000006870 function Effects 0.000 claims description 38
- 239000007788 liquid Substances 0.000 claims description 30
- 239000007787 solid Substances 0.000 claims description 14
- 239000000463 material Substances 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 10
- 239000007790 solid phase Substances 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 7
- 238000009792 diffusion process Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 7
- 230000008719 thickening Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000000518 rheometry Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000010008 shearing Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000903 blocking effect Effects 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000006071 cream Substances 0.000 description 2
- 230000003278 mimic effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 239000006260 foam Substances 0.000 description 1
- 235000008960 ketchup Nutrition 0.000 description 1
- 239000000314 lubricant Substances 0.000 description 1
- 239000011859 microparticle Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 235000015067 sauces Nutrition 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a non-Newtonian fluid simulation method and a non-Newtonian fluid simulation device, wherein the method comprises the following steps: respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation; coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspended body model with the property between the liquid phase and the suspended particle phase; constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model; and constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastic-plastic constitutive model for simulating the rheological property of the non-Newtonian fluid. The invention realizes the simulation of non-Newtonian property by modeling the real physical mechanism of the bottom layer, provides physical interpretable parameter control, thereby uniformly simulating different types of non-Newtonian fluids and accurately depicting the rheological property of the non-Newtonian fluids.
Description
Technical Field
The invention relates to the technical field of computer graphics fluid simulation and rendering, in particular to a non-Newtonian fluid simulation method and device.
Background
Non-newtonian fluids differ from general fluids in that their viscosity changes with shear flow and are widely used in everyday life and industrial applications, such as ketchup, cream and various lubricants. Non-newtonian fluids are also gaining increasing attention in the graphics arts due to their wide use and display of rich visual effects. The complex rheological properties of non-newtonian fluids are not fully understood compared to general newtonian fluids, and their highly nonlinear behavior also presents significant challenges to simulation.
It is now common practice in the industry to describe non-newtonian fluids using an exponential model, for example, the Herschel-Bulkley model,wherein τ is shear stress, τ0The yield stress, k is the consistency index and n is the fluidity index. When n is more than 1, the fluid is shear thickening fluid; when n is<1, shear-thinning fluid. Based on this unique image description, there have been many efforts in graphics to successfully simulate various non-Newtonian fluids, such as, for example, the method of codimensional simulation proposed by Zhu et al, which can simulate a non-Newtonian fluid that is spread into a thin layer, and the method of simulating substances such as shaving foam that have a viscosity that varies with shear rate proposed by Yue et al.
However, these prior art methods treat non-newtonian fluids as homogeneous single fluids, ignoring the internal microstructure. However, most non-Newtonian fluids are suspensions, i.e., a complex system of small particles suspended in Newtonian fluids, such as various sauces, creams, blood, etc., are suspensions that do not truly mimic the properties of non-Newtonian fluids. Therefore, there is a need for a non-Newtonian fluid simulation method and apparatus to solve the above problems.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a non-Newtonian fluid simulation method and a non-Newtonian fluid simulation device.
The invention provides a non-Newtonian fluid simulation method, which comprises the following steps:
respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation;
coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase;
constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model;
and constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastic-plastic constitutive model for simulating the rheological property of the non-Newtonian fluid.
According to the non-Newtonian fluid simulation method provided by the invention, the fluid resistance formula is as follows:
wherein k isdDenotes the resistance constant, usDenotes the suspended particle phase velocity ufDenotes the liquid phase velocity, d denotes the suspended particle diameter, eta denotes the fluid viscosity, phisRepresents the suspended particle phase volume fraction;
the fluid pressure equation is:
wherein p isfDenotes the liquid phase density, pf0Denotes the liquid phase static density and κ denotes the stiffness coefficient.
According to the non-Newtonian fluid simulation method provided by the invention, the elastic-plastic constitutive model based on the suspended particle phase is constructed according to the plastic shear stress and the volume change of the two-phase suspension model, and the method comprises the following steps:
from the plastic shear stress, a first yield condition is established:
τ≤(μ+β)p;
wherein τ represents shear stress, μ represents friction coefficient, IvDenotes the viscosity number, phimRepresents the maximum suspended particle phase volume fraction, beta represents the expansion angle, p represents the pressure,representing the plastic shear rate, eta the liquid viscosity coefficient, mu1Denotes the viscosity number IvCoefficient of friction, mu, to 02Denotes the viscosity number IvCoefficient of friction to ∞, I0Is a constant parameter;
from the plastic stretching, a second yield condition is established:
p+∈≥0;
wherein epsilon represents a preset viscosity threshold value;
from the plastic compression, a third yield condition is established:
wherein φ represents the suspended particle phase volume fraction, and g (φ) represents a coefficient that varies with φ;
and constructing and obtaining an elastic-plastic constitutive model based on the suspended particle phase according to the yield surfaces corresponding to the first yield condition, the second yield condition and the third yield condition respectively.
According to the non-Newtonian fluid simulation method provided by the invention, the formula of the friction function is as follows:
φm=φj+(φc-φj)Ψ;
wherein the content of the first and second substances,representing the friction function of the particle after a change in the skeletal structure, cfA constant indicating a control change rate, H indicates a hardening rate,representing the plastic shear rate, ΨmRepresenting the maximum value of the friction function, psi representing the friction function before the change of the particle skeleton structure, xi representing the buckling coefficient, tau representing the shear stress, tau*Denotes the repulsion constant between particles, cbAnd l is a constant, phi, controlling the buckling effectmDenotes the maximum suspended particle phase volume fraction, phi denotes the suspended particle phase volume fraction, phicIs indicative of phimLower bound of (phi)jIs indicative of phimThe upper bound of (c).
According to the invention, the non-Newtonian fluid simulation method further comprises the following steps:
and carrying out discrete simulation solution on the elastic-plastic constitutive model by a material point method to obtain the rheological property of the non-Newtonian fluid.
According to the non-newtonian fluid simulation method provided by the invention, the elastic-plastic constitutive model is solved through discrete simulation by a material point method, and the method further comprises the following steps:
simulating the permeation between the liquid and the dry particles based on a liquid drift velocity model, the liquid drift velocity model being:
wherein, cdDenotes the diffusion constant,. phifIndicating the liquid phase volume fraction.
According to the non-Newtonian fluid simulation method provided by the invention, the liquid phase and the suspended particle phase of the non-Newtonian fluid are respectively simulated based on the continuity equation and the momentum equation, and the method comprises the following steps:
simulating a suspended particle phase in the non-Newtonian fluid based on a solid continuity equation and a solid momentum equation, the solid continuity equation being:
the solid momentum equation is as follows:
wherein the content of the first and second substances,represents the solid-phase equivalent density, ρsDenotes the solid phase density, t denotes the time, usRepresenting the suspended particle phase velocity, g representing gravity, fdRepresenting the fluid resistance, σsDenotes shear stress,. phisRepresents the solid phase volume fraction;
simulating a liquid phase in the non-Newtonian fluid based on a liquid continuity equation and a liquid momentum equation, wherein the liquid continuity equation is as follows:
the liquid momentum equation is as follows:
wherein phi isfRepresenting the liquid phase volume fraction, pfDenotes the liquid phase density, ufWhich is indicative of the velocity of the liquid phase,denotes the liquid phase equivalent density, σfIndicating viscous stress.
The present invention also provides a non-Newtonian fluid simulation device, comprising:
the first simulation module is used for respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation;
the coupling module is used for coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase;
the model construction module is used for constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model;
and the second simulation module is used for constructing a friction function according to the particle skeletal structure change of the two-phase suspension model and adding the friction function into the elastic-plastic constitutive model so as to simulate the rheological property of the non-Newtonian fluid.
The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor executes the program to implement the steps of any one of the above-described non-newtonian fluid simulation methods.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the non-newtonian fluid simulation method as described in any one of the above.
The non-Newtonian fluid simulation method and device provided by the invention provide a non-Newtonian fluid simulation method with higher expansibility, realize the simulation of non-Newtonian properties by modeling the real physical mechanism of the bottom layer of the non-Newtonian fluid simulation method, and provide parameter control which can be interpreted physically, thereby uniformly simulating different types of non-Newtonian fluids and accurately describing the rheological property of the non-Newtonian fluids.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of a non-Newtonian fluid simulation method provided by the present invention;
FIG. 2 is a schematic illustration of a yielding surface provided by the present invention;
FIG. 3 is a schematic representation of three shear thickening rheology mechanisms provided by the present invention;
FIG. 4 is a schematic diagram of a non-Newtonian fluid simulation device in accordance with the present invention;
fig. 5 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Although most non-Newtonian fluids belong to suspension, the existing simulation method for non-Newtonian fluids ignores the internal microstructure and treats the non-Newtonian fluids as a uniform single fluid, so that the simulation results in that the non-Newtonian fluids can not completely describe the real characteristics. Based on the above, the invention uses two phases to simulate the non-Newtonian fluid, can simulate the non-uniform non-Newtonian fluid, and can even directly start from raw materials (namely water and powder) to simulate the formation of the non-Newtonian fluid, thereby greatly expanding the phenomenon which can be simulated. Particularly for a discontinuous type of shear thickening phenomenon, in materials such as the European Union Jack (ovleck), which may exhibit rigid-like properties when sheared, but may return to fluid when left undisturbed, this is closely related to the structural changes within.
Fig. 1 is a schematic flow chart of a non-newtonian fluid simulation method provided by the present invention, and as shown in fig. 1, the present invention provides a non-newtonian fluid simulation method, including:
According to the method, a liquid phase and a suspended particle phase of the non-Newtonian fluid are obtained through simulation by taking two independent continuous medium phases as targets through a continuity equation and a momentum equation corresponding to liquid and suspended particles in the non-Newtonian fluid, wherein the suspended particle phase is a micro particle suspended in the non-Newtonian fluid and is described through the continuity equation and the momentum equation. Specifically, continuity equations are used to update particle density, while kinetic equations, including gravity, viscous forces, pressure, and shear stress, are used to update particle velocity.
And 102, coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase.
In the invention, fluid resistance and fluid pressure are applied to the two phases to generate the effect of solid-liquid strong coupling, and the property between the solid and the liquid is displayed, so that a two-phase suspension model of the property between the liquid phase and the suspended particle phase is obtained.
And 103, constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model.
According to the method, 3 plastic yield conditions are designed according to the characteristics of the two-phase suspension model, so that the elastoplasticity constitutive model based on the suspended particle phase is constructed and obtained.
And 104, constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastoplasticity constitutive model to simulate the rheological property of the non-Newtonian fluid.
In the invention, when a liquid phase and a suspended particle phase of a non-Newtonian fluid are respectively simulated, firstly, model parameters suitable for a target non-Newtonian fluid are preset, wherein the model parameters comprise yield condition parameters, updating parameters of a friction function, suspended particle diameter, liquid viscosity coefficient and the like, and the friction function is used for describing the change of a microscopic particle bone structure in a shear flow, so that the phenomenon of discontinuous shear thickening is simulated; further, the invention adopts a material point method, the material to be simulated (namely a liquid phase and a suspended particle phase) is converted into discrete material point particles, fluid resistance and fluid pressure are applied to the two phases, the effect of solid-liquid strong coupling is generated, each particle has proper density, mass and volume fraction, and a grid with proper size and resolution is generated for a simulation area for the calculation of the material point method; then, calculating respective volume fractions of two phases at any position in space through an elastic-plastic constitutive model; after the stress tensor of the suspended particles is updated, projecting the updated stress tensor of the suspended particle phase particles to yield surfaces corresponding to 3 yield conditions in an elastic-plastic constitutive model, so as to generate the characteristics of the non-Newtonian fluid and simulate the motion of the non-Newtonian fluid; finally, the invention draws each time step information obtained in the simulation on a display terminal to obtain a simulation image of the non-Newtonian fluid, and the simulation image is exported to be a data file for subsequent processing.
The non-Newtonian fluid simulation method provided by the invention provides a non-Newtonian fluid simulation method with higher expansibility, the simulation of non-Newtonian properties is realized by modeling a real physical mechanism of the bottom layer of the non-Newtonian fluid simulation method, and parameter control which can be interpreted physically is provided, so that different types of non-Newtonian fluids are simulated uniformly, and the rheological property of the non-Newtonian fluid is accurately described.
On the basis of the above embodiment, the fluid resistance formula is:
wherein k isdDenotes the resistance constant, usDenotes the suspended particle phase velocity ufDenotes the liquid phase velocity, d denotes the suspended particle diameter, eta denotes the fluid viscosity, phisRepresents the suspended particle phase volume fraction;
the fluid pressure equation is:
wherein p isfDenotes the liquid phase density, pf0Denotes the liquid phase static density and κ denotes the stiffness coefficient.
In the present invention, the two phases obtained by simulation of the above embodiments are coupled by fluid resistance and fluid pressure, and in the fluid resistance formula, the resistance is proportional to the velocity difference of the two phases and inversely proportional to the square of the suspended particle diameter. In the case of small suspended particles, the resistance becomes very large, and the present invention calculates it by using the reverse euler method. The present invention couples two phases via resistance and pressure to simulate the complex rheology of non-newtonian fluids between fluids and solids.
On the basis of the above embodiment, the constructing an elastoplastic constitutive model based on suspended particle phases according to the plastic shear stress and the volume change of the two-phase suspension model includes:
from the plastic shear stress, a first yield condition is established:
τ≤(μ+β)p;
wherein τ represents shear stress; μ denotes the coefficient of friction, IvRepresenting viscosity number, coefficient of friction with viscosity number Iv(ii) a change; phi is amRepresents the maximum suspended particle phase volume fraction; beta represents the expansion angle, beta ═ ca(φ-φeq),caIs a constant number of times, and is,is the equilibrium volume fraction; p represents the pressure of the gas to be pressurized,representing the plastic shear rate, eta the liquid viscosity coefficient, mu1Denotes the viscosity number IvCoefficient of friction, mu, to 02Denotes the viscosity number IvCoefficient of friction to ∞, I0Is a constant parameter;
from the plastic stretching, a second yield condition is established:
p+∈≥0;
wherein epsilon represents a preset viscosity threshold value;
from the plastic compression, a third yield condition is established:
wherein φ represents the suspended particle phase volume fraction, and g (φ) represents a coefficient that varies with φ;
and constructing and obtaining an elastic-plastic constitutive model based on the suspended particle phase according to the yield surfaces corresponding to the first yield condition, the second yield condition and the third yield condition respectively. FIG. 2 is a schematic representation of the yield surface provided by the present invention, the yield surface corresponding to each of the three yield conditions (i.e., f)1、f2And f3) As can be seen with reference to fig. 2, in the first yield condition, the factor before pressure is a function of shear rateAnd an increase in the suspended particle volume fraction phi; the pressure is limited by the third yield condition, when phi is less than or equal to phimWhile, the same is true with shear rateIs increased with an increase in; at phi > phimThe pressure is no longer limited and can be increased to any value, which means that during shearing, an unlimited increase in pressure causes an unlimited increase in the shear stress τ, so that a blocking phenomenon occurs.
According to the invention, an elastoplastic constitutive model which is correspondingly obeyed by suspended particles is designed according to the plastic shear stress and the volume change of the non-Newtonian fluid, and in the model, the internal stress of the suspended particle phase is limited by the product of pressure and friction coefficient, and is simultaneously influenced by the shear rate, the volume fraction of the suspended particle phase and the like, thereby forming a relatively complete description of the physical properties of the suspended body.
On the basis of the above embodiment, the formula of the friction function is:
φm=φj+(φc-φj)Ψ;
wherein the content of the first and second substances,representing the friction function of the particle after a change in the skeletal structure, cfA constant indicating a control change rate, H indicates a hardening rate,represents the plastic shear rate; ΨmRepresenting the maximum value of the friction function, psi representing the friction function before the change of the particle bone structure, the friction function being defined in the interval 0, psim]To (c) to (d); ξ denotes the buckling coefficient, τ denotes the shear stress, τ*Denotes the repulsion constant between particles, cbAnd l is a constant, phi, controlling the buckling effectmDenotes the maximum suspended particle phase volume fraction, phi denotes the suspended particle phase volume fraction, phicIs indicative of phimLower bound of (phi)jIs indicative of phimThe upper bound of (c). In the present invention, the friction function directly affects the maximum allowable volume fraction φmAnd in turn, the rheological behavior of the material.
In the present invention, psi is usedmWhen Ψ increases, meaning that the friction between suspended particles increases, the maximum volume fraction allowed decreases, going towards φcClose, and phi does not change. When phi ismDecreasing to near φ, a large increase in viscosity occurs, causing discontinuous shear thickening. FIG. 3 is a schematic diagram of three shear thickening rheology mechanisms provided by the present invention, as shown in FIG. 3, there are 3 rheology zones when φ<φcAny shear rate is permissible; when phi isc<φ<φjWhen the shear rate increases, the friction function increases, resulting in phimIf is less than phi, the blockage will be caused, therefore the shearing rate must be less than a certain threshold value; when phi is more than phijWhen the fluid is blocked, any shear flow is unlikely to occur.
The invention describes the proportion of the suspended particle phase and the liquid phase in the space by volume fraction, the volume fraction has key effect on the characteristics of the suspension, wherein, the volume fraction of the suspended particle phase is marked as phi, which describes the density degree of the space occupied by the suspended particle phase, and the volume fraction of the liquid phase at the same position is 1-phi according to the definition. In the present invention, the equivalent viscosity of the suspension increases with increasing φ and approaches φmWhen phi is larger than or equal to phi, the temperature is greatly increasedmA blocking phenomenon occurs and the material cannot flow in shear.
On the basis of the above embodiment, the method further includes:
and carrying out discrete simulation solution on the elastic-plastic constitutive model by a material point method to obtain the rheological property of the non-Newtonian fluid.
In the invention, a material dot method is adopted for simulation, and the specific steps are as follows:
step 201, interpolating the mass and the speed of the particles on a grid, and calculating the pressure gradient and the internal stress;
step 202, updating grid speed and applying boundary conditions;
step 203, interpolating the updated grid speed back to the particle, updating the position of the particle, and calculating the speed gradient of the position of the particle;
step 204, updating the stress tensor carried by the suspension phase particles, and calculating the pressure of the liquid phase particles;
step 205, obtaining the shear rate according to the elastic-plastic yield condition of the suspension phaseShear stress tau and pressure p, from which the friction function psi and the maximum volume fraction phi of the suspended phase particles are updatedm;
And step 206, repeating the steps 201 to 205 until the algorithm is ended.
On the basis of the above embodiment, the discrete simulation solution is performed on the elastic-plastic constitutive model by a material point method, and the method further includes:
simulating the permeation between the liquid and the dry particles based on a liquid drift velocity model, the liquid drift velocity model being:
wherein, cdDenotes the diffusion constant,. phifIndicating the liquid phase volume fraction.
In the present invention, in order to simulate the penetration between the liquid and the dry particles, the present application introduces diffusion based on a drift velocity model when interpolating the updated mesh velocity back to the particles (i.e., step 203 in the above-described embodiment), thereby simulating the phenomenon of penetration diffusion between two phases.
On the basis of the above embodiment, the simulating the liquid phase and the suspended particle phase of the non-newtonian fluid based on the continuity equation and the momentum equation respectively includes:
simulating a suspended particle phase in the non-Newtonian fluid based on a solid continuity equation and a solid momentum equation, the solid continuity equation being:
the solid momentum equation is as follows:
wherein the content of the first and second substances,represents the solid-phase equivalent density, ρsDenotes the solid phase density, t denotes the time, usRepresenting the suspended particle phase velocity, g representing gravity, fdRepresenting the fluid resistance, σsDenotes shear stress,. phisRepresents the solid phase volume fraction;
simulating a liquid phase in the non-Newtonian fluid based on a liquid continuity equation and a liquid momentum equation, wherein the liquid continuity equation is as follows:
the liquid momentum equation is as follows:
wherein phi isfRepresenting the liquid phase volume fraction, pfDenotes the liquid phase density, ufWhich is indicative of the velocity of the liquid phase,denotes the liquid phase equivalent density, σfIndicating viscous stress. In the invention, the pressure of the fluid is fixed as follows from the momentum equation
The invention constructs a unified framework for simulating non-Newtonian fluids, and simulates various non-Newtonian fluids through physically interpretable parameters. When simulating a suspension of arbitrary volume fraction, the coupling between them can be easily simulated. The yield criterion and the friction function adopted by the invention provide an accurate physical tool and can simulate the discontinuous shearing thickening phenomenon. Finally, the present invention encompasses an osmotic diffusion mechanism that can mimic the mixing process between separate phases in their initial state, thereby mimicking the formation of non-Newtonian fluids.
Fig. 4 is a schematic structural diagram of a non-newtonian fluid simulation apparatus provided by the present invention, and as shown in fig. 4, the present invention provides a non-newtonian fluid simulation apparatus, which includes a first simulation module 401, a coupling module 402, a model building module 403, and a second simulation module 404, where the first simulation module 401 is configured to respectively simulate a liquid phase and a suspended particle phase of a non-newtonian fluid based on a continuity equation and a momentum equation; the coupling module 402 is configured to couple the liquid phase and the suspended particle phase according to a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of a property between the liquid phase and the suspended particle phase; the model construction module 403 is configured to construct an elasto-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model; the second simulation module 404 is configured to construct a friction function according to the particle skeletal structure variation of the two-phase suspension model, and add the friction function to the elasto-plastic constitutive model for simulating the rheological properties of the non-newtonian fluid.
The non-Newtonian fluid simulation device provided by the invention provides a non-Newtonian fluid simulation method with higher expansibility, the simulation of non-Newtonian properties is realized by modeling a real physical mechanism of the bottom layer of the non-Newtonian fluid simulation device, and parameter control which can be interpreted physically is provided, so that different types of non-Newtonian fluids are simulated uniformly, and the rheological property of the non-Newtonian fluids is accurately described.
The apparatus provided by the present invention is used for executing the above method embodiments, and for details and flow, reference is made to the above embodiments, which are not described herein again.
Fig. 5 is a schematic structural diagram of an electronic device provided in the present invention, and as shown in fig. 5, the electronic device may include: a processor (processor)501, a communication interface (communication interface)502, a memory (memory)503 and a communication bus 504, wherein the processor 501, the communication interface 502 and the memory 503 are communicated with each other through the communication bus 504. Processor 501 may invoke logic instructions in memory 503 to perform a non-newtonian fluid simulation method, the method comprising: respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation; coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase; constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model; and constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastic-plastic constitutive model for simulating the rheological property of the non-Newtonian fluid.
In addition, the logic instructions in the memory 503 may be implemented in the form of software functional units and stored in a computer readable storage medium when the logic instructions are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-only memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
In another aspect, the present invention also provides a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, enable the computer to perform the non-newtonian fluid simulation method provided by the above methods, the method comprising: respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation; coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase; constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model; and constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastic-plastic constitutive model for simulating the rheological property of the non-Newtonian fluid.
In yet another aspect, the present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program, which when executed by a processor is implemented to perform the non-newtonian fluid simulation method provided in the above embodiments, the method comprising: respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation; coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase; constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model; and constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastic-plastic constitutive model for simulating the rheological property of the non-Newtonian fluid.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method of simulating a non-newtonian fluid, comprising:
respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation;
coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase;
constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model;
and constructing a friction function according to the particle skeletal structure change of the two-phase suspension model, and adding the friction function into the elastic-plastic constitutive model for simulating the rheological property of the non-Newtonian fluid.
2. A non-newtonian fluid simulation method according to claim 1, wherein the fluid resistance formula is:
wherein k isdDenotes the resistance constant, usDenotes the suspended particle phase velocity ufDenotes the liquid phase velocity, d denotes the suspended particle diameter, eta denotes the fluid viscosity, phisRepresents the suspended particle phase volume fraction;
the fluid pressure equation is:
wherein p isfDenotes the liquid phase density, pf0Denotes the liquid phase static density and κ denotes the stiffness coefficient.
3. The non-Newtonian fluid simulation method of claim 1, wherein the constructing an elasto-plastic constitutive model based on suspended particle phase according to the plastic shear stress and volume change of the two-phase suspension model comprises:
from the plastic shear stress, a first yield condition is established:
τ≤(μ+β)p;
wherein τ represents shear stress, μ represents friction coefficient, IvDenotes the viscosity number, phimRepresents the maximum suspended particle phase volume fraction, beta represents the expansion angle, p represents the pressure,representing the plastic shear rate, eta the liquid viscosity coefficient, mu1Denotes the viscosity number IvCoefficient of friction, mu, to 02Denotes the viscosity number IvCoefficient of friction to ∞, I0Is a constant parameter;
from the plastic stretching, a second yield condition is established:
p+∈≥0;
wherein epsilon represents a preset viscosity threshold value;
from the plastic compression, a third yield condition is established:
wherein φ represents the suspended particle phase volume fraction, and g (φ) represents a coefficient that varies with φ;
and constructing and obtaining an elastic-plastic constitutive model based on the suspended particle phase according to the yield surfaces corresponding to the first yield condition, the second yield condition and the third yield condition respectively.
4. A non-newtonian fluid simulation method according to claim 1, wherein the friction function is formulated as:
φm=φj+(φc-φj)Ψ;
wherein the content of the first and second substances,representing the friction function of the particle after a change in the skeletal structure, cfA constant indicating a control change rate, H indicates a hardening rate,representing the plastic shear rate, ΨmRepresenting the maximum value of the friction function, psi representing the friction function before the change of the particle skeleton structure, xi representing the buckling coefficient, tau representing the shear stress, tau*Denotes the repulsion constant between particles, cbAnd l is a constant, phi, controlling the buckling effectmDenotes the maximum suspended particle phase volume fraction, phi denotes the suspended particle phase volume fraction, phicIs indicative of phimLower bound of (phi)jIs indicative of phimThe upper bound of (c).
5. A non-newtonian fluid simulation method according to claim 1, further comprising:
and carrying out discrete simulation solution on the elastic-plastic constitutive model by a material point method to obtain the rheological property of the non-Newtonian fluid.
6. A non-Newtonian fluid simulation method according to claim 5, wherein the elasto-plastic constitutive model is solved by discrete simulation by a particle method, the method further comprising:
simulating the permeation between the liquid and the dry particles based on a liquid drift velocity model, the liquid drift velocity model being:
wherein, cdDenotes the diffusion constant,. phifIndicating the liquid phase volume fraction.
7. A method for modeling a non-newtonian fluid according to claim 1, wherein the modeling the liquid phase and the suspended particle phase of the non-newtonian fluid based on the continuity equation and the momentum equation, respectively, comprises:
simulating a suspended particle phase in the non-Newtonian fluid based on a solid continuity equation and a solid momentum equation, the solid continuity equation being:
the solid momentum equation is as follows:
wherein the content of the first and second substances,represents the solid-phase equivalent density, ρsDenotes the solid phase density, t denotes the time, usRepresenting the suspended particle phase velocity, g representing gravity, fdRepresenting the fluid resistance, σsDenotes shear stress,. phisRepresents the solid phase volume fraction;
simulating a liquid phase in the non-Newtonian fluid based on a liquid continuity equation and a liquid momentum equation, wherein the liquid continuity equation is as follows:
the liquid momentum equation is as follows:
8. A non-newtonian fluid simulation device, comprising:
the first simulation module is used for respectively simulating a liquid phase and a suspended particle phase of the non-Newtonian fluid based on a continuity equation and a momentum equation;
the coupling module is used for coupling the liquid phase and the suspended particle phase through a fluid resistance formula and a fluid pressure formula to obtain a two-phase suspension model of the property between the liquid phase and the suspended particle phase;
the model construction module is used for constructing an elastic-plastic constitutive model based on the suspended particle phase according to the plastic shear stress and the volume change of the two-phase suspension model;
and the second simulation module is used for constructing a friction function according to the particle skeletal structure change of the two-phase suspension model and adding the friction function into the elastic-plastic constitutive model so as to simulate the rheological property of the non-Newtonian fluid.
9. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor when executing the computer program implements the steps of the non-newtonian fluid simulation method of any of claims 1-7.
10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements the steps of the non-newtonian fluid simulation method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110246996.4A CN113051842B (en) | 2021-03-05 | 2021-03-05 | Non-Newtonian fluid simulation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110246996.4A CN113051842B (en) | 2021-03-05 | 2021-03-05 | Non-Newtonian fluid simulation method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113051842A true CN113051842A (en) | 2021-06-29 |
CN113051842B CN113051842B (en) | 2024-06-04 |
Family
ID=76510650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110246996.4A Active CN113051842B (en) | 2021-03-05 | 2021-03-05 | Non-Newtonian fluid simulation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113051842B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114510858A (en) * | 2022-01-06 | 2022-05-17 | 河北乐凯化工工程设计有限公司 | Flow field determination method, equipment and medium based on non-Newtonian fluid |
CN115270651A (en) * | 2022-06-20 | 2022-11-01 | 北京科技大学 | Monocular video oriented non-Newtonian fluid simulation reconstruction method |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4637250A (en) * | 1985-01-25 | 1987-01-20 | State University Of New York | Apparatus and method for viscosity measurements for Newtonian and non-Newtonian fluids |
AU2001274674A1 (en) * | 2000-06-30 | 2002-01-08 | Cicero Dental Systems B.V. | Dental prosthesis and process for the making thereof |
CN1450345A (en) * | 2003-03-27 | 2003-10-22 | 北京科技大学 | Instrument for determining rheological properties of Non-Newtonian viscosity fluid |
WO2010140957A1 (en) * | 2009-06-03 | 2010-12-09 | Algoryx Simulation Ab | A method, an apparatus and computer program product for simulating dynamic fluids |
US20110153299A1 (en) * | 2009-12-18 | 2011-06-23 | Young Hee Kim | Method and apparatus for simulating viscoelastic fluid in smoothed particle hydrodynamics based fluid simulation |
US20130332089A1 (en) * | 2012-06-10 | 2013-12-12 | Halliburton Energy Services, Inc. | Methods for Predicting Dynamic Sag Using Viscometer/Rheometer Data |
JP2019125102A (en) * | 2018-01-15 | 2019-07-25 | 株式会社豊田中央研究所 | Fluid analysis device, fluid analysis method, and fluid analysis program |
CN110532727A (en) * | 2019-09-09 | 2019-12-03 | 扬州大学 | It can be used for the method for numerical simulation of common non-newtonian fluid |
CN112100943A (en) * | 2020-09-21 | 2020-12-18 | 湖南御家化妆品制造有限公司 | Method for predicting mixing performance of high-viscosity fluid |
-
2021
- 2021-03-05 CN CN202110246996.4A patent/CN113051842B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4637250A (en) * | 1985-01-25 | 1987-01-20 | State University Of New York | Apparatus and method for viscosity measurements for Newtonian and non-Newtonian fluids |
AU2001274674A1 (en) * | 2000-06-30 | 2002-01-08 | Cicero Dental Systems B.V. | Dental prosthesis and process for the making thereof |
CN1450345A (en) * | 2003-03-27 | 2003-10-22 | 北京科技大学 | Instrument for determining rheological properties of Non-Newtonian viscosity fluid |
WO2010140957A1 (en) * | 2009-06-03 | 2010-12-09 | Algoryx Simulation Ab | A method, an apparatus and computer program product for simulating dynamic fluids |
US20120123754A1 (en) * | 2009-06-03 | 2012-05-17 | Algoryx Simulation Ab | Method, an apparatus and computer program product for simulating dynamic fluids |
US20110153299A1 (en) * | 2009-12-18 | 2011-06-23 | Young Hee Kim | Method and apparatus for simulating viscoelastic fluid in smoothed particle hydrodynamics based fluid simulation |
US20130332089A1 (en) * | 2012-06-10 | 2013-12-12 | Halliburton Energy Services, Inc. | Methods for Predicting Dynamic Sag Using Viscometer/Rheometer Data |
JP2019125102A (en) * | 2018-01-15 | 2019-07-25 | 株式会社豊田中央研究所 | Fluid analysis device, fluid analysis method, and fluid analysis program |
CN110532727A (en) * | 2019-09-09 | 2019-12-03 | 扬州大学 | It can be used for the method for numerical simulation of common non-newtonian fluid |
CN112100943A (en) * | 2020-09-21 | 2020-12-18 | 湖南御家化妆品制造有限公司 | Method for predicting mixing performance of high-viscosity fluid |
Non-Patent Citations (1)
Title |
---|
孙洪全;韩纪庆;: "基于粒子的非牛顿流体模拟统一模型", 哈尔滨工业大学学报, no. 07, 15 July 2010 (2010-07-15), pages 101 - 104 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114510858A (en) * | 2022-01-06 | 2022-05-17 | 河北乐凯化工工程设计有限公司 | Flow field determination method, equipment and medium based on non-Newtonian fluid |
CN114510858B (en) * | 2022-01-06 | 2023-06-30 | 河北乐凯化工工程设计有限公司 | Flow field determining method, device and medium based on non-Newtonian fluid |
CN115270651A (en) * | 2022-06-20 | 2022-11-01 | 北京科技大学 | Monocular video oriented non-Newtonian fluid simulation reconstruction method |
CN115270651B (en) * | 2022-06-20 | 2024-03-15 | 北京科技大学 | Monocular video-oriented non-Newtonian fluid simulation reconstruction method |
Also Published As
Publication number | Publication date |
---|---|
CN113051842B (en) | 2024-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fraggedakis et al. | Yielding the yield-stress analysis: a study focused on the effects of elasticity on the settling of a single spherical particle in simple yield-stress fluids | |
McKee et al. | The MAC method | |
Stroeven et al. | Numerical determination of representative volumes for granular materials | |
Xie et al. | Lattice Boltzmann modeling for multiphase viscoplastic fluid flow | |
Vananroye et al. | Effect of confinement on the steady-state behavior of single droplets during shear flow | |
Flenner et al. | Kinetic Monte Carlo and cellular particle dynamics simulations of multicellular systems | |
Wallevik et al. | Analysis of shear rate inside a concrete truck mixer | |
CN113051842A (en) | non-Newtonian fluid simulation method and device | |
Pan et al. | Dissipative particle dynamics simulation of droplet suspension in shear flow at low Capillary number | |
Zeidan | Assessment of mixture two-phase flow equations for volcanic flows using Godunov-type methods | |
CN111475937B (en) | Simulation method for flow-solid two-phase flow fluidization landslide | |
Alexandrou et al. | Flow instabilities of Herschel–Bulkley fluids | |
Kefayati et al. | Lattice Boltzmann method for the simulation of the steady flow of a Bingham fluid in a pipe of square cross-section | |
Richefeu et al. | Lattice Boltzmann modelling of liquid distribution in unsaturated granular media | |
Phan-Thien et al. | A spring model for suspended particles in dissipative particle dynamics | |
Moschopoulos et al. | The concept of elasto-visco-plasticity and its application to a bubble rising in yield stress fluids | |
Pan et al. | Numerical study of spheres settling in Oldroyd-B fluids | |
Verhulst et al. | Influence of viscoelasticity on drop deformation and orientation in shear flow. Part 2: Dynamics | |
Le-Cao et al. | A dissipative particle dynamics model for thixotropic materials exhibiting pseudo-yield stress behaviour | |
Dzwinel et al. | A two-level, discrete particle approach for large-scale simulation of colloidal aggregates | |
Gac et al. | A two-dimensional modeling of binary coalescence time using the two-color lattice-Boltzmann method | |
Vahabi | The effect of thixotropy on deformation of a single droplet under simple shear flow | |
Derksen | Solid particle mobility in agitated Bingham liquids | |
Hooper et al. | Assessing a flow-based finite element model for the sintering of viscoelastic particles | |
Bašić | An Overview of Numerical Methods for Non-Newtonian Flows |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |