Droplet Trajectory Movement Modeling Using a Drop-on-Demand Inkjet Printhead Simulations

Abstract

Previous studies of the authors were focused on the vertical movement of the jet print when the printed head was stationary. In this work, the following study was presented, in which the movement of droplets is achieved using a moving horizontal print head. The printed head moves at various velocities, which affects the time of printing and deposition accuracy. This study provides a 3D numerical model with a complete turnover/interchange of the droplet shape at different time steps during the formation and movement process. By considering the dynamics of a droplet surrounded by air, we modeled them using the two-phase flow coupling and level set function from the computational fluid dynamics module by COMSOL Multiphysics. The trajectory shifts of the inkjet droplet are considered from its ejection to its impact on the surface at each time step. The conclusions summarize all the factors responsible for the trajectory shift of the droplet during vertical fall.

1. Introduction

Ink trajectory. With the modernization of modern digital printing, research on the trajectory of the ink droplet movement remains relevant. The trajectory of droplets of inkjet, experiment, and modeling is given by Mohebi and Evans [1]. They used a curve-fitting method to determine the drag coefficient for droplets (small size) in inkjet printing and to evaluate the velocity at exit. Considering the use of the application of inkjet printing, the analysis of the flying droplet and the prognosis of its landing on a non-flat substrate is given by Arango et al. [2]. They mentioned that the density of the liquid affects the loss of the rate of falling droplets. In addition, the functional fluid’s use is limited due to physicochemical conditions forcing the Weber number (the range is not large) [2]. Wang [3] studied the aerodynamic influence on the location of the inkjet droplet and satellite. He mentioned that numerical modeling is a powerful tool in the study and development of an inkjet product. A study on the distribution of an electric field and the trajectory of droplets during electrohydrodynamic jet printing was given by Wu et al. [4]. During research, based on the study of the distribution of the electric field, the analysis of the speed and direction of the ejected droplets was confirmed by plotting their trajectories. The movement of the droplets of ink in wide-format printers was investigated by Link et al. [5]. They mentioned that, in accordance with the results, it is possible to significantly increase the relative movement of the head media without causing significant turbulence with the subsequent flow interference, which can worsen quality imprints. To examine the droplet ejection process for a piezoelectric (PZT) inkjet printhead, a numerical simulation model is presented by Yang [6]. His presented theoretical model takes account of a set of three-dimensional, time-dependent conservation equations of mass and momentum, with the incorporation of the continuous surface force model for treating the interfacial surface tension effect. Inkjet printing technologies are presented by Hudd [7]. He mentioned that the accuracy of placing the droplet depends on the exact position of the landing. Placement accuracy is uncertain due to various parameters, such as the nozzle and moistening of the surface, variations of the jet-to-jet, sensitivity to nozzle straightness, single jet-over-time, contamination of the nozzle plate, the condition and formulation, as well as the rate of the droplet. This problem is worse for longer flight paths or “throwing distance”.

Implementation/application of inkjet printing.

Manufacturers are increasingly utilizing inkjet printing technology due to its flexibility. Increasing demand for personalized production, such as small print run papers, souvenir production, textile printing, and even microchip manufacturing, meeting high-quality standards. Formalization and improvement of the design rules for inkjet printing are given by Mashayekhi et al. [8]. For improving design rules and scaling down device dimensions, they present a design approach that combines pre-patterned, high-resolution substrates with digital inkjet fabrication as a demonstration of the capabilities of combining inkjet with other fabrication technologies. The cons and pros of inkjet technology in industrial inkjet printing are outlined in Zapka’s work [9]. He also gave information on the comparison of inkjet and screen technologies.

Other implementation/application. There are many other implementations of the presented technology where ink can play an important role. Nowadays, trends in pharmaceutical and medical applications of inkjet printing technology are outlined by Skoutaris et al. [10]. They noted that inkjet printing could help create drugs targeted to specific patient groups and could also be used in bioprinting in regenerative medicine to repair and replace damaged tissues as well as cells. Priyadarshini et al. [11] presented 3D printing of biocompatible materials using Multi Jet Fusion technology for application in bioreactors. They noted that when using the concept as a bioreactor for Multi Jet Fusion (MJF) printed materials, it is necessary to support (given the MJF technology) not only cell proliferation and growth but also other highly regulated and complex processes, as an example osteogenesis.

