ACU-T: 4000 Transient Dam Break Simulation

Prerequisites

This tutorial provides instructions for setting, solving and viewing results for a transient dam break simulation using the Level Set method. Prior to starting this tutorial, you should have already run through the introductory HyperWorks tutorial, ACU-T: 1000 HyperWorks UI Introduction, and have a basic understanding of HyperMesh, AcuSolve, and HyperView. To run this simulation, you will need access to a licensed version of HyperMesh and AcuSolve.

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

Since the HyperMesh database (.hm file) contains meshed geometry, this tutorial does not include steps related to geometry import and mesh generation.

Problem Description

The problem to be addressed in this tutorial is shown schematically in the figure below. It consists of a square water column held in place by the reservoir walls. At time t=0, the walls are removed and the water column is now free to flow out. The simulation can be used to visualize and study the surge patterns as the column of water rushes out, as in a dam wall break.



Figure 1.

Open the HyperMesh Model Database

  1. Start HyperMesh and load the AcuSolve user profile.
    Refer to the HM introductory tutorial, ACU-T: 1000 HyperWorks UI Introduction, to learn how to select AcuSolve from User Profiles.
  2. Click the Open Model icon located on the standard toolbar.
    The Open Model dialog opens.
  3. Browse to the directory where you saved the model file. Select the HyperMesh file ACU-T4000_DamBreak.hm and click Open.
  4. Click File > Save As.
    The Save Model As dialog opens.
  5. Create a new directory named DamBreak 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 DamBreak as the file name for the database, or choose any name of your preference.
  7. Click Save to create the database.

Set the Transient Simulation Parameters

Set the Analysis Parameters

  1. Go to the Solver Browser, expand 01.Global, then click PROBLEM_DESCRIPTION.
  2. Change the Analysis type to Transient.
  3. Set the Multiphase equation to Level Set.


    Figure 2.

Specify the Solver Settings

  1. In the Solver Browser, click 02.SOLVER_SETTINGS under 01.Global.
  2. Set the Max time steps to 0 and press Enter.
  3. Set the Final time to 1.0 sec.
  4. Set the Initial time increment to 0.002498 sec.
  5. Change the Max stagger iterations to 5 and press Enter.
  6. Change the Relaxation factor to 0.0.
  7. Check that the Flow and Field options are turned On.


    Figure 3.

Define the Nodal Outputs

  1. In the Solver Browser, expand 17.Output then click NODAL_OUTPUT.
  2. In the Entity Editor, set the Time step frequency to 1 and toggle on the Output initial condition field.


    Figure 4.

Create a Multiphase Model and Set the Body Force

Create a Multiphase Material Model

  1. In the Solver Browser, expand 02.Materials then right-click on MULTIPHASE and select Create.
  2. Enter Air-Water as the name of the model.
  3. In the Entity Editor, verify that the Field Interaction Type is set to Level Set.
  4. Select Air_HM as the Field 1 Material.
  5. Enter Air as the Field 1 name.
  6. Select Water_HM as the Field 2 Material.
  7. Enter Water as the Field 2 name.


    Figure 5.

Set the Body Force

  1. In the Solver Browser, expand 03.Body_Force then expand BODY_FORCE and click Gravity_HM.
  2. In the Entity Editor, set the Y gravity to -9.81 m/sec2 and Z gravity to 0.0.


    Figure 6.

Set the Nodal Initial Conditions and Boundary Conditions

In this step, you will start by creating a node-set for assigning nodal initial conditions for volume fraction. All the nodes in the set are assigned a volume fraction of 1 for water, thereby creating a water column.

Create a Node Set

  1. Go to the Model Browser, right-click on empty space in the browser area, and select Create > Set.
  2. In the Entity Editor, rename the set to Water_Column.

Set the Boundary Conditions

By default, all components are assigned to the wall boundary condition. In this step, you will change them to the appropriate boundary conditions and assign material properties to the fluid volumes.
  1. In the Solver Browser, expand 12.Surfaces > WALL.
  2. Click z_pos. In the Entity Editor, change the Type to SYMMETRY.


    Figure 7.
  3. Similarly, change the z_neg component Type to SYMMETRY.
  4. Click sides. In the Entity Editor, verify that the Type is set to WALL.


    Figure 8.
  5. Click Fluid. In the Entity Editor,
    1. Change the Type to MULTIPHASE.
    2. Select Air-Water as the Material.
    3. Select Gravity_HM as the Body Force.


    Figure 9.

