mlfem_nac issueshttps://gitlab.tue.nl/stem/mlfem_nac/-/issues2020-10-19T17:27:23Zhttps://gitlab.tue.nl/stem/mlfem_nac/-/issues/8calculate flux in femlin_cd / elcd2020-10-19T17:27:23ZTurnhout, M.C. vancalculate flux in femlin_cd / elcdfemlin_cd and elcd do not calculate flux/sigma (as e.g. elm1d & ele).
-> write elcd_d (as e.g. elm1d_d & ele_d)femlin_cd and elcd do not calculate flux/sigma (as e.g. elm1d & ele).
-> write elcd_d (as e.g. elm1d_d & ele_d)https://gitlab.tue.nl/stem/mlfem_nac/-/issues/7Column with shape functions2020-09-19T15:07:10ZTurnhout, M.C. vanColumn with shape functionsThe column N in the book is a row n in mlfem, and this seems clumsy.
But `fixing' this might involve a lot of functions...The column N in the book is a row n in mlfem, and this seems clumsy.
But `fixing' this might involve a lot of functions...https://gitlab.tue.nl/stem/mlfem_nac/-/issues/6The element 'mass' matrix of elme2020-08-19T11:21:47ZTurnhout, M.C. vanThe element 'mass' matrix of elmeThe element elme (also) produces the element `mass' matrix, but there does not seem to be a (linear solid ) engine that does anything with that mass matrix.The element elme (also) produces the element `mass' matrix, but there does not seem to be a (linear solid ) engine that does anything with that mass matrix.https://gitlab.tue.nl/stem/mlfem_nac/-/issues/5steady-state fluids: femlin_e or femnl?2021-03-30T13:36:24ZTurnhout, M.C. vansteady-state fluids: femlin_e or femnl?Some demos use rho > 0 and femnl for steady-state fluid problems, some use rho = 0 and femnl, and some use femlin_e.
E.g. demo_2df_axi_Poiseuille used rho > 0 & femnl, but these results differ from rho = 0 & femnl and femlin_e. rho = 0 and femlin_e appear to give the correct (analytical) solution for the outflow.
Which engine do you use, when and how?
What is the added value for rho > 0 and/or femnl??
It further turns out that results with femnl/elup only match femlin_e for nincr = 1, and depend on nincr otherwise (see also issue #3)Some demos use rho > 0 and femnl for steady-state fluid problems, some use rho = 0 and femnl, and some use femlin_e.
E.g. demo_2df_axi_Poiseuille used rho > 0 & femnl, but these results differ from rho = 0 & femnl and femlin_e. rho = 0 and femlin_e appear to give the correct (analytical) solution for the outflow.
Which engine do you use, when and how?
What is the added value for rho > 0 and/or femnl??
It further turns out that results with femnl/elup only match femlin_e for nincr = 1, and depend on nincr otherwise (see also issue #3)https://gitlab.tue.nl/stem/mlfem_nac/-/issues/4check implementation of time-dependent boundary conditions2020-08-03T19:01:49ZTurnhout, M.C. vancheck implementation of time-dependent boundary conditionsTime-dependent boundary conditions (`pulsatile flow') are now hard-coded in femnlt, and the implementation looks ... weird.
Validate and move to a decicated function (adaptable in /userfiles or work-directory)Time-dependent boundary conditions (`pulsatile flow') are now hard-coded in femnlt, and the implementation looks ... weird.
Validate and move to a decicated function (adaptable in /userfiles or work-directory)https://gitlab.tue.nl/stem/mlfem_nac/-/issues/3check femnl and/or cneo-elements2020-08-18T16:01:20ZTurnhout, M.C. vancheck femnl and/or cneo-elementsObservations from demo_2ds_barcneo_bending.m:
- the total Lagrange element elcneotot does not converge in the demo
- results with the `normal' updated Lagrange element elcneo depend on the number of iterations for femnl (?!)
- a sensible comparison with the demo_2ds_bar_bending seems impossible...
- results with elup only match femlin_e for nincr = 1, and (also) depend on nincr otherwiseObservations from demo_2ds_barcneo_bending.m:
- the total Lagrange element elcneotot does not converge in the demo
- results with the `normal' updated Lagrange element elcneo depend on the number of iterations for femnl (?!)
- a sensible comparison with the demo_2ds_bar_bending seems impossible...
- results with elup only match femlin_e for nincr = 1, and (also) depend on nincr otherwisehttps://gitlab.tue.nl/stem/mlfem_nac/-/issues/2femnlt & istat, mat.mat2020-08-01T14:32:15ZTurnhout, M.C. vanfemnlt & istat, mat.matFor femnlt, we store dt and theta (for istat == 2) in mat.mat, but not ntime.
Check what is going on and whether that makes (enough) sense.For femnlt, we store dt and theta (for istat == 2) in mat.mat, but not ntime.
Check what is going on and whether that makes (enough) sense.https://gitlab.tue.nl/stem/mlfem_nac/-/issues/1Check consistency of comments/help in element (support) functions2021-03-30T16:52:00ZTurnhout, M.C. vanCheck consistency of comments/help in element (support) functionsThe element library contains al lot of remnants from ad-hoc copy-paste actions.
Matlab help and code comment is inconsistent and sometimes just wrong (in particular on mat.mat, in particular on mat.mat(11))The element library contains al lot of remnants from ad-hoc copy-paste actions.
Matlab help and code comment is inconsistent and sometimes just wrong (in particular on mat.mat, in particular on mat.mat(11))