ACU-T: 6500 Flow Through Porous Medium

Prerequisites

This tutorial provides the instructions for setting up, solving and viewing results for a simulation of a flow through porous medium. 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 and AcuSolve. To run this simulation, you will need access to a licensed version of HyperWorks 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 the figure below. It consists of a cylindrical channel with a porous medium in the flow section. As the flow passes through this section, a pressure drop is observed. In this simulation, an inlet velocity will be assigned to the flow and pressure drop across the porous medium will be calculated. The length of the porous section is 0.06 m and the fluid is air with a density of 1.225 kg/m3 and a molecular viscosity of 1.781e-5 kg/m-s. The inlet velocity of the flow is 0.2 m/s.



Figure 1.

In AcuSolve, a porous media problem can be solved using one of two methods: the superficial velocity method and physical velocity method. In this tutorial, we adopt the superficial velocity method, which provides good representation of the pressure drop through a porous media. However, it cannot predict the velocity increase in the porous region since the velocity in the porous region remains the same as those outside of the porous zones. The more accurate representation of velocity inside the porous region can be obtained using the physical velocity method, which solves the continuity and momentum equation inside the porous media using the intrinsic averaging method.

The pressure loss across the porous region using the superficial velocity method is modelled as:

ΔP Δ x i = C D μ k V i + C p ρ k V m V i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq qHuoarcaWGqbaabaGaeuiLdqKaamiEamaaBaaaleaacaWGPbaabeaa aaGccqGH9aqpcaWGdbWaaSbaaSqaaiaadseaaeqaaOWaaSaaaeaacq aH8oqBaeaacaWGRbaaaiaadAfadaWgaaWcbaGaamyAaaqabaGccqGH RaWkcaWGdbWaaSbaaSqaaiaadchaaeqaaOWaaSaaaeaacqaHbpGCae aadaGcaaqaaiaadUgaaSqabaaaaOGaamOvamaaBaaaleaacaWGTbaa beaakiaadAfadaWgaaWcbaGaamyAaaqabaaaaa@4CFC@

where Δ P MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iuaaaa@3831@ is the pressure drop (Pa) across the porous region, Δ x i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iEamaaBaaaleaacaWGPbaabeaaaaa@3973@ is the porous region length (m), k MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Aaaaa@36E6@ is the permeability, C D MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGebaabeaaaaa@37B3@ is the Darcy coefficient, C p MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGWbaabeaaaaa@37DF@ is the Forchheimer coefficient, μ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiVd0gaaa@37AC@ is the viscosity (kg/m-sec), ρ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdihaaa@37B6@ is density (kg/m3), V m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOvamaaBa aaleaacaWGTbaabeaaaaa@37EF@ is velocity magnitude (m/sec), and V i MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOvamaaBa aaleaacaWGPbaabeaaaaa@37EB@ is velocity (m/sec).

In this tutorial, you will utilize the porous coefficient calculator to calculate the Darcy coefficient and the Forchheimer coefficient after curve fitting of the array data of pressure drop vs velocity. The pressure drop across the porous region can be simplified as the permeability is not included in the curve fitting.

ΔP Δx = C D μ V x + C p ρ V x 2 B V x +A V x 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq qHuoarcaWGqbaabaGaeuiLdqKaamiEaaaacqGH9aqpcaWGdbWaaSba aSqaaiaadseaaeqaaOGaeqiVd0MaamOvamaaBaaaleaacaWG4baabe aakiabgUcaRiaadoeadaWgaaWcbaGaamiCaaqabaGccqaHbpGCcaWG wbWaa0baaSqaaiaadIhaaeaacaaIYaaaaOGaeSipIOJaamOqaiaadA fadaWgaaWcbaGaamiEaaqabaGccqGHRaWkcaWGbbGaamOvamaaDaaa leaacaWG4baabaGaaGOmaaaaaaa@50FC@

The figure below shows the Darcy and Forchheimer coefficient calculator after providing the velocity (m/s) and pressure drop (Pa) data in addition to the thickness (m) of the porous medium zone, density (kg/m3) and viscosity (kg/m-s) of air.


Figure 2.

The Darcy coefficient C D MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGebaabeaaaaa@37B3@ and Forchheimer coefficient C p MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGWbaabeaaaaa@37DF@ are then calculated based on the curve fit coefficients

