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 Tree90Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 18, 2014, 21:09
Default
  #801
New Member
 
yalong cai
Join Date: Feb 2014
Location: New York
Posts: 13
Rep Power: 3
pizicai is on a distinguished road
Hello, Niels,

I have some confusions about the term 'initializationName' suggested in waveProperties.input.
this term appears twice in openfoamWiki
1. The entry initializationName in waveProperties defines which sub-dictionary is used for initialising the wave field.
2. //Used by setWaveField and it is the sub-dictionary,
// which defines which wave theory to be used for initialization


what confuses me is why the term 'initializationName' is set to 'outlet' in each 'waveProperties.input' file rather than 'inlet'. actually, after I change it to 'inlet'
error comes out:

FOAM Warning :
From function setWaveField::setWaveField( const fvMesh& mesh, volVectorField& U, volScalarField& alpha, volScalarField& p)
in file preProcessing/setWaveField/setWaveField.C at line 77

The specified value of Tsoft is non-zero in the waveType: `stokesFirst'
specified in the sub-dictionary waveProperties::inletCoeffs

Consequently, the initialised `wave field' is set to a horizontal free surface with zero velocity.

thanks for your attention.
pizicai is offline   Reply With Quote

Old   June 20, 2014, 12:45
Default porouswavefoam/gambit mesh
  #802
New Member
 
Join Date: Apr 2014
Posts: 6
Rep Power: 3
da14 is on a distinguished road
Hii

Im try to simulate a porous breakwater with wave generation on waves2foam using porousWaveFoam with Gambi mesh.

My problem is that when i run fluentMeshToFoam to convert the msh file the porosity zone is not defined at cellZones. So when i run the solver the porosity information is not found so i get error.

Can you help me?

Regards

Diogo Afonso
da14 is offline   Reply With Quote

Old   June 21, 2014, 05:57
Default
  #803
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,660
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Good morning,

@Yalong: The program does the right thing. The choice of using the outlet for initialisation is merely random. When Tsoft > 0 it does not matter. If is only if you want a non-flat initial field that you have to think about the zone to use and the magnitude of Tsoft.

@DA14: This is a pre-processing/meshing question, so please refer it to an appropriate thread.

Kind regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is online now   Reply With Quote

Old   June 27, 2014, 05:46
Default $wave_dir
  #804
New Member
 
Dmitrijs Gavrilovs-Stepanovs
Join Date: Jun 2014
Posts: 15
Rep Power: 3
Dmitrjs is on a distinguished road
Hi everyone,
I am new to Linux and OpenFOAM, so it is extremely painful for me to configure it correctly. I am trying to install waves2Foam toolbox at the moment, I followed the instructions, but when I try to execute Allwmake file, it gives me an error message:
~~~~~~~~
FATAL ERROR.
The directory path WAVES_DIR=:/home/waves2Foam does not exist.
Correct the path in bin/bashrc

This means that you should correct the environmental variable $WAVES_DIR
such that it is pointed (in absolute sense) to the path, where waves2Foam
is located.

Please note that once bin/bashrc is created, bin/bashrc.org is only an
inactive file. The latter is also the only of the two files, which is
updated through the SVN-repository.

EXITING
~~~~~~~~~~~~
I know that I have to change the path in the bashrc for the WAVES_DIR variable, my question is what is the right path if my waves2Foam folder is just in my home directory? I have tried /home/waves2Foam , but it did not work.

Thank you very much in advance!
Dmitrjs is offline   Reply With Quote

Old   June 27, 2014, 06:36
Default
  #805
New Member
 
Dmitrijs Gavrilovs-Stepanovs
Join Date: Jun 2014
Posts: 15
Rep Power: 3
Dmitrjs is on a distinguished road
Hi again!
Please ignore the previous post since I figured it out!
Dmitrjs is offline   Reply With Quote

Old   June 27, 2014, 08:40
Default
  #806
Senior Member
 
Paritosh Vasava
Join Date: Oct 2012
Location: Lappeenranta, Finland
Posts: 553
Rep Power: 13
vasava will become famous soon enough
Dmitrijs, could you please explain how you did it?
vasava is offline   Reply With Quote

