Saturday, December 22, 2018

Reynolds stress transport models notes

Reynolds stress transport models

https://www.openfoam.com/documentation/cpp-guide/html/guide-turbulence-ras-reynolds-stress-transport-models.html


CFD Online - Reynolds stress models

https://www.cfd-online.com/Wiki/Reynolds_stress_model_(RSM)

RANS turbulence modeling of a square channel in openFOAM (part II)


https://www.youtube.com/watch?v=j3U7-17cZ6E

Open Foam Tutorials:::

https://vdocuments.site/basic-open-foam-tutorials-guide.html

======================

To work on convergence:::

https://cfd.direct/openfoam/user-guide/v6-controlDict/ 


Mod-06 Lec-38 Two -equation model for turbulent flow; Numerical calculation of turbulent


https://www.youtube.com/watch?v=yGUg3WV3QLE



k-epsilon Turbulence Model

https://www.youtube.com/watch?v=Rq0AaSnZ1GM

https://www.cfd-online.com/Wiki/V2-f_models

RSM Models:

https://www.openfoam.com/documentation/cpp-guide/html/guide-turbulence-ras-reynolds-stress-transport-models.html

About LRR/SSG:

https://www.cfd-online.com/Forums/openfoam/82493-lrr-launder-reece-rodi-rstm-turbulent-model.html

https://www.cfd-online.com/Forums/openfoam/82493-lrr-launder-reece-rodi-rstm-turbulent-model-2.html

http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/Berlin2008/Papers/OSCIC-08_KarvinenAku.pdf

https://www.cfd-online.com/Forums/openfoam-solving/212081-les-models-reynolds-stress.html


https://www.slideshare.net/MarinaFilonovich/phdthesismfilonovichapril2015-49741564


https://github.com/OpenFOAM/OpenFOAM-dev/tree/76213ac30147e5d63d63354b41d9e4ffcdef2392/src/TurbulenceModels/turbulenceModels/RAS/SSG

https://github.com/OpenFOAM/OpenFOAM-dev/tree/76213ac30147e5d63d63354b41d9e4ffcdef2392/src/TurbulenceModels/turbulenceModels/RAS/LRR

https://www.cfd-online.com/Forums/openfoam/122190-simulation-axisymmetric-free-jet-using-lrr-reynolds-stress-model.html


Also useful for convergence:

https://cfd.direct/openfoam/user-guide/v6-controlDict/

More convergence links:

https://www.cfd-online.com/Forums/openfoam/82493-lrr-launder-reece-rodi-rstm-turbulent-model.html


https://www.cfd-online.com/Forums/openfoam-solving/58365-simplefoam-convergence-problems.html


https://www.cfd-online.com/Forums/openfoam/82183-where-log-file-created.html

Friday, December 21, 2018

CFD slot die results for kinematic viscosity of 0.2 (turbulence turned off)

in case/system/fvSolution

........

relaxationFactors
{
    equations
    {
        U               0.85; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.85; // 0.9 is more stable but 0.95 more convergent
    }
}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}

........

The program crashes with the turbulence turned off



===================================================
 relaxationFactors
{
    equations
    {
        U               0.9; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.9; // 0.9 is more stable but 0.95 more convergent
    }
}





===================================================

relaxationFactors
{
    equations
    {
        U               0.95; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.95; // 0.9 is more stable but 0.95 more convergent
    }
}


===================================================

relaxationFactors
{
    equations
    {
        U               0.98; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.98; // 0.9 is more stable but 0.95 more convergent
    }
}








Hmm...I wonder if I can run things longer for the higher relaxation factors.


====================================
fvSolution:::

relaxationFactors
{
    equations
    {
        U               0.98; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.98; // 0.9 is more stable but 0.95 more convergent
    }
}

controlDict::

stopAt          endTime;

endTime         20;    //// this is double the previous 0.98 relaxation run



==================================================

fvSolution:::

relaxationFactors
{
    equations
    {
        U               0.98; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.98; // 0.9 is more stable but 0.95 more convergent
    }
}

controlDict::

stopAt          endTime;

endTime         30;    //// this is double the previous 0.98 relaxation run


===================================================

/constant/turbulenceProperties/

simulationType RAS;

RAS
{
    // Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, v2f,
    // ShihQuadraticKE, LienCubicKE.
    RASModel        kEpsilon;

    turbulence      on;

    printCoeffs     on;
}



// ************************************************************************* //

with everything else (i.e. relaxationFactors, ControlDict) like the previous run!


ahhh.... crashing...

==================================================

fvSolution::

.......

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent      yes;

    residualControl
    {
        p               1e-1;
        U               1e-1;
        "(k|epsilon|omega|f|v2)" 1e-3;
    }
}

relaxationFactors
{
    equations
    {
        U               0.98; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.98; // 0.9 is more stable but 0.95 more convergent
    }
}

............

controlDict::

startTime       0;

stopAt          endTime;

endTime         30;


