 |
 |
 |
|
|
|
[Sponsors] | ||||
compressibleInterDyMFoam and the troubles of the right thermophysical Parameters |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 6, 2017, 10:27
|
compressibleInterDyMFoam and the troubles of the right thermophysical Parameters
|
#1 |
|
New Member
Oliver K
Join Date: May 2017
Posts: 15
Rep Power: 8  |
Hey there,
I'm solving right now my first compressible Case and got a little bit of a hard time with the right thermophysical parameters. Actually I've run the case with interDyMFoam without any errors, so I don't think that the error could be provoked of the normal boundary conditions. That's how the error message looks like: Code:
Create mesh for time = 0
Selecting dynamicFvMesh solidBodyMotionFvMesh
AMI: Creating addressing and weights between 123 source faces and 2304 target faces
AMI: Patch source sum(weights) min/max/average = 0, 0.0644961, 0.00100825
AMI: Patch target sum(weights) min/max/average = 0, 1, 0.00392754
Selecting solid-body motion function tabulated6DoFMotion
Applying solid body motion to cellZone movingZone
Reading g
PIMPLE: max iterations = 100
field p_rgh : relTol 0, tolerance 5e-05
field U : relTol 0, tolerance 0.0001
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Constructing twoPhaseMixtureThermo
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectFluid;
specie specie;
energy sensibleInternalEnergy;
}
#0 Foam::error::printStack(Foam::Ostream&) in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 ? in "/lib64/libc.so.6"
#3 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::he(Foam::Field<double> const&, Foam::Field<double> const&, int) const in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::init() in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"#5 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heRhoThermo(Foam::fvMesh const&, Foam::word const&) in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#7 Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/projects2/OF_bin/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#8 Foam::autoPtr<Foam::rhoThermo> [...]
Floating point exception (core dumped)
Code:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectFluid;
specie specie;
energy sensibleInternalEnergy;
}
mixture
{
nMoles 1;
molWeight 18.0;
}
equationOfState
{
R 3000;
rho0 1027;
}
thermodynamics
{
Cp 4195;
Hf 0;
}
transport
{
mu 3.645e-4;
Pr 2.289;
}
}
Code:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;//was sensibleInternalEnergy latest sensibleEnthalpy
}
mixture
{
specie
{
nMoles 1;
molWeight 28.96;
}
thermodynamics
{
Cp 1004.5;//was 1007
Hf 2.544e+06;//was 0
}
transport
{
mu 1.8e-05;//was 1.84e-05
Pr 0.7;
}
}
//Edit: So for I got the first part to run. The issue results of incorrect T file for my ACMI BC. Actually the first part is loading but the error still exists for the air Part. Any ideas which BC is incorrect? Code:
internalField uniform 293.15;
boundaryField
{
atmosphere
{
type calculated;
value $internalField;
}
kaverneWall
{
type zeroGradient;
}
outlet
{
type inletOutlet;//was zeroGradient
inletValue uniform 293.15;
value $internalField;
}
pipeWall
{
type zeroGradient;
}
AMI
{
type zeroGradient;
}
slide_interface
{
type zeroGradient;
}
ACMI1_couple
{
type cyclicACMI;
value $internalField;
}
ACMI1_blockage
{
type zeroGradient;
// type calculated;
// value uniform 0;
}
ACMI2_couple
{
type cyclicACMI;
value $internalField;
}
ACMI2_blockage
{
type zeroGradient;
//type calculated;
//value uniform 0;
}
}
I want to simulate a the outflow of a tank through a pipe which is closed after some time with a valve and want to calculate the water hammer effect. The patch atmosphere is the upper part of the tank and the rest explains itself i hope. Cheers, Oli Last edited by silencebreak; September 8, 2017 at 08:20. |
|
|
|
|
|
September 8, 2017, 08:19
|
|
#2 |
|
New Member
Oliver K
Join Date: May 2017
Posts: 15
Rep Power: 8 |
Problem is solved. For everyone other expecting similar errors: You've got to check your T and especially for the ACMI-Patches you need to specify the value as internalFields.
|
|
|
|
|
|
|
August 26, 2020, 04:01
|
|
#3 |
|
New Member
Prasad ADHAV
Join Date: Apr 2020
Location: Belval, Luxembourg
Posts: 10
Rep Power: 6 |
I have the same issue.
I am trying to simulate a water jet through nozzle. But I am unsure of the boundary conditions, since no tutorial really covers a good case. I have the following 0/T: Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
wall
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
inlet2
{
type zeroGradient;
}
inlet3
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
I get the following error: Code:
Create time
Create mesh for time = 0
Selecting dynamicFvMesh staticFvMesh
PIMPLE: No convergence criteria found
PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
Constructing twoPhaseMixtureThermo
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo eConst;
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::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7 Foam::twoPhaseMixtureThermo::twoPhaseMixtureThermo(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) at ??:?
#8 ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/compressibleInterFoam"
Floating point exception (core dumped)
Thank you for the help. |
|
|
|
|
|
|
November 19, 2020, 10:26
|
May just set everything to 300 Kelvin
|
#4 |
|
New Member
Goeran Gabriel
Join Date: Oct 2016
Posts: 5
Rep Power: 9 |
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2006 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { piston { type fixedValue; value uniform 300; } fixedWall { type fixedValue; value uniform 300; } walls { type fixedValue; value uniform 300; } inlet { type fixedValue; value uniform 300; } outlet { type inletOutlet; inletValue uniform 300; value uniform 300; } } // ************************************************** *********************** // |
|
|
|
|
|
|
January 27, 2021, 01:44
|
|
#5 | |
|
Member
Hector
Join Date: Jul 2010
Location: Barcelona
Posts: 30
Rep Power: 15 |
Quote:
Do you mean specifying value just for ACMI1_couple and ACMI2_couple? Or for ACMI1_blockage and ACMI2_blockage as well? |
||
|
|
|
||
|
March 22, 2023, 08:51
|
|
#6 | |
|
New Member
Join Date: Sep 2021
Posts: 6
Rep Power: 4 |
Quote:
Thank you |
||
|
|
|
||
|
|
|
Linear Mode
