ACU-T: 5200 Rigid-Body Dynamics of a Check Valve

This tutorial provides the instructions for setting up, solving, and viewing results for a simulation of the opening of a pressure check valve. In this simulation, AcuSolve is used to compute the forces on the valve due to the time-varying inlet flow field and to compute the motion of the valve that results from these flow forces. This tutorial is designed to introduce you to a number of modeling concepts necessary to perform simulations of rigid-body dynamics.

Prerequisites

Prior to starting this tutorial, you should have already run through the introductory tutorial, ACU-T: 1000 Basic Flow Set Up, and have a basic understanding of HyperMesh CFD and AcuSolve. 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 to be addressed in this tutorial is shown schematically in Figure 1. It consists of a cylindrical pipe containing water that flows past a check valve with a shutter attached to a virtual spring (not included in the geometry). The inlet pressure varies over time and the movement of the shutter will be determined as a function of the balance of the fluid forces against the reactive force of the spring. The problem is rotationally periodic at 30° increments about the longitudinal axis, and it is assumed that the resulting flow is also rotationally periodic, allowing for modeling with the use of a wedge-shaped section. For this tutorial, a 30° section of the geometry is modeled, as shown in the figure. Modeling a portion of an rotationally periodic geometry leads to reduced computation time while still providing an accurate solution.



Figure 1. Schematic of Check Valve with Spring-Loaded Shutter

The pipe has an inlet diameter of 0.08 m, and is 0.3 m long. The check-valve assembly is 0.085 m downstream of the inlet. It consists of a plate 0.005 m thick with a centered orifice 0.044 m in diameter and a shutter with an initial position 0.005 m from the opening, simulating a nearly closed condition. The shutter plate is 0.05 m in diameter and 0.005 m thick. The shutter plate is attached to a stem 0.03 m long and 0.01 m in diameter. The mass of the shutter and stem is 0.2 kg and its motion is affected by a virtual spring with a stiffness of 2162 N/m. The motion of the valve shutter is limited by a stop mounted on a perforated plate downstream of the shutter.

Note that AcuSolve's internal rigid-body-dynamics solver is not able to simulate contact. Therefore, this problem is formulated to avoid contact between the valve and the stop.

Modeling the geometry as a 30° section requires that the fluid model is set up to be consistent with the rigid-body model. Since only 1/12 of the rigid body is modeled, the forces computed by AcuSolve that act on the valve shutter represent 1/12 of the actual force on the device. Therefore, it is also necessary to account for this in the simulation. There are two methods that can be used to accomplish this:
  1. Scale up the fluid forces calculated by AcuSolve by a factor of 12 to represent the full load on the device when the displacement of the body is computed.

    Using this approach, the full stiffness of the valve spring is used in the rigid-body solution, and the full mass of the valve is used.

  2. Scale down the mass of the valve and the stiffness of the spring to by a factor of 12 to match the fraction of the valve geometry to be modeled.

    Using this approach, the loading passed to the rigid-body solver is not scaled.

This second approach is used in this tutorial; the scaled mass of 0.0167 kg and the scaled stiffness of 180.1667 N/m will be used .

The fluid in this problem is water, which has a density (ρ) of 1000 kg/m3 and a molecular viscosity (μ) of 1 X 10-3 kg/m-sec.

At the start of the simulation the flow field is stationary. Flow is driven by the pressure at the inlet, which varies over time as a piecewise linear function shown in Figure 2. As the pressure at the inlet rises, the flow will accelerate as the valve opens. The turbulence viscosity ratio is assumed to be 10.

The initial inlet pressure is 0 Pa. At 0.002 s, the pressure begins to ramp up, and reaches 26,500 Pa at 0.05 s. The pressure is held at 26,500 Pa and begins to ramp back down starting at 0.2 s, reaching the initial pressure at 0.25 s where it remains for the rest of the simulation.


Figure 2. Transient Pressure at the Inlet

Prior simulations of this geometry indicate that the average velocity at the inlet reaches a maximum of 0.9 m/s. At this velocity, the Reynolds number for the flow is 72,000. When the Reynolds number is above 4,000, it is generally accepted that flow should be modeled as turbulent.

Note that the initial conditions of the flow are actually laminar, however, the increase in flow velocity and flow around the valve shutter is expected to cause a rapid transition to turbulent conditions. Therefore, the simulation will be set up to model transient, turbulent flow. When performing a transient analysis, convergence is achieved at every time step based on the defined stagger criteria. Mesh motion will be modeled using arbitrary mesh movement (arbitrary Lagrangian-Eulerian mesh motion).