Assign Reference Pressure

  1. In the Solver Browser, right-click on 15.Nodal_Boundary_Condition and select Create.
  2. Set the Name to Fixed Pressure Node and change the Definition to Nodes.
  3. Set the Number of Nodes to 1.
  4. Click on the Node collector. Then, in the graphics window, click on a node in the Fluid volume.


    Figure 10.
  5. Click proceed in the panel area.
  6. Change the Boundary condition variable to Pressure and leave the remaining settings as default.


    Figure 11.

Set the Nodal Initial Conditions

  1. In the Solver Browser, expand 01.Global then click 03.NODAL_INITIAL_CONDITION.
  2. In the Entity Editor, under the Field tab, set the Type to Nodal Values.
  3. Set the Multiphase Field to Water.
  4. Set the Select nodes by option to Node Set.
  5. With the Node Set Count set to 1, select Water_Column as the Node Set.
  6. Set the Initial Volume Fraction to 1.0.


    Figure 12.

Assign Nodes to the Node Set

  1. Go to the Model Browser, right-click on empty space in the browser area, and select Create > Block.
  2. In the Entity Editor, rename the block to Water_Column.
  3. Set the Min and Max coordinates to (0, 0, 0) and (0.05715, 0.05715, 0.002748), respectively.
  4. In the Model Browser, click Sets > Water_Column. In the Entity Editor, click on the Nodes collector in the Entity IDs field.
  5. In the panel area, click the nodes collector and select the by block option.
  6. Toggle on the Water_Column block then click select.
    All the nodes in the Water column block are highlighted in the graphics area.


    Figure 13.
  7. Click proceed.
  8. Optional: Turn off the display of the Water_Column block.
  9. Save the model.

Compute the Solution

In this step, you will launch AcuSolve directly from HyperMesh and compute the solution.

Run AcuSolve

  1. Turn on the visibility of all mesh components.
    For the analysis to run, the mesh for all active components must be visible.
  2. Click on the ACU toolbar.
    The Solver job Launcher dialog opens.
  3. Optional: For a faster solution time, set the number of processors to a higher number (4 or 8) based on availability.
  4. The Output time steps can be set to All or Final. Since this is a Transient analysis, set it to All.
  5. Leave the remaining options as default and click Launch to start the solution process.


    Figure 14.

Post-Process the Results with HyperView

Once the solution has converged, close the AcuProbe and AcuTail windows. Go to the HyperMesh window and close the AcuSolve Control tab.

Open HyperView and Load the Model and Results

  1. In the HyperMesh main menu area, click Applications > 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 DamBreak.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 the Water Flow Animation

In this step, you will create an animation of the water flow as it surges once the walls restricting the water column are removed.

  1. Orient the display to the xy-plane by clicking on the Standard Views toolbar.
  2. Click on the Results toolbar to open the Contour panel.
  3. Select Volume_fraction-2-Water (s) as the Result type.
  4. Click Apply to display the volume fraction contour at the first time step.
  5. Click the Legend tab then click Edit Legend.
  6. In the Edit Legend dialog, change the Number of levels to 2 and the Numeric format to Fixed.


    Figure 15.
  7. On the Animation toolbar, click the Animation Controls icon .
  8. Drag the Max frame Rate slider to 50 fps.
  9. Click the Start/Pause Animation icon to play the animation in the graphics area.

Save the Animation

  1. In the menu area, select Preferences > Export Settings > AVI.
  2. In the Export Settings AVI dialog, set the Frame rate to 50 fps and click OK.
  3. On the ImageCapture toolbar, make sure that the Save Image to File option is On.


  4. Click the Capture Graphics Area Video icon .
    The Save Graphics Area Video As dialog opens.
  5. Navigate to the location where you want to save the file, enter a name of your choice, and click Save.

Summary

In this tutorial, you successfully learned how to set up and solve a multiphase flow problem using HyperMesh and AcuSolve. You also learned how to create a multiphase model using the Level Set method. Once the solution was computed, you post-processed the results in HyperView where you generated an animation of the water flow as it surged once the dam walls were removed.