philipmorrisintl / aerosolved

AeroSolved is a Computational Fluid Dynamics code, based on the OpenFOAM software package, for simulation of the generation, transport, evolution and deposition of multispecies aerosol mixtures.
GNU General Public License v3.0
55 stars 32 forks source link

'bentPipe' case failing #8

Open fotisAnagnostopoulos opened 1 year ago

fotisAnagnostopoulos commented 1 year ago

Hi, Trying to run 'bentPipe' case using OF2006 (and pull request #3), i encountered very large Reynolds numbers for "V.0i" and segmentation fault. """"""" fullStokes: Solving for V.01, max(Re) = 64934176509.411056519, Initial residual = 0.99999988820481799934, Final residual = 0.16363358688145657194, No Iterations 3

fullStokes: Solving for V.02, max(Re) = 87787037316.140594482, Initial residual = 0.99999987303583404152, Final residual = 0.14972580643221766739, No Iterations 3

fullStokes: Solving for V.03, max(Re) = 106094341466.79673767, Initial residual = 0.99999986628891768259, Final residual = 0.14400644272231580278, No Iterations 3

fullStokes: Solving for V.04, max(Re) = 121744486520.0188446, Initial residual = 0.9999998623041878032, Final residual = 0.14100586997151468016, No Iterations 3

fullStokes: Solving for V.05, max(Re) = 135618136862.88098145, Initial residual = 0.99999985970297511351, Final residual = 0.13916990372873153103, No Iterations 3

fullStokes: Solving for V.06, max(Re) = 148167653095.57357788, Initial residual = 0.99999985786913936359, Final residual = 0.13792752319661194527, No Iterations 3

fullStokes: Solving for V.07, max(Re) = 159764261683.24920654, Initial residual = 0.99999985649623557293, Final residual = 0.13703875455545583284, No Iterations 3

fullStokes: Solving for V.08, max(Re) = 170583779486.52407837, Initial residual = 0.99999985543825653522, Final residual = 0.13637232286533981784, No Iterations 3

fullStokes: Solving for V.09, max(Re) = 180760247145.01428223, Initial residual = 0.99999985460199858345, Final residual = 0.13585436214190790838, No Iterations 3

fullStokes: Solving for V.10, max(Re) = 190396175564.88552856, Initial residual = 0.99999985391475088026, Final residual = 0.13543892275081365062, No Iterations 3

fullStokes: Solving for V.11, max(Re) = 199569202903.19943237, Initial residual = 0.99999985334088925537, Final residual = 0.13509846824438262458, No Iterations 3

fullStokes: Solving for V.12, max(Re) = 208337211493.70684814, Initial residual = 0.99999985285414705327, Final residual = 0.13481554956720773197, No Iterations 3

fullStokes: Solving for V.13, max(Re) = 216754016370.39968872, Initial residual = 0.9999998524351539908, Final residual = 0.13457650466825513136, No Iterations 3

fullStokes: Solving for V.14, max(Re) = 224856779774.82620239, Initial residual = 0.99999985207023056599, Final residual = 0.13437195775841456746, No Iterations 3

fullStokes: Solving for V.15, max(Re) = 232675248324.25796509, Initial residual = 0.9999998517523506214, Final residual = 0.13419506131861025366, No Iterations 3

fullStokes: Solving for V.16, max(Re) = 240241693675.44580078, Initial residual = 0.99999985147149528242, Final residual = 0.1340404677976244352, No Iterations 3

[5] '#'0 Foam::error::printStack(Foam::Ostream&) at ??:? [5]' #'1 Foam::sigFpe::sigHandler(int) at ??:? [5] '#'2 ? in /lib/x86_64-linux-gnu/libc.so.6 [5] '#'3 ? in /lib/x86_64-linux-gnu/libm.so.6 [5] '#'4 exp in /lib/x86_64-linux-gnu/libm.so.6 [5] '#'5 Foam::StokesEinstein::diffusivity(double const&, double const&, double const&) const at ??:? [5] '#'6 Foam::BrownianModel::D(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:? [5] '#'7 Foam::aerosolModels::fixedSectional::updateDrift() at ??:? [5] '#'8 Foam::aerosolModel::correct() at ??:? [5] '#'9 ? in ~/OpenFOAM/fa-v2006/platforms/linux64GccDPInt32Opt/bin/aerosolEulerFoam [5] '#'10 ? in /lib/x86_64-linux-gnu/libc.so.6 [5] '#'11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 [5] '#'12 ? in ~/OpenFOAM/fa-v2006/platforms/linux64GccDPInt32Opt/bin/aerosolEulerFoam """""""

Any ideas/suggestions could be valuable.

Thanks, FA

fort-NA commented 2 months ago

I also met same problem with you. I also do not know how to solve it. Does anybody have some suggestion? Thank you very much. log.txt

fotisAnagnostopoulos commented 2 months ago

Dear @fort-NA, i was able to perform the needed simulations using a standard solver and cloud functions. In particular, i used "buoyantReactingFoam".