Next Article in Journal
Effects of Coating Parameters of Hot Filament Chemical Vapour Deposition on Tool Wear in Micro-Drilling of High-Frequency Printed Circuit Board
Next Article in Special Issue
Design and Fabrication Process Optimization of Silver-Based Inkjet-Printed Microheater
Previous Article in Journal
A Study on the Application of Discrete Wavelet Decomposition for Fault Diagnosis on a Ship Oil Purifier
Previous Article in Special Issue
Suppression and Utilization of Satellite Droplets for Inkjet Printing: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling 3D Droplet Movement Using a Drop-on-Demand Inkjet Printhead Model

by
Tim Tofan
1,
Sergejus Borodinas
2,3,
Rimantas Kačianauskas
2,3 and
Raimondas Jasevičius
1,3,*
1
Department of Mechanics and Materials Engineering, Vilnius Gediminas Technical University, J. Basanavicius Street 28, 03224 Vilnius, Lithuania
2
Department of Applied Mechanics, Vilnius Gediminas Technical University, Saulėtekio al. 11, 10223 Vilnius, Lithuania
3
Institute of Mechanical Science, Vilnius Gediminas Technical University, J. Basanavicius Street 28, 03224 Vilnius, Lithuania
*
Author to whom correspondence should be addressed.
Processes 2022, 10(8), 1467; https://doi.org/10.3390/pr10081467
Submission received: 10 July 2022 / Revised: 20 July 2022 / Accepted: 25 July 2022 / Published: 27 July 2022
(This article belongs to the Special Issue Advanced Research and Applications of Inkjet Printing (IJP) Technique)

Abstract

:
This article presents a numerical simulation of a printhead model for drop-on-demand (DoD) inkjet printers. A three-dimensional droplet model is provided for the numerical study of inks, ejection parameters, droplet movement, and the analysis of droplet impacts on the surface. This work is devoted to the analysis of different droplet ejection settings during the printing process, when the behavior of the droplet directly affects the accuracy of the printing process itself. A numerical model was also developed to investigate the effect of various settings on droplet stability, including printhead size and nozzle orifice, motion parameters (pulse strength and droplet ejection amplitude) and fluid properties. The results reflect the behavior of the ink droplet over time. The behavior of the drop was tested at different waveform ejection parameters and a mass turnover was observed.

1. Introduction

Nowadays, inkjet printing technology has become more widely adopted by manufacturers due to its simplicity and flexibility. There is increasing demand for paper-based products, including in such areas as T-shirt printing, microchip production or personalized production, while meeting increasing quality requirements.
An image analysis study of ink–paper interactions was conducted by Fiadeiro et al. [1]. It is important to identify variations in the interaction of ink and paper, manifested in the machine and cross directions of the paper (Mendes et al. [2]). The influence of paper properties on the print quality of an inkjet printer is described by Bandyopadhyay [3]. Skowronski [4] analyzed the complexity of the interaction between ink and paper in inkjet printing and their influence on the properties of printed paper [4]. The initial coating of the paper remains important. The interaction between offset ink and coated paper is described by Ström [5]. Liu and Derby [6] analyzed the droplet formation of inkjet printers using different printing ink properties and actuation strength (voltage). A physical experiment was carried out to determine the possibility of ejection of well-shaped drops without the formation of satellite drops. The influence of voltage parameters on the piezoelectric ejection process of ink droplets with different rheological properties was investigated and analyzed by Jiao et al. [7].
The droplet ejection process in inkjet printers is modeled by a single dimension presented by Jiang and Tan [8]. Droplet stability after ejection was analyzed by Zhong et al. [9]. To control the formation of droplets, Wang and Chiu [10] investigated data-driven models of droplet formation. To improve the print quality by increasing the dots per inch (DPI) of the piezoelectric inkjet printhead, the width and length of the nozzle chamber become smaller. In this situation, the formation of droplets with different waveform times was performed by Wei et al. [11]. The printing quality is directly related to the print excitation parameters (Qiu-min et al. [12]). They mentioned that the piezoelectric (PZT) driving system is an important part in inkjet printing equipment. The main function is to eject droplets of an ideal shape and volume at the correct velocity and without satellites to reach the surface at the intended location. A machine learning approach for predicting the maximum spreading factor of droplets upon impact on surfaces with various wettability’s were investigated by Tembely et al. [13]. They found that machine learning algorithms can more accurately predict the experimental results of the maximum spreading factor. Increasing the quality of DoD inkjet drops was analyzed by Oktavianty et al. [14]. They mentioned that it is necessary to establish a flexible waveform design for different fluid compositions. The influence of excitation on dynamic characteristics of piezoelectric micro-jets was investigated by Li et al. [15]. They demonstrated that the direct coupling method proves to be an effective and feasible method to quantitatively simulate the fluid dynamic characteristics of the micro-jet. An evaluation of the relationship between interdigital geometry quality and inkjet parameters was analyzed by Bertolucci et al. [16]. They developed a methodology for assessing the relationships between geometrical quantities and printing parameters. The effect of ink supply pressure on the piezoelectric inkjet was investigated by Kim et al. [17]. They mentioned that decreasing the ink supply pressure can be expected to improve the stability and productivity of inkjet printing. Simulations and sensitivity analysis of piezoelectric inkjet printheads were provided by Nguyen et al. [18]. They developed a lumped element model (LEM) for an inkjet head to simulate its droplet generation in inkjet printing. The inkjet printing of PEDOT:PSS-based conductive patterns for 3D forming applications was investigated by Basak et al. [19]. They developed a polymer PEDOT: PSS-based ink formulation that was inkjet printable through a 21 µm orifice. The effect of process parameters on the performance of drop-on-demand 3D inkjet printing was analyzed by Elkaseer et al. [20]. They developed an optimization-oriented simulation tool of droplet behavior during the drop-on-demand 3D inkjet printing process.
In this paper, we perform numerical simulation using the COMSOL multiphysics finite element method and analyze the ejection of droplet through a nozzle into the air, the shape turnover of droplet, and the distribution of satellites during the fall. Inkjet printing inks are calculated as Newtonian fluid. The flow of ink is modeled by the Navier–Stokes equations for an incompressible fluid. The simulation model presented is solved using the level set method to determine the ink–air interface.

