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 Tree54Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   November 20, 2012, 17:22
Default
  #221
Senior Member
 
Kevin Smith
Join Date: Mar 2009
Posts: 103
Rep Power: 8
kev4573 is on a distinguished road
Niels,

I'm personally okay with manually copying the waveProperties.org file, since I already do this sort of thing with the fields in the zero directory. I think your proposed approach makes sense anyway.

Kevin
kev4573 is offline   Reply With Quote

Old   November 20, 2012, 17:50
Default
  #222
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Kevin

Good, I will put it on my waves2Foam-TODO.

With respect the need of two files during the running of setWaveParameters, then the reason is very simple.

Previously, I wrote generated wave parameters of the type Field<Type> to a separate file using the IOField<Type> method in OpenFoam. However, it turned out that it worked perfectly in serial, but for versions >= 2.0, the data files could suddenly not be read in parallel. In order to overcome this problem, those of my colleagues who are running irregular waves, had to run setWaveParameters and manually copy the data into waveProperties. The need of IOField<Type> is due to a restriction in the IOdictionary format, which does not support adding a nonuniform List<Type> to the dict. It was added per say, however, as a raw list without any formatting. For some reason the data was then read in the waveTheory as a bunch of zeros.

As the practical work-around was too cumbersome and error-prone, I am now writing a completely fresh waveProperties file, where all the data are either read from the source file or computed in setWaveProperties. This gives me complete control over the output format, hence it is also possible to have nonuniform fields in the output - so parallel runs can be executed.

The drawback is that the final waveProperties file is stripped of all comments, etc, since each derived class from setWaveProperties know exactly, which bits and piece are needed for say stokesSecond - everything else is disregarded.

I hope that someone liked this long explanation

Have a nice evening,

Niels
kev4573 likes this.
ngj is offline   Reply With Quote

Old   November 21, 2012, 12:09
Default
  #223
New Member
 
niall o sullivan
Join Date: Jan 2011
Posts: 2
Rep Power: 0
niallosully is on a distinguished road
Hi Niels,

I am preforming pimplefoam LES simulations of airflow over research vessels and wanted to implement the wave2foam at the air-sea interface. my question is it possible to apply a uniforminletvelocity at the air section of the wave2foam simulation using this type of simulation with a K-epsilon model in place and would you have any pointers on where i would start changing the waves2foam tutorials to achieve this.

Niall
niallosully is offline   Reply With Quote

Old   November 22, 2012, 05:39
Default
  #224
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Niall

It is possible to add a constant wind velocity in the air in the relaxation zones. Others using waves2Foam have tried this, however, I am unaware of their results.

Please consult the wiki, section 5.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   November 25, 2012, 15:26
Default
  #225
mf.
New Member
 
Marco Fitzner
Join Date: Nov 2012
Posts: 3
Rep Power: 4
mf. is on a distinguished road
Hi Nils,

after a few days of interrupt using OpenFOAM now my waveTheory is running. I am using the OpenFOAM version 2.1 so I had to modify the interFoam solver on my own to get wave2Foam as it's discribed on your wiki. During this the variable DOFVERSION was to set by hand in the Make/options file of the solver. To compile a new waveTheory I had to do this again for the Make/options file of the waveTheories but I didn't. What a bad mistake.

Thank you again for helping me to localize my problem.
Marco
mf. is offline   Reply With Quote

Old   November 26, 2012, 03:11
Default
  #226
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Good morning Marco,

I am a bit puzzled with your post (but congrats that it is now working):

1. waveFoam for 2.1 is supported directly without your need to do anything. I tried to compile it on my own 2.1. distribution, and everything went smoothly. If you do not have the directory waves2Foam/applications/solvers/solvers210, then you should do an update from svn.

2. I do not understand what you mean about setting OFVERSION for waveTheories. If you execute the Allwmake script in either waves2Foam or waves2Foam/src, then the OFVERSION is set automatically; thus no need to manual editing.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   November 28, 2012, 08:22
Default
  #227
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
Hello Niels,

I hope you are doing well. If you might recall my previous posts, I asked you about the possibility of simulating waves in 3D around an offshore structure. To that you replied that it is indeed possible to do so and that your colleagues run such simulations on a regular basis. Now I have been trying to do the same thing for the past one week but the simulations are taking A LOT of time. I wanted to ask you if this is something normal or do I need to make any changes to my case files. My mesh has a total number of 125,526 points and I'm using a time step of 0.01s for an end time equal to 10s. I'm running these simulations on a Intel i5 Quad-Core 3.20 GHz machine with 4 GB RAM.

Thanks,
Sagun
Sagun is offline   Reply With Quote

Old   November 28, 2012, 11:50
Default
  #228
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Well, it sounds like a small problem with merely a 1000 time steps (are you sure it stays there this large?) Furthermore, the mesh sounds like a small of the kind, so I would not be scared at all. For 2D problems some hundred thousands are not uncommon for my simulations.

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   November 28, 2012, 11:56
Default
  #229
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
Well, it has been almost 4 hours now since I started the last simulation and it has only reached 1.16s as of right now. This is the kind of processing speed that I am having to face.
Sagun is offline   Reply With Quote

