OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Component lable delphi<1e-4

Browse Source
This commit is contained in:
James E McClure
2015-08-08 21:02:19 -04:00
parent 8c3ec59e92
commit fca04fe4f3
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
common/TwoPhase.h
View File

@@ -634,7 +634,7 @@ void TwoPhase::ComponentAverages(){
ComponentAverages_NWP(VSQ,LabelNWP) += 0.125*(Vel_x(n)*Vel_x(n)+Vel_y(n)*Vel_y(n)+Vel_z(n)*Vel_z(n));
// volume the for pressure averaging excludes the interfacial region
if (DelPhi(n) < 1e-2 ){
if (DelPhi(n) < 1e-4 ){
ComponentAverages_NWP(TRIMVOL,LabelNWP) += 0.125;
ComponentAverages_NWP(PRS,LabelNWP) += 0.125*Press(n);
}
@@ -649,7 +649,7 @@ void TwoPhase::ComponentAverages(){
ComponentAverages_WP(VSQ,LabelWP) += 0.125*(Vel_x(n)*Vel_x(n)+Vel_y(n)*Vel_y(n)+Vel_z(n)*Vel_z(n));
// volume the for pressure averaging excludes the interfacial region
if (DelPhi(n) < 1e-2){
if (DelPhi(n) < 1e-4){
ComponentAverages_WP(TRIMVOL,LabelWP) += 0.125;
ComponentAverages_WP(PRS,LabelWP) += 0.125*Press(n);
}