2. Problem Formulation

There are many inkjet printers that use inkjet printheads that move in one direction, making the printers easier to use. The print material is moving in the opposite direction of the printhead.
It can also be noted that there are many articles (Introduction section) on similar topics related to inkjet printing, and these topics remain relevant, leading not only to improvements in the printing process itself, but also to such fundamental things as the fundamental behavior of a liquid drop during its formation, movement into air and interaction with a solid surface.
Our study is characterized by the principles of numerical simulation and aspects of their optimization; when using finite elements, we study the chamber of the nozzle of the inkjet printhead, the ejection parameters, the fall of the droplet on the surface and the impact on it. In addition, similar inkjet printing processes are investigated by Elkaseer et al. [20] and Wilson et al. [21].
Elkaseer et al. [20] are investigating inkjet printing technology as an alternative manufacturing method for producing functional parts from multi-material, such as micro-schemes.
They simulated and investigated the printing process, which takes into account the effects of droplet volume, droplet center-to-center distance, jetted droplets coverage percentage, ink–solid substrate contact angle, and overlapping droplets coalescence performance, as well as the number of printed layers. The simulation tool was created using the MATLAB programming language. The simulation results were experimentally validated, and the results showed good agreement with the maximum deviation for the horizontal features. The method they developed was used to find out the effect of processed TIFF image resolution and droplet diameter on the dimensional accuracy of inkjet 3D printed parts [20].
For our further research, it is useful to use a geometric-based approach to simulate the coalescence of droplets on a substrate, taking into account the resolution and percentage of coverage of the printed image file and the droplet diameter.
Wilson et al. [21] studied inkjet process optimization for a laboratory inkjet printing system designed to print with polymer inks, with a focus on transparent and electrically conductive polymers. In this study, the authors investigated the speed of droplets from a nozzle as a function of the input voltage waveform and dwelling time difference. With the correct ejection parameters, the authors obtained a continuous line from droplets on the polymer surface. This study [21] presents a new model incorporating the effects of inertia, surface tension, and ink viscoelasticity. Due to the increased viscosity of the ink at low rates, the ink becomes more shear thinning. These studies can be used for our further research on the similar effect of inks on the inkjet process.
It should be noted that when using COMSOL Multiphysics to conduct a numerical experiment, it is necessary to determine the detailed paint parameters. The corresponding geometric parameters of the head are also important. Several of the referenced publications [20,21] do not use COMSOL Multiphysics, so comparing the results may be one of the future tasks.