Old   June 27, 2014, 08:53
Default
  #807
New Member
 
Dmitrijs Gavrilovs-Stepanovs
Join Date: Jun 2014
Posts: 15
Rep Power: 3
Dmitrjs is on a distinguished road
Quote:
Originally Posted by vasava View Post
Dmitrijs, could you please explain how you did it?
OK, so you just go to the wave2Foam directory which should be in your HOME directory. Then you go to /bin and select bashrc file. In that file find label ### USER DEFINED ENVIRONMENTAL VARIABLES and find 'export WAVES_DIR=<some location>' just under it. This is the path you need to correct. I changed it to:
export WAVES_DIR=$HOME/waves2Foam
and it worked.
Hope it will work for you too!
Dmitrjs is offline   Reply With Quote

Old   June 27, 2014, 10:35
Default Tutorials
  #808
New Member
 
Dmitrijs Gavrilovs-Stepanovs
Join Date: Jun 2014
Posts: 15
Rep Power: 3
Dmitrjs is on a distinguished road
Hi guys!
It might be a silly question, but are there any instructions or additional explanations for the 6 tutorial cases that come with waves2Foam?
Thank you in advance!
Dmitrjs is offline   Reply With Quote

Old   June 30, 2014, 08:22
Default
  #809
New Member
 
andrew odonnell
Join Date: Jun 2014
Posts: 8
Rep Power: 3
andr5w5 is on a distinguished road
Hi there,

I am very new to OpenFoam and Linux and I have trying to install the waves2Foam toolbox. I am using OpenFoam 2.3.0 and so I have followed the steps on the wiki page to copy the "interFoam" folder into the waves2foam solvers folder and to change the waveFoam.C, createFields.H, the make/files and the make/options files. I have also changed the Path in the bashrc file to the correct location of the toolbox.

However when I try to run the waveFlume tutorial via "./Allrun" I get the following error message in the log.postProcessWaves2Foam file:

" /opt/openfoam230/bin/tools/RunFunctions: line 52: postProcessWaves2Foam: command not found "

I also get similar error message in the other log files except from the blockMesh file!

Can anyone tell me where I have went wrong?

Many thanks in advance.
andr5w5 is offline   Reply With Quote

Old   June 30, 2014, 16:36
Default
  #810
New Member
 
yalong cai
Join Date: Feb 2014
Location: New York
Posts: 13
Rep Power: 3
pizicai is on a distinguished road
hello, Niels,

I want to model a solitary wave impacts on a fixed plate above SWL. in this model, the wave is generated by Wave2Foam and the mesh is created by snappyHexMesh, that means the plate is just digged out from a plate (this is a 2D model), whose surface is defined as walls. the model could be seen in photo named "model description"

my problem is: before the wave comes through the plate, the water close to the plate produce velocity, you can see the details in the photo named "velocity_graph_before_wavecomes", as a rusult, the plate carry some force from water before wave comes, which could been seen in the photo named "plate_genwave_z"

i am suspicious that there is some problems with "waveProperties.input". you can find it in the zip-file named "plate_solitary_model"

I wonder if it is obligatory to set outlet BC as potential wave when I just want the water flow out smoothly, could I just set a relaxationzone without setting potential wave? what is interesting is water in the whole domain rotate with a small velocity from start, which i assume is caused by potential wave?

(all files could be found in dropbox, link has been attached)
https://www.dropbox.com/sh/duwtl5pym...JzR62FjTZJFB4a
pizicai is offline   Reply With Quote

Old   June 30, 2014, 17:07
Default
  #811
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,660
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Good evening,

@Dmitrjs: To the best of my knowledge, there is only the documentation on the wiki, except for the porous dambreak, see our paper Jensen et. al (2014), and the Beji Battjes case, see e.g. my thesis or the recent publication by Cadelho.

@Andrew: Either your compilation was unsuccessful for some reason, or the solvers and utilities are not in your search path. You do not need to make waveFoam yourself, as it is already there for 2.3.0. Furthermore, note that the guidelines are outdated, as is also explicitly stated on the wiki, but they do show the general idea.

