ACU-T: 2000 Turbulent Flow in a Mixing Elbow

Prerequisites

Prior to starting this tutorial, you should have already run through the introductory tutorial, ACU-T: 1000 UI Introduction. To run this simulation, you will need access to a licensed version of HyperMesh CFD and AcuSolve.

Before you begin, copy the file(s) used in this tutorial to your working directory.

Problem Description

The problem characteristics shown here determine if the flow is laminar or turbulent by calculating the Reynolds number in the pipe. The diameter of the large inlet is 0.1 m, and the inlet velocity (v) is 0.4 m/s. The diameter of the small inlet is 0.025 m, and the inlet velocity is 1.2 m/s.



Figure 1. Schematic of Mixing Elbow

The fluid in this problem is water with the following properties that do not change with temperature: a density (ρ) of 1000 kg/m3 and a molecular viscosity (μ) of 1 X 10-3 kg/m-sec.

Based on mass conservation, the combined flow rate (Q) yields a pipe velocity of 0.475 m/s downstream of the small inlet. The following equations are used to calculate the pipe velocity:(1) Q L arg e I n l e t = A L arg e I n l e t V L arg e I n l e t Q S m a l l I n l e t = A S m a l l I n l e t V S m a l l I n l e t Q T o t a l = Q L arg e I n l e t + Q S m a l l I n l e t V p i p e = Q t o t a l A L arg e I n l e t = 0.475 m / s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacaWGrb WaaSbaaSqaaiaadYeaciGGHbGaaiOCaiaacEgacaWGLbGaamysaiaa d6gacaWGSbGaamyzaiaadshaaeqaaOGaeyypa0JaamyqamaaBaaale aacaWGmbGaciyyaiaackhacaGGNbGaamyzaiaadMeacaWGUbGaamiB aiaadwgacaWG0baabeaakiaadAfadaWgaaWcbaGaamitaiGacggaca GGYbGaai4zaiaadwgacaWGjbGaamOBaiaadYgacaWGLbGaamiDaaqa baaakeaacaWGrbWaaSbaaSqaaiaadofacaWGTbGaamyyaiaadYgaca WGSbGaamysaiaad6gacaWGSbGaamyzaiaadshaaeqaaOGaeyypa0Ja amyqamaaBaaaleaacaWGtbGaamyBaiaadggacaWGSbGaamiBaiaadM eacaWGUbGaamiBaiaadwgacaWG0baabeaakiaadAfadaWgaaWcbaGa am4uaiaad2gacaWGHbGaamiBaiaadYgacaWGjbGaamOBaiaadYgaca WGLbGaamiDaaqabaaakeaacaWGrbWaaSbaaSqaaiaadsfacaWGVbGa amiDaiaadggacaWGSbaabeaakiabg2da9iaadgfadaWgaaWcbaGaam itaiGacggacaGGYbGaai4zaiaadwgacaWGjbGaamOBaiaadYgacaWG LbGaamiDaaqabaGccqGHRaWkcaWGrbWaaSbaaSqaaiaadofacaWGTb GaamyyaiaadYgacaWGSbGaamysaiaad6gacaWGSbGaamyzaiaadsha aeqaaaGcbaaabaGaamOvamaaBaaaleaacaWGWbGaamyAaiaadchaca WGLbaabeaakiabg2da9maalaaabaGaamyuamaaBaaaleaacaWG0bGa am4BaiaadshacaWGHbGaamiBaaqabaaakeaacaWGbbWaaSbaaSqaai aadYeaciGGHbGaaiOCaiaacEgacaWGLbGaamysaiaad6gacaWGSbGa amyzaiaadshaaeqaaaaakiabg2da9iaaicdacaGGUaGaaGinaiaaiE dacaaI1aGaaGjbVlaad2gacaGGVaGaam4Caaaaaa@AF83@

This value is useful in determining the Reynolds number, which in turn can be used to determine if the flow should be modeled as turbulent, or if it should be modeled as laminar.