For this case, the transient behavior of interest occurs in the time it takes for the pressure to ramp up and ramp back down, which is given by the transient pressure profile. To allow time for the spring to recover, additional time will be simulated. For this tutorial, 0.1 s is added after the pressure drops back to initial conditions, for a total duration of 0.35 s.

Another critical decision in a transient simulation is choosing the time increment. The time increment is the change in time during a given time step of the simulation. It is important to choose a time increment that is short enough to capture the changes in flow properties of interest, but does not require unnecessary computation time.

The change in inlet pressure from initial conditions to maximum occurs over 0.048 s. A time increment of 0.002 s would allow for excellent resolution of the transient changes, without requiring excessive computational time.

Start HyperMesh CFD and Create the HyperMesh Model Database

  1. Start HyperMesh CFD from the Windows Start menu by clicking Start > Altair <version> > HyperMesh CFD.
    When HyperMesh CFD is loaded, the Geometry ribbon is open by default.
  2. Create a new .hm database in one of the following ways:
    • From the menu bar, click File > Save.
    • From the Home tools, Files tool group, click the Save As tool.


      Figure 3.
  3. In the Save File As dialog, navigate to the directory where you would like to save the database.
  4. Enter CheckValve_RBD as the name for the database then click Save.
    This will be your problem directory and all the files related to the simulation will be stored in this location.

Import and Validate the Geometry

Import the Geometry

  1. From the menu bar, click File > Import > Geometry Model.
  2. In the Import File dialog, browse to your working directory then select ACU-T5200_pressureCheckValve.x_t and click Open.
  3. In the Geometry Import Options dialog, leave all the default options unchanged then click Import.


    Figure 4.


    Figure 5.

Validate the Geometry

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


    Figure 6.
    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 7.
  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 Flow

Set Up the Simulation Parameters and Solver Settings

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


    Figure 8.
    The Setup dialog opens.
  2. Under the Physics models setting:
    1. Set Time marching to Transient.
    2. Set the Time step size to 0.002 and the Final time to 0.35.
    3. Select Spalart-Allmaras as the Turbulence model.


    Figure 9.
  3. Click the Solver controls setting and set the Maximum stagger iterations to 3.
    Figure 10.
  4. Close the dialog and save the model.

Create a Multiplier Function for the Inlet Pressure

  1. From the Flow ribbon, click the arrow next to the Setup tool set, then select Multipliers.


    Figure 11.
  2. Click in the Multiplier Library dialog.
  3. In the multiplier creation dialog, change the name of the multiplier function to Inlet Pressure by clicking on the top-left corner.
  4. Set the Type to Piecewise Linear.
  5. Verify that the Variable is set to Time and Evaluation is set to Time Step.
  6. Click four times to add four rows to the bottom of the table.
  7. Enter the table values according to the image below.


    Figure 12.

Assign Material Properties

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


    Figure 13.
  2. Click anywhere on the pipe geometry.
    The entire geometry is highlighted.


    Figure 14.
  3. In the microdialog, select Water from the Material drop-down menu.
  4. On the guide bar, click to execute the command and exit the tool.
  5. Save the model.

Define Flow Boundary Conditions

In this step, you will assign the boundary conditions for the inlet, outlet and the symmetry faces.
  1. From the Flow ribbon, Pressure tool group, click the Stagnation Pressure tool.


    Figure 15.
  2. In the modeling window, click on the inlet surface, as shown in Figure 16.
  3. In the microdialog, enter a value of 1 N/m2 for the Stagnation pressure.
  4. Click the Multiplier function drop-down and select the function Inlet Pressure from the list.


    Figure 16.
  5. Click the Turbulence icon the microdialog, set the Turbulence input type to Viscosity Ratio, and the enter a value of 10 for the Turbulence viscosity ratio.


    Figure 17.
  6. On the guide bar, click to execute the command and remain in the tool.
  7. In the Boundaries legend, right-click on Stagnation pressure and rename it to Inlet then press Enter.
  8. Click the Outlet tool.


    Figure 18.
  9. Click the outlet surface shown in the figure below.


    Figure 19.
  10. On the guide bar, click to execute the command and exit the tool.
  11. Click the Symmetry tool.


    Figure 20.
  12. Select the two surfaces shown in the figure below.


    Figure 21.
  13. On the guide bar, click to execute the command and remain in the tool.
  14. In the Boundaries legend, rename Symmetry to Front_symmetry.
  15. Rotate the model and select the other two symmetry faces.


    Figure 22.
  16. On the guide bar, click to execute the command and remain in the tool.
  17. In the Boundaries legend, rename Symmetry to Back_symmetry.
  18. Save the model.

