====== Commit 2019-09-24 ====== ===== Description ===== **Lobatto integration on boundary elements:** In commit [[commit:2019:01_17|]] I had interfaced the Lobatto Integration for volume Elements, this commit extends it to boundary condition elements. **Fix inconsistency in some of the boundary elements:** In some of the boundary elements i.e. Tm2Convection3DElement.cpp, the function "getNumberOfTherNod()" returned the total number of nodes instead of the number of thermal nodes. **Bug fix:** OnEdgeGpHexa::fillLobattoPoints was mistakenly using SimpsonIntegration instead of Lobatto integration. ===== Tests Lobatto ===== I added some AM tests comparing the use of lobatto integration for volume elements, boundary elements, and both. {{:commit:2019:09_24:complobatto.png?600|}} One can see that using the Lobatto integration on the boundary elements has little influence on the results while using it on the volume element has a great impact, at least in simulations with coarse meshes. {{:commit:2019:09_24:geomlobatto.jpg?600|}} ===== Addition of battery tests ===== I added battery tests in monosThermoMeca to test lobatto integration on Volume/boundary/both. I added battery tests for AM using lobatto integration on Volume/boundary/both. ===== Added [a] / deleted [d] / modified [m] / renamed [m] files ====== ====Code:==== [m] mtElements/boundaries/Tm2Convection2DElement.cpp [m] mtElements/boundaries/Tm2Convection3DElement.cpp [m] mtElements/boundaries/Tm2ConvectionElement.h [m] mtElements/boundaries/Tm2ConvectionElement.inl [m] mtElements/boundaries/Tm2ConvectionElShcuts.h [m] mtElements/boundaries/Tm2ConvectionElShcuts.inl [m] mtElements/boundaries/Tm2Flux3DElement.cpp [m] mtElements/boundaries/Tm2FluxElement.h [m] mtElements/boundaries/Tm2FluxElement.inl [m] mtElements/boundaries/Tm2FluxElShcuts.h [m] mtElements/boundaries/Tm2FluxElShcuts.inl [m] mtElements/boundaries/Tm2Rayonnement2DElement.cpp [m] mtElements/boundaries/Tm2Rayonnement3DElement.cpp [m] mtElements/boundaries/Tm2RayonnementElement.h [m] mtElements/boundaries/Tm2RayonnementElement.inl [m] mtElements/boundaries/Tm2RayonnementElShcuts.h [m] mtElements/boundaries/Tm2RayonnementElShcuts.inl [m] mtElements/boundaries/Tm2VolFlux2DElement.cpp [m] mtElements/boundaries/Tm2VolFlux3DElement.cpp [m] mtElements/boundaries/TmConvection2DElement.cpp [m] mtElements/boundaries/TmConvection3DElement.cpp [m] mtElements/boundaries/TmConvectionElement.h [m] mtElements/boundaries/TmConvectionElement.inl [m] mtElements/boundaries/TmConvectionElShcuts.h [m] mtElements/boundaries/TmConvectionElShcuts.inl [m] mtElements/boundaries/TmFlux2DElement.cpp [m] mtElements/boundaries/TmFlux3DElement.cpp [m] mtElements/boundaries/TmFluxElement.h [m] mtElements/boundaries/TmFluxElement.inl [m] mtElements/boundaries/TmFluxElShcuts.h [m] mtElements/boundaries/TmFluxElShcuts.inl [m] mtElements/boundaries/TmRayonnement2DElement.cpp [m] mtElements/boundaries/TmRayonnement3DElement.cpp [m] mtElements/boundaries/TmRayonnementElement.h [m] mtElements/boundaries/TmRayonnementElement.inl [m] mtElements/boundaries/TmRayonnementElShcuts.h [m] mtElements/boundaries/TmRayonnementElShcuts.inl [m] mtElements/boundaries/TmVolFlux2DElement.cpp [m] mtElements/boundaries/TmVolFlux3DElement.cpp [m] mtElements/boundaries/Traction2DElement.cpp [m] mtElements/boundaries/Traction3DElement.cpp [m] mtElements/boundaries/TractionElement.h [m] mtElements/boundaries/TractionElement.inl [m] mtElements/boundaries/TractionElShcuts.h [m] mtElements/boundaries/TractionElShcuts.inl [m] mtElements/boundaries/TriangleTraction3DElement.cpp [m] mtElements/boundaries/TriangleTractionElement.h [m] mtElements/boundaries/TriangleTractionElement.inl [m] mtElements/boundaries/TriangleTractionElShcuts.h [m] mtElements/boundaries/TriangleTractionElShcuts.inl [m] mtShapeFunctions/OnEdgeGpHexa.cpp ====Tests:==== [a] additiveM/battery/withAM/Chiumenti2LobattoBc.py [a] additiveM/battery/withAM/Chiumenti2LobattoFull.py [a] additiveM/battery/withAM/ChiumentiLobattoBc.py [a] additiveM/battery/withAM/ChiumentiLobattoFull.py [a] apps/monosThermoMeca/convection3dLobattoBc.py [a] apps/monosThermoMeca/convection3dLobattoFull.py [a] apps/monosThermoMeca/convection3dLobattoVol.py [a] apps/monosThermoMeca/flux3dLobattoBc.py [a] apps/monosThermoMeca/flux3dLobattoFull.py [a] apps/monosThermoMeca/flux3dLobattoVol.py [r] additiveM/battery/withAM/Chiumenti2LobattoVol.py (previously "Chiumenti2Lobatto") [r] additiveM/battery/withAM/ChiumentiLobattoVol.py (previously "ChiumentiLobatto") --- //[[Cedric.Laruelle@ULiege.be|Cédric Laruelle]] 2019/09/24//