|
[Sponsors] |
Dimensions issue using EBRSM turbulence model with chtMultiRegionSimpleFoam |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Lorenzo Faedo
Join Date: Aug 2020
Location: Bologna
Posts: 1
Rep Power: 0 ![]() |
Hello everyone,
I am trying to run a simulation with chtMultiRegionSimpleFoam using the Ellipting Blending Reynolds Stress Model (EBRSM). I am using OpenFOAM-v2306. I started from the case in the tutorial section (EBRSM turbulence model with simpleFoam solver: $FOAM_TUTORIALS/verificationAndValidation/turbulenceModels/planeChannel) and I modified this case to solve it with chtMultiRegionSimpleFoam. The geometry is a rectangular domain with a volumetric momentum source, the BC are two walls at the same temperature in the upper and lower sides and a cyclic condition in the left and right side. When I run the simulation with chtMultiRegionSimpleFoam it crashes during the first time step. here's the output: Code:
Time = 1 Solving for fluid region region1 smoothSolver: Solving for Ux, Initial residual = 0.9999999966456584, Final residual = 1.092316707381895e-11, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 3.15409836055221e-13, No Iterations 11 Pressure gradient source: uncorrected Ubar = 17.54514880973776, pressure gradient = 0.04728518145885395 smoothSolver: Solving for h, Initial residual = 0.9999999417030885, Final residual = 5.504674013964634e-08, No Iterations 1000 Min/max T:373.1499998029909 373.1500011392363 GAMG: Solving for p_rgh, Initial residual = 0.8116908744296613, Final residual = 9.470812287532965e-07, No Iterations 16 GAMG: Solving for p_rgh, Initial residual = 5.8644174102195e-07, Final residual = 6.048763452135292e-08, No Iterations 3 Pressure gradient source: uncorrected Ubar = 17.5451678969471, pressure gradient = 0.02837495316405573 time step continuity errors : sum local = 1.822343471660741e-08, global = -7.188472161815702e-17, cumulative = -7.188472161815702e-17 Min/max rho:1 1 DICPBiCGStab: Solving for f, Initial residual = 1, Final residual = 6.806232488646261e-20, No Iterations 3 --> FOAM FATAL ERROR: (openfoam-2306) Incompatible dimensions for operation [epsilon[1 -1 -4 0 0 0 0] ] + [((Cstability*epsilon)|k)[0 2 -4 0 0 0 0] ] From void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::DimensionedField<Type, Foam::volMesh>&, const char*) [with Type = double] in file /usr/src/packages/BUILD/src/finiteVolume/lnInclude/fvMatrix.C at line 1897. FOAM aborting [stack trace] ============= #1 Foam::error::simpleExit(int, bool) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #2 Foam::error::exiting(int, bool) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so #3 void Foam::checkMethod<double>(Foam::fvMatrix<double> const&, Foam::DimensionedField<double, Foam::volMesh> const&, char const*) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam #4 Foam::tmp<Foam::fvMatrix<double> > Foam::operator+<double>(Foam::tmp<Foam::fvMatrix<double> > const&, Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so #5 Foam::RASModels::EBRSM<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so #6 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam #7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #8 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionSimpleFoam ============= Aborted I have already tried to switch to another turbulence model (realizableKE) and with that the simulation runs perfectly. I've also tried to run another case with EBRSM turbulence using rhoSimpleFoam as a solver and the error is the same. I am struggling to find out how to solve this problem, maybe the EBRSM turbulence model works only with incompressible solvers? I have attached the case, any help is appreciated! Regards Lorenzo |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Engineer
Join Date: Mar 2013
Posts: 27
Rep Power: 13 ![]() |
Hi Lorenzo,
Were you able to find the root cause of the error. I have been trying to implement the same turbulence model and end up having the same error. The EBRSM model is compatible with compressible and incompressible flows. Regards Nagaraj |
|
![]() |
![]() |
![]() |
Tags |
chtmultiregionsimplefoam, ebrsm, turbulence |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
PANS Turbulence Model Compilation Issue - OpenFOAM v-11 | Ishak.Mazing | OpenFOAM Programming & Development | 0 | December 6, 2023 10:30 |
New turbulence model in OpenFOAM 5 with porosity scalar | shinjanghosh | OpenFOAM Programming & Development | 0 | January 15, 2020 00:06 |
Only two turbulence options available in CFX Pre | Jack001 | CFX | 5 | March 30, 2016 02:47 |
Overflow Error in Multiphase Modelling with Two Continuous Fluids | ashtonJ | CFX | 6 | August 11, 2014 14:32 |
Compressible turbulence model issues | 351Cleveland | OpenFOAM | 5 | October 24, 2013 15:41 |