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:

Friday, November 30, 2018

The earth rings

https://www.nationalgeographic.com/science/2018/11/strange-earthquake-waves-rippled-around-world-earth-geology/
Posted by at No comments:

More Scale the Mesh 3


After: transformPoints -scale '(0.001 0.001 0.001)'    (Scale by mm)

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(0.02 0 0)
(0 0 0.2)
(0.02 0 0.2)
(0 0.2 0)
(0.02 0.2 0)
(0 0.2 0.2)
(0.02 0.2 0.2)
(0 0 0.01)
(0 0 0.02)
(0 0 0.03)
(0 0 0.04)
(0 0 0.05)
(0 0 0.06)
(0 0 0.07)
(0 0 0.08)
(0 0 0.09)
(0 0 0.1)
(0 0 0.11)
(0 0 0.12)
(0 0 0.13)
(0 0 0.14)
(0 0 0.15)
(0 0 0.16)
(0 0 0.17)
(0 0 0.18)
(0 0 0.19)


-----------------------

Before:


/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(20 0 0)
(0 0 200)
(20 0 200)
(0 200 0)
(20 200 0)
(0 200 200)
(20 200 200)
(0 0 10)
(0 0 20)
(0 0 30)
(0 0 40)
(0 0 50)
(0 0 60)
(0 0 70)
(0 0 80)
(0 0 90)
(0 0 100)
(0 0 110)
(0 0 120)
(0 0 130)
(0 0 140)
(0 0 150)
(0 0 160)
(0 0 170)
(0 0 180)
(0 0 190)

-------------------------

transformPoints -scale '(0.01 0.01 0.01)'      scale by cm
Posted by at No comments:

Thursday, November 29, 2018

More Scale the Mesh 2

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.1;

vertices
(
    (0 0 0)
    (1 0 0)
    (1 1 0)
    (0 1 0)
    (0 0 0.1)
    (1 0 0.1)
    (1 1 0.1)
    (0 1 0.1)
);
-----------------------------------------------------------------------------------

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(0.005 0 0)
(0.01 0 0)
(0.015 0 0)
(0.02 0 0)
(0.025 0 0)
(0.03 0 0)
(0.035 0 0)
(0.04 0 0)
(0.045 0 0)
(0.05 0 0)
(0.055 0 0)
(0.06 0 0)
(0.065 0 0)
(0.07 0 0)
(0.075 0 0)
(0.08 0 0)
(0.085 0 0)
(0.09 0 0)
(0.095 0 0)
(0.1 0 0)
----------------------------------------------------------------------------------------------------

The step size is 0.1 / 20 or the "convertToMeters 0.1" divided by the number of patches "20".

What number of patches would one use for an irregular shape? Would it be based on the number of patches in the shortest face?

----------------------------------------------------------------------------------------------------

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

3
(
    movingWall
    {
        type            wall;
        inGroups        1(wall);
        nFaces          20;
        startFace       760;
    }
    fixedWalls
    {
        type            wall;
        inGroups        1(wall);
        nFaces          60;
        startFace       780;
    }
    frontAndBack
    {
        type            empty;
        inGroups        1(empty);
        nFaces          800;
        startFace       840;
    }
)
Posted by at No comments:

More Scale the MEsh

incompressible/icoFoam/Mesh_1
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Processing tag:164
Starting reading units at line 3.
l:1
units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15


Processing tag:2420
Skipping tag 2420 on line 9
Skipping section at line 9.

Processing tag:2411
Starting reading points at line 20.
Read 882 points.

Processing tag:2412
Starting reading cells at line 1787.
First occurrence of element type 11 for cell 1 at line 1788
First occurrence of element type 44 for cell 85 at line 2280
First occurrence of element type 115 for cell 1045 at line 4040
Read 400 cells and 880 boundary faces.

