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 HyperWorks CFD, AcuSolve, and HyperView. To run this simulation, you will need access to a licensed version of HyperWorks CFD and AcuSolve.

Prior to running through this tutorial, copy HyperWorksCFD_tutorial_inputs.zip from <Altair_installation_directory>\hwcfdsolvers\acusolve\win64\model_files\tutorials\AcuSolve to a local directory. Extract ACU-T5200_pressureCheckValve.x_t from HyperWorksCFD_tutorial_inputs.zip.

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 HyperWorks CFD and Create the HyperMesh Model Database

  1. Start HyperWorks CFD from the Windows Start menu by clicking Start > Altair <version> > HyperWorks CFD.
    When HyperWorks 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 doesn’t 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 select Exit from the context menu to exit the tool.
  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. Click the Time setting.
  3. Set the radio button to Transient, the Time step size to 0.002 sec, and the Final time to 0.35 sec.


    Figure 9.
  4. Click the Flow setting.
  5. Set the radio button to Turbulent and set the Turbulence model to Spalart-Allmaras..


    Figure 10.
  6. Click the Solver controls setting and set the Maximum stagger iterations to 3.


    Figure 11.
  7. Exit the Setup 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 12.
  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 13.

Assign Material Properties

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


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


    Figure 15.
  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 16.
  2. In the modeling window, click on the inlet surface, as shown in Figure 17.
  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 17.
  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 18.
  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 19.
  9. Click the outlet surface shown in the figure below.


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


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


    Figure 22.
  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 23.
  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 24.
  2. In the Setup dialog, change the Mesh motion to Arbitrary.


    Figure 25.

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 26.
  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 27.
  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 28.
  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 29.
  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 30.
  9. Click on the X-rotation arrow and de-activate it.


    Figure 31.
  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 32.
  2. Select the two front symmetry surfaces shown in Figure 33.
  3. In the Mesh Motion legend, rename Planar Slip to Front_symmetry by double-clicking on the name.


    Figure 33.
  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 34.
  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 35.
  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 36.
  4. Click the center of the manipulator and enter the following coordinates for the center (-0.06, 0.02, 0)


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


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


    Figure 39.
  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 40.
  10. Save the model.

Define the Boundary Layer Controls

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


    Figure 41.
  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 42.
  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 43.
  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 44.
  12. In the Advanced mesh settings dialog, change the Boundary layer element type to Tetra then close the dialog.


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

Generate the Mesh

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


    Figure 46.
  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 47.

    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 48.
  2. In the Field Outputs dialog, set the Time step interval to 3 for the Solution variables.


    Figure 49.

Launch AcuSolve

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


    Figure 50.
  2. In the Launch AcuSolve dialog, enter the following text in the Additional arguments box: -tlog -lprobe.
    This will launch the AcuProbe and AcuTail utilities automatically once the run is launched.
  3. Set the Parallel processing option to Intel MPI.
  4. Optional: Set the number of processors to 4 or 8 based on availability.
  5. Deactivate the Automatically define pressure reference option.
  6. Leave the remaining options as default and click Run to launch AcuSolve.


    Figure 51.
    The Run Status dialog opens. Once the run is complete, the status is updated and you can close the dialog.

Monitor the Solution with AcuProbe

While AcuSolve is running, you can monitor flow characteristics such as inlet pressure, displacement of the valve, and velocity of the valve using AcuProbe.

  1. Go to the AcuProbe window. In the Data Tree, expand Surface Output > Rigid Body – surface output > Geometry.
  2. Right-click on mesh_x_displacement and select Plot.
  3. Similarly, expand Inlet - Output > Pressure under Surface Output.
  4. Right-click on total_pressure and select Plot.
    Note: You might need to click on the toolbar in order to properly display the plot.


    Figure 52.

Post-Process the Results with HyperView

In this step, you will create an animation of the valve motion as the water flows across the valve.

Open HyperView and Load the Model and Results

  1. Start HyperView from the Windows Start menu by clicking Start > All Programs > Altair <version> > HyperView.
    Once the HyperView window is loaded, the Load model and results panel should be open by default. If you do not see the panel, click File > Open > Model.
  2. In the Load model and results panel, click next to Load model.
  3. In the Load Model File dialog, navigate to your working directory and select the AcuSolve .Log file for the solution run that you want to post-process. In this example, the file to be selected is CheckValve_RBD.1.Log.
  4. Click Open.
  5. Click Apply in the panel area to load the model and results.
    The model is colored by geometry after loading.

Create an Animation of Velocity Magnitude

  1. In Results Browser, expand the list of Components and then click on the Isolate Shown icon .
  2. Click the Front_symmetry - Output component to turn off the visibility of all the components except the front symmetry surface.


    Figure 53.
  3. Orient the display to the xy-plane by clicking on the Standard Views toolbar.
  4. Click on the Results toolbar to open the Contour panel.
  5. In the panel area, set the Result type to Velocity (v) and verify that the drop-down below is set to Mag (magnitude).
  6. Click the Components entity selector.


    Figure 54.
  7. In the Extended Entity Selection dialog, select Displayed.
  8. Click Apply in the panel area.
  9. Go to the Legend tab and then click Edit Legend.
  10. In the Edit Legend dialog, change the Type to Dynamic scale and the Numeric format to Fixed, then click OK.
  11. Click on the Animation toolbar to play the velocity magnitude animation on the front symmetry plane.


    Figure 55.
  12. On the ImageCapture toolbar, click on the Capture Graphics Area Video icon .
  13. In the Save Graphics Area Video As dialog, browse to the directory where you want to save the animation and give a name to the video, ex: velocity mag, then click Save.

Display Pressure Contours and Velocity Vectors on a Section Cut

In the next step, you will create a section cut on the mid-z plane and then display the pressure contours and velocity vectors on that cross section.

  1. In the Results Browser, turn off the display of all the Components except the Fluid component (Fluid_Main_SolidBody_2_1).


    Figure 56.
  2. Click on the Results toolbar to open the Contour panel.
  3. In the panel area, set the Result type to Pressure (s).
  4. Click the Components entity selector. In the Extended Entity Selection dialog, select Displayed.
  5. Click Apply to create the contour plot of pressure.
  6. In the panel area, under the Result tab, activate the Overlay result display check box (if not set already).


    Figure 57.
  7. Click the Section cut icon on the HV-Display toolbar.
  8. In the panel area, click Add to create a new section cut named Section 1.
  9. In the Define plane section, set the axis to Z Axis and then click Apply.
  10. Set the Z Base coordinate to 0 and press Enter.
  11. Change the Display options from Clipping plane to Cross section.
  12. Click Gridline. In the Gridline Options dialog, deactivate the Show check box under Grid line then click OK.
  13. Click the Vector icon on the Results toolbar to open the Vector panel.
  14. In the panel area, set the Result type to Velocity (v).
  15. Click the Selection drop-down and select Sections from the list of options.


    Figure 58.
  16. Click the Sections entity selector then select All.
  17. In the panel area, activate the Overlay result display check box (if not set already).
  18. Click Apply.
  19. Under the Plot tab, verify that only the X+Y+Z Resultant option is selected.
  20. Go to the Display tab, set the Size scaling option to Uniform, and enter a value of 0.002 in the size field.
  21. Set the Color by option to Direction and set the X+Y+Z color to White.


    Figure 59.
  22. Go to the Section tab, activate the Projected check box, then click Apply.
    The vector plot should look like the one shown in the figure below. The result at 0.156 sec is shown.


    Figure 60.

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 AcuProbe and HyperView.