Ink as material. A guide to printing inks is provided by Leach et al. [12]. This guide contains all the information needed for the day-to-day development of printing inks and an understanding of this technique. The influence of ink types and printing processes on the efficiency of flotation deinking of recycled waste paper is presented by Nie et al. [13]. They highlight the need for research to better understand the surface physical and chemical reactions of the ink/fiber system, particularly during the printing process for office waste, in order to optimize the flotation process for this material’s efficient recycling. The weight change in the inkjet droplets for the two-phase flow was investigated by Turkan [14]. He performed Taguchi analysis and observed that the most effective parameter on ink weight was the inlet velocity. The fluid generation mechanism in a modified Continuous Inkjet Print (CIJ) method was investigated by Desai and Lovell [15]. Given the complex drop formation mechanism, a computational fluid dynamics (CFD) model is developed that is further validated using an ultrahigh-speed photography experimental setup. The influence of solvent volatility on the reproducibility of droplet formation in inkjet printing in the pharmaceutical industry was described by Mau and Seitz [16]. They noted that drop-on-demand (DoD) inkjet printing technology enables precise dosing and positioning of individual picoliter droplets.

Ink interaction. A study of the interaction of ink-paper using image analysis, volumetric inspection, and surface is presented by Fiadeiro et al. [17]. Two optical systems and their implementation were demonstrated, which were used to perform bulk and surface monitoring of ink-paper interactions through image analysis.

Future. Future challenges and opportunities of inkjet technologies are given by Castrejón-Pita et al. [18]. They noted that the technology needs consistent development to become widespread in new application areas such as additive manufacturing (3D printing).

This paper performs numerical simulation to analyze the different printhead horizontal moving velocities impact on the droplet form turnover during vertical fall, droplet shift, and droplet vertical movement velocities. Many known studies are mainly focused on the conditions of fixed printheads, Jiang and Tan [19], Cadarso et al. [20], and Wu et al. [21]. With this condition, the movement of the drop does not depend on the moving print head. In this work, the movement of the print head is taken into account; this allows us to analyze the movement of the droplets not only in the vertical direction but also in the horizontal.

2. Problem Formulation

Printing velocity plays an important role in ensuring print quality. This velocity is determined by the movement of the printing head, which has a direct effect on the behavior of the ink droplets during printing (spray and fall of a droplet). Printing velocity and droplet deposition precision are the main objectives that every printer manufacturer tries to reach. During the contemporary challenges to high-quality printing without losing time, it is necessary to investigate moving inkjet printhead models numerically. The main challenge is to ensure the accuracy of precipitation/deposition of droplets on the surface as a fundamental factor of print quality.

During a numerical experiment with COMSOL software (software version 6.2; developer location: 100 District Avenue, Burlington, MA 01803, USA), it is important to accurately reproduce the printing head parameters, its movements, and the printing velocity. This must be conducted in the appropriate three-dimensional space. Only by following appropriate steps in a numerical experiment can one expect the results of a droplet trajectory, which is needed to determine the numerical results close to reality (physically). On the one hand, 3D modeling gives the ability to ensure a more realistic experiment, but on the other hand, it requires large computer resources. One of the tasks of the work was to choose the appropriate option/parameter for the mesh refinement, which was needed to properly determine and set the air–liquid interface. This is the essential task of the printing head movement during a droplet spray. The numerical experiment requires a higher number of finite elements and computer resources to calculate/achieve a printing head motion using a considered/appropriate theoretical model. It was also important to model/simulate the droplet spray/fall/deposition itself, which is important because if the droplet is not properly formed, it can form aerosol-type droplets or not spray at all. Therefore, it was important to determine and establish the correct initial parameters during the numerical experiment, which are necessary to reproduce physical printed droplets (size, movement, and behavior). However, despite the practical significance of numerical results in this work, important droplet dynamical motion parameters such as falling velocity, droplet shift displacement, and moving trajectory will be presented. The next section will present the methodology of a numerical experiment.

