


REDUCEVECTORSTOKES - compute a reduced vector to get rid of the bubble grid
This function return the reduced elementary load vector, so that we ca
solve the problem without having a grid in the middle of the pentaelem
[ ] [ ]
[ ] [ ]
Kebubble=[ Kib ] and Pe_temp=[ Pi ]
[ ] [ ]
[-----] [-------]
[ Kbb ] [ Pb ]
indice b represents the bubble function (node in the middle of the penta
indice i represents the other nodes (the six principal nodes of the penta)
So that the elementary reduced load vector is pe_reduced=Pi-Kib*Kbb^-1*Pb
See theoryguide for more informations
Usage:
pe_reduced=ReduceVectorStokes(pentaelem,Kebubble_temp,pe_temp)
See also REDUCEMATRIXSTOKES, CREATEPVECTOR

0001 function pe_reduced=ReduceVectorStokes(pentaelem,Kebubble_temp,pe_temp) 0002 %REDUCEVECTORSTOKES - compute a reduced vector to get rid of the bubble grid 0003 % 0004 % This function return the reduced elementary load vector, so that we ca 0005 % solve the problem without having a grid in the middle of the pentaelem 0006 % [ ] [ ] 0007 % [ ] [ ] 0008 % Kebubble=[ Kib ] and Pe_temp=[ Pi ] 0009 % [ ] [ ] 0010 % [-----] [-------] 0011 % [ Kbb ] [ Pb ] 0012 % indice b represents the bubble function (node in the middle of the penta 0013 % indice i represents the other nodes (the six principal nodes of the penta) 0014 % So that the elementary reduced load vector is pe_reduced=Pi-Kib*Kbb^-1*Pb 0015 % See theoryguide for more informations 0016 % 0017 % Usage: 0018 % pe_reduced=ReduceVectorStokes(pentaelem,Kebubble_temp,pe_temp) 0019 % 0020 % See also REDUCEMATRIXSTOKES, CREATEPVECTOR 0021 0022 %Create load vector whithout using the bubble function 0023 Kib=Kebubble_temp(1:24,:); 0024 Kbb=Kebubble_temp(25:27,:); 0025 0026 Pi=pe_temp(1:24,:); 0027 Pb=pe_temp(25:27,:); 0028 0029 pe_reduced=Pi-Kib*Kbb^-1*Pb;