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

Waves2Foam Related Topics

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

Like Tree44Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   September 23, 2012, 14:59
Default
  #141
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Kevin,

Are you sure that the files you have sent me are from a 2.1.0 distribution? The reason I am asking is that they do not compile on a freshly installed 2.1.0 distribution from the Ubuntu repository.

The compile error is:

Code:
Making dependency list for source file waveFoam.C
SOURCE=waveFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam210/src/transportModels -I/opt/openfoam210/src/transportModels/incompressible/lnInclude -I/opt/openfoam210/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam210/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam210/src/finiteVolume/lnInclude -DOFVERSION=21 -I./../../../../src/lnInclude -IlnInclude -I. -I/opt/openfoam210/src/OpenFOAM/lnInclude -I/opt/openfoam210/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/waveFoam.o
In file included from waveFoam.C:63:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:35: error: ‘class Foam::twoPhaseMixture’ has no member named ‘alpha1’
In file included from waveFoam.C:93:
UEqn.H:5: error: ‘class Foam::incompressible::turbulenceModel’ has no member named ‘divDevRhoReff’
/opt/openfoam210/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/waveFoam.o] Error 1
I will have another look at the update to 2.1.0 in waves2Foam one of the following days.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   September 24, 2012, 10:40
Default
  #142
Member
 
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 7
dkingsley is on a distinguished road
Niels,

The version of waves2Foam that Kevin sent to you was compiled for the bleeding edge 2.1.1 git repository.

If you need a 2.1.0 version, he can do that when he gets back from travel next week.

Dennis
dkingsley is offline   Reply With Quote

Old   September 24, 2012, 10:57
Default
  #143
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Dennis

So that would mean that the solvers compile differently from 2.1.0 to 2.1.1? In that case I am very reluctant to distribute a solver in waveFoam, because it will definitely create even more confusion in my perspective. Here the guidelines on the wiki will do better.

Do you agree?

However, I will update the tutorials with the missing files, so there will not be any problems running them.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   September 28, 2012, 07:06
Default
  #144
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 5
JanL is on a distinguished road
Dear All,

I'm still looking for a way to simulate a ship with forward speed in waves. As I'd like to use waveFoam for that purpose, I was wondering if anybody could give me some advice how to realize that type of simulations. I thought about combining potential current with a regular wave theory using combinedWaves. Has anybody experience doing that, or would this actually be a solution to my problem?
Apart from that I would like to vary the incident angle of the waves (e.g. still forward speed of the ship but waves from behind or beam seas) which leads me to the problem, how the INLET and OUTLET relaxationZones could be realised in that case for each waveType.

I'd very much appreciate any advice!

Best regards
Jan
JanL is offline   Reply With Quote

Old   September 28, 2012, 08:19
Default
  #145
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Jan,

Maybe you should consider streamFunction waves instead. The wave-current interaction is a natural output from the estimation of the harmonic coefficients.

What do you mean with INLET/OUTLET. They are only directional flags and has nothing to do with the behaviour of the specific relaxation zone, but can help to easily change the relaxation direction. Try to use relaxationZoneLayout utility and switch between INLET/OUTLET and you will understand the difference.

They are planned to be removed in a future release, as they do cause confusion.

Have a nice weekend,

Niels
ngj is offline   Reply With Quote

Old   September 28, 2012, 08:34
Default
  #146
New Member
 
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 5
JanL is on a distinguished road
Dear Niels,

thank's for your quick answer and your advice. I'll have a look into the streamFunction waveType.

What I meant with INLET/OUTLET, was that if I have a regular wave coming from one direction and a potential current from the opposite (as it would be the case if I'm simulating a ship with forward speed and waves coming from behind, if I'm not wrong!?), than I'd have to apply on both ends an INLET, as well as an OUTLET condition for one respectively the other waveType, right? Or is there another way around?

Have a nice weekend,
Jan
JanL is offline   Reply With Quote

Old   September 28, 2012, 08:56
Default
  #147
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
No, they are only used to change the direction in which the relaxation zone works. As stated, play around with the relaxationZoneLayout utility, and you will understand.

- Niels
ngj is offline   Reply With Quote

Old   September 30, 2012, 23:53
Default waveDyMFoam problem with floating object
  #148
New Member
 
David Hickerson
Join Date: Aug 2009
Location: Monument, CO
Posts: 17
Rep Power: 6
dahicke is on a distinguished road
Send a message via AIM to dahicke
Hi all,

