[evxta3]

Hi IHFOAM team,

Is there any compressible solver (e.g. compressibleinterfoam) already implemented in IHFOAM?

Thanks,
Tommaso

[IHFOAM Team]

Hi Tommaso,

Yes, you can use the waves library in a compressible solver (compressibleInterFoam).

And yes, we have our own version of a porous medium solver for compressible cases. Please, check this presentation:
https://www.researchgate.net/publication/361998233_Novel_Methodology_for_a_Fast_3D_Numerica l_Analysis_of_The_PTO_Damping_Force_on_a_Dual-Chamber_OWC


Best,
IHFOAM Team

[IHFOAM Team]

Hi AlxB,

Where could these come from ?

I recommend you to first run the cnoidal case (laminar) provided in the tutorials and check if if works correctly; if so, there is something wrong with the numerical setup of your turbulence model.

Please, check carefully that the initial seeding of the turbulence is correct:
https://www.cfd-online.com/Wiki/Turbulence_free-stream_boundary_conditions

Best,
IHFOAM Team

[evxta3]

Hi,

I modified the wave/solitary tutorial using the solver compressibleInterFoam, however it gives the following error:

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting surfaceTensionModel constant
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectFluid;
specie specie;
energy sensibleInternalEnergy;
}

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib64/libpthread.so.0
#3 Foam::heThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::he(Foam::Field<double> const&, Foam::Field<double> const&, int) const at ??:?
#4 Foam::heThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::init(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#5 Foam::heThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam::rhoThermo::addfvMeshConstructorToTable<Foam: :heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7 Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMe sh const&, Foam::word const&) at ??:?
#8 Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#9 Foam::twoPhaseMixtureThermo::twoPhaseMixtureThermo (Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#10 ? at ??:?
#11 __libc_start_main in /lib64/libc.so.6
#12 ? at ??:?



Any ideas where this comes from? Thanks!

[IHFOAM Team]

Hi Tommaso

Originally Posted by evxta3

Any ideas where this comes from? Thanks!

I am afraid we cannot understand or fix the error with just the log file.

Sorry,
IHFOAM Team

p.s: Just curious, does your surname start with an A?

[evxta3]

Hi IHFOAM Team,

Thanks for your reply. We menaged to fix the error and to run IHFOAM with a compressible solver.

Originally Posted by IHFOAM Team

p.s: Just curious, does your surname start with an A?

Yes, it does start with an A .

Thanks,
Tommaso

[IHFOAM Team]

Great, well done !!!

( Good luck with the air pockets )

IHFOAM Team

[lucia_995]

Hi, I have problem to generate irregular waves with one direction. Can you help me?
Thanks

[IHFOAM Team]

Hola Lucia,

Originally Posted by lucia_995

I have problem to generate irregular waves with one direction. Can you help me?
Thanks

After the summer nothing works as it should do, but we keep on foaming . Can you give us more info about your problem?

Best,
IHFOAM Team

[lucia_995]

I try to reproduce a channel with a slope and a rubble mound breakwater, regular and irregular waves conditions.
With regular waves I haven't problems but with irregular waves I don't understand gow to setup the wavePropeerties file.
I see the case irregularMultiDirectional in interFoam and I notice that there are a "series" of periods, heights, etc
I don't understand how to derive this series and I want know if there is a way to setup only Tp and Hs like in waves2Foam.

Thank you very much!

[lucia_995]

Originally Posted by IHFOAM Team

Hola Lucia,



After the summer nothing works as it should do, but we keep on foaming . Can you give us more info about your problem?

Best,
IHFOAM Team

I try to reproduce a channel with a slope and a rubble mound breakwater, regular and irregular waves conditions.
With regular waves I haven't problems but with irregular waves I don't understand gow to setup the wavePropeerties file.
I see the case irregularMultiDirectional in interFoam and I notice that there are a "series" of periods, heights, etc
I don't understand how to derive this series and I want know if there is a way to setup only Tp and Hs like in waves2Foam.

Thank you very much!

[IHFOAM Team]

Hi there again,

Originally Posted by lucia_995

I see the case irregularMultiDirectional in interFoam and I notice that there are a "series" of periods, heights, etc

You can find the answer in the great book "Wavess in oceanic and coastal waters" " by Leo H. Holthuijsen (I believe is in page 44 a clear explanation).

Best,
IHFOAM TEAM

p.d: Sorry, you can not just define Tp or Hs in our model, we believe that following our choice is much better to generate "fast & correct" irregular cases.

[bknorris]

Originally Posted by lucia_995

I try to reproduce a channel with a slope and a rubble mound breakwater, regular and irregular waves conditions.
With regular waves I haven't problems but with irregular waves I don't understand gow to setup the wavePropeerties file.
I see the case irregularMultiDirectional in interFoam and I notice that there are a "series" of periods, heights, etc
I don't understand how to derive this series and I want know if there is a way to setup only Tp and Hs like in waves2Foam.

Thank you very much!

Hi Lucia,

IHFOAM does contain a tutorial case for irregular waves if you haven't found it already,

<OpenFOAM directory>/tutorials/multiphase/interFoam/laminar/waves/irregularMultiDirection

Note this is for multi-directional irregular waves and so to make the waves constant from one direction you would just set 'waveDirs' equal to 0 for each period, height, and phase.

To estimate wave period, height, and phase, well it depends on what data you're working with. For example, when I was using the irregular wave model, I computed the fft of a sea surface pressure time series recorded by an instrument, and used the complex spectrum to calculate wave height (= 2*amplitude), period (1/frequency), and wave phase for each frequency component. If you lack instrument data, you could simply use a program (e.g., https://www.mathworks.com/matlabcent...2-wavemodeling) to generate a spectrum and then calculate these components yourself.

Hope that helps,

Ben