Old   November 28, 2012, 13:45
Default
  #230
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
That sounds uncommonly slow, but it is hard to say anymore with such limited information.
A colleagues of mine has been running entire irregular sea-states (3 hours), and I have not heard him complain.

/Niels
ngj is offline   Reply With Quote

Old   November 28, 2012, 13:52
Default
  #231
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
Any other information I can provide which can help you suggest something?
Sagun is offline   Reply With Quote

Old   November 28, 2012, 14:22
Default
  #232
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
The full - ready to run - test case would be a nice start.

/ Niels
ngj is offline   Reply With Quote

Old   November 28, 2012, 14:33
Default
  #233
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
These are my entire case files:

http://www.4shared.com/archive/XbKTs...Objecttar.html

I would be extremely grateful if you could take a look and suggest something.

Regards,
Sagun
Sagun is offline   Reply With Quote

Old   November 28, 2012, 14:41
Default
  #234
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Is it correct that I need to create an account to download the files? If yes, please find another way to send the files.

/ Niels
ngj is offline   Reply With Quote

Old   November 28, 2012, 14:53
Default
  #235
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
Ohh I didn't know that. I realized that the file size is only 3.7 MB. So may be I can simply mail it to you?
Sagun is offline   Reply With Quote

Old   November 29, 2012, 03:26
Default 6DoF + waveDyMFoam crashing
  #236
New Member
 
Nima
Join Date: Feb 2012
Location: Perth, Western Australia
Posts: 13
Rep Power: 5
nima3906m is on a distinguished road
Quote:
Originally Posted by dahicke View Post
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
Hi Dave

I just saw your thread today , was wondering if you have solved your problem or not? I have been struggling to find a solution for long time

Regards
Nima
nima3906m is offline   Reply With Quote

Old   November 29, 2012, 05:24
Default
  #237
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Good morning Sagun,

You can find my contact on http://www.dtu.dk/Service/Telefonbog.aspx

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   November 29, 2012, 10:45
Default
  #238
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
Good afternoon Niels,

Please follow this link. And you don't need to create an account on this one.

https://www.yousendit.com/download/W...NEhTRTcwZXNUQw

Thanks once again,
Sagun
Sagun is offline   Reply With Quote

Old   November 29, 2012, 11:44
Default
  #239
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Rotterdam, The Netherlands
Posts: 1,562
Rep Power: 24
ngj will become famous soon enoughngj will become famous soon enough
Hi Sagun

You need to use the correct wave properties. Those defined in your file are merely the standard from the waveFlume tutorial. Please run setWaveParameters. The solution you got is actually the correct one, but based on incorrect input parameters; e.g. the velocities are 200 m/s after the very first time step.

Furthermore, your wave gauges (see controlDict) are partly defined outside the computational domain - leading to segmentation fault after Time = 2 s.

Correcting your setup makes it possible to complete the simulation in the time is takes to go for a cup of coffee

Kind regards,

Niels
ngj is offline   Reply With Quote

Old   November 30, 2012, 10:04
Default
  #240
Member
 
Sagun Tripathi
Join Date: Aug 2012
Location: Amherst, USA
Posts: 78
Rep Power: 4
Sagun is on a distinguished road
Hello Niels,

Thank you so much for taking a look at my case files. Since there was a problem with some of the wave gauges, I removed all of them for the time being. I kept the wave parameters same as those in the waveFlume tutorial since I just wanted to experiment with simple Airy wave theory first. On running setWaveParameters, the simulation does end in about 10 mins but when I try to visualize the flow in paraView there is no flow at all. There seems to be some problem with the generation of waves at the inlet. I have looked hard at the waveProperties file of the 3Dwaves tutorial but I don't seem to be able to understand what I'm doing wrong. I believe I have not understood how to assign wave parameters (esp. direction to waves and relaxation zones) properly. I have copied my waveProperties file here. I'll be grateful if you could take another look.

timeShift 0;

seaLevel 8;

relaxationNames ( inlet outlet );

initializationName outlet;

pName pd;

inletCoeffs
{
waveType stokesFirst;
depth 8;
period 2;
direction ( 1 0 0 );
phi 0;
height 0.2;
waveNumber (1.00608 0 0);
omega 3.14159;

relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;
relaxType INLET;
startX ( -40 8 0 );
endX ( -20 8 3 );
orientation ( 1 0 0 );
}
}

outletCoeffs
{
waveType potentialCurrent;
U ( 0 0 0 );
Tsoft 2;

relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;
relaxType OUTLET;
startX ( 20 8 0 );
endX ( 40 8 3 );
orientation ( 1 0 0 );
}
}

Thanks,
Sagun
Sagun 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 11: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 22:06.