I have been playing with the wave solvers for a while and am running into a problem with a dynamic mesh version using 2.1.x, though I don't think it is version related. I have a 3D wave flume that appears to work OK without floating object (a box). However, when I add the box in one case with the mass calculated to the submerged section, it slowly picks up velocity and sinks reaching over -60 m/sec. In the other case with weight less than calculated, it starts to rise until it reaches an upward velocity of over 60 m/sec. Both of these reach local high velocities that case the Courant number is over 1 and the case crashes.

I have tried turning on momentum prediction, increasing nCorrectors to 3 and nNonOrthogonalCorrectors to 2, which did not seem to change this run away z velocity.

With this new version I had to add a definition for div((muEff*dev(T(grad(U))))) which I set to Gauss linear;

Has anyone encountered this kind of velocity run away?

Dave Hickerson
dahicke is offline   Reply With Quote

Old   October 1, 2012, 02:01
Default
  #149
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Dave,

What happens, if you have the box, but no waves?

BTW: You say it is a dynamic mesh, however, it is pure motion and no refinement, correct?

Kind regards

Niels
ngj is offline   Reply With Quote

Old   October 2, 2012, 02:28
Default
  #150
New Member
 
David Hickerson
Join Date: Aug 2009
Location: Monument, CO
Posts: 17
Rep Power: 6
dahicke is on a distinguished road
Send a message via AIM to dahicke
I ran it without the waves, basically setting the inletCoeffs to long periods and zero height:
inletCoeffs
{
waveType stokesFirst;
period 100000000000;
direction ( 1 0 0 );
//Tsoft 2;
depth 2.0;
omega 1;
phi 0;
waveNumber ( 0.1047 0 0 );
height 0;
relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;
relaxType INLET;
startX ( 0 -1.5 0 );
endX ( 2 1.5 0 );
orientation ( 1 0 0 );
}
}

It appears to work OK. There is a slight oscillation in the z direction as it finds buoyancy center.

Have you noticed problems using the sixDoFRigidBodyDisplacement with the wave BC?

Other than the water height changing local to the rigid body, I don't see how the these two would interact in the code. In my case, the wave has not reached the box yet, so it shouldn't be effected at the time of the crash.

I played with the mass of the body a bit, which affected how fast in the z direction it would get prior to crashing. Closer to the calculated mass the faster the speed, but longer the calculation time. The crash always came around 2.1 seconds simulation time.

Any ideas?

Dave
dahicke is offline   Reply With Quote

Old   October 2, 2012, 03:43
Default
  #151
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Dave,

I am sorry that I cannot help you, as I have never played with the 6DOF in OF - and I was merely asking to find out, if there was a coupling between the waves and the 6DOF, which there appears to be

Good luck,

Niels
ngj is offline   Reply With Quote

Old   October 2, 2012, 19:31
Default
  #152
New Member
 
David Hickerson
Join Date: Aug 2009
Location: Monument, CO
Posts: 17
Rep Power: 6
dahicke is on a distinguished road
Send a message via AIM to dahicke
Niels,

I have a thought. If I run the case on multiple cores with the inlet and outlet sections in different areas than the test section, then I could separate out this behavior.

Is there issues with the wave code on multiple cpu's?

PS: I still want to find the coupling, though.


Dave
dahicke is offline   Reply With Quote

Old   October 2, 2012, 22:57
Default
  #153
Senior Member
 
Kevin Smith
Join Date: Mar 2009
Posts: 103
Rep Power: 7
kev4573 is on a distinguished road
Quote:
Originally Posted by ngj View Post
Hi Kevin,

Are you sure that the files you have sent me are from a 2.1.0 distribution? The reason I am asking is that they do not compile on a freshly installed 2.1.0 distribution from the Ubuntu repository.

The compile error is:

Code:
Making dependency list for source file waveFoam.C
SOURCE=waveFoam.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam210/src/transportModels -I/opt/openfoam210/src/transportModels/incompressible/lnInclude -I/opt/openfoam210/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam210/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam210/src/finiteVolume/lnInclude -DOFVERSION=21 -I./../../../../src/lnInclude -IlnInclude -I. -I/opt/openfoam210/src/OpenFOAM/lnInclude -I/opt/openfoam210/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/waveFoam.o
In file included from waveFoam.C:63:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:35: error: ‘class Foam::twoPhaseMixture’ has no member named ‘alpha1’
In file included from waveFoam.C:93:
UEqn.H:5: error: ‘class Foam::incompressible::turbulenceModel’ has no member named ‘divDevRhoReff’
/opt/openfoam210/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/waveFoam.o] Error 1
I will have another look at the update to 2.1.0 in waves2Foam one of the following days.