3. Methodology

The single printhead ink chamber was modeled using COMSOL Multiphysics software (software version 6.2; developer location: 100 District Avenue, Burlington, MA 01803, USA). Inkjet printhead nozzle geometry dimensions are presented in Figure 1.

A droplet is generated by applying fluid intake velocity to the inlet cross-sectional ${\mathrm{r}}_1=0.1 \mathrm{m}\mathrm{m}$ to generate pressure inside the nozzle chamber (in Figure 1 presented in dark blue color) and push inks through the nozzle cross-section ${\mathrm{r}}_2=12 \mu \mathrm{m}$ to the air-filled chamber (in Figure 1 presented in light blue). The modeled air-filled cylindrical chamber has a radius $\mathrm{R}=0.1 \mathrm{m}\mathrm{m}$ and distance to surface $\mathrm{L}=1.1 \mathrm{m}\mathrm{m}.$

Constant pressure on the outlet has been set. The value of the pressure here is not directly important because the velocity depends only on the pressure gradient. On all other boundaries except the impact surface, there are set no-slip boundary conditions. The wetted wall boundary condition on the surface where the droplet falls off. Ink intake velocity and waveform dwelling time at the inlet are used to generate the pressure at the nozzle chamber. At the inlet, ink velocity increases from 0 to 0.5 m/s by the smooth parabolic profile actuated by an optimized waveform with a dwelling time of $2 \mathrm{\mu }\mathrm{s}.$ The time-dependent velocity profile is defined in seconds and is a very smooth rectangular pulse function with the transition points at 3 μs and 5 μs with a 2 μs dwelling period at the peak (Figure 2).

Additional ink is injected through the inlet during a period of 2 μs, and it consequently forces the ink to push out of the nozzle, forming the jet and transforming to the droplet furthermore. When the injection stops, a droplet of ink snaps off and continues to travel until it hits the surface.

The integration of the continuity equation and the Navier–Stokes equation in numerical simulations provides a coherent representation of incompressible substances two-phase flows. The formula used to calculate the continuity equation for an incompressible fluid is as follows:

