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 Tree76Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2014, 14:17
Default
  #701
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Nick,

Please do read the error message carefully. I cannot say it more clearly myself.
If still in doubt, please address the posts above (e.g. use google), where the same question has been asked and resovled a number of times.

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 offline   Reply With Quote

Old   January 30, 2014, 19:34
Default
  #702
New Member
 
Nick Krgs
Join Date: Jan 2014
Posts: 11
Rep Power: 3
Nick_civ is on a distinguished road
Hi Niels,

There was some confusion about directories. I finally compiled it succesfully.
Thank you for your time.

NK
Nick_civ is offline   Reply With Quote

Old   February 6, 2014, 10:41
Default potentialCurrent
  #703
New Member
 
nikita tryaskin
Join Date: Oct 2009
Location: Helsinki/St. Petersburg
Posts: 22
Rep Power: 7
eternityboy is on a distinguished road
hello, i would like to set potentialCurrent on both sides. from the inlet side and the outlet

but when I do so, i got strange freeSurface behavior after the relaxation zone

inlet2Coeffs
{
waveType potentialCurrent;
U (2 0 0);
Tsoft 0;
relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;

relaxType INLET;
startX (-8 2.7 0.1868);
endX (-1.5 -2.7 0.1868);
orientation (1.0 0.0 0.0);
}
};

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

relaxationZone
{
relaxationScheme Spatial;
relaxationShape Rectangular;
beachType Empty;

relaxType OUTLET;
startX (8 2.7 0.1868);
endX (10 -2.7 0.1868);
orientation (1.0 0.0 0.0);
}
};
eternityboy is offline   Reply With Quote

Old   February 6, 2014, 11:55
Default
  #704
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
It is very difficult to help, when you are not stating what type of weird behaviour you see. But, I will make a guess that one problem could be that you impose 2.0 m/s from one side and have 0.0 m/s on the other.

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 offline   Reply With Quote

Old   February 8, 2014, 15:14
Default
  #705
New Member
 
behrang chenari
Join Date: Nov 2013
Location: Coimbra, Portugal
Posts: 11
Rep Power: 3
Behrang is on a distinguished road
Send a message via Skype™ to Behrang
Hi Neils,

In waveFlume tutorial, I modified the wave theory in waveProperties file, in order to simulate wave within the wave flume based on Stokes second and fifth order wave theory. By visualizing the results, I can see that waves are generated, but How can I support it? (I mean, is there any experimental data available to compare the results and validate them?)
Please advise.

Best,
Behrang
Behrang is offline   Reply With Quote

Old   February 9, 2014, 15:21
Default
  #706
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Good evening Behrang,

The literature is very rich on analytical expressions and experimental (field and laboratory) data, so I would recommend you to find data sets, which are closely related to your topic at hand.

Another source of information could be any of the articles, which have already been applying waves2Foam: http://openfoamwiki.net/index.php/Co...in_These_Works

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 offline   Reply With Quote

Old   February 11, 2014, 10:06
Default New wave theory
  #707
New Member
 
kian
Join Date: Jan 2014
Posts: 3
Rep Power: 3
kiang is on a distinguished road
Hi Dear Niels,

Thanks a lot for your waveFoam solver. I am going to make a new waveTheory but I did not find a tutorial about it. Would you guide me how can I make a new wave based on existing wave theories in waves2Foam?

Regards,

Kiang
kiang is offline   Reply With Quote

Old   February 11, 2014, 10:47
Default
  #708
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hallo Kiang,

Have a look in the folder:

Code:
waves2Foam/doc/templateWaveTheory
If you execute the script "makeNewWaveTheory", it should create the necessary folder for you, so all you need to modify is the algebraic expressions. The script is heavily commented, so it should be possible to read through it step-by-step.

Example of execution:

Code:
./makeNewWaveTheory myNewWaveTheory
The output to screen is then:

Code:
A directory called 'myNewWaveTheory' with [C,H]-files is created.

Steps to finalise:
1. Copy the directory to the correct wave type in waves2Foam/src/waveTheories
2. Insert the needed variables and expressions for eta, U, etc.
3. Add a line in waves2Foam/src/Make/files to include the new theory in the compilation.
4. Done.
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 offline   Reply With Quote

Old   February 12, 2014, 08:33
Default
  #709
New Member
 
Jorge Gadelho
Join Date: Feb 2013
Posts: 11
Rep Power: 4
JGadelho is on a distinguished road
Hello Niels,
first of all, sorry for the newbie question, but I don't have the time to explore more and I guess it may help in the future other users.