3. Methodology

Based on standard PZT printhead models, a single inkjet printhead nozzle chamber was modeled using COMSOL Multiphysics software. The geometry of the three-dimensional numerical model of the inkjet printhead nozzle is shown in Figure 1.
A droplet is formed by jetting a fixed amount of ink into an air-filled cylindrical printing chamber with a length L = 1   mm and a radius R = 0.1   mm . The radius of the jetting channel is equal to r 1 = 0.1   mm . The initial circular shape and the cross-sectional radius r 2 = 12   μ m of the droplet are predefined by the channel shape and size.
The waveforms are generated by the top inlet of the nozzle by applying the fluid intake velocity and various waveform dwelling times. Different waveforms dwelling times are shown at Figure 2.
Here, the printhead was actuated by an optimized waveform with different dwelling times t 1 = 1   µ s , t 2 = 3   µ s , and t 3 = 12   µ s . The considered droplet ejection times were set to demonstrate differences in droplet behavior. Typically, PZT inkjet printheads use different dwelling times to control droplet size on the surface. To better understand the droplet behavior and impact differences, we chose the different dwelling times mentioned above. As the results will show, these are the time intervals when the droplet formed on the printed surface has a surface close to the shape of a hemisphere, with a dwelling time of 1 µs, and the spherical surface of the droplet is completely lost with a dwelling time of 12 µs; an additional intermediate dwelling time interval was chosen to be 3 µs. The ink inlet velocity was V in = 0.5   m s . When the pressure waveform appears, the liquid begins to push out the ink through the printhead nozzle, forming the jet. Then the jet pinches off from the nozzle, the main droplet begins to form, and peripheral satellites may appear.
Numerical simulation based on the continuity equation and the N–S equation reconcile two-phase flows of incompressible substances. The continuity equation for an incompressible fluid is calculated using the following equation:
· u = u x x + u y y + u z z = 0
Here, ∇ is the vector differential operator, and u is the fluid velocity. From the definitions of calculus, we know that the gradient operator (∇) is calculated by the equation
= x i + y j + z k
where i, j, k are the unit vectors in the directions of the x, y and z coordinates.
The N–S equation describes the momentum conservation equation for a viscous incompressible fluid in motion. This equation is written as follows:
ρ ( u · ) u = μ ( ( · u + ( · u T ) ) ) p + ρ g
Here, ρ is the substance density; μ is the substance viscosity; p is the pressure; and g is the acceleration due to gravity.
The following equation describes the interface between the ink and air phases. The convection of the reinitialized level set function is written as the following:
ϕ t + u · ϕ = γ · ( ε ϕ ϕ ( 1 ϕ ) ϕ | ϕ | )
Here, the coefficient ϕ represents the level set interface between ink and air. Accordingly,   ϕ = 0 represents the air environment, and ϕ = 1 represents ink. In the transition layer close to the interface, ϕ smoothly changes from 0 to 1. The fluid velocity, u , represents the motion of interface during the environments interplay. The parameter γ determines the amount of stabilization of the level set function. A suitable value for γ is the maximum magnitude occurring in the velocity field.
By default, this reinitialization parameter γ for the current simulation is set to γ = 10   m s . Here, ε refers to the thickness of the transition layer. For this simulation, ε is taken as half of the representative mesh size in the area that the falling droplet passes through. For this simulation, the interface thickness is controlled by the default and set to 1.25 µm.
By defining the interface of different substrates, the level set function is used to stabilize the density and viscosity calculations through the interface. In numerical simulation, it is necessary to take into account the dynamic viscosity (shear viscosity) µ . It describes the interface between the shear rate and shear stresses in the substrate. The substrate density and dynamic viscosity in different finite elements depend on the interface coefficient ϕ . Density and dynamic viscosity are written as follows:
ρ = ρ 1 + ( ρ 2 ρ 1 ) ϕ
µ = µ 1 + ( µ 2 µ 1 ) ϕ
Here, ρ 1 indicates the air density and ρ 2 indicates the ink density; µ 1 indicates the air viscosity and µ 2 indicates the ink viscosity.
In numerical simulation, ink and air are investigated as a substrate, from the flow of which viscous stresses arise at each point. The N–S equations describe the mass and momentum transfer of an incompressible fluid, including the surface tension and gravity forces. The motion of the droplet, taking into account the conservation of mass, is calculated from the following differential equations:
( · u ) = 0
ρ ( u t + u · u g ) · ( μ ( u + u T ) ) + p = F st
The surface tension force is calculated by the equation
F st = σ δ ĸ n
Here, δ is the Dirac delta function; σ is the surface tension coefficient; ĸ is the curvature; and n is the unit vector (normal direction). The unit vector is calculated from the following equation:
n = ϕ | ϕ |
The delta function is approximated by the following equation:
Δ = 6 | ϕ ( 1 ϕ ) | | ϕ |
To calculate the droplet motion, including mass conservation, the adaptive mesh refinement function is used to locally improve the mesh at the ink–air interface. For this simulation, the finite elements sequence and turnover are shown at Figure 3.
This functionality will essentially divide the simulation progress into several time intervals and locally refine the mesh in the region, where the phase interface is present in each interval to improve the accuracy of the calculation. To begin with, we have about 220,000 finite elements, and during the ejection, the mesh of the droplet changes 12 times. At the second mesh refinement, the number of finite elements reaches 800,000. The average quality of the elements decreases with an increase in the number of finite elements.
The final droplet mass can be calculated using the following COMSOL Multiphysics numerical simulation equation:
m d = ρ ( V · ϕ )
There are several dimensionless parameters that are very important to the inkjet printhead working process. The Reynolds number Re determines the ratio of a fluid’s inertia to its viscosity, while the Weber number We determines the ratio of inertia to its surface tension (Seipel et al. [22]). These dimensionless numbers are calculated using the following equations:
Re = ρ dv d μ
  We = ρ dv d 2 σ