@Pizicai: I am unfortunately not in a position, where I can help you with resolving your problems. It sounds like a classical generation of energy on refined meshes due to non-balanced momentum equation in interFoam.

Kind regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is online now   Reply With Quote

Old   June 30, 2014, 20:12
Default
  #812
New Member
 
andrew odonnell
Join Date: Jun 2014
Posts: 8
Rep Power: 3
andr5w5 is on a distinguished road
Hi Niels,

So the interFoam folder does not need to be copied over into a new solvers "solvers21" folder ?

I will try to play around with this more this morning and hopefully have more success than Ive had the last few days.

Thanks for your reply.
andr5w5 is offline   Reply With Quote

Old   July 1, 2014, 03:27
Default Wavemaker simulation
  #813
New Member
 
David AKNIN
Join Date: Jun 2014
Posts: 2
Rep Power: 0
davidicl is on a distinguished road
Hello everyone,

I'm trying to simulate a wavemaker with the interFoam solver (air phase + water phase) considering a laminar flow.
I have implemented my own input kinematics in the waves2foam package based on the uniform velocity profile of a piston-type wavemaker.
The velocity BC at the wavemaker is a no-slip.
My mesh is a 3m long wave flume, with a water depth of 0.6m.
The resolution is 5cm, except around the free surface where I had to refine the mesh to capture the free surface position accurately (resolution of 2mm).
I also had to refine the mesh at the wavemaker to account for the local high velocity gradients due to the no-slip boundary condition an the paddle (resolution of 5mm).

My simulation always crashes at the same time (t = 0.40s).
I didn't use adaptative time-stepping, and reducing the time step doesn't avoid the crash.
Also, the courant number doesn't go up before it crashes.

My boundary conditions are:

U:

inlet
{
type waveVelocity;
refValue uniform ( 0 0 0 );
refGradient uniform ( 0 0 0 );
valueFraction uniform 1;
value uniform ( 0 0 0 );
}
bottom
{
type fixedValue;
value uniform ( 0 0 0 );
}
outlet
{
type fixedValue;
value uniform ( 0 0 0 );
}
atmosphere
{
type pressureInletOutletVelocity;
value uniform ( 0 0 0 );
}
front
{
type slip;
}
back
{
type slip;
}


p_rgh:

front
{
type zeroGradient;
}
back
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type totalPressure;
rho none;
psi none;
gamma 1;
p0 uniform 0;
value uniform 0;
}


alpha1:

inlet
{
type waveAlpha;
refValue uniform 0;
refGrad uniform 0;
valueFraction uniform 1;
value uniform 0;
}
bottom
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
front
{
type zeroGradient;
}
back
{
type zeroGradient;
}


Here are the last 3 iterations before it crashes, could anyone help ?


MULES: Solving for alpha1
Phase-1 volume fraction = 0.884502 Min(alpha1) = -0.000918702 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884501 Min(alpha1) = -0.000975817 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884501 Min(alpha1) = -0.00102995 Max(alpha1) = 1
Relaxing time: 0.04 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000459341, Final residual = 9.65422e-15, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000712516, Final residual = 2.7008e-14, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.00157154, Final residual = 3.37197e-14, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00131087, Final residual = 8.49745e-13, No Iterations 26
GAMG: Solving for p_rgh, Initial residual = 5.6643e-05, Final residual = 8.0606e-13, No Iterations 19
GAMG: Solving for p_rgh, Initial residual = 5.58024e-06, Final residual = 9.06201e-13, No Iterations 15
GAMG: Solving for p_rgh, Initial residual = 3.8308e-07, Final residual = 3.82669e-13, No Iterations 12
time step continuity errors : sum local = 2.65472e-16, global = -1.89997e-17, cumulative = -1.20753e-15
GAMG: Solving for p_rgh, Initial residual = 2.30848e-05, Final residual = 9.23408e-13, No Iterations 18
GAMG: Solving for p_rgh, Initial residual = 2.16917e-06, Final residual = 6.50879e-13, No Iterations 13
GAMG: Solving for p_rgh, Initial residual = 2.36853e-07, Final residual = 9.31239e-13, No Iterations 11
GAMG: Solving for p_rgh, Initial residual = 1.9956e-08, Final residual = 7.22653e-13, No Iterations 8
time step continuity errors : sum local = 5.01315e-16, global = -7.7352e-17, cumulative = -1.28489e-15
GAMG: Solving for p_rgh, Initial residual = 8.72431e-07, Final residual = 9.6311e-13, No Iterations 12
GAMG: Solving for p_rgh, Initial residual = 5.31826e-08, Final residual = 8.48754e-13, No Iterations 9
GAMG: Solving for p_rgh, Initial residual = 5.95305e-09, Final residual = 7.37499e-13, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 5.83901e-10, Final residual = 6.19992e-13, No Iterations 4
time step continuity errors : sum local = 4.301e-16, global = -2.95754e-17, cumulative = -1.31446e-15
ExecutionTime = 905.95 s ClockTime = 911 s