Set Up Motion

In this step, you will activate the mesh motion and define the rigid body motion for the valve wall. Then, you will define the mesh displacement boundary conditions for the symmetry surfaces.

Define the Mesh Motion Type

  1. From the Motion ribbon, click the Settings tool.


    Figure 23.
  2. In the Setup dialog, change the Mesh motion to Arbitrary.


    Figure 24.

Define Rigid Body Motion

In this step, you will define the rigid body mesh motion on the valve walls.
  1. From the Motion ribbon, click the Rigid Body tool.


    Figure 25.
  2. On the guide bar, activate the Surfaces selector by clicking on it.
  3. For convenience, hide the Front_symmetry and Back_symmetry faces in the modeling window by selecting the four surface and then pressing H on the keyboard or by right-click and selecting Hide from the context menu.


    Figure 26.
  4. Select the valve wall surfaces, shown in the figure below, by using box selection (hold and drag the left mouse button).
    The number of surfaces selected should be 8, which will show on the guide bar.


    Figure 27.
  5. On the guide bar, click on DOF then click anywhere on the valve wall surfaces in the modeling window.
  6. In the microdialog that appears, enter the following values:
    • Mass: 0.0167 kg
    • Center: (0, 0, 0)


    Figure 28.
  7. Click on the X-displacement arrow in the modeling window.
  8. In the microdialog, enter a value of 180.1667 kg/sec2 for Kxx, the stiffness of the spring.


    Figure 29.
  9. Click on the X-rotation arrow and de-activate it.


    Figure 30.
  10. Similarly, click on the remaining four arrows (Y-displacement, Y-rotation, Z-displacement, and Z-rotation) and de-activate them.
    This will ensure that the valve body translates only in the x-direction.
  11. On the guide bar, click to execute the command and exit the tool.
  12. Turn on the display of all the surfaces and save the model.

Define the Mesh Displacement Boundary Conditions

  1. From the Motion ribbon, click the Planar Slip tool.


    Figure 31.
  2. Select the two front symmetry surfaces shown in Figure 32.
  3. In the Mesh Motion legend, rename Planar Slip to Front_symmetry by double-clicking on the name.


    Figure 32.
  4. On the guide bar, click to execute the command and remain in the tool.
  5. Rotate the model and select the two back symmetry faces.
  6. In the Mesh Motion legend, rename Planar Slip to Back_symmetry.
  7. On the guide bar, click to execute the command and exit the tool.
  8. Save the model.

Generate the Mesh

In this step, you will define the mesh controls and then generate the mesh.

Define the Zone Mesh Controls

The zone mesh controls allow you to define a local volume mesh size for specific regions using standard shaped zones like cylindrical, spherical, and cuboidal zones.
  1. From the Mesh ribbon, Zones tool group, click the Cylinder tool.


    Figure 33.
  2. In the modeling window, hover the mouse around the point shown in the figure below. When the preview cylinder zone is parallel to the axis of the pipe, click on the model near the point shown below.
    This point will be the center of the front face of the cylinder. In the next few steps you will edit the co-ordinates of this point manually.


    Figure 34.
  3. Move the mouse cursor away from the center and then click again.
    A preview of the zone should be displayed on the screen along with the manipulator. The manipulator allows you to modify the location and orientation of the zone.


    Figure 35.
  4. Click the center of the manipulator and enter the following coordinates for the center (-0.06, 0.02, 0)


    Figure 36.
  5. Double-click on the cylindrical surface of the zone and enter 0.022 m as the radius of the cylinder.


    Figure 37.
  6. Double-click on the base of the cylinder and enter 0.1 m as the height of the cylinder.


    Figure 38.
  7. Click anywhere in the empty space in the modeling window.
  8. Enter a value of 0.0015 m for the mesh size in the microdialog.
  9. In the Refinement zones legend, rename the zone to Valve region mesh refinement.


    Figure 39.
  10. Save the model.