I have OF 2.2.1 with Waves2Foam (probably R1984) working.
I need to use the surfaceElevation function with InterDyMFoam.

I added this code to the ControlDict file:

Code:
libs
(
    "libwaves2Foam.so"
    "libwaves2FoamSampling.so"
    "libOpenFOAM.so"
    "libincompressibleRASModels.so"
    "libfvMotionSolvers.so"
    "libforces.so"
    
);

functions
{
    surfaceElevation
    {
        type            surfaceElevation;
        functionObjectLibs ( "libwaves2Foam.so" );
        
        outputControl   timeStep; // Alternative: outputTime
        outputInterval  1;
        
        // Additional output controls in waves2Foam
        // Set outputInterval = 1 and timeStep output control if surfaceSampleDeltaT is used.
        samplingStartTime 0; 
        surfaceSampleDeltaT 0.025; // Only an approximate output interval! No feedback on deltaT in the simulation
        
        setFormat raw;
        interpolationScheme cellPointFace;
        fields (alpha1);
        sets
        (
            gauge0
            {
                type        face;
                axis        y;

                start       (0.6 -0.3 0.005);
                end         (0.6 0.2 0.005);
                nPoints     100;
            }  
        
            
        );
    }
}
The surfaceElevation function writes the first time step, but then the simulation crashes with this error message:

Code:
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::valuePointPatchField<double>::updateCoeffs() at ??:?
#4  Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>::GeometricBoundaryField::updateCoeffs() at ??:?
#5  Foam::velocityComponentLaplacianFvMotionSolver::solve() at ??:?
#6  Foam::motionSolver::newPoints() at ??:?
#7  Foam::dynamicMotionSolverFvMesh::update() at ??:?
#8  
 at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  
 at ??:?
Am I missing something?

thank you!
JGadelho is offline   Reply With Quote

Old   February 12, 2014, 13:08
Default WaveFoam error
  #710
New Member
 
Luis Miguel
Join Date: Apr 2013
Location: Colombia
Posts: 13
Rep Power: 4
luigi21 is on a distinguished road
Hello everybody..

I'm trying to run a waveflume case in order to compare some results with experimental data, I've created a flume with (40mx0.5mx0.1m) (x,y,z), and the resolution in each direction is (0.02, 0.0025, 1), I've got the error below:

ecopetrol@ecopetrol-OptiPlex-755:~/OpenFOAM/OpenFOAM_simulations/coimbra_cluster_ext$ [0] [1] ##0 0 Foam::error:rintStack(Foam::Ostream&)Foam::error :rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 at ??:?
[1] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) at ??:?
[1] #4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) at ??:?
[0] #4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) at ??:?
[1] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) at ??:?
[0] #5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) at ??:?
[1] #6 at ??:?
[0] #6

[1] at ??:?
[1] #7 __libc_start_main[0] at ??:?
[0] #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #8 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #8