Kind regards,

Niels
Hi Niels,

Like Dennis said I am compiling against 2.1.x. I agree it would cause confusion if this is not compiling properly with 2.1.0. A possible solution might be to add the minor version number to the solver directory, something like 'solvers210', and only target stable releases of openfoam (not the git repo version).

Kevin
kev4573 is offline   Reply With Quote

Old   October 3, 2012, 04:03
Default
  #154
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Kevin,

That could be a possible solution. Can I make a version of waveFoam based on my OF-2.1 and send it to you, so you can test if it compiles? If it does, then we will be going with that solution.

The installation is on another computer, so you should have it in a couple of days.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   October 3, 2012, 22:39
Default
  #155
Senior Member
 
Kevin Smith
Join Date: Mar 2009
Posts: 103
Rep Power: 7
kev4573 is on a distinguished road
Niels,

Sure, I'd be happy to test it for you.

Kevin
kev4573 is offline   Reply With Quote

Old   October 6, 2012, 08:17
Default Progressive 2D wave
  #156
New Member
 
ernest
Join Date: Jun 2010
Posts: 21
Rep Power: 5
erncyc is an unknown quantity at this point
Hi Niels,

Thank you for the tremendous job you have done with waves2Foam. I have successfully installed it and am experimenting with it. I was wondering which tutorial addresses a traveling (non-reflected) 2D wave. I tried the standing wave tutorial (wave type ....1st order) but the animation appears to show (as the name suggests) a standing wave. I am a bit confused because in the wiki...I read that there's another wave type called First order standing wave. Or do I need to change the outlet boundary condition? Please help. Thanks.



Ernest
erncyc is offline   Reply With Quote

Old   October 8, 2012, 03:33
Default
  #157
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Good morning,

@Kevin: Please find the waveFoam for 2.1.0 here: http://www.student.dtu.dk/~ngja/waveFoam.tar

@Ernest: Well, if you merely want a propagating wave, then you should look into the "waveFlume" tutorial. The standing wave tutorial shows that you can filter the out-going wave component so well that you get a correct standing wave pattern in the computational domain.
The waveType "standingWave" can be used, if you e.g. want to initialise a standing wave in a closed domain.

Good luck,

Niels
ngj is offline   Reply With Quote

Old   October 8, 2012, 05:54
Default waveFlume
  #158
New Member
 
ernest
Join Date: Jun 2010
Posts: 21
Rep Power: 5
erncyc is an unknown quantity at this point
Thank you Neil for your response. I am now running waveFlume and hopefully results will be positive. I was just wondering if there is some manual I could read to help me understand some of the terms (coefficients) in the waveProperties file. Thank you.

Ernest
erncyc is offline   Reply With Quote

Old   October 8, 2012, 06:33
Default
  #159
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,510
Rep Power: 23
ngj will become famous soon enoughngj will become famous soon enough
Hi Ernest,

You could either read the Wiki or look into the header files of the waveTheories. Most of them have a reference to the literature, which was used to implement the different wave theories. Essentially, the coefficients originate from potential wave theory, hence it is where you can find more information.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   October 8, 2012, 20:31
Default
  #160
Senior Member
 
Kevin Smith
Join Date: Mar 2009
Posts: 103
Rep Power: 7
kev4573 is on a distinguished road
Niels,

Thanks, I'll test this out.

Is there any reason not to make a version for the latest stable release (2.1.1) as well?

Cheers,
Kevin
kev4573 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
Other Topics at OpenFOAM Workshop Milan 2008 hjasak OpenFOAM 2 October 26, 2013 04:33
Sections / Topics in CFD Wiki Roberthealy1 CFD-Wiki 6 August 23, 2007 17:58
CFD Related Educational Programmes Jonas Larsson Main CFD Forum 3 February 9, 2007 10:11
project topics vivekanand CFX 0 October 27, 2004 05:17
Advanced Topics in Aerodynamics Antonio Filippone Main CFD Forum 0 August 28, 1999 12:16


All times are GMT -4. The time now is 21:14.