|
[Sponsors] |
November 12, 2012, 19:46 |
"solidDisplacementFoam" solver usage
|
#1 |
New Member
Vlad Miryaha
Join Date: Nov 2012
Posts: 2
Rep Power: 0 |
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 |
|
November 12, 2012, 20:43 |
|
#2 |
Senior Member
Hisham Elsafti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 257
Blog Entries: 10
Rep Power: 17 |
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 |
|
November 13, 2012, 07:21 |
|
#3 |
New Member
Vlad Miryaha
Join Date: Nov 2012
Posts: 2
Rep Power: 0 |
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. |
|
November 13, 2012, 09:21 |
|
#4 |
Senior Member
Hisham Elsafti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 257
Blog Entries: 10
Rep Power: 17 |
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 |
|
Tags |
initial condition, soliddisplacmentfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 16: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 05:48 |
why the solver reject it? Anyone with experience? | bearcat | CFX | 6 | April 28, 2008 15:08 |