In order to determine whether the modeled flow would be turbulent or whether it would be laminar, the Reynolds number (Re) should be calculated. The Reynolds number is given by:(2) Re = ρ V D μ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOuaiaacw gacqGH9aqpdaWcaaqaaiabeg8aYjaadAfacaWGebaabaGaeqiVd0ga aaaa@3DE7@

where ρ is the fluid density, V is the fluid velocity, D is the diameter of the flow region, and μ is the molecular viscosity of the fluid. When the Reynolds number is above 4,000, it is generally accepted that flow should be modeled as turbulent.

The Reynolds numbers of 40,000 at the large inlet, 30,000 at the small inlet, and 47,500 for the combined flow indicate that the flow is turbulent throughout the flow domain. The simulation will be set up to model steady state, turbulent flow.

Start HyperMesh CFD and Open the HyperMesh Database

  1. Start HyperMesh CFD from the Windows Start menu by clicking Start > Altair <version> > HyperMesh CFD.
  2. From the Home tools, Files tool group, click the Open Model tool.


    Figure 2.
    The Open File dialog opens.
  3. Browse to the directory where you saved the model file. Select the HyperMesh file ACU-T2000_MixingElbow.hm and click Open.
  4. Click File > Save As.
  5. Create a new directory named MixingElbow_Turbulent and navigate into this directory.
    This will be the working directory and all the files related to the simulation will be stored in this location.
  6. Enter MixingElbow as the file name for the database, or choose any name of your preference.
  7. Click Save to create the database.

Validate the Geometry

  1. From the Geometry ribbon, click the Validate tool.


    Figure 3.
    The Validate tool scans through the entire model, performs checks on the surfaces and solids, and flags any defects in the geometry, such as free edges, closed shells, intersections, duplicates, and slivers.

    The current model does not have any of the issues mentioned above. Alternatively, if any issues are found, they are indicated by the number in the brackets adjacent to the tool name.

    Observe that a blue check mark appears on the top-left corner of the Validate icon. This indicates that the tool found no issues with the geometry model.


    Figure 4.
  2. Press Esc or right-click in the modeling window and swipe the cursor over the green check mark from right to left.
  3. Save the database.

Set Up the Problem

Set Up the Simulation Parameters and Solver Settings

  1. From the Flow ribbon, click the Physics tool.


    Figure 5.
    The Setup dialog opens.
  2. Under the Physics models setting:
    1. Verify that the Incompressible option is selected under Single phase flow.
    2. Set the Time marching to Steady.
    3. Select Spalart-Allmaras as the Turbulence model.
      The robustness and accuracy of the Spalart Allmaras turbulence model makes it an excellent choice for simulation of steady state flows.


    Figure 6.
  3. Click the Solver controls setting and verify that the parameters are set as shown in the figure below.


    Figure 7.
  4. Close the dialog and save the model.

Assign Material Properties

  1. From the Flow ribbon, click the Material tool.


    Figure 8.
  2. Select the model body.
    The entire solid is highlighted.


    Figure 9.
  3. In the microdialog, click the drop-down menu next to Material and select Water.
  4. On the guide bar, click to execute the command and exit the tool.

Assign Flow Boundary Conditions

Set Boundary Conditions for the Large Inlet

The large inlet of the mixing elbow will be defined as a profiled inlet based on average velocity.
  1. From the Flow ribbon, Profiled tool group, click the Profiled Inlet tool.


    Figure 10.
  2. Click the face of the large inlet.


    Figure 11.
  3. In the microdialog, enter a value of 0.4 m/s for Average velocity.


    Figure 12.
  4. Rename the inlet.
    1. From the legend on the left side of the modeling window, double-click on Inlet
    2. Type Large_Inlet and press Enter.
  5. On the guide bar, click to execute the command and remain in the tool.
    Note: The number of inlets created appears in parenthesis on the top-right of the Profiled tool icon.

