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

"solidDisplacementFoam" solver usage

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 12, 2012, 19:46
Default "solidDisplacementFoam" solver usage
  #1
New Member
 
Vlad Miryaha
Join Date: Nov 2012
Posts: 2
Rep Power: 0
Vlad Miryaha is on a distinguished road
Hello,
in the paper http://powerlab.fsb.hr/ped/kturbo/op...tresspaper.pdf
it`s written, that
"The initial condition consists of specified distribution of u and du/dt at time zero."
I think that the paper authors use "solidDisplacementFoam" solver.

Does anybody know how to set initial du/dt distribution???

TIA
Vlad Miryaha is offline   Reply With Quote

Old   November 12, 2012, 20:43
Default
  #2
Senior Member
 
Hisham's Avatar
 
Hisham El Safti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 247
Blog Entries: 10
Rep Power: 8
Hisham is on a distinguished road
Hi Vlad,

The solver in question is to my understanding solidDisplacementFoam as you mentioned. Currently, I think it is not possible to introduce initial velocity. I am not aware of the initial version where that was possible. Nevertheless, you can introduce a new volVectorField called initialVelocity that is read from a dictionary and add a term to the equilibrium equation that is fvc::ddt(initialVelocity) ... I do not know if this is a proper way for your application but if you can explain more about your problem maybe you can get better help!

Best regards,
Hisham
Hisham is offline   Reply With Quote

Old   November 13, 2012, 07:21
Default
  #3
New Member
 
Vlad Miryaha
Join Date: Nov 2012
Posts: 2
Rep Power: 0
Vlad Miryaha is on a distinguished road
Thank you for your answer, Hisham!
At this moment I`m doing some test computations of transient wave propagation in a cube (1m*1m*1m).
At time zero initial displacement at one side is set. (I want to see a compression plane wave propagetion and a reflection from opposite side).
Wave propagates and it`s reflected, but it becomes smearing.
Wave becomes smearing even before it reachs opposite side.
I had rectangle mesh 40*40*40.
Boundary conditions are set to zeroGradient (I tried symmetryPlane too).
The value of "relTol" and "nCorrectors" in fvSolution file are set to 0.01 and 10 respectively.
Approximation orders in fvScheme file are set to cubic/fourth for spatial derivatives.

I can`t understand the reason of such significant smearing.

I hope that if I set initial velocity instead of displacement result becomes better.
Vlad Miryaha is offline   Reply With Quote

Old   November 13, 2012, 09:21
Default
  #4
Senior Member
 
Hisham's Avatar
 
Hisham El Safti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 247
Blog Entries: 10
Rep Power: 8
Hisham is on a distinguished road
Dear Vlad

I do not think you need to change the way of input. What I think could be the problem is the mesh size ... if the mesh is relatively coarse, it could cause dampening of waves. Try to increase the density of your mesh to the point where it does not significantly affect wave reflection as you expect it.

Best regards
Hisham
Hisham is offline   Reply With Quote

Reply

Tags
initial condition, soliddisplacmentfoam

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
Unexplained Error during Solver Runs cfb CFX 6 November 9, 2012 16:42
Strange residuals of the Density Based Solver Pat84 FLUENT 0 October 22, 2012 15:59
Suggestion for the usage of appropriate solver gandesk OpenFOAM 5 November 10, 2010 13:00
OpenFOAM Solver/BC usage description murrayjc OpenFOAM 3 August 25, 2009 04:48
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


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