Sunday, December 30, 2018

For When It Burns (preventing the Magic Smoke in Electronics)

https://en.wikipedia.org/wiki/Electronic_circuit_simulation

https://en.wikipedia.org/wiki/Quite_Universal_Circuit_Simulator
Posted by at No comments:

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
Posted by at No comments:

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;


Posted by at No comments:

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
Posted by at No comments:

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

Posted by at No comments:

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
    }


Posted by at No comments:

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
Posted by at No comments:

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




 Solver: SimpleFoam

File: B2CoarseMesh2
https://github.com/bshambaugh/openfoam-experiments/tree/master/B2_CoarseMesh2

In  B2CoarseMesh2/constant/transportProperties:

nu [0 2 -1 0 0 0 0] 0.2;
transportModel Newtonian; 
In B2_CoarseMesh2/constant/turbulenceProperties:


simulationType RAS;

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

    turbulence      off;

    printCoeffs     on;
}



For turbulence turned on, check out:

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

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
    }

For vector field plots in ParaView See:

[ParaView] Creating Vector Fields
https://www.youtube.com/watch?v=cygVdhn-kG0



Convergence (Residual Consideration):
How To Plot Residuals in OpenFOAM
 https://www.youtube.com/watch?v=ngyBbYn5-V0

How to Plot and Control simpleFoam Residuals in OpenFoam
 https://www.youtube.com/watch?v=diCbJMUjZlk


How to plot variables:

OpenFOAM v4.0 Data Analysis: Quick Post-Processing of a Graph
https://www.youtube.com/watch?v=IYBLfmFA7Gc

Variable Extraction from results:
Extracting variables from any case in OpenFoam (use of foamLog command)   https://www.youtube.com/watch?v=HK3xutgOyI0


Calculation and Plotting in Paraview:
Flow and pressure integrals:
https://www.youtube.com/watch?v=ZcfJzbd6mwQ

Paraview CFD Tutorial - Advanced Postprocessing in Paraview
 https://www.youtube.com/watch?v=5M6VgXyHB60

Bonus:
How to calculate aerodynamic forces with Paraview
 https://www.youtube.com/watch?v=J944HOj_4b0

Across the jet region:



 
Cross section and the midpoint of the y-axis (die face axis) from jet:


Posted by at No comments:

Wednesday, December 12, 2018

Bond Graphs Links

AUTOMATED MODELING AND SIMULATION USING THE BOND GRAPH METHOD FOR THE AEROSPACE INDUSTRY 
http://people.idsia.ch/~andrea/sim/simbond.html   >>> http://mtt.sourceforge.net/   (Transform tool for Octave)

Object–oriented modeling with bond graphs and Modelica
https://www.modelica.org/publications/papers/icbgm99bnkw.pdf

Introduction to Physical Systems Modelling with Bond Graphs
http://www.me.umn.edu/~lixxx099/ME8287_S11/BondGraphsV2.pdf

About Bond graphs
-The System Modeling World
https://groups.csail.mit.edu/drl/journal_club/papers/Samantaray__2001__www.bondgraphs.com_about.pdf


Bond Graph Modelling of Engineering Systems
https://www3.math.tu-berlin.de/numerik/MoRTransPhen/materials/Breedveld/BreedveldChapter.pdf
Posted by at No comments:

Tuesday, December 11, 2018

A plan to get pitzDaily to work for what I am doing????

grr...I have trouble with the PitzDaily tutorial... maybe I should run the simplest simulation possible.... Just the inlet mesh.

Then I will add the FineMesh to the Inlet Mesh and run with the simplest boundary conditions possible....


useful links:

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

https://www.cfd-online.com/Forums/openfoam-solving/68072-pimplefoam-vs-simplefoam-vs-pisofoam-vs-icofoam.html


https://openfoamwiki.net/index.php/PotentialFoam

Posted by at No comments:

Sunday, December 9, 2018

owh....how do I run to prevent floating point exceptions???? Here is an example here....

For running the case:

https://github.com/bshambaugh/openfoam-experiments/tree/master/bkup_BS_JG_pitzDaily

Time = 28

smoothSolver:  Solving for Ux, Initial residual = 0.223626, Final residual = 0.00841321, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.123804, Final residual = 0.0119879, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0302533, Final residual = 0.00182838, No Iterations 2
#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::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::fvMatrix<double>::solve() at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? at ??:?
Floating point exception (core dumped)

 See error message:
https://www.cfd-online.com/Forums/openfoam-solving/122962-floating-point-exception-core-dumped.html

See Also::

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

Try initializing with OpenFoam:
https://www.cfd-online.com/Forums/openfoam-solving/107559-initialization-potentialfoam.html
Posted by at No comments:

KiCAD Links

https://www.youtube.com/watch?v=wuek-NyDaWo   [ KiCad - From Design to Manufacture - Part 1 ]
https://www.youtube.com/watch?v=UBIhCX_ET4M   [ KiCad - From Design to Manufacture - Part 2 ]
https://www.youtube.com/watch?v=pwZ3zbqxIo4    [ KiCad - From Design to Manufacture - Part 3 ]
https://www.youtube.com/watch?v=lYQqY5knVRs    [  https://www.youtube.com/watch?v=lYQqY5knVRs ]


Because there is not a push button with an LED, I will have to make one:
http://kicadhowto.org/MakeCompFootP.htm

How To Use KiCad to Make a Printed Circuit Board: Part 1


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


How To Use KiCad To Make a Printed Circuit Board: Part 2

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


How To Use KiCad to Make a Printed Circuit Board: Part 3

https://www.youtube.com/watch?v=mHx-fnpNc7Q



https://oshpark.com/#services   [ Print your board ]
JLCPCB 


Edit (2/15/2021):

EEVblog #244 - How To Lay Out A PCB - PSU Design Part 9
https://www.youtube.com/watch?v=2b1UdOmxVrw

KiCAD PCB Design

Shawn Hymel KiCaD PCB Tutorial (probably the best one)
https://www.youtube.com/watch?v=vaCVh2SAZY4&list=PL3bNyZYHcRSUhUXUt51W6nKvxx2ORvUQB


How to install & use KiCad FreeRouting/FreeRoute auto-router (incl. manual cleanup & DRC) in KiCad 5
https://www.youtube.com/watch?v=ESrHILARL7Y


KiCad 3.0 - Creating Power and Ground Planes In KiCad
https://www.youtube.com/watch?v=WcdJ7FAmD7k


EEVblog #975 - Human vs Autorouter
https://www.youtube.com/watch?v=sffuvnGhano


file:///usr/share/doc/kicad/help/en/pcbnew.html#_adjustment_of_the_zoom_level


Posted by at No comments:

Matlab Drone Simulation and Control

https://www.youtube.com/playlist?list=PLn8PRpmsu08oOLBVYYIwwN_nvuyUqEjrj
Posted by at No comments:

Saturday, December 8, 2018

boundary conditions for pitzDaily (use as template)

inlet : boundary : patch
inlet : p : zeroGradient
inlet : U : fixedValue : uniform (10 0 0)
inlet : k : fixedvalue : uniform 0.375
inlet : nut : calculated : uniform 0
inlet : omega : fixedValue : uniform 440.15
inlet : nuTilda : fixedValue : uniform 0
inlet : v2 : zeroGradient
inlet : epsilon : fixedValue : uniform 14.855

outlet : boundary : patch
outlet : p : fixedValue : uniform 0
outlet : U : zeroGradient
outlet : k : zeroGradient
outlet : nut : calculated : uniform 0
outlet : omega : zeroGradient
outlet : nuTilda : zeroGradient
outlet : v2 :  zeroGradient
outlet : epsilon : zeroGradient


upperWall : boundary : wall
upperWall : p : zeroGradient
upperWall : U : noSlip
upperWall : k : kqRWallFunction : uniform 0.375
upperWall : nut : nutkWallFunction : uniform 0
upperWall : omega : omegaWallFunction : uniform 440.15
upperWall : nuTilda : zeroGradient
upperWall : v2 : fWallFunction : uniform 0
upperWall : epsilon : epsilonWallFunction : uniform 14.855

LOOk at Josef Nagy's videos to kappa epsilon videos linked to earlier
Posted by at No comments:

Advanced Simulation Library for simulation with CPU and GPU

http://asl.org.il/

Advanced Simulation Library
Posted by at No comments:

Saturday, December 1, 2018

Steps to follow to run die case

Steps to follow:

0) Prepare the environment:
https://raptorlicious.blogspot.com/2018/08/explorations-with-openfoam.html

1) https://raptorlicious.blogspot.com/2018/12/mesh-parameters-for-holly5bkuphdf.html

2) https://raptorlicious.blogspot.com/2018/11/method-of-creation.html

3) Run using https://raptorlicious.blogspot.com/2018/11/more-ras-model-links-trying-to-find-out.html 
learning from https://raptorlicious.blogspot.com/2018/11/to-do-nov-24th.html

