mhamadmahdialloush / uFVM

Learning the Finite Volume Method in CFD with MATLAB Programming
https://www.aub.edu.lb/msfea/research/Pages/cfd.aspx
139 stars 44 forks source link

Non-orthogonal part and velocity gradient transpose part at boundary while assembling stress term #11

Open zhanggh97 opened 3 years ago

zhanggh97 commented 3 years ago

Hi, I have a question about assembing stress term in uFVM. When we assemble stress term at interior face,we need to add non-orthogonal part and velocity gradient transpose part to FluxVf, which is as followed: % Non-orthogonal part local_FluxVf = - mu_f.*dot(UGrad_f(:,:,iComponent)',Tf')'; % Add transpose term local_FluxVf = local_FluxVf + mu_f.*dot(tUGrad_f(:,:,iComponent)',Sf')';

But we didn't do the same thing at boundary: % non-linear fluxes local_FluxVb = zeros(size(local_FluxFb)); Why?

mhamadmahdialloush commented 3 years ago

The non-ortho stress term at boundaries is negligible, we skip it usually.

On Mon, Nov 8, 2021 at 8:45 AM zhanggh97 @.***> wrote:

Hi, I have a question about assembing stress term in uFVM. When we assemble stress term at interior face,we need to add non-orthogonal part and velocity gradient transpose part to FluxVf, which is as followed: % Non-orthogonal part local_FluxVf = - mu_f.dot(UGrad_f(:,:,iComponent)',Tf')'; % Add transpose term local_FluxVf = local_FluxVf + mu_f.dot(tUGrad_f(:,:,iComponent)',Sf')';

But we didn't do the same thing at boundary: % non-linear fluxes local_FluxVb = zeros(size(local_FluxFb)); Why?

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/mhamadAlloush/uFVM/issues/11, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFLISFM6CVXCU6G2PJUKN6LUK552HANCNFSM5HR63RRQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

zhanggh97 commented 3 years ago

How about velocity gradient transpose part at boundary?

mhamadmahdialloush commented 3 years ago

Also not of a significant order

On Mon, Nov 8, 2021 at 11:19 AM zhanggh97 @.***> wrote:

How about velocity gradient transpose part at boundary?

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/mhamadAlloush/uFVM/issues/11#issuecomment-963006076, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFLISFPGBDCUQCLEBSJ7M73UK6P4HANCNFSM5HR63RRQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

zhanggh97 commented 3 years ago

Got it, thanks