Fluent Tutorial Guide Printer Ink
Fluent Tutorial Guide Printer Ink
17.1. Introduction
This tutorial examines the flow of ink as it is ejected from the nozzle of a printhead in an inkjet printer.
Using ANSYS Fluent’s volume of fluid (VOF) multiphase modeling capability, you will be able to predict
the shape and motion of the resulting droplets in an air chamber.
• Set up and solve a transient problem using the pressure-based solver and VOF model.
• Examine the flow and interface of the two fluids using volume fraction contours.
17.2. Prerequisites
This tutorial is written with the assumption that you have completed one or more of the introductory
tutorials found in this manual:
• Introduction to Using ANSYS Fluent in ANSYS Workbench: Fluid Flow and Heat Transfer in a Mixing El-
bow (p. 1)
• Introduction to Using ANSYS Fluent: Fluid Flow and Heat Transfer in a Mixing Elbow (p. 121)
and that you are familiar with the ANSYS Fluent tree and ribbon structure. Some steps in the setup and
solution procedure will not be shown explicitly.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 715
Using the VOF Model
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
716 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
• At time zero, the nozzle is filled with ink, while the rest of the domain is filled with air. Both fluids are assumed
to be at rest. To initiate the ejection, the ink velocity at the inlet boundary (which is modeled in this simulation
by a user-defined function) suddenly increases from 0 to 3.58 m/s and then decreases according to a cosine
law.
The calculation is run for 30 microseconds overall, that is, three times longer than the duration of the
initial impulse.
Because the dimensions are small, the double-precision version of ANSYS Fluent will be used. Air will
be designated as the primary phase, and ink (which will be modeled with the properties of liquid water)
will be designated as the secondary phase. Patching will be required to fill the ink chamber with the
secondary phase. Gravity will not be included in the simulation. To capture the capillary effect of the
ejected ink, the surface tension and prescription of the wetting angle will be specified. The surface inside
the nozzle will be modeled as neutrally wettable, while the surface surrounding the nozzle orifice will
be non-wettable.
17.4.1. Preparation
To prepare for running this tutorial:
Note
If you do not have a login, you can request one by clicking Customer Registration on
the log in page.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 717
Using the VOF Model
4. Narrow the results by using the filter on the left side of the page.
The files inkjet.msh and inlet1.c can be found in the vof directory created on unzipping the
file.
Fluent Launcher displays your Display Options preferences from the previous session.
9. Ensure that the Display Mesh After Reading and Workbench Color Scheme options are enabled.
For more information about the Fluent Launcher, see starting ANSYS Fluent using the Fluent Launcher
in the Fluent Getting Started Guide.
Note
The double precision solver is recommended for modeling multiphase flows simulation.
A warning message will be displayed twice in the console. You need not take any action at this point,
as the issue will be resolved when you define the solver settings in General Settings (p. 723).
2. Examine the mesh (Figure 17.2: Default Display of the Nozzle Mesh (p. 719)).
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
718 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Tip
By zooming in with the middle mouse button, you can see that the interior of the model
is composed of a fine mesh of quadrilateral cells (see Figure 17.3: The Quadrilateral
Mesh (p. 720)).
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 719
Using the VOF Model
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
720 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Ensure that All is selected from the Animation Option drop-down list.
Selecting All will allow you to see the movement of the entire mesh as you manipulate the Camera
view in the next step.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 721
Using the VOF Model
c. Click Apply.
d. Click the Camera... button to open the Camera Parameters dialog box.
Note
You may notice that the scale of the dimensions in the Camera Parameters dialog
box appear very large given the problem dimensions. This is because you have not
yet scaled the mesh to the correct units. You will do this in a later step.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
722 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Drag the indicator of the dial with the left mouse button in the clockwise direction until the upright
view is displayed (Figure 17.4: Mesh Display of the Nozzle Mirrored and Upright (p. 723)).
ANSYS Fluent will perform various checks on the mesh and report the progress in the console. Make sure
that the reported minimum volume is a positive number.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 723
Using the VOF Model
Note
It is a good idea to check the mesh after you manipulate it (that is, scale, convert to
polyhedra, merge, separate, fuse, add zones, or smooth and swap.) This will ensure that
the quality of the mesh has not been compromised.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
724 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
5. Retain the default setting of Pressure-Based in the Solver group box of the Setting Up Physics ribbon
tab.
7. Select Axisymmetric from the drop-down list in the Solver group box (below Reference Values).
17.4.4. Models
1. Enable the Volume of Fluid multiphase model.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 725
Using the VOF Model
b. Retain the default settings and click OK to close the Multiphase Model dialog box.
17.4.5. Materials
The default properties of air and water defined in ANSYS Fluent are suitable for this problem. In this step,
you will make sure that both materials are available for selection in later steps.
1. Add water to the list of fluid materials by copying it from the ANSYS Fluent materials database.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
726 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Click Fluent Database... in the Create/Edit Materials dialog box to open the Fluent Database Mater-
ials dialog box.
i. Select water-liquid (h2o < l >) from the Fluent Fluid Materials selection list.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 727
Using the VOF Model
Scroll down the Fluent Fluid Materials list to locate water-liquid (h2o < l >).
ii. Click Copy to copy the information for water to your list of fluid materials.
17.4.6. Phases
In the following steps, you will define water as the secondary phase. When you define the initial solution,
you will patch water in the nozzle region. In general, you can specify the primary and secondary phases
whichever way you prefer. It is a good idea to consider how your choice will affect the ease of problem setup,
especially with more complicated problems.
In the Phases dialog box, select phase 1 – Primary Phase and click Edit... to open the Primary
Phase dialog box.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
728 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
b. Retain the default selection of air in the Phase Material drop-down list.
In the Phases dialog box, select phase 2 – Secondary Phase and click Edit... to open the Secondary
Phase dialog box.
In the Phases dialog box, click Interaction... to open the Phase Interaction dialog box.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 729
Using the VOF Model
The surface tension inputs is displayed and the Continuum Surface Force model is set as the default.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
730 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
You will set the Reference Pressure Location to be a point where the fluid will always be 100 air.
If the UDF source file is not in your working directory, then you must enter the entire directory path
for Source File Name instead of just entering the file name. Alternatively, click the Browse... button
and select inlet1.c in the vof directory that was created after you unzipped the original file.
b. Click Interpret.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 731
Using the VOF Model
The UDF defined in inlet1.c is now visible and available for selection as udf membrane_speed in the
drop-down lists of relevant graphical user interface dialog boxes.
2. Set the boundary conditions at the inlet (inlet) for the secondary phase by selecting water-liquid from
the Phase drop-down list in the Boundary Conditions task page.
a. Click the Multiphase tab and enter 1 for the Volume Fraction.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
732 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
3. Set the boundary conditions at the outlet (outlet) for the secondary phase by selecting water-liquid from
the Phase drop-down list in the Boundary Conditions task page.
a. Click the Multiphase tab and retain the default setting of 0 for the Backflow Volume Fraction.
4. Set the conditions at the top wall of the air chamber (wall_no_wet) for the mixture by selecting mixture
from the Phase drop-down list in the Boundary Conditions task page.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 733
Using the VOF Model
Note
This angle affects the dynamics of droplet formation. You can repeat this simulation
to find out how the result changes when the wall is hydrophilic (that is, using a small
contact angle, say 10 degrees).
5. Set the conditions at the side wall of the ink chamber (wall_wet) for the mixture.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
734 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
17.4.10. Solution
1. Set the solution methods.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 735
Using the VOF Model
The non-iterative time advancement (NITA) scheme is often advantageous compared to the iterative
schemes as it is less CPU intensive. Although smaller time steps must be used with NITA compared
to the iterative schemes, the total CPU expense is often smaller. If the NITA scheme leads to convergence
difficulties, then the iterative schemes (for example, PISO, SIMPLE) should be used instead.
b. Select Fractional Step from the Scheme drop-down list in the Pressure-Velocity Coupling group
box.
c. Retain the default selection of Least Squares Cell Based from the Gradient drop-down list in the
Spatial Discretization group box.
d. Retain the default selection of PRESTO! from the Pressure drop-down list.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
736 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 737
Using the VOF Model
a. Retain the default settings for all the parameters and click Initialize (either in the ribbon or in the
Solution Initialization task page.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
738 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Retain the default setting of 0 mm for X Min and Y Min in the Input Coordinates group box.
d. Click Mark.
ANSYS Fluent will report in the console that 1500 cells were marked for refinement while zero cells
were marked for coarsening.
Tip
You can display and manipulate adaption registers, which are generated using the
Mark command, using the Manage Adaption Registers dialog box. Click the Man-
age... button in the Region Adaption dialog box to open the Manage Adaption
Registers dialog box. For details, see Adapting the Mesh (p. 174).
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 739
Using the VOF Model
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
740 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
b. Ensure that time-step is selected from the Append File Name with drop-down list.
ANSYS Fluent will append the time step value to the file name prefix (inkjet). The standard .dat
extension will also be appended. This will yield file names of the form inkjet-1-00200.dat,
where 200 is the time step number.
Optionally, you can add the extension .gz to the end of the file name (for example, inkjet.gz),
which instructs ANSYS Fluent to save the data files in a compressed format, yielding file names of the
form inkjet-1-00200.dat.gz.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 741
Using the VOF Model
Note
Small time steps are required to capture the oscillation of the droplet interface and
the associated high velocities. Failure to use sufficiently small time steps may cause
differences in the results between platforms.
c. Click Calculate.
17.4.11. Postprocessing
1. Read the data file for the solution after 6 microseconds (inkjet-1-00600.dat.gz).
2. Create and display a filled contour of water volume fraction after 6 microseconds (Figure 17.5: Contours
of Water Volume Fraction After 6 μs (p. 744)).
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
742 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
c. Select Phases... and Volume fraction from the Contours of drop-down lists.
e. Click Save/Display.
Tip
In order to display the contour plot in the graphics window, you may need to click the
3. Display contours of water volume fraction after 12, 18, 24, and 30 microseconds (Figure 17.6: Contours of
Water Volume Fraction After 12 μs (p. 744) — Figure 17.9: Contours of Water Volume Fraction After 30
μs (p. 746)).
a. Read the data file for the solution after 12 microseconds (inkjet-1-01200.dat.gz).
c. Repeat these steps for the 18, 24, and 30 microseconds files.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 743
Using the VOF Model
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
744 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 745
Using the VOF Model
17.5. Summary
This tutorial demonstrated the application of the volume of fluid method with surface tension effects.
The problem involved the 2D axisymmetric modeling of a transient liquid-gas interface, and postpro-
cessing showed how the position and shape of the surface between the two immiscible fluids changed
over time.
For additional details about VOF multiphase flow modeling, see the Fluent Theory Guide.
Release 18.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
746 of ANSYS, Inc. and its subsidiaries and affiliates.