[0] at ??:?
[ecopetrol-OptiPlex-755:07467] *** Process received signal ***
[ecopetrol-OptiPlex-755:07467] Signal: Floating point exception (8)
[ecopetrol-OptiPlex-755:07467] Signal code: (-6)
[ecopetrol-OptiPlex-755:07467] Failing at address: 0x3e800001d2b
[1] at ??:?
[ecopetrol-OptiPlex-755:07468] *** Process received signal ***
[ecopetrol-OptiPlex-755:07468] Signal: Floating point exception (8)
[ecopetrol-OptiPlex-755:07468] Signal code: (-6)
[ecopetrol-OptiPlex-755:07468] Failing at address: 0x3e800001d2c
[ecopetrol-OptiPlex-755:07467] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f882b4fb0b0]
[ecopetrol-OptiPlex-755:07467] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f882b4fb037]
[ecopetrol-OptiPlex-755:07467] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f882b4fb0b0]
[ecopetrol-OptiPlex-755:07467] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geo metricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEER KT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7vo lMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMes hEEESK_RKT0_RKT1_ddi+0xe3c) [0x7f882d66be1c]
[ecopetrol-OptiPlex-755:07467] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES5limitINS_17geome tricOneFieldENS_9zeroFieldES3_EEvRKT_RKNS_14Geomet ricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS7_Id NS_13fvsPatchFieldENS_11surfaceMeshEEERSF_RKT0_RKT 1_ddib+0x244) [0x7f882d66fad4]
[ecopetrol-OptiPlex-755:07467] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveER NS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMesh EEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEER S8_dd+0x79) [0x7f882d660ba9]
[ecopetrol-OptiPlex-755:07467] [ 6] waveFoam() [0x42fd9b]
[ecopetrol-OptiPlex-755:07467] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f882b4e5ea5]
[ecopetrol-OptiPlex-755:07468] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f37f121d0b0]
[ecopetrol-OptiPlex-755:07467] [ 8] waveFoam() [0x430f35]
[ecopetrol-OptiPlex-755:07467] *** End of error message ***
[ecopetrol-OptiPlex-755:07468] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f37f121d037]
[ecopetrol-OptiPlex-755:07468] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x370b0) [0x7f37f121d0b0]
[ecopetrol-OptiPlex-755:07468] [ 3] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geo metricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEER KT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7vo lMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMes hEEESK_RKT0_RKT1_ddi+0xe6e) [0x7f37f338de4e]
[ecopetrol-OptiPlex-755:07468] [ 4] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES5limitINS_17geome tricOneFieldENS_9zeroFieldES3_EEvRKT_RKNS_14Geomet ricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS7_Id NS_13fvsPatchFieldENS_11surfaceMeshEEERSF_RKT0_RKT 1_ddib+0x244) [0x7f37f3391ad4]
[ecopetrol-OptiPlex-755:07468] [ 5] /opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveER NS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMesh EEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEER S8_dd+0x79) [0x7f37f3382ba9]
[ecopetrol-OptiPlex-755:07468] [ 6] waveFoam() [0x42fd9b]
[ecopetrol-OptiPlex-755:07468] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f37f1207ea5]
[ecopetrol-OptiPlex-755:07468] [ 8] waveFoam() [0x430f35]
[ecopetrol-OptiPlex-755:07468] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 7467 on node ecopetrol-OptiPlex-755 exited on signal 8 (Floating point exception).



I attached the whole case looking forward someone could catch some error, I'd appreciate it so much, thank you guys.

coimbra_case.zip
luigi21 is offline   Reply With Quote

Old   February 12, 2014, 13:21
Default
  #711
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Jorge,

It will be quite difficult for anyone running in OF2.2.1 to help you as long as you do not deliver a running case.
Also note, that I do not have the time to resolve this in the code without an actual solution/understanding of the issue.

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 offline   Reply With Quote

Old   February 13, 2014, 04:01
Default
  #712
New Member
 
Join Date: Oct 2013
Posts: 4
Rep Power: 3
weezer579 is on a distinguished road
Helllo Niels,

Thank you very much for this great toolbox. It really saves me tons of time in my research.

I am trying to simulate a 2nd order Stokes wave (using waveFlume from tutorial) with the following properties:

stokesSecond
wave height: 0.3
wave number: 1.70048
omega: 3.1415
wave depth: 2
relaxationZone: default

The result looks very nice in the relaxationZone (Inlet). However, right after passing the this region, the wave crest seems to be pushed further to the right and making the wavelength longer, wave height shorter. I am not sure if this could be due to the my wrong set up with relaxationZone or something else. I have attached the surface elevation to this post. I'd really appreciate if you could take a quick look.

Thank you again.

Last edited by weezer579; February 14, 2014 at 00:32.
weezer579 is offline   Reply With Quote

Old   February 13, 2014, 04:16
Default
  #713
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Good morning Weezer,

Your specifications of water depth, cyclic frequency and wave number does not match the linear dispersion relation. Use the setWaveParameters utility to correct this issue.

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 offline   Reply With Quote

Old   February 13, 2014, 06:30
Default
  #714
New Member
 
Jorge Gadelho
Join Date: Feb 2013
Posts: 11
Rep Power: 4
JGadelho is on a distinguished road
Quote:
Originally Posted by ngj View Post
Hi Jorge,

It will be quite difficult for anyone running in OF2.2.1 to help you as long as you do not deliver a running case.
Also note, that I do not have the time to resolve this in the code without an actual solution/understanding of the issue.

Kind regards,

Niels
Hello Niels,
my initial question was:

if it is possible to couple surfaceElevation with interDyMFoam or waveDyMFoam?

I've made several tests, and surfaceElevation only works on static mesh solvers, like waveFoam or interFoam. I don't know if I'm making something wrong.

Can you suggest an alternative to read the free surface elevation time series on dynamic mesh solvers?

Thank you again,
Jorge.
JGadelho is offline   Reply With Quote

Old   February 13, 2014, 14:07
Default
  #715
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,608
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Jorge,

Yes, I have used the surfaceElevation utility (runTime and post-processing) successfully on moving meshes in 1.6-ext. It should be noted that my motion is vertical, i.e. no change of the location of the grid lines in the hortizontal direction.

Good luck,

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

Old   February 14, 2014, 00:31
Default
  #716
New Member
 
Join Date: Oct 2013
Posts: 4
Rep Power: 3
weezer579 is on a distinguished road
Good morning Neils,

Thank you very much for your reply. I did try your suggestion and changed my omega value to 4.2129 to satisfy the dispersion relation. The result looks a lot better and matches the theoretical shape very well. However, the amplitude is still smaller than prediction. I've attached the results to this post. Would you mind taking a look again?

Thank you again for your time and help.

weezer579 is offline   Reply With Quote

Old   February 15, 2014, 18:42
Default
  #717
New Member
 
Ji
Join Date: Oct 2013
Posts: 3
Rep Power: 3
aisy is on a distinguished road
Hi Neils,
Can you tell me how to set the debug in waveProperties for stokesSecond ?
Is it right if I set debug as on ?
Thanks!
aisy is offline   Reply With Quote

Old   February 18, 2014, 10:45
Default FreeSurface sampling on a plane
  #718
New Member
 
Ed Ransley
Join Date: Jul 2012
Posts: 29
Rep Power: 5
Ed R is on a distinguished road
Dear all,
I've been struggling recently. I want to plot/output the location of the freesurface on a slice/plane of a 3D case. The function object (or sampling method) for the isoSurfaceCell recordeds the isoValue 0.5 for the entire domain (which is typically very large). Has anyone managed to record the position of the freesurface on a plane dissecting the domain?

Thanks in advance.

Ed
Ed R is offline   Reply With Quote

Old   February 19, 2014, 00:16
Default Draft sensitive solution
  #719
New Member
 
Nima
Join Date: Feb 2012
Location: Perth, Western Australia
Posts: 13
Rep Power: 5
nima3906m is on a distinguished road
Dear all

I am trying to simulate a floating body in 2D. The body is a simple rectangular box and the wave is airy. To my understanding the boundary condition for the floating object should be movingWallVelocity and the mass of the object should be equal to the mass of the displaced water. However my problem starts here! when I try to change the body draft and make it lower say half of the height of the body in the water the simulation crashes. While when it is 0.8 or higher the simulation runs perfectly. I can't understand what is wrong the only thing I see is that there exist undesired pressure region under the body for low draft values! and before it crashes it seems that the body has buoyancy problem . ( Do I calculate the mass in a wrong way ? ) Does it have a meaning for the centre of mass to be at sea level?
nima3906m is offline   Reply With Quote

Old   February 22, 2014, 07:40
Default
  #720
New Member
 
Nick Krgs
Join Date: Jan 2014
Posts: 11
Rep Power: 3
Nick_civ is on a distinguished road
P { margin-bottom: 0.08in; } Hi Neils,

I am trying to simulate a case in which an object (e.g. breakwater) is inserted into the waveFlume just before the relaxation zone. I created an STL file and i executed
1. blockMesh
2. Snappyhexmesh - overwrite
3. Create the new boundary in 0 dir
3. Extrudemesh
4. Delete some files from Polymesh
5. waveGaugesNProbes
6. setWaveField
but when i run it, it doesn't work - here s the output

No finite volume options present

time step continuity errors : sum local = 0, global = 0, cumulative = 0
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.0011976
Time = 0.0011976

MULES: Solving for alpha1
Phase-1 volume fraction = 0.217392 Min(alpha1) = 0 Max(alpha1) = 1
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::relaxationSchemes::relaxationSchemeSpatial:: correct() in "/home/nikos/OpenFOAM/nikos-2.2.2/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#4 Foam::relaxationZone::correct() in "/home/nikos/OpenFOAM/nikos-2.2.2/platforms/linux64GccDPOpt/lib/libwaves2Foam.so"
#5
in "/home/nikos/OpenFOAM/nikos-2.2.2/platforms/linux64GccDPOpt/bin/waveFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/home/nikos/OpenFOAM/nikos-2.2.2/platforms/linux64GccDPOpt/bin/waveFoam"
Floating point exception (core dumped)

I attached the case and i would be grateful if you took a look. (I didn't attach Postprocessed2foam, surfaceElevationAnyname and wavegaugesnprobes. Also, 0 directory's files include the new boundary condition of breakwater and they are just before running and the stl file in the triface is txt due to upload limit here)

Thank you in advance.
Nick
Attached Files
File Type: gz test7.tar.gz (22.9 KB, 8 views)
Nick_civ 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 15:20.