Edit: It looks good ( https://github.com/bshambaugh/openfoam-experiments/tree/master/B1_CoarseMesh  ) but there is an error with matching the geometry.


https://gist.github.com/bshambaugh/60ded31417283c9cedc4e2c3e4fca212

 Maybe I can see what the points are, so that I can change them. Rather than going through the effort to write a parser,

cat points | grep -o "^([0-9].* [0-9].* [0-9].*)"

cat points | grep -o "^([0-9].* [0-9].* [0-9].*)" | sort -rn

Trying the FineMesh and Inlet Mesh together with mergeMeshes, while shrinking before hand, leads to a misaligned mesh. Maybe the mesh merging should be done beforehand, and then the scaling of the meshes should follow.


Okay, I did this and then ran:

transformPoints -scale '(0.00005 0.00001 0.000008571428571428571)'


Now moving to solvers involving pitzDaily, which may or may not be for turbulent flow:

https://www.cfd-online.com/Forums/openfoam-solving/68072-pimplefoam-vs-simplefoam-vs-pisofoam-vs-icofoam.html

tutorials/incompressible/pisoFoam/les/pitzDailyMapped
tutorials/incompressible/pisoFoam/les/pitzDaily

tutorials/incompressible/simpleFoam/les/pitzDaily
tutorials/incompressible/simpleFoam/les/pitzDailyMapped
tutorials/incompressible/simpleFoam/les/pitzDailyExptInlet

tutorials/incompressible/pimpleFoam/pitzDaily
tutorials/incompressible/adjointShapeOptimizationFoam/pitzDaily/

tutorials/compressible/rhoPimpleFoam/les/pitzDaily

basic/scalarTransportFoam/pitzDaily
basic/potentialFoam/pitzDaily

Which one should I chose?? I would like a turbulent flow solver.

I like  this pose from the simplefoam-vs-pisofoam-vs-icofoam thread above: "Hi
icoFoam is transient solver for incompressible, laminar flow of Newtonian fluids.

pimpleFoam is large time-step transient solver for incompressible, flow using the PIMPLE(merged PISO-SIMPLE) algorithm.

pisoFoam is transient solver for incompressible flow.
Turbulence modelling is generic, i.e. laminar, RAS or LES may be selected.

simpleFoam is steady-state solver for compressible, turbulent flow

Regards

Ata"

Do I want the flow to be compressible or incompressible?
https://www.cfd-online.com/Forums/main/8344-incompressible-compressible-flow.html . I found from this link that the density varies.

Follow Josef Nagy's tutorial, and modify it:
https://www.youtube.com/watch?v=-46pgYQYER8&list=PLcOe4WUSsMkGPdwCpZfKcpn7w-EkgAMB8

Edit:: For convergence:

Change the Solver
rhoPimpleFoam of compressible fluids for HVAC and similar applications, with optional
mesh motion and mesh topology change
rhoSimpleFoam Steady-state solver for turbulent flow of compressible fluids.
rhoCentralFoam Density-based compressible flow solver based on central-upwind schemes of Kurganov and Tadmor.

Plot Residuals
Alter Residuals
See albcem's excellent comment at:
https://www.cfd-online.com/Forums/openfoam-solving/58365-simplefoam-convergence-problems.html
Change the grid size
Change the time steps

 Change Model
kEpsilon Standard k-epsilon turbulence model for incompressible flows.
LRR Launder, Reece and Rodi Reynolds-stress turbulence model for compressible flows.
kEpsilon Standard k-epsilon turbulence model for compressible flows including rapid distor-
tion theory (RDT) based compression term.

Posted by at No comments:

mesh parameters for holly5_bkup.hdf

 This is for FineMesh.unv , InletMesh.unv , CoarseMesh.unv :
================================================


ideasUnvToFoam FineMesh.unv

...

units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15

...

Constructing mesh with non-default patches of size:
    fineSymmetryWall    80
    fineWalls    160
    fineCyclicBoundary    160
    fineplug    80
    faceFine    13120

...

transformPoints -scale '(0.00005 0.00001 0.00000857142857)'


gives    max( X Y Z) === >  (0.001  0.006  0.03)

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

ideasUnvToFoam CourseMesh.unv

...

unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15
....

Constructing mesh with non-default patches of size:
    courseFace    40000
    courseCyclicBoundary    200
    courseWalls    100
    freeBoundaryCourse    200
    courseSymmetryWalls    100
....

transformPoints -scale '(0.00005 0.00001 0.000009523809523809525)'

gives   max(X Y Z) ===>  (0.001 0.03 0.1)

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

 ideasUnvToFoam InletMesh.unv 

...

units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15

...

Constructing mesh with non-default patches of size:
    inletFace    450
    inlet    15
    inletWalls    30
    outletInlet    15

...

 transformPoints -scale '(0.00005 0.000005598092074237316 0.00001)'

gives max(X Y Z) ===> (0.001 0.002886751346 0.005)
Posted by at No comments:
Subscribe to: Posts (Atom)