Here, v d is the droplet velocity and d is the length of the nozzle diameter.
Several effects are observed during the motion of the ink droplet. A typical picture of a droplet, a fragment of the mesh and a contour plot of the velocity magnitude accompanied with the Cell Reynold’s number are shown in Figure 4.
The influence of droplet velocity in the inkjet printing process can be eliminated by forming the Ohnesorge number Oh . There is a solution based on the N–S equations for expressing droplet ejection limitations considering the interfacial, viscous and inertial properties of the fluid (Seipel et al. [22]):
Z = σ ρ d μ = 1 Oh = Re We  
The Oh value can only be calculated from the physical properties of the inks and the nozzle size scale. It is independent of the waveform extraction impulse that controls the velocity. Oh closely affects the behavior of an expanding droplet from a printhead nozzle.
The initial ink printability specification is required for ejection of stable droplets and should be in the acceptable range of 1 < Z < 10, (Fromm [23]; Derby and Reis [24]). If the number Z > 10, then during the printing process, satellite droplets will form, which degrading the print quality. If the number Z < 1, then the viscous forces of the fluid will prevent the expansion of the droplet.
Proper operation of the printing process in a DoD inkjet printhead requires a propriate combination of physical parameters, which will also depend on the droplet size and velocity (via the Reynolds or Weber number value). The results entail empirical data that are embedded in numerical simulation calculations based on the dimensionless numbers We, Re and Z.

4. Experiment and Results

The physical parameters of ink and air used in numerical simulation of the inkjet printhead nozzle are presented in Table 1.
The droplet ejection progress was taken from numerical simulation using COMSOL Multiphysics software. The momentum strength of the droplet ejection is predicted from the different velocity at the fluid inlet. For these simulations, the inlet velocity was set to 0.5 m/s. The presented graphs with different waveform dwelling times are shown in Figure 5.
When the dwelling time is 1 µs, the final mass of the droplet is 2 × 10−13 kg. The droplet forms a geometrical surface close to a hemisphere after the pitch off the nozzle and reaches the surface approximately 40 µs from the start of actuation. After increasing the dwelling period to 3 µs, the mass of the droplet increases to 3 × 10−12 kg. Larger droplets also form a geometrical surface close to a hemisphere and reach the surface faster by 35 µs. After increasing the formatting time to 12 µs, the droplet is ejected with a large tail and does not form a hemisphere when falling. After the impact with a surface, the droplet splashes and does not form the final round shape on the surface. The final mass before droplet impact with a surface and splashing is 1.4 × 10−10 kg.
The final mass of the droplet depends on the different dwelling time period. The mass turnover of droplets is shown in Figure 6. Here, us and ng stand for microseconds (µs) and nanograms (ng), respectively.
For this research, we investigate the droplet mass formatting with different ejection times. During the droplet formatting periods, the droplet mass increases in accordance with the effect of the inlet velocity. After the waveform decreases, the droplet continues to move forward toward the surface and pitch off from the nozzle.
If the inlet velocity was too low, the droplet would shrink to the nozzle and not move forward to the surface, and the kinetic energy of the droplet would be too low. After the droplet reaches the surface, the droplet splashes and we can see a sudden change in the graph. For a large formatting time of 12 µs, we see that a large mass grows in after 35–40 µs. This is caused by splashing and the droplet jet replacing the main droplet.

