CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Problems in calculating the fluid traction on the current structure frame in 3D models

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

LinkBack Thread Tools Display Modes
Old   August 6, 2008, 12:04
Default Hi everyone, I am working o
feng wang
Join Date: Mar 2009
Posts: 67
Rep Power: 8
fw407 is on a distinguished road
Hi everyone,

I am working on an FSI project, the solver is based on the icoFsiFoam but the structure code can deal with geometrically nonlinear structure problems and works in the initial structure frame. the solver is working properly in 2-D and I have tried 2-D axially symmetric cases and the results are satisfactory, but I get problems in running a 3-D model say a complete pipe,

I have debuged the solver, and the problem lies in the "setPressure.H" which looks like:


tensorField gradD = fvc::grad(D);

vectorField nc = stressMesh.Sf().boundaryField()[solidPatchID] & inv (I + gradD.T()) * det (I + gradD);

vectorField n = nc / mag(nc);

symmTensorField solidPatchTraction = sigmaFluid.boundaryField()[fluidPatchID] ;

solidPatchTraction *= rhoFluid.value();

tForce.traction() = solidPatchTraction & (n);


for a 3-D model, the equation "vectorField n = nc / mag(nc);" can not be calculated correctly and the computation diverges at the beginning:

snail@snail-laptop:~/OpenFOAM/snail-1.4.1-dev/run/tutorials/icoStVkFsiFoam$ icoStVkFsiFoam . halfVine
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1-dev |
| \ / A nd | Web: |
| \/ M anipulation | |

Exec : icoStVkFsiFoam . halfVine
Date : Aug 06 2008
Time : 16:49:22
Host : snail-laptop
PID : 11528
Root : /home/snail/OpenFOAM/snail-1.4.1-dev/run/tutorials/icoStVkFsiFoam
Case : halfVine
Nprocs : 1
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: laplaceFaceDecomposition
Selecting motion diffusivity: quadratic
Reading transportProperties

Reading field p

Reading field U

Calculating stress field sigmaFluid

Reading/calculating face flux field phi

Reading mechanical properties

Calculating Lame's coefficients

mu = 46153.8 Pa
lambda = 69230.8 Pa
threeK = 300000 Pa
Reading field D

Calculating field V

Calculating stress field gradD

Calculating stress field sigma

Calculating explicit part of div(sigma) Q

Reading coupling properties
the relaxation factor is

Starting time loop

Courant Number mean: 0 max: 0
Time = 1e-07

Setting traction on solid patch
Floating point exception (core dumped)

from the debug results, the "gradD" and "nc" are calculated correctly but the calculation diverges when trying to calculate "vectorField n = nc / mag(nc)". so far,I haven't figured out why, could anyone give me some ideas?

Kind regards

fw407 is offline   Reply With Quote


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
Fluid Structure Interaction using icoFsiFoam Problems lr103476 OpenFOAM Running, Solving & CFD 79 August 7, 2014 09:30
Two-fluid models vs mixture models for bubbly flows Hansong Tang Main CFD Forum 6 December 8, 2009 04:21
Problems calculating field gh with interFoam cricke OpenFOAM Running, Solving & CFD 0 December 10, 2007 08:17
OpenFOAM Debian packaging current status problems and TODOs oseen OpenFOAM Installation 9 August 26, 2007 13:50
More Rotating reference frame problems DaveG CFX 2 October 26, 2005 16:47

All times are GMT -4. The time now is 06:25.