Processing tag:2467
Starting reading patches at line 4842.
For group 1 named sourceFaceleft trying to read 400 patch face indices.
For group 2 named destinationFaceleft trying to read 400 patch face indices.
For group 3 named channelWallsleft trying to read 40 patch face indices.
For group 4 named leftRightWallleft trying to read 20 patch face indices.
For group 5 named inlet trying to read 20 patch face indices.

Sorting boundary faces according to group (patch)
0: sourceFaceleft is patch
1: destinationFaceleft is patch
2: channelWallsleft is patch
3: leftRightWallleft is patch
4: inlet is patch

Constructing mesh with non-default patches of size:
    sourceFaceleft    400
    destinationFaceleft    400
    channelWallsleft    40
    leftRightWallleft    20
    inlet    20

End
---------------------------------------------
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


(0 0 0)
(0 0 10)
(0 0 20)
(0 0 30)
(0 0 40)
(0 0 50)
(0 0 60)
(0 0 70)
(0 0 80)
(0 0 90)
(0 0 100)
(0 0 110)
(0 0 120)
(0 0 130)
(0 0 140)
(0 0 150)
(0 0 160)
(0 0 170)
(0 0 180)
(0 0 190)
(0 0 200)



Posted by at No comments:

Saturday, November 24, 2018

Scale the Mesh

Scale the mesh:




By looking at:

https://www.cfd-online.com/Forums/openfoam-meshing/198268-how-scale-stl-when-importing-salome-any-cad-tool-openfoam.html

try:

transformPoints -scale '(0.001 0.001 0.001)'

https://cfd.direct/openfoam/user-guide/v6-basic-file-format/

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

What is the default mesh size???

http://www.salome-platform.org/forum/forum_10/549857363

-------------------------------------------------------------------------------------------------------
  -1
   164
         1  SI: Meter (newton)         2
    1.0000000000000000E+0    1.0000000000000000E+0    1.0000000000000000E+0
    2.7314999999999998E+2
    -1
    -1
  2420
         1
SMESH_Mesh
         1         0         0
Global Cartesian Coordinate System
    1.0000000000000000E+0    0.0000000000000000E+0    0.0000000000000000E+0
    0.0000000000000000E+0    1.0000000000000000E+0    0.0000000000000000E+0
    0.0000000000000000E+0    0.0000000000000000E+0    1.0000000000000000E+0
    0.0000000000000000E+0    0.0000000000000000E+0    0.0000000000000000E+0
    -1
    -1
  2411
         1         1         1        11
   2.0000000000000000E+01   4.5500000000000000E+02   0.0000000000000000E+00
         2         1         1        11
   0.0000000000000000E+00   4.5500000000000000E+02   0.0000000000000000E+00
         3         1         1        11
   2.0000000000000000E+01   5.2000000000000000E+02   0.0000000000000000E+00
         4         1         1        11
   0.0000000000000000E+00   5.2000000000000000E+02   0.0000000000000000E+00
         5         1         1        11
   2.0000000000000000E+01   1.6600000000000000E+02   5.0000000000000000E+02

-----------------------------------------------------------------------------------------------------
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : ideasUnvToFoam InletMesh.unv
Date   : Nov 26 2018
Time   : 18:32:30
Host   : "caelinux"
PID    : 22307
Case   : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/test-InletMesh
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Processing tag:164
Starting reading units at line 3.
l:1
units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15


-------------------------------------------------------------------------------
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(0.005 0 0)
(0.01 0 0)
(0.015 0 0)
(0.02 0 0)
(0.025 0 0)
(0.03 0 0)
(0.035 0 0)
(0.04 0 0)
(0.045 0 0)
(0.05 0 0)
(0.055 0 0)
(0.06 0 0)
(0.065 0 0)
(0.07 0 0)
(0.075 0 0)
(0.08 0 0)
(0.085 0 0)
(0.09 0 0)
(0.095 0 0)
(0.1 0 0)

---------------------------------

See: https://cfd.direct/openfoam/user-guide/v6-cavity/   ... length is 0.1 meters

Posted by at No comments:
Subscribe to: Posts (Atom)