Courant Number mean: 0.00122961 max: 0.218001
Interface Courant Number mean: 8.42295e-05 max: 0.113623
Time = 0.4055

MULES: Solving for alpha1
Phase-1 volume fraction = 0.884501 Min(alpha1) = -0.00108301 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.8845 Min(alpha1) = -0.00113319 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.8845 Min(alpha1) = -0.00118055 Max(alpha1) = 1
Relaxing time: 0.04 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000458296, Final residual = 1.93453e-14, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000708513, Final residual = 3.09945e-14, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.000680929, Final residual = 2.04159e-13, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00145775, Final residual = 7.57413e-13, No Iterations 105
GAMG: Solving for p_rgh, Initial residual = 0.000112549, Final residual = 7.91032e-13, No Iterations 86
GAMG: Solving for p_rgh, Initial residual = 1.2736e-05, Final residual = 8.64368e-13, No Iterations 87
GAMG: Solving for p_rgh, Initial residual = 1.70598e-06, Final residual = 8.5032e-13, No Iterations 73
time step continuity errors : sum local = 5.88191e-16, global = -9.4146e-17, cumulative = -1.40861e-15
GAMG: Solving for p_rgh, Initial residual = 4.04471e-05, Final residual = 7.85174e-13, No Iterations 86
GAMG: Solving for p_rgh, Initial residual = 6.2229e-06, Final residual = 8.98764e-13, No Iterations 82
GAMG: Solving for p_rgh, Initial residual = 7.91941e-07, Final residual = 9.64218e-13, No Iterations 70
GAMG: Solving for p_rgh, Initial residual = 1.15067e-07, Final residual = 8.55931e-13, No Iterations 60
time step continuity errors : sum local = 5.92062e-16, global = 8.55776e-17, cumulative = -1.32303e-15
GAMG: Solving for p_rgh, Initial residual = 6.83839e-07, Final residual = 8.45255e-13, No Iterations 70
GAMG: Solving for p_rgh, Initial residual = 7.13413e-08, Final residual = 9.04201e-13, No Iterations 52
GAMG: Solving for p_rgh, Initial residual = 3.75992e-09, Final residual = 8.70047e-13, No Iterations 23
GAMG: Solving for p_rgh, Initial residual = 5.88567e-10, Final residual = 8.6016e-13, No Iterations 29
time step continuity errors : sum local = 5.94978e-16, global = -4.85994e-17, cumulative = -1.37163e-15
ExecutionTime = 908.83 s ClockTime = 914 s

Courant Number mean: 0.00122992 max: 0.21805
Interface Courant Number mean: 8.40171e-05 max: 0.113785
Time = 0.406

MULES: Solving for alpha1
Phase-1 volume fraction = 0.884499 Min(alpha1) = -0.00123552 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884499 Min(alpha1) = -0.00128781 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884499 Min(alpha1) = -0.00133746 Max(alpha1) = 1
Relaxing time: 0.04 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000459398, Final residual = 1.06474e-14, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000707475, Final residual = 2.65673e-14, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.00125732, Final residual = 1.15904e-13, No Iterations 4
[3] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/david/OpenFOAM/david-2.2.1/platforms/linux64GccDPOpt/bin/waveFoam"
[3] #8
[3] in "/home/david/OpenFOAM/david-2.2.1/platforms/linux64GccDPOpt/bin/waveFoam"
[3] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #10
[3] in "/home/david/OpenFOAM/david-2.2.1/platforms/linux64GccDPOpt/bin/waveFoam"