5. Conclusions

At this research, we investigate the droplet ejection turnover over time. The main aim is to understand the behavior of the droplet at different waveform ejection parameters. Mass turnover is also observed. It is important to investigate the trajectory of the fall of the droplet and the mass of the droplet to understand the final place and diameter of the droplet on the surface. These aspects affect the final printing quality. Theoretical studies of three-dimensional numerical simulation help to evaluate such a component as the waveform dwelling time. The different dwelling times indicate a direct effect on the impact of the droplet on the surface. The masses of the ejected droplet before impact with the surface at different dwelling periods are as follows: 1   µ s     2.5 × 10 3   µ g ; 3   µ s     16 × 10 3   µ g ; and 12   µ s   65 × 10 3   µ g . This numerical simulation is an important alternative to quality improvement analysis; further research should be focused on the impact of droplets on different surfaces.

Author Contributions

Conceptualization, methodology, investigation, writing—original draft preparation T.T., S.B., R.K. and R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of UNIVERSITAT POLITÈCNICADE VALÈNCIA (date of approval: 23 October 2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Data Availability Statement: The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Data Availability Statement

The data presented in this article is available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

DoDdrop-on-demand;
N-SNavier–Stokes;
PZTpiezoelectric;
TIFFtagged image file format.
Roman
dis the diameter of the nozzle (mm);
F st is the surface tension force (N);
gis the acceleration due to gravity (m/s²);
L is the air length between nozzle and surface (mm);
m d is the droplet mass (kg);
nis the unit vector in the normal direction;
Ohis the Ohnesorge number, a dimensionless coefficient;
pis the pressure (Pa);
Ris the air chamber radius (mm);
r 1 is the inlet channel radius (mm);
r 2 is   the   jetting   nozzle   radius   ( μ m ) ;
Reis the Reynolds number, a dimensionless coefficient;
tis time (s);
V d is the droplet volume (pL);
V in is the inlet ink velocity (m/s);
Weis the Weber number, a dimensionless coefficient;
Zis the ink printability range, a dimensionless number.
Greek
βis the maximum spreading factor;
γis the parameter which determines the repetition of initiations;
Δis the delta function;
is the vector differential operator;
δis the Dirac delta function;
εis the representative mesh size in the area passed by the falling droplet;
𝜂is fluid viscosity (m/s);
ĸ is the curvature;
μis the dynamic viscosity (N·s/m²);
ρis the density (kg/m3);
σis the surface tension coefficient (mN/m);
ϕis the coefficient of level set interface between air and ink;
v d is the droplet velocity (m/s).

References

  1. Fiadeiro, P.T.; Mendes, A.O.; Ramos, A.M.M.; Sousa, S.C.L. Study of the ink-paper interaction by image analysis: Surface and bulk inspection. In Proceedings of the 8th Iberoamerican Optics Meeting and 11th Latin American Meeting on Optics, Lasers, and Applications, Porto, Portugal, 22–26 July 2013. [Google Scholar] [CrossRef]
  2. Mendes, A.O.; Fiadeiro, P.T.; Ramos, A.M.M.; de Sousa, S.C.L. Development of an optical system for analysis of the ink–paper interaction. Mach. Vis. Appl. 2013, 24, 1733–1750. [Google Scholar] [CrossRef]
  3. Bandyopadhyay, S. Effect of paper properties on print quality of ink jet printer. In Proceedings of the IS&T NIP17: International Conference on Digital Printing Technologies, Fort Lauderdale, FL, USA, 30 September–5 October 2001; pp. 141–143. [Google Scholar]
  4. Skowronski, J.W.A. Clarification of the Complexity of Ink–Paper Interactions in Inkjet Printing and their Impact on Printed Paper Properties. Zeszyty Naukowe. Rozprawy Naukowe/Politechnika Łódzka 2009, 374, 110. [Google Scholar]
  5. Ström, G. Review: Interaction between offset ink and coated paper–a review of the present understanding. In Advances in Paper Science and Technology, Transactions of the XIIIth Fundamental Research Symposium, Cambridge, 2005; I’Anson, S.J., Ed.; FRC: Manchester, UK, 2018; pp. 1101–1137. [Google Scholar] [CrossRef]
  6. Liu, Y.; Derby, B. Experimental study of the parameters for stable drop-on-demand inkjet performance. Phys. Fluids 2019, 31, 032004. [Google Scholar] [CrossRef] [Green Version]
  7. Jiao, T.; Lian, Q.; Zhao, T.; Wang, H. Influence of ink properties and voltage parameters on piezoelectric inkjet droplet formation. Appl. Phys. 2021, 127, 11. [Google Scholar] [CrossRef]
  8. Jiang, H.; Tan, H. One dimensional model for droplet ejection process in inkjet devices. Fluids 2018, 3, 28. [Google Scholar] [CrossRef] [Green Version]
  9. Zhong, Y.; Fang, H.; Ma, Q.; Dong, X. Analysis of droplet stability after ejection from an inkjet nozzle. J. Fluid Mech. 2018, 845, 378–391. [Google Scholar] [CrossRef]
  10. Wang, J.; Chiu, G.T.C. Data-driven drop formation modeling in nanoliter drop-on-demand inkjet printing. In Proceedings of the ASME 2020 Dynamic Systems and Control Conference, Online, 5–7 October 2020; Volume 84287, p. V002T28A002. [Google Scholar] [CrossRef]
  11. Wei, H.; Xiao, X.; Yin, Z.; Yi, M.; Zou, H. A waveform design method for high DPI piezoelectric inkjet print-head based on numerical simulation. Microsyst. Technol. 2017, 23, 5365–5373. [Google Scholar] [CrossRef]
  12. Qiu-Min, W.U.; Xin-Yu, C.U.I.; Lei, X.U.; Xin, S.U. Research on driving waveform and residual oscillation in piezo inkjet. Res. Sq. 2022. [Google Scholar] [CrossRef]
  13. Tembely, M.; Vadillo, D.C.; Dolatabadi, A.; Soucemarianadin, A. A machine learning approach for predicting the maximum spreading factor of droplets upon impact on surfaces with various wettabilities. Processes 2022, 10, 1141. [Google Scholar] [CrossRef]
  14. Oktavianty, O.; Haruyama, S.; Ishii, Y. Enhancing droplet quality of edible ink in single and multi-drop methods by optimization the waveform design of DoD inkjet printer. Processes 2022, 10, 91. [Google Scholar] [CrossRef]
  15. Li, K.; Liu, J.-K.; Chen, W.-S.; Zhang, L. Influences of excitation on dynamic characteristics of piezoelectric micro-jets. Micromachines 2017, 8, 213. [Google Scholar] [CrossRef] [PubMed]
  16. Bertolucci, F.; Berdozzi, N.; Rebaioli, L.; Patil, T.; Vertechy, R.; Fassi, I. Assessing the relationships between interdigital geometry quality and inkjet printing parameters. Micromachines 2022, 13, 57. [Google Scholar] [CrossRef] [PubMed]
  17. Kim, S.; Choi, J.H.; Sohn, D.K.; Ko, H.S. The effect of ink supply pressure on piezoelectric inkjet. Micromachines 2022, 13, 615. [Google Scholar] [CrossRef] [PubMed]
  18. Nguyen, V.-T.; Leong, J.Y.C.; Watanabe, S.; Morooka, T.; Shimizu, T. A multi-fidelity model for simulations and sensitivity analysis of piezoelectric inkjet printheads. Micromachine 2021, 12, 1038. [Google Scholar] [CrossRef] [PubMed]
  19. Basak, I.; Nowicki, G.; Ruttens, B.; Desta, D.; Prooth, J.; Jose, M.; Nagels, S.; Boyen, H.G.; D’Haen, J.; Buntinx, M.; et al. Inkjet printing of PEDOT:PSS based conductive patterns for 3D forming applications. Polymers 2020, 12, 2915. [Google Scholar] [CrossRef] [PubMed]
  20. Elkaseer, A.; Schneider, S.; Deng, Y.; Scholz, S.G. Effect of process parameters on the performance of drop-on-demand 3D inkjet printing: Geometrical-based modeling and experimental validation. Polymers 2022, 14, 2557. [Google Scholar] [CrossRef] [PubMed]
  21. Wilson, P.; Lekakou, C.; Watts, J.F. System design and process optimization for the inkjet printing of PEDOT: Poly (styrenesulfonate). J. Micro Nano-Manuf. 2014, 2. [Google Scholar] [CrossRef]
  22. Seipel, S.; Yu, J.; Nierstrasz, V.A. Effect of physical parameters and temperature on the piezo-electric jetting behaviour of UV-curable photochromic inks. Sci. Rep. 2020, 10, 18841. [Google Scholar] [CrossRef] [PubMed]
  23. Fromm, J.E. Numerical calculation of the fluid dynamics of drop-on-demand jets. IBM J. Res. Dev. 1984, 28, 322–333. [Google Scholar] [CrossRef]
  24. Derby, B.; Reis, N. Inkjet printing of highly loaded particulate suspensions. MRS Bull. 2003, 28, 815–818. [Google Scholar] [CrossRef]
Figure 1. Inkjet printhead nozzle geometry.
Figure 1. Inkjet printhead nozzle geometry.
Processes 10 01467 g001
Figure 2. The inlet velocity for inkjet printhead nozzle (determined by waveform functions).
Figure 2. The inlet velocity for inkjet printhead nozzle (determined by waveform functions).
Processes 10 01467 g002
Figure 3. Representation of an adaptive mesh refinement sequence.
Figure 3. Representation of an adaptive mesh refinement sequence.
Processes 10 01467 g003
Figure 4. Illustration of the shape of a droplet in motion with a mesh and distribution of velocities.
Figure 4. Illustration of the shape of a droplet in motion with a mesh and distribution of velocities.
Processes 10 01467 g004
Figure 5. Representation of the movement of ejected droplet with different ejection parameters. Waveform dwelling time: 1 µs (a); 3 µs (b); 12 µs (c). Inlet velocity V in = 0.5   m s .
Figure 5. Representation of the movement of ejected droplet with different ejection parameters. Waveform dwelling time: 1 µs (a); 3 µs (b); 12 µs (c). Inlet velocity V in = 0.5   m s .
Processes 10 01467 g005aProcesses 10 01467 g005b
Figure 6. Representation of the droplet mass turnover. Mass versus formatting time.
Figure 6. Representation of the droplet mass turnover. Mass versus formatting time.
Processes 10 01467 g006
Table 1. Physical parameters of the substances used in the numerical modeling of the inkjet printhead nozzle.
Table 1. Physical parameters of the substances used in the numerical modeling of the inkjet printhead nozzle.
Substanceρ
kg/m³
μ
mN·s/m²
σ
mN/m
OhZ
Air1.2251.789 × 10−2---
Ink 111014720.323.13
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tofan, T.; Borodinas, S.; Kačianauskas, R.; Jasevičius, R. Modeling 3D Droplet Movement Using a Drop-on-Demand Inkjet Printhead Model. Processes 2022, 10, 1467. https://doi.org/10.3390/pr10081467

AMA Style

Tofan T, Borodinas S, Kačianauskas R, Jasevičius R. Modeling 3D Droplet Movement Using a Drop-on-Demand Inkjet Printhead Model. Processes. 2022; 10(8):1467. https://doi.org/10.3390/pr10081467

Chicago/Turabian Style

Tofan, Tim, Sergejus Borodinas, Rimantas Kačianauskas, and Raimondas Jasevičius. 2022. "Modeling 3D Droplet Movement Using a Drop-on-Demand Inkjet Printhead Model" Processes 10, no. 8: 1467. https://doi.org/10.3390/pr10081467

APA Style

Tofan, T., Borodinas, S., Kačianauskas, R., & Jasevičius, R. (2022). Modeling 3D Droplet Movement Using a Drop-on-Demand Inkjet Printhead Model. Processes, 10(8), 1467. https://doi.org/10.3390/pr10081467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop