Ansys Polyflow Tutorial Guide 2024 R1
Ansys Polyflow Tutorial Guide 2024 R1
ANSYS, Ansys Workbench, AUTODYN, CFX, FLUENT and any and all ANSYS, Inc. brand, product, service and feature
names, logos and slogans are registered trademarks or trademarks of ANSYS, Inc. or its subsidiaries located in the
United States or other countries. ICEM CFD is a trademark used by ANSYS, Inc. under license. CFX is a trademark
of Sony Corporation in Japan. All other brand, product, service and feature names or trademarks are the property
of their respective owners. FLEXlm and FLEXnet are trademarks of Flexera Software LLC.
Disclaimer Notice
THIS ANSYS SOFTWARE PRODUCT AND PROGRAM DOCUMENTATION INCLUDE TRADE SECRETS AND ARE CONFID-
ENTIAL AND PROPRIETARY PRODUCTS OF ANSYS, INC., ITS SUBSIDIARIES, OR LICENSORS. The software products
and documentation are furnished by ANSYS, Inc., its subsidiaries, or affiliates under a software license agreement
that contains provisions concerning non-disclosure, copying, length and nature of use, compliance with exporting
laws, warranties, disclaimers, limitations of liability, and remedies, and other provisions. The software products
and documentation may be used, disclosed, transferred, or copied only in accordance with the terms and conditions
of that software license agreement.
ANSYS, Inc. and ANSYS Europe, Ltd. are UL registered ISO 9001: 2015 companies.
For U.S. Government users, except as specifically granted by the ANSYS, Inc. software license agreement, the use,
duplication, or disclosure by the United States Government is subject to restrictions stated in the ANSYS, Inc.
software license agreement and FAR 12.212 (for non-DOD licenses).
Third-Party Software
See the legal information in the product help files for the complete Legal Notice for ANSYS proprietary software
and third-party software. If you are unable to access the Legal Notice, contact ANSYS, Inc.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. iii
Tutorial Guide
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
iv of ANSYS, Inc. and its subsidiaries and affiliates.
Tutorial Guide
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. v
Tutorial Guide
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
vi of ANSYS, Inc. and its subsidiaries and affiliates.
Tutorial Guide
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. vii
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
viii of ANSYS, Inc. and its subsidiaries and affiliates.
List of Figures
1.1. Problem Description ............................................................................................................................. 18
1.2. The Boundary Set for the Problem ......................................................................................................... 19
1.3. Selecting the Fluid Flow (Polyflow) Analysis System in Ansys Workbench ............................................... 21
1.4. Ansys Workbench with a New Ansys Polyflow-Based Fluid Flow Analysis System .................................... 22
1.5. Displaying the Files View after Adding an Ansys Polyflow-Based Fluid Flow Analysis System ................... 24
1.6. The Imported Geometry in the Ansys DesignModeler Application .......................................................... 26
1.7. The Ansys Meshing Application with the Extrusion Geometry Loaded .................................................... 27
1.8. Rotated View ........................................................................................................................................ 28
1.9. Selecting the Inlet Face ......................................................................................................................... 29
1.10. The Zero Velocity Faces Selected ......................................................................................................... 30
1.11. The Free Surface Faces Selected .......................................................................................................... 31
1.12. Rotated View ...................................................................................................................................... 32
1.13. The First Pair of Symmetry Faces Selected ............................................................................................ 33
1.14. The Second Pair of Symmetry Faces Selected ....................................................................................... 34
1.15. The Flow Exit Face Selected ................................................................................................................. 35
1.16. The Computational Mesh for the Extrusion Geometry .......................................................................... 38
1.17. The Ansys Polydata Application ........................................................................................................... 41
1.18. The Extrusion Geometry Loaded into Ansys CFD-Post .......................................................................... 46
1.19. Rotating the View ............................................................................................................................... 47
1.20. Contours of Velocity Magnitude .......................................................................................................... 49
1.21. Velocity Profiles at cross-sections ......................................................................................................... 54
1.22. Creating an Expression for an Output Parameter .................................................................................. 59
1.23. Displaying the Files View after Viewing Results in Ansys CFD-Post ........................................................ 60
1.24. The Parameters and Design Points View ............................................................................................... 62
1.25.The Chart of the Current Design Point .................................................................................................. 64
1.26. Displaying Values for All of the Design Points ....................................................................................... 65
1.27. The Chart of All of the Design Points .................................................................................................... 65
1.28. Displaying the Files View after Exploring Solutions ............................................................................... 66
1.1. Problem Schematic ............................................................................................................................... 70
1.2. Subdomains and Boundary Sets for the Problem ................................................................................... 71
1.3. Contours of Static Pressure .................................................................................................................... 95
1.4. Contours of Static Pressure after Applying Rotation ............................................................................... 98
1.5. Velocity Vectors ................................................................................................................................... 100
1.6. Velocity Vectors with Mesh .................................................................................................................. 104
2.1. A Schematic Diagram of the Fluid and the Circular Die ......................................................................... 106
2.2. Boundaries and Sub-domains .............................................................................................................. 107
2.3. Pressure Contours ............................................................................................................................... 123
2.4. Velocity Distribution ............................................................................................................................ 124
2.5. Velocity Vectors ................................................................................................................................... 125
2.6. Temperature Distribution (Celsius) ....................................................................................................... 126
2.7. Temperature Distribution (Kelvin) ........................................................................................................ 127
2.8. Temperature Profile Along the Line Y = 0.006 m ................................................................................... 130
3.1. Problem Description ........................................................................................................................... 132
3.2. Boundaries and Subdomains ............................................................................................................... 134
3.3. Contours of Pressure ........................................................................................................................... 154
3.4. Velocity Profile .................................................................................................................................... 155
3.5. Velocity Vectors ................................................................................................................................... 157
3.6. Temperature Profile ............................................................................................................................. 159
3.7. Temperature Profile Near the Die Exit .................................................................................................. 160
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. ix
Tutorial Guide
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
x of ANSYS, Inc. and its subsidiaries and affiliates.
Tutorial Guide
3.1. Plug-Assisted Thermoforming of a Blister: Plug (Orange), Sheet (Blue), and Mold (Green) in the Initial
Configuration ........................................................................................................................................... 354
3.2. Finite Element Mesh, Subdomains and Boundary Sets .......................................................................... 355
3.3. Ramp Function for Mold Velocity ......................................................................................................... 361
3.4. Multi-Ramp Function for Plug Velocity ................................................................................................. 365
3.5. Ramp Function for Pressure ................................................................................................................. 368
3.6. Contours of Thickness 1/4 Geometry ................................................................................................... 378
3.7. Contours of Thickness on the Whole Thermoformed Blister .................................................................. 380
3.8. Thickness of the Film at Time = 0.0001 s ............................................................................................... 382
3.9. Contours of Thickness at Time = 0.1208 s ............................................................................................. 383
3.10. Contours of Thickness at Time = 0.1291 s ........................................................................................... 384
3.11. Contours of Thickness at Time = 1.0 s ................................................................................................. 385
4.1. Blow Molding Initial Configuration ...................................................................................................... 390
4.2. Finite Element Mesh, Subdomains, and Boundary Sets ......................................................................... 391
4.3. Ramp Function for Right Mold Velocity ................................................................................................ 397
4.4. Ramp Function for Left Mold Velocity .................................................................................................. 400
4.5. Ramp Function for Pressure ................................................................................................................. 404
4.6. Contours of Thickness on the Parison .................................................................................................. 413
4.7. Thickness of the Film at Time = 0.01 s .................................................................................................. 415
4.8. Contours of Thickness at Time = 0.09253 s ........................................................................................... 416
4.9. Contours of Thickness at Time = 0.23409 s ........................................................................................... 417
4.10. Contours of Thickness at Time = 2.0 s ................................................................................................. 418
4.11. Effect of Adaption on Final Mesh and Thickness Variation ................................................................... 421
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. xi
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
xii of ANSYS, Inc. and its subsidiaries and affiliates.
Using This Manual
1. The Contents of This Manual
The Polyflow Tutorial Guide contains a number of example problems with complete detailed instructions,
commentary, and postprocessing of results.
• The Polyflow User's Guide explains how to use Ansys Polydata and Ansys Polyflow to set up and
solve a problem.
• The Polyflow Tutorial Guide contains a number of example problems with complete detailed
instructions, commentary, and postprocessing of results.
Tutorials for release 2024 R1 are available on the Ansys Help site.
• The Polyflow in Workbench User's Guide explains how to use the Ansys Polyflow application
within Ansys Workbench.
• The Polymat User's Guide explains how to use the Ansys Polymat module for material property
evaluation.
• The Polystat User's Guide explains how to set up a MIXING task in Ansys Polydata and how to
use the Ansys Polystat module for statistical postprocessing of results.
• The Polyflow Examples Manual provides overviews of solutions to a variety of problem types
and is available on the Ansys Help site.
• The CFD-Post User's Guide explains how to use CFD-Post to examine your results.
For details on how to access the Ansys Polyflow manuals, see the Ansys Polyflow documentation section
of the separate Polyflow User's Guide.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. xiii
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
xiv of ANSYS, Inc. and its subsidiaries and affiliates.
Part 1: Introduction to Using Polyflow in Workbench
The following Workbench tutorial is available:
• Create a computational mesh for the geometry using the Ansys Meshing application.
• Import a data file, and modify it using Ansys Polydata to include a user-defined template for the die
inlet flow rate.
• View the initial results and create an output parameter for the maximum velocity of the extrudate in
CFD-Post.
• Generate results for multiple design points using the Parameter and Design Points view, and chart
how the outflow velocity varies with the inlet flow rate.
1.2. Prerequisites
This tutorial assumes that you have little to no experience with Ansys DesignModeler, Ansys Meshing,
Ansys Polyflow, CFD-Post, or the Parameter and Design Points view of Ansys Workbench, and so each
step will be explicitly described.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 17
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
The melt enters the die as shown in Figure 1.1: Problem Description (p. 18) at an initial flow rate of
cm3/s (this flow rate is a quarter of that for the complete physical system) and the extrudate is
obtained at the exit. It is assumed that the extrudate is fully deformed at the end of the computational
domain, and that it will not deform any further (that is, subdomain 2 is long enough to account for all
the deformation of the extrudate).
The incompressibility and momentum equations are solved over the computational domain. The domain
for the problem is divided into two subdomains (as shown in Figure 1.1: Problem Description (p. 18))
so that a remeshing algorithm can be applied only to the portion of the mesh that will be deformed.
Subdomain 1 represents the die where the fluid is confined. Subdomain 2 corresponds to the extrudate
that is in contact with the air and can deform freely. The calculation will determine the location of the
free surface (the skin of the extrudate), as well as the velocity of the extrudate at the exit.
The boundary set for the problem is shown in Figure 1.2: The Boundary Set for the Problem (p. 19), and
the conditions at the boundaries of the domains are:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
18 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the ext3_workbench.zip file you have downloaded to your working folder.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 19
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
The extracted files include the geometry file ext3d.x_t , the data file polyflow.dat , and a
solution folder that contains the solution files created during the preparation of the tutorial.
Note:
This tutorial is prepared using Ansys Polyflow on a Windows system. The screen shots
and graphic images that follow may be slightly different than the appearance on your
system, depending on the operating system or graphics card.
The Ansys Workbench application window will open, containing the Toolbox on the left and the
Project Schematic on the right. The Toolbox lists the various supported analyses and applications,
and the Project Schematic provides a space to display the components of the analysis systems you
select.
Note:
When you first start Ansys Workbench, the Getting Started message window is dis-
played, offering assistance through the online help for using the application. You can
keep the window open, or close it by clicking the red ‘X’ icon in the upper right corner.
If you need to access the online help at any time, use the Help menu, or press the
F1 key.
2. Create a new fluid flow analysis system by double-clicking the Fluid Flow (Polyflow) option
under Analysis Systems in the Toolbox.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
20 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Figure 1.3: Selecting the Fluid Flow (Polyflow) Analysis System in Ansys Workbench
Extra:
You can also create a new fluid flow analysis system by dragging-and-dropping the
analysis system into the Project Schematic: a green dotted outline will indicate a po-
tential location in the Project Schematic for the new system, which will turn into a
red box when you attempt to drop it.
A new Ansys Polyflow-based fluid flow analysis system will be displayed in the Project Schematic.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 21
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
Figure 1.4: Ansys Workbench with a New Ansys Polyflow-Based Fluid Flow Analysis System
Note:
The Ansys Polyflow-based fluid flow analysis system is composed of various cells
(Geometry, Mesh, and so on) that represent the work flow for performing the analysis.
Ansys Workbench is composed of multiple data-integrated (for example, Ansys Polyflow)
and native applications into a single, seamless project flow, where individual cells can
obtain data from and provide data to other cells. Ansys Workbench provides visual
indications of a cell’s state at any given time via icons on the right side of each cell.
Because of the constant flow of data, a cell’s state can quickly change. Brief descriptions
of the various states are provided below. For more information about cell states, see
the Ansys Workbench online help.
• Unfulfilled ( ) indicates that required upstream data does not exist. For example, when you
first create a new Fluid Flow (Polyflow) analysis system, all cells downstream of the Geometry
cell appear as Unfulfilled because you have not yet specified a geometry for the system.
• Refresh Required ( ) indicates that upstream data has changed since the last refresh or up-
date. For example, after you assign a geometry to the Geometry cell in your new Fluid Flow
(Polyflow) analysis system, the Mesh cell appears as Refresh Required since the geometry
data has not yet been passed from the Geometry cell to the Mesh cell.
• Attention Required ( ) indicates that the current upstream data has been passed to the cell,
however, you must take some action to proceed. For example, after you launch Ansys Poly-
data from the Setup cell in a Fluid Flow (Polyflow) analysis system that has a valid mesh, the
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
22 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Setup cell appears as Attention Required because additional data must be entered in Ansys
Polydata before you can calculate a solution.
• Update Required ( ) indicates that local data has changed and the output of the cell must
be regenerated. For example, after you launch Ansys Meshing from the Mesh cell in a Fluid
Flow (Polyflow) analysis system that has a valid geometry, the Mesh cell appears as Update
Required because the Mesh cell has all the data it requires to generate an Ansys Polyflow mesh
file, but the Ansys Polyflow mesh file has not yet been generated.
• Up-to-Date ( ) indicates that an update has been performed on the cell and no failures have
occurred (or an interactive calculation has been completed successfully). For example, after
Ansys Polyflow finishes performing the number of necessary solver iterations, the Solution
cell appears as Up-to-Date.
• Interrupted ( ) indicates that you have interrupted an update (or stopped an interactive
calculation that is in progress). For example, if you select the stop button ( ) in the Progress
Monitor of Ansys Workbench at a point where Ansys Polyflow has generated results but has
not yet completed the calculation (such as during a transient simulation), then verify the action
in the dialog box that opens, Ansys Polyflow is immediately stopped and the Solution cell
appears as Interrupted.
• Input Changes Pending ( ) indicates that the cell is locally up-to-date, but may change
when next updated as a result of changes made to upstream cells. For example, if you change
the Mesh in an Up-to-Date Fluid Flow (Polyflow) analysis system, the Setup cell appears as
Refresh Required, and the Solution and Results cells appear as Input Changes Pending.
• Pending ( ) indicates that a batch or asynchronous solution is in progress. This icon will only
appear when the Solution cell is in background mode.
• Refresh Failed, Refresh Required ( ) indicates that the last attempt to refresh cell input
data failed, and so the cell must be refreshed.
• Update Failed, Update Required ( ) indicates that the last attempt to update the cell and
calculate output data failed, and so the cell must be updated. For example, if you update the
Solution cell and the solver diverges during the calculation, the Solution cell appears as Update
Failed, Update Required.
• Update Failed, Attention Required ( ) indicates that the last attempt to update the cell
and calculate output data failed, and so the cell requires attention.
a. Double-click the Fluid Flow (Polyflow) label underneath the analysis system.
a. Select the Save option under the File menu in Ansys Workbench.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 23
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
File → Save
The Save As dialog will open, where you can browse to a specific directory and enter a specific
name for the Ansys Workbench project.
b. In your working directory, enter ext3d-wb as the project File name and click the Save
button to save the project. Ansys Workbench saves the project with a .wbpj extension, as
well as supporting files for the project.
5. View the files generated by Ansys Workbench, by enabling the Files option under the View
menu.
View → Files
Figure 1.5: Displaying the Files View after Adding an Ansys Polyflow-Based Fluid Flow
Analysis System
Ansys Workbench allows you to easily view the files associated with your project using the Files view.
You can see the name and type of file, the ID of the cell the file is associated with, the size of the file,
the location of the file, and other information. For more information about the Files view, see the
separate Polyflow in Workbench User's Guide and the Ansys Workbench online help. To access docu-
mentation files, go to the documentation area of the customer site.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
24 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Note:
Ansys DesignModeler is licensed separately from Ansys Polyflow. If you do not have access
to Ansys DesignModeler, you can instead import a geometry file that does not need to be
modified, as noted in step 1.c.
a. Right-click the Geometry cell in the ext3d fluid flow analysis system (cell A2 in the Ansys
Workbench Project Schematic).
b. Move your pointer over Import Geometry in the context menu that opens, and click
Browse....
c. Use the Open dialog box to browse to the folder you unzipped in a previous step, select
ext3d.x_t, and click Open.
Note:
If you do not have access to Ansys DesignModeler, select PFL.agdb in the Open
dialog box instead, then skip to Meshing the Geometry in the Ansys Meshing Ap-
plication (p. 27).
The state of the Geometry cell becomes Up-to-Date, indicating that there is a geometry now asso-
ciated with the fluid flow analysis system.
Right-click the Geometry cell in the ext3d fluid flow analysis system, and select Edit Geometry
in DesignModeler..., to launch the Ansys DesignModeler application.
3. Finish importing the geometry file by clicking Generate in the Ansys DesignModeler toolbar.
The geometry will be displayed in the Graphics window.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 25
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
Note that the scale at the bottom of the Graphics window can be used to confirm that the overall
length of the domain is 0.6 m.
4. Modify the geometry so that the separate domains ("bodies") are treated as a single entity (a
"part"), by performing the following actions in the Tree Outline.
By uniting the multiple bodies of the geometry into a single part, you will create a conformal mesh
between the separate domains of the bodies.
d. Right-click the highlighted objects and click Form New Part in the menu that opens.
You can simply close the Ansys DesignModeler application. Ansys Workbench automatically saves
the geometry and updates the Project Schematic accordingly.
6. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
26 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Note the addition of the geometry file (PFL.agdb, where PFL indicates a Polyflow-based fluid flow
system) to the list of files.
Double-click the Mesh cell in the ext3d fluid flow analysis system (cell A3) to launch the Ansys
Meshing application with the extrusion geometry already loaded.
Extra:
You can also right-click the Mesh cell to display the context menu where you can
select the Edit... option.
Figure 1.7: The Ansys Meshing Application with the Extrusion Geometry Loaded
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 27
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
2. Group the faces and create named selections to match the boundary set shown in Figure 1.2: The
Boundary Set for the Problem (p. 19).
a. Rotate the view to get your display similar to that shown in Figure 1.8: Rotated View (p. 28),
by holding the center mouse button and moving your pointer in the geometry window.
(You can also manipulate the view by clicking in the Ansys Meshing toolbar and dragging
the model).
Tip:
Look at the orientation of the axis triad, , to assist when rotating the
view.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
28 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Information will be displayed about the mesh in the Details view below the Outline tree view.
Note:
Note that since the Ansys Meshing application automatically detects that you are
going to perform a CFD fluid flow analysis, CFD is selected from the Physics
Preference drop-down list.
c. Select the face that will represent the inlet, as shown highlighted in green in Figure 1.9: Se-
lecting the Inlet Face (p. 29).
d. Right-click and select the Create Named Selection option (from the menu that opens) to
open the Selection Name dialog box.
e. Enter inlet for the name in the Selection Name dialog box, and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 29
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
f. Hold down the Ctrl key, select the 2 faces that will represent the zero velocity boundary (as
highlighted in green in Figure 1.10: The Zero Velocity Faces Selected (p. 30)), then create a
selection named die wall in a manner similar to the previous steps.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
30 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
g. Hold down the Ctrl key, select the 2 faces that will represent the free surface boundary (as
highlighted in green in Figure 1.11: The Free Surface Faces Selected (p. 31)), and create a
selection named free surface in a manner similar to the previous steps.
h. Rotate the view to get your display to be similar to that shown in Figure 1.12: Rotated
View (p. 32), by holding the center mouse button and moving your pointer in the geometry
window.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 31
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
i. Hold down the Ctrl key, select the 2 faces that will represent one of the symmetry bound-
aries (as highlighted in green in Figure 1.13: The First Pair of Symmetry Faces Selected (p. 33)),
and create a selection named symmetry 1 in a manner similar to the previous steps.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
32 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
j. Hold down the Ctrl key, select the 2 faces that will represent the other of the symmetry
boundaries (as highlighted in green in Figure 1.14: The Second Pair of Symmetry Faces Se-
lected (p. 34)), and create a selection named symmetry 2 in a manner similar to the pre-
vious steps.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 33
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
k. Select the face that will represent the flow exit boundary (as highlighted in green in Fig-
ure 1.15: The Flow Exit Face Selected (p. 35)), and create a selection named outlet in a
manner similar to the previous steps.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
34 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
3. Set the appropriate meshing parameters for the Ansys Meshing application in the Details view.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 35
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
a. Confirm that Polyflow is selected from the Solver Preference drop-down list under Defaults.
a. Right-click Mesh in the Outline tree view, and select Update in the context menu.
Note:
Using the Generate Mesh option from the Mesh context menu creates the mesh,
but does not actually create the relevant mesh files for the project and is optional
if you already know that the mesh is acceptable. Using the Update option auto-
matically generates the mesh and creates the relevant mesh files for your project
and updates the Ansys Workbench cell that references this mesh.
i. Enter 0.008m for Element Size under Defaults in the Details view.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
36 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Right-click Mesh in the Outline tree view, and select Update in the context menu.
Extra:
After the mesh is generated, you can view the mesh statistics by expanding the
Statistics node in the Details view to reveal information about the number of nodes,
the number of elements, and other details.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 37
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
When you close the Ansys Meshing application, Ansys Workbench automatically saves the mesh
and updates the Project Schematic accordingly (the state of the Mesh cell changes from Refresh
Required to Up-to-Date, indicating that there is a mesh now associated with the fluid flow
analysis system).
6. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
38 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Note the addition of the mesh files (PFL.1.poly and PFL.mshdb) to the list of files. The
PFL.1.poly file was created when you updated the mesh, and the PFL.mshdb file was generated
when you closed the Ansys Meshing application.
The data file you will import has already been set up for a 3D extrusion simulation with a single inlet
flow rate. For details on how to set up a similar data file in Ansys Polydata, see 3D Extrusion (p. 165).
a. Right-click the Setup cell in the ext3d fluid flow analysis system, and click Import Polyflow
Dat ... in the context menu that opens.
b. Use the Open dialog box to browse to the folder you unzipped in a previous step, select
polyflow.dat, and click Open.
The state of the Setup cell remains Refresh Required, indicating that even though there is a
data file now associated with the fluid flow analysis system, you still must perform an update
for the cell.
c. Right-click the Setup cell and click Update in the context menu that opens.
After Ansys Polydata checks for coherence between the mesh and data files, the state for the
Setup cell becomes Up-to-Date. At this point it would be possible to run the Ansys Polyflow solver
for your simulation; however, for this tutorial you will first modify the data file.
2. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 39
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
Note the addition of the data file (polyflow.dat) to the list of files.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
40 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Double-click the Setup cell in the ext3d fluid flow analysis system.
Extra:
You can also launch Ansys Polydata by right-clicking the Setup cell and clicking Edit...
in the context menu that opens.
Note:
The mesh is automatically loaded and displayed in the Graphics Display window by
default.
4. View the mesh information, in order to verify the unit that should be used for length when de-
fining your inputs for the simulation.
It is a good practice to always perform this step with new meshes. Polydata and Polyflow do not
consider units when calculating a solution, so it is your responsibility to ensure that you enter values
for the velocity, material data, and other settings that are consistent with each other and with the
mesh.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 41
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
b. Click the Info button to open a panel that displays the mesh information.
c. Verify that the Length Unit used to create the mesh was meters and that the dimensions
of the Bounding Box (which surrounds the mesh) are 0.1 x 0.1 x 0.6.
It is therefore recommended that you use meters for the length unit when specifying the inputs
for the simulation. Note that if you decided that you would rather work with a different length
unit, you could scale the mesh using Polyfuse, as described in the Polyflow User's Guide.
d. Close the panel and click the Menus tab at the bottom of the Polydata window.
5. Modify the data file so that the inlet flow rate is flagged as modifiable in a user-defined template
(UDT).
Note:
F.E.M. Task 1
3D die swell
3D die swell is the name that was given to the sub-task for the flow problem when the data
file was created.
d. Select Inflow along INLET in the Flow boundary conditions menu and click Modify.
e. Click the UPDT button at the top of the Ansys Polydata application window, to enable
template inputs.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
42 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
g. Retain the selections of Automatic and Volumetric flow rate in the Inflow calculation on
INLET menu, and note that the flow rate is already set to 1 x 10-5 (which, since you are using
meters for your length unit, is equal to 10 cm3/s—that is, the initial flow rate proposed in
Problem Description (p. 18)). Then click Upper level menu.
h. Click Create a new template entry in the Create template entry menu.
i. Click the UPDT button again at the top of the Ansys Polydata application window, to disable
template inputs.
j. Click Upper level menu four times to return to the main Polydata menu.
6. Verify the system of units that will be passed to CFD-Post for postprocessing.
It is a good practice to always perform this step before running a new simulation, to ensure consistency
with the mesh and the units you used when defining the velocity, material data, and other settings.
Outputs
a. Note that the Current output(s) in the Outputs menu indicate that CFD-Post is currently
selected as the intended postprocessor.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 43
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
c. Note that meter, kilogram, and second are currently selected for Length, Mass, and Time,
respectively (which is consistent with the values used in setting up the data file).
d. Click Upper level menu twice to return to the main Polydata menu.
A Parameters cell will be added to the ext3d fluid flow analysis system in the Ansys Workbench
Project Schematic (cell A7). Also, a Parameter Set bar will be added below the system with an
inbound arrow, indicating that an input parameter has been created.
8. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Note the addition of the template file (templates.upd) to the list of files.
1.4.6. Solution
1. Start Ansys Polyflow.
In the Ansys Workbench Project Schematic, right-click the Solution cell in the ext3d fluid flow
analysis system (cell A5), and click Update in the context menu that opens.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
44 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The Ansys Polyflow solver will begin running. When the calculation is complete, the state for the
Solution cell becomes Up-to-Date.
2. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Note the addition of the listing file (polyflow.lst), the Ansys Polyflow results file (res), the output
mesh file (res.msh), the CFD-Post file (cfx.res), and the automatically generated probe files
(.prb) to the list of files. For more information about Ansys Polyflow (and the files associated with
it), see files written and read by AnsysPolydata in the Polyflow User's Guide. To access documentation
files, go to the documentation area of the customer site.
1.4.7. Postprocessing
In the following steps you will use Ansys CFD-Post to view the results of your initial simulation, create
an expression that can be used as an output parameter for Ansys Workbench, then review the list of
files generated by Ansys Workbench.
In the Ansys Workbench Project Schematic, double-click the Results cell in the ext3d fluid
flow analysis system (cell A6).
Extra:
You can also start Ansys CFD-Post by right-clicking the Results cell and selecting the
Edit... option in the context menu that opens.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 45
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
The Ansys CFD-Post application will launch with the extrusion geometry already loaded (displayed in
outline mode). Note that Ansys Polyflow results are also automatically loaded into Ansys CFD-Post.
2. Obtain the view shown in Figure 1.19: Rotating the View (p. 47).
a. Rotate the view, by holding the center mouse button and moving your pointer in the
viewer area.
b. Reduce the magnification of the view by clicking the Zoom icon at the top of the viewer
area ( ), holding the left mouse button, and moving your pointer in the viewer area.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
46 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
3. Display contours of velocity magnitude on the boundaries (Figure 1.20: Contours of Velocity
Magnitude (p. 49)).
Insert → Contour
b. Retain the default entry of Contour 1 for Name and click OK to close the dialog box.
Information about Contour 1 will be displayed in the Details view below the Tree view in
Ansys CFD-Post. The Details view contains all of the settings for a contour object.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 47
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
c. Open the Location Selector dialog box by clicking the location editor button ( ) next to
the Locations drop-down list in the Geometry tab.
i. Select all of the boundaries listed under ext3d by clicking the first one in the list
(PART_1_1_SOL_DIE_WALL), holding the Shift key, and clicking the last one in the
list (PART_1_2_SOL_SYMMETRY_2).
e. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
48 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The velocity is 0 along the die wall (as expected) and there is a fully developed profile at the inlet of
the die. At the die outlet, the velocity profile changes to become constant throughout the extrudate
cross-section. The transition between these two states can be seen in the first third of the extrudate.
4. Display contours of velocity in cross-sections (Figure 1.21: Velocity Profiles at cross-sections (p. 54)).
a. Disable Contour 1 under User Locations and Plots in the Outline tab of the Tree view.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 49
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
i. Select Plane from the Location drop-down menu, located in the toolbar.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
50 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Retain the default entry of Plane 1 for Name in the Insert Plane dialog box that opens,
and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 51
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
iii. Retain the default selection of XY Plane for Method in the Geometry tab of the Details
view for Plane 1.
v. Click Apply
c. In a similar manner, create cross-section planes at =0.08 m, 0.15 m, and 0.45 m named
Plane 2, Plane 3, and Plane 4 respectively. Note that you will retain the default selection
of XY Plane for Method and enter appropriate values for Z in the Details view.
d. Disable Plane 1, Plane 2, Plane 3, and Plane 4 under User Locations and Plots in the
Outline tab of the Tree view, so that the planes are no longer colored gray in the viewer
area.
Insert → Contour
f. Retain the default entry of Contour 2 for Name and click OK to close the dialog box.
Information about Contour 2 will be displayed in the Details view below the Tree view.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
52 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
g. Open the Location Selector dialog box by clicking the location editor button ( ) next to
the Locations drop-down list in the Geometry tab.
i. Select all of the planes listed under User Locations and Plots by clicking Plane 1,
holding the Shift key, and clicking Plane 4.
i. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 53
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
Velocity profiles at the flow inlet, the flow outlet, and planes just before and just after the die exit are
displayed. Compare the velocity profile within the die to the velocity profile just after the die exit at
the end of the computational domain. In the die the flow is fully developed. The velocity profile is flat
(that is, all the particles in the cross-section are at the same velocity) in the extrudate, far away from
the die exit. In the transitional zone just beyond the die exit, the velocity profile is reorganized. The
velocity profile on the plane =0.15 m is no longer fully developed, but it is not yet flat either. The
velocity rearrangement is the source of the deformation of the extrudate.
5. Create an expression for the maximum velocity at the flow exit, which can be used as an output
parameter in Ansys Workbench.
b. Right-click anywhere in the Expressions tab and click New in the menu that opens to create
a new expression.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
54 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
c. Right-click in the Definition tab of the Details view, move your pointer over Functions,
move your pointer over CFD-Post, and click maxVal, to specify that the function in the ex-
pression obtains the maximum value.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 55
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
d. Make sure that the cursor is between the parentheses of maxVal()@, right-click in the Details
view again, move your pointer over Variables, and click VELOCITIES, to specify that the
variables obtained in the expression are velocities.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
56 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
e. Move the cursor so that it is after the @ symbol of maxVal (VELOCITIES)@, right-click in the
Details view again, move your pointer over Locations, and click PART_1_2_SOL_OUTLET,
to specify that the variables are obtained for the expression at the flow exit.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 57
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
f. Click Apply.
The expression in the Definition tab of the Details view will be defined as maxVal (VELOCITIES)@
PART_1_2_SOL_OUTLET with a Value of approximately 7.8 x 10-4m/s, and maxvelocity will
be added to the list in the Expressions tab of the Tree view, as shown in Figure 1.22: Creating
an Expression for an Output Parameter (p. 59).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
58 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
g. Right-click maxvelocity in the Expressions tab of the Tree view and select Use as Work-
bench Output Parameter in the context menu that opens.
An outbound arrow will be added from the Parameters cell to the Parameter Set bar in
the Project Schematic , indicating that an output parameter has been created.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 59
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
Note:
Note that the Ansys CFD-Post state files are automatically saved when you exit Ansys
CFD-Post and return to Ansys Workbench.
File → Save
8. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Figure 1.23: Displaying the Files View after Viewing Results in Ansys CFD-Post
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
60 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Note the addition of the Ansys CFD-Post state file (ext3d.cst) to the list of files. For more inform-
ation about Ansys CFD-Post (and the files associated with it), see the Ansys CFD-Post documentation.
Note:
Ansys DesignXplorer is licensed separately from Ansys Polyflow. If you do not have access
to Ansys DesignXplorer, you will not be able to perform some of the steps that follow,
such as computing multiple design points or plotting results in a chart.
1. Open the Parameters Set tab, which contains the Parameters and Design Points view (Fig-
ure 1.24: The Parameters and Design Points View (p. 62)).
In the Ansys Workbench Project Schematic, double-click the Parameter Set bar below the
ext3d fluid flow analysis system.
Extra:
You can also open the Parameters and Design Points view by right-clicking the
Parameter Set bar and selecting the Edit... option in the context menu that opens.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 61
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
If you do not see the panes shown in the previous figure, make them visible by enabling Outline,
Properties, Table, and Chart from the View menu.
2. Run the calculation again with a new inlet flow rate for the current design point.
a. Enter 8E-6 under P1 - volumetric flow rate for the DP0 (Current) design point (cell B3)
in the Table of Design Points.
An Update Required icon will be added to the cell under P2- maxvelocity for the DP0
(Current) design point (cell C3 ).
b. Right-click the cell under P2 - maxvelocity for the DP0 (Current) design point and select
Update Selected Design Points in the context menu that opens, to generate the maximum
velocity at the flow exit with the revised inlet flow rate.
Extra:
You can also update the design point by clicking Update All Design Points in
the Ansys Workbench toolbar.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
62 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
A dialog box may open to inform you that some open editors may close during this process.
Click OK to proceed.
Ansys Polydata will update the data file based on the revised inlet flow rate and Ansys Polyflow will
run again. When the calculation is complete, the Table of Design Points will display a new value
of approximately 6.2 x 10-4 m/s under P2 - maxvelocity for the DP0 (Current) design point.
a. Click P1 under Input Parameters (cell A4) in the Outline of All Parameters.
The Ansys Workbench Toolbox will display options for Parameter Charts.
The Properties of Outline A11:0 window will display an initial setup for Parameter Chart 0
, in which P1 - volumetric flow rate is selected from the X-Axis (Bottom) drop-down list.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 63
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
c. Select P2 - maxvelocity from the Y-Axis (Left) drop-down list in the Properties of Outline
A11:0 window.
The current design point will be plotted in Parameter Chart 0 (Figure 1.25: The Chart of the Current
Design Point (p. 64)).
a. Enter 1E-5 for P1 - volumetric flow rate in the row beneath the DP0 (Current) design
point (cell B*) in the Table of Design Points, so that a new row is added (4) with DP 1 as
the Name.
Extra:
By default, Workbench only saves the calculated data for the design point in the
row labeled Current. You can specify that the data generated for any other design
points is saved within the project by enabling the Retain option in column D.
After the design points are updated, you can then right-click a design point in
the Table of Design Points and select Set as Current to access the data.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
64 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
5. Generate the values for the maximum velocity at the flow exit for all of the new design points.
Ansys Polydata will update and Ansys Polyflow will run repeatedly to solve for each of the design
points. As each calculation completes, the Table of Design Points (Figure 1.26: Displaying Values
for All of the Design Points (p. 65)) and Parameter Chart 0 (Figure 1.27: The Chart of All of the
Design Points (p. 65)) will be updated.
File → Save
7. Return to the Project Schematic view by clicking the Project tab above the Ansys Workbench tool-
bar.
8. View the files generated by Ansys Workbench, as displayed in the Project Schematic.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 65
Ansys Polyflow in Ansys Workbench Tutorial: 3D Extrusion
Note that the list of files shows that the design point file (designPoint.wbdp) was updated. For
more information about the files associated with Ansys Workbench, see the Ansys Workbench docu-
mentation.
1.5. Summary
In this tutorial, portions of Ansys Workbench were used to simulate a 3D extrusion and to compare the
flow exit velocities associated with a range of inlet flow rates.
Ansys DesignModeler was used to prepare the geometry, Ansys Meshing was used to create a compu-
tational mesh, Ansys Polydata was used to set up the simulation, Ansys Polyflow was used to calculate
the fluid flow throughout the geometry using the computational mesh, and CFD-Post was used to
analyze the results. In addition, the Parameters and Design Points view of Ansys Workbench was used
to add additional design points and compare their associated flow exit velocities on a chart.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
66 of ANSYS, Inc. and its subsidiaries and affiliates.
Part 2: Extrusion
The following extrusion tutorials are available:
1.1. Introduction
This tutorial illustrates the setup and solution of a 2.5D axisymmetric extrusion problem. The problem
corresponds to a simplified 2D simulation of a swirling flow that occurs around the head of an extrusion
screw. The fluid is forced through the die and exits the extruder after a short die land. The model involves
a free surface, the position of which is unknown.
• Create a sub-task.
• Set material properties and boundary conditions for a 2.5D axisymmetric extrusion problem.
1.2. Prerequisites
This tutorial assumes that you have little experience with Polyflow and its associated modules.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 69
2.5D Axisymmetric Extrusion
die. The position of the free surface is unknown. A portion of the mesh is affected by this unknown
boundary. A remeshing technique will be applied on this part of the mesh.
Since the problem involves a free surface, the domain is divided into two subdomains: one for the region
near the free surface and the other for the rest of the domain, as shown in Figure 1.2: Subdomains and
Boundary Sets for the Problem (p. 71)
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
70 of ANSYS, Inc. and its subsidiaries and affiliates.
Problem Description
The boundary sets for the problem are also shown in Figure 1.2: Subdomains and Boundary Sets for the
Problem (p. 71), and the conditions at the boundaries of the domains are:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 71
2.5D Axisymmetric Extrusion
1.4. Preparation
To prepare for running this tutorial:
3. Unzip the 25_axi_extrusion.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
Note:
If you create the mesh in GAMBIT or a third-party CAD package, you need to convert it
before you read it into Polydata. In this tutorial, the mesh file has already been converted.
So you can read the mesh file directly into Polydata.
1. Create a Fluid Flow (Polyflow) analysis system by drag and drop in Ansys Workbench.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
72 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Rename the project name to Tutorial 1 by double-clicking and editing the text Fluid Flow
(Polyflow).
Right-click the Mesh cell, hover over Import Mesh File… and click Browse....
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 73
2.5D Axisymmetric Extrusion
a. Select ext2d.msh.
b. Click Open.
3. Double-click the Setup cell to start Polydata and read in the mesh. When Polydata starts, the
Create a new task menu item appears in bold text, and the geometry for the problem is displayed
in the Graphics Display window.
Note:
At this point (when Create a new task appears in bold text) if you realize that you
have read the wrong mesh file, click STOP at the top of the Polydata menu and repeat
the process to access the correct mesh file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
74 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
• F.E.M. task
• Steady-state problem(s)
The Current setup (above the selected options) will be updated to reflect the selection. In
any problem solved using Polyflow, first an F.E.M. task is defined to calculate the flow field. If
information regarding the trajectories is necessary, specify a MIXING task after solving the
problem with the F.E.M. task specification and obtaining the results file. Then solve the problem
once again. 3D velocity components (u,v,w) are prescribed in a 2D cylindrical reference frame
(r,z), so 2D 1/2 axisymmetric geometry has been chosen. A steady-state condition is assumed
for this problem.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 75
2.5D Axisymmetric Extrusion
Note:
At this point (when Create a sub-task appears in bold text) if you realize that you
have made a mistake in the creation of the task and you need to return to that
menu, do the following:
ii. Select Redefine global parameters of a task and make the necessary
changes.
iii. Click Accept the current setup when you are satisfied with the corrected
settings.
Create a sub-task
A small dialog box appears asking for the title of the problem.
Note:
At this point (when Domain of the sub-task appears in bold text) if you realize
that you have made a mistake in the creation of the sub-task and you need to return
to that menu, do the following:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
76 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Select Redefine global parameters of a sub-task and make the necessary
changes.
Since this flow involves a free surface, the domain is divided into two subdomains: one for the
region near the free surface and the other for the rest of the domain. Here, the sub-task applies
to both subdomains (the default condition).
Accept the default selection of both subdomains by clicking Upper level menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 77
2.5D Axisymmetric Extrusion
when inertia, gravity, heat convection, or natural convection is taken into account. Since gravitational
effects are not included in the model, the default value of zero is retained for gravity.
Material data
(1.1)
where:
= shear rate
3. Specify the value , referred to as “fac” in the graphical user interface (compare the equation at
the top of the Cross law menu to Equation 1.1 (p. 78)).
Modify fac
Enter 85000 [units: poise] as the New value and click OK.
4. Specify the value for , referred to as “tnat” in the graphical user interface.
Modify tnat
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
78 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
5. Specify the value for , referred to as “expom” in the graphical user interface.
Modify expom
6. Check whether the values of the constants are correct, and repeat the previous steps if you need
to modify the constants again.
7. Click Upper level menu three times to leave the Material Data specification.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 79
2.5D Axisymmetric Extrusion
a. Select Zero wall velocity (vn=vs=0) along BS_1 and click Modify.
b. Click Inflow.
When the Automatic option is selected, Polydata chooses the most appropriate method to compute
the inflow. In this case, Polydata will use a 1D finite-element technique to compute a 1D fully-de-
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
80 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
veloped velocity profile, based on the specified material properties and flow rate. Moreover, the
inflow boundary condition requires that the computational domain be built in such a way that the
basic assumptions of fully-developed flow are satisfied. In axisymmetric geometries, the inflow
section must be perpendicular to the axial direction.
The fluid is assumed to stick to the wall, since at a solid-liquid interface the velocity of the liquid is that
of the solid surface. This is commonly known as the no-slip assumption because the liquid is assumed
to adhere to the wall, and therefore has no velocity relative to the wall.
Retain the default condition Zero wall velocity (vn=vs=0) along BS_2.
In a steady-state problem, the velocity field must be tangential to a free surface, since no fluid
particles go out of the domain through the free surface. This constraint is called the kinematic
condition, = 0. This equation requires an initial condition at the starting point of the free surface,
which in this case is located at the intersection of BS_2 and BS_3.
a. Select Zero wall velocity (vn=vs=0) along BS_3 and click Modify.
Note:
Do not select the Outlet option. It is only applicable for die design problems.
As mentioned above, the starting point of the free surface is at the intersection of BS_2 and BS_3.
h. Retain the default settings for the Normal force and Direction of motion.
j. Click Upper level menu to return to the Flow boundary conditions menu.
It is reasonable to consider that a uniform velocity profile is obtained at the exit. In most cases, a bulk
flow is obtained and therefore no force is acting, so the selection of zero normal and tangential forces
is appropriate. In situations involving pulling velocity or force or gravity, the corresponding boundary
condition should be selected.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 81
2.5D Axisymmetric Extrusion
a. Select Zero wall velocity (vn=vs=0) along BS_4 and click Modify.
c. Accept the default value of 0 for the normal force by clicking Upper level menu.
d. Accept the default value of 0 for the tangential force by clicking Upper level menu.
5. Retain the default condition at the symmetry axis (Axis of symmetry along BS_5).
For axisymmetric models, the axis of symmetry is always the y axis. Polydata determines the axis of
symmetry from the mesh file, and automatically imposes the symmetry condition along the line r=0
(x=0).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
82 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Since the screw is rotating with angular velocity = 2 = 6.2832 rad/s, the rotational velocity along
this boundary is prescribed to increase linearly with ( = 6.2832 ). In the equation for , X denotes
the direction and Y denotes the direction. Since the fluid sticks to the wall, = 0 = .
a. Select Zero wall velocity (vn=vs=0) along BS_6 and click Modify.
c. Accept the default value of 0 for the normal velocity ( ) and tangential velocity ( ) by selecting
Upper level menu twice.
e. Click Velocity w imposed and select 'w' velocity = linear function of coordinates.
f. Accept the default value of 0 [units: cm/s] for the constant A by clicking OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 83
2.5D Axisymmetric Extrusion
g. Enter 6.2832 [units: rad/s] as the New value for the constant B and click OK.
h. Accept the default value of 0 [units: rad/s] for the constant C by clicking OK.
j. Click Upper level menu at the top of the Flow boundary conditions menu.
1.5.5. Remeshing
This model involves a free surface for which the position is unknown. A portion of the mesh is affected by
this unknown boundary. Hence a remeshing technique is applied on this part of the mesh. The free surface
is entirely contained within subdomain 2, and hence only subdomain 2 will be affected by the relocation
of the free surface.
Global remeshing
In some cases, when the mesh is geometrically complex, it may be necessary to split it into additional
subdomains in order to define a specific remeshing method on each of them. For this purpose, Polydata
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
84 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
allows you to create several local remeshings. For the current problem, a single local remeshing is
sufficient.
SD_1 is moved from the top list to the bottom list, indicating that only SD_2 will be remeshed.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 85
2.5D Axisymmetric Extrusion
If you accidentally remove the wrong subdomain, select it and click Add to restore it. Then, follow
the instructions to remove the correct subdomain.
The purpose of the remeshing technique is to relocate internal nodes according to the displacement
of boundary nodes due to the motion of the free surface. Mesh nodes are organized along lines of
remeshing (spines), which are collections of nodes logically arranged in a one-dimensional manner.
This technique is most suited for 2D extrusion problems. Polydata requires the specification of the first
and last spines that the fluid encounters (inlet of spines and outlet of spines, respectively).
In this case, the inlet of spines is the intersection of subdomain 2 with subdomain 1, and the outlet of
spines is the intersection of subdomain 2 with the flow exit (boundary 4).
Method of Spines
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
86 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
At this point, if you realize that you have made a mistake in global remeshing, click die swell at
the bottom of the menu and perform this Step again.
1. Select Condition on the stream function for field 1. Click No in the window that pops up.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 87
2.5D Axisymmetric Extrusion
If you have made a mistake in assigning the stream function, click F.E.M. Task 1 to get into that
menu and then repeat this Step.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
88 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1.5.7. Outputs
After Polyflow calculates a solution, it can save the results in several different formats. Choose the format
that is appropriate for your postprocessor. In this case, save the outputs in the default format for Ansys
CFD-Post.
Outputs
2. Accept the default output option for CFD-Post by clicking Upper level menu.
When exiting the menu, Polydata asks you to confirm the current system units and fields that are to
be saved to the results file for postprocessing.
If you do not enter the menu Outputs, Polydata will ask you to confirm the current system units at
the end of the session, if it is a new session.
If this was not yet done before (see above), Polydata asks you to confirm convergence strategies.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 89
2.5D Axisymmetric Extrusion
Polydata also asks you to confirm the current system units. It will also ask to confirm fields that are
to be saved to the results file for postprocessing.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
90 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
2. Click Accept.
3. Click Continue.
This accepts the default names for graphical output files ( cfx.res ) that are to be saved for post-
processing, and for the Polyflow format results file ( res ).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 91
2.5D Axisymmetric Extrusion
1.5.9. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update from the
shortcut menu.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Ansys Workbench opens the View listing file dialog box, which displays the listing file.
b. In the View listing file dialog box, disable the Show only last 200 lines option and click the
Manual refresh button. Then find the SOLVER section that relates to F.E.M. Task 1; at
the end of this section, a message Convergence assumed is displayed. This indicates that
the solution has converged. See the Polyflow User's Guide for more information on convergence.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
92 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1.5.10. Postprocessing
Ansys CFD-Post has similar interfaces for UNIX and Windows, the postprocessing steps are illustrated for
Windows.
1. Double-click the Results cell in the Ansys Workbench analysis and read the results files saved by
Polyflow.
Ansys CFD-Post reads the solution fields that were saved to the results file.
a. Right-click a blank area anywhere in the graphical window, hover over Predefined Camera
in the context menu, and select View From +Z.
The central mouse button allows you to zoom in and zoom out. The left mouse button allows ro-
tating the image. The right mouse button allows you zoom to an area.
b. Also, right-click a blank area anywhere in the graphical window and deselect Ruler, if needed.
a. Click the Insert menu and select Contour or click the Contour button ( ).
b. Click OK to accept the default name (Contour 1) and open the details view below the Outline
tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 93
2.5D Axisymmetric Extrusion
c. In the details view for Contour 1, specify the following settings under Geometry:
i. Next to Locations, click the ellipsis button ( ) on the right and select SD_1_surf and
SD_2_surf (use Ctrl to select multiple items).
ii. Select PRESSURE from the Variable drop-down list, or click the ellipsis button ( ) on
the right and select PRESSURE.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
94 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 95
2.5D Axisymmetric Extrusion
The details view of Default Transform will open below the Outline tab.
ii. In the details view of Default Transform, disable Instancing Info From Domain.
iv. Select Y from the Axis drop-down list in the Axis Definition group box.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
96 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Click the Insert menu and select Text or click the button.
ii. Click OK to accept the default name (Text 1) and open the details view below the Outline
tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 97
2.5D Axisymmetric Extrusion
iii. Enter Polyflow Results for Text String in the Definition tab of the details view.
a. Deselect Contour 1.
b. Click the Insert menu and select Vector or click the Vector button ( ).
c. Click OK to accept the default name (Vector 1) and open the details view below the Outline
tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
98 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. In the Geometry tab, click the ellipsis button ( ) next to Locations, select SD_1_surf
and SD_2_surf (use Ctrl to select multiple items), and click OK.
ii. Select VELOCITIES from the Variable drop-down list, or click the ellipsis button ( ) on
the right and select VELOCITIES.
iv. Define the attributes of vectors: in the Symbol tab, set Symbol to Arrow3D and Symbol
Size to 3.
v. Click Apply.
i. Deselect Text 1 in the Outline tab, under User locations and plots.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 99
2.5D Axisymmetric Extrusion
i. Rotate to the isometric view by clicking the cyan-blue dot in the axis triad (bottom right
of the graphics window).
This allows you to better visualize the magnitude of the velocity vectors.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
100 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Enable Normalize Symbols in the Symbol tab of the details view for Vector 1.
This allows you to better visualize the direction of the velocity vectors.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 101
2.5D Axisymmetric Extrusion
The velocity vectors take all components of the velocity into account. Along the screw tip, the
rotational component is important, leading to long vectors that are not in the xy plane. After
the die exit, a rearrangement of the velocity field takes place. The flow slows down along the
axis of symmetry and accelerates on the outside. This makes the particles go toward the free
surface, creating the swelling.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
102 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Under Render tab, deselect Show Faces and select Show Mesh Lines.
iii. Under Color tab, click the ellipsis ( ), select the color white and click OK.
i. Double-click Vector 1 under User Locations and Plots in the Outline tab.
i. Right-click a blank area in the graphics window, hover over Predefined Camera in the context
menu, and select View from +Z.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 103
2.5D Axisymmetric Extrusion
i. Move the mouse to the left-hand border of the graphic window until the cursor suggests
a rotation along a vertical line.
ii. Click and move the mouse slowly to the right-hand side.
1.6. Summary
This tutorial demonstrated how to set up and solve a 2.5D axisymmetric extrusion problem. It showed
how to set up a free surface problem and the associated remeshing, and demonstrated the use of CFD-
Post to examine the flow behavior associated with the problem.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
104 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 2: Fluid Flow and Conjugate Heat Transfer
This tutorial is divided into the following sections:
2.1. Introduction
2.2. Prerequisites
2.3. Problem Description
2.4. Setup and Solution
2.5. Summary
2.1. Introduction
This tutorial illustrates the setup and solution of a problem involving heat transfer between a Newtonian
fluid and a cooled circular die. Along with a good die design, rheological and thermo physical properties
of the melt and the thermal settings in the die are very important in obtaining a geometrically well-
defined polymer product. The heat transfer calculation is important when temperature-sensitive polymers
are shaped and when product surface qualities are of critical importance. The temperature field at the
die exit influences the swelling and drawing behavior of the product.
In this tutorial, you will solve the non-isothermal flow problem for the fluid and the heat conduction
in the die, making some assumptions regarding the rheological and thermo physical properties of the
melt.
• Set material properties and boundary conditions for a fluid-solid heat conduction and flow problem.
2.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the set up procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 105
Fluid Flow and Conjugate Heat Transfer
To solve the coupled problem, two sub-tasks are defined: one for the fluid (sub-task 1) and the other
for the solid (sub-task 2). Each sub-task will contain a particular model, domain of definition, material
properties, and boundary conditions, including interface conditions with the other sub-task. The sub-
tasks are coupled, because the global solution of the problem depends on the values of the solution
variables at the intersection of the fluid and solid domains.
Figure 2.1: A Schematic Diagram of the Fluid and the Circular Die
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
106 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Viscous heating is taken into account. For the solid region, the thermal conductivity is 35 W/m-˚C.
The boundary sets for the problem are shown in Figure 2.2: Boundaries and Sub-domains (p. 107), and
the conditions at the boundaries of the domains are as follows:
• boundary 2: insulated
• boundary 3: = 100°C
• boundary 4: insulated
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 107
Fluid Flow and Conjugate Heat Transfer
2.4.8. Postprocessing
2.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the fluid_solid.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering fluid-solid as the name of
the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Steady-state problem(s)
• 2D axisymmetric geometry
Since the problem involves an axisymmetric steel die, the computational domain for the problem
is chosen to be a 2D cylindrical reference frame (r,z) with r=0 as the axis of symmetry, and involves
two velocity components (u,v); hence 2D axisymmetric geometry has been chosen. A Steady-state
condition is assumed for the problem.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
108 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Create a sub-task
Note:
To solve the coupled problem, the computational domain is divided into two sub-domains with
a common intersection. A sub-task with its own model, material properties, and boundary condi-
tions is defined on each of the non-overlapping subdomains. Sub-task 1 is defined for SUBDO-
MAIN_1, since SUBDOMAIN_1 represents the fluid (as shown in Figure 2.2: Boundaries and Sub-
domains (p. 107)).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 109
Fluid Flow and Conjugate Heat Transfer
SUBDOMAIN_2 is moved from the top list to the bottom list, indicating that subtask 1 is
defined on SUBDOMAIN_1.
Polydata indicates which material properties are relevant for your sub-task by graying out the irrelevant
properties. In this sub-task, Polyflow solves energy, incompressibility and momentum equations, so
you have to define viscosity, density, thermal conductivity, heat capacity per unit mass, and viscous
heating. For a non-isothermal flow problem, the viscosity can depend on both shear rate and temper-
ature. In this case, the viscosity is constant, so it depends on neither of them.
Material Data
Since the fluid flow is Newtonian, specify a constant value for the viscosity.
ii. Specify the value for , referred to as “fac” in the graphical user interface.
Modify fac
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
110 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. Enter 2500 [units: Pa•s] as the New value and click OK.
iv. Click Upper level menu two times to continue the Material Data specification.
Polydata displays the following message, confirming that there is no temperature dependence
for the viscosity.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 111
Fluid Flow and Conjugate Heat Transfer
iii. Click Upper level menu to continue the Material Data specification.
c. Click Density.
Modification of density
i. Enter 950 [units: kg/m3] as the New value and click OK.
ii. Click Upper level menu to continue the Material Data specification.
As shown at the top of the menu, the thermal conductivity is defined as a nonlinear function of
the temperature:
(2.1)
In this problem, the thermal conductivity is assumed to be a constant for the fluid so only
the constant coefficient is modified.
Modify a
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
112 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Enter 0.5 [units: W/m-˚C] as the New value and click OK.
ii. Click Upper level menu to continue with the Material Data specification.
As shown at the top of the menu, the heat capacity per unit mass is defined as a nonlinear
function of temperature:
(2.2)
The temperature variation of differs with the nature of the polymer melts. In this problem,
is assumed to be constant, so only the constant coefficient is modified.
Modify a
i. Enter 2300 [units: J/kg-˚C] as the New value and click OK.
ii. Click Upper level menu to continue with the Material Data specification.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 113
Fluid Flow and Conjugate Heat Transfer
When shearing occurs in a flow, the friction of the different fluid layers generates heat. When
the fluid is highly viscous and/or the shear rate is high, the heating of the fluid caused due to
this phenomenon must be taken into account.
ii. Click Upper level menu to return to the Material Data specification.
The intersection acts as a wall for the fluid, and since the fluid is assumed to stick to the wall,
zero normal and tangential velocities is imposed along this boundary.
i. Retain the default condition Zero wall velocity (vn=vs=0) along SUBDOMAIN_2.
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY_1 and click Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
114 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. Retain the default settings, Automatic and Volumetric flow rate.
v. Enter 0.6e-06 [units: m3/s] as the New value in the dialog box that appears, and click
OK.
The flow rate of the melt is very low due to the highly viscous nature of the melt.
When the Automatic option is selected, Polydata automatically chooses the most appro-
priate method to compute the inflow condition.
It is assumed that a fully developed velocity profile is reached at the exit, so the outflow condition
is appropriate. This condition imposes a zero normal force, (which includes a pressure term),
and zero tangential velocity, .
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY_5 and click Modify.
For axisymmetric models, Polydata recognizes the axis of symmetry from the mesh file, and
automatically imposes the symmetry condition along the line r=0. This condition imposes a zero
normal velocity and zero tangential force along this boundary.
e. Click Upper level menu at the top of the Flow boundary conditions menu to return to
the fluid menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 115
Fluid Flow and Conjugate Heat Transfer
For non-isothermal problems, specify either the temperature or the heat flux on each boundary seg-
ment. The temperature along a given boundary can be a constant or a prescribed function of coordin-
ates.
An interface condition is set at the intersection of subdomain 1 and subdomain 2. This condition
ensures continuity of the temperature field and of the heat flux along the interface. Since you
are solving a coupled problem, this condition of continuity is essential for the global solution of
the temperature and heat flux variables.
iii. Click Upper level menu to accept the default setting (continuous heat flux along the
interface).
In the case of an interface condition, both the heat flux and temperature are usually con-
tinuous along the interface. It is possible to specify a nonzero value for the heat flux jump
( ), but this is mainly used in problems where internal radiation is simulated. Here, accept
the default value for the definition of heat flux discontinuity, = 0.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
116 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Polydata prompts you for the new value of the constant temperature.
iv. Enter 180 [units: ˚C] as the New value and click OK.
v. Click Upper level menu to return to the Thermal boundary conditions menu.
e. Click Upper level menu twice to return to the F.E.M. Task 1 menu.
Create a sub-task
a. Polydata warns you that a thermal interface has been defined, and that no other sub-task shares
the interface. Click OK to close the warning.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 117
Fluid Flow and Conjugate Heat Transfer
c. Click No, since this sub-task has different parameters associated with it.
A small dialog box appears asking for the title of the problem.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
118 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Material Data
In this problem, thermal conductivity is assumed to be a constant, so only the constant coefficient
is modified.
i. Select Modify a.
ii. Enter 35 [units: W/m-˚C] as the New value and click OK.
b. Click Upper level menu two times to return to the solid menu.
In this step, set the conditions at each of the boundaries of the domain. When a boundary set
is selected, it is highlighted in red in the graphics window.
iii. Click Upper level menu to accept the default option for continuity of temperature and
heat flux.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 119
Fluid Flow and Conjugate Heat Transfer
Polydata prompts you for the new value of the constant temperature.
iv. Enter 100 [units: ˚C] as the New value and click OK.
v. Click Upper level menu to return to the Thermal boundary conditions menu.
e. Click Upper level menu three times to return to the top-level Polydata menu.
Polydata asks you to confirm the current system units and fields that are to be saved to the results
file for postprocessing.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
120 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
6. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and for the Polyflow format results file (res).
2.4.7. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
2.4.8. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
1. Double-click the Results cell in the Workbench analysis and read the results files saved by Polyflow.
CFD-Post reads the solution fields that were saved to the results file.
Right-click the Graphics window and select View from +Z under Predefined Camera.
(Or you can click +Z on the axis triad in the graphic window.)
a. Click the Insert menu and select Contour or click the Contour button ( ).
b. Click OK to accept the default name (Contour 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 121
Fluid Flow and Conjugate Heat Transfer
c. In the details view for Contour 1, specify the following settings under Geometry:
i. Next to Locations, click the ellipsis button ( ) on the right and select SUBDO-
MAIN_1_surf and SUBDOMAIN_2_surf (use Ctrl to select multiple items).
ii. Select PRESSURE from the Variable drop-down list, or click the ellipsis button ( )
on the right and select PRESSURE.
Most of the pressure drop occurs in the upper part of the die where the cross-section is smallest (Fig-
ure 2.3: Pressure Contours (p. 123)). The pressure is linear except in the contraction zone. The isobars
are perpendicular to the flow direction, as expected for the fully developed flow that occurs in the
second part of the die.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
122 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. In the details of Contour 1, select VELOCITIES from the Variable drop-down list.
b. Click Apply.
The velocity is higher in the second part of the die where the cross-section is smaller (Figure 2.4: Velocity
Distribution (p. 124)). It reaches a maximum value in the center of the thin tube.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 123
Fluid Flow and Conjugate Heat Transfer
a. In the Outline tab, under User Locations and Plots, deselect Contour 1.
i. Click the Insert menu and select Vector or click the button.
iii. In the Geometry tab of the details view of Vector 1, click the button next to Loc-
ations.
iv. Select the location SUBDOMAIN_1 and click OK to close the Location Selector dialog
box.
v. In the Symbol tab, select Arrow 3D and increase the Symbol Size to 3.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
124 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The flow is fully developed in the downstream part of the die (Figure 2.5: Velocity Vectors (p. 125)).
Observe the classical parabolic velocity profile. The Poiseuille flow is rapidly reached after the contrac-
tion because inertia is not taken into account here.
6. Display of the temperature distribution in the solid and the fluid regions.
a. Deselect Vector 1 in the Outline tab under User Locations and plots.
b. Select, then double-click Contour 1 in the Outline tab under User Locations and plots.
c. In the details view for Contour 1, specify the following settings under Geometry:
ii. Select Temperature from the Variable drop-down list, or click the ellipsis button ( )
on the right and select Temperature.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 125
Fluid Flow and Conjugate Heat Transfer
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
126 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
As expected, the temperature gradients are larger in the fluid than in the die, (the isolines are closer
together in the fluid than in the die—see Figure 2.6: Temperature Distribution (Celsius) (p. 126)). This
is because the conductivity of the fluid is much lower than the conductivity of the solid. The temper-
ature isolines are perpendicular to the boundaries where the (normal) flux becomes zero. The heating
of the fluid due to viscous dissipation can be clearly seen. In order to visualize contours in Kelvin, select
the Edit/Options... menu item, click Units under Common, select K from the Temperature drop-
down menu, and click OK.
ii. Click OK to accept the default name (Line 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 127
Fluid Flow and Conjugate Heat Transfer
iv. Enter 0, 0.006, and 0 for Point 1 and enter 0.008, 0.006, and 0 for Point 2.
Note:
You will need to ensure that your unit of length is set to meter in CFD-Post.
v. Click Apply.
b. Create a plot.
i. Click the Insert menu and select Chart or click the Chart button ( ).
ii. Click OK to accept the default name (Chart 1) and display the details view below the
Outline tab tree.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
128 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. In the General tab, ensure XY is selected for the Chart Type and enter Temperature
Profile for the Title.
iv. In the Data Series tab, select Line 1 from the Location drop-down list.
vi. In the Y Axis tab, select Temperature from the Variable drop-down list.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 129
Fluid Flow and Conjugate Heat Transfer
The thermal boundary layer located along the die wall is clearly visible. This boundary layer is the
result of the low thermal conductivity and high heat capacity of the fluid. The temperature of the
fluid at the center is not affected by the low temperature of the solid. The heat does not diffuse quickly
enough through the fluid layer to reach the axis of symmetry, before the fluid exits the die.
2.5. Summary
This tutorial introduced the coupling of sub-tasks of different types: a non-isothermal flow problem and
a heat conduction problem in a solid. Coupled calculations like this are very useful in polymer processing
applications where thermal effects are critical (for example: extrusion, coating, fiber spinning). Coupling
can also be applied through fields other than temperature (for example: electrical potential and pressure
in porous media).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
130 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 3: Non-Isothermal Flow Through a Cooled
Die
This tutorial is divided into the following sections:
3.1. Introduction
3.2. Prerequisites
3.3. Problem Description
3.4. Setup and Solution
3.5. Summary
3.1. Introduction
This tutorial examines the flow of a polymer melt through a die. The temperature of the melt increases
due to viscous dissipation caused by the shearing taking place in the die. The temperature of the fluid
is critical for the process. The viscosity of the fluid changes with temperature, which leads to the
modification of the shape of the extrudate. The polymer might degrade if the temperature is too high,
so a numerical simulation is of great interest to optimize the operating conditions.
• Set material properties and boundary conditions for the contraction flow problem.
3.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the set up procedure
will not be shown explicitly.
In solving for the free surface location, the position variables are also coupled to the temperature, velocity,
and pressure fields. To solve the coupled problem, you will define two sub-tasks: one each for the fluid
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 131
Non-Isothermal Flow Through a Cooled Die
(sub-task 1) and the solid (sub-task 2). Each sub-task contains a particular model, domain of definition,
material properties, and boundary conditions, including interface conditions with the other sub-task.
The sub-tasks are coupled because the global solution of the problem depends on the values of the
solution variables at the intersection of the fluid and solid domains.
The high flow rate introduces strong nonlinearity in the problem, which can lead to a loss of convergence
in the iterative scheme. In Polyflow, convergence strategies are available to solve such highly nonlinear
problems.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
132 of ANSYS, Inc. and its subsidiaries and affiliates.
Problem Description
Viscous heating is taken into account and the shear-rate dependence of viscosity obeys the Bird-Carreau
law. For the solid region, the thermal conductivity ( ) is 30 W/m-°C.
The boundary sets for the problem are shown in Figure 3.2: Boundaries and Subdomains (p. 134), and
the flow and thermal conditions for the fluid and the die at the boundaries of the domains are:
• boundary 4: free surface with convective heat transfer to surroundings ( = 20 W/m2-°C, α= 20°C)
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 133
Non-Isothermal Flow Through a Cooled Die
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
134 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
3.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the non_iso_flow.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering non-iso-flow as the name
of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• 2D axisymmetric geometry
The Current setup (above the selected options) is updated to reflect your selections. Since the
problem involves an axisymmetric die, Polyflow uses a 2D cylindrical reference frame (r,z) with r=0
as the axis of symmetry.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 135
Non-Isothermal Flow Through a Cooled Die
Create a sub-task
Note:
To solve the coupled problem, the computational domain is divided into three subdomains. There
are two sub-tasks in this problem. Define a sub-task with its own model, material properties, and
boundary conditions for the fluid region. Since this problem involves a free surface, the domain for
sub-task 1 is divided into two subdomains: one for the region near the free surface (SUBDOMAIN_2)
and the other for the rest of the fluid domain (SUBDOMAIN_1). In this problem, sub-task 1 applies to
SUBDOMAIN_1 and SUBDOMAIN_2.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
136 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
SUBDOMAIN_3 is moved from the top list to the bottom list, indicating that subtask 1 is defined
on SUBDOMAIN_1 and SUBDOMAIN_2.
b. Click Upper level menu at the top of the Domain of the sub-task menu.
Polydata indicates the material properties that are relevant for your sub-task by graying out the irrel-
evant properties. In this sub-task, Polyflow solves energy, incompressibility, and momentum equations.
Hence, define viscosity, density, thermal conductivity, and heat capacity per unit mass. For a non-
isothermal generalized Newtonian fluid, the viscosity depends on the shear rate and the temperature.
Hence, define the shear-rate dependence of viscosity and the temperature dependence of viscosity.
Material data
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 137
Non-Isothermal Flow Through a Cooled Die
where is the viscosity at zero shear rate, is the shear rate, is the Bird-Carreau
law index, and is the natural time.
ii. Specify the value , referred to as “fac” in the graphical user interface (compare the
equation at the top of the Bird-Carreau law menu to Equation 3.1 (p. 137)).
Modify fac
Enter 5000 [units: Pa•s] as the New value and click OK.
iii. Specify the value , referred to as “tnat” in the graphical user interface.
Modify tnat
iv. Specify the value for , referred to as “expo” in the graphical user interface.
Modify expo
When you click Upper level menu , Polydata displays the following warning message:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
138 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
For this tutorial, you will use the convergence strategy for viscosity/slip.
vii. Click Upper level menu again to continue with the Material Data specification.
For this problem, assume that the dependence of viscosity on temperature follows the Arrhenius
law.
where is the ratio of the activation energy to the thermodynamic constant and is a
reference temperature for which = 1. The parameter denotes the absolute 0 tem-
perature in your selected temperature scale. It is set to 0, when and are absolute
temperatures. In this example, specify the temperatures in Celsius, so enter a value of -273
for .
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 139
Non-Isothermal Flow Through a Cooled Die
ii. Specify the value for , referred to as “alfa” in the graphical user interface (compare
the equation at the top of the Temperature dependence of viscosity menu to Equa-
tion 3.2 (p. 139)).
Modify alfa
Enter 2300 [units: 1/°C] as the New Value and click OK.
iii. Specify the value for , referred to as “talfa” by the graphical user interface.
Modify talfa
Enter 200 [units: °C] as the New Value and click OK.
iv. Specify the value for , referred to as “t0” by the graphical user interface.
Modify t0
Enter -273 [units: °C] as the New Value and click OK.
v. Click Upper level menu two times to continue with the Material Data specification.
c. Click Density.
Modification of density
i. Enter 950 [units: kg/m3] as the New value and click OK.
ii. Click Upper level menu to continue with the Material Data specification.
For this problem, the thermal conductivity of the fluid is assumed to be a constant. So only the
constant coefficient is modified.
Modify a
i. Enter 0.5 [units: W/m-°C] as the New Value and click OK.
ii. Click Upper level menu to continue with the Material Data specification.
The heat capacity per unit mass is defined as a nonlinear function of temperature:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
140 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
(3.4)
The temperature variation of depends on the nature of the polymer melt. For this problem,
is assumed to be constant, so only the constant coefficient is modified.
Modify a
i. Enter 2300 [units: J/kg-°C] as the New value and click OK.
ii. Click Upper level menu to continue with the Material Data specification.
When shearing occurs in a flow, the friction of the different fluid layers generates heat. When
the fluid is highly viscous and/or the shear rate is high, the heating of the fluid caused by this
phenomenon must be taken into account.
ii. Click Upper level menu to return to the Material Data specification.
a. Retain the default condition Zero wall velocity (vn=vs=0) along SUBDOMAIN_3 at the
intersection of SUBDOMAIN_1 and SUBDOMAIN_3.
The liquid is assumed to stick to the wall, since at a solid-liquid interface the velocity of the liquid
is that of the solid surface. This is known as the no-slip assumption because the liquid is assumed
to adhere to the wall, and hence, has no velocity relative to the wall.
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY_1 and click Modify.
Polydata prompts for the new value of the volumetric /mass flow rate.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 141
Non-Isothermal Flow Through a Cooled Die
Enter 5e-06 [units: m3/s] as the New Value and click OK.
When the Automatic option is selected, Polydata automatically chooses the most appro-
priate method to compute the inflow condition.
vi. Click Upper level menu to return to the Flow boundary conditions menu.
For axisymmetric models, Polydata recognizes the axis of symmetry from the mesh file, and
automatically imposes the symmetry condition along the line = 0. This condition imposes a
zero surface normal velocity ( ) and zero tangential force ( ) along this boundary.
It is assumed that a uniform velocity profile is reached at the exit. The melt is not subjected to
any externally applied stress at the exit, so the condition of zero normal and tangential forces
is selected.
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY_3 and click Modify.
iii. Click Upper level menu to accept the default value of 0 [units: Pa] for .
iv. Click Upper level menu to accept the default value of 0 [units: Pa] for .
In a steady-state problem, the velocity field must be tangential to a free surface, since no fluid
particles leave the domain through the free surface. This constraint is called the kinematic con-
dition, = 0. This equation requires an initial condition, which is the starting line of the free
surface. In this problem, the starting line of the free surface is the intersection of BOUNDARY_4
and SUBDOMAIN_3 (see Figure 3.2: Boundaries and Subdomains (p. 134) ).
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY_4 and click Modify.
vii. Click Upper level menu to return to the Kinematic condition menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
142 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Polydata prompts for the new value of the Y-component of the direction-of-displacement
vector.
xiv. Click Upper level menu to return to the Kinematic condition menu.
xv. Click Upper level menu to return to the Flow boundary condition menu.
For non-isothermal problems, specify either the temperature or the heat flux on each boundary set.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 143
Non-Isothermal Flow Through a Cooled Die
iii. Click Upper level menu to accept the default setting (continuous heat flux along the
interface).
For an interface condition, both the heat flux and temperature are usually continuous along
the interface. It is possible to specify a nonzero value for the heat flux jump ( ), but this
is mainly used in problems where internal radiation is simulated. Accept the default value
for the definition of heat flux discontinuity ( =0).
iv. Enter 200 [units: °C] as the New Value and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
144 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
v. Click Upper level menu to return to the Thermal boundary conditions menu.
If the heat transfer from radiation is neglected, the heat flux can be written as
(3.5)
where is the heat convection coefficient and is the reference temperature (in this case,
the temperature of the air surrounding the extrudate).
Modify alpha
Modify Talpha
v. Click Upper level menu to return to the Thermal boundary conditions menu.
This model involves a free surface for which the position is unknown. A portion of the mesh is affected
by the relocation of this boundary. Hence, a remeshing technique is applied on this part of the mesh.
The free surface is entirely contained within SUBDOMAIN_2 and hence, only SUBDOMAIN_2 is affected
by the relocation of the free surface.
Global remeshing
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 145
Non-Isothermal Flow Through a Cooled Die
SUBDOMAIN_1 is moved from the top list to the bottom list, indicating that only SUBDO-
MAIN_2 will be remeshed.
The purpose of the remeshing technique is to relocate internal nodes according to the displace-
ment of boundary nodes due to the motion of the free surface. Mesh nodes are organized along
lines of remeshing (spines), which are collections of nodes logically arranged in a one-dimensional
manner.
Polydata requires the specification of the first and last spines (inlet and outlet) that the fluid
encounters. In this case, the inlet of spines is the intersection of SUBDOMAIN_2 with SUBDO-
MAIN_1, and the outlet of spines is the intersection of SUBDOMAIN_2 with the flow exit
(BOUNDARY_3).
Method of Spines
i. Specify the inlet for the system of spines by selecting SUBDOMAIN_1 and clicking
Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
146 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Specify the outlet for the system of spines by selecting BOUNDARY_3 and clicking
Modify.
iii. Click Accept the current setup in the Element distortion check menu.
The finite-element mesh can undergo great deformations. The Element distortion
check menu deals with the detection of all possible distortions of the elements.
For this problem, accept the default options and proceed to the next step.
Interpolation
You can expect important temperature gradients in the calculation. Therefore, you can retain the
quadratic interpolation (9 unknowns per element) for velocity and the linear interpolation (4 unknowns
per element) for pressure, but it is recommended that you select the 4x4 interpolation for temperature.
In the 4x4 discretization scheme, each finite element is divided into 16 sub-elements, with the temper-
ature being linearly interpolated over each sub-element. This leads to 25 temperature unknowns per
element.
a. Scroll down to select 4x4 element for temperature in the Interpolation menu.
The Current setup (at the top of the menu) is updated to reflect your selection.
b. Click Upper level menu two times to return to the F.E.M. Task 1 menu.
Create a sub-task
a. Polydata warns you that a thermal interface has been defined, and that no other sub-task shares
the interface. Click OK to close the warning.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 147
Non-Isothermal Flow Through a Cooled Die
c. Click No, since this sub-task has different parameters associated with it.
c. Click Upper level menu at the top of the Domain of the sub-task menu.
For this problem, specify a constant value for the thermal conductivity .
Material data
For this problem, thermal conductivity is assumed to be a constant, so only the constant coefficient
is modified.
b. Select Modify a.
d. Click Upper level menu two times to return to the solid menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
148 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Set the conditions at each of the boundaries of the domain. The selected boundary set will be high-
lighted (in red) in the graphics window as you select them..
iii. Click Upper level menu to accept the default option for continuity of temperature and
heat flux.
Take only the heat convection into account: see Equation 3.5 (p. 145) .
Modifiy alpha
Modifiy Talpha
v. Click Upper level menu to return to the Thermal boundary conditions menu.
Modify alpha
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 149
Non-Isothermal Flow Through a Cooled Die
Modifiy Talpha
v. Click Upper level menu to return to the Thermal boundary conditions menu.
5. Click Upper level menu twice to return to the F.E.M. Task 1 menu.
3.4.6. Outputs
Outputs
Set to metric_MKSA+Celsius
2. Click Upper level menu three times to return to the top-level Polydata menu.
1. Select Enable convergence strategy for rheology / slip in the Convergence strategies menu,
and click Accept.
3. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and for the Polyflow format results file (res).
3.4.8. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
150 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file panel, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
3.4.9. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
1. Double-click the Results cell in the Workbench analysis and read the results files saved by Polyflow.
CFD-Post reads the solution fields that were saved to the results file.
In the graphical window, right-click, and select the option Predefined Camera.
a. Right-click in the graphical window and select View from +Z under Predefined Camera.
b. To remove the ruler right-click in the graphical window, select Viewer Options, and disable
Ruler Visibility.
a. Click the Insert menu and select Contour or click the button.
b. In the panel that opens, click OK to accept the default name (Contour 1) display the details
view below the Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 151
Non-Isothermal Flow Through a Cooled Die
c. Perform the following steps In the Geometry tab of the details view for Contour 1:
i. Next to Locations, click the ellipsis button ( ) on the right and select SUBDO-
MAIN_1_surf and SUBDOMAIN_2_surf (use Ctrl to select multiple items).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
152 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Select PRESSURE from the Variable drop-down list, or click the ellipsis button ( )
on the right and select PRESSURE.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 153
Non-Isothermal Flow Through a Cooled Die
a. In the details view of Contour 1, select VELOCITIES from the Variable drop-down list.
b. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
154 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The fluid experiences high velocity gradients in the narrow section of the die. This leads to important
viscous dissipation effects that cause the temperature of the melt to increase.
a. In the Outline tab under User Locations and Plots, disable Contour 1.
b. Click the Insert menu and select Vector or click the button.
c. Click OK to accept the default name (Vector 1) and open the details view below the Outline
tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 155
Non-Isothermal Flow Through a Cooled Die
i. In the Geometry tab, click the button next to Locations to open the Location
Selector dialog box.
ii. Select SUBDOMAIN_1_surf and SUBDOMAIN_2_surf (use Ctrl to select multiple items).
iii. In the Symbol tab, select Arrow3D and retain the default Symbol Size of 1.0.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
156 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The velocity vectors in the wide section of the die are very small compared to those in the narrow
section of the die (Figure 3.5: Velocity Vectors (p. 157)). Also, the important velocity re-arrangement
takes place at the die exit. This leads to the swelling of the extrudate.
6. Display the temperature distribution in the solid and the fluid regions.
a. In the Outline tab, under User Locations and Plots, disable Vector 1, and enable and
double-click Contour 1.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 157
Non-Isothermal Flow Through a Cooled Die
i. Next to Locations, click the ellipsis button ( ) on the right and select SUBDO-
MAIN_1_surf, SUBDOMAIN_2_surf and SUBDOMAIN_3_surf (use Ctrl to select multiple
items).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
158 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 159
Non-Isothermal Flow Through a Cooled Die
Figure 3.7: Temperature Profile Near the Die Exit (p. 160) shows a magnified view of the temperature
contours near the die exit. The high velocity gradients near the die exit lead to an important viscous
dissipation effect. The temperature of the polymer melt increases from the converging zone to the
die lip. This increase in temperature must be monitored to avoid melt degradation. The simulation
helps optimize the geometry of the die, the flow section for the cooling fluid, and other conditions in
order to maximize the flow rate and the extrudate speed.
a. Verify that you have millimeters selected as your units for length in CFD-Post.
b. Define the line for the plot with the points (0, 1, 0) and (15, 1, 0).
ii. Click OK to accept the default name (Line 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
160 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
v. Click Apply.
c. Create a plot.
ii. Click OK to accept the default name (Chart 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 161
Non-Isothermal Flow Through a Cooled Die
iii. In the General tab of the details view, ensure XY is selected for the chart Type and
disable Display Title.
iv. In the Data Series tab, select Line 1 from the Locations drop-down list for Series 1.
vi. In the Y Axis tab, select TEMPERATURE from the Variable drop-down list.
vii. With Series 1 (Line 1) enabled under the Line Display tab, select Rectangle from the
Symbols drop-down list. Retain the default setting forAutomatically generate Symbol
Color.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
162 of ANSYS, Inc. and its subsidiaries and affiliates.
Summary
3.5. Summary
In this tutorial, you solved the non-isothermal flow of a polymer melt through a cooled die. You set the
material properties for the melt and supplied suitable boundary conditions. A specific interpolation
scheme was used for the temperature in order to cope with the important gradients.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 163
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
164 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 4: 3D Extrusion
This tutorial is divided into the following sections:
4.1. Introduction
4.2. Prerequisites
4.3. Problem Description
4.4. Preparation
4.5. Setup and Solution
4.6. Summary
4.1. Introduction
This tutorial illustrates the simulation of a 3D extrusion process. Due to the velocity rearrangement that
occurs at the die exit, the shape of the extrudate is usually different from the die lip cross-section.
Polyflow is capable of handling 3D free surfaces, so it can predict the extrudate shape that corresponds
to a given die geometry under prescribed operating conditions.
4.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the setup procedure
will not be shown explicitly.
The melt enters the die as shown in Figure 4.1: Problem Description (p. 166) at a flow rate of = 10
3
cm /s (a quarter of the actual flow rate) and the extrudate is obtained at the exit. At the end of the
computational domain, it is assumed that the extrudate is fully deformed and that it will not deform
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 165
3D Extrusion
any further. It is assumed that subdomain 2 is long enough to account for all the deformation of the
extrudate.
The incompressibility and momentum equations are solved over the computational domain. The domain
for the problem is divided into two subdomains (as shown in Figure 4.1: Problem Description (p. 166))
so that the remeshing algorithm can be applied only to the portion of the mesh that will be deformed.
The subdomain 1 represents the die where the fluid is confined. The subdomain 2 corresponds to the
extrudate that is in contact with the air and can deform freely. The main aim of the calculation is to
find the location of the free surface (the skin of the extrudate).
The boundary sets for the problem are shown in Figure 4.2: Boundary Sets for the Problem (p. 167), and
the conditions at the boundaries of the domains are as follows.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
166 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
4.4. Preparation
To prepare for running this tutorial:
3. Unzip the 3d_extrusion.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 167
3D Extrusion
4.5.7. Solution
4.5.8. Postprocessing
2. Save the Ansys Workbench project using File → Save, entering 3D-extrusion as the name
of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Steady-state problem(s)
This problem is a 3D simulation of the extrusion process, that is, a three-dimensional geometry
is assumed for the die. A Cartesian (x,y,z) reference frame is used for the 3D calculations. A steady-
state condition is assumed for the problem.
Create a sub-task
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
168 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Since this problem involves a free surface, the domain is divided into two subdomains; one for the
region near the free surface (subdomain 2) and the other for the rest of the domain (subdomain 1).
In this problem, the sub-task applies to both subdomains, which is the default condition.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 169
3D Extrusion
Material Data
The Cross law exhibits shear-thinning (the decrease in viscosity as the shear rate increases) that is a
characteristic of many polymers. The viscosity in this tutorial is given by the Cross law:
(4.1)
where:
= shear rate
3. Specify the value , referred to as “fac” in the graphical user interface (compare the equation
at the top of the Cross law menu to Equation 4.1 (p. 170)).
Modify fac
Enter 85000 [units: poise] as the New value and click OK.
4. Specify the value for , referred to as “tnat” in the graphical user interface.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
170 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Modify tnat
5. Specify the value for , referred to as “expom” in the graphical user interface.
Modify expom
6. Select Upper level menu three times to return to the 3D die swell menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 171
3D Extrusion
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_1 and click Modify.
b. Click Inflow.
When the Automatic option is selected, Polydata automatically chooses the most appropriate
method to compute the inflow condition.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
172 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
2. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY_2 at the wall of
SUBDOMAIN_1 (BOUNDARY_2).
At a solid-liquid interface, the velocity of the liquid is that of the solid surface. Hence the fluid is as-
sumed to stick to the wall. This is known as the no-slip condition because the liquid is assumed to
adhere to the wall, and hence, has no velocity relative to the wall.
By default, Polydata imposes = = 0 along all boundaries. No action is required to accept the
default condition.
In 2D axisymmetric problems, Polydata automatically identifies the axis of symmetry, but for 3D flows,
you must manually identify a plane of symmetry.
The normal velocity ( ) and the tangential force ( ) are set to zero on a symmetry plane. A particle
cannot cross the plane ( = 0) due to the symmetry, so the particles flow at the same velocity on
both sides of the symmetry plane, leading to a zero tangential force.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_3 and click Modify.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_4 and click Modify.
In a steady-state problem, the velocity field must be tangential to a free surface, since no fluid particles
go out of the domain through the free surface. This constraint is called the kinematic condition,
. This equation requires an initial condition, that is, the starting line of the free surface. In the
current problem, the starting line of the free surface is the intersection of boundary 2 and boundary
5 (see Figure 4.2: Boundary Sets for the Problem (p. 167)).
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_5 and click Modify.
h. Click Upper level menu to return to the Flow boundary conditions panel.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 173
3D Extrusion
It is assumed that a uniform velocity profile is reached at the exit. The melt is not subjected to any
externally applied stress at the exit, so the condition of zero normal and tangential forces is selected.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_6 and click Modify.
c. Accept the default value of 0 for the normal force, , by clicking Upper level menu.
d. Accept the default value of 0 for the tangential force, , by clicking Upper level menu.
e. Click Upper level menu at the top of the Flow boundary conditions panel.
4.5.5. Remeshing
This model involves a free surface, whose shape is unknown a priori, which will be calculated together
with the flow equations. A portion of the mesh is affected by the relocation of this boundary, so a
remeshing technique is applied on this part of the mesh. The free surface is entirely contained within
subdomain 2, therefore only subdomain 2 is affected by the relocation of the free surface.
Global remeshing
In general, only one local remeshing is required for direct extrusion simulations. It becomes necessary
to define multiple local remeshings for inverse extrusion simulations. A single local remeshing is suf-
ficient for this case.
SUBDOMAIN_1 is moved from the top list to the bottom list, indicating that only SUBDOMAIN_2
will be remeshed.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
174 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The purpose of the remeshing technique is to relocate internal nodes according to the displacement
of boundary nodes due to the motion of the free surface, since a part of the mesh is deformed. For
3D extrusion problems where large deformations of the extrudate are expected, the optimesh
remeshing technique is recommended
The optimesh remeshing technique requires the direction of extrusion to be parallel to the , , or
axis, and all slices into which the remeshing domain is cut must be perpendicular to the extrusion
axis.
The domain to be remeshed is cut into a series of 2D slices (planes) in a direction perpendicular to
the direction of extrusion, and each plane is remeshed independently. For this process, Polyflow requires
the selection of the initial plane and the final plane. In this problem, the initial plane is the intersection
of subdomain 2 with subdomain 1, and the final plane is the intersection of subdomain 2 with the
flow exit (boundary 6).
a. Specify the initial plane for the optimesh remeshing technique, by selecting SUBDOMAIN_1
and clicking Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 175
3D Extrusion
b. Specify the final plane for the optimesh remeshing technique, by selecting BOUNDARY_6
and clicking Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
176 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Polydata asks if you want to change from the surface kinematic condition to the line kinematic
condition.
The line kinematic condition is recommended for extrusion problems, and should be used in
combination with the optimesh remeshing technique.
Polydata asks you to confirm the current system units and fields that are to be saved to the results file for
postprocessing.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 177
3D Extrusion
4. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and for the Polyflow format results file (res).
4.5.7. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
4.5.8. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
178 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1. Double-click the Results cell in the Workbench analysis and read the results files saved by Polyflow.
CFD-Post reads the solution fields that were saved to the results file.
a. Click the Insert menu and select Contour or click the button.
b. Click OK to accept the default name (Contour 1) and display the details view below the
Outline tab.
c. In the Outline tree tab, under User Locations and Plots, deselect Wireframe.
ii. Select all topological entities under SUBDOMAIN_1 and SUBDOMAIN_2 in the Location
Selector dialog box (use Shift) and click OK.
iii. Select VELOCITIES from the Variable drop-down list (or by clicking ).
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 179
3D Extrusion
In Figure 4.3: Contours of Velocity Magnitude (p. 180), the velocity is zero along the die wall (as expected)
and there is a fully developed profile at the inlet of the die. At the die outlet, the velocity profile
changes to become constant throughout the extrudate cross-section. The transition between these
two states can be seen in the beginning section of the extrudate.
a. Deselect Contour 1 in the Outline tree tab under User Locations and Plots.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
180 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Click OK to accept the default name (Plane 1) and display the details view below the
Outline tab.
iii. In the Geometry tab of the details view of Plane 1, ensure XY Plane is selected from
the Method drop-down list.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 181
3D Extrusion
vi. Repeat steps 3.b.i.–v. for the other planes, at Z = 0.08, 0.15, and 0.45 m.
vii. In the Outline tree tab, under User Locations and Plots, deselect Plane 1, Plane 2,
Plane 3, and Plane 4.
c. Click the Insert menu and select Contour or click the button.
d. Click OK to accept the default name (Contour 2) and display the details view below the
Outline tab.
e. In the Outline tree tab under User Locations and Plots, select Wireframe.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
182 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Select all planes under User Locations and Plots (use Shift for multiple selection).
iv. Select VELOCITIES from the Variable drop-down list (or click ).
Velocity profiles at the flow inlet, the flow outlet, and planes just before and just after the die exit are
displayed (Figure 4.4: Velocity Profiles at Cross-Sections (p. 184)). Compare the velocity profile within
the die to the velocity profile just after the die exit at the end of the computational domain. In the
die the flow is fully developed. The velocity profile is flat in the extrudate, far away from the die exit;
all the particles in the cross-section plane are at the same velocity. Just beyond the die exit, in the
transitional zone, the velocity profile is reorganized. The velocity profile on the plane = 15 is no
longer fully developed, but it is not yet flat either. The velocity rearrangement is the source of the
deformation of the extrudate.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 183
3D Extrusion
In the Outline tree tab, under User Locations and Plots, deselect Contour 2 and Wireframe.
ii. Click OK to accept the default name (Polyline 1) and display the details view below
the Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
184 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. In the Geometry tab of the details view, select Boundary Intersection from the
Method drop-down list.
vi. Under the Color tab, click next to Color and select dark blue.
ii. Click OK to accept the default name (Polyline 2) and display the details view below
the Outline tab.
iii. In the Geometry tab of the details view, select Boundary Intersection from the
Method drop-down list.
d. On the axis triad in the graphics window click +Z to view from Z-direction.
This allows you to compare the size and shape of the flow inlet with that of the flow outlet
without distortion due to perspective.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 185
3D Extrusion
Since the model involves a generalized Newtonian fluid, there are no viscoelastic effects. The swelling
(Figure 4.5: Swelling of the Extrudate (p. 186)) is only due to reorganization of the velocity profile at
the die exit. Fluid from the high-speed region moves to the low-speed region and pushes the free
surface to the exterior.
a. Define the line of the plot with the points (0.0, 0.1, 0.1) and (0.1, 0.0, 0.1).
ii. Click OK to accept the default name (Line 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
186 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. Enter 0, 0.1, 0.1 for Point 1 and 0.1, 0, 0.1 for Point 2.
b. Create a plot.
ii. Click OK to accept the default name (Chart 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 187
3D Extrusion
iii. In the General tab of the details view, ensure XY - Line is selected for Type, and enter
3D Extrusion for Title.
iv. In the Data Series tab, for Series 1, select Line 1 from the Locations drop-down list
vi. In the Y Axis tab, select VELOCITIES from the Variable drop-down list.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
188 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The shear-thinning introduced by the Cross law is not clearly visible in Figure 4.6: Velocity Magnitude
Along a Diagonal of Die Exit Section (p. 188) due to the large finite elements along the die wall. The
mesh should be refined in that zone.
a. Define the line of the plot with the points (0.08, 0.02, 0.00) and (0.08, 0.02, 0.50).
ii. Click OK to accept the default name (Line 2) and display the details view below the
Outline tab.
iii. Enter 0.08, 0.02, 0 for Point 1 and 0.08, 0.02, 0.5 for Point 2.
b. Create a plot.
ii. Click OK to accept the default name (Chart 2) and display the details view below the
Outline tab.
iii. In the General tab of the details view, ensure XY is selected for Type, and disable
Display Title.
iv. In the Data Series tab, select Line 2 from the Locations drop-down list for Series 1.
vi. In the Y Axis tab, select VELOCITIES X from the Variable drop-down list.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 189
3D Extrusion
Figure 4.7: X-Velocities Along a Line Close to the Center of the Die
a. In the Y Axis tab of the details of Chart 2, select VELOCITIES Y from the Variable drop-
down list.
b. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
190 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Figure 4.8: Y-Velocities Along a Line Close to the Center of the Die
a. In the Y Axis tab of the details view of Chart 2, select VELOCITIES Z from the Variable
drop-down list.
b. Click Apply.
Figure 4.9: Velocities Along a Line Close to the Center of the Die (p. 192) shows that the flow
slows down ( decreases) after the die exit. Meanwhile, particles travel from the center of
the extrudate toward the edge, creating the swelling of the extrudate. Figure 4.7: X-Velocities
Along a Line Close to the Center of the Die (p. 190) and Figure 4.8: Y-Velocities Along a Line
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 191
3D Extrusion
Close to the Center of the Die (p. 191) show that the peak values of and are located
at the very beginning of the extrudate, and vanish at the end of the free jet.
Figure 4.9: Velocities Along a Line Close to the Center of the Die
4.6. Summary
This tutorial introduced the concept of a 3D extrusion problem. You solved the problem using a specific
3D geometry for the die and made suitable assumptions about the physics of the problem. You analyzed
the factors affecting the extrudate shape. In Polydata you learned how to use the optimesh remeshing
method, which is recommended for 3D extrusion problems.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
192 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 5: Direct Extrusion
This tutorial is divided into the following sections:
5.1. Introduction
5.2. Prerequisites
5.3. Problem Description
5.4. Setup and Solution
5.5. Summary
5.6. Appendix
5.1. Introduction
This tutorial is similar to the 3D extrusion problem solved in 3D Extrusion (p. 165), where the shape of
the extrudate was computed from the die geometry. In this tutorial, a complex geometry (free surface)
is associated with the exit section of the die and undergoes large deformations during the extrusion
process. Consequently, the problem becomes highly nonlinear and special convergence techniques are
required to obtain a solution. This tutorial introduces the evolution procedure in Polyflow that is used
to handle nonlinear problems.
• Set material properties and boundary conditions for a direct extrusion problem.
5.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the set up procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 193
Direct Extrusion
The melt enters the die as shown in Figure 5.1: Problem Description (p. 194) at a flow rate = 10 cm3/s
(a quarter of the actual flow rate) and the extrudate is obtained at the exit. It is assumed that subdomain
2 is long enough to account for all the deformation of the extrudate.
The incompressibility and momentum equations are solved over the computational domain. The domain
for the problem is divided into two subdomains (as shown in Figure 5.1: Problem Description (p. 194))
so that the remeshing algorithm can be applied only to the portion of the mesh that is deformed.
Subdomain 1 represents the fluid as it enters and is confined by the die. Subdomain 2 corresponds to
the extrudate that is in contact with the air (and can deform freely). The main aim of the calculation is
to find the location of the free surface (the skin of the extrudate).
The boundary set for the problem are shown in Figure 5.2: Boundary Set for the Problem (p. 195), and
the conditions at the boundaries of the domains are:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
194 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 195
Direct Extrusion
5.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the direct_extrusion.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering direct-extrusion as the
name of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Evolution problem(s)
The complex geometry associated with the free surface of the extrudate introduces nonlinear
terms into the kinematic condition equation used to find its location. An evolution scheme
is used to handle the nonlinear problem.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
196 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Create a sub-task
Since this problem involves a free surface, the domain is divided into two subdomains; one for
the region near the free surface (SUBDOMAIN_2) and the other for the rest of the domain (SUB-
DOMAIN_1). In this problem, the sub-task applies to both subdomains, which is the default
condition.
Accept the default selection of both subdomains by clicking Upper level menu.
Material Data
The viscosity in this tutorial is given by the power law. For information on power law, see Ap-
pendix (p. 214) .
3. Specify the value for , referred to as “fac” in the graphical user interface (compare the equation
at the top of the Power law menu to the equation shown in the Appendix (p. 214)).
Modify fac
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 197
Direct Extrusion
Enter 3e+05 [units: poise] as the New value and click OK.
4. Retain the default value for , referred to as “tnat” in the graphical user interface.
Modify tnat
5. Specify the value for , referred to as “expo” in the graphical user interface.
Modify expo
6. Click Upper level menu three times to return to the Direct extrusion menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
198 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_1 and click Modify.
b. Click Inflow.
c. Ensure Volumetric flow rate is selected and click Modify volumetric flow rate.
When the Automatic option is selected, Polydata chooses the most appropriate method to
compute the inflow condition.
In 2D axisymmetric problems, the axis of symmetry is automatically identified by Polydata, but for
3D flows, you must manually identify a plane of symmetry. The normal velocity ( ) and the tangential
force ( ) are set to zero on a symmetry plane. A particle cannot cross the plane ( = 0) due to the
symmetry, so the particles flow at the same velocity on both sides of the symmetry plane, leading to
a zero tangential force.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_2 and click Modify.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_3 and click Modify.
4. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY_4 at the wall of
SUBDOMAIN_1 (BOUNDARY_4).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 199
Direct Extrusion
At a solid-liquid interface, the velocity of the liquid is that of the solid surface. Hence the velocity the
fluid is assumed to stick to the wall. This is known as the no-slip assumption because the liquid is
assumed to adhere to the wall, and hence, has no velocity relative to the wall.
By default, Polydata imposes = = 0 along all boundaries. No action is required to accept the
default condition.
In a steady-state problem, the velocity field must be tangential to a free surface, since no fluid particles
go out of the domain through the free surface. This constraint is called the kinematic condition,
= 0. This equation requires an initial condition, that is, the starting line of the free surface. In this
problem, the starting line of the free surface is the intersection of boundary 4 and boundary 5 (see
Figure 5.2: Boundary Set for the Problem (p. 195)).
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_5 and click Modify.
d. Select No condition along BOUNDARY_4 (the boundary where the free surface starts) and
click Modify.
f. Click Upper level menu to return to the Boundary conditions on the moving surface
menu.
g. Click Upper level menu two times to return to the Flow boundary conditions menu.
It is assumed that a uniform velocity profile is reached at the exit. The melt is not subjected to any
externally applied stress at the exit, so the condition of zero normal and tangential forces is selected.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_6 and click Modify.
c. Accept the default value of 0 for the normal force, , by clicking Upper level menu.
d. Accept the default value of 0 for the tangential force, , by clicking Upper level menu.
5.4.6. Remeshing
This model involves a free surface, whose shape is unknown a priori, which will be calculated together
with the flow equations. A portion of the mesh is affected by the relocation of this boundary, so a
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
200 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
remeshing technique is applied on this part of the mesh. The free surface is entirely contained within
subdomain 2, therefore only subdomain 2 is affected by the relocation of the free surface.
Global remeshing
In general, only one local remeshing is required for direct extrusion simulations. It becomes necessary
to define multiple local remeshings for inverse extrusion simulations. A single local remeshing is suf-
ficient for this case.
SUBDOMAIN_1 is moved from the top list to the bottom list, indicating that only SUBDOMAIN_2
will be remeshed.
The purpose of the remeshing technique is to relocate internal nodes according to the displacement
of boundary nodes due to the motion of the free surface, since a part of the mesh is deformed. For
3D extrusion problems where large deformations of the extrudate are expected, the optimesh
remeshing technique is recommended. For information on optimesh remeshing technique, refer to
the Appendix (p. 214).
a. Specify the initial plane for the optimesh remeshing technique, by selecting SUBDOMAIN_1
and clicking Modify.
b. Specify the final plane for remeshing technique, by selecting BOUNDARY_6 and clicking
Modify.
Polydata asks if you want to change from the surface kinematic condition to the line kinematic
condition.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 201
Direct Extrusion
The line kinematic condition is recommended for extrusion problems, and must be used in
combination with the optimesh remeshing technique.
d. Click Accept the current setup in the Element distortion check menu.
In complex extrusion simulations, the finite element mesh can undergo great deformations. The
Element distortion check menu deals with the detection of all possible distortions of the elements.
Numerical parameters
a. Define the starting solution for the iterative scheme in the calculation of the free surface
location.
The first calculation is performed at . Increase the value of the initial increment of ( ) to
reduce the number of evolution steps and to speed up the calculation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
202 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
3. Click Upper level menu two times to return to the top-level Polydata menu.
5.4.8. Outputs
After Polyflow calculates a solution, it can save the results in several different formats. Choose the one
that is appropriate for your postprocessor. In this case, save the outputs in IGES format, as well as the
default format for CFD-Post.
Outputs
1. Retain the default output (CFD-Post) and click Enable Iges file output.
The default CFD-Post output is used for postprocessing with CFD-Post. The IGES output contains the
modified geometry of the extrudate (after remeshing) calculated at every step of the evolution pro-
cedure. For information on IGES output, see Appendix (p. 214).
Polydata asks you to confirm the current system units and fields that are to be saved to the results
file for postprocessing.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 203
Direct Extrusion
Polydata asks you to confirm fields that are to be saved to the results file for postprocessing.
3. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and for the Polyflow format results file (res).
5.4.10. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
5.4.11. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
1. Double-click the Results tab in the Polyflow analysis system. This will start CFD-Post and read
the results files saved by Polyflow. CFD-Post reads the mesh information and the solution fields
that were saved to the results file.
Deselect Wireframe in the Outline tree tab, under User Locations and Plots.
a. Click the Insert menu and select Contour or click the button.
b. Click OK to accept the default name (Contour 1) and display the details view below the
Outline tab.
c. Perform the following steps in the Geometry tab of the details view of Contour 1:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
204 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Select all topological entities under SUBDOMAIN_1 and SUBDOMAIN_2 in the Location
Selector dialog box (use Shift for multiple selection) and click OK.
iii. Select VELOCITIES from the Variable drop-down list (or by clicking ).
d. Rotate the image so that you can see the fluid at the inlet of the die, as shown in Fig-
ure 5.3: Contours of Velocity Magnitude (p. 206).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 205
Direct Extrusion
Observe that the velocity is zero along the die wall, as expected, and there is a fully developed profile
at the inlet of the die. At the die outlet, the velocity profile changes to become constant throughout
the extrudate cross-section. The transition between these two states can be seen in the first third of
the extrudate.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
206 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In the Outline tree tab, under User Locations and Plots, deselect Contour 1.
ii. Click OK to accept the default name (Plane 1) and display the details view below the
Outline tab.
iii. In the Geometry tab of the details view, ensure XY Plane is selected from the Method
drop-down list.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 207
Direct Extrusion
v. Click Apply.
vi. Repeat steps 3.b.i.–v. for the other planes, at Z = 0.03, 0.07, and 0.1999 m.
vii. In the Outline tree tab, under User Locations and Plots, deselect Plane 1, Plane 2,
Plane 3, and Plane 4.
i. Click the Insert menu and select Contour or click the button.
ii. Click OK to accept the default name (Contour 2) and display the details view below
the Outline tab.
iii. In the Outline tree tab under User Locations and Plots, select Wireframe.
iv. In the Geometry tab of the details view of Contour 2, click the button next to
Locations.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
208 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
v. Select all planes under User Locations and Plots (use Shift for multiple selection).
vii. Select VELOCITIES from the Variable drop-down list (or click ).
Figure 5.4: Velocity Profiles at Cross-Sections (p. 210) shows the velocity profiles at the flow inlet, the
flow outlet, and at the planes just before and just after the die exit.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 209
Direct Extrusion
Compare the velocity profile within the die to the velocity profile just after the die exit at the end of
the computational domain. In the die the flow is fully developed. In the extrudate, far away from the
die exit, the velocity profile is flat. That is, all the particles in a cross-sectional plane are at the same
speed. Just after the die exit, there is a transitional zone where the velocity profile is reorganized. The
velocity profile on the plane Z = 7 cm is no longer fully developed, but it is not yet flat either. The
velocity rearrangement is the source of the deformation of the extrudate.
In the Outline tree tab, under User Locations and Plots, deselect Contour 2 and Wireframe.
ii. Click OK to accept the default name (Polyline 1) and display the details view below
the Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
210 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. In Geometry tab of the details view, select Boundary Intersection from the Method
drop-down list.
vi. In the Color tab, click next to Color and select dark blue.
ii. Click OK to accept the default name (Polyline 2) and display the details view below
the Outline tab.
iii. In Geometry tab of the details view, select Boundary Intersection from the Method
drop-down list.
d. Right-click in the graphic window and select View From +Z under Predefined Camera.
This allows you to compare the size and shape of the flow inlet with that of the flow outlet
without distortion due to perspective.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 211
Direct Extrusion
a. Click the Insert menu and select Instance Transform, or click the button.
b. Click OK to accept the default name (Instance Transform 1) and display the details view
below the Outline tab.
iv. Ensure Method is set to Principal Axis and Z is selected from the Axis drop-down list.
d. In the Outline tree tab, under User Locations and Plots, right-click Polyline 1 and click
Edit (or double-click Polyline 1).
e. In the View tab, scroll down and enable Apply Instancing Transform.
g. Click Apply.
h. In the Outline tree tab, under User Locations and Plots, right-click Polyline 2 and click
Edit (or double-click Polyline 2).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
212 of ANSYS, Inc. and its subsidiaries and affiliates.
Summary
You can use the central-mouse button to zoom in and out. This allows you to compare the size and
shape of the flow inlet with that of the flow outlet without distortion due to perspective.
You can also click the fit view button ( ) to properly center the image.
The deformations come from the rearrangement of the velocity profile. Particles coming from
high-speed regions in the die must slow down, while particles coming from low-speed regions
must accelerate. Observe that the central part of the cross, where the fluid easily flows in the
die, is enlarged in the extrudate, while the extremities of the branches are smaller in the extrudate.
Since the combined effect of cross-sectional expansions and reductions is very difficult to guess,
a numerical simulation is necessary for a moderate to high complexity die.
5.5. Summary
This tutorial introduced the concept of a direct extrusion problem. You solved the problem using a
specific 3D geometry for the die, made suitable assumptions about the physics of the problem, and
analyzed the factors affecting the extrudate shape. The nonlinear problem was solved using an evolution
technique to reach the convergence.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 213
Direct Extrusion
5.6. Appendix
The appendix contains the following sections:
5.6.1. Power Law
5.6.2. Optimesh Remeshing Technique
5.6.3. Evolution Scheme
5.6.4. IGES Output
where:
= consistency factor
= power-law index
= natural time
The domain to be remeshed is cut into a series of 2D slices (planes) in a direction perpendicular to
the direction of extrusion, and each plane will be remeshed independently. For this process, Polyflow
requires the selection of the initial plane and the final plane. In this problem, the initial plane is the
intersection of SUBDOMAIN_2 with SUBDOMAIN_1, and the final plane is the intersection of SUBDO-
MAIN_2 with the flow exit (boundary 6).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
214 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 215
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
216 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 6: Inverse Extrusion
This tutorial is divided into the following sections:
6.1. Introduction
6.2. Prerequisites
6.3. Problem Description
6.4. Setup and Solution
6.5. Summary
6.6. Appendix
6.1. Introduction
Inverse extrusion deals with the computation of the shape of a die that produces an extrudate of the
desired shape. This tutorial illustrates how to handle a complex inverse extrusion problem. In this tutorial,
slip conditions along the die walls are considered and evolution on the slip coefficient is enabled to
aid convergence.
• Set material properties and boundary conditions for an inverse extrusion problem.
6.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the set up procedure
will not be shown explicitly.
The melt enters the die as shown in Figure 6.1: Problem Description (p. 218) at a flow rate = 10 cm3/s
(a quarter of the actual flow rate) and the extrudate is obtained at the exit.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 217
Inverse Extrusion
The incompressibility and momentum equations are solved over the computational domain. The domain
for the problem is divided into two sub-domains (as shown in Figure 6.1: Problem Description (p. 218))
so that specific remeshing algorithms can be applied in each sub-domain to accurately predict the die
profile. Subdomain 1 represents the fluid as it enters and is confined by the die. Subdomain 2 corresponds
to the extrudate that is in contact with the air (and can deform freely). The main aim of the calculation
is to compute the geometry of the die to obtain the desired extrudate.
The boundary set for the problem are shown in Figure 6.2: Boundary Set for the Problem (p. 218), and
the conditions at the boundaries of the domains are given below
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
218 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
6.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the inverse_extrusion.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 219
Inverse Extrusion
2. Save the Ansys Workbench project using File → Save, entering inverse-extrusion as the
name of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Evolution problem(s)
Apply the evolution scheme on the slip coefficient along the outer wall of the die (boundary 4)
when you define the slip boundary conditions.
Create a sub-task
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
220 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
This problem involves a free surface, so the domain is divided into two sub-domains; one for the region
near the free surface (SUBDOMAIN_2) and the other for the rest of the domain (SUBDOMAIN_1). In
this problem, the sub-task applies to both sub-domains, which is the default condition.
Accept the default selection of both subdomains by clicking Upper level menu.
Material Data
The viscosity in this tutorial is given by the power law. For information on power law, see Power
Law (p. 240) .
3. Specify the value of , referred to as “fac” in the graphical user interface (compare the equation
at the top of the Power law menu to the equation shown in the Power Law (p. 240)).
Modify fac
Enter 300000 [units: poise] as the New value and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 221
Inverse Extrusion
4. Retain the default value for , referred to as “tnat” in the graphical user interface.
Modify tnat
5. Specify the value for , referred to as “expo” in the graphical user interface.
Modify expo
6. Click Upper level menu three times to return to the Inverse Extrusion menu.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_1 and click Modify.
b. Click Inflow.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
222 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
When the Automatic option is selected, Polydata chooses the most appropriate method to
compute the inflow condition.
In 2D axisymmetric problems, the axis of symmetry is automatically identified by Polydata, but for
3D flows, you must manually identify a plane of symmetry.
The normal velocity ( ) and the tangential force ( ) are set to zero on a symmetry plane. A particle
cannot cross the plane ( = 0) due to the symmetry, so the particles flow at the same velocity on
both sides of the symmetry plane, leading to a zero tangential force.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_2 and click Modify.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_3 and click Modify.
4. Set the conditions along the outer wall of the die (BOUNDARY_4).
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_4 and click Modify.
The evolution scheme is used to aid with convergence by starting with a low value for the slip
coefficient and slowly increasing the value of the coefficient to reach a no-slip condition. With
a low value for the slip coefficient there is no swelling of the extrudate, simplifying the calculation.
As the slip coefficient increases, the extrudate begins to swell because the fluid in contact with
the wall slows down, which increases the velocity of the fluid in the center of the die. For more
information on the evolution scheme, see Evolution Scheme (p. 240).
i. Click EVOL at the top of the Polydata menu to enable the evolution inputs for the slip
coefficient.
Modify k
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 223
Inverse Extrusion
You will retain the default values for a and b, and will modify the values for c and d.
xi. Click EVOL at the top of the Polydata menu to disable the evolution inputs.
xiv. Click Upper level menu two times to return to the Flow boundary conditions menu.
In a steady-state problem, the velocity field must be tangential to a free surface, since no fluid particles
go out of the domain through the free surface. This constraint is called the kinematic condition,
= 0. This equation requires an initial condition (the starting line of the free surface). In this problem,
the starting line of the free surface is the intersection of boundary 4 and boundary 5 (see Fig-
ure 6.2: Boundary Set for the Problem (p. 218)).
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_5 and click Modify.
d. Select No condition along BOUNDARY_4 (the boundary where the free surface starts) and
click Modify.
f. Click Upper level menu to return to the Boundary conditions on the moving surface
menu.
h. Click Outlet (Inv. prediction) to define the outlet of the moving surface.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
224 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In inverse extrusion problems, you have to predict the appropriate die exit cross-section to obtain
a given extrudate cross-section. By defining the outlet of a free surface, inform Polyflow the desired
extrudate shape. Hence, you impose the outlet of the moving surface along the last section of
the free jet as the outlet of the free surface. This section will not be modified during the calculation.
i. Select BOUNDARY_6 as the outlet of the moving surface and click Confirm.
j. Click Upper level menu to return to the Flow boundary conditions menu.
It is assumed that a uniform velocity profile is reached at the exit. The melt is not subjected to any
externally applied stress at the exit, so the condition of zero normal and tangential forces is selected.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY_6 and click Modify.
c. Accept the default value of 0 for the normal force, , by clicking Upper level menu.
d. Accept the default value of 0 for the tangential force, , by clicking Upper level menu.
7. Click Upper level menu at the top of the Flow boundary conditions menu.
6.4.6. Remeshing
The purpose of the remeshing technique is to relocate internal nodes according to the displacement of
boundary nodes due to the motion of the free surface, since a part of the mesh is deformed. For information
on remeshing technique, see Appendix (p. 240).
Global remeshing
SUBDOMAIN_1 is moved from the top list to the bottom list, indicating that only SUBDOMAIN_2
will be remeshed.
For 3D extrusion problems where large deformations of the extrudate are expected, the optimesh
remeshing technique is recommended. For information on optimesh remeshing technique see
Appendix (p. 240).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 225
Inverse Extrusion
a. Specify the initial plane for the optimesh remeshing technique, by selecting SUBDOMAIN_1
and clicking Modify.
b. Specify the final plane for the remeshing technique, by selecting BOUNDARY_6 and clicking
Modify.
c. Polydata asks if you want to change from the surface kinematic condition to the line kinematic
condition.
The line kinematic condition is recommended for extrusion problems, and must be used in
combination with the optimesh remeshing technique.
e. Click Accept the current setup in the Element distortion check menu.
In complex extrusion simulations, the finite element mesh can undergo great deformations. The
Element distortion check menu deals with the detection of all possible distortions of the elements.
Accept default options.
The technique of inverse prediction is selected to calculate the profile for the “constant section”
region of the die.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
226 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Specify the inlet plane of the system of planes, by selecting BOUNDARY_1 and clicking
Modify.
b. Specify the outlet plane for the system of planes, by selecting SUBDOMAIN_2 and
clicking Modify.
Numerical parameters
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 227
Inverse Extrusion
Note:
Setting the final value of equal to 20 creates a large enough slip coefficient
that it is equivalent to a no-slip condition at the die wall (BOUNDARY_4, as dis-
cussed in a previous step).
2. Click Upper level menu three times to return to the top-level Polydata menu.
6.4.8. Outputs
After Polyflow calculates a solution, it can save the results in several different formats. Choose the one
that is appropriate for your postprocessor. In this case, save the outputs in IGES format, as well as the
default format for CFD-Post.
Outputs
1. Retain the default output (CFD-Post) and click Enable Iges file output.
The default CFD-Post output is used for postprocessing with CFD-Post. The IGES output contains the
modified geometry of the extrudate (after remeshing) calculated at every step of the evolution pro-
cedure. For information on IGES output, see Appendix (p. 240).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
228 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Polydata asks you to confirm the current system units and fields that are to be saved to the results
file for postprocessing.
3. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and for the Polyflow format results file (res).
6.4.10. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
A cfx.res file (corresponding to the eight evolution steps of the flow case) will be created.
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 229
Inverse Extrusion
6.4.11. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
1. Double-click the Results tab in the Workbench analysis and read the results files saved by Polyflow.
CFD-Post reads the solution fields that were saved to the results file.
a. Deselect Wireframe in the Outline tree tab, under User Locations and Plots.
b. Click the Insert menu and select Contour or click the button.
c. Click OK to accept the default name (Contour 1) and open the details view below the
Outline tab.
d. Perform the following steps In the Geometry tab of the details view.
ii. Select all topological entities under SUBDOMAIN_1 and SUBDOMAIN_2 in the Location
Selector dialog box (use Shift for multiple selection) and click OK.
iii. Select VELOCITIES from the Variable drop-down list (or by clicking ).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
230 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
You can see in Figure 6.3: Contours of Velocity Magnitude (p. 231) that the velocity is zero along the
die wall, as expected, and there is a fully developed profile at the inlet of the die. At the die outlet,
the velocity profile changes to become constant throughout the extrudate cross-section. The transition
between these two states can be seen at the beginning of extrudate section.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 231
Inverse Extrusion
In the Outline tree tab, under User Locations and Plots deselect Contour 1.
ii. Click OK to accept the default name (Plane 1) and display the details view below the
Outline tab.
iii. In the Geometry tab of the details view, select XY Plane from the Method drop-down
list.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
232 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
vi. Repeat steps 3.b.i.–v. to create the other planes at Z = 0.03, 0.07, and 0.1999 m.
vii. In the Outline tree tab, under User Locations and Plots, deselect Plane 1, Plane 2,
Plane 3, and Plane 4.
c. Click the Insert menu and select Contour or click the button.
d. Click OK to accept the default name (Contour 2) and open the details view below the
Outline tree tab.
e. In the Outline tree tab under User Locations and Plots, select Wireframe.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 233
Inverse Extrusion
ii. Select all planes under User Locations and Plots (use Shift for multiple selection).
iv. Select VELOCITIES from the Variable drop-down list (or click ).
The velocity profiles planes are located at the flow inlet, the flow outlet, and planes just before and
after the die exit as shown in Figure 6.4: Velocity Profile Planes (p. 235).
Compare the velocity profile within the die to the velocity profile just after the die exit at the end of
the computational domain.
• The velocity profile is flat in the extrudate, far away from the die exit. All particles in the cross-section
plane are at the same velocity.
• Just after the die exit, there is a transitional zone where the velocity profile is reorganized.
• The velocity profile on the plane = 7 cm is not fully developed, but it is not flat either.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
234 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In the Outline tree tab, under User Locations and Plots, deselect Contour 2 and Wireframe.
ii. Click OK to accept the default name (Polyline 1) and display the details view below
the Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 235
Inverse Extrusion
iii. In the Geometry tab of the details view, select Boundary Intersection from the
Method drop-down list.
vi. In the Color tab, click next to Color and select dark blue.
ii. Click OK to accept the default name (Polyline 2) and display the details view below
the Outline tab.
iii. In the Geometry tab of the details view, select Boundary Intersection from the
Method drop-down list.
d. Right-click in the graphic window and select View From +Z under Predefined Camera.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
236 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Click the Insert menu and select Instance Transform, or click the button.
b. Click OK to accept the default name (Instance Transform 1) and display the details view
below the Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 237
Inverse Extrusion
c. Perform the following steps in the Definition tab of the details view of Instance Transform
1:
iv. Ensure Principal Axis and Z are selected for Method and Axis in the Axis Definition
group box.
d. In the Outline tree tab, under User Locations and Plots, right-click Polyline 1 and click
Edit (or double-click Polyline 1).
e. In the View tab of the details view, scroll down and enable Apply Instancing Transform.
g. Click Apply.
h. In the Outline tree tab, under User Locations and Plots, right-click Polyline 2 and click
Edit (or double-click Polyline 2).
You can use the central-mouse button to zoom in and out. This allows you to compare the size and
shape of the flow inlet with that of the flow outlet without distortion due to perspective.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
238 of ANSYS, Inc. and its subsidiaries and affiliates.
Summary
You can also click the fit view button ( ) to properly center the image.
The deformation of the extrudate is the result of the rearrangement taking place at the die exit.
Particles coming from high-speed regions in the die must slow down, while particles coming from
low-speed regions must accelerate. You can change the speed by enlarging the flowing section. A
tube of fluid at high speed in the die will enlarge its cross-section in the extrudate to decrease its av-
erage velocity. A tube of fluid at low speed in the die will reduce its cross-section in the extrudate in
order to increase average speed. In Figure 6.5: Swelling of the Extrudate (p. 239), you can see that the
die design tool compensated for these effects. The central part of the cross-section where the fluid
easily flowed in the original die has been reduced by the die design tool, while the extremities of the
branches were enlarged because the flow was much slower in the original die. Since the combined
effects of the cross-sectional enlargements and reductions are very difficult to guess, the numerical
simulation is necessary to help the die designer reduce the number of trial-and-error iterations.
6.5. Summary
This tutorial introduced the concept of an inverse extrusion problem. You solved the problem assuming
suitable conditions for the physics of the problem and Polyflow predicted the shape of the die. You
used multiple domain calculations with remeshing methods most suited to 3D inverse extrusion problems.
The nonlinear problem was solved using an evolution technique to aid convergence.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 239
Inverse Extrusion
6.6. Appendix
The appendix contains the following sections:
6.6.1. Power Law
6.6.2. Evolution Scheme
6.6.3. Remeshing Technique
6.6.4. Optimesh Remeshing Technique
6.6.5. IGES Output
where:
= consistency factor
= power-law index
= natural time
where:
= material parameters
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
240 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
= material parameters
This technique modifies the location of the section where the boundary conditions on the kinematic
condition apply (the die-lip region). Apply local remeshing technique to the region between the entry
section and the die-lip area (subdomain 1). In this tutorial, you define a "constant section" on this
subdomain. This means that the die cross-section is constant from the die entry to the die exit (a
parallel die). Using this two-stage remeshing technique, Polyflow calculates the die profile that produces
an extrudate of the desired shape. More complex deformations of the die are available via the
definition of different local remeshings within the die.
The domain to be remeshed is cut into a series of 2D slices (planes) in a direction perpendicular to
the direction of extrusion, and each plane will be remeshed independently. For this process, Polyflow
requires the selection of the initial plane and the final plane. In this problem, the initial plane is the
intersection of subdomain 2 with subdomain 1, and the final plane is the intersection of subdomain
2 with the flow exit (boundary 6).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 241
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
242 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 7: Flow of Two Immiscible Fluids
This tutorial is divided into the following sections:
7.1. Introduction
7.2. Prerequisites
7.3. Problem Description
7.4. Setup and Solution
7.5. Summary
7.1. Introduction
This tutorial examines the flow of two fluids in a single die. Two polymer melts with distinct physical
properties are fed through different channels into a die. The aim of the calculation is to predict the
location of the interface between the two fluids.
• Set material properties and boundary conditions for a moving interface problem.
7.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the set up procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 243
Flow of Two Immiscible Fluids
Incompressibility and momentum equations are solved in the fluid domains. To solve the fully coupled
problem, two sub-tasks are defined one each for fluid 1 (sub-task 1) and fluid 2 (sub-task 2). Each sub-
task will contain a particular model, domain of definition, material properties, and boundary conditions,
including the moving interface along the intersection of the two sub-tasks.
The domain of definition for the problem is divided into four subdomains: sub-task 1 is defined on
subdomain 1 and subdomain 2, and sub-task 2 is defined on subdomain 3 and subdomain 4 (see Fig-
ure 7.2: Boundary Sets and Subdomains for the Problem (p. 245)). Each sub-task is defined over two
subdomains to allow for the definition of the remeshing method only where it is necessary, (in the area
near the moving interface).
Fluid 1 has a viscosity of = 10000 poise, and fluid 2 has a viscosity of = 5000 poise.
The boundary sets for the problem are shown in Figure 7.2: Boundary Sets and Subdomains for the
Problem (p. 245).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
244 of ANSYS, Inc. and its subsidiaries and affiliates.
Problem Description
An interface is defined at the intersection of subdomain 2 and subdomain 4. In this problem, the interface
is a moving one, since the exact line of separation between the fluids is unknown.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 245
Flow of Two Immiscible Fluids
7.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the two_fluids.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering two-fluids as the name of
the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
246 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
• F.E.M. task
• Steady-state problem(s)
• 2D axisymmetric geometry
Create a sub-task
A small dialog box appears asking for the title of the problem.
Sub-task 1 is defined for SUBDOMAIN2 (the region of fluid 1 near the moving interface) and SUBDO-
MAIN1 (the rest of fluid 1), as shown in Figure 7.2: Boundary Sets and Subdomains for the Prob-
lem (p. 245).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 247
Flow of Two Immiscible Fluids
Material Data
Polydata indicates which material properties are relevant for your sub-task by graying out the irrelevant
properties. In this case, viscosity, density, inertia terms, and gravity are available for specification. For
this model, define only the viscosity of the material.
5. Specify the value for , referred to as “fac” in the graphical user interface.
Modify fac
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
248 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Enter 10000 [units: poise] as the New value and click OK.
6. Click Upper level menu three times to return to the fluid 1 menu.
7. Specify the flow boundary conditions for fluid 1 (SUBDOMAIN1 and SUBDOMAIN2).
The boundary uses the interface condition, which is the standard boundary condition between
two adjacent fluids. This condition establishes the continuity of the velocity field and the contact
forces in the momentum equation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 249
Flow of Two Immiscible Fluids
The position of the line separating the two fluids is unknown at the start of the problem and is
calculated as part of the solution, so the intersection will be defined as a moving interface. In
steady flows and problems involving immiscible fluids, the interface must be a streamline. To
satisfy this condition and to obtain the exact location of the line of separation, an additional
equation, the kinematic condition, ( = 0), is added to the system. This guarantees that the
material points do not cross the interface.
i. Select Zero wall velocity (vn=vs=0) along SUBDOMAIN4 and click Modify.
Polydata asks you to select the boundary or subdomain on which the position of the moving
surface is to be imposed.
viii. Click Upper level menu twice to return to the Kinematic condition menu.
ix. Select Upwinding in the kinematic equation. Click Upper level menu.
x. Click Accept the current setup to return to the Flow boundary conditions menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
250 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY1 and click Modify.
iv. Enter 3 [units: cm3/s] as the New value and click OK.
When this option is selected, Polydata automatically chooses the most appropriate method
to compute the inflow condition.
vi. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY2 at the
outer wall of SUBDOMAIN1 (BOUNDARY2).
At a solid-liquid interface, the velocity of the liquid is that of the solid surface. Hence the
fluid is assumed to stick to the wall. This is known as the no-slip assumption because the
liquid is assumed to adhere to the wall, and therefore has no velocity relative to the wall.
It is assumed that a fully developed velocity profile is reached at the exit, so the outflow condition
is appropriate. This condition essentially imposes a zero normal force ( ) that includes a pressure
term, and zero tangential velocity ( ).
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY5 and click Modify.
d. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY6 at the outer
wall common to SUBDOMAIN1 and SUBDOMAIN2 (BOUNDARY6).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 251
Flow of Two Immiscible Fluids
The fluid is assumed to stick to the wall, since at a solid-liquid interface the velocity of the liquid
is that of the solid surface.
This model involves a free surface, whose shape is unknown a priori, which will be calculated together
with the flow equations. A portion of the mesh is affected by the relocation of this boundary. Hence
a remeshing technique is applied on this part of the mesh. The moving interface is entirely contained
within SUBDOMAIN2, and hence only SUBDOMAIN2 will be affected by the relocation of the moving
interface.
Global remeshing
If you have a complex geometry, it may be necessary to split it into additional subdomains in
order to define a specific remeshing method on each of them.
For this purpose, Polydata allows you to create several local remeshings. For this problem, a
single local remeshing is sufficient.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
252 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The purpose of the remeshing technique is to relocate internal nodes according to the displace-
ment of boundary nodes due to the motion of the interface. Mesh nodes must be organized
along lines of remeshing (spines), which are collections of nodes logically arranged in a one-di-
mensional manner. Polydata requires the specification of the first and last spines that the fluid
encounters (inlet of spines and outlet of spines, respectively). In this case, the inlet of spines is
the intersection of SUBDOMAIN2 with SUBDOMAIN1, and the outlet of spines is the intersection
of SUBDOMAIN2 with the flow exit (BOUNDARY5).
Method of Spines
i. To specify the inlet for the system of spines, select SUBDOMAIN1 and click Modify.
ii. Specify the outlet for the system of spines, select BOUNDARY5 and click Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 253
Flow of Two Immiscible Fluids
Create a sub-task
a. Polydata asks you if you want to copy data from an existing sub-task.
b. Click No, since this sub-task has different parameters associated with it.
A small dialog box appears asking for the title of the problem.
Material Data
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
254 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
c. Specify the value for , referred to as “fac” in the graphical user interface.
Modify fac
Enter 5000 [units: poise] as the New value and click OK.
d. Click Upper level menu three times to return to the fluid 2 menu.
4. Specify the flow boundary conditions for fluid 2 (SUBDOMAIN3 and SUBDOMAIN4).
i. Select Zero wall velocity (vn=vs=0) along SUBDOMAIN2 and click Modify.
The interface condition was defined as a moving interface when setting the boundary
conditions for fluid 1. So further inputs are not required to define the moving interface for
fluid 2. Surface tension effects are neglected in this problem.
b. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY2 at the outer
wall of SUBDOMAIN3 (BOUNDARY2).
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY3 and click Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 255
Flow of Two Immiscible Fluids
iv. Accept the default value of 1 [units: cm3/s] for the flow rate by clicking OK.
d. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY4 at the outer
wall common to SUBDOMAIN3 and SUBDOMAIN4 (BOUNDARY4).
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY5 and click Modify.
f. Click Upper level menu to return to the Flow boundary conditions menu.
Global remeshing
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
256 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In this case, the inlet of spines is the intersection of SUBDOMAIN3 with SUBDOMAIN4, and the
outlet of spines is the intersection of SUBDOMAIN4 with the flow exit (BOUNDARY5).
Method of Spines
i. Specify the inlet for the system of spines by selecting SUBDOMAIN3 and click Modify.
ii. Specify the outlet for the system of spines by selecting BOUNDARY5 and click Modify.
Polydata asks you to confirm the current system units and fields that are to be saved to the results file for
postprocessing.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 257
Flow of Two Immiscible Fluids
4. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and for the Polyflow format results file (res).
7.4.7. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
258 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
7.4.8. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
1. Double-click the Results tab in the Workbench analysis and read the results files saved by Polyflow.
CFD-Post reads the solution fields that were saved to the results file.
a. In the graphical window, right-click, and select the option Predefined Camera.
a. Click the Insert menu and select Contour or click the button.
b. Click OK to accept the default name (Contour 1) and display the details view below the
Outline tab.
ii. In the Location Selector dialog box that opens, select SUBDOMAIN1_surf, SUBDO-
MAIN2_surf, SUBDOMAIN3_surf, and SUBDOMAIN4_surf (use Ctrl for multiple selec-
tion) and then click OK.
iii. Select VELOCITIES from the Variable drop-down list (or by clicking ).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 259
Flow of Two Immiscible Fluids
The velocity is much larger at the inlet of fluid 1 than at the inlet of fluid 2. There are two reasons for
this:
• The flow rate is three times larger for fluid 1 than for fluid 2.
• You are modeling an annular die. Hence the flow section is smaller for the interior channel than
for the exterior channel.
When the two fluids come into contact with each other, the interface between the two fluids is pushed
towards the exterior of the annular die.
• The die is annular, so even identical flow rates cause the interface to move in order to equilibrate
the flow sections.
• The viscosity of fluid 1 is higher than the viscosity of fluid 2. In the process of giving more room to
the most viscous fluid, its shearing decreases. This leads to a smaller global dissipation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
260 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In the Outline tree tab, under User Locations and Plots deselect Contour 1.
b. Click the Insert menu and select Vector or click the button.
c. Click OK to accept the default name (Vector 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 261
Flow of Two Immiscible Fluids
ii. In the Location Selector dialog box that opens, select SUBDOMAIN1_surf, SUBDO-
MAIN2_surf, SUBDOMAIN3_surf, and SUBDOMAIN4_surf (use Ctrl for multiple selec-
tion) and click OK.
iv. In the Symbol tab, set Symbol to Arrow3D and increase the Symbol Size to 3.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
262 of ANSYS, Inc. and its subsidiaries and affiliates.
Summary
You can see that the velocity is continuous across the interface. As both the fluids are Newtonian, the
velocity profile is a parabola on both sides of the interface. Since the force must be continuous across
the interface, the shear stress generated within fluid 1 is equal to the shear stress generated within
fluid 2 along the interface.
7.5. Summary
This tutorial introduced the concept of fluid layers flowing in the same duct. In Polydata, you learned
how to set up a multiple-domain calculation, including the definition of a moving interface and associated
remeshing methods.
The location of the interface depends largely on the physical properties of the fluids involved, the
geometry of the channels, and the operating conditions (for example: flow rates of the fluids). A CFD
simulation with Polyflow allows you to test different setups (for example: in order to optimize the
feeding of a coextrusion die).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 263
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
264 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 8: Flow of Two Immiscible Fluids by Species
Method
This tutorial is divided into the following sections:
8.1. Introduction
8.2. Prerequisites
8.3. Problem Description
8.4. Setup and Solution
8.5. Summary
8.1. Introduction
This tutorial examines the flow of two fluids in a single die. Two polymer melts with distinct physical
properties are fed through different channels into a die. The aim of the calculation is to predict the
location of the interface between the two fluids.
• Define a species.
8.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the set up procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 265
Flow of Two Immiscible Fluids by Species Method
Polyflow. The location of the interface depends on the physical properties of the fluids, the flow rates
of the fluids, and the geometry of the die.
Incompressibility and momentum equations are solved in the fluid domains. To determine the interface,
an extra scalar transport equation is solved and material properties are made functions of this scalar
using PMAT. If the scalar value is greater than 0.5, material properties of first fluid are used and if the
scalar value is less than 0.5, material properties of second species are used. A Scalar value of 0.5 determ-
ines the location of interface.
Note that the same problem has been solved using the interface tracking method (see Flow of Two
Immiscible Fluids (p. 243)).
The advantage of this method over the interface tracking method is that it can be used for more complex
geometries, and it is less computationally expensive than the interface tracking method (no remeshing
method must be defined). However this comes at a loss of accuracy. The interface tracking method
gives a very accurate position of interface, whereas the species method produces a blurred interface.
The geometry and mesh from Flow of Two Immiscible Fluids (p. 243) is used.
Fluid 1 has a viscosity of = 10000 poise, and fluid 2 has a viscosity of = 5000 poise.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
266 of ANSYS, Inc. and its subsidiaries and affiliates.
Problem Description
The boundary sets for the problem are shown in Figure 8.2: Boundary Sets and Subdomains for the
Problem (p. 268).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 267
Flow of Two Immiscible Fluids by Species Method
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
268 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Note that when using this method for a sharp interface, you should ensure that the scalar doesn't diffuse
much into the domain. To ensure this, an evolution is applied on scalar diffusivity starting from a large
value and gradually decreasing it to a very small number.
8.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the two_fluids_species.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering two-fluids-species as the
name of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 269
Flow of Two Immiscible Fluids by Species Method
• F.E.M. task
• Evolution problem(s)
• 2D axisymmetric geometry
1. Create a species A.
Define species
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
270 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Create a sub-task
b. Click SpeciesA.
Material data
Polydata indicates the material properties that are relevant for the sub-task by graying out the irrel-
evant properties. For this model, you will only define the diffusivity of the species. The evolution will
be applied on diffusivity with an initial high value (1) and decreases it to a small value (1e-9).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 271
Flow of Two Immiscible Fluids by Species Method
a. Click Diffusivity.
b. Click EVOL button at the top of Polydata menu to enable evolution inputs.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
272 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Polydata will take you to evolution panel. Here you will make species diffusivity a function of
evolution parameter . Since the diffusivity must be decreased by several orders of magnitude,
is selected.
h. Click the EVOL button at the top of the Polydata menu to disable evolution inputs.
i. Click Upper level menu twice to return to the Transport of SpeciesA menu.
6. Specify the concentration boundary conditions for fluid 1 (SUBDOMAIN1 and SUBDOMAIN2).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 273
Flow of Two Immiscible Fluids by Species Method
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
274 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. Repeat step (i) and step (ii) for BOUNDARY5 and BOUNDARY6.
Create a sub-task
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 275
Flow of Two Immiscible Fluids by Species Method
Material Data
Polydata indicates the material properties that are relevant for the sub-task by graying out the irrel-
evant properties. For this model, you will define only the viscosity of the material. The viscosity of
material 1 will be used if the species concentration is greater than 0.5, otherwise the viscosity of ma-
terial 2 will be used. This can be achieved by the use of PMAT.
d. Specify the value for , referred to as “fac” in the graphical user interface.
Modify fac
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
276 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
g. Click f1(...).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 277
Flow of Two Immiscible Fluids by Species Method
Polydata will ask to define pairs of values. A minimum of two pairs must be defined. Here you
will define (0.495, 5000) and (0.505, 10000), where the first index stands for species concentration
and the second index stands for viscosity [units: poise] value.
Enter 0.495 for X1( 1), and 5000 for f(X1)( 1).
n. Select SpeciesA.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
278 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
p. Click Upper level menu six times to return to the fluid 1 and 2 menu.
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY1 and click Modify.
Polydata prompts for the new value of the volumetric flow rate.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 279
Flow of Two Immiscible Fluids by Species Method
iv. Enter 3 [units: cm3/s] as the New value and click OK.
When the Automatic option is selected, Polydata automatically chooses the most appro-
priate method to compute the inflow condition.
b. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY2.
At a solid-liquid interface, the velocity of the liquid is that of the solid surface. Hence the fluid is
assumed to stick to the wall. This is known as the no-slip assumption because the liquid is assumed
to adhere to the wall, and so has no velocity relative to the wall.
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY3 and click Modify.
Polydata prompts for the new value of the volumetric flow rate.
iv. Click OK to accept the default value of 1 [units: cm3/s] for New Value.
d. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY4.
It is assumed that a fully developed velocity profile is reached at the exit, so the outflow condition
is appropriate. This condition essentially imposes a zero normal force ( ) that includes a pressure
term, and a zero tangential velocity ( ).
i. Select Zero wall velocity (vn=vs=0) along BOUNDARY5 and click Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
280 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
f. Retain the default condition Zero wall velocity (vn=vs=0) along BOUNDARY6.
The fluid is assumed to stick to the wall, since at a solid-liquid interface the velocity of the liquid
is that of the solid surface.
g. Click Upper level menu three times to return to the top-level Polydata menu.
Polydata asks you to confirm the current system units and fields that are to be saved to the results file for
postprocessing.
2. Click Accept.
3. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and the Polyflow format results file (res).
8.4.7. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens theView listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 281
Flow of Two Immiscible Fluids by Species Method
8.4.8. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
1. Double-click the Results tab in the Workbench analysis and read the results files saved by Polyflow.
CFD-Post reads the solution fields that were saved to the results file.
In the graphical window, right-click, and select View from +Z under Predefined Camera.
a. Click the Insert menu and select Contour or click the Contour button ( ).
b. Click OK to accept the default name (Contour 1) and display the details view below the
Outline tab.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
282 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. In the Location Selector dialog box that opens, select SUBDOMAIN1_surf, SUBDO-
MAIN2_surf, SUBDOMAIN3_surf, and SUBDOMAIN4_surf (use Ctrl for multiple selec-
tion) and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 283
Flow of Two Immiscible Fluids by Species Method
The velocity is much larger at the inlet of fluid 1 than at the inlet of fluid 2. There are two reasons for
this:
• The flow rate is three times larger for fluid 1 than for fluid 2.
• You are modeling an annular die. Hence the flow section is smaller for the interior channel than
for the exterior channel.
When the two fluids come into contact with each other, the interface between the two fluids is pushed
towards the exterior of the annular die.
• The die is annular, so even identical flow rates cause the interface to move in order to equilibrate
the flow sections.
• The viscosity of fluid 1 is higher than the viscosity of fluid 2. In the process of giving more room to
the most viscous fluid, its shearing decreases. This leads to a smaller global dissipation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
284 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. In Outline tree tab, under User Locations and Plots, deselect Contour 1.
b. Click the Insert menu and select Vector or click the button.
c. Click OK to accept the default name (Vector 1) and display the details view below the
Outline tab.
ii. In the Location Selector dialog box that opens, select the locations SUBDOMAIN1_surf,
SUBDOMAIN2_surf, SUBDOMAIN3_surf and SUBDOMAIN4_surf (use ctrl for multiple
selections) and click OK.
iii. In the Symbol tab, select Arrow 3D from the Symbol drop-down list.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 285
Flow of Two Immiscible Fluids by Species Method
You can see that the velocity is continuous across the interface. As both the fluids are Newtonian, the
velocity profile is a parabola on both sides of the interface. Since the force must be continuous across
the interface, the shear stress generated within fluid 1 is equal to the shear stress generated within
fluid 2 along the interface.
a. In Outline tree tab, under User Locations and Plots, deselect Vector 1.
b. Click the Insert menu and select Contour or click the Contour button ( ).
c. Click OK to accept the default name (Contour 2) and display the details view below the
Outline tab.
ii. In the Location Selector dialog box that opens, select SUBDOMAIN1_surf, SUBDO-
MAIN2_surf, SUBDOMAIN3_surf, and SUBDOMAIN4_surf (use Ctrl for multiple selec-
tion), and then click OK.
iv. In the Variable Selector dialog box that opens, select SpeciesA, and then click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
286 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. In Outline tree tab, under User Locations and Plots, deselect Contour 2.
c. Click OK to accept the default name (Isosurface 1) and display the details view below the
Outline tab.
i. In the Geometry tab, click the button next to Variable and select SpeciesA.
ii. Enter 0.5 for Value in order to locate the interface line.
iii. In the Color tab, select Constant from the Mode drop-down list and select pink by
iv. In the Render tab, select Draw As Lines from the Draw Mode drop-down list.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 287
Flow of Two Immiscible Fluids by Species Method
8.5. Summary
This tutorial introduced the concept of fluid layers flowing in the same duct. In Polydata, you learned
how to set up a species transport equation, a PMAT function, and you learned how to define a coextrusion
problem using species transport and PMAT functions. This method avoids the use of the remeshing
method, which is computationally expensive.
The species method, although less accurate, can help in quickly finding a solution when the die has a
complex shape. For more accurate results, the interface tracking method, as demonstrated in Flow of
Two Immiscible Fluids (p. 243), should be used. Generating a mesh for a complex die may be an issue
with the interface tracking method.
The location of the interface depends largely on the physical properties of the fluids involved, the
geometry of the channels, and the operating conditions (for example: flow rates of the fluids). A CFD
simulation with Polyflow allows you to test different setups (for example: in order to optimize the
feeding of a coextrusion die).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
288 of ANSYS, Inc. and its subsidiaries and affiliates.
Part 3: Blow Molding
The following blow molding tutorials are available:
1.1. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the setup procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 291
3D Thermoforming of a Blister
To reduce the computational run time, and utilizing the symmetric nature of the blister, only one quarter
of the blister/mold is modeled, Figure 1.1: Thermoforming of a Blister, Sheet (blue) and Mold (red) (p. 292).
From a geometric point of view, the initial (1/4) film has the following dimensions:
• Length = 15 mm
• Width = 5 mm
Dimensions are intentionally given in millimeters due to the small size of the object. The simulation will be
built around the system of units consisting of millimeters, grams and seconds.
The thickness compared to the length/width of the blister is rather small. This allows for the use of the
membrane (shell) element, which is suited for the analysis of 3D blow molding and thermoforming
simulations. The use of the membrane element is presently restricted to time-dependent flows and is
combined with Lagrangian representation (each mesh node is a material point). Node displacement
results from the time integration of nodal velocity.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
292 of ANSYS, Inc. and its subsidiaries and affiliates.
Problem Description
The finite element mesh and the boundary conditions are displayed in Figure 1.2: Finite Element Mesh,
Subdomains and Boundary Sets (p. 293). A 3D surface mesh has been generated for both the mold and
the film. The most important aspect is the proper description of the inner mold surfaces that will shape
the blister.
• Viscosity = 105
As seen in Figure 1.2: Finite Element Mesh, Subdomains and Boundary Sets (p. 293), the topology involves
two subdomains:
• Subdomain 1 = film
• Subdomain 2 = mold
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 293
3D Thermoforming of a Blister
The inflation pressure will be defined on the subdomain representing the film (Subdomain 1).
An important new concept is introduced in this tutorial: contact with a mold. Typically, two cases may
be encountered:
• The moving mold comes in contact with the shell and the shell acquires the mold velocity.
• The shell is inflated according to a certain rate and eventually comes into contact with the mold,
acquiring its shape.
1.3.1. Preparation
To prepare for running this tutorial:
3. Unzip the 3d_thermo_blister.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
294 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
2. Save the Ansys Workbench project using File → Save, entering blister as the name of the
project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Time-dependent problem(s)
• 2D shell geometry
Define molds
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 295
3D Thermoforming of a Blister
b. Click Upper level menu at the top of the Domain of the mold menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
296 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Contact conditions
Mold motion
A dialog box opens, asking you to specify the type of mold motion.
Enter 1 as the New value to impose a translation velocity, and click OK.
b. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 297
3D Thermoforming of a Blister
Click OK to accept the default value of 0 for the New value of velocity-x.
Polydata asks you to the specify the time dependence of the x-velocity. Click Upper
level menu, as there is no velocity in the x-direction.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
298 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Click OK to accept the default value of 0 for the New value of velocity-y, and click
Upper level menu, as there is no velocity in the y-direction.
Enter 10 [units: mm/s] for the New value of velocity-z and click OK.
v. Define the coordinate pairs (a,b) and (c,d) for the points that define the ramp function.
vi. In a similar manner, set the values for b, c, and d to 1.0, 0.103, and 0, respectively.
Figure 1.3: Ramp Function for Mold Velocity (p. 299) shows the ramp function you just defined.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 299
3D Thermoforming of a Blister
vii. Click the EVOL button at the top of the menu to disable evolution inputs.
viii. Click Upper level menu three times to return to the F.E.M. Task 1 menu.
Create a sub-task
b. Click Upper level menu button at the top of the Domain of the sub-task menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
300 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
b. Select Zero wall velocity (vn=vs=0) along Boundary 3 and click Modify.
iii. Click Upper level menu to continue specifying flow boundary conditions.
c. Select Zero wall velocity (vn=vs=0) along Boundary 4 and click Modify.
iii. Click Upper level menu to return to the Flow boundary conditions menu.
d. Click Inflation pressure imposed at the bottom of the Flow boundary conditions menu.
A dialog box opens, asking for the new value of the constant.
Enter 100000 [units: Pa] as the New value and click OK.
f. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
ii. Click Modify the value of a, and enter 0.1 as the New value.
iii. In a similar manner, set constants b, c, and d to 0, 0.11, and 1.0, respectively.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 301
3D Thermoforming of a Blister
Figure 1.4: Ramp Function for Pressure (p. 302) shows the ramp function you just defined.
Click the EVOL button at the top of the Polydata menu to disable evolution inputs.
Define contacts
The Modification of a contact problem menu will open with the Select a contact wall menu
item highlighted.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
302 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Magnify the view of the mold to ensure that you can see the darts that will be displayed.
Darts will be displayed in the Graphics Display window, as shown in the following figure.
iii. Click No in the dialog box that opens, to specify that the darts are not pointing towards
the mold body.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 303
3D Thermoforming of a Blister
If the direction of the darts is not clear to you, you can close the dialog box, rotate the view
and/or change the magnification, click Specify mold side / cavity side again, and then
answer the question appropriately.
A warning dialog box opens, saying that velocity prediction must be disabled, and that the
modification has automatically been done.
Click OK.
Define layers
In the dialog box that opens, enter blister as the New value.
The blister menu will open with the Material data menu item highlighted.
Material Data
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
304 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iv. Click Upper level menu twice to continue with material data specification.
v. Click Density.
Enter 0.001 [units: g/mm3] as the New value and click OK.
vii. Click Upper level menu to continue with the material data specification.
Select Inertia will be taken into account and click Upper level menu to continue
with material data specification.
Enter 5e-12 [units: g-mm/s/mm2] as the New value and click OK.
xi. Click Upper level menu twice to return to the blister menu.
Initial thickness
i. Click Constant.
Enter 0.35 [units: mm] as the New value and click OK.
ii. Click Upper level menu four times to return to the F.E.M. Task 1 menu.
Create a sub-task
Click No in the dialog box that opens asking if you want to copy the data of an existing sub-task.
a. Click Postprocessor.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 305
3D Thermoforming of a Blister
Enter Mass of product as the New value in the dialog box that asks for the title of the
sub-task and click OK.
Click OK in the dialog box that opens saying that the calculation will be performed on
Subdomain 1.
c. Confirm that the Density of layer [blister] is set to 1.00E-03, and click Upper level menu.
The plane is calculated by the following equation (which is displayed at the top of the Re-
striction of Layers by cutting planes menu).
(1.1)
Create a sub-task
Click No in the dialog box that opens asking if you want to copy the data of an existing sub-task.
a. Click Postprocessor.
Enter Permeability of product as the New value in the dialog box that asks for the
title of the sub-task and click OK.
Click OK in the dialog box that opens saying that the calculation will be performed on
Subdomain 1.
c. Confirm that the Permeability of layer [blister] is set to 5.00E-12, and click Upper level
menu.
See Equation 1.1 (p. 306) for more information on the coefficients.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
306 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Create a sub-task
Click No in the dialog box that opens asking if you want to copy the data of an existing sub-task.
a. Click Postprocessor.
Enter Volume of product as the New value in the dialog box that asks for the title of
the sub-task and click OK.
Click OK in the dialog box that opens saying that the calculation will be performed on
Subdomain 1.
See Equation 1.1 (p. 306) for more information on the coefficients.
f. Click Upper level menu to ignore contact with Subdomain 2 and return to F.E.M. Task 1.
Numerical parameters
Click OK to retain the default value of 0.0 [units: s] as the New value.
Click OK to retain the default value of 0.0001 [units: s] as the New Value.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 307
3D Thermoforming of a Blister
9. Click Upper level menu three times to return to the top-level Polydata menu.
1.3.7. Outputs
Outputs
Set to metric_mm/g/s/mA+Celsius
2. Click Upper level menu two times to return to the Outputs menu.
Output Triggering
3. Click Upper level menu twice to return to the top-level Polydata menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
308 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1. Click Accept.
2. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and the Polyflow format results file (res).
1.3.9. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 309
3D Thermoforming of a Blister
Here you can see the results of the postprocessor sub-tasks you created in Polydata. For additional
information on postprocessing sub-tasks, see Postprocessing Sub-Tasks (p. 305) and derived quantities
in the Polyflow User's Guide.
1.3.10. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
CFD-Post reads the solution fields that were saved to the results file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
310 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
b. In the dialog box that opens, click OK to accept the default name (Contour 1) and display
the details of Contour 1 below the Outline tree.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 311
3D Thermoforming of a Blister
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
312 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Double-click Default Transform in the Outline tree tab, under User Locations and Plots
(or right-click Default Transform and select Edit).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 313
3D Thermoforming of a Blister
b. Disable Instancing Info From Domain under the Definition tab in the details of Default
Transform.
d. Select Value from the Determine Angle From drop-down list in the Instance Definition
group box.
f. Enable Apply Reflection, and select ZX Plane from the Method drop-down list.
h. Click Apply and revise the magnification of the view to show the whole mold.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
314 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 315
3D Thermoforming of a Blister
c. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
316 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 317
3D Thermoforming of a Blister
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
318 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 319
3D Thermoforming of a Blister
c. To save the animation, expand the dialog box by clicking the button at the lower-right.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
320 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Click the file icon to the right and enter the path where you would like to save the
animation.
Enter Thickness.wmv for the name of the file and click Save to close the Save Movie
dialog box.
iv. Click the play button, to play the animation and save it as a file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 321
3D Thermoforming of a Blister
1.4. Summary
This tutorial introduced the concept of a blow molding problem. The mold moved into contact with
the film, where a constant pressure was applied to the film. This blew the film into the mold where it
assumed the shape of the mold.
You represented the film by a shell geometry under the valid assumption that the thickness of the film
was much smaller than the other two dimensions. Polyflow linearly interpolated the process vari-
ables—thickness, velocity and position. By reporting the individual time steps to CFD-Post you were
able to view the thickness of the product as a function of time.
F.E.M. Task 1
• Numerical parameters
– Adaptive meshing
• Modify size_min = .1
• Modify size_max = 1
• Modify dist_crit = .5
→ Modify Nstep = 4
For additional information on adaptive meshing, see adaptive meshing in the Polyflow User's Guide.
The results are shown in Figure 1.11: Effect of Adaption on Final Mesh and Thickness Variation (p. 323).
Note how the mesh changes as the Maxdiv value increases. The results do not change very much as
a result of the adaption, which indicates that the original solution was already mesh independent.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
322 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
1.6. Appendix
The appendix contains the following topics:
1.6.1. Contact Boundary Conditions
1.6.2. Remark on the Penalty Coefficient
1.6.3. Remeshing
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 323
3D Thermoforming of a Blister
Next to these operating attributes, some numerical parameters must be specified. A geometrical al-
gorithm is applied for detecting the occurrence of contact, while a penalty formation is used for the
treatment of contact. A penalty coefficient ensures that a geometrical contact is detected. It should
not be too small. A coefficient is also specified in the tangential direction. If the fluid sticks along the
wall, this tangential coefficient should preferably receive the same value as the penalty coefficient.
Two additional coefficients are also needed; a tolerance on penetration accuracy and an element
dilatation.
Presently, the penalty coefficient has been set to 1012, while the same value has been selected for
the coefficient along the tangential direction (slipping coefficient). The tolerance on penetration and
element dilatation equal 0.05 and 0.05 mm, respectively.
In classical thermoforming applications, such residual velocity will not produce any significant numer-
ical penetration of the fluid film/sheet through the mold in view of the short times involved (physically,
the thermoforming process is very fast). However, some situations may involve longer time scales
such as in the glass industry.
The question that is now raised concerns the best evaluation of the penalty coefficient. The penalty
formulation mainly establishes a balance between a force (for example: the inflation pressure, ) and
a penalty force because of contact. The penalty force is simply the product of the penalty coefficient,
, and the residual velocity of the film/sheet upon contact. The other elements of the momentum
equation can be ignored for the present consideration. Assuming a typical time scale (for example:
the simulation time), and a maximum penetration depth , a good penalty coefficient can be selected
as:
1.6.3. Remeshing
No remeshing must be specified. In the context of the membrane element, a Lagrangian representation
is applied where all mesh nodes are considered material points.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
324 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 2: 2D Axisymmetric Blow Molding
This tutorial contains the following sections:
2.1. Introduction
2.2. Prerequisites
2.3. Problem Description
2.4. Setup and Solution
2.5. Summary
2.6. Appendix
2.1. Introduction
Molding is a process of forcing a preform or a parison (preshaped sleeve) into a mold cavity so that the
preform assumes the shape of the cavity. There are numerous molding methods, including blow
molding, compression and transfer molding, and slush and rotational molding. These methods differ
in the formation of the preform and the filling of the mold cavity. Also, each processing method is
suitable for a specific class of polymers.
Blow molding is an important processing method for molding hollow articles such as bottles. The preform
is usually made by extrusion and forced between the mold halves by pressurization (blowing air). The
polymer solidifies upon contact with the cold mold and the finished product is then ejected. The homo-
geneity and rheological properties of the preform along with the operating conditions (temperature
and pressure variations) are crucial in this step and will affect the design of the processing machinery.
This process reflects all facets of polymer processing— the isothermal and transient flow of Newtonian
fluids in complex geometries with simultaneous structuring and solidification.
• Set material properties and boundary conditions for a 2D axisymmetric blow molding problem.
2.2. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the setup procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 325
2D Axisymmetric Blow Molding
A large pressure is applied to the preform, which enters the mold and eventually takes its shape. The
operating conditions must account for a low pressure drop at the entrance, low material waste, and
slow cooling to avoid premature solidification of the preform.
The cylindrical geometry of the preform (Figure 2.1: Problem Description (p. 326)) has an internal radius
of 2 cm and external radius of 3 cm (the initial thickness of the preform is 1 cm). The height of the
preform is 10 cm.
The domain for the problem is divided into two subdomains: one for the fluid preform (subdomain 1)
and the other for the mold (subdomain 2). Incompressibility and momentum equations are solved in
subdomain 1 (the fluid preform). The problem involves two free surfaces (boundary 2 and boundary 4,
shown in Figure 2). boundary 2 will eventually come into contact with the mold, and its position is
calculated as a part of the solution.
The fluid preform has a density of =1 g/cm3 and a viscosity of = 100000 poise. Inertia terms and
the effects of gravity will be included in the calculation.
The boundary sets for the problem are shown in Figure 2.2: Boundary Set for the Problem (p. 327), and
the conditions at the boundaries of the domain (for the preform) are:
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
326 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 327
2D Axisymmetric Blow Molding
2.4.1. Preparation
To prepare for running this tutorial:
3. Unzip the 2d_axi_blow_molding.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering Final-blow-mold as the name
of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Time-dependent problem(s)
• 2D axisymmetric geometry
The Current setup is updated to reflect the selected options. This example is a simulation of
blow molding for a 2D axisymmetric bottle and the mold is two-dimensional. The problem is
assumed to be time-dependent.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
328 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Define molds
a. Create a mold.
Adiabatic mold
c. When prompted, click OK to retain the default name for the mold (Mold 1).
SUBDOMAIN_1 is moved from the top list to the bottom list, indicating that the mold is
defined as SUBDOMAIN_2 .
e. Specify the boundary that represents the part of the mold that comes into contact with the
fluid.
Polyflow uses this information to determine the penetration distance (into the mold) of
every point of the free surface (BOUNDARY2).
Contact conditions
The free surface, BOUNDARY2 of the preform comes into contact with the mold wall,
BOUNDARY5 (as shown in Figure 2.2: Boundary Set for the Problem (p. 327)).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 329
2D Axisymmetric Blow Molding
f. Click Upper level menu four times to return to the F.E.M. Task 1 menu.
Create a sub-task
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
330 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The domain is divided into two subdomains, one for the fluid preform (SUBDOMAIN_1) and the other
for the mold (SUBDOMAIN_2). In this problem, the sub-task applies only to the preform.
SUBDOMAIN_2 is moved from the top list to the bottom list, indicating that the sub-task is
defined on SUBDOMAIN_1 .
Material data
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 331
2D Axisymmetric Blow Molding
c. Click Modify fac to specify the value of , which is referred to as “fac” in the graphical user
interface.
d. Enter 100000 [units: poise] as the New value and click OK.
e. Click Upper level menu two times to continue with the Material Data specification.
a. Click Density.
d. Click Upper level menu to continue with the Material Data specification.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
332 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In this problem inertia plays an important role. When internal pressure is applied, the preform expands,
and the fluid accelerates towards the mold. In order to obtain a realistic blowing time, inertia must
be taken into account.
c. Click Upper level menu to continue with the Material Data specification.
The fluid preform flows in the negative y direction under gravity, so specify the component of gravity
along the y direction ( ).
a. Click Gravity.
Modify gy
Polydata prompts for the new value of the gravity along the y-axis.
c. Enter -981 [units: cm/s2] as the New value and click OK.
d. Click Upper level menu two times to return to the blow molding menu.
No action is required to accept the default value. You can simply proceed to the next step. For 2D
axisymmetric models, Polydata recognizes the axis of symmetry from the mesh file and automatically
imposes the symmetry condition along the line .
The free surface boundary condition in contact detection problems is different from their simulations
in Polyflow. In blow molding problems, a free surface comes into contact with a solid mold. Polyflow
applies a contact detection algorithm at each location of the surface to detect the occurrence of the
contact.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 333
2D Axisymmetric Blow Molding
You need to specify the following for the free surface on BOUNDARY2:
• the contact wall (the boundary of the mold along which the contact is detected)
This determines the accuracy of the contact; the smaller its value, the deeper the contact is.
The fluid may slip along the contact wall, so to take this factor into account, a slipping coefficient
must be specified along the tangential direction.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY2 and click Modify.
iii. Specify where the free surface will contact the mold.
Polyflow uses the definition of the contact wall in the determination of the penetration
distance (into the mold) of every point of the free surface (BOUNDARY2).
As shown in Figure 2.2: Boundary Set for the Problem (p. 327), the free surface (BOUNDARY2)
of the preform comes into contact with the mold (BOUNDARY5).
With such a high value of the slipping coefficient, the fluid will stick to the contact wall.
vii. Click Upper level menu two times to return to the Kinematic condition menu.
d. Click Upper level menu to return to the Flow boundary conditions panel.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
334 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In contact detection problems, abrupt changes in the velocity field occur at the contact points
between the fluid preform and the mold. Polydata gives the warning message shown below.
Since the prediction of the velocity field in such cases destroys the prediction scheme, you can
continue by clicking OK.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY3 and click Modify.
c. Click Upper level menu to accept the default value of 0 for the normal velocity, .
d. Click Upper level menu to accept the default value of 0 for the tangential force, .
This boundary of the preform is subjected to pressure by the application of a normal force, so specify
a normal force along this boundary.
a. Select Zero wall velocity (vn=vs=0) along BOUNDARY4 and click Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 335
2D Axisymmetric Blow Molding
d. Click Upper level menu to return to the Flow boundary conditions menu.
2.4.6. Remeshing
This model involves free surfaces for which the positions are unknown. A portion of the mesh is affected
by the relocation of these boundaries. Hence a remeshing technique is applied on this part of the mesh.
The free surfaces are entirely contained within SUBDOMAIN_1, and hence only SUBDOMAIN_1 is affected
by the relocation of the free surfaces.
Global remeshing
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
336 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Click Accept the current setup in the Element distortion check menu.
In blow molding simulations, the finite-element mesh can undergo great deformations. The Ele-
ment distortion check menu deals with the detection of all possible distortions of the elements.
In this problem, you can accept the default options and proceed to the next step.
b. Click Upper level menu two times to return to the F.E.M. Task 1 menu.
Numerical parameters
1. Specify the parameters for the iterative scheme in the calculation of the free surface.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 337
2D Axisymmetric Blow Molding
This option specifies the time at which the solution procedure stops.
This option is used to define the initial time step, which is used for the calculation of the next
two time steps. After that, the step size is automatically calculated by Polyflow. This first time
step should be set according to the characteristic time scale of the process considered.
If a calculated value for the time step falls below the minimum for the time step at any point in
the calculation, the iterative scheme stops since this might be a symptom of calculation difficulties.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
338 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
In order to guarantee accuracy of the time-marching scheme and to avoid useless calculations
(rejection of inaccurate time steps), you can limit the growth of the time increment.
The tolerance is the admissible error between the predicted solution and the exact solution at a
particular time step. A very small value of the tolerance can result in large computational costs
and a very large value can result in wrong solution.
This option is used to select the maximum number of converged steps. If this value is reached,
the calculation stops, even if the upper time limit has not been reached.
2. Click Upper level menu three times to return to the top-level Polydata menu.
2.4.8. Outputs
You can specify how often Polyflow saves the solution data when it calculates a solution. In this tutorial,
save the results at every 4 time steps.
Outputs
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 339
2D Axisymmetric Blow Molding
4. Click Upper level menu twice to return to the top-level Polydata menu.
F.E.M. Task 1
Create a sub-task
b. Click Postprocessor.
c. Enter parison thickness as the New value for the title and click OK.
6. Specify the boundary sets representing the starting and ending borders to be used in the thickness
calculation.
Polyflow evaluates the distance between these borders at a point between them to determine the
thickness at that location.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
340 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
7. Click Upper level menu five times to return to the top-level Polydata menu.
Polydata asks you to confirm fields that are to be saved to the results file for postprocessing.
1. Click Accept.
2. Click Continue.
This accepts the default names for graphical output files (cfx.res) that are to be saved for postpro-
cessing, and the Polyflow format results file is (res).
2.4.11. Solution
In the following steps you will run Polyflow to calculate a solution for the model you just defined using
Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
2.4.12. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
CFD-Post reads the solution fields that were saved to the results file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 341
2D Axisymmetric Blow Molding
a. Click the Insert menu and select Contour or click the button.
b. In the box that opens, click OK to accept the default name (Contour 1) and display the
details below the Outline tree.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
342 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Select estim. THICKNESS from the Variable drop-down list, or click the ellipsis button
i. Click the Insert menu and select Text or click the button.
ii. Click OK to accept the default name (Text 1) and display the details view below the
Outline tab.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 343
2D Axisymmetric Blow Molding
a. Double-click Default Transform in the Outline tree tab, under User Locations and Plots
(or right-click Default Transform and select Edit).
b. Disable Instancing Info From Domain under the Definition tab in the details of Default
Transform.
c. Enable Apply Reflection, and select YZ Plane from the Method drop-down list.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
344 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
e. Click Apply.
Display the results at several time steps to see the shape and thickness of the parison during the blow
molding process.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 345
2D Axisymmetric Blow Molding
th
a. Select the 20 timestep and click Apply.
th
b. Select the 40 timestep and click Apply.
th
c. Select the 60 timestep and click Apply.
The thickness decreases as the parison inflates. At the final time step, the thickness is smallest
where the parison has been the most extended, (in the corner of the bottle). It is largest at the
top where the deformation was much less important due to the small diameter here.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
346 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 347
2D Axisymmetric Blow Molding
a. Click the Tools menu and select Animation or click the button.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
348 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
a. Double-click Contour 1 under the Outline tab to display the details view.
b. Select VELOCITIES from the Variable drop-down list and click Apply.
There is zero velocity at the contact between the parison and the mold, but the velocity magnitude
is still important where the fluid does not yet touch the mold. At the final time of the simulation,
the velocity is near zero, which indicates that the contact is completed. The residual value origin-
ates from the penalty formulation used for the contact, as explained in 3D Thermoforming of a
Blister (p. 291).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 349
2D Axisymmetric Blow Molding
2.5. Summary
This tutorial introduced a time-dependent problem with a 2D axisymmetric geometry for the mold.
Suitable assumptions were made regarding the nature of the preform and the operating conditions.
You analyzed the factors affecting the process in the postprocessing section. An optimization of the
preform shape could be performed in order to minimize the weight of the bottle while avoiding weak
(too thin) bottle walls.
You used a remeshing method that is most suited for contact detection problems. This problem also
introduced the concept of the calculation of free surfaces for contact detection problems. You used
efficient numerical techniques to more accurately solve a time-dependent problem.
2.6. Appendix
The appendix covers the following topics:
2.6.1. Remeshing Technique
2.6.2.Time Marching Scheme
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
350 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
Because this application involves contact occurring over time, a Lagrangian representation is used
for the free surface that undergoes the contact; this improves the robustness of the contact algorithm.
The Lagrangian on the border only technique remeshes based on the combination of a Lagrangian
representation on the border of the fluid domain and a minimum-pseudo-energy representation for
the inner mesh nodes. For additional information on this technique, see Lagrangian method on borders
in the Polyflow User's Guide.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 351
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
352 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 3: Plug-Assisted Thermoforming of a Blister
This tutorial is divided into the following sections:
3.1. Prerequisites
3.2. Problem Description
3.3. Setup and Solution
3.4. Summary
3.5. Appendix
3.1. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the setup procedure
will not be shown explicitly.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 353
Plug-Assisted Thermoforming of a Blister
Figure 3.1: Plug-Assisted Thermoforming of a Blister: Plug (Orange), Sheet (Blue), and Mold (Green)
in the Initial Configuration
To reduce the computational run time, and utilizing the symmetric nature of the blister, only one quarter
of the blister/plug/mold is modeled, Figure 3.1: Plug-Assisted Thermoforming of a Blister: Plug (Orange),
Sheet (Blue), and Mold (Green) in the Initial Configuration (p. 354). From a geometric point of view, the
initial (1/4) film has the following dimensions:
• length = 15 mm
• width = 5 mm
Dimensions are intentionally given in millimeters due to the small size of the object. The simulation will be
built around the system of units consisting of millimeters, grams and seconds.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
354 of ANSYS, Inc. and its subsidiaries and affiliates.
Problem Description
The thickness compared to the length/width of the blister is rather small. This allows the use of the
membrane (shell) element, which is suited for the analysis of 3D blow molding and thermoforming
simulations. The use of the membrane element is presently restricted to time-dependant flows and is
combined with Lagrangian representation. That is, each mesh node is a material point.
The finite element mesh and the boundary conditions are displayed in Figure 3.2: Finite Element Mesh,
Subdomains and Boundary Sets (p. 355). As shown, a full 3D finite element is built for the mold, the plug
and the film. Only a surface mesh is required for the three subdomains, but the most important aspect
remains the proper description of the inner mold surfaces which will shape the blister.
• viscosity = 105
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 355
Plug-Assisted Thermoforming of a Blister
As seen in Figure 3.2: Finite Element Mesh, Subdomains and Boundary Sets (p. 355), the mesh topology
involves three subdomains:
• Subdomain 1 = film
• Subdomain 2 = mold
• Subdomain 3 = plug
The inflation pressure will be defined on the subdomain representing the film (Subdomain 1).
An important new concept is introduced in this tutorial: plug-assisted contact with a mold. The film
acquires the mold velocity then the plug guides the film into the mold. Once inside the mold, the film
is inflated according to a certain rate where it eventually comes into contact with the mold, finally ac-
quiring its shape.
3.3.1. Preparation
To prepare for running this tutorial:
3. Unzip the plug_thermo_blister.zip file you have downloaded to your working folder.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
356 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering Blister-plug-assist as
the name of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Time-dependent problem(s)
• 2D shell geometry
Define molds
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 357
Plug-Assisted Thermoforming of a Blister
c. Click Upper level menu at the top of the Domain of the mold menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
358 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Contact conditions
Mold motion
A dialog box opens, asking you to specify the type of mold motion.
Enter 1 as the New value, to impose a translation velocity, and click OK.
b. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 359
Plug-Assisted Thermoforming of a Blister
Click OK to accept the default value of 0 [units: mm/s] for the New value of velocity-
x.
Polydata asks you to the specify the time dependence of the x-velocity. Click Upper
level menu, as there is no velocity in the x direction.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
360 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Click OK to accept the default value of 0 [units: mm/s] for the New value of velocity-
y, and click Upper level menu, as there is no velocity in the y direction.
Enter 10 [units: mm/s] for the New value of velocity-z and click OK.
v. Define the coordinate pairs (a,b) and (c,d) for the points that define the ramp function.
vi. In a similar manner, set the values for b, c, and d to 1.0, 0.103, and 0, respectively.
Figure 3.3: Ramp Function for Mold Velocity (p. 361) shows the ramp function you just defined.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 361
Plug-Assisted Thermoforming of a Blister
vii. Click the EVOL button at the top of the menu to disable evolution inputs.
viii. Click Upper level menu two times to return to the Define molds menu.
A dialog box opens, asking if you want to copy the data of an existing mold.
Click No.
c. Click Upper level menu at the top of the Domain of the mold menu.
Contact conditions
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
362 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Mold motion
A dialog box opens, asking you to specify the type of mold motion.
Enter 1 as the New value, to impose a translation velocity, and click OK.
b. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
Click OK to accept the default value of 0 [units: mm/s] for the New value of velocity-
x.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 363
Plug-Assisted Thermoforming of a Blister
Polydata asks you to the specify the time dependence of the x-velocity. Click Upper
level menu, as there is no velocity in the x direction.
Click OK to accept the default value of 0 [units: mm/s] for the New value of velocity-
y, and click Upper level menu, as there is no velocity in the y direction.
Enter -500 [units: mm/s] for the New value of velocity-z and click OK.
Enter 0.11 as the New value for time( 1) and click OK.
Retain the default, 0 as the New value for f(time)( 1) and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
364 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
vi. Click Insert new pair, and in a similar manner, define the following three pairs: (0.12,
1), (0.13, 1), and (0.14, 0).
Figure 3.4: Multi-Ramp Function for Plug Velocity (p. 365) shows the multi-ramp function
you just defined.
Click Upper level menu three times to return the Mold motion menu.
vii. Click the EVOL button at the top of the menu to disable evolution inputs.
viii. Click Upper level menu three times to return to the F.E.M. Task 1 menu.
Create a sub-task
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 365
Plug-Assisted Thermoforming of a Blister
c. Click Upper level menu button at the top of the Domain of the sub-task menu.
b. Select Zero wall velocity (vn=vs=0) along Boundary 3 and click Modify.
iii. Click Upper level menu to continue specifying flow boundary conditions.
c. Select Zero wall velocity (vn=vs=0) along Boundary 4 and click Modify.
iii. Click Upper level menu to return to the Flow boundary conditions menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
366 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
d. Click Inflation pressure imposed at the bottom of the Flow boundary conditions menu.
A dialog box opens, asking for the new value of the constant.
Enter 1e05 [units: Pa] as the New value and click OK.
f. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
ii. Click Modify the value of a, and enter 0.14 as the New value.
iii. In a similar manner, set constants b, c, and d to 0, 0.15, and 1.0 respectively.
Figure 3.5: Ramp Function for Pressure (p. 368) shows the ramp function you just defined.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 367
Plug-Assisted Thermoforming of a Blister
Click the EVOL button at the top of the Polydata menu to disable evolution inputs.
Define contacts
The Modification of a contact problem menu will open with the Select a contact wall menu
item highlighted.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
368 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Increase the size of the darts that will be used to display the orientation, to ensure that
they are visible.
ii. Rotate the view to an oblique angle and zoom in on the mold.
Darts will be displayed in the Graphics Display window, as shown in the following figure.
iv. Click No in the dialog box that opens, to specify that the darts are not pointing towards
the mold body.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 369
Plug-Assisted Thermoforming of a Blister
If the direction of the darts is not clear to you, you can close the dialog box, rotate the view
and/or change the magnification, click Specify mold side / cavity side again, and then
answer the question appropriately.
The Modification of a contact problem menu will open with the Select a contact wall menu
item highlighted.
Click Yes in the dialog box that asks if you want to activate contact release.
i. Rotate the view and change the magnification so that you can see the plug.
Darts will be displayed in the Graphics Display window, as shown in the following figure.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
370 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. Click Yes in the dialog box that opens, to verify that the darts are pointing toward the
plug body (away from contact with the film).
If the direction of the darts is not clear to you, you can close the dialog box, rotate the view
and/or change the magnification, click Specify mold side / cavity side again, and then
answer the question appropriately.
iv. Click Upper level menu to return to the Define contacts menu.
e. Click OK in the warning box that opens, saying that velocity prediction must be disabled.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 371
Plug-Assisted Thermoforming of a Blister
Define layers
Enter blister as the New value in the dialog box that opens and click OK.
The blister menu will open with the Material data menu item highlighted.
Material Data
iv. Click Upper level menu twice to continue with material data specification.
v. Click Density.
Enter 0.001 [units: g/mm3] as the New value and click OK.
vii. Click Upper level menu to continue with the material data specification.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
372 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Initial thickness
i. Click Constant.
Enter 0.35 [units: mm] as the New value and click OK.
ii. Click Upper level menu four times to return to the F.E.M. Task 1 menu.
Numerical parameters
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 373
Plug-Assisted Thermoforming of a Blister
The maximum number of steps must be increased due to the case containing two contacts and “sharp”
changes in the kinematics. This higher number was determined by running the problem with a max-
imum of 200 steps and observing that more steps were required to reach the final time.
It is recommended that you specify a convergence criterion of 0.0001 or lower when contact release
is activated.
12. Click Upper level menu three times to return to the top-level Polydata menu.
3.3.7. Outputs
Outputs
Set to metric_mm/g/s/mA+Celsius
4. Click Upper level menu three times to return to the top-level Polydata menu.
1. Click Accept.
2. Click Continue.
This accepts the default names for the graphical output files (cfx.res) that are to be saved for
postprocessing, and the Polyflow format results file (res).
3.3.9. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
374 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
3.3.10. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
CFD-Post reads the solution fields that were saved to the results file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 375
Plug-Assisted Thermoforming of a Blister
a. Click the Insert menu and select Contour or click the button.
b. In the box that opens, click OK to accept the default name (Contour 1) and display the
details view below the Outline tree.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
376 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Select THICKNESS from the Variable drop-down list, or click the ellipsis button ( )
on the right and select THICKNESS.
v. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 377
Plug-Assisted Thermoforming of a Blister
a. Double-click Default Transform in the Outline tree tab, under User Locations and Plots
(or right-click Default Transform and select Edit).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
378 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
The details view of Default Transform will be displayed below the Outline tab.
b. Perform the following steps in the Definition tab of the details view.
iii. Select Value from the Determine Angle From drop-down list in the Instance Definition
group box.
v. Enable Apply Reflection, and select ZX Plane from the Method drop-down list.
c. Rotate the view and change the magnification, as shown in Figure 3.7: Contours of Thickness
on the Whole Thermoformed Blister (p. 380).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 379
Plug-Assisted Thermoforming of a Blister
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
380 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
c. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 381
Plug-Assisted Thermoforming of a Blister
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
382 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 383
Plug-Assisted Thermoforming of a Blister
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
384 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
c. To save the animation, expand the dialog box by clicking the button at the lower-right.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 385
Plug-Assisted Thermoforming of a Blister
ii. Click the file icon to the right and enter the path where you would like to save the
animation.
Enter Thickness.wmv for the name of the file and click Save to close the Save Movie
dialog box.
iv. Click the play button, to play the animation and save it as a file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
386 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
3.4. Summary
This tutorial introduced the concept of a plug-assisted blow molding problem. The mold moved into
contact with the film, where a plug guided the film into the mold, and a constant pressure was applied
to the film. This blew the film into the mold where it assumed the shape of the mold.
You represented the film by a shell geometry under the valid assumption that the thickness of the film
was much smaller than the other two dimensions. Polyflow linearly interpolated the process vari-
ables—thickness, velocity and position. By reporting the individual time steps to CFD-Post you were
able to view the thickness of the product as a function of time.
3.5. Appendix
The appendix covers the following topics:
3.5.1. Contact Boundary Conditions
3.5.2. Remark on the Penalty Coefficient
3.5.3. Remeshing
Next to these operating attributes, some numerical parameters have to be specified. A geometrical
algorithm is applied for detecting the occurrence of contact, while a penalty formation is used for
the treatment of contact. A penalty coefficient makes sure that a geometrical contact is detected. It
should not be too small. A coefficient is also to be specified in the tangential direction. If the fluid
sticks along the wall, this tangential coefficient should preferably receive the same value as the penalty
coefficient. Two additional coefficients are also needed; a tolerance on penetration accuracy and an
element dilatation.
Presently, the penalty coefficient has been set to 1012, while the same value has been selected for
the coefficient along the tangential direction (slipping coefficient). The tolerance on penetration and
element dilatation equal 0.05 and 0.05 mm respectively.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 387
Plug-Assisted Thermoforming of a Blister
In classical thermoforming applications, such residual velocity will not produce any significant numer-
ical penetration of the fluid film/sheet through the mold in view of the short times involved (physically,
the thermoforming process is very fast). However, some situations may involve longer time scales
such as in the glass industry.
The question that is now raised concerns the best evaluation of the penalty coefficient. The penalty
formulation mainly establishes a balance between a force (for example, the inflation pressure, ) and
a penalty force because of contact. The penalty force is simply the product of the penalty coefficient,
, and the residual velocity of the film/sheet upon contact. The other elements of the momentum
equation can be ignored for the present consideration. Assuming a typical time scale (for example,
the simulation time), and a maximum penetration depth , a good penalty coefficient can be selected
as:
3.5.3. Remeshing
The results of this tutorial could be refined and improved with the use of adaptive meshing.
In the context of the membrane element, a Lagrangian representation is applied where all mesh
nodes are considered material points. Therefore, the only available remeshing technique is Lagrangian
and is the one that should be specified for this case.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
388 of ANSYS, Inc. and its subsidiaries and affiliates.
Chapter 4: 3D Blow Molding of a Bottle
This tutorial is divided into the following sections:
4.1. Prerequisites
4.2. Description
4.3. Setup and Solution
4.4. Summary
4.5. Further Improvements
4.6. Appendix
4.1. Prerequisites
This tutorial assumes that you are familiar with the menu structure in Polydata and Workbench and
that you have solved or read 2.5D Axisymmetric Extrusion (p. 69). Some steps in the setup procedure
will not be shown explicitly.
4.2. Description
This tutorial simulates a typical blow molding situation for a bottle. In the present case, it is assumed
that a cylindrical parison with uniform thickness distribution has been extruded. The present calculation
involves two major steps; parison pinch-off due to mold closing, and inflation. Figure 4.1: Blow Molding
Initial Configuration (p. 390) shows a sketch of the process in the initial configuration, before the pinch-
off and parison inflation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 389
3D Blow Molding of a Bottle
From a geometric point of view, the initial parison has the following dimensions:
• height = 0.276 m
• radius = 0.0225 m
The thickness of the fluid parison is much smaller than the other two dimensions of the bottle, which
allows for the use of the membrane (shell) element, which is suited for the analysis of 3D blow molding
simulations. It is important to remember when preparing the surface mesh, that the mesh elements on
the mold should not be the same order of magnitude as the expected final local thickness. The use of
the membrane element is presently restricted to time-dependant flows and is combined with Lagrangian
representation. That is, each mesh node is a material point.
The finite element mesh and the boundary conditions are displayed in Figure 4.2: Finite Element Mesh,
Subdomains, and Boundary Sets (p. 391). As shown, a full 3D finite element is built for both the mold
and the parison. Only a surface mesh is needed for both the mold and the parison, but the most im-
portant aspect remains the proper description of the inner mold surfaces that will shape the bottle.
• viscosity = 104
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
390 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
As seen in Figure 4.2: Finite Element Mesh, Subdomains, and Boundary Sets (p. 391), the mesh topology
involves three subdomains (MoldLeft, parison, and Moldright) and two boundary sets (TopEdge and
BottomEdge). The fluid parison is covered by the subdomain named parison while MoldLeft and
MoldRight will be defined as molds. Along boundary sets TopEdge and BottomEdge, a symmetry
boundary condition will be imposed. The inflation pressure will be defined on the subdomain repres-
enting the parison.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 391
3D Blow Molding of a Bottle
4.3.10. Postprocessing
4.3.1. Preparation
To prepare for running this tutorial:
3. Unzip the 3d_blow_molding_bottle.zip file you have downloaded to your working folder.
4. Start Workbench from Start > Ansys 2024 R1 > Workbench 2024 R1.
2. Save the Ansys Workbench project using File → Save, entering instanet-PF-only as the
name of the project.
When Polydata starts, the Create a new task menu item is highlighted, and the geometry for the
problem is displayed in the Graphics Display window.
• F.E.M. task
• Time-dependent problem(s)
• 2D shell geometry
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
392 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Define molds
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 393
3D Blow Molding of a Bottle
c. Click Upper level menu at the top of the Domain of the mold menu.
Contact conditions
Mold motion
A small dialog box opens, asking you to specify the type of mold motion.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
394 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Enter 1 as the New value, to impose a translation velocity, and click OK.
b. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 395
3D Blow Molding of a Bottle
Click OK to accept the default value of 0 for the New value of velocity-x.
Polydata asks you to the specify the time dependence of the x-velocity. Click Upper
level menu, as there is no velocity in the x direction.
Click OK to accept the default value of 0 for the New value of velocity-y, and click
Upper level menu, as there is no velocity in the y direction.
Enter 0.736842 [units: m/s] for the New value of velocity-z and click OK.
v. Define the coordinate pairs (a,b) and (c,d) for the points that define the ramp function.
vi. In a similar manner, set the values for b, c, and d to -1.0, 0.1, and 0, respectively.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
396 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Figure 4.3: Ramp Function for Right Mold Velocity (p. 397) shows the ramp function you just
defined.
vii. Click the EVOL button at the top of the menu to disable evolution inputs.
viii. Click Upper level menu two times to return to the Define molds menu.
A dialog box opens, asking if you want to copy the data of an existing mold.
Click No.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 397
3D Blow Molding of a Bottle
c. Click Upper level menu at the top of the Domain of the mold menu.
Contact conditions
Mold motion
A small dialog box opens, asking you to specify the type of mold motion.
Enter 1 as the New value, to impose a translation velocity, and click OK.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
398 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
b. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
Click OK to accept the default value of 0 for the New value of velocity-x.
Polydata asks you to the specify the time dependence of the x-velocity. Click Upper
level menu, as there is no velocity in the x direction.
Click OK to accept the default value of 0 for the New value of velocity-y, and click
Upper level menu, as there is no velocity in the y direction.
Enter 0.736842 [units: m/s] for the New value of velocity-z and click OK.
v. Define the coordinate pairs (a,b) and (c,d) for the points that define the ramp function.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 399
3D Blow Molding of a Bottle
vi. In a similar manner, set the values for b, c, and d to 1.0, 0.1, and 0, respectively.
Figure 4.4: Ramp Function for Left Mold Velocity (p. 400) shows the ramp function you just
defined.
vii. Click Upper level menu to return the Mold motion menu.
d. Click the EVOL button at the top of the menu to disable evolution inputs.
e. Click Upper level menu three times to return to the F.E.M. Task 1 menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
400 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Create a sub-task
A small dialog box opens, asking for the title of the problem.
c. Click Upper level menu button at the top of the Domain of the sub-task menu.
a. Select Zero wall velocity (vn=vs=0) along BOTTOMEDGE and click Modify.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 401
3D Blow Molding of a Bottle
iii. Click Upper level menu to continue specifying flow boundary conditions.
b. Select Zero wall velocity (vn=vs=0) along TOPEDGE and click Modify.
iii. Click Upper level menu to return to the Flow boundary conditions menu.
c. Click Inflation pressure imposed at the bottom of the Flow boundary conditions menu.
Darts will be displayed in the Graphics Display window, to indicate the orientation of the pressure
on the parison.
As shown in the figure that follows, the darts point into the center of the parison.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
402 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
A dialog box opens, asking for the new value of the constant.
Enter -1e4 [units: Pa] as the New value and click OK.
f. Click the EVOL button at the top of the Polydata menu to enable evolution inputs.
h. Define the settings in the Time dependence of inflation pressure menu that opens.
ii. Click Modify the value of a, enter 0.1 as the New value, and click OK.
iii. In a similar manner, set the values for b, c, and d to 0, 0.2, and 1.0, respectively.
Figure 4.5: Ramp Function for Pressure (p. 404) shows the ramp function you just defined.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 403
3D Blow Molding of a Bottle
j. Click the EVOL button at the top of the Polydata menu to disable evolution inputs.
Define contacts
The Modification of a contact problem menu will open with the Select a contact wall menu
item highlighted.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
404 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
i. Increase the size of the darts that will be used to display the orientation, to ensure that
they are visible.
Darts will be displayed in the Graphics Display window, as shown in the following figure.
iv. Click No in the dialog box that opens, to specify that the darts are not pointing towards
the mold body.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 405
3D Blow Molding of a Bottle
If the direction of the darts is not clear to you, you can close the dialog box, rotate the view
and/or change the magnification, click Specify mold side / cavity side again, and then
answer the question appropriately.
The Modification of a contact problem menu will open with the Select a contact wall menu
item highlighted.
i. Rotate the view slightly to display the bottle-shaped cavity from an oblique angle.
Darts will be displayed in the Graphics Display window, as shown in the following figure.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
406 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
iii. Click No in the dialog box that opens, to specify that the darts are not pointing towards
the mold body.
If the direction of the darts is not clear to you, you can close the dialog box, rotate the view
and/or change the magnification, click Specify mold side / cavity side again, and then
answer the question appropriately.
iv. Click Upper level menu to return to the Define contacts menu.
e. Click OK in the warning box that opens, to acknowledge that the velocity prediction must
be disabled.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 407
3D Blow Molding of a Bottle
Define layers
In the dialog box that opens, enter parison as the New value.
The parison menu will open with the Material data menu item highlighted.
Material Data
iv. Click Upper level menu twice to continue with material data specification.
v. Click Density.
Enter 900 [units: kg/m3] as the New value and click OK.
vii. Click Upper level menu to continue with the material data specification.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
408 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ix. Select Inertia will be taken into account and click Upper level menu twice to return
to the parison menu.
Initial thickness
i. Click Constant.
ii. Click Upper level menu four times to return to the F.E.M. Task 1 menu.
Numerical parameters
9. Click Upper level menu three times to return to the top-level Polydata menu.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 409
3D Blow Molding of a Bottle
4.3.7. Outputs
Outputs
The Current system is shown at the top of the Change system of Units for specific outputs
menu.
3. Click Upper level menu twice to return to the top-level Polydata menu.
1. Click Accept.
2. Click Continue.
This accepts the default names for the graphical output files (cfx.res) that are to be saved for
postprocessing, and the Polyflow format results file (res).
4.3.9. Solution
Run Polyflow to calculate a solution for the model you just defined using Polydata.
1. Run Polyflow by right-clicking the Solution cell of the simulation and selecting Update.
This executes Polyflow using the data file as standard input, and writes information about the problem
description, calculations, and convergence to a listing file (polyflow.lst).
Workbench opens the View listing file dialog box, which displays the listing file.
b. It is a common practice to confirm that the solution proceeded as expected by looking for
the following printed at the bottom of the listing file:
The computation succeeded.
4.3.10. Postprocessing
Use CFD-Post to view the results of the Polyflow simulation.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
410 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
CFD-Post reads the solution fields that were saved to the results file.
2. Change the view in the Graphics Display window as shown in the figure that follows.
a. Click the Insert menu and select Contour or click the button.
b. In the dialog box that opens, click OK to accept the default name (Contour 1) and display
the details view below the Outline tree.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 411
3D Blow Molding of a Bottle
v. Click Apply.
d. Disable the Wireframe in the Outline tree tab, under User Locations and Plots.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
412 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 413
3D Blow Molding of a Bottle
b. Enter 0.7 for Transparency in the Render tab in the details view of MOLDLEFT_surf.
The contours of thickness on the parison would not be visible without increasing the transparency
of the mold.
c. Click Apply.
c. Click Apply.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
414 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 415
3D Blow Molding of a Bottle
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
416 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 417
3D Blow Molding of a Bottle
c. To save the animation, expand the dialog box by clicking the button at the lower-right.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
418 of ANSYS, Inc. and its subsidiaries and affiliates.
Setup and Solution
ii. Click the file icon to the right and enter the path where you would like to save the
animation.
Enter Thickness.wmv for the name of the file and click Save to close the Save Movie
dialog box.
iv. Click the play button, to play the animation and save it as a file.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 419
3D Blow Molding of a Bottle
4.4. Summary
This tutorial introduced the concept of a parison blow molding problem. The two halves of the mold
moved into contact with the parison, where it became pinched, and a vacuum was applied to the par-
ison. This blew the parison into the mold where it assumed the shape of the mold, which was a bottle
in this case.
You represented the parison by a shell geometry under the valid assumption that the thickness of the
parison was much smaller than the other two dimensions (diameter and height). Polyflow linearly inter-
polated the process variables—thickness, velocity and position. By reporting the individual time steps
to CFD-Post you were able to view the thickness of the product as a function of time.
F.E.M. Task 1
• Numerical parameters
– Adaptive meshing
→ Modify Nstep = 5
→ Modify Maxdiv = 1
For additional information on adaptive meshing, see adaptive meshing in the Polyflow User's Guide.
The results are shown in Figure 4.11: Effect of Adaption on Final Mesh and Thickness Variation (p. 421).
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
420 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
4.6. Appendix
The appendix contains the following topics:
4.6.1. Contact Boundary Conditions
4.6.2. Remark on the Penalty Coefficient
4.6.3. Remeshing
4.6.4. Evolutions
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 421
3D Blow Molding of a Bottle
Along with these operating attributes, some numerical parameters must be specified. A geometrical
algorithm is applied for detecting the contact, while a penalty formation is used for the treatment of
the contact. A penalty coefficient ensures that a geometrical contact is detected. It should not be too
small. A coefficient is also specified in the tangential direction. If the fluid sticks along the wall, this
tangential coefficient should preferably receive the same value as the penalty coefficient. Two addi-
tional coefficients are also needed; a tolerance on penetration accuracy and an element dilatation.
Presently, the penalty coefficient and the slipping coefficient (tangential direction) are both set to
1012. The tolerance on penetration and element dilatation are equal to 0.001 and 0.002 m respectively.
In classical blow molding applications, such residual velocity will not produce any significant numer-
ical penetration of the fluid parison through the mold in view of the short times involved (physically,
blow molding process is very fast). However, some situations may involve longer time scales such as
in the glass industry.
The question that is now raised concerns the best evaluation of the penalty coefficient. The penalty
formulation mainly establishes a balance between a force (for example: the inflation pressure, ) and
a penalty force because of contact. The penalty force is simply the product of the penalty coefficient,
, and the residual velocity of the parison upon contact. The other variables of the momentum
equation can be omitted for the present problem. Considering a typical time scale, (for example,
the simulation time), and a maximum penetration depth that can be practically accepted, , a good
penalty coefficient can be selected as:
4.6.3. Remeshing
In the context of the membrane element, a Lagrangian representation is applied where all mesh
nodes are considered material points. Therefore, the only available remeshing technique is Lagrangian
and is the one that should be specified for this case.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
422 of ANSYS, Inc. and its subsidiaries and affiliates.
Appendix
4.6.4. Evolutions
The present case involves a mold motion followed by inflation.
For the mold motion, the x and y-components are zeros. The two mold halves move only in the z
direction at the same speed but in opposite directions. The two mold halves move at 0.736842 m/s
in the z direction. To control the duration and the direction of the motion, a simple ramp function is
applied on the mold speed. The ramp function is multiplied by the z-velocity component to give each
half of the mold the proper speed in the appropriate direction.
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
of ANSYS, Inc. and its subsidiaries and affiliates. 423
Release 2024 R1 - © ANSYS, Inc. All rights reserved. - Contains proprietary and confidential information
424 of ANSYS, Inc. and its subsidiaries and affiliates.