Mixing tools such as mixers or blenders for adhesives and sealants are widely used. Mixers are commercially available for two-part adhesives. A static or dynamic mixer is a precision engineered device for the continuous mixing of fluid materials, without or with moving components. Normally the fluids to be mixed are liquid, but mixers can also be used to mix gas streams, disperse gas into liquid or blend immiscible liquids. The energy needed for mixing comes from a loss in pressure as fluids flow through the static mixer. One design of mixer is a plate-type mixer and another common device type consists of mixer elements (e.g., a rotating screw) contained in a cylindrical (tube) or squared housing.
There is a desire to optimize mixing tools to improve the mixing performance and quality, as well as minimizing waste. The present disclosure provides systems and methods for optimizing and designing mixing tools using modeling and visualization.
In one aspect, the present disclosure describes a method of optimizing a mixing tool to mix a plurality of fluid materials to obtain a mixture. The method includes representing a first geometry of the mixing tool with a first digital three-dimensional (3D) model; providing the first digital 3D model to a computational fluid dynamics (CFD) simulation model implemented by a processor to simulate a mixing process of the plurality of fluids to generate a particle density distribution of the mixture inside the mixing tool; converting the particle density distribution of the mixture to a first spatial distribution of fluid concentration of the mixture; measuring a second spatial distribution of fluid concentration of the mixture inside the mixing tool when mixing the plurality of fluids using the mixing tool; and comparing the first and second spatial distributions of fluid concentration to determine whether the first and second spatial distributions match with each other.
In another aspect, the present disclosure describes a computer-implemented method to design and optimize a tool. The method includes representing a first geometry of the tool with a first digital 3D model; providing the first digital 3D model to a computational fluid dynamics (CFD) simulation model implemented by a processor to simulate a fluid flow inside the tool to generate a particle density distribution of one or more fluids inside the tool; converting the particle density distribution to a spatial distribution of fluid concentration of the mixture; and visualizing the spatial distribution of fluid concentration and a measured spatial distribution of fluid concentration in a graphic user interface (GUI)
In another aspect, the present disclosure describes a computer-implemented system to design and optimize a tool. The system includes a module to represent a first geometry of the tool with a first digital 3D model; a module to provide the first digital 3D model to a computational fluid dynamics (CFD) simulation model implemented by a processor to simulate a fluid flow inside the tool to generate a particle density distribution of one or more fluids inside the tool; a module to convert the particle density distribution to a spatial distribution of fluid concentration of the mixture; and a module to visualize the spatial distribution of fluid concentration and a measured spatial distribution of fluid concentration in a graphic user interface (GUI).
Various unexpected results and advantages are obtained in exemplary embodiments of the disclosure. One such advantage of exemplary embodiments of the present disclosure is that the systems and methods can carry out multiple comparisons for tool design and optimization through high-definition visualizations of simulation data.
The performance of adhesives depends on the mixing quality in the dispensing process of multiple parts. The mixer tip design is critical for achieving a targeted mixing performance as well as minimizing both waste and the pressure drop during dispensing. The methods and systems described herein can also apply to other mixing systems with relatively high viscosity such as sealant mixing, dental adhesive, etc.
Various aspects and advantages of exemplary embodiments of the disclosure have been summarized. The above Summary is not intended to describe each illustrated embodiment or every implementation of the present certain exemplary embodiments of the present disclosure. The Drawings and the Detailed Description that follow more particularly exemplify certain preferred embodiments using the principles disclosed herein.
The disclosure may be more completely understood in consideration of the following detailed description of various embodiments of the disclosure in connection with the accompanying figures, in which:
In the drawings, like reference numerals indicate like elements. While the above-identified drawing, which may not be drawn to scale, sets forth various embodiments of the present disclosure, other embodiments are also contemplated, as noted in the Detailed Description. In all cases, this disclosure describes the presently disclosed disclosure by way of representation of exemplary embodiments and not by express limitations. It should be understood that numerous other modifications and embodiments can be devised by those skilled in the art, which fall within the scope and spirit of this disclosure.
Systems and methods for designing and optimizing tools are provided. A computational fluid dynamics (CFD) simulation model implemented by a processor is provided to simulate a fluid flow inside the tool to generate a particle density distribution of one or more fluids being mixed inside the tool. The particle density distribution is converted to a spatial distribution of fluid concentration of the mixture, and visualized in a graphic user interface (GUI) to allow for high resolution visualization of any spatially resolved discontinuous variable of fluid flow in the tool.
At 110, a computation fluid dynamics (CFD) simulation model can be implemented by a computing device such as, for example, a processor, to simulate a mixing process of multiple fluids inside any given mixing tool. A geometry of the mixing tool can be represented by a digital three-dimensional (3D) model. The digital 3D model can be in the form of an electronic file for computer-aided design (CAD), computer-aided manufacturing (CAM), computer-aided engineering (CAE), or other suitable applications. One exemplary digital 3D model of a mixing tool is visualized in a graphic user interface shown in
The digital 3D model of the mixing tool can be provided, along with other simulation model parameters, to the computer-implemented simulation method to simulate the process of mixing fluid materials to form a mixture by the mixing tool. The simulation results can be further processed to generate a spatial distribution of fluid concentration of the mixture. A spatial distribution of fluid concentration described herein refers to the quantity of a given type of fluid occupying a unit volume. The simulation model parameters for the CFD simulation model may include, for example, mesh density, time step, etc.
At 120, a physical mixing tool can be provided having the same geometry as that for the CFD simulation model at 110. Fluid materials can be mixed by the mixing tool under the same operation conditions for the simulation at 110 including, for example, the same rotational speed (RPM) of an internal rotating screw, the same fluid flow rates, etc. The spatial distribution of fluid concentration of the mixture can be physically measured for the mixture inside the mixing tool. One exemplary method is to use particle laden experimental CT-scan to measure fluid concentrations. The method 100 then proceeds to 130.
At 130, the simulated spatial distribution of fluid concentration from 110 and the measured spatial distribution of fluid concentration from 120 can be compared to determine whether the two spatial distributions match with each other at 135. When the two spatial distributions match with each other, the CFD simulation model at 110 can be validated and the method 100 proceeds to 160. When the two spatial distributions do not match with each other, the computation fluid dynamics (CFD) simulation model is modified and the method 100 proceeds back to 110. For example, the simulated and measured spatial distributions can be respectively visualized, e.g., being represented by iso-surfaces which can be overlaid in a graphic user interface (GUI). When the respective shapes and distributions of the iso-surfaces match with each other, the CFD simulation model at 110 can be validated. When the respective shapes and distributions of the iso-surfaces do not match with each other, the CFD simulation model can be modified at 140 and the method 100 proceeds back to 110 until the simulated and measured spatial distributions match with each other and the CFD simulation model is validated. Then the method 100 proceeds to 160.
At 160, a tool can be designed, and its geometry can be represented by a digital 3D model. Then, at 170, the validated CFD simulation model can be applied to the designed tool, and the process of mixing fluid materials to form a mixture can be simulated via a computing device. The simulation results can be further processed to generate data related to a fluid concentration distribution of the mixture, which can be visualized at 180. By reviewing the visualized data, a user can determine at 185 whether the designed tool meet the target performance, e.g., whether the fluid materials are uniformly mixed by the designed cool according to the simulation results. When the target performance is reached, the design is optimized and the method 100 proceeds to 190. When the target performance is not reached, the method 100 proceeds back to 160 to keep optimizing the tool design. At 190, a physical tool according to the optimized design is made, and its physical performance can be measured by, e.g., standard peel test, shear test, etc., to physically validate the tool.
The processor 12 can be included in any computing device. The processor 12 may include, for example, one or more general-purpose microprocessors, specially designed processors, application specific integrated circuits (ASIC), field programmable gate arrays (FPGA), a collection of discrete logic, and/or any type of processing device capable of executing the techniques described herein. In some embodiments, the processor (or any other processors described herein) may be described as a computing device. In some embodiments, the memory may be configured to store program instructions (e.g., software instructions) that are executed by the processor to carry out the processes or methods described herein. In other embodiments, the processes or methods described herein may be executed by specifically programmed circuitry of the processor. The processor may thus be configured to execute the techniques for any simulation and analysis described herein. The processor (or any other processors described herein) may include one or more processors.
At 210, a tool design with a given geometry is provided. The tool geometry can be represented by a digital 3D model. The tool geometry can be provided to a computation fluid dynamics (CFD) simulation model at 220 to simulate the process of mixing multiple fluid materials or species to generate a particle density distribution of the mixture at 230.
The CFD model described herein utilizes a particle mixing model that tracks weight-less particles in different components of the mixture of multiple fluids. Input for the particle mixing model includes the spatial coordinates of discrete tracer particles (or fluid representing particles) of the respective fluid species (e.g., a mixture including fluid species or parts A and B). In other words, a space distribution of a finite number of fluid representing particles (or tracer particles) is utilized to represent the continuous fluids to be mixed.
Numerical diffusion can be avoided in the particle mixing model as compared to a scalar mixing model. A typical scalar mixing model was described in, for example, Alméras, E., Plais, C., Euzenat, F., Risso, F., Roig, V. and Augier, F., 2016, Scalar mixing in bubbly flows: Experimental investigation and diffusivity modelling, Chemical Engineering Science, 140, pp. 114-122. The scalar mixing model solves passive scalar transport equation(s) for each type of fluid to obtain the concentration of each fluid. One problem accompanying the scalar mixing model is that the numerical method to solve such equation(s) will introduce error called numerical diffusion which results in overestimation of mixing quality compared to the reality. The particle mixing model applied herein can overcome the challenges by tracking the coordinates of fluid representing particles according to the fluid velocity field. In general, the particle mixing model is more accurate than the scalar mixing model by avoiding the problem of numerical diffusion.
A particle density distribution described herein refers to the distributions or locations of each fluid representing particle inside the mixing tool at a given time. The number of fluid representing particles depends on the mesh density (number of cells) of the CFD model. Usually the number of fluid representing particles is one order of magnitude lower than the number of cells used for the CFD model. The typical number of cells can be, for example, about 2 million or more in a 3D configuration of a tool. The computational time increases exponentially with the number of fluid representing particles used for a simulation. Thus, the number of fluid representing particles is limited by computational resources.
Referring again to
One exemplary conversion method uses a spatial coarse-graining of the particle density field. In other words, the method takes a set of points and provides them with some volume to better approximate the fluid elements they are supposed to represent. This is achieved through convolution of a particle's position with a reasonably chosen coarse graining function. The local concentration ca(r) can then be calculated through Equation (1) below:
c
α(r)=Σi=1N
Where ri is the position of the i-th particle of specie a and Na is the number of such particles. One embodiment employs a Gaussian as the coarse-graining function as shown in Equation (2) below:
In this case, a is a coarse-graining length scale, which sets how much volume the particle is intended to occupy. In addition, to enable more efficient computation, the Gaussian is truncated and shifted to zero at a specified cut-off rc and then normalized.
While a Gaussian coarse-graining function is used in Equation (2) to weight the contribution from nearby particles while converting a particle density distribution of the mixture to a spatial distribution of fluid concentration of the mixture, it is to be understood that any suitable decay function can be used to weight the contribution from nearby particles, including, for example, linear functions, polynomial functions, etc.
Referring again to
In one example, copper metal flakes with particle size 325 Mesh (44 microns and finer) was premixed with one part of the adhesive with volume fraction of less than 2%. Then, the mixing experiment of the two-part adhesive using different mixer tip was carried out. After the experiment, the filled mixer tip is CT-scanned. Due to the high density of the copper comparing to the adhesive material, the part with the copper particle can show higher density in the CT-scan. By analyzing the density intensity of the CT-scan image, the detailed mixing process can be obtained. For well mixed regions, the density intensity is similar throughout the space. For example, at the exit of the mixing tool where the species are well mixed, a substantially uniform density intensity is shown. For poorly mixed regions, the density intensity has a large spatial variation. For example, near the inlet of the mixing tool, a substantially non-uniform density intensity can be shown.
Referring again to
When the two spatial distributions match with each other, the CFD simulation model at 220 can be validated at 290. When the two spatial distributions do not match with each other, the computation fluid dynamics (CFD) simulation model is modified at 280 and the method 200 proceeds back to 220. Modifying a CFD simulation model may include, for example, varying simulation parameters including mesh density, time steps, etc.
The obtained local fluid concentration ca(r) at 240 can be further processed to calculate a spatially resolved mixing index MP(r) that can easily be visualized to assess mixing performance. The local mixing quality depends on the respective fluid concentration distributions of the components (species) of the mixture. For example, for a binary mixture with a concentration ratio between the two components or species of N1:N2, a mixing index can be defined by Equation (3):
where f1t is the target well mixing fraction of specie 1. The local weight fraction of specie 1 is defined as f1=c1(c2+c1), where c1 and c2 being the respective concentrations of species 1 and 2. The value of mixing index Mp varies between 0 and 1, where the value of 1 refers to a state of mixture being mixed to target proportions and the value of 0 refers to pure species.
Computer-implemented computational methods and systems are developed herein to determine the mixing quality as represented by, e.g., the mixing index MP(r), at desired locations within a mixing tool based on the obtained local fluid concentration ca(r). It is to be understood that the mixing index MP(r) and the local fluid concentration ca(r) can be derived from each other for one or more specifies of the mixture. The spatial distribution of mixing index can be visualized in the graphic user interface (GUI), along with a spatial distribution of fluid concentration.
A static mixer was used to validate a computational fluid dynamics (CFD) simulation model. The geometry of the static mixer was created via a CAD file. The CAD file was input to the CFD model for simulation of mixing fluid A and fluid B to form a mixture. Fluid A is 3M Scotch-Weld low odor acrylic adhesive accelerator part (reference DP 8810NS Green) with a viscosity of 35000 cP and density of 1.08 g/cc. Fluid B is 3M Scotch-Weld low odor acrylic adhesive base part (reference DP8810NS Green) with a viscosity of 90000 cP and density of 1.08 g/cc. The mix ratio by volume and weight is 10 part Fluid A: 1 part Fluid B. The CFD model is a particle mixing model with suitable simulation parameters including a mesh density of 2,377,000 cells, a time step of 0.005 seconds, a finite number of 18,700 fluid representing particles for fluid A, a finite number of 187,000 fluid representing particles for fluid B. Particle density distributions were calculated for fluid A and fluid B after mixing. The particle density distributions for fluid A and fluid B were converted to the respective spatial distributions of fluid concentration by using the equations (2) and (3) above with α=0.6 mm and rc=1.2 mm. A mixing index for the mixture was calculated based on the spatial distributions for fluid concentration.
A static mixer prototype was made by 3D printing by using the same CAD file. The particle laden experimental CT-scan after using mixing tip was developed to quantify the detailed mixing process. Copper metal flakes with particle size 325 Mesh (44 microns and finer) was premixed with 3M Scotch-Weld low odor acrylic adhesive DP8810 accelerator part with volume fraction of less than 2%. Then, the mixing experiment of the two-part adhesive using different mixer tips was carried out. After the experiment, the filled mixer tip is CT-scanned using a CT-scan machine commercially available from North Star Imaging (Rogers, Minn.) under the trade designation of X25 Industrial 3D X-Ray Inspection System. Due to the high density of the copper comparing to the adhesive material, the part with the copper particle will show higher density in the CT-scan.
A dynamic mixer was used to validate a computational fluid dynamics (CFD) simulation model. The geometry of the dynamic mixer was created via a CAD file. The CAD file was input to the CFD model for simulation of mixing fluid A and fluid B to form a mixture. Fluids A and B are the same as in Example 1. The CFD model is a particle mixing model with suitable simulation parameters including a mesh density of 2,817,314, a time step of 0.005 seconds, a finite number of 13,800 fluid representing particles for Fluid A, a finite number of 138,000 fluid representing particles for Fluid B. Particle density distributions were calculated for Fluid A and Fluid B after mixing. The particle density distributions for Fluid A and Fluid B were converted to the respective spatial distributions of fluid concentration by using the equations (2) and (3) above with σ=0.6 mm and rc=1.2 mm. A mixing index for the mixture was calculated based on the spatial distributions for fluid concentration.
A dynamic mixer prototype was made by 3D printing by using the same CAD file as for the simulation. The particle laden experimental CT-scan after using mixing tip was developed to quantify the detailed mixing process. Copper metal flakes with particle size 325 Mesh (44 microns and finer) was premixed with one part of the adhesive with volume fraction of less than 2%. Then, the mixing experiment of the two-part adhesive using different mixer tip was carried out. After the experiment, the filled mixer tip is CT-scanned. Due to the high density of the copper comparing to the adhesive material, the part with the copper particle will show higher density in the CT-scan.
The simulated fluid concentration as visualized in
The validated computational fluid dynamics (CFD) simulation model was used to optimize the operation and design of a dynamic mixer. Simulations were performed to output the mixing index along the length of the dynamic mixer for different design and operation conditions. For mixing fluids A and part B for an adhesive,
Unless otherwise indicated, all numbers expressing quantities or ingredients, measurement of properties and so forth used in the specification and embodiments are to be understood as being modified in all instances by the term “about.” Accordingly, unless indicated to the contrary, the numerical parameters set forth in the foregoing specification and attached listing of embodiments can vary depending upon the desired properties sought to be obtained by those skilled in the art utilizing the teachings of the present disclosure. At the very least, and not as an attempt to limit the application of the doctrine of equivalents to the scope of the claimed embodiments, each numerical parameter should at least be construed in light of the number of reported significant digits and by applying ordinary rounding techniques.
Exemplary embodiments of the present disclosure may take on various modifications and alterations without departing from the spirit and scope of the present disclosure. Accordingly, it is to be understood that the embodiments of the present disclosure are not to be limited to the following described exemplary embodiments, but is to be controlled by the limitations set forth in the claims and any equivalents thereof.
Exemplary embodiments are listed below. It is to be understood that any one of embodiments 1-19, 20-21 and 22-24 can be combined.
Embodiment 1 is a method of optimizing a mixing tool to mix a plurality of fluid materials to obtain a mixture, the method comprising:
representing a first geometry of the mixing tool with a first digital 3D model;
providing the first digital three-dimensional (3D) model to a computational fluid dynamics (CFD) simulation model implemented by a processor to simulate a mixing process of the plurality of fluids to generate a particle density distribution of the mixture inside the mixing tool;
converting the particle density distribution of the mixture to a first spatial distribution of fluid concentration of the mixture;
measuring a second spatial distribution of fluid concentration of the mixture inside the mixing tool when mixing the plurality of fluids using the mixing tool; and
comparing the first and second spatial distributions of fluid concentration to determine whether the first and second spatial distributions match with each other.
Embodiment 2 is the method of embodiment 1, wherein when the first and second spatial distributions match with each other, the CFD simulation model is validated.
Embodiment 3 is the method of embodiment 2, wherein the first and second spatial distributions are respectively represented by first and second sets of iso-surfaces, and when the respective shapes and distributions of the first and second sets of iso-surfaces match with each other, the CFD simulation model is validated.
Embodiment 4 is the method of embodiment 2 or 3, further comprising optimizing the first geometry of the mixing tool to a second geometry using the validated CFD simulation model.
Embodiment 5 is the method of embodiment 4, further comprising representing the second geometry with a second digital 3D model.
Embodiment 6 is the method of embodiment 5, further comprising providing the second digital 3D model to the validated CFD simulation model to simulate the mixing of the plurality of fluids, and implementing, via the processor, the CFD simulation model to generate an updated particle density distribution of the mixture inside the mixing tool.
Embodiment 7 is the method of embodiment 6, further comprising converting the updated particle density distribution of the mixture to an updated spatial distribution of fluid concentration of the mixture.
Embodiment 8 is the method of embodiment 7, further comprising determining a mixing spatial distribution of index based on the updated spatial distribution of fluid concentration of the mixture.
Embodiment 9 is the method of embodiment 8, further comprising visualizing at least one of the spatial distributions of fluid concentration and mixing index in a graphic user interface to guide a user to determine whether the plurality of fluids is uniformly mixed.
Embodiment 10 is the method of any one of embodiments 1-9, wherein when the first and second spatial distributions do not match with each other, adjusting the CFD simulation model to generate an updated particle density distribution of the mixture.
Embodiment 11 is the method of any one of embodiments 1-10, wherein converting the particle density distribution of the mixture to the first spatial distribution of fluid comprises calculating a fluid concentration at a given point by weighting adjacent discrete fluid representing particles of the particle density distribution.
Embodiment 12 is the method of any one of embodiments 1-11, wherein the comparison comprises visualizing the first and second spatial distributions of fluid concentration, optionally comprising overlaying digital representations of the first and second first and second spatial distributions in a graphic user interface.
Embodiment 13 is the method of embodiment 12, wherein overlaying the digital representations of the first and second first and second spatial distributions comprises importing the corresponding polygon surfaces to the same coordinate system in the graphic user interface.
Embodiment 14 is the method of any one of embodiments 1-13, wherein the mixing of the plurality of fluids via the simulation with the CFD simulation model and via the mixing tool is under the same operation conditions.
Embodiment 15 is the method of embodiment 14, wherein the same operation conditions include the same flow parameters.
Embodiment 16 is the method of any one of embodiments 1-15, wherein the plurality of fluids includes two or more specifies of an adhesive.
Embodiment 17 is the method of any one of embodiments 1-16, wherein the second spatial distribution of fluid concentration is measured via an X-ray scan, optionally including a particle laden experimental CT-scan of the mixing tool.
Embodiment 18 is the method of embodiment 17, wherein copper metal flakes are used for the particle laden experimental CT-scan.
Embodiment 19 is the method of any one of embodiments 1-18, further comprising physically testing the mixture.
Embodiment 20 is a computer-implemented method to visualize a spatial distribution of fluid concentration inside a tool, the method comprising:
representing a first geometry of the tool with a first digital 3D model;
providing the first digital 3D model to a computational fluid dynamics (CFD) simulation model implemented by a processor to simulate a fluid flow inside the tool to generate a particle density distribution of one or more fluids inside the tool;
converting the particle density distribution to a spatial distribution of fluid concentration of the mixture; and
visualizing the spatial distribution of fluid concentration in a graphic user interface (GUI).
Embodiment 21 is the method of embodiment 20, further comprising determining a spatial distribution of mixing index based on the spatial distribution of fluid concentration, and optionally visualizing the spatial distribution of mixing index in a graphic user interface (GUI).
Embodiment 22 is a computer-implemented system to design and optimize a tool, comprising:
a module to represent a first geometry of the tool with a first digital 3D model;
a module to provide the first digital 3D model to a computational fluid dynamics (CFD) simulation model implemented by a processor to simulate a fluid flow inside the tool to generate a particle density distribution of one or more fluids inside the tool;
a module to convert the particle density distribution to a spatial distribution of fluid concentration of the mixture; and
a module to visualize the spatial distribution of fluid concentration and a measured spatial distribution of fluid concentration in a graphic user interface (GUI).
Embodiment 23 is the computer-implemented system of embodiment 22, wherein the module to module to convert the particle density distribution is configured further to determine a spatial distribution of mixing index based on the spatial distribution of fluid concentration.
Embodiment 24 is the computer-implemented system of embodiment 22 or 23, wherein the module to visualize the spatial distribution of fluid concentration is configured further to visualize the spatial distribution of mixing index in a graphic user interface (GUI).
Reference throughout this specification to “one embodiment,” “certain embodiments,” “one or more embodiments,” or “an embodiment,” whether or not including the term “exemplary” preceding the term “embodiment,” means that a particular feature, structure, material, or characteristic described in connection with the embodiment is included in at least one embodiment of the certain exemplary embodiments of the present disclosure. Thus, the appearances of the phrases such as “in one or more embodiments,” “in certain embodiments,” “in one embodiment,” or “in an embodiment” in various places throughout this specification are not necessarily referring to the same embodiment of the certain exemplary embodiments of the present disclosure. Furthermore, the particular features, structures, materials, or characteristics may be combined in any suitable manner in one or more embodiments. While the specification has described in detail certain exemplary embodiments, it will be appreciated that those skilled in the art, upon attaining an understanding of the foregoing, may readily conceive of alterations to, variations of, and equivalents to these embodiments. Accordingly, it should be understood that this disclosure is not to be unduly limited to the illustrative embodiments set forth hereinabove. Furthermore, various exemplary embodiments have been described. These and other embodiments are within the scope of the following claims.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2020/061085 | 11/24/2020 | WO |
Number | Date | Country | |
---|---|---|---|
62940613 | Nov 2019 | US |