CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[IHFOAM] The IHFOAM Thread

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

Like Tree57Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 7, 2023, 14:49
Default Compressible fluid solver
  #381
New Member
 
Tommaso
Join Date: Jan 2020
Posts: 16
Rep Power: 6
evxta3 is on a distinguished road
Hi IHFOAM team,

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

Thanks,
Tommaso
evxta3 is offline   Reply With Quote

Old   March 8, 2023, 03:02
Default
  #382
Senior Member
 
IHFOAM Team's Avatar
 
IHFOAM The Environmental Hydraulics Institute "IHCantabria"
Join Date: Sep 2017
Location: Santander (Spain)
Posts: 120
Rep Power: 9
IHFOAM Team is on a distinguished road
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
evxta3 likes this.
__________________
http://ihfoam.ihcantabria.com/
IHFOAM Team is offline   Reply With Quote

Old   March 8, 2023, 03:08
Default
  #383
Senior Member
 
IHFOAM Team's Avatar
 
IHFOAM The Environmental Hydraulics Institute "IHCantabria"
Join Date: Sep 2017
Location: Santander (Spain)
Posts: 120
Rep Power: 9
IHFOAM Team is on a distinguished road
Hi AlxB,

Quote:
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
__________________
http://ihfoam.ihcantabria.com/
IHFOAM Team is offline   Reply With Quote

Old   March 15, 2023, 06:39
Default compressibleInterFoam
  #384
New Member
 
Tommaso
Join Date: Jan 2020
Posts: 16
Rep Power: 6
evxta3 is on a distinguished road
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!
evxta3 is offline   Reply With Quote

Old   March 30, 2023, 06:28
Default
  #385
Senior Member
 
IHFOAM Team's Avatar
 
IHFOAM The Environmental Hydraulics Institute "IHCantabria"
Join Date: Sep 2017
Location: Santander (Spain)
Posts: 120
Rep Power: 9
IHFOAM Team is on a distinguished road
Hi Tommaso

Quote:
Originally Posted by evxta3 View Post
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?
__________________
http://ihfoam.ihcantabria.com/
IHFOAM Team is offline   Reply With Quote

Old   March 30, 2023, 06:52
Default
  #386
New Member
 
Tommaso
Join Date: Jan 2020
Posts: 16
Rep Power: 6
evxta3 is on a distinguished road
Hi IHFOAM Team,

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

Quote:
Originally Posted by IHFOAM Team View Post
p.s: Just curious, does your surname start with an A?
Yes, it does start with an A .

Thanks,
Tommaso
evxta3 is offline   Reply With Quote

Old   March 30, 2023, 07:13
Default
  #387
Senior Member
 
IHFOAM Team's Avatar
 
IHFOAM The Environmental Hydraulics Institute "IHCantabria"
Join Date: Sep 2017
Location: Santander (Spain)
Posts: 120
Rep Power: 9
IHFOAM Team is on a distinguished road
Great, well done !!!

( Good luck with the air pockets )

IHFOAM Team
evxta3 likes this.
__________________
http://ihfoam.ihcantabria.com/
IHFOAM Team is offline   Reply With Quote

Old   September 4, 2023, 09:58
Default
  #388
New Member
 
Lucia Ligorio
Join Date: Feb 2022
Posts: 14
Rep Power: 4
lucia_995 is on a distinguished road
Hi, I have problem to generate irregular waves with one direction. Can you help me?
Thanks
lucia_995 is offline   Reply With Quote

Old   September 4, 2023, 11:12
Default
  #389
Senior Member
 
IHFOAM Team's Avatar
 
IHFOAM The Environmental Hydraulics Institute "IHCantabria"
Join Date: Sep 2017
Location: Santander (Spain)
Posts: 120
Rep Power: 9
IHFOAM Team is on a distinguished road
Hola Lucia,

Quote:
Originally Posted by lucia_995 View Post
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
__________________
http://ihfoam.ihcantabria.com/
IHFOAM Team is offline   Reply With Quote

Old   September 4, 2023, 11:38
Default
  #390
New Member
 
Lucia Ligorio
Join Date: Feb 2022
Posts: 14
Rep Power: 4
lucia_995 is on a distinguished road
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 is offline   Reply With Quote

Old   September 4, 2023, 11:39
Default
  #391
New Member
 
Lucia Ligorio
Join Date: Feb 2022
Posts: 14
Rep Power: 4
lucia_995 is on a distinguished road
Quote:
Originally Posted by IHFOAM Team View Post
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!
lucia_995 is offline   Reply With Quote

Old   September 5, 2023, 07:46
Default
  #392
Senior Member
 
IHFOAM Team's Avatar
 
IHFOAM The Environmental Hydraulics Institute "IHCantabria"
Join Date: Sep 2017
Location: Santander (Spain)
Posts: 120
Rep Power: 9
IHFOAM Team is on a distinguished road
Hi there again,

Quote:
Originally Posted by lucia_995 View Post
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.
__________________
http://ihfoam.ihcantabria.com/
IHFOAM Team is offline   Reply With Quote

Old   September 8, 2023, 18:10
Default Irregular waves model
  #393
New Member
 
Benjamin Norris
Join Date: Oct 2020
Location: California
Posts: 17
Rep Power: 5
bknorris is on a distinguished road
Quote:
Originally Posted by lucia_995 View Post
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,

Code:
<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
bknorris is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Divergence detected in AMG solver: k when udf loaded google9002 Fluent UDF and Scheme Programming 3 November 7, 2019 23:34
udf problem jane Fluent UDF and Scheme Programming 37 February 20, 2018 04:17
UDF velocity profile willroca Fluent UDF and Scheme Programming 2 January 10, 2016 03:13
Error messages atg enGrid 7 August 30, 2013 11:16
Phase locked average in run time panara OpenFOAM 2 February 20, 2008 14:37


All times are GMT -4. The time now is 22:07.