Chapter 3: Fluent Material Processing Tutorials

This chapter include several Material Processing Tutorials for Ansys Fluent.

3.1. 3D Polymer Extrusion

3.1.1. Introduction

This tutorial illustrates the simulation of a 3D extrusion process using the Fluent Materials Processing workspace (Fluent Materials Processing Workspace). The workspace will allow you to set up and solve a polymer extrusion problem so you can easily obtain an accurate prediction of the extrudate shape for a given die geometry under prescribed operating conditions.

In this tutorial you will learn how to use the Fluent Materials Processing workspace to:

  • Easily set up your extrusion simulation

  • Calculate a solution

  • Analyze the results.

3.1.2. Problem Description

This problem deals with the flow of a Newtonian fluid through a three-dimensional die. Due to the symmetry of the problem (the cross-section of the die is a square), the computational domain is defined for a quarter of the geometry and two planes of symmetry are defined.

The melt enters the die as shown in Figure 3.1: Problem Description at a flow rate of Q = 10 cm3/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 any further. It is assumed that subdomain 2 is long enough to account for all the deformation of the extrudate.

Figure 3.1: Problem Description

Problem Description

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 3.1: Problem Description) 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 3.2: Boundary Sets for the Problem, and the conditions at the boundaries of the domains are as follows.

  • boundary 1: flow inlet, volumetric flow rate Q = 10 cm3/s

  • boundary 2: zero velocity

  • boundary 3: symmetry plane

  • boundary 4: symmetry plane

  • boundary 5: free surface

  • boundary 6: flow exit

Figure 3.2: Boundary Sets for the Problem

Boundary Sets for the Problem

3.1.3. Setup and Solution

3.1.3.1. Preparation

To prepare for running this tutorial:

  1. Prepare a working folder for your simulation.

  2. Download the 3d_extrusion.zip file here .

  3. Unzip the 3d_extrusion.zip file you have downloaded to your working folder.

  4. The mesh file ext3d.msh can be found in the unzipped folder.

3.1.3.2. Launching Ansys Fluent

  1. Use the Fluent Launcher to start Ansys Fluent.

  2. Select Materials Processing in the list of Fluent workspaces.

  3. Click Start.

3.1.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

The workspace is a version of Ansys Fluent that utilizes the power of the Ansys Polyflow solver to simulate polymer flows such as extrusion, blow molding, pressing, and so forth.

  1. Open an extrusion geometry by reading in a Polyflow mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the mesh file from your working folder (ext3d.msh)

    3. Set the Mesh Length Unit to cm, and close the dialog box.

  2. Set up an extrusion simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Template...

    This will display the Simulation Template dialog.

    1. Set the Type to Extrusion.

    2. Set the Goal to be Predict extrudate shape.

    3. Keep the Number of fluids as 1.

    4. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.1.3.4. General Properties

Click the General node in the Outline View to display general settings for the problem setup.

  1. Select General in the Outline View to review general properties of the simulation.

  2. Review the settings in the properties page.

  3. Make sure that Steady is selected for the Calculation Type.

  4. Enter 3D die swell for the Task Name.

3.1.3.5. Material Properties

Fluent indicates which material properties are relevant by graying out the irrelevant properties. For this model you will only define the viscosity of the material.

Even though the default material's (fluid) properties are up-to-date, you will change how that material's viscosity is defined.

  1. Select fluid under Materials in the Outline View to review material properties of the simulation.

  2. In the Properties - fluid panel, click the icon to the left of the Viscosity Law category to define the viscosity properties of the fluid.

  3. For Shear Rate Dependence, select Cross.

    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:

    (3–1)

    where:

    = zero shear-rate viscosity = 85000 poise

    = natural time = 0.2 s

    = Cross law index = 0.3

    = shear rate

  4. For Zero Shear Viscosity [Pa s], enter 85000.

  5. For Time Constant [s], enter 0.2.

  6. For Cross Law Index, enter 0.3.

3.1.3.6. Cell Zone Properties

In this portion of the tutorial, you will assign your fluid cell zones.

  1. Select fluid-zone under Cell Zones in the Outline View to review cell zone properties of the simulation.

  2. For Zones, select the field to open the ZoneId dialog.

  3. In the ZoneId dialog, select both subdomain1 and subdomain2 and click OK to keep your selections and close the dialog.

3.1.3.7. Boundary Condition Properties

In the following steps you will set the conditions at each of the boundaries of the simulation.

  1. Set inlet conditions.

    1. Select inlet under Fluid Boundary Conditions in the Outline View to review inlet boundary conditions for the simulation.

    2. Select boundary1 as the Boundary Zone.

    3. Select Volume Flow as the Flow Specification.

    4. Select fluid as the Incoming material.

    5. Enter 0.00001 m^3/s (10 cm^3/s) as the Volume Flow Rate.

  2. Set outlet conditions.

    1. Select extrudate-exit under Fluid Boundary Conditions in the Outline View to review exit boundary conditions for the simulation.

    2. Select boundary6 as the Boundary Zone.

  3. Set symmetry conditions.

    1. Select Fluid Boundary Conditions in the Outline View, and click the New... button in its properties window, or right-click in the tree and add a new condition to the tree.

    2. Rename the new condition to Symmetry.

    3. Set the Type to Symmetry.

    4. Select boundary3-subdomain1, boundary3-subdomain2, boundary4-subdomain1, and boundary4-subdomain2 as the Boundary Zone.

  4. Set free surface conditions.

    1. Select free-surface under Fluid Boundary Conditions in the Outline View to review free surface boundary conditions for the simulation.

    2. Select boundary5 as the Boundary Zone.

    3. Select boundary2 as the Fixed Part.

  5. Set wall conditions.

    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 condition because the liquid is assumed to adhere to the wall, and hence, has no velocity relative to the wall.

    1. Select wall under Fluid Boundary Conditions in the Outline View to review wall boundary conditions for the simulation.

    2. Select boundary2 as the Boundary Zone.

3.1.3.8. Mesh Deformation Properties

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.

  1. Under Mesh Deformations, select extrudate to edit the default properties of the extrudate exit for your simulation.

  2. Select subdomain2 for Zones.

  3. Select Optimesh-3D for the Extrudate Deformation Method.

  4. Select subdomain2-subdomain1 for the Inlet Section.

  5. Select boundary6 for the Outlet Section.

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).

3.1.3.9. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and other solution properties. Most items indicate that current default values are appropriate.

  1. Review your problem setup.

    1. Select Run Calculation to see available options.

    2. Click Check to review your simulation settings. Fluent will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  2. Calculate a solution.

    1. Click Calculate to begin the computing a solution.

    2. Once the calculations are complete, check the solution's listing in the Transcript tab, at the bottom of the Graphics window, to confirm its success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  3. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering extrusion-3D as the name of the session.

3.1.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can setup contours, vectors, and so on.

  1. Display the velocity distribution on the boundaries.

     ResultsGraphics Contour New

    Figure 3.3: Contour Properties

    Contour Properties

    1. Keep the default Name (contour-1).

    2. For the Field, select VELOCITIES[].

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Once your selections are complete, click OK to close the dialog.

      Figure 3.4: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. Keep the remaining defaults and click Display to see the velocity magnitude contour plot in the Graphics window.

      Use the toolbars in the Graphics window to adjust the display so that it is in an isometric orientation, and fits the window.

      Figure 3.5: Contours of Velocity Magnitude

      Contours of Velocity Magnitude


  2. Display contours of velocity in cross-sectional planes, at Z = 0, 0.08, 0.15, and 0.45 m.

    1. Create a cross-section planes at Z = 0 m.

       Results Surfaces Plane New

      In the Properties - plane-1 panel, under Plane Settings, select XY Plane for the Creation Mode, enter 0 for Z, and click Display.

    2. Create additional cross-section planes for Z = 0.08, 0.15, and 0.45 m.

      Once complete, you should have four separate plane surfaces available.

      You could also have used the Create Multiple Plane dialog to easily create multiple surfaces.

    3. Display contours of velocity in the cross-sectional planes.

       ResultsGraphics Contour New

      In the properties of contour-2, select VELOCITIES[] for the Field, and for the Surfaces, select the four plane surfaces you just created, keep the remaining defaults, and click Display..

      Figure 3.6: Velocity Profiles at Cross-Sections

      Velocity Profiles at Cross-Sections

3.1.5. Summary

This tutorial introduced the concept of a 3D extrusion problem using the Fluent Materials Processing workspace. 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 addition you learned how to use the optimesh remeshing method, which is recommended for 3D extrusion problems.

3.2. Multiple Material Coextruded Tubing

3.2.1. Introduction

This tutorial demonstrates the simulation of a multiple material coextrusion process using the Fluent Materials Processing workspace. A coextrusion process allows for two or more different polymer materials to be combined into a single product. Depending on the materials used, the process allows for the creation of products with unique properties, such as increased strength, resistance to corrosion, or insulation. The versatility of coextruded plastic profiles and tubing makes coextrusion an ideal manufacturing process for a wide range of products.

In this tutorial you will learn how to use the Fluent Materials Processing workspace to set up a coextrusion simulation using multiple materials, calculate the solution, and analyze the results.

3.2.2. Problem Description

This tutorial considers the flow of two different isothermal polymer materials through a three-dimensional die. The problem consists of two inlets, one for each material, which are coextruded through the die to produce a multi-layered tubing structure at the exit.

Figure 3.7: Problem Description Inlet

Problem Description Inlet

The domain for the problem is divided into two parts, the first part is the die where the two fluid materials are confined between the inner wall and outer wall, as can be seen in Figure 3.7: Problem Description Inlet. The second part corresponds to the multi-material extrudate tubing which is in contact with the air and can deform freely along the inner and outer surfaces, as can be seen in Figure 3.8: Problem Description Exit.

The first inlet (inlet1) is positioned at the top of the die where the first material, high density polyethylene, enters into the domain at a flow rate of Q = 175 * 10 -6 m 3 / s. The second inlet (inlet 2) is positioned at the bottom of the die where the second material, linear low density polyethylene, enters into the domain at a flow rate of Q = 25 * 10 -6 m 3 / s.

Figure 3.8: Problem Description Exit

Problem Description Exit

At the end of the computational domain, it is assumed that the motion of the multi-material extruded tubing is fully developed and that it will not deform any further. It is assumed that the extrusion domain is long enough to account for the entire deformation of the extrudate.

3.2.3. Setup and Solution

3.2.3.1. Preparation