hm-trp-tool configuration notes

 1985  cd Downloads/
 1986  ls
 1987  cd hm-trp-tool/
 1988  ls
 1989  cat main.cpp
 1990  nano main.cpp
 1991  qmake Configurator.pro && make
 1992  nano main.cpp
 1993  qmake Configurator.pro && make
 1994  ls
 1995  cat configurator.h
 1996  ls
 1997  cat configurator.h | less
 1998  locate qmainwindow
 1999  nano configurator.h
 2000  qmake Configurator.pro && make
 2001  nano configurator.h
 2002  qmake Configurator.pro && make
 2003  locate ui_creator.h
 2004  locate uicreator
 2005  locate uicreator.h
 2006  locate creator
 2007  locate creator,ui
 2008  locate creator.ui
 2009  qmake Configurator.pro && make
 2010  ls
 2011  cat configurator.cpp
 2012  cat configurator.cpp | less
 2013  nano configurator.cpp
 2014  qmake Configurator.pro && make
 2015  ls
 2016  history


in main.cpp

#include <QtWidgets/qapplication.h>

/usr/include/x86_64-linux-gnu/qt5/QtWidgets/qapplication.h

in configurator.h

#include <QtWidgets/qmainwindow.h>



http://doc.qt.io/qt-5/designer-using-a-ui-file.html

http://doc.qt.io/qt-5/qtwidgets-index.html

NMEA data

root@brent-P65-67HSHP:/home/brent# cat /dev/ttyUSB0
$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$PSRF150,1*3E

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,01,00,000,,02,00,000,,03,00,000,,04,00,000,*7C

$GPGSV,3,2,12,05,00,000,,06,00,000,,07,00,000,,08,00,000,*77

$GPGSV,3,3,12,09,00,000,,10,00,000,,11,00,000,,12,00,000,*71

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,21,00,000,20,01,00,000,,02,00,000,,03,00,000,*79

$GPGSV,3,2,12,04,00,000,,05,00,000,,06,00,000,,07,00,000,*7B

$GPGSV,3,3,12,08,00,000,,09,00,000,,10,00,000,,11,00,000,*7A

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,21,00,000,20,01,00,000,,02,00,000,,03,00,000,*79

$GPGSV,3,2,12,04,00,000,,05,00,000,,06,00,000,,07,00,000,*7B

$GPGSV,3,3,12,08,00,000,,09,00,000,,10,00,000,,11,00,000,*7A

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPRMC,,V,,,,,,,,,,N*53

$GPGGA,,,,,,0,00,,,M,0.0,M,,0000*48

$GPGSA,A,1,,,,,,,,,,,,,,,*1E

$GPGSV,3,1,12,21,00,000,21,01,00,000,,02,00,000,,03,00,000,*78

$GPGSV,3,2,12,04,00,000,,05,00,000,,06,00,000,,07,00,000,*7B

$GPGSV,3,3,12,08,00,000,,09,00,000,,10,00,000,,11,00,000,*7A

$GPRMC,,V,,,,,,,,,,N*53

^C

https://www.eecis.udel.edu/~mills/ntp/html/drivers/driver20.html

https://unix.stackexchange.com/questions/361534/how-to-read-nmea-data-from-the-holux-m-1200-gps-receiver

https://www.systutorials.com/docs/linux/man/n-nmea/

https://www.nmea.org/content/nmea_standards/nmea_2000_ed3_10.asp

https://www.linuxquestions.org/questions/linux-hardware-18/reading-nmea-from-gps-and-output-it-in-the-usb-port-727476/

https://www.gpsinformation.org/dale/nmea.htm?fbclid=IwAR3aSp8MVAVGuWSQpPolshKg-GYBTG4_p0FuaYPOEddNr4Tpxzqaj2QQtJ8

Parse NMEA data with:

https://github.com/Knio/pynmea2

Wednesday, December 19, 2018

CFD slot die results for kinematic viscosity of 0.2 (turbulence turned on)

relaxationFactors
{
    equations
    {
        U               0.9; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.9; // 0.9 is more stable but 0.95 more convergent
    }

https://www.cfd-online.com/Forums/openfoam-solving/58365-simplefoam-convergence-problems.html


https://www.cfd-online.com/Forums/openfoam-solving/107559-initialization-potentialfoam.html

relaxationFactors
{
    equations
    {
        U               0.85; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.85; // 0.9 is more stable but 0.95 more convergent
    }


relaxationFactors
{
    equations
    {
        U               0.95; // 0.95 is more stable but 0.95 more convergent
        ".*"            0.95; // 0.95 is more stable but 0.95 more convergent
    }


relaxationFactors
{
    equations
    {
        U               0.98; // 0.95 is more stable but 0.95 more convergent
        ".*"            0.98; // 0.95 is more stable but 0.95 more convergent
    }


Monday, December 17, 2018

Apogee Components::: Building the Level-2 Fiberglass Rocket Kit - 13: Mount Ejection Canisters
https://www.youtube.com/watch?v=KQ1gGyLwB2c

Apogee Torrent Instructions:

https://www.youtube.com/watch?v=5w0LJ7Hw72w&list=PLmJQzoMsBYwPal0Dce79sJz9qa2XA8gBU