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

dynamic mesh problem with Oscillating B.C.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 25, 2014, 06:47
Question dynamic mesh problem with Oscillating B.C.
  #1
New Member
 
h. lu
Join Date: Oct 2009
Location: TX, USA
Posts: 1
Rep Power: 0
Haier is on a distinguished road
Dear Foamers,

I have been struggling with a problem (similar to wave tank/basin) related to dynamic mesh where I am trying to apply an oscillating boundary. The following are the basic settings:

  • OF v2.3.0 on Suse 13.1
  • solver: interDyMFoam

In constant/dynamicMeshDict
Code:
dynamicFvMesh       dynamicMotionSolverFvMesh;
motionSolverLibs    ("libfvMotionSolvers.so");
solver              displacementLaplacian;
displacementLaplacianCoeffs
{
    diffusivity  inverseDistance (movingBody);
}
In 0/pointDisplacement
Code:
boundaryField
{
...
    sidePlane
    {   
        type            slip;
    }
...
    movingBody
    {
        type            oscillatingDisplacement;
        amplitude    (0 0 1.0);
        omega        1.0;
        value           uniform (0 0 0);
    }
}
also, in system/fvSolution.solvers:
Code:
    cellDisplacement
    {   
        solver          GAMG;
        tolerance       1e-8;
        relTol          0;  
        smoother        GaussSeidel;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10; 
        agglomerator    faceAreaPair;
        mergeLevels     1;  
    }   
    cellDisplacementFinal
    {   
        $cellDisplacement;
        tolerance       1e-8;
        relTol          0;  
    }
Here is part of the error message I got:
Code:
[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2   in "/lib64/libc.so.6"
[0] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:? 
[0] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? 
[0] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[0] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[0] #8  
[0]  at ??:?
[0] #9  __libc_start_main in "/lib64/libc.so.6"
[0] #10  
[0]  at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125
The error message seems to indicate something went wrong while using GAMG solver. When I visualized the mesh movement, weird "ripple" near moving boundary were observed on the "sidePlane" which is supposed to be flat, see attachment. "Ripple" occurred around t=0.699 and the simulation crashed shortly after (around t=0.793).

Right now I am doubting if the "slip" type BC is the cause. Anyone has experience to share or bother to give any comments?

Thanks,
H.
Attached Images
File Type: png mesh_at_0.002400.png (46.8 KB, 16 views)
File Type: png mesh_at_0.606726.png (47.6 KB, 16 views)
File Type: png mesh_at_0.699011.png (54.9 KB, 15 views)
File Type: png mesh_at_0.792755.png (69.9 KB, 16 views)
Haier is offline   Reply With Quote

Reply

Tags
dynamic mesh, oscillatingdisplacement

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
dynamic mesh for multi-region problem alundilong OpenFOAM Programming & Development 4 October 22, 2015 19:43
Update of the variables after dynamic mesh motion. gtg258f OpenFOAM Programming & Development 9 January 18, 2014 11:08
Dynamic mesh simulation (pimpleDyMFoam) problem sidlof OpenFOAM Running, Solving & CFD 6 April 12, 2013 11:24
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
Dynamic mesh problem Karthick FLUENT 1 October 11, 2004 05:47


All times are GMT -4. The time now is 02:55.