To prepare for running this tutorial:

  1. Prepare a working folder for your simulation.

  2. Download the coextrusion.zip file here .

  3. Unzip the coextrusion.zip file you have downloaded to your working folder.

  4. The mesh file coextruded_tubing.msh can be found in the unzipped folder.

3.2.3.2. Launching Ansys Fluent

  1. Use the Fluent Launcher to start Ansys Fluent.

  2. Select Materials Processing in the list of Fluent workspaces.

  3. Click Start.

3.2.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

  1. Open the geometry by reading in a Polyflow mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the mesh file from your working folder (coextruded_tubing.msh)

    3. Set the Mesh Length Unit to m, and close the dialog box.

  2. Set up an extrusion simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Template...

    This will display the Simulation Template dialog.

    1. Set the Type to Extrusion.

    2. Set the Goal to be Predict extrudate shape.

    3. Enter 2 for the Number of fluids.

    4. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.2.3.4. General Properties

Display and review the general settings for the problem setup.

  1. Select General in the Outline View to review general properties of the simulation.

  2. Review the settings in the properties page.

  3. Make sure that Steady is selected for the Calculation Type.

3.2.3.5. Material Properties

For this tutorial, you will replace the two fluids that were created for you with two materials imported from an installed material library.

  1. Import a new material from the materials library.

    Click the Import from library button in the Materials property page. Alternatively, you can right-click on the Materials node, and select Import from library from the context menu.

     Materials  Import from library

    This displays the Library of Materials dialog box.

    1. Select Extrusion for the Family.

    2. Select Extrusion_HDPE_isothermal_463K for the Material.

    3. Click OK to close the dialog.

      The new material will be added to the list of materials in the Outline View.

  2. Import another material from the library by repeating the same procedure.

    1. Select Extrusion for the Family

    2. Select Extrusion_LLDPE_isothermal_463K for the Material.

    3. Click OK to close the dialog.

      The new material will be added to the list of materials in the Outline View.

  3. Remove the original fluid materials (fluid-1 and fluid-2).

    1. Select fluid-1 in the Outline View, right-click and select Delete from the context menu.

       Setup Materials fluid-1  Delete

      You will be prompted to make sure you want to remove the material. Click Yes to remove the material.

    2. Perform the same operation for fluid-2 to remove it from your setup.

  4. Review material properties.

    1. Select Extrusion_HDPE_isothermal_463K under Materials in the Outline View to review material properties of the simulation.

    2. Review the density and the viscosity settings for the first material.

    3. Review the density and the viscosity settings for the Extrusion_LLDPE_isothermal_463K material.

3.2.3.6. Cell Zone Properties

In this portion of the tutorial, you will assign your fluid cell zones.

  1. Select fluid-zone under Cell Zones in the Outline View to review cell zone properties of the simulation.

  2. For Zones, select the field to open a selection dialog.

  3. In the selection dialog, select both available zones: die and extrudate and click OK to keep your selections and close the dialog.

  4. For Fluid Material(s), select the field to open a selection dialog.

  5. In the selection dialog, select the two materials previously created: Extrusion_HDPE_isothermal_463K and Extrusion_LLDPE_isothermal_463K and click OK to keep your selections and close the dialog.

3.2.3.7. Boundary Condition Properties

In the following steps you will set the conditions at each of the boundaries of the simulation.

  1. Set outlet conditions.

    1. Select extrudate-exit under Fluid Boundary Conditions in the Outline View to review the exit boundary conditions for the simulation.

    2. Select exit as the Boundary Zone.

    3. Select Take-up Velocity as the Flow Specification.

    4. Enter 0 m/s for Vx.

    5. Enter 0 m/s for Vy.

    6. Enter 0.5 m/s for Vz.

  2. Set free surface conditions.

    1. Select free-surface under Fluid Boundary Conditions in the Outline View to review free surface boundary conditions for the simulation.

    2. Select innersurface as the Boundary Zone.

    3. Select innerwall as the Fixed Part.

  3. Define a new free surface boundary.

    Click the New... button in the Fluid Boundary Conditions property page. Alternatively, you can right-click on the Fluid Boundary Conditions node, and select New... from the context menu.

     Fluid Boundary Conditions  New...

    1. Enter outer-free-surface for the Name.

    2. Select Free surface for the Type.

    3. Select outersurface for the Boundary Zone.

    4. Select outerwall for the Fixed Part.

  4. Set inlet conditions for the Extrusion_HDPE_isothermal_463K material.

    1. Select inlet-for-fluid-1 under Fluid Boundary Conditions in the Outline View to review inlet boundary conditions for the simulation.

    2. Select inlet1 for the Boundary Zone.

    3. Select Volume Flow for the Flow Specification.

    4. For the Incoming Material, select Extrusion_HDPE_isothermal_463K.

    5. Enter 175e-6 m^3/s for the Volume Flow Rate.

  5. Set inlet conditions for the Extrusion_LLDPE_isothermal_463K material.

    1. Select inlet-for-fluid-2 under Fluid Boundary Conditions in the Outline View to review inlet boundary conditions for the simulation.

    2. Select inlet2 for the Boundary Zone.

    3. Select Volume Flow for the Flow Specification.

    4. For the Incoming Material, select Extrusion_LLDPE_isothermal_463K.

    5. Enter 25e-6 m^3/s for the Volume Flow Rate.

  6. Set wall surface conditions.

    1. Select wall under Fluid Boundary Conditions in the Outline View to review free surface boundary conditions for the simulation.

    2. Select innerwall as the Boundary Zone.

    3. Select Partial slip for the Slip Specification.

    4. Retain the selection of Navier law for the Slip Model.

    5. Enter 1e6 kg/(m^2 s) for the Friction Coefficient.

  7. Define a new wall boundary.

    Click the New... button in the Fluid Boundary Conditions property page. Alternatively, you can right-click on the Fluid Boundary Conditions node, and select New... from the context menu.

     Fluid Boundary Conditions  New...

    1. Enter outer-wall for the Name.

    2. Select Wall for the Type.

    3. Select outerwall for the Boundary Zone.

    4. Select Partial slip for the Slip Specification.

    5. Retain the selection of Navier law for the Slip Model.

    6. Enter 1e6 kg/(m^2 s) for the Friction Coefficient.

3.2.3.8. Mesh Deformation Properties

This model involves two free surfaces, 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 these boundaries, so a mesh deformation technique is applied on this part of the mesh. The extrudate is entirely contained between the inner and outer surfaces, therefore this will be the only zone affected by the mesh deformation algorithm.

  1. Under Mesh Deformations, select extrudate to edit the default properties of the extrudate exit for your simulation.

  2. Select extrudate for Zones.

  3. Select Streamwise for the Extrudate Deformation Method.

  4. Select extrudate-die for the Inlet Section.

  5. Select exit for the Outlet Section.

3.2.3.9. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and specify other solution properties.

  1. Define the Calculation Activities.

    1. The number of processors is set to 4 by default, which is sufficient for most cases. However, for this problem we will increase the number of processors to 7 which will moderately increase the marginal speedup. Note that the value specified for the processor count does not need to be a power of 2.

      1. Enter 7 for the Number of Processors.

    2. The primary source of non-linearity in a coextrusion flow originates from the transport of immiscible materials and their respective properties. At each location in the calculation domain, equations must be inspected for the appropriate properties assignment. For facilitating the calculation, convergence strategies are available. Next to convergence strategy for multiple materials, which is often recommended for coextrusion flow calculations, we also enable convergence strategies for free surfaces and moving interfaces as well as for viscosity and slip.

      1. Enable the Multiple Materials option.

      2. Enable the Free Surfaces and Moving interfaces option.

      3. Enable the Viscosity and Slip option.

  2. Review your problem setup.

    1. Select Run Calculation to see available options.

    2. Click Check to review your simulation settings. The Fluent Materials Processing workspace will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  3. Calculate the solution.

    1. Click Calculate to begin the computation of the solution.

      Depending on your machine, the calculation should take approximately 20 to 25 minutes to complete.

    2. Once the calculations are complete, check the solution's listing in the Transcript tab to confirm its success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  4. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering coextruded-tubing as the name of the session.

3.2.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can setup contours, vectors, and so on.

  1. Display the fluid fraction distribution on the boundaries.

     ResultsGraphics Contour New

    Figure 3.9: Contour Properties

    Contour Properties

    The fluid fraction distibution (FLUIDFRACTION1) represents the proportion or fraction of fluid material within the multi-layered tubing.

    1. For the Name, enter coextrusion-fluid-fraction.

    2. For the Field, select FLUIDFRACTION1.

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Once your selections are complete, click OK to close the dialog.

      Figure 3.10: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. Keep the remaining defaults and click Display to see the fluid fraction contour plot in the graphics window.

      Figure 3.11: Contours of FLUIDFRACTION1

      Contours of FLUIDFRACTION1

In Figure 3.11: Contours of FLUIDFRACTION1 we can see the high density polyethylene in blue, and the linear low density polyethylene in red. The interface between these two polymers is shown in green. A more refined mesh in the zone where the fluid fraction varies will improve the interface location.

3.2.5. Summary

This tutorial introduced the concept of a 3D coextrusion problem with multiple materials using the Fluent Materials Processing workspace. You solved the problem using a specific 3D geometry for the die and made suitable assumptions about the physics of the problem.

3.3. 3D Polymer Blow Molding & Thermoforming

3.3.1. Introduction

This tutorial illustrates the simulation of a 3D blow molding process using the Fluent Materials Processing workspace (Fluent Materials Processing Workspace). The workspace will allow you to set up and solve a polymer blow molding problem so you can easily obtain an accurate prediction of the shaped object for a given mold geometry under prescribed operating conditions.

In this tutorial you will learn how to use the Fluent Materials Processing workspace to:

  • Easily set up your blow molding / thermoforming simulation

  • Calculate a solution

  • Analyze the results, such as generating thickness-based contour plots and running transient animations.

3.3.2. Problem Description

This tutorial simulates a typical thermoforming situation for a blister. Figure 3.12: Thermoforming of a Blister, Sheet (blue) and Mold (red) shows a view of the process in the initial configuration.

Figure 3.12: Thermoforming of a Blister, Sheet (blue) and Mold (red)

Thermoforming of a Blister, Sheet (blue) and Mold (red)

To reduce the computational run time, and utilizing the symmetric nature of the blister, only one quarter of the blister/mold is modeled, Figure 3.12: Thermoforming of a Blister, Sheet (blue) and Mold (red). From a geometric point of view, the initial (1/4) film has the following dimensions:

  • Length = 0.015 m

  • Width = 0.05 m

  • Initial thickness = 0.0005 m