Start HyperWorks CFD and Open the HyperMesh Database

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


    Figure 3.
    The Open File dialog opens.
  3. Browse to the directory where you saved the model file. Select the HyperMesh file ACU-T6500_PorousMedia.hm and click Open.
  4. Click File > Save As.
  5. Create a new directory named PorousMedia 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 PorousMedia as the file name for the database, or choose any name of your preference.
  7. Click Save to create the database.

Validate the Geometry

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.

To focus on the physics part of the simulation, this tutorial input file contains geometry which has already been validated. Observe that a blue check mark appears on the top-left corner of the Validate icon on the Geometry ribbon. This indicates that the geometry is valid, and you can go to the flow set up.


Figure 4.

Set Up Flow

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. Set Time marching to Steady.
    2. Select Laminar as the Turbulence model.


    Figure 6.
  3. Click the Solver controls setting.
  4. Confirm that the Steady update factor and the Steady maximum steps are set to 0.6 and 100, respectively.


    Figure 7.

Assign Material Properties

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


    Figure 8.
  2. Verify that the material Air is assigned to the model's three solids.


    Figure 9.

Define the Porous Medium

  1. From the Flow ribbon, Porous tool group, click the Cartesian Porous Media tool.


    Figure 10.
  2. Select the middle solid on the model.


    Figure 11.
  3. On the guide bar, click Orientation.
  4. Left-click to place a point anywhere on the selected solid.
  5. In the microdialog, click the coefficient calculator and enter the following values for the coefficient calculator.
    • porous zone thickness = 0.06 (m)
    • air density = 1.225 kg/m3
    • air viscosity = 1.781e-05 kg/m-sec


    Figure 12.
  6. Close the coefficient calculator and confirm the values of the Darcy coefficient and Forchheimer coefficient. Enter permeability values for the three directions as shown in the figure below


    Figure 13.
  7. In the microdialog, click to open the Orient tool then verify that the direction is aligned to the global x-axis.


    Figure 14.
  8. On the guide bar, click to execute the command and exit the tool.

Assign the Flow Boundary Conditions

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


    Figure 15.
  2. Select the inlet face.


    Figure 16.
  3. In the microdialog, set the velocity parameters as shown below.


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


    Figure 18.
  6. Select the outlet face.


    Figure 19.
  7. Accept the default parameters then click on the guide bar.

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 20.
    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.01.
  3. Accept all other default parameters.


    Figure 21.
  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.

Run AcuSolve

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


    Figure 22.
  2. Set the Parallel processing option to Intel MPI.
  3. Optional: Set the number of processors to 4 or 8 based on availability.
  4. Leave the remaining options as default and click Run to launch AcuSolve.


    Figure 23.

Post-Process with the Plot Tool

  1. In the Run Status dialog, right-click on AcuSolve run and select Plot time history to launch the Plot Manager.


    Figure 24.
  2. Click to create a new plot.
  3. Click under the Model panel to create a user defined function.
  4. In the Data Tree, expand Integrated Surface Output > Pressure then select pressure.
  5. In the Name field, enter dP.
  6. In the Value field, type value=.
  7. Under Parts, choose AUTO Porous-1 SolidBody_2_1 internal then click Insert Field to insert the field as part of the value, as shown below.


    Figure 25.
  8. Type - at the end of the user defined function value.
  9. Under Parts, choose AUTO Porous-3 SolidBody_4_3 internal then click Insert Field to insert the field as part of the value.
  10. Click to add the user defined function.


    Figure 26.
  11. Click the x axis to switch to Time Steps.


    Figure 27.

    The pressure drop across the porous inlet (AUTO Porous-1 SolidBody_2_1 internal) and porous outlet (AUTO Porous-3 SolidBody_4_3 internal) surfaces is 0.0223 Pa.

Summary

In this tutorial, you learned how to set up and solve a flow simulation with porous medium. You started by importing the HyperWorks CFD input database and then you defined the porous medium. You also learned how to compute the Darcy and Forchheimer coefficients using the calculator available in HyperWorks CFD. Next, you assigned the flow boundary conditions and generated the mesh. Once the solution was computed, you created a plot of the pressure drop across the porous section using the Plot tool.