CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Help needed with pimpleDyMFoam case!

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 21, 2017, 16:10
Default Help needed with pimpleDyMFoam case!
  #1
New Member
 
Jeremy
Join Date: May 2016
Posts: 17
Rep Power: 3
j_moulton is on a distinguished road
Hi everyone,

I am attempting to build a case that can model some airflow into an out of an idealized bronchiole-alveolus. Some background: inhalation is the process of your diaphragm contracting, which causes a pressure drop between the alveolar sacs and the atmosphere, creating a rush in of air; exhalation is the same but works backwards (diaphragm expansion causes the alveoli to contract, forcing higher pressure which is equilibrated by air leaving the body).

To do this, I defined my dynamicMeshDict as
Code:
dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver          velocityLaplacian;

velocityLaplacianCoeffs
{
    diffusivity inverseDistance (alv1);
    diffusivity inverseDistance (alv2);
    diffusivity inverseDistance (alv3);
    diffusivity inverseDistance (alv4);
}
The boundary conditions are:

p:

Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
      type         uniformFixedValue;
      uniformValue sine;
      uniformValueCoeffs
      {
          frequency    1;
          amplitude    1;
          scale        1;  // Scale factor for wave
          level        0;  // Offset
      }
    }
    pipeWall
    {
    type           zeroGradient;
    }
    alv1
    {
      type         zeroGradient;
    }
    alv2
    {
      type         zeroGradient;
    }
    alv3
    {
      type         zeroGradient;
    }
    alv4
    {
      type         zeroGradient;
    }
    frontAndBack
    {
        type            empty;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
U:

Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
      type         uniformFixedValue;
      uniformValue sine;
      uniformValueCoeffs
      {
          frequency    1;
          amplitude    1;
          scale        (1 0 0);  // Scale factor for wave
          level        (0 0 0);  // Offset
      }
    }
    pipeWall
    {
        type            zeroGradient;
    }
    alv1
    {
      type              zeroGradient;
    }
    alv2
    {
      type              zeroGradient;
    }
    alv3
    {
      type              zeroGradient;
    }
    alv4
    {
      type              zeroGradient;
    }
    frontAndBack
    {
        type            empty;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
and

pointMotionU:

Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    alv1
    {
      type oscillatingVelocity;
      omega 0.5;
      amplitude (-.25 .25 0);
      value uniform (0 0 0);
    }
    alv2
    {
      type oscillatingVelocity;
      omega 0.5;
      amplitude (.25 .25 0);
      value uniform (0 0 0);
    }
    alv3
    {
      type oscillatingVelocity;
      omega 0.5;
      amplitude (.25 -.25 0);
      value uniform (0 0 0);
    }
    alv4
    {
      type oscillatingVelocity;
      omega 0.5;
      amplitude (-.25 -.25 0);
      value uniform (0 0 0);
    }
    pipeWall
    {
        type            uniformFixedValue;
        uniformValue    constant (0 0 0);
    }
    inlet
    {
        type            uniformFixedValue;
        uniformValue    constant (0 0 0);
    }
    frontAndBack
    {
        type            empty;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
I am stuck at making this work the way I want it to and I need some help and input from those that know OF (much) better than me. There should be a pressure gradient from the expanding walls out to the inlet, and there should be a gradient of velocity entering the inlet. I tried to force the flow by oscillating the velocity and pressure at the inlet, but it's still not coming up right.

Can someone please help me sort this out? Or explain to me where I'm failing. All case files are zipped & attached, a short .avi of the expanding geometry is included for reference.

Case pimpleDyMFoam: https://www.dropbox.com/s/cn1tt7i6ch...MFoam.zip?dl=0

Video of expanding geometry:
https://www.dropbox.com/s/9yeorfujey...nding.avi?dl=0 ( I have a low frequency set here, but you get the idea).
j_moulton is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFD Post - How to check for case in an expression pilakin ANSYS 0 September 26, 2014 04:44
OpenFOAM Case Setup Needed mpl CFD Freelancers 4 October 10, 2013 16:14
Error reading new case montag dp FLUENT 5 September 15, 2011 06:00
Pressure instabilities with interDyMFoam for the floatingObject case nbadano OpenFOAM Running, Solving & CFD 14 March 8, 2011 09:19
MFX FSI for linear case. Help needed. peiji1984 CFX 4 February 22, 2011 13:11


All times are GMT -4. The time now is 12:43.