This tutorial demonstrates the modeling of a polymer sheet made up of several layers, as it is often the case for packaging food or medication, where one layer provides mechanical resistance and another layer provides a humidity barrier.

The thickness of the layer, when 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.

The finite element mesh and the boundary conditions are displayed in Figure 3.13: Finite Element Mesh, Subdomains and Boundary Sets. 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.

The film has the following material properties:

  • First layer: Viscosity = 5500 Pa -s

    Second layer: Viscosity = 4000 Pa -s

  • First layer: Density = 975 kg/m3

    Second layer: Density = 1100 kg/m3

  • First layer: Initial thickness = 0.0002 m

    Second layer: Initial thickness = 0.0003 m

Figure 3.13: Finite Element Mesh, Subdomains and Boundary Sets

Finite Element Mesh, Subdomains and Boundary Sets

As seen in Figure 3.13: Finite Element Mesh, Subdomains and Boundary Sets, the topology involves two subdomains:

  • Subdomain 1 = film

  • Subdomain 2 = mold

and four boundary sets:

  • Boundary 1: will be fixed (clamped boundary)

  • Boundary 2: will be fixed (clamped boundary)

  • Boundary 3: symmetry boundary condition with respect to the x-axis

  • Boundary 4: symmetry boundary condition with respect to the y-axis

The inflation pressure will be defined on the subdomain representing the film (Subdomain 1).

When considering 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.

Often, both types of contact are encountered in a given application.

3.3.3. Setup and Solution

3.3.3.1. Preparation

To prepare for running this tutorial:

  1. Prepare a working folder for your simulation.

  2. Download the 3d_thermo_blister.zip file here .

  3. Unzip the 3d_thermo_blister.zip file you have downloaded to your working folder.

    The mesh file blister.msh can be found in the unzipped folder.

3.3.3.2. Launching Ansys Fluent

  1. Use the Fluent Launcher to start Ansys Fluent.

  2. Select Materials Processing in the list of Fluent workspaces.

  3. Click Start.

3.3.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

The workspace is a version of Ansys Fluent that utilizes the power of the Ansys Polyflow solver to simulate polymer flows such as extrusion, blow molding, pressing, and so forth.

  1. Open a blow molding geometry by reading in a Polyflow mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the mesh file from your working folder (blister.msh)

    3. In this case, the mesh file does not contain any unit information, so you are prompted to provide the units through the Define Mesh Unit dialog. Set the Mesh Length Unit to mm, and close the dialog box.


      Note:  The Fluent Materials Processing workspace always performs the setup and the calculations in MKS. If the mesh is provided in other units, then the coordinates are converted to meters prior to the calculation.


  2. Set up a blow molding simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Template...

    This will display the Simulation Template dialog.

    1. Set the Type to Blow Molding / Thermoforming.

    2. Keep the Number of Fixed Molds as 0.

    3. Set the Number of Moving Molds to 1.

    4. Set the Number of Layers to 2.

    5. Set the Duration to 0.5 seconds.

    6. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.3.3.4. General Properties

Click the General node in the Outline View to display general settings for the problem setup. For this tutorial, you can keep the default settings.

3.3.3.5. Material Properties

For this tutorial, you will replace the two fluids that were created for you with two special materials imported from an installed material library.

  1. Import a new material from the materials library.

    Click the Import from library button in the Materials property page. Alternatively, you can right-click on the Materials node, and select Import from library from the context menu.

     Materials  Import from library

    This displays the Library of Materials dialog box.

    1. Select Blow molding for the Family

    2. Select BlowMolding_HDPE_isothermal_463K for the Material.

    3. Click OK to close the dialog.

      The new material will be added to the list of materials in the Outline View.

  2. Import another material from the library by repeating the same procedure.

    1. Select Blow molding for the Family

    2. Select BlowMolding_PP_isothermal_493K for the Material.

    3. Click OK to close the dialog.

      The new material will be added to the list of materials in the Outline View.

  3. Remove the original fluid materials (fluid-1 and fluid-2).

    1. Select fluid-1 in the Outline View, right-click and select Delete from the context menu.

       Setup Materials fluid-1  Delete

      You will be prompted to make sure you want to remove the material. Click Yes to remove the material.

    2. Perform the same operation for fluid-2 to remove it from your setup.

  4. Review material properties.

    1. Select BlowMolding_HDPE_isothermal_463K under Materials in the Outline View to review material properties of the simulation.

      When you select an item in the Outline View, you can see and edit its properties and settings in the separate Properties View window. There are numerous material properties available for your materials using the workspace. Material properties are listed by category and you can expand the properties of the categories you are interested in by clicking the icon to the left of the category.

    2. Review the density and the viscosity settings for the first material.

    3. Review the density and the viscosity settings for the second material.

3.3.3.6. Cell Zone Properties

In this portion of the tutorial, you will assign which cell zone is the fluid which is the moving mold, and prescribe the motion of the mold using an expression.

  1. Set the cell zone conditions for the fluid zone.

    1. Select fluid-zone under Cell Zones in the Outline View to review cell zone properties of the simulation.

    2. For Zones, select the field to open a selection dialog. In the selection dialog, select subdomain1 and click OK to keep your selections and close the dialog.

    3. Set the Inflation Pressure to 10000 Pa.

      Click Display to visualize arrows on the fluid zone that demonstrate the orientation of the pressure.

    4. Select expression from the Time Dependency drop-down list, then enter IF(time <= 0.1, 0.0, IF(time >= 0.11, 1.0, (time-0.1)/(0.11-0.1))) for the time dependency.

  2. Set the cell zone conditions for the moving mold zone.

    1. Select moving-mold under Cell Zones in the Outline View to review cell zone properties of the simulation.

    2. For Zones, select the field to open a selection dialog. In the selection dialog, select subdomain2 and click OK to keep your selections and close the dialog.

    3. Check the orientation.

      Click Display to visualize arrows on the moving mold zone that demonstrate the orientation, and to confirm that the arrows are pointing to the mold body.

    4. Under Translation Velocity, select expression from the Vz drop-down list and enter 0.01 * IF(time <= 0.097, 1.0, IF(time >= 0.103, 0.0, 17.167 - time * 166.67)) for Vz.

3.3.3.7. Layer Properties