(1)$\nabla ·\mathbf{u}={\displaystyle \frac{\partial {\mathbf{u}}_{\mathrm{x}}}{\partial \mathrm{x}}}+{\displaystyle \frac{\partial {\mathbf{u}}_{\mathrm{y}}}{\partial \mathrm{y}}}+{\displaystyle \frac{\partial {\mathbf{u}}_{\mathrm{z}}}{\partial \mathrm{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:

(2)$\nabla ={\displaystyle \frac{\partial }{\partial \mathrm{x}}}\mathbf{i}+{\displaystyle \frac{\partial }{\partial \mathrm{y}}}\mathbf{j}+{\displaystyle \frac{\partial }{\partial \mathrm{z}}}\mathbf{k},$

where i, j, and 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:

(3)$\mathsf{\rho }\left(\mathbf{u}·\nabla \right)\mathbf{u}=\mu \left(\nabla \left(\nabla ·\mathbf{u}+(\nabla ·{\mathbf{u}}^{\mathrm{T}}\right))\right)-\nabla \mathrm{p}+\mathsf{\rho }\mathrm{g}.$

Here, $\mathsf{\rho }$ is the substance density; $\mu$ is the substance viscosity; $\mathrm{p}$ is the pressure; $\mathrm{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 follows:

(4)${\displaystyle \frac{\partial \varphi }{\partial \mathrm{t}}}+\mathbf{u}·\nabla \varphi =\gamma \nabla ·\left(\epsilon \nabla \varphi -\varphi \left(1-\varphi \right){\displaystyle \frac{\nabla \varphi }{\left|\nabla \varphi \right|}}\right).$

Here, the coefficient $\varphi$ represents the level set interface between ink and air, according to the COMSOL Module User’s Guide [22]. Accordingly,$\varphi =0$ represents the air environment, and $\varphi =1$ represents ink. In the transition layer close to the interface, $\varphi$ smoothly changes from 0 to 1. Fluid velocity, $\mathbf{u}$ together represent the motion of the interface during the environment’s interplay. The parameter $\gamma$ determines the amount of stabilization of the level set function. A suitable value for $\gamma$ is the maximum magnitude occurring in the velocity field. Here, a time-dependent phase initialization approach is used. A numerical model for time-dependent two-phase flow requires an initialization of a level set function. By default, this reinitialization parameter $\gamma$ for the current simulation is set to $\gamma =10 \mathrm{m}/\mathrm{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 default and set to 1.25 μm.

The Two-Phase Flow, Level Set Multiphysics coupling feature defines the density and dynamic viscosity of the fluid used in the Laminar Flow or Turbulent Flow interface, and it defines the surface tension on the interface in the form of a volume force used in the momentum equation. It also enables the Level Set interface to use the velocity field calculated from the Laminar Flow interface to transport the interface, according to the COMSOL Module User’s Guide [22].

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 consider the dynamic viscosity (shear viscosity) $\mu$. It describes the interface between shear rate and shear stresses in the substrate. The substrate density and dynamic viscosity in different finite elements depend on the interface coefficient $\varphi$. Density and dynamic viscosity are written as follows, Olsson and Kreiss [23]:

(5)$\mathsf{\rho }={\mathsf{\rho }}_1+\left({\mathsf{\rho }}_2-{\mathsf{\rho }}_1\right)\varphi ;$

(6)$\mu ={\mu }_1+({\mu }_2-{\mu }_1)\varphi .$

In this instance, ${\mathsf{\rho }}_1$ represents the density of air, while ${\mathsf{\rho }}_2$ denotes the density of ink. Additionally, ${\mu }_1$ signifies the viscosity of air, and ${\mu }_2$ refers to the viscosity of ink.

The Navier–Stokes equations characterize the transfer of mass and momentum in an incompressible fluid, taking into account the effects of surface tension and gravitational forces. This research explores the ink and air substrates, analyzing how viscous stress develops at every point in the flow. The dynamics of the droplet, based on the conservation of mass, are derived from the subsequent differential equations, Sohr [24]:

(7)$(\nabla ·\mathbf{u})=0;$

(8)$\mathsf{\rho }\left({\displaystyle \frac{\partial \mathbf{u}}{\partial \mathrm{t}}}+\mathbf{u}·\nabla \mathbf{u}-\mathrm{g}\right)-\nabla ·\left(\mu \left(\nabla \mathbf{u}+\nabla {\mathbf{u}}^{\mathrm{T}}\right)\right)+\nabla \mathrm{p}={\mathbf{F}}_{\mathrm{s}\mathrm{t}}.$

The force of surface tension can be determined using the equation provided below [22]:

(9)${\mathbf{F}}_{\mathrm{s}\mathrm{t}}=\sigma \delta \mathrm{ĸ}\mathbf{n}.$

Here, δ is the Dirac delta function, σ is the surface tension coefficient, ĸ is the curvature, and n is the unit vector (normal direction). The calculation of the unit vector is derived from the equation below [22]:

(10)$\mathbf{n}={\displaystyle \frac{\nabla \varphi }{\left|\nabla \varphi \right|}}.$

The delta function is approximated by the following equation [22]:

(11)$\Delta =6\left|\varphi (1-\varphi )\right|\left|\nabla \varphi \right|.$

This functionality will partition the simulation progress into several time intervals, facilitating localized mesh refinement in the areas where the ink and air interface is observed during each interval, thus increasing the accuracy of the computations. At the outset, our model consists of around 220,000 finite elements, and the mesh is refined 12 times during the droplet movement at the air substrate. During the second mesh refinement, the total count of finite elements increases to 800,000. As the number of finite elements increases, the average quality of the elements declines to 0.64. The mass of the final droplet can be determined using the equation provided below [22]:

(12)${\mathrm{m}}_{\mathrm{d}}=\mathsf{\rho }\left(\mathrm{V}·\varphi \right).$

Several key dimensionless parameters significantly influence the working mechanism of the inkjet printhead. The Weber number (We) measures the ratio of inertia to surface tension, in contrast to the Reynolds number (Re), which serves to measure the ratio of a fluid’s inertia to its viscosity, Seipel et al. [25]. The calculation of these dimensionless numbers is performed using the following equations:

(13)$\mathrm{W}\mathrm{e}={\displaystyle \frac{\mathsf{\rho }{{\mathrm{d}\mathrm{v}}_{\mathrm{d}}}^2}{\sigma }}.$

(14)$\mathrm{R}\mathrm{e}={\displaystyle \frac{\mathsf{\rho }{\mathrm{d}\mathrm{v}}_{\mathrm{d}}}{\mu }};$

Here, $\mathrm{d}$ is the nozzle diameter, and ${\mathrm{v}}_{\mathrm{d}}$ is the droplet velocity.

The role of droplet velocity in the inkjet printing process can be predicted by the Ohnesorge number $\mathrm{O}\mathrm{h}$. A solution derived from the Navier–Stokes equations exists for articulating the limitations of droplet ejection, taking into account the interfacial, viscous, and inertial characteristics of the fluid [25]:

(15)$\mathrm{Z}={\displaystyle \frac{\sqrt{\sigma \mathsf{\rho }\mathrm{d}}}{\mu }}=\mathrm{}{\displaystyle \frac{1}{\mathrm{O}\mathrm{h}}}.$

The calculation of the Oh value is based on the physical properties of the inks and the dimensions of the printhead nozzle. The dynamics of an expanding droplet emitted from a printhead nozzle are closely impacted by Oh value.

It is necessary to establish an initial specification for ink printability to ensure the ejection of stable form droplets, which should be maintained within the acceptable limits of 1 < Z < 10. (Fromm [26]; Derby and Reis [27]). If Z is greater than 10, the printing operation will generate satellite droplets, which will diminish the overall quality of the print. When Z is less than 1, the fluid’s viscous forces will obstruct the expansion of the droplet. The results comprise empirical data that are contained within numerical simulation calculations, which rely on the dimensionless numbers We, Re, and Z.

4. Experiment and Results

An appropriate combination of physical parameters is necessary for the printing process in a DoD inkjet printhead to work properly, which will also depend on the size and velocity of the droplets. The physical characteristics of the ink and air utilized in the numerical simulation of the inkjet printhead nozzle are detailed in Table 1. Printing inks are described as Newtonian fluids for numerical calculations.

The form turnover in droplet ejection was achieved based on numerical simulations performed using COMSOL Multiphysics software. The presented graphs represent different horizontal printheads moving to the right side, with the different velocities (0 m/s, 1 m/s, 2 m/s, and 3 m/s) of the moving printhead shown in Figure 3. Here, μs represent time from the ejection process start. The velocity u increases to a maximum at the beginning of the ejection process because droplet diameter increases immediately while ejecting, and droplet fluid takes the same velocity as printhead moving velocity. After the droplet pinches off from the printhead, it moves forward to the surface and changes form during free falling. The parasite satellites are formed from the droplet tail while the droplets form the spherical form. Then, the printhead is not moving, and the satellites are deposited in the center of the droplet impact place. When printhead motion in the horizontal direction is added, the droplet tail formation is different.

To measure the droplet, the shift from the start ejecting position to the final form on the surface is shown in Figure 4. The investigation of droplet shifts by different printhead-moving velocities is represented. Simulation results represent how printhead moving velocity affects the deposition of the formatted droplet. The printhead horizontal move for 1 m/s impacts a 12 μm distance shift, 2 m/s–40 μm, 3 m/s–60 μm. All four ejected droplets reach the surface almost at the same time (±1 μs). There is visual graphic curvature turnover from 15 μs to 25 μs; at this period, the droplet jet tail unites with the main droplet mass and creates the sphere and some small malign satellites. Droplet vertical fall velocity also affects printing quality and droplet impact on the surface. Droplet movement velocity turnover during time is represented in Figure 5.

The investigation of droplet vertical movement velocity is not affected much by different printhead moving velocities. That proves that different printhead moving velocities affected droplets reaching the surface almost at the same time (32–34 μs).

5. Conclusions

In this study, we investigated the effects of different moving printhead velocities on droplet formation processes and droplet falling trajectories. A 3D numerical model has been investigated that allows a detailed study of the trajectories of a droplet falling on a surface when spraying ink with different spray parameters and at different velocities of movement of the print head.

Figure 3 shows the turnover in the form of a droplet during vertical fall. The tail progress of the droplet is observed. The conversion of droplet satellites and the curvature of the droplet tail can be investigated. The droplet shift increases due to the increase in the movement velocity of the print head.

Figure 4 shows the curvature of the droplet shift trajectory in more detail. The main changes in the shift are observed with a time step from 15 μs to 25 μs. During this period, the droplet tail shrinks to the main body of the droplet, and it forms a sphere. The droplet deposition shift increases after this formation.

Figure 5 shows a stable increase in the droplet falling velocity till the 15 μs time step. After this step, the droplet tail starts to tighten up, which reduces the droplet’s falling velocity. The velocity of movement of the print head does not affect the droplet velocity in the falling phase. It is observed that the droplet falling motion at different printhead movement velocities does not have a stable trajectory. It depends on many parameters, such as the physical parameters of the ink, the geometric parameters of the printhead, the droplet ejection parameters, and the print head movement velocity, but the main changes are related to the droplet formation phases during the falling process.

The main result is the analysis of the trajectory and velocity of the droplet at each time step in the case when the printing head moves horizontally. Certain shifts affect the quality of the print when the droplets begin the ejection process and reach the surface. This shift increased due to an increase in the velocity of the inkjet printhead. As a rule, an inkjet printhead prints and moves in two directions (moving forward and backward), and in each print cycle, droplets should reach the same initially predicted place. The main advantage is to identify the effect of a shift on the movement of droplets, the accuracy of such a movement, and the numerical investigation of such a dynamic process itself. Many studies are limited by the vertical movement of droplets. Our investigation considers not only the vertical (normal direction) movement of droplets but also the horizontal (tangential direction) movement, taking into account the velocity and movement of a horizontally moving known physical printing head.

Considering what is mentioned above, further analysis is still needed, in which various printhead actuation/ejection parameters will be investigated.

Author Contributions

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

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

DoD	drop-on-demand;
CIJ	Continuous Inkjet Print;
CFD	Computational fluid dynamics;
N-S	Navier–Stokes;
PZT	piezoelectric;
MJF	Multi Jet Fusion.
Greek
β	is the maximum spreading factor;
$\gamma$	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/m2);
$\mathsf{\rho }$	is the density (kg/m3);
$\sigma$	is the surface tension coefficient (mN/m);
ϕ	is the coefficient of level set interface between air and ink.
Roman
${\mathbf{F}}_{\mathrm{s}\mathrm{t}}$	is the surface tension force (N);
g	is the acceleration due to gravity (m/s2);
$\mathrm{L}$	is the air length between nozzle and surface (mm);
n	is the unit vector in the normal direction;
Oh	is the Ohnesorge number, a dimensionless coefficient;
p	is the pressure (Pa);
R	is the air chamber radius (mm);
${\mathrm{r}}_1$	is the inlet channel radius (mm);
${\mathrm{r}}_2$	$\mathrm{is} \mathrm{the} \mathrm{jetting} \mathrm{nozzle} \mathrm{radius} (\mu \mathrm{m})$;
Re	is the Reynolds number, a dimensionless coefficient;
t	is time (s);
u	is the fluid velocity (m/s);
${\mathrm{V}}_{\mathrm{i}\mathrm{n}}$	is the inlet ink velocity (m/s);
${\mathrm{v}}_{\mathrm{d}}$	is the droplet velocity at vertical direction (m/s).
We	is the Weber number, a dimensionless coefficient;
Z	is the ink printability range, a dimensionless number.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Table

Figure 1. Inkjet printhead nozzle geometry.

Figure 2. Fluid intake velocity and waveform dwelling time velocity intake to inlet.

Figure 3. Representation of the droplet form turnover during vertical fall. The presented graphs represent different horizontal printheads moving to the right side by different printhead moving velocities of 0 m/s (black color), 1 m/s (blue color), 2 m/s (green color), and 3 m/s (red color).

Figure 4. Representation of the droplet shift by different printhead moving velocities.

Figure 5. Representation of the droplet vertical movement velocity turnover during time.

Table 1

Physical parameters of the substances used in the numerical modelling of the inkjet printhead nozzle.

Substance	μmN·s/m2	ρkg/m3	Oh	σ mN/m	Z
Ink	14	1110	0.32	72	3.13
Air	1.789·10−2	1.225	-	-	-

References

1. Mohebi, M.M.; Evans, J.R.G. The trajectory of ink-jet droplets: Modelling and experiment. Chem. Eng. Sci.; 2005; 60, pp. 3469-3476. [DOI: https://dx.doi.org/10.1016/j.ces.2005.01.036]

2. Arango, I.; Bonil, L.; Posada, D.; Arcila, J. Prediction of a flying droplet landing over a non-flat substrate for ink-jet applications. Int. J. Interact. Des. Manuf. (IJIDeM); 2019; 13, pp. 967-980. [DOI: https://dx.doi.org/10.1007/s12008-019-00547-w]

3. Wang, S. Aerodynamic effect on inkjet main drop and satelite dot placement. Proceedings of the IS&Ts NIP 14: 1998 International Conference on Digital Printing Technologies; Toronto, ON, Canada, 18–23 October 1998; pp. 5-8.

4. Wu, W.Z.; Yang, X.; Zhang, B.Y.; Yin, Z.F.; Jia, B.Q. The study on electric field distribution and droplet trajectory during electrohydrodynamic jet printing. Microsyst. Technol.; 2021; 27, pp. 2745-2750. [DOI: https://dx.doi.org/10.1007/s00542-020-05040-z]

5. Link, N.; Lampert, S.; Gurka, R.; Liberzon, A.; Hetsroni, G.; Semiat, R. Ink drop motion in wide-format printers: II. Airflow investigation. Chem. Eng. Proc. Proc. Intensif.; 2009; 48, pp. 84-91. [DOI: https://dx.doi.org/10.1016/j.cep.2008.02.006]

6. Yang, A.S.; Cheng, C.H.; Lin, C.T. Investigation of droplet-ejection characteristics for a piezoelectric inkjet printhead. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci.; 2006; 220, pp. 435-445. [DOI: https://dx.doi.org/10.1243/09544062C04305]

7. Hudd, A. Inkjet printing technologies. The Chemistry of Inkjet Inks; Magdassi, S. World Scientific Publishing: Hackensack, NJ, USA, 2010; pp. 3-18.

8. Mashayekhi, M.; Conde, A.; Ng, T.N.; Mei, P.; Ramon, E.; Martinez-Domingo, C.; Alcalde, A.; Terés, L.; Carrabina Bordoll, J. Inkjet Printing Design Rules Formalization and Improvement. J. Disp. Technol.; 2015; 11, pp. 658-665. [DOI: https://dx.doi.org/10.1109/JDT.2015.2405135]

9. Zapka, W. Pros and Cons of Inkjet Technology in Industrial Inkjet Printing. Handbook of Industrial Inkjet Printing; Wiley: Weinheim, Germany, 2017; pp. 1-6.

10. Scoutaris, N.; Ross, S.; Douroumis, D. Current Trends on Medical and Pharmaceutical Applications of Inkjet Printing Technology. Pharm. Res.; 2016; 33, pp. 1799-1816. [DOI: https://dx.doi.org/10.1007/s11095-016-1931-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27174300]

11. Priyadarshini, B.M.; Kok, W.K.; Dikshit, V.; Feng, S.; Li, K.H.H.; Zhang, Y. 3D printing biocompatible materials with Multi Jet Fusion for bioreactor applications. Int. J. Bioprint.; 2023; 9, pp. 14-35. [DOI: https://dx.doi.org/10.18063/ijb.v9i1.623] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36636131]

12. Leach, R.H.; Pierce, R.J.; Hickman, E.P.; Mackenzie, M.J.; Smith, H.G. The Printing Ink Manual; 5th ed. Springer: Dordrecht, The Netherlands, 1993; [DOI: https://dx.doi.org/10.1007/978-1-4020-6187-5]

13. Nie, X.; Miller, J.D.; Yeboah, Y.D. The effect of ink types and printing processes on flotation deinking efficiency of wastepaper recycling. Environ. Eng. Policy; 1998; 1, pp. 47-58. [DOI: https://dx.doi.org/10.1007/s100220050005]

14. Türkan, B. Optimization analysis of the inkjet using the Taguchi method. Pamukkale Üniversitesi Mühendislik Bilim. Derg.; 2023; 29, pp. 281-288. [DOI: https://dx.doi.org/10.5505/pajes.2022.97957]

15. Desai, S.; Lovell, M. CFD analysis of a continuous inkjet print head for direct write fabrication. Proceedings of the ASME International Mechanical Engineering Congress and Exposition 2007; Seattle, WA, USA, 11–15 November 2007; Volume 43076, pp. 209-213. [DOI: https://dx.doi.org/10.1115/IMECE2007-43692]

16. Mau, R.; Seitz, H. Influence of the volatility of solvent on the reproducibility of droplet formation in pharmaceutical inkjet printing. Pharmaceutics; 2023; 15, 367. [DOI: https://dx.doi.org/10.3390/pharmaceutics15020367] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36839689]

17. 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. Proceedings of the SPIE 8785, 8th Iberoamerican Optics Meeting and 11th Latin American Meeting on Optics, Lasers, and Applications; Porto, Portugal, 22–26 July 2013; 8785BV. [DOI: https://dx.doi.org/10.1117/12.2024991]

18. Castrejón-Pita, J.R.; Baxter, W.R.S.; Morgan, J.; Temple, S.; Martin, G.D.; Hutchings, I.M. Future, opportunities and challenges of inkjet technologies. At. Sprays; 2013; 23, pp. 541-565. [DOI: https://dx.doi.org/10.1615/AtomizSpr.2013007653]

19. Jiang, H.; Tan, H. One dimensional model for droplet ejection process in inkjet devices. Fluids; 2018; 3, 28. [DOI: https://dx.doi.org/10.3390/fluids3020028]

20. Cadarso, V.J.; Perera-Nuñez, J.; Mendez-Vilas, A.; Labajos-Broncano, L.; González-Martín, M.-L.; Brugger, J. Microdrop generation and deposition of ionic liquids. J. Mater. Res.; 2014; 29, pp. 2100-2107. [DOI: https://dx.doi.org/10.1557/jmr.2014.162]

21. Wu, H.-C.; Shan, T.-R.; Hwang, W.-S.; Lin, H.-J. Study of micro-droplet behavior for a piezoelectric inkjet printing device using a single pulse voltage pattern. Mater. Trans.; 2004; 45, pp. 1794-1801. [DOI: https://dx.doi.org/10.2320/matertrans.45.1794]

22. CFD Module User’s Guide; COMSOL: Burlington, MA, USA, 2018; Available online: https://doc.comsol.com/5.4/doc/com.comsol.help.cfd/CFDModuleUsersGuide.pdf (accessed on 31 October 2024).

23. Olsson, E.; Kreiss, G. A conservative level set method for two phase flow. J. Comput. Phys.; 2005; 21, pp. 225-246. [DOI: https://dx.doi.org/10.1016/j.jcp.2005.04.007]

24. Sohr, H. The Navier-Stokes Equations: An Elementary Functional Analytic Approach; Springer Science & Business Media: Berlin, Germany, 2012; 367. [DOI: https://dx.doi.org/10.1007/978-3-0348-8255-2]

25. 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. [DOI: https://dx.doi.org/10.1038/s41598-020-75449-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33139808]

26. Fromm, J.E. Numerical calculation of the fluid dynamics of drop-on-demand jets. IBM J. Res. Dev.; 1984; 28, pp. 322-333. [DOI: https://dx.doi.org/10.1147/rd.283.0322]

27. Derby, B.; Reis, N. Inkjet printing of highly loaded particulate suspensions. MRS Bull.; 2003; 28, pp. 815-818. [DOI: https://dx.doi.org/10.1557/mrs2003.230]