Thanks !
David
davidicl is offline   Reply With Quote

Old   July 1, 2014, 08:36
Default
  #814
New Member
 
andrew odonnell
Join Date: Jun 2014
Posts: 8
Rep Power: 3
andr5w5 is on a distinguished road
Sorry to bother you again with what I'm sure is a silly mistake on my part but I still cant seem to get the compilation to work without errors.

It seems the problem lies in the "COMPILE LIBRARY" part of the code where I continue to get an error message stating:

" /bin/sh: 1: g+: not found
make: *** [/home/andrew/OpenFOAM/andrew-2.3.0/platforms/linuxGccDPOpt/lib/libwaves2Foam.so] Error 127"

I have made sure that g++ is installed, is there an error that the program is looking for "g+" and not "g++" ?

Thanks.
andr5w5 is offline   Reply With Quote

Old   July 1, 2014, 16:52
Default
  #815
New Member
 
yalong cai
Join Date: Feb 2014
Location: New York
Posts: 13
Rep Power: 3
pizicai is on a distinguished road
Hello, Niels,

Thanks for your prompt reply. Actually, I found the mistake afterward. This is supposed be a 2D model, however, after snappyHexMesh, the mesh around the plate has more than one grids, after I run extrudeMesh, the chaos around plate before wave comes disappears. sorry to ask this question without a careful check.

By the way, my last question still exists, I wonder if it is obligatory to set outlet BC as potential wave when I just want the water flow out smoothly, could I just set a relaxationzone without setting potential wave? I tried delete 'wave type' part and kept
relaxationzone definition part of the "outletcoefficients section', and the result is it cannot run.
pizicai is offline   Reply With Quote

Old   July 5, 2014, 09:47
Default
  #816
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,660
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Good afternoon,

@Andrew: I have no idea, why you are getting that error. Best of luck.

@Pizicai: The use of relaxation zones is the way to have waves smoothly leaving the domain in waves2Foam. There are other methods in the literature, but they are not implemented in waves2Foam. You cannot define a relaxation zone without defining all the required properties.

Kind regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is online now   Reply With Quote

Old   July 9, 2014, 02:50
Default
  #817
Member
 
maryam morta
Join Date: Sep 2013
Posts: 54
Rep Power: 4
mary mor is on a distinguished road
Quote:
Originally Posted by mary mor View Post
Hi,
I used waves2Foam solver successfully and got the surfaceElevations right for waveflume case.
Then I tried to compile waveDyMFoam which worked successfully. In order to do this I recompiled waves2Foam and other softwares. Then I tried to compile the waveFlume case again. But the error below is shown in log.postProcessWaves2Foam file:

Code:
Create time


Reading g
Processing: regularSpectrum
    Processing sub-action: interpolateSurfaceElevation
        - Interpolation range: [10; 0.969063]


--> FOAM FATAL ERROR: 
bad size -903

    From function List<T>::List(const label size, const T&)
    in file lnInclude/List.C at line 69.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Foam::List<double>::List(int, double const&) in "/opt/openfoam210/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#3  Foam::postProcessingWaves::equidistantTime(Foam::List<std::pair<double, int> > const&, Foam::dictionary const&) in "/home/maryam/OpenFOAM/maryam-2.1.0/platforms/linuxGccDPOpt/lib/libwaves2FoamProcessing.so"
#4  Foam::interpolateSurfaceElevation::evaluate() in "/home/maryam/OpenFOAM/maryam-2.1.0/platforms/linuxGccDPOpt/lib/libwaves2FoamProcessing.so"
#5  
 in "/home/maryam/OpenFOAM/maryam-2.1.0/platforms/linuxGccDPOpt/bin/postProcessWaves2Foam"
#6  __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7  
 in "/home/maryam/OpenFOAM/maryam-2.1.0/platforms/linuxGccDPOpt/bin/postProcessWaves2Foam"