Using the blow molding simulation template, this tutorial defined two layers. In the following steps you will further define the various layers of the simulation, each having its own material properties and initial thickness (which can be defined as a constant, or a linear function of X, Y, Z coordinates. Layers can overlap (partly or entirely), but they are not required to do so. Here, each layer is defined on the same fluid-zone. The solver will model the deformation of the polymer sheet that combines (the resistance of) the two layers. The layers are meant to stick together, so the deformation of each layer is identical, since the area stretch can only be computed once.

  1. Set conditions for the first layer.

    1. Select layer-of-fluid-1 under Layers in the Outline View to review this layer's settings for the simulation.

    2. Select BlowMolding_HDPE_isothermal_463K as the Material.

    3. Select subdomain1 as the Zones.

    4. Select f(X,Y,Z) under Option.

    5. The layer's thickness H is defined by the polynomial: A*(1+B*x+C*y+D*z). All that is required for this tutorial is a value for the first coefficient, A, so enter 0.0002 m.

  2. Set conditions for the second layer.

    1. Select layer-of-fluid-2 under Layers in the Outline View to review this layer's settings for the simulation.

    2. Select BlowMolding_PP_isothermal_493K as the Material.

    3. Select subdomain1 as the Zones.

    4. Select f(X,Y,Z) under Option.

    5. Enter 0.0003 m for the thickness coefficient A.

3.3.3.8. Fluid Boundary Condition Properties

In the following steps you will set the conditions at each of the boundaries of the simulation.

  1. Set conditions for the fixed edges.

    1. Select fixed-edges under Fluid Boundary Conditions in the Outline View to review this boundary condition for the simulation.

    2. Select boundary1 and boundary2 as the Boundary Zone.

  2. Set two symmetry conditions.

    1. Select Fluid Boundary Conditions in the Outline View, and click the New... button in its properties window, or right-click in the tree and add a new condition to the tree.

    2. Rename the new condition to symmetry-1.

    3. Set the Type to Symmetry.

    4. Select boundary3 as the Boundary Zone.

    5. Select Normal direction along X axis as the Plane of Symmetry.

    6. Repeat the process and assign a symmetry boundary called symmetry-2 to boundary4, assigning its plane of symmetry to be Normal direction along Y axis.

3.3.3.9. Contact Boundary Conditions

Contact between the fluid film and the solid mold can be accounted for using contact conditions.

  1. Under Contact Boundary Conditions, select contact-with-moving-mold to edit the default properties of this contact condition for your simulation.

  2. Select subdomain1 for the Fluid Zone of Contact.

  3. Select subdomain2 for the Mold Zone of Contact.

3.3.3.10. Mesh Deformation Properties

Here, we'll assign the film cell zone to be a deforming fluid.

  1. Under Mesh Deformations, select fluid-deformation to edit the default properties of the fluid deformation zone for your simulation.

  2. Select subdomain1 for Zones.

3.3.3.11. Adaptive Mesh Settings

In many practical cases, the use of adaptive meshing based on contact, remeshing, or both may be useful to selectively and automatically refine the mesh during the solution. The blow molding simulation template automatically includes adaptive mesh refinement, however, it is not necessary for this tutorial.

Select Adaptive Meshing in the Outline View and disable it in the properties window.

3.3.3.12. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and other solution properties. Most items indicate that current default values are appropriate.

  1. Review the derived quantities.

    1. Select Derived Quantities to see available options.

    2. Under Extension, leave the Area Stretch Ratio enabled and retain all other default settings.

  2. Review your problem setup.

    1. Select Run Calculation to see available options.

    2. Click Check to review your simulation settings. Fluent will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  3. Calculate a solution.

    1. Click Calculate to start the calculation.

    2. Once the calculations are complete, check the solution's listing in the Transcript tab to confirm it's success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  4. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering blister as the name of the session

3.3.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can setup contours, vectors, and so on.

  1. Improve the view of the geometry using mirror planes.

    1. Open the Results tab of the Ribbon and, under Graphics, select Mirror Planes... to open the Mirror Planes dialog.

       ResultsGraphics Mirror Planes...

    2. In the Mirror Planes dialog, select X and Y and click Close to close the dialog.

  2. Display contours of thickness (THICKNESS_1) in the fluid region (subdomain1).

     ResultsGraphics Contour New

    Figure 3.14: Thickness Contour Properties

    Thickness Contour Properties

    1. Specify a new Name (Thickness1).

    2. For the Field, select THICKNESS_1.

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Select subdomain1 from the list, and click OK to close the dialog.

      Figure 3.15: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. For the contour properties, keep the remaining defaults and click Display to see the thickness contour plot in the Graphics window.

      Use the toolbars in the Graphics window to adjust the display so that it is in an isometric orientation, and fits the window.

      Figure 3.16: Contours of Thickness (THICKNESS_1)

      Contours of Thickness (THICKNESS_1)


    5. Review an animation of the contour.

       ResultsTimestep

      1. To reset the timesteps, change the Current animation time step value to 1. Alternatively, you can also move the slider all the way to the left, or click the Jump to Start button ( ).

      2. To start the animation., click the Play button ( ) to start the animation.

      3. To stop the animation., click the Stop button ( ).

        Figure 3.17: Contours of Thickness (THICKNESS_1): Timestep =1

        Contours of Thickness (THICKNESS_1): Timestep =1

        Figure 3.18: Contours of Thickness (THICKNESS_1): Timestep =25

        Contours of Thickness (THICKNESS_1): Timestep =25

        Figure 3.19: Contours of Thickness (THICKNESS_1): Timestep=50

        Contours of Thickness (THICKNESS_1): Timestep=50

  3. Display contours of thickness (THICKNESS_2) in the fluid region (subdomain1).

     ResultsGraphics Contour New

    Repeat the previous steps to create a new contour over the subdomain1 surface with the Name of Thickness2 with the Field set to THICKNESS_2.

    Figure 3.20: Contours of Thickness (THICKNESS_2)

    Contours of Thickness (THICKNESS_2)

    Repeat the previous steps to review an animation of the contour.

  4. Display contours of thickness (THICKNESS_TOT) in the fluid region (subdomain1).

     ResultsGraphics Contour New

    Repeat the previous steps to create a new contour over the subdomain1 surface with the Name of Thickness with the Field set to THICKNESS_TOT.

    Figure 3.21: Contours of Thickness (THICKNESS_TOT)

    Contours of Thickness (THICKNESS_TOT)

    Repeat the previous steps to review an animation of the contour.

  5. Display contours of area stretch (AREA_STRETCH) in the fluid region (subdomain1).

     ResultsGraphics Contour New

    Repeat the previous steps to create a new contour over the subdomain1 surface with the Name of AreaStretch with the Field set to AREA_STRETCH.

    Figure 3.22: Contours of Area Stretch (AREA_STRETCH)

    Contours of Area Stretch (AREA_STRETCH)

    Repeat the previous steps to review an animation of the contour.

3.3.5. Summary

This tutorial introduced the concept of a 3D blow molding problem using the Fluent Materials Processing workspace. You solved the problem using a 3D shell geometry for the mold and the film and made suitable assumptions about the physics of the 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.

3.4. Injection Stretch Blow Molding (ISBM)

3.4.1. Introduction

This tutorial illustrates the simulation of a 3D injection stretch blow molding (ISBM) process using the Fluent Materials Processing workspace (Fluent Materials Processing Workspace). The workspace will allow you to set up and solve a polymer ISBM problem so you can easily obtain an accurate prediction of the shaped object for a given mold geometry under prescribed operating conditions.

In this tutorial you will learn how to use the Fluent Materials Processing workspace to:

  • Easily set up your injection stretch blow molding (ISBM) simulation.

  • Calculate a solution.

  • Analyze the results, such as generating thickness-based contour plots and running transient animations.

3.4.2. Problem Description

This tutorial covers the injection stretch blow molding (ISBM) process for a 3D bottle. A preform is first set into a mold and stretched by means of a moving rod. Once the rod has reached the desired position (near the bottom of the mold), an internal pressure is applied to the preform to obtain the desired bottle shape. The initial configuration of the model is shown in Figure 3.23: Problem Description and consists of the following three domains:

  • preform (displayed in orange)

  • mold (displayed in light blue)

  • rod (displayed in light blue)

Figure 3.23: Problem Description

Problem Description

The dimensions of the domains listed above are as follows:

  • Rod diameter = 0.02 m

  • Rod height = 0.1 m

  • Mold diameter = 0.097 m

  • Mold height = 0.3 m

  • Preform initial thickness = 0.005 m

The grooves which are visible on the mold suggest that the model does not involve any geometric symmetry.

For this problem, contact release between the preform and the rod is assumed upon application of the shaping pressure. During the tutorial setup, the contact release will be defined using contact boundary conditions and cell zone and fluid boundary conditions will be specified to account for the following:

  1. The velocity-driven motion of the rod during a prescribed time interval (until the rod has reached the specified displacement).

  2. The applied shaping pressure.

Additionally, the following properties will be specified during the tutorial setup:

  • Downward velocity of the rod = 0.246 m/s

    This velocity ensures that the rod reaches a position of 0.2 m into the mold within 0.8 s, prior to applying the internal pressure to the preform.

  • Preform material density = 1000 kg/m3

  • Applied shaping pressure = 106 Pa

    The shaping pressure will be defined to apply after 0.86 seconds, once the rod reaches the desired position within the preform.

  • Adhesion Force = 10 Pa

    This value will be used to express the contact release of the preform from the rod.

  • Viscosity = 105 Pa s

3.4.3. Setup and Solution

3.4.3.1. Preparation

To prepare for running this tutorial:

  1. Download the ISBM.zip file here .

  2. Unzip ISBM.zip to your working directory.

  3. The file, ISBM.msh can be found in the folder.

3.4.3.2. Launching Ansys Fluent

  1. Use the Fluent Launcher to start Ansys Fluent.

  2. Select Materials Processing in the list of Fluent workspaces.

  3. Click Start.

3.4.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

The workspace is a version of Ansys Fluent that utilizes the power of the Ansys Polyflow solver to simulate polymer flows such as extrusion, blow molding, pressing, and so forth.

  1. Load the provided mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the mesh file from your working folder (ISBM.msh).

    3. Set the Mesh Length Unit to m, and click Accept to close the dialog box.

  2. Set up a blow molding simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Template...

    This will display the Simulation Template dialog.

    1. Set the Type to Blow Molding / Thermoforming.

    2. Set the Number of Fixed Molds to 1.

    3. Set the Number of Moving Molds to 1.

    4. Keep the Number of Layers as 1.

    5. Set the Duration to 1.05 seconds.

    6. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.4.3.4. General Properties

Click the General node in the Outline View to display general settings for the problem setup. For this tutorial, you can keep the default settings.

3.4.3.5. Material Properties

  1. Edit the density and viscosity settings for the fluid material.

     Setup Materials fluid

    1. Expand the Density Law and Viscosity Law categories.

    2. In the Density Law category, enter 1000 kg/m3 for the Density.

    3. In the Viscosity Law category, enter 1e5 Pa s for the Viscosity.

3.4.3.6. Cell Zone Properties

In this portion of the tutorial, you will assign your fluid cell zones.

  1. Define the cell zone properties for the fluid-zone.

     Setup Cell Zones fluid-zone

    1. For Zones, select the field to open a selection dialog.

    2. In the selection dialog, select preform and click OK to keep the selection and close the dialog.

    3. Retain the default selection of Generalized newtonian for Fluid Model .

    4. Enter -1e6 Pa for the Inflation Pressure.

    5. Click the arrow next to the Time Dependency entry and select Expression from the drop-down list, then enter the following expression for the time dependency:

      IF(time <= 0.86, 0., IF(time >= 0.9, 1., 0.+(1.-0.)/(0.9-0.86)*(time-0.86)))

      This expression corresponds to a ramp function which increases from 0 to 1 within the time interval of 0.86 to 0.9.

    6. Click Display to visualize arrows on the fluid zone.

      The arrows displayed on the fluid-zone show the orientation of the pressure.

  2. Define the cell zone properties for the fixed-mold.

     Setup Cell Zones fixed-mold

    1. For Zones, select mold.

    2. Click Display to visualize arrows on the fixed mold.

      The arrows displayed on the fixed-mold show that the direction is pointing to the mold body side.

  3. Define the cell zone properties for the moving-mold.

     Setup Cell Zones moving-mold

    1. For Zones, select rod.

    2. Click Display to visualize the arrows on the rod (moving-mold).

    3. In the Graphics window, rotate the rod to view the arrows from the +Y direction.

      The arrows displayed on the moving-mold point away from the rod, they should point to the rod axis.

    4. Enable the Flip Darts Orientation option to make the arrows point towards the rod.

      In the above image, the darts may appear to still be facing away from the rod. However, the darts are pointing towards the rod axis, but are overlapping with the rod geometry due to scaling.

    5. Click the arrow next to the Vy entry and select Expression from the drop-down list, then enter the following expression for Vy:

      IF(time <= 0.8, -0.246, IF(time >= 0.85, 0., -0.246+(0.+0.246)/(0.85-0.8)*(time-0.8)))

      This expression corresponds to a ramp function which increases from -0.246 to 0 within the time interval of 0.8 to 0.85.

3.4.3.7. Layer Properties

  1. Define the layer-of-fluid properties.

     Setup Layers layer-of-fluid

    1. For Zones, select the field to open a selection dialog.

    2. In the selection dialog, select preform and click OK to keep the selection and close the dialog.

    3. Enter 0.005 m as the thickness for A.

3.4.3.8. Boundary Condition Properties

  1. Set the conditions for the fixed-edges boundary.

    1. Select fixed-edges under Fluid Boundary Conditions in the Outline View to review this boundary condition for the simulation.

    2. Select preformedge as the Boundary Zone.

3.4.3.9. Contact Boundary Conditions

Contact between the preform and the solid and moving molds can be accounted for using contact conditions.

  1. Define the contact between the preform and the fixed mold.

    1. Under Contact Boundary Conditions, select contact-with-fixed-mold to edit the default properties of this contact condition for your simulation.

    2. Select preform for the Fluid Zone of Contact.

    3. Select mold for the Mold Zone of Contact.

    4. Enter 1e-4 m for the Penetration Accuracy.

  2. Define the contact between the preform and the moving mold.

    1. Under Contact Boundary Conditions, select contact-with-moving-mold to edit the default properties of this contact condition for your simulation.

    2. Select preform for the Fluid Zone of Contact.

    3. Select rod for the Mold Zone of Contact.

    4. Enter 1e-4 m for the Penetration Accuracy.

    5. Enable the Allow Contact Release option.

    6. Enter 10 Pa for Adhesion Force per Unit Area.

3.4.3.10. Mesh Deformation Properties

Here, we will assign the preform cell zone to be a deforming fluid.

  1. Under Mesh Deformations, select fluid-deformation to edit the default properties of the fluid deformation zone for your simulation.

  2. Retain the selection of Lagrangian as the deformation Type.

  3. Select preform for Zones.

3.4.3.11. Adaptive Mesh Settings

Define the conditions and settings for adaptive meshing that will be used to refine the mesh during the solution.

  1. Define the adaptive meshing condition (condition-1).

     Setup Adaptive Meshing Conditions condition-1

    1. Select preform for Zones.

    2. Select mold for Mold Surfaces of Contact.

    3. Retain the selection of angle and curvature for Method.

    4. Enter 0.005 m for Minimum Size.

    5. Enter 0.2 m for Maximum Size.

    6. Enter 0.001 m for Tolerance.

    7. Enter 0.02 m for Critical Distance.

  2. Edit the adaptive meshing properties.

     Setup Adaptive Meshing

    1. Enter 5 for the Number of Steps.

3.4.3.12. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and other solution properties. Most items indicate that current default values are appropriate.

  1. Edit the Derived Quantities.

     Solution Derived Quantities

    1. Under Extension, enable the Track Vectors option to expose the additional options.

    2. For Initial Orientation, select initially parallel to D and initially perpendicular to D.

    3. Retain the entry of 0 m for reference points Px, Py, and Pz.

    4. Enter 0 for reference direction Dir-x.

    5. Enter 1 for reference direction Dir-y.

    6. Enter 0 for reference direction Dir-z.

    7. Enable Self-Contact and retain the selection of Warning for Action if Self-Contact Detected

  2. Edit the Calculation Activities.

     Solution Calculation Activities

    1. Select Calculation Activities to see available options.

    2. From the Convergence Criterion drop-down list, select User specified.

    3. Enter 0.0001 for the Convergence Tolerance.

    4. Enter 1e-4 seconds for the Initial Time Step.

    5. Enter 1e-6 seconds for the Minimum Time Step.

    6. Enter 1e-2 seconds for the Maximum Time Step.

    7. Enter 0.001 for the Tolerance.

  3. Review your problem setup.

     Solution Run Calculation

    1. Click Check to review your simulation settings. Fluent will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  4. Calculate a solution.

    1. Click Calculate to start the calculation.

    2. Click the Plots tab to view the plot of the convergence monitors.

    3. Once the calculations are complete, check the solution's listing in the Transcript tab to confirm it's success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  5. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering ISBM as the name of the session

3.4.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can setup contours, vectors, and so on.

  1. Display contours of thickness (THICKNESS) for the preform surface.

     ResultsGraphics Contour New

    Figure 3.24: Thickness Contour Properties

    Thickness Contour Properties

    1. Specify a new Name (Thickness).

    2. For the Field, select THICKNESS.

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Select preform from the list, and click OK to close the dialog.

      Figure 3.25: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. For the contour properties, keep the remaining defaults and click Display to see the thickness contour plot in the graphics window.

      Figure 3.26: Contours of Thickness (THICKNESS)

      Contours of Thickness (THICKNESS)

  2. Review an animation of the contour of THICKNESS.

     ResultsTimestep

    1. To reset the timesteps, change the Current animation time step value to 0. Alternatively, you can also move the slider all the way to the left, or click the Jump to Start button ( ).

    2. To start the animation, click the Play button ( ).

    3. To stop the animation, click the Stop button ( ).

      Figure 3.27: Contours of Thickness (THICKNESS): Timestep =0

      Contours of Thickness (THICKNESS): Timestep =0

      Figure 3.28: Contours of Thickness (THICKNESS): Timestep =208

      Contours of Thickness (THICKNESS): Timestep =208

      Figure 3.29: Contours of Thickness (THICKNESS): Timestep=417

      Contours of Thickness (THICKNESS): Timestep=417

      Figure 3.30: Contours of Thickness (THICKNESS): Timestep=625

      Contours of Thickness (THICKNESS): Timestep=625

  3. Display contours of area stretch (AREA_STRETCH) for the preform surface.

     ResultsGraphics Contour New...

    Repeat the previous steps to create a new contour over the preform surface with the Name of AreaStretch with the Field set to AREA_STRETCH.

    Figure 3.31: Contours of Area Stretch (AREA_STRETCH)

    Contours of Area Stretch (AREA_STRETCH)

    Repeat the previous steps to review an animation of the contour.

  4. Display contours of self-contact (SELF_CONTACT) for the preform surface.

     ResultsGraphics Contour New...

    Repeat the previous steps to create a new contour over the preform surface with the Name of Self-contact with the Field set to SELF_CONTACT.

    Figure 3.32: Contours of Self-contact (SELF_CONTACT)

    Contours of Self-contact (SELF_CONTACT)

    Repeat the previous steps to review an animation of the contour.

3.4.5. Summary

This tutorial introduced the concept of a 3D injection stretch blow molding problem using the Fluent Materials Processing workspace. You solved the problem using a 3D geometry for the mold, the preform, and the moving rod and made suitable assumptions about the physics of the problem. The preform was lowered into a 3D mold using a moving rod, then an internal pressure was applied to obtain the desired bottle shape.

3.5. Glass Pressing

3.5.1. Introduction

This tutorial illustrates the simulation of a glass pressing process using the Fluent Materials Processing workspace (Fluent Materials Processing Workspace). The workspace will allow you to set up and solve a glass pressing problem so you can easily obtain an accurate prediction of the shaped glass object for a given mold geometry under prescribed operating conditions.

In this tutorial you will learn how to use the Fluent Materials Processing workspace to:

  • Easily set up your pressing simulation.

  • Calculate a solution.

  • Analyze the results, such as generating contour plots and running transient animations.

3.5.2. Problem Description

In this tutorial, we simulate a simple isothermal glass pressing process. In particular, deformations undergone by the fluid as well as contact evolution are investigated. Due to the large deformations involved in the pressing process, adaptive remeshing will also be applied for the fluid zone.

The number of steps between each mesh adaption will be set to 5. This value is chosen to allow the time step adaption algorithm to reduce or increase the time step according to the accuracy of the transient scheme.

Figure 3.33: Problem Description shows a sketch of the geometry in the initial configuration, before pressing the fluid. A gob of material is placed onto a mold and at a specified time, the plunger moves downwards and presses the fluid volume. Due to the symmetry of the molds and the fluid gob, only a quarter of the geometry is modelled.

Figure 3.33: Problem Description

Problem Description

From a geometrical point of view, the mesh within three volumes is made of tetrahedra. In the initial configuration, the plunger is at rest and has a mass of 100 kg. At the start of the process, a constant downward force of 10 N is applied. As the process progresses, the downward velocity of the plunger varies to maintain the applied force against the viscous forces of the fluid, which resist deformation. Additionally, the plunger has a maximum displacement of 0.027 m.

The fluid has a constant viscosity of 500 Pa s and density of 2500 kg/m3. It is assumed that the fluid material sticks on walls while establishing contact with the plunger and the fixed mold. The contact condition between the fluid and both the plunger and fixed mold will be defined with a penetration accuracy of 1e-4 m.

3.5.3. Setup and Solution

3.5.3.1. Preparation

To prepare for running this tutorial:

  1. Download the glass_pressing.zip file here .

  2. Unzip glass_pressing.zip to your working directory.

  3. The file, glass_pressing.msh can be found in the folder.

3.5.3.2. Launching Ansys Fluent

  1. Use the Fluent Launcher to start Ansys Fluent.

  2. Select Materials Processing in the list of Fluent workspaces.

  3. Click Start.

3.5.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

The workspace is a version of Ansys Fluent that utilizes the power of the Ansys Polyflow solver to simulate flows of complex rheology fluids and covers processes such as extrusion, blow molding, pressing, and so forth.

  1. Load the provided mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the mesh file from your working folder (glass_pressing.msh).

    3. Set the Mesh Length Unit to cm, and click Accept to close the dialog box.

  2. Set up a pressing simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Template...

    This will display the Simulation Template dialog.

    1. Set the Type to Pressing.

    2. Set the Number of Fixed Molds to 1.

    3. Set the Number of Moving Molds to 1.

    4. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.5.3.4. General Properties

Click the General node in the Outline View to display general settings for the problem setup. For this tutorial, you can keep the default settings.

3.5.3.5. Material Properties

  1. Edit the density and viscosity settings for the fluid material.

     Setup Materials fluid

    1. Enter glass for the Name.

    2. Expand the Density Law and Viscosity Law categories.

    3. In the Density Law category, enter 2500 kg/m3 for the Density.

    4. In the Viscosity Law category, enter 500 Pa s for the Viscosity.

3.5.3.6. Cell Zone Properties

  1. Define the cell zone properties for the fluid-zone.

     Setup Cell Zones fluid-zone

    1. For Zones, select the field to open a selection dialog.

    2. In the selection dialog, select glass and click OK to keep the selection and close the dialog.

    3. Retain the default selection of Generalized newtonian for Fluid Model .

  2. Define the cell zone properties for the fixed-mold.

     Setup Cell Zones fixed-mold

    1. For Zones, select the field to open a selection dialog.

    2. In the selection dialog, select mold and click OK to keep the selection and close the dialog.

  3. Define the cell zone properties for the moving-mold.

     Setup Cell Zones moving-mold

    1. For Zones, select the field to open a selection dialog.

    2. In the selection dialog, select plunger and click OK to keep the selection and close the dialog.

    3. From the Motion Type drop-down list, select Translation force imposed.

    4. Enter 100 kg for Mass of the Mold.

    5. Enable the Limit Mold Displacement option.

    6. Enter 0.027 m for Maximum Displacement.

    7. Enter -10 N for Fy.

3.5.3.7. Boundary Condition Properties

  1. Set the conditions for the free-surface boundary.

    1. Select free-surface under Fluid Boundary Conditions in the Outline View to review this boundary condition for the simulation.

    2. Select freesurface as the Boundary Zone.

  2. Click Fluid Boundary Conditions in the Outline View, then click New... in the properties window.

    1. Enter Symmetry-1 for Name.

    2. From the Type drop-down list, select Symmetry.

    3. Select symmetry1 as the Boundary Zone.

  3. Click Fluid Boundary Conditions in the outline view, then click New... in the properties window.

    1. Enter Symmetry-2 for Name.

    2. From the Type drop-down list, select Symmetry.

    3. Select symmetry2 as the Boundary Zone.

3.5.3.8. Contact Boundary Conditions

Contact between the glass and the solid and moving molds can be accounted for using contact conditions.

  1. Define the contact between the fluid zone and the fixed mold.

    1. Under Contact Boundary Conditions, select contact-with-fixed-mold to edit the default properties of this contact condition for your simulation.

    2. Select freesurface for the Fluid Zone of Contact.

    3. Select moldcontact for the Mold Zone of Contact.

    4. Enter 1e-4 m for the Penetration Accuracy.

  2. Define the contact between the fluid zone and the moving mold.

    1. Under Contact Boundary Conditions, select contact-with-moving-mold to edit the default properties of this contact condition for your simulation.

    2. Select freesurface for the Fluid Zone of Contact.

    3. Select plungercontact for the Mold Zone of Contact.

    4. Enter 1e-4 m for the Penetration Accuracy.

3.5.3.9. Mesh Deformation Properties

Here, we will assign the glass to be a deforming fluid.

  1. Under Mesh Deformations, select fluid-deformation to edit the default properties of the fluid deformation zone for your simulation.

  2. Retain the selection of Lagrangian as the deformation Type.

  3. Select glass for Zones.

3.5.3.10. Adaptive Mesh Settings

Define the conditions and settings for adaptive meshing that will be used to refine the mesh during the solution.

  1. Define the adaptive meshing condition (condition-1).

     Setup Adaptive Meshing Conditions condition-1

    1. Select glass for Zones.

    2. Select moldcontact and plungercontact for Mold Surfaces of Contact.

    3. Select distance from the Method drop-down list.

    4. Enter 0.001 m for Minimum Size.

    5. Enter 0.001 m for Minimum Distance.

    6. Enter 0.002 m for Maximum Size.

    7. Enter 0.002 m for Maximum Distance.

  2. Edit the adaptive meshing properties.

     Setup Adaptive Meshing

    1. Enter 5 for the Number of Steps.

3.5.3.11. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and other solution properties. Most items indicate that current default values are appropriate.

  1. Edit the Calculation Activities.

     Solution Calculation Activities

    1. From the Convergence Criterion drop-down list, select User specified.

    2. Enter 0.0001 for the Convergence Tolerance.

  2. Review your problem setup.

     Solution Run Calculation

    1. Click Check to review your simulation settings. Fluent will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  3. Calculate a solution.

    1. Click Calculate to start the calculation.

    2. Click the Plots tab to view the plot of the convergence monitors.

    3. Once the calculations are complete, check the solution's listing in the Transcript tab to confirm it's success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  4. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering glass_pressing as the name of the session

3.5.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can setup contours, vectors, and so on.

  1. Display contours of contact time (CONTACT_TIME_2) for the freesurface, symmetry1, and symmetry2 surfaces.

    The CONTACT_TIME_2 field corresponds with the length of time that the plunger is in contact with the fluid.

     ResultsGraphics Contour New

    Figure 3.34: Contact Time Contour Properties

    Contact Time Contour Properties

    1. Specify a new Name (contact-time-2).

    2. For the Field, select CONTACT_TIME_2.

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Select freesurface, symmetry1, and symmetry2 from the list, and click OK to close the dialog.

      Figure 3.35: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. For the contour properties, keep the remaining defaults and click Display to see the thickness contour plot in the graphics window.

      Figure 3.36: Contours of Contact Time (CONTACT_TIME_2)

      Contours of Contact Time (CONTACT_TIME_2)

  2. Review an animation of the contour of CONTACT_TIME_2.

     ResultsTimestep

    1. To reset the timesteps, change the Current animation time step value to 1. Alternatively, you can also move the slider all the way to the left, or click the Jump to Start button ( ).

    2. To start the animation, click the Play button ( ) to start the animation.

    3. To stop the animation, click the Stop button ( ).

      Figure 3.37: Contours of Contact Time (CONTACT_TIME_2): Timestep =1

      Contours of Contact Time (CONTACT_TIME_2): Timestep =1

      Figure 3.38: Contours of Contact Time (CONTACT_TIME_2): Timestep =42

      Contours of Contact Time (CONTACT_TIME_2): Timestep =42

      Figure 3.39: Contours of Contact Time (CONTACT_TIME_2): Timestep=84

      Contours of Contact Time (CONTACT_TIME_2): Timestep=84

      Figure 3.40: Contours of Contact Time (CONTACT_TIME_2): Timestep=126

      Contours of Contact Time (CONTACT_TIME_2): Timestep=126

  3. Display contours of initial coordinate y (COORDINI[Y]).

     ResultsGraphics Contour New

    Repeat the previous steps to create a new contour over the freesurface, symmetry1, and symmetry2 surfaces with the Name of initial-coordinate-y and the Field set to COORDINI[Y].

    Figure 3.41: Contours of Initial Coordinate Y (COORDINI[Y])

    Contours of Initial Coordinate Y (COORDINI[Y])

    Repeat the previous steps to review an animation of the contour.

  4. Display contours of initial coordinate x (COORDINI[X]).

     ResultsGraphics Contour New

    Repeat the previous steps to create a new contour over the freesurface, symmetry1, and symmetry2 surfaces, with Name of initial-coordinate-x, Field set to COORDINI[X], and the Use Global Range option disabled.

    Figure 3.42: Contours of Initial Coordinate X (COORDINI[X])

    Contours of Initial Coordinate X (COORDINI[X])

    Repeat the previous steps to review an animation of the contour.

3.5.5. Summary

This tutorial introduced the concept of a 3D pressing problem using the Fluent Materials Processing workspace. You solved the problem using a 3D geometry for the fixed mold and the plunger and made suitable assumptions about the physics of the problem. The plunger moves downward and presses the fluid against a fixed 3D mold to obtain the desired shape.

3.6. 3D Inverse Extrusion

3.6.1. Introduction

This tutorial illustrates the simulation of an extrusion process with adaption of a die geometry to obtain the desired extrudate shape, using the Fluent Materials Processing workspace (Fluent Materials Processing Workspace).

In this tutorial you will learn how to use the Fluent Materials Processing workspace to:

  • Easily set up your inverse extrusion simulation.

  • Calculate a solution.

  • Analyze the results.

3.6.2. Problem Description

This tutorial covers the simulation of a non-isothermal extrusion process to predict the shape of a die for a given extrudate shape. The initial configuration of the die and extrudate geometry is shown in Figure 3.43: Problem Description, which consists of four distinct cell zones:

  • extrudate

    Cell zone for the extrudate geometry.

  • constant

    Cell zone for the constant portion of the die.

  • adaptive

    Cell zone for the adaptive portion of the die.

  • fixed

    Cell zone for the fixed portion of the die.

Figure 3.43: Problem Description

Problem Description

The extrudate cell zone contains the following boundaries:

The constant cell zone contains the following boundaries:

The adaptive cell zone contains the following boundaries:

The fixed cell zone contains the following boundaries:

The figures below show the relevant boundaries for the extrudate, constant, adaptive, and fixed cell zones.

Figure 3.44: Boundaries

Boundaries

Figure 3.44: Boundaries shows the free-surface displayed in pink and the extrudate-exit and inlet displayed with red arrows and blue arrows, respectively. Note that walls (displayed in gray) consists of the wall-constant, wall-adaptive, and wall-fixed boundaries.

Figure 3.45: Constant, Adaptive, and Fixed Walls

Constant, Adaptive, and Fixed Walls

Figure 3.45: Constant, Adaptive, and Fixed Walls shows boundaries wall-constant, wall-adaptive, and wall-fixed.

Figure 3.46: Extrudate Guide

Extrudate Guide

Figure 3.46: Extrudate Guide shows the guide within the extrudate cell zone.

Figure 3.47: Symmetries

Symmetries

Figure 3.47: Symmetries shows the boundaries that support a symmetry condition: symmetry-extrudate, symmetry-constant, symmetry-adaptive, and symmetry-fixed.

Convergence can be difficult to achieve when computing inverse extrusion problems, especially for non-isothermal cases. However, the Fluent Materials Processing workspace is equipped with convergence strategies for handling such difficulties. The appropriate convergence strategies are activated automatically by the wizard, used in this tutorial.

3.6.3. Setup and Solution

3.6.3.1. Preparation

To prepare for running this tutorial:

  1. Prepare a working folder for your simulation.

  2. Download the inverse_extrusion.zip file here .

  3. Unzip the inverse_extrusion.zip file you have downloaded to your working folder.

  4. The polyflow mesh file invext.poly can be found in the unzipped folder.

3.6.3.2. Launching Ansys Fluent

  1. Use the Fluent Launcher to start Ansys Fluent.

  2. Select Materials Processing in the list of Fluent workspaces.

  3. Click Start.

3.6.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

The workspace is a version of Ansys Fluent that utilizes the power of the Ansys Polyflow solver to simulate polymer flows such as extrusion, blow molding, pressing, and so forth.

  1. Open an extrusion geometry by reading in a Polyflow mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the polyflow mesh file from your working folder (invext.poly).

  2. Use the wizard to set up an inverse extrusion simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Wizard...

    This will display the Wizard dialog.

    1. Set the Type to Simple inverse extrusion.

    2. For Fluid Material select Extrusion_PS_thermal_453K.

    3. Enter 0.085 kg/s for the Mass Flow Rate.

    4. For Die Type select Fixed/Adaptive/Constant Sections.

    5. Enable the Thermal option.

    6. Enter 423 K and 373 K for the Inlet Temperature and Wall Temperature, respectively.

    7. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.6.3.4. General Properties

Click the General node in the Outline View to display general settings for the problem setup.

  1. Select General in the Outline View to review general properties of the simulation.

  2. Review the settings in the properties page.

3.6.3.5. Boundary Condition Properties

In the following steps you will set the conditions at each of the boundaries of the simulation.

  1. Set the free surface conditions.

    1. Select free-surface under Fluid Boundary Conditions in the Outline View to review free surface boundary conditions for the simulation.

    2. Select Convection for the thermal condition Option.

    3. Enter 5 W/(m2K) for the Heat Transfer Coefficient.

    4. Enter 303 K for the Convection Temperature.

  2. Set the wall conditions.

    1. Select wall under Fluid Boundary Conditions in the Outline View to review wall boundary conditions for the simulation.

    2. Select Partial Slip for Slip Specification as the Boundary Zone.

    3. Enter 7e6 kg/(m2s) for the Friction Coefficient.

    4. Select Convection for the thermal condition Option.

    5. Enter 15 W/(m2K) for the Heat Transfer Coefficient.

    6. Enter 373 K for the Convection Temperature.

  3. Create the boundary condition for the guide.

    1. Select Fluid Boundary Conditions in the Outline View, and click the New... button in the properties window, or right-click in the tree and add a new condition to the tree.

    2. Rename the new condition to guide.

    3. Set the Type to Wall.

    4. Select guide as the Boundary Zone.

    5. Select Free slip for Slip Specification.

    6. Select Convection for the thermal condition Option.

    7. Enter 10 W/(m2K) for the Heat Transfer Coefficient.

    8. Enter 323 K for the Convection Temperature.

3.6.3.6. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and other solution properties. Most items indicate that current default values are appropriate.

  1. Review your problem setup.

    1. Select Run Calculation to see available options.

    2. Click Check to review your simulation settings. The workspace will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  2. Calculate a solution.

    1. Click Calculate to start the calculation.

    2. Click the Plots tab to view the plot of the convergence monitors.

    3. Once the calculations are complete, check the solution's listing in the Transcript tab to confirm its success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  3. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering inverse_extrusion as the name of the session.

3.6.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can setup contours, vectors, and so on.

  1. Click the Graphics tab to display the graphics window.

  2. Display an isometric view of the geometry by clicking the   button in the graphics toolbar, then select Isometric from the drop-down.

  3. Display the velocity distribution for the fluid zone.

     ResultsQuick-View Velocity fluid-zone

    Figure 3.48: Contours of Normal Component of Velocity

    Contours of Normal Component of Velocity

  4. Repeat the previous step to create a Pressure contour for the fluid zone using the Quick-View in the ribbon.

    Figure 3.49: Contours of Pressure

    Contours of Pressure

  5. Repeat the previous steps to create a Temperature contour for the fluid zone using the Quick-View in the ribbon.

    Figure 3.50: Contours of Temperature

    Contours of Temperature

    The temperature plot shows that the temperature of the fluid has increased due to viscous dissipation. Additionally, heat exchange with the wall has caused the fluid to cool in the corners of the channel.

  6. Display the geometry from the +Y direction by clicking the   button in the graphics toolbar, then select From +Y Direction from the drop-down.

    Figure 3.51: Contours of Temperature From +Y Direction

    Contours of Temperature From +Y Direction

    From the +Y direction, the deformation of the die and the free jet can be seen.

  7. Display the geometry from the +Z direction by clicking the   button in the graphics toolbar, then select From +Z Direction from the drop-down.

    Figure 3.52: Contours of Temperature From +Z Direction

    Contours of Temperature From +Z Direction

  8. Display contours of velocity magnitude for the die lip.

     ResultsGraphics Contour New

    Figure 3.53: Contour Properties

    Contour Properties

    1. Enter Die_lip for Name.

    2. For the Field, select VELOCITIES[].

    3. For Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Select extrudate-constant from the list and click OK to close the dialog.

      Figure 3.54: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. Click the Display button in the properties panel to view the die lip.

      Figure 3.55: Die Lip Contours of Velocity Magnitude

      Die Lip Contours of Velocity Magnitude

    5. In the properties panel, enable Draw Mesh, then select Mesh Outline from the Overlayed Mesh drop-down.

    6. Click Display.

    7. Use the Graphics window to adjust the view as shown in the following figure.

      Figure 3.56: Contours of Velocity Magnitude

      Contours of Velocity Magnitude

      The plot shows that the left part of the die has been enlarged and the right part has been reduced, relative to the extrudate shape. This is due to the slower fluid velocity in the left portion of the die and higher fluid velocity in the right portion. Additionally, the fluid in the left portion of the die is accelerated leading to the geometry reduction, while the fluid in the right portion of the die is slowed leading to geometry enlargement.

  9. Create a YZ plane for observing the fluid velocity in the left portion of the die.

     Results Surfaces Plane New

    1. Select YZ Plane for Creation Mode.

    2. Enter 0.01 m for X.

  10. Create an additional YZ plane for observing the fluid velocity in the right portion of the die.

    1. Select YZ Plane for Creation Mode.

    2. Enter 0.075 m for X.

  11. Display contours of velocity to observe the fluid velocity in the left portion of the die.

     ResultsGraphics Contour New

    Figure 3.57: Contour Properties

    Contour Properties

    1. Enter velocity-left-block for Name.

    2. For the Field, select VELOCITIES[Z].

    3. For the Surfaces, select plane-1.

    4. Click Display.

    5. View the geometry from the +X direction by clicking the   button in the graphics toolbar, then select From +X Direction from the drop-down.

      Figure 3.58: Contours of Velocity-Z From +X Direction

      Contours of Velocity-Z From +X Direction

  12. Display contours of velocity to observe the fluid velocity in the right portion of the die.

     ResultsGraphics Contour New

    Figure 3.59: Contour Properties

    Contour Properties

    1. Enter velocity-right-block for Name.

    2. For the Field, select VELOCITIES[Z].

    3. For the Surfaces, select plane-2.

    4. Click Display.

    Figure 3.60: Contours of Velocity-Z From +X Direction

    Contours of Velocity-Z From +X Direction

  13. Create a line surface for the left portion of the die.

     ResultsSurfaces Line New

    Figure 3.61: Line Properties

    Line Properties

    1. Set First Points as 0.01, 0.0125, and -0.2 m for X, Y, and Z, respectively.

    2. Set Second Points as 0.01, 0.0125, and 0.35 m for X, Y, and Z, respectively.

  14. Repeat the previous procedure to create a line surface for the right portion of the die, by entering 0.075, 0.005, -0.2, 0.075, 0.005, and 0.35 for First Point X, First Point Y, First Point Z, Second Point X, Second Point Y, and Second Point Z, respectively.

  15. Create an XY plot to observe the fluid acceleration in the left and right portions of the die.

     ResultsPlots XY Plot New

    Figure 3.62: XY Plot Properties

    XY Plot Properties

    1. Select line-1 and line-2 for Lines.

    2. Select Velocities[Z] for the Field.

    3. Set the Plot Direction as 0, 0, and 1 for X Component, Y Component, and Z Component, respectively.

    4. Expand the Curves category and select - - - - from the Pattern drop-down list.

    5. Select * from the Symbol drop-down list, then enter 0.4 for the Size.

    6. Click Plot.

      Figure 3.63: XY Plot

      XY Plot

      From the plot it can be seen that the fluid accelerates in the extrudate along line-1 (corresponding to the left block of the die), while the fluid slows down in the extrudate along line-2 (corresponding to the right block of the die).

3.6.5. Summary

This tutorial introduced the concept of a 3D inverse extrusion problem using the Fluent Materials Processing workspace. You solved the problem using a simple geometry for the extrudate, made suitable assumptions about the physics of the problem, and analyzed the factors affecting the predicted die shape.

3.7. Mesh Superposition Technique

3.7.1. Introduction

This tutorial describes how to use the mesh superposition technique within the Fluent Materials Processing workspace (Fluent Materials Processing Workspace).

In this tutorial you will learn how to use the Fluent Materials Processing workspace to:

  • Easily set up your simulation using the mesh superposition technique.

  • Calculate the solution.

  • Analyze the results and generate a scene of the fluid temperature distribution.

3.7.2. Problem Description

This tutorial demonstrates the use of the mesh superposition technique in the Fluent Materials Processing workspace to analyze the flow of a non-isothermal fluid through a single rotating screw extrusion device.

Figure 3.64: Problem Description

Problem Description

The flow domain consists of an extruder wall, screw, screw shaft, one inlet, and one outlet. These components are shown below in Figure 3.65: Flow Domain Components. For the purpose of this description, a clipped surface of the extruder wall is displayed.

Figure 3.65: Flow Domain Components

Flow Domain Components

The fluid enters into the extruder through an inlet shown below in Figure 3.66: Extruder Boundaries at a volume flow rate of Q = 3 * 10-5 m3 / s, before encountering the solid steel rotating screw. The fluid has the following material properties:

  • Viscosity: Bird-Carreau

    with,

  • Density = 1200 kg / m3

  • Thermal conductivity = 1 W / m / K

  • Heat capacity (Cp) = 1700 J / kg / K

Figure 3.66: Extruder Boundaries

Extruder Boundaries

For simplicity, the screw is defined as an overlapping part by use of the mesh superposition technique. The mesh of the flow domain can be seen below in Figure 3.67: Mesh of the Flow Domain.

Figure 3.67: Mesh of the Flow Domain

Mesh of the Flow Domain

All nodes of the flow domain inside the screw will be constrained and their velocity will be imposed. The density, thermal conductivity, and heat capacity for solid steel will be applied on nodes of the flow domain that overlap with the screw, which is shown below in Figure 3.68: Screw Mesh.

Figure 3.68: Screw Mesh

Screw Mesh

3.7.3. Setup and Solution

3.7.3.1. Preparation

To prepare for running this tutorial:

  1. Download the SSE_thermal.zip file here .

  2. Unzip SSE_thermal.zip to your working directory.

  3. The file, sse_thermal.msh can be found in the folder.

3.7.3.2. Launching Ansys Fluent

  1. Use the to start Ansys Fluent.

  2. In the Fluent Launcher, set the Capability Level to Enterprise.

    The mesh superposition technique (MST) in the Fluent Materials Processing workspace is only available with the Enterprise licensing level.

  3. Select Materials Processing in the list of Fluent workspaces.

  4. Click Start.

3.7.3.3. Setup Your Simulation

Clicking Start in the Fluent Launcher opens the Fluent Materials Processing workspace.

  1. Load the provided mesh file.

    1. In the Setup properties, under Read, click the Mesh... button.

      You can also use the File menu, and choose Read > Polyflow Mesh....

    2. Locate and select the mesh file from your working folder (sse_thermal.msh).

    3. Set the Mesh Length Unit to mm, and click Accept to close the dialog box.

  2. Set up an extrusion simulation (using the Simulation category of the Setup Ribbon).

     Setup Simulation Use Template...

    This will display the Simulation Template dialog.

    1. Set the Type to Compounding / Mixing.

    2. Set the Mode to Extruder.

    3. Set the Number of Moving Parts to 1.

    4. Set the Duration to 3 seconds.

    5. Enable the Thermal option.

    6. Click Apply.

      This will instruct the workspace to set up the appropriate objects and settings based on your selections.


      Note:  Notice the Outline View's use of status icons. A green check mark indicates the properties of that object are satisfactory. A red 'x' indicates that attention is required for that object.


      You can progress down the Outline View (or across the Ribbon) to complete your simulation settings using each object's property pages.

3.7.3.4. General Properties

Click the General node in the Outline View to display general settings for the problem setup. For this tutorial, you can keep the default settings.

3.7.3.5. Material Properties

  1. Edit the density and viscosity settings for the fluid material.

     Setup Materials fluid

    1. Expand the Density Law category.

      1. Enter 1200 kg/m3 for the Density.

    2. Expand the Thermal Properties category.

      1. Under Thermal Conductivity, enter 1 W/(m K) for K.

      2. Under Heat Capacity per Unit Mass, enter 1700 J/(kg K) for Cp.

    3. Expand the Viscosity Law category.

      1. Select Bird-Carreau from the Shear Rate Dependence drop-down list.

      2. Under Bird-Carreau Law, enter 100000 Pa s for Zero Shear Viscosity, 1 s for Time Constant, and 0.7 for the Power Law Index.

  2. Import a new material for the screw.

    Click the Import from library button in the Materials property page. Alternatively, you can right-click on the Materials node, and select Import from library from the context menu.

     Materials  Import from library

    1. Select Miscellaneous for the Family

    2. Select Steel_Solid_300K for the Material.

    3. Click OK to close the dialog.

      The new material will be added to the list of materials in the Outline View.

  3. Remove the original solid material solid.

    1. Select solid in the Outline View, right-click and select Delete from the context menu.

       Setup Materials solid  Delete

      You will be prompted to make sure you want to remove the material. Click Yes to remove the material.

3.7.3.6. Cell Zone Properties

In this portion of the tutorial, you will assign the fluid and moving part cell zones.

  1. Define the cell zone properties for the fluid-zone.

     Setup Cell Zones fluid-zone

    1. For Zones, select the field to open a selection dialog.

    2. In the selection dialog, select fluid and click OK to keep the selection and close the dialog.

    3. Retain the default selection of Generalized newtonian for Fluid Model .

  2. Define the cell zone properties for the moving-part.

     Setup Cell Zones moving-part

    1. For Zones, select screw.

    2. For the Solid Material, select Solid Steel.

    3. Under Angular Velocity, enter -60 rpm for Omega.

3.7.3.7. Boundary Condition Properties

  1. Set the conditions for the inlet boundary.

    1. Select inlet under Fluid Boundary Conditions in the Outline View.

    2. Select Volume flow for the Flow Specification.

    3. For the Method, select Computed dynamically.

      Due to the rotation of the extruder wall, computed dynamically is specified for the inflow condition, which will search the normal force to the inlet that must be set in order to obtain the prescribed volumetric flow rate.

    4. For the Volume Flow Rate, enter 3e-5 m3/s.

    5. Under Thermal Condition, enter 353.15 K for the Temperature.

  2. Set the conditions for the outlet boundary.

    1. Select outlet under Fluid Boundary Conditions in the Outline View.

    2. Select outlet for the Boundary Zone.

  3. Set the conditions for the wall boundary.

    1. Select wall under Fluid Boundary Conditions in the Outline View.

    2. Select wall for the Boundary Zone.

    3. Under Thermal Condition, select Convection for the Option.

    4. Enter 300 W/(m2 K) for the Heat Transfer Coefficient.

    5. Enter 383.15 K for the Convection Temperature.

  4. Create a new boundary condition for the screw shaft.

    1. Click Fluid Boundary Conditions in the outline view, then click New... in the properties window.

    2. Enter shaft for the Name.

    3. From the Type drop-down list, select Wall.

    4. Select shaft as the Boundary Zone.

    5. Select Moving for the Wall Velocity.

    6. Under First Point of Axis, enter 0 m for X1, 0 m for Y1, and 0 m for Z1.

    7. Under Second Point of Axis, enter 0 m for X2, 0 m for Y2, and 1 m for Z2.

    8. Under Angular Velocity, enter -6.2832 rad/s for Omega.

      Although the shaft is rotating at the same angular velocity as the screw, the units of angular velocity for rotating walls is specified in rpm while rad/s is used for moving parts.

    9. Under Translation Velocity, enter 0 m/s for Vx, 0 m/s for Vy, and 0 m/s for Vz.

    10. Under Thermal Condition, select Convection for the Option.

    11. Enter 300 W/(m2 K) for the Heat Transfer Coefficient.

    12. Enter 353.15 K for the Convection Temperature.

3.7.3.8. Solution

Open the Solution branch of the Outline View (or use the Ribbon). Here, you can review problem setup and other solution properties. Most items indicate that current default values are appropriate.

  1. Set the Temperature Initializations.

    1. Click temperature-in-fluid-zone under Temperature Initialization in the Outline View.

    2. Enter 353.15 K for A.

  2. Edit the solution Methods.

     Solution Methods

    1. Enable the Decouple Computation of Temperature option.

    2. Select Hybrid for Integration Rule for Transient Terms.

  3. Edit the Calculation Activities.

     Solution Calculation Activities

    1. Under Launch Options, enter 5 for the Number of Processors.

    2. Under Transient Controls, enter 0.05 s for the Initial Time Step.

  4. Edit the solution Outputs.

     Solution Outputs

    1. Enter 0.05 s for the Time Interval.

  5. Review your problem setup.

     Solution Run Calculation

    1. Click Check to review your simulation settings. The workspace will inform you as to whether or not your setup contains any problems and will provide guidance to solve any issues.

  6. Calculate a solution.

    1. Click Calculate to start the calculation.

    2. Click the Plots tab to view the plot of the convergence monitors.

    3. Once the calculations are complete, check the solution's listing in the Transcript tab to confirm it's success. You can confirm that the solution proceeded as expected by looking for the following printed at the bottom of the listing file:

      The computation succeeded.
  7. Save the work.

    Your work can be preserved using a "session" file - a Materials Processing workspace-specific case file (*.mprcas) that contains your settings and your results.

    File Write Session..., entering sse_thermal as the name of the session

3.7.4. Results

Review results using the various tools in the Results section of the Outline View or the Ribbon where you can create new surfaces, setup contours, and create scenes.

  1. Display the temperature distribution on the inlet, outlet, and extruder wall.

     ResultsGraphics Contour New

    Figure 3.69: Temperature Contour Properties

    Temperature Contour Properties

    1. Enter temperature for the Name.

    2. For the Field, select TEMPERATURE.

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Select inlet, outlet, and wall from the list, then click OK to close the dialog.

      Figure 3.70: Selecting Boundaries in the Surfaces Dialog

      Selecting Boundaries in the Surfaces Dialog

    4. For the contour properties, keep the remaining defaults and click Display to see the temperature contour plot in the graphics window.

      Figure 3.71: Contours of Temperature

      Contours of Temperature

    In the following steps, you will create new surfaces for displaying the temperature field, which will then be used to create a scene.

  2. Create a new iso-surface on the domain.

     Results Sufaces   New Iso-Surface...

    Figure 3.72: Iso-Surface Properties

    Iso-Surface Properties

    1. Enter x=0 for the Name.

    2. For the Field, select Coordinates[X].

    3. For Iso-Value, enter 0 m.

  3. Create a new iso-clip for one half of the extruder wall.

     Results Sufaces   New Iso-Clip...

    Figure 3.73: Iso-Clip x<0 Properties

    Iso-Clip x<0 Properties

    1. Enter x<0 for the Name.

    2. For the Field, select Coordinates[X].

    3. From the Surfaces list, select wall.

    4. Enter -0.021 m for the Minimum.

    5. Enter 0 m for the Maximum.

  4. Repeat the previous step to create another iso-clip for the second half of the wall surface, with the Name of x>0, the Minimum set to 0 m, and the Maximum set to 0.021 m.

    Figure 3.74: Iso-Clip x>0 Properties

    Iso-Clip x>0 Properties

  5. Display the temperature distribution on the second half of the flow domain.

     ResultsGraphics Contour New

    Figure 3.75: Temperature Contour Properties

    Temperature Contour Properties

    1. Enter temperature-x-greater-than-0 for the Name.

    2. For the Field, select TEMPERATURE.

    3. For the Surfaces, select the field to open the Surfaces dialog where you can select the relevant boundaries. Select x>0, x=0, inlet, and outlet from the list, then click OK to close the dialog.

    4. Deselect the Auto-Compute Range option.

    5. Enter 350 K and 440 K for the Minimum Value and Maximum Value, respectively.

    6. Expand the Color Map category.

    7. Enter 100 for the Color Map Size.

    8. From the Type list, select float.

    9. Enter 0 for the Precision.

    10. Disable the Automatically Skip Labels option, then enter 11 for Skip.

    11. Click Display to see the temperature contour plot in the graphics window.

      Figure 3.76: Contours of Temperature for Second Half of Flow Domain

      Contours of Temperature for Second Half of Flow Domain

  6. Repeat the previous step to create an additional contour plot for the other half of the extruder wall, with the Name of temperature-x-less-than-0 and only x<0 selected for Surfaces.

    Figure 3.77: Contours of Temperature for the First Half of the Extruder Wall

    Contours of Temperature for the First Half of the Extruder Wall

  7. Create a new mesh object for the screw.

     Results Meshes   New...

    Figure 3.78: Screw Mesh Properties

    Screw Mesh Properties

    1. Enter screw for the Name.

    2. From the Surfaces list, select screwsurf and shaft.

    3. Disable the Edges option.

  8. Using the newly created graphics objects, create a scene of the temperature field.

     Results Scenes   New...

    Figure 3.79: Scene Properties

    Scene Properties

    1. Deselect any objects that are enabled by default within the list of Graphics Objects.

    2. Under the deformed category in the list of Graphics Objects, enable screw, temperature-x-greater-than-0, and temperature-x-less-than-0.

    3. For temperature-x-less-than-0, set the Transparency to 80.

    4. Click Display to view the scene in the Graphics window.

  9. Play the animation for the new scene.

     ResultsTimestep

    1. To reset the timesteps, change the Current animation timestep value to 1. Alternatively, you can also move the slider all the way to the left, or click the Jump to Start button ( ).

    2. To start the animation, click the Play button ( ).

    3. To stop the animation, click the Stop button ( ).

      Figure 3.80: Contours of Temperature: Timestep = 1

      Contours of Temperature: Timestep = 1

      Figure 3.81: Contours of Temperature: Timestep = 20

      Contours of Temperature: Timestep = 20

      Figure 3.82: Contours of Temperature: Timestep= 40

      Contours of Temperature: Timestep= 40

      Figure 3.83: Contours of Temperature: Timestep = 60

      Contours of Temperature: Timestep = 60

    The results displayed above show that the extruder wall and the screw experience the greatest temperature rise and the temperature within the flow domain increases from inlet to outlet due to the generated heat. The heat is generated almost exclusively within the thin layer adjacent to the barrel wall, and it appears that the screw is the primary source heating the main stream of fluid.

3.7.5. Summary

This tutorial introduced the use of the mesh superposition technique when using the Fluent Materials Processing workspace. You solved the problem using a 3D geometry of a single screw extrusion device and made suitable assumptions about the physics of the problem.