Set Boundary Conditions for the Small Inlet

The small inlet of the mixing elbow will also be defined using average velocity. The guide bar for the Profiled tool should still be open from the last step.
  1. Click the face of the small inlet.


    Figure 13.
  2. In the microdialog, enter a value of 1.2 m/s for Average velocity.


    Figure 14.
  3. Rename the inlet.
    1. From the legend on the left side of the modeling window, double-click on Inlet
    2. Type Small_Inlet and press Enter.
  4. On the guide bar, click to execute the command and exit the tool.

Set Boundary Conditions for the Outlet

  1. From the Flow ribbon, click the Outlet tool.


    Figure 15.
  2. Click the face of the outlet.


    Figure 16.
  3. In the microdialog, make sure the Static pressure is 0 Pa and the Pressure loss factor is also 0.


    Figure 17.
  4. On the guide bar, click to execute the command and exit the tool.

Set Boundary Conditions for the Symmetry Plane

This geometry is symmetric about the XY midplane, and can therefore be modeled with half of the geometry. In order to take advantage of this, the midplane needs to be identified as a symmetry plane. The symmetry boundary condition enforces constraints such that the flow field from one side of the plane is a mirror image of that on the other side.

  1. From the Flow ribbon, click the Symmetry tool.


    Figure 18.
  2. Select the face of the symmetry plane.


    Figure 19.
  3. In the microdialog, accept the default symmetry conditions.


    Figure 20.
  4. On the guide bar, click to execute the command and exit the tool.
  5. Save the database.

Generate the Mesh

The meshing parameters for this tutorial are already set in the input file.
  1. From the Mesh ribbon, click the Volume tool.


    Figure 21.
    The Meshing Operations dialog opens.
    Note: If the model has not been validated, you are prompted to create the simulation model before running the batch mesh.
  2. Check that the Average element size is 0.025 m.
  3. Accept all other default parameters.


    Figure 22.
  4. Click Mesh.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.
    Tip: Right-click on the mesh job and select View log file to view a summary of the meshing process.
  5. From the back side of the model, observe the refined mesh around the small inlet.


    Figure 23.

Run AcuSolve

  1. From the Solution ribbon, click the Run tool.


    Figure 24.
    The Launch AcuSolve dialog opens.
  2. Set the Parallel processing option to Intel MPI.
  3. Optional: Set the number of processors to 4 or 8 based on availability.
  4. Deactivate the Automatically define pressure reference option.
  5. Leave the remaining options as default and click Run to launch AcuSolve.


    Figure 25.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.
    Tip: While AcuSolve is running, right-click on the AcuSolve job in the Run Status dialog and select View Log File to monitor the solution process.

Post-Process the Results with HM-CFD Post

  1. Once the solution is completed, navigate to the Post ribbon.
  2. From the menu bar, click File > Open > Results.
  3. Select the AcuSolve log file in your problem directory to load the results for post-processing.
    The solid and all the surfaces are loaded in the Post Browser.
  4. Isolate the Symmetry flow boundary in the Post Browser.
  5. Right-click on Symmetry and select Edit.
  6. In the display properties microdialog, set the display to velocity and activate the Legend toggle.
  7. Click and set the Colormap Name to Rainbow Uniform.


    Figure 26.
  8. Click on the guide bar to show the velocity magnitude (m/s) contour.


    Figure 27.
  9. Again, right-click on the Symmetry and select Edit.
  10. Display the pressure (Pa).
  11. Click to reset the range.


    Figure 28.
  12. Click on the guide bar to show the pressure (Pa) contours.


    Figure 29.

Summary

In this tutorial, you worked through a basic workflow to set up a CFD model, carry out a CFD simulation, and post-process the results. You started by importing the model in HyperMesh CFD. Then, you defined the simulation parameters and launched AcuSolve directly from within HyperMesh CFD. Upon completion of the solution by AcuSolve, you used the Post ribbon to create contour plots.