How can I fix this? I think there were some problems with the recomplition of the solver.
Besides is there any problem with the working of waveDyMFoam and postProcessWaves2Foam?
Thanks in advanced
Hi
Regarding to my previous post, I found out the problem was that endTime in controlDict was less than the tMin in postProcessing file.

Now I've got another question. The distribution of gauges in probeDefinition is as lineDistribution with a constant distance between them.
How can I add gauges in the specified points in the tank. I've got a floating object in the domain, so I can not specify them as lineDistribution, because some of them will be in the field of floating object.
I saw there was a pointDistributions.C file in the source folder. Is it possible to use this kind of distribution.

I would be thankful Niel if you answer me soon because I'm in a hurry
Thanks so much for your attention
mary mor is offline   Reply With Quote

Old   July 9, 2014, 08:06
Default
  #818
New Member
 
Dmitrijs Gavrilovs-Stepanovs
Join Date: Jun 2014
Posts: 15
Rep Power: 3
Dmitrjs is on a distinguished road
Good afternoon everyone,
Hopefully you all are having a good summer and I was wondering if you could help me with the following problem.
I am trying to modify the waveFlume tutorial by changing the geometry (increasing the domain) and wave inlet conditions. Everything runs OK, but when I view it in paraFoam, the wave flume stops at certain point in the tank as if there is a relaxation zone. If I change the period to larger value, then the wavelength becomes longer and waves reach the relaxation zone I set up at the outlet, however, I need to use a lower value of period. Here is my waveProperties.input file:

seaLevel 0.00;

// A list of the relaxation zones in the simulation. The parameters are given
// in <name>Coeffs below.
relaxationNames (inlet outlet);

initializationName outlet;

inletCoeffs
{
// Wave type to be used at boundary "inlet" and in relaxation zone "inlet"
waveType stokesFirst;

// Ramp time of 2 s
Tsoft 2;

// Water depth at the boundary and in the relaxation zone
depth 1.0;

// Wave period
period 1.0;

// Phase shift in the wave
phi 0.000000;

// Wave number vector, k.
direction (1.0 0.0 0.0);

// Wave height
height 0.033;

// Specifications on the relaxation zone shape and relaxation scheme
relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;

relaxType INLET;
startX (0 0.0 -1);
endX (3 0.0 1);
orientation (1.0 0.0 0.0);
}
};

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

relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;

relaxType OUTLET;
startX (14.408 0.0 -1);
endX (17.965 0.0 1);
orientation (1.0 0.0 0.0);
}
};

I do not know whether there is something wrong with the relaxation zones (the dimensions of them) or the value of wave height and period are too small for this domain, so I thought someone could help me to understand what is wrong.
Thank you in advance!
Dmitrjs is offline   Reply With Quote

Old   July 9, 2014, 16:53
Default
  #819
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,660
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Good evening,

@Maryam: If you make a spelling error in lineDistribution you are given all the possible options when the pre-processing program fails - otherwise it should be straight forward for you to code your own stuff.

@Dimitrijs: Please read this thread carefully: How to give enough info to get help (Hint: You have by far given enough information for anyone to help you).

Kind regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is online now   Reply With Quote

Old   July 17, 2014, 01:22
Default
  #820
New Member
 
Betsy
Join Date: Jul 2014
Location: Honolulu, HI
Posts: 5
Rep Power: 3
betsybrite is on a distinguished road
Hi Niels,
First of all, I want to echo many others in this post by saying that I am so grateful for your waves2Foam toolbox and your attention to this thread! I am working on adding the waves2Foam toolbox for use with compressibleInterFoam instead of interFoam. I have been able to generate waves with the same amplitude as in interFoam. The problem I am running into is that in order to use the compressibility factor psi in compressibleInterFoam I need to be solving for absolute pressure instead of gauge pressure. It appears in the compressible tutorials, this is accomplished by setting:

defaultFieldValues
(
volScalarFieldValue alpha1 1
volScalarFieldValue p_rgh 1e5
);

in setFieldsDict. Is there a different way to do this through waves2Foam since we don't use setFieldsDict?

Thank you so much!
Betsy
betsybrite 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 17:54.