Define the Boundary Layer Controls

  1. From the Mesh ribbon, click the Boundary Layer tool.


    Figure 40.
  2. Right-click in the modeling window and go to Select > Advanced Select > By Boundaries > Default Wall.
    All the wall surfaces should be highlighted and a dialog for BL specification appears.
  3. Enter the following values in the dialog:
    1. First layer thickness: 0.00035
    2. Total number of layers: 3
    3. Growth method: Constant
    4. Initial growth rate: 1.2
    5. Termination policy: Truncate
    6. Activate Enable surface mesh modification.


    Figure 41.
  4. In the BL Controls legend, rename the BL control to Pipe walls BL.
  5. On the guide bar, click to execute the command and remain in the tool.
  6. In the next few steps, you will define the boundary layer mesh control for the valve walls. For convenience, hide the front and back symmetry faces.
  7. Select the valve wall surfaces using the box selection method.
  8. In the dialog that appears, enter the following values for the BL specification:
    1. First layer thickness: 0.00015
    2. Total number of layers: 3
    3. Growth method: Constant
    4. Initial growth rate: 1.2
    5. Termination policy: Truncate
    6. Activate Enable surface mesh modification.


    Figure 42.
  9. In the BL Controls legend, rename the BL control to Valve walls BL.
  10. On the guide bar, click to execute the command and exit the tool.
  11. Click the drop-down next to the Mesh Controls tool set and select Advanced.


    Figure 43.
  12. In the Advanced mesh settings dialog, change the Boundary layer element type to Tetra then close the dialog.


    Figure 44.
  13. Turn on the display of all the surfaces and save the model.

Generate the Mesh

  1. From the Mesh ribbon, click the Volume tool.


    Figure 45.
  2. In the Meshing Operations dialog, set the Mesh size option to Maximum size and set the Maximum element size to 0.005.
  3. Deactivate Curvature based surface refinement.
  4. Click Mesh to generate the mesh.


    Figure 46.

    Once the meshing process has started, the Run Status dialog appears and the application moves to the Solution ribbon.

  5. Once the meshing is completed, close the Run Status dialog.

Compute the Solution

Define Nodal Outputs

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


    Figure 47.
  2. In the Field Outputs dialog, set the Time step interval to 3 for the Solution variables.


    Figure 48.

Launch AcuSolve

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


    Figure 49.
  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 50.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.

Monitor Results

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


    Figure 51.
  2. Click next to the Library tab.
  3. Keep the X-axis set to Step. For the Y-axis, click the tree browser and select mesh_displacement under Surface Output.
    mesh_x_displacement is selected by default for the coordinate direction.
  4. Check the box for Rigid Body - surface output under Selector.
  5. Check the box for Twin Y-Axis.
  6. Click the tree browser for the second Y-axis and select pressure under Surface Output.
  7. Check the box for Inlet - Output under Selector.


    Figure 52.
  8. Click Create.


    Figure 53.

Post-Process the Results with HM-CFD Post

Create an Animation of Velocity Magnitude

  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. Click the Top face on the View Cube to align the model.


    Figure 54.
  5. In the Post Browser, right-click on Front_symmetry and select Isolate.
  6. Right-click on Front_symmetry once more and select Edit.
  7. In the display properties microdialog, set the display to velocity.
  8. Activate the Legend toggle and set the max velocity to 6.0.
  9. Click and set the Colormap Name to Rainbow Uniform.


    Figure 55.
  10. Click on the guide bar.
  11. Click at the bottom of the modeling window to view a live animation of the flow.


    Figure 56.
  12. Save the animation.
    1. Go to File > Screen Capture > Advanced Capture.
    2. Click on the toolbar.
    3. Uncheck Include mouse cursor.
    4. Set the frame rate to 50.
    5. Click on the toolbar then drag over the area you want to record.
    6. Click to start recording and the same button to stop recording.
    7. Name the file and save it.

Display Pressure Contours and Velocity Vectors on a Section Cut

  1. Right-click in the modeling window and select Show All.
  2. Click the Slice Planes tool.


    Figure 57.
  3. Select the plane perpendicular to the z-axis.
  4. Make sure the View Cube is aligned with the Top face.


    Figure 58.
  5. Click in the slice plane microdialog.
  6. Turn off the Display in the Contour tab of the display properties microdialog.


    Figure 59.
  7. Click the Vector Display tab.
  8. Turn the Display on, set the variable to velocity, and the value to 0.001.
  9. Set the Coloring to pressure.
  10. Activate the Legend toggle and click to reset the range.
  11. Click , set the Legend location to Upper Center, the Legend orientation to Horizontal, and the Colormap Name to Rainbow Uniform.


    Figure 60.
  12. Click on the guide bar.
  13. In the Post Browser, right-click on Slice Plane 1 and select Isolate.


    Figure 61.

Summary

In this tutorial, you learned how to set up and solve a CFD simulation involving rigid body mesh motion. You started by importing geometry, then defined the flow set up. Next, you defined the mesh motion set up and then generated the mesh. Once the solution was computed, you post-processed the results using the Plot tool and HyperMesh CFD Post.