Skip to content

Commit

Permalink
Surface Flux, Const. Massflow (AdaptiveType=4): Correct weighting bas…
Browse files Browse the repository at this point in the history 
…ed on the volume flow rate
  • Loading branch information
@pnizenkov
pnizenkov committed Dec 5, 2022
1 parent 4ee759d commit 6882692
Show file tree
Hide file tree
Showing 5 changed files with 59 additions and 96 deletions.
80 changes: 30 additions & 50 deletions regressioncheck/WEK_DSMC/ChannelFlow_AdaptiveBoundary_ConstMassflow/PartAnalyze_ref.csv
View file
Original file line number Diff line number Diff line change
@@ -1,52 +1,32 @@
001-TIME,002-Massflow-Spec-001-SF-001,003-Pressure-Spec-001-SF-001,004-Massflow-Spec-001-SF-002,005-Pressure-Spec-001-SF-002
0.0000000000000000E+000,0.0000000000000000E+000,0.0000000000000000E+000,0.0000000000000000E+000,0.0000000000000000E+000
0.1000000000000000E-005,-.1974504031999998E-014,0.5101742792810926E+000,-.4952200640000046E-015,0.2542441898736914E+001
0.2000000000000000E-005,0.6335841247999995E-014,0.1291460470962966E+001,0.4804484784000007E-014,0.2503815298384615E+001
0.3000000000000000E-005,0.1270993984000001E-013,0.1597378707630111E+001,0.7324155967999996E-014,0.2487501534161306E+001
0.4000000000000000E-005,0.1620092248000001E-013,0.1829005527225455E+001,0.5510120240000000E-014,0.2470963365440040E+001
0.5000000000000000E-005,0.1939115988799999E-013,0.2004734853455359E+001,0.3579187072000004E-014,0.2475567562968181E+001
0.6000000000000001E-005,0.2113505715200000E-013,0.2334225246046347E+001,0.7120116799999953E-015,0.2468310809669783E+001
0.7000000000000002E-005,0.2172060705600000E-013,0.2621163565453777E+001,-.8586648319999970E-015,0.2471276432351922E+001
0.8000000000000001E-005,0.2461009923200001E-013,0.2854472800398014E+001,-.4278446304000005E-014,0.2483880659976196E+001
0.9000000000000002E-005,0.2519139831999998E-013,0.3054888687339466E+001,-.6989404207999997E-014,0.2505084554352881E+001
0.1000000000000000E-004,0.2494697639999999E-013,0.3243656393703143E+001,-.8843822687999999E-014,0.2512998912435173E+001
0.1100000000000000E-004,0.2566961512000000E-013,0.3400256975212518E+001,-.1026465793600000E-013,0.2528877222346422E+001
0.1200000000000000E-004,0.2611382539199998E-013,0.3542974078331929E+001,-.1324873076800001E-013,0.2534329297962050E+001
0.1300000000000000E-004,0.2852935158399998E-013,0.3692760963488928E+001,-.1606383366400001E-013,0.2534805308582880E+001
0.1400000000000000E-004,0.2845389960000001E-013,0.3833606765181604E+001,-.1748679432000001E-013,0.2541050397792644E+001
0.1500000000000001E-004,0.2856335811200004E-013,0.3952073704962856E+001,-.1825831742400000E-013,0.2542386532895929E+001
0.1600000000000000E-004,0.2954954742400001E-013,0.4056500091185906E+001,-.2075142100800000E-013,0.2546735137098510E+001
0.1700000000000000E-004,0.3060906331200000E-013,0.4172502615535060E+001,-.2171848164800000E-013,0.2555341398416259E+001
0.1800000000000000E-004,0.3044753230400001E-013,0.4248860828199861E+001,-.2319138939200000E-013,0.2555577958645423E+001
0.1899999999999999E-004,0.3108302929599998E-013,0.4316843991001847E+001,-.2368660945600000E-013,0.2543009040842799E+001
0.1999999999999999E-004,0.3158037476800001E-013,0.4399639988275036E+001,-.2502349108799999E-013,0.2542432078855079E+001
0.2099999999999999E-004,0.3089599339200001E-013,0.4480677699463522E+001,-.2526153678400001E-013,0.2531285348388810E+001
0.2199999999999999E-004,0.3160481695999998E-013,0.4539417552997691E+001,-.2708513684800000E-013,0.2519695077855990E+001
0.2299999999999998E-004,0.3215317222400000E-013,0.4605838469402263E+001,-.2769300353599999E-013,0.2514827387846171E+001
0.2399999999999998E-004,0.3345179651200004E-013,0.4657787191766029E+001,-.2897568726400000E-013,0.2522566007539098E+001
0.2500000000000000E-004,0.3358675991999849E-013,0.4702201070712960E+001,-.2811595972799875E-013,0.2521226242151291E+001
0.2600000000000000E-004,0.3344435758400001E-013,0.4737025488934150E+001,-.2891511313600000E-013,0.2525362602392541E+001
0.2700000000000000E-004,0.3251767969599998E-013,0.4773666162726126E+001,-.2957824043200000E-013,0.2526377647463425E+001
0.2799999999999999E-004,0.3406497672000001E-013,0.4799012181701939E+001,-.3048897776000001E-013,0.2530070078036206E+001
0.2899999999999999E-004,0.3243053796800000E-013,0.4833292514342970E+001,-.3138164912000001E-013,0.2531712279864621E+001
0.2999999999999999E-004,0.3293851048000004E-013,0.4853118254077882E+001,-.3110215796800001E-013,0.2536822428623428E+001
0.3099999999999999E-004,0.3404159723199999E-013,0.4890066723804255E+001,-.3238165358400000E-013,0.2545403574506630E+001
0.3199999999999998E-004,0.3375360444800000E-013,0.4898869745732007E+001,-.3117229643200000E-013,0.2550700924993627E+001
0.3299999999999998E-004,0.3380461423999999E-013,0.4929836993123001E+001,-.3313298531200000E-013,0.2558883387327509E+001
0.3399999999999998E-004,0.3396614524799998E-013,0.4958391882818118E+001,-.3343479324800000E-013,0.2543203058492540E+001
0.3499999999999998E-004,0.3498209027200003E-013,0.4970918732750711E+001,-.3257400300800001E-013,0.2542691294072281E+001
0.3599999999999997E-004,0.3440716740800003E-013,0.4980224573400316E+001,-.3298208134400001E-013,0.2526933389266032E+001
0.3699999999999997E-004,0.3439122684799996E-013,0.5000689291186351E+001,-.3357719558400000E-013,0.2525941194659184E+001
0.3799999999999997E-004,0.3545393084800000E-013,0.5004695893758515E+001,-.3316274102400001E-013,0.2514265108821918E+001
0.3899999999999997E-004,0.3309579067200000E-013,0.5010430661751538E+001,-.3319143403200000E-013,0.2517133827088020E+001
0.3999999999999996E-004,0.3398208580800000E-013,0.5034034268675430E+001,-.3406391401600000E-013,0.2515266646462496E+001
0.4099999999999996E-004,0.3374616552000000E-013,0.5052493362512732E+001,-.3382374291200000E-013,0.2514718077625061E+001
0.4199999999999996E-004,0.3299058297599997E-013,0.5061493186021456E+001,-.3283861630399999E-013,0.2515582798087293E+001
0.4299999999999995E-004,0.3324882004799997E-013,0.5081563781841816E+001,-.3448580750400001E-013,0.2509421205593210E+001
0.4399999999999995E-004,0.3368877950400001E-013,0.5076224818850775E+001,-.3358463451200000E-013,0.2511144614887840E+001
0.4499999999999995E-004,0.3331789580800002E-013,0.5077865736917214E+001,-.3481205763200000E-013,0.2497934464939569E+001
0.4599999999999995E-004,0.3470153641600000E-013,0.5067855302804303E+001,-.3403628371200000E-013,0.2499355259397003E+001
0.4699999999999994E-004,0.3622757935999999E-013,0.5079292771797843E+001,-.3583544158400000E-013,0.2495381069439121E+001
0.4799999999999994E-004,0.3432427649600003E-013,0.5083461994026425E+001,-.3579187072000000E-013,0.2502324728408006E+001
0.4899999999999994E-004,0.3704373603200001E-013,0.5094953026699098E+001,-.3547837304000001E-013,0.2501689704129798E+001
0.5000000000000000E-004,0.3505966766399535E-013,0.5101040982105030E+001,-.3431258675199544E-013,0.2508849092952619E+001
0.1000000000000000E-005,-.2160477232000007E-014,0.5438100504172863E+000,0.1301281047999999E-013,0.2551613014067587E+001
0.2000000000000000E-005,0.4997896912000020E-014,0.1294771852625714E+001,0.1821793467200001E-013,0.2522460339588515E+001
0.3000000000000000E-005,0.1148783024000000E-013,0.1607708488706000E+001,0.1194479295999999E-013,0.2507726067154505E+001
0.4000000000000000E-005,0.1461430540800002E-013,0.1844670173859858E+001,0.6447425167999980E-014,0.2508081298992916E+001
0.5000000000000000E-005,0.1741240503999998E-013,0.2042083618751872E+001,-.8395361599999189E-016,0.2510614047512956E+001
0.6000000000000001E-005,0.2088000819200000E-013,0.2370947583486195E+001,-.4208307840000003E-014,0.2515927055831678E+001
0.7000000000000002E-005,0.2256651944000001E-013,0.2679163296266943E+001,-.9919279135999999E-014,0.2522177171671231E+001
0.8000000000000001E-005,0.2453783536000000E-013,0.2937511632231268E+001,-.1417647136000000E-013,0.2538381520496905E+001
0.9000000000000002E-005,0.2645176526400002E-013,0.3153553886098444E+001,-.1623386630400000E-013,0.2553399419209018E+001
0.1000000000000000E-004,0.2653784428800007E-013,0.3334749837248620E+001,-.1990657132800005E-013,0.2566844253875795E+001
0.1100000000000000E-004,0.2777270633600001E-013,0.3492953624889923E+001,-.2192783433600001E-013,0.2570335420450902E+001
0.1200000000000000E-004,0.3024243043200002E-013,0.3621778685508302E+001,-.2422327497599998E-013,0.2574217960123442E+001
0.1300000000000000E-004,0.3107771577599999E-013,0.3738757438283882E+001,-.2563242047999999E-013,0.2579682290620827E+001
0.1400000000000000E-004,0.3111809852799999E-013,0.3823138284903422E+001,-.2775251496000002E-013,0.2575956093329860E+001
0.1500000000000000E-004,0.3034976353599999E-013,0.3902495665276909E+001,-.2966963297600000E-013,0.2570122033964163E+001
0.1600000000000000E-004,0.3368559139199998E-013,0.3971010192852202E+001,-.3042627822400002E-013,0.2571869474677951E+001
0.1700000000000000E-004,0.3243903960000001E-013,0.4015182312569924E+001,-.3090662043200000E-013,0.2567865329490435E+001
0.1800000000000000E-004,0.3333277366400000E-013,0.4054365668882586E+001,-.3153149038400000E-013,0.2557568567319311E+001
0.1899999999999999E-004,0.3436040843200001E-013,0.4105523487124288E+001,-.3189706056000001E-013,0.2551679261347710E+001
0.2000000000000000E-004,0.3273447131199942E-013,0.4131746825633499E+001,-.3220843283199944E-013,0.2552036487020768E+001
0.2100000000000000E-004,0.3345923544000000E-013,0.4157128781561378E+001,-.3276210161600001E-013,0.2539856945002305E+001
0.2200000000000000E-004,0.3470366182399999E-013,0.4187447617109554E+001,-.3400440259199999E-013,0.2537024662230533E+001
0.2299999999999999E-004,0.3539335672000000E-013,0.4204019117093607E+001,-.3398421121600001E-013,0.2538983258948380E+001
0.2399999999999999E-004,0.3449005831999997E-013,0.4220405331369403E+001,-.3314148694400001E-013,0.2537614088417073E+001
0.2499999999999999E-004,0.3453788000000000E-013,0.4241800629420518E+001,-.3421906880000000E-013,0.2535444492139198E+001
0.2599999999999998E-004,0.3390769652799999E-013,0.4254092336491262E+001,-.3445711449600000E-013,0.2544904136886347E+001
0.2699999999999998E-004,0.3445286368000001E-013,0.4264381532485316E+001,-.3507242011200000E-013,0.2547404515918131E+001
0.2799999999999998E-004,0.3477805110400002E-013,0.4267891937629085E+001,-.3342097809600000E-013,0.2549471886043859E+001
0.2899999999999998E-004,0.3523076300799998E-013,0.4263639231454835E+001,-.3468878396800000E-013,0.2553729581195537E+001
0.3000000000000000E-004,0.3573873551999816E-013,0.4262598101513217E+001,-.3429664619199823E-013,0.2546734480158401E+001
Binary file modified regressioncheck/WEK_DSMC/ChannelFlow_AdaptiveBoundary_ConstMassflow/cube_mesh.h5
View file
Binary file not shown.
4 changes: 2 additions & 2 deletions regressioncheck/WEK_DSMC/ChannelFlow_AdaptiveBoundary_ConstMassflow/hopr.ini
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4,8 +4,8 @@ DebugVisuLevel=2
NVisu =1
Mode =1

Corner =(/0.0,0.0,0.0 ,, 15.0,0.0,0.0 ,, 15.0,1.0,0.0 ,, 0.0,1.0,0.0 ,, 0.0,0.0,1.0 ,, 15.0,0.0,1.0 ,, 15.0,1.0,1.0 ,, 0.0,1.0,1.0 /)
nElems =(/50,4,4/)
Corner =(/0.0,0.0,0.0 ,, 10.0,0.0,0.0 ,, 10.0,1.0,0.0 ,, 0.0,1.0,0.0 ,, 0.0,0.0,1.0 ,, 10.0,0.0,1.0 ,, 10.0,1.0,1.0 ,, 0.0,1.0,1.0 /)
nElems =(/40,4,4/)
elemtype =108

BCIndex =(/6 ,4 ,1 ,3 ,2 ,5/)
Expand Down
7 changes: 4 additions & 3 deletions regressioncheck/WEK_DSMC/ChannelFlow_AdaptiveBoundary_ConstMassflow/parameter.ini
View file
Original file line number Diff line number Diff line change
Expand Up @@ -31,8 +31,8 @@ TimeStampLength = 14
! =============================================================================== !
! CALCULATION
! =============================================================================== !
tend = 5.0E-5
Analyze_dt = 2.5E-5
tend = 3.0E-5
Analyze_dt = 1.0E-5
CFLscale = 0.2 ! Scaling of theoretical CFL number
! Load balancing
DoLoadBalance = T
Expand Down Expand Up @@ -82,7 +82,7 @@ Part-FIBGMdeltas=(/2e-6,2e-6,2e-6/)
! Species1 - O2
! =============================================================================== !
Part-Species1-MassIC=5.31352E-26
Part-Species1-MacroParticleFactor=2E2
Part-Species1-MacroParticleFactor=1E2

Part-Species1-nInits=1
Part-Species1-Init1-SpaceIC=cell_local
Expand Down Expand Up @@ -123,6 +123,7 @@ Part-Species1-Surfaceflux2-Adaptive-Pressure=2.5
AdaptiveBC-RelaxationFactor = 0.1
! Sample is reset every given number of iterations
AdaptiveBC-SamplingIteration = 0,10,10,10
!AdaptiveBC-SamplingIteration = 10
! Sample is not reset, current value replaces the oldest (DEFAULT: TRUE)
AdaptiveBC-TruncateRunningAverage = F,F,T,T
! Average the sample over the whole BC
Expand Down
64 changes: 23 additions & 41 deletions src/particles/emission/particle_surface_flux.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -100,7 +100,7 @@ SUBROUTINE ParticleSurfaceflux()
#if USE_MPI
CALL MPI_ALLREDUCE(MPI_IN_PLACE,AdaptBCPartNumOut(iSpec,iSF),1,MPI_INTEGER,MPI_SUM,PartMPI%COMM,IERROR)
#endif
IF(.NOT.ALMOSTEQUAL(SF%AdaptiveMassflow,0.)) CALL AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)
IF(.NOT.ALMOSTEQUAL(SF%AdaptiveMassflow,0.)) CALL CalcConstMassflowWeight(iSpec,iSF)
END IF
!Calc Particles for insertion in standard case
IF ((.NOT.DoPoissonRounding).AND.(.NOT. DoTimeDepInflow).AND.(.NOT.RadialWeighting%DoRadialWeighting) &
Expand Down Expand Up @@ -304,6 +304,7 @@ SUBROUTINE ParticleSurfaceflux()
! instead of an UpdateNextfreePosition we update the particleVecLength only - enough ?!?
PDM%CurrentNextFreePosition = PDM%CurrentNextFreePosition + NbrOfParticle
PDM%ParticleVecLength = PDM%ParticleVecLength + NbrOfParticle
IF(PDM%ParticleVecLength.GT.PDM%maxParticleNumber) CALL abort(__STAMP__, 'Error ParticleSurfaceflux: Maximum number of particles reached!')
#if USE_LOADBALANCE
CALL LBPauseTime(LB_SURFFLUX,tLBStart)
#endif /*USE_LOADBALANCE*/
Expand Down Expand Up @@ -1015,17 +1016,16 @@ END SUBROUTINE CalcPartInsAdaptive


!===================================================================================================================================
!> Routine calculates the weights of the triangles for AdaptiveType=4 to scale up the number of particles to be inserted
!> Routine calculates the weights of the triangles for AdaptiveType=4 depending on the side-specific volume flow rate
!===================================================================================================================================
SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)
SUBROUTINE CalcConstMassflowWeight(iSpec,iSF)
! MODULES !
!----------------------------------------------------------------------------------------------------------------------------------!
USE MOD_Globals
USE MOD_Globals_Vars ,ONLY: BoltzmannConst, Pi
USE MOD_Particle_Vars ,ONLY: Species
USE MOD_Particle_Sampling_Vars ,ONLY: AdaptBCMacroVal, AdaptBCMapElemToSample, AdaptBCBackupVelocity
USE MOD_Particle_Surfaces_Vars ,ONLY: SurfMeshSubSideData, BCdata_auxSF, SurfFluxSideSize
USE MOD_TimeDisc_Vars ,ONLY: dt, RKdtFrac
USE MOD_Mesh_Vars ,ONLY: SideToElem, offsetElem
USE MOD_Particle_Mesh_Tools ,ONLY: GetGlobalNonUniqueSideID
#if USE_MPI
Expand All @@ -1040,21 +1040,17 @@ SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)
! OUTPUT VARIABLES
!-----------------------------------------------------------------------------------------------------------------------------------
! LOCAL VARIABLES
INTEGER :: iSide, BCSideID, ElemID, iLocSide, SideID, currentBC, PartInsSubSum, iSample, jSample
INTEGER :: SampleElemID
INTEGER, ALLOCATABLE :: PartInsSubSidesAdapt(:,:,:)
REAL :: VeloVec(1:3), vec_nIn(1:3), ElemPartDensity, VeloIC, VeloVecIC(1:3), projFak
REAL :: v_thermal, a, vSF, nVFR, RandVal1, area
INTEGER :: iSide, BCSideID, ElemID, SideID, currentBC, iSample, jSample, SampleElemID
REAL :: VeloVec(1:3), vec_nIn(1:3), nVFRTotal, VeloIC, VeloVecIC(1:3), projFak
REAL :: v_thermal, a, vSF, nVFR, area
!===================================================================================================================================

ASSOCIATE(SF => Species(iSpec)%Surfaceflux(iSF))

currentBC = SF%BC
SDEALLOCATE(PartInsSubSidesAdapt)
ALLOCATE(PartInsSubSidesAdapt(1:SurfFluxSideSize(1),1:SurfFluxSideSize(2),1:BCdata_auxSF(currentBC)%SideNumber))
PartInsSubSidesAdapt = 0

PartInsSubSum = 0
nVFRTotal = 0.
SF%ConstMassflowWeight = 0.

DO iSide=1,BCdata_auxSF(currentBC)%SideNumber
! Skip sides outside of the circular inflow region
Expand All @@ -1064,8 +1060,7 @@ SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)
BCSideID=BCdata_auxSF(currentBC)%SideList(iSide)
ElemID = SideToElem(S2E_ELEM_ID,BCSideID)
SampleElemID = AdaptBCMapElemToSample(ElemID)
iLocSide = SideToElem(S2E_LOC_SIDE_ID,BCSideID)
SideID=GetGlobalNonUniqueSideID(offsetElem+ElemID,iLocSide)
SideID=GetGlobalNonUniqueSideID(offsetElem+ElemID,SideToElem(S2E_LOC_SIDE_ID,BCSideID))
! Get the sampled velocity vector
VeloVec(1:3) = AdaptBCMacroVal(1:3,SampleElemID,iSpec)
! Determine the velocity magnitude
Expand All @@ -1089,7 +1084,6 @@ SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)
END IF
! Loop over the triangles
DO jSample=1,SurfFluxSideSize(2); DO iSample=1,SurfFluxSideSize(1)
ElemPartDensity = 0.
! Set the area of the side, different area for circular inflow
IF(SF%CircularInflow) THEN
area = SF%CircleAreaPerTriaSide(iSample,jSample,iSide)
Expand All @@ -1114,35 +1108,25 @@ SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)
vSF = v_thermal / (2.0*SQRT(PI)) * ( EXP(-(a*a)) + a*SQRT(PI)*(1+ERF(a)) ) !mean flux velocity through normal sub-face
nVFR = area * vSF !VFR projected to inwards normal of sub-side
CASE DEFAULT
CALL abort(__STAMP__,'ERROR in AdaptiveBoundary_ConstMassflow_Weight: Wrong velocity distribution!')
CALL abort(__STAMP__,'ERROR in CalcConstMassflowWeight: Wrong velocity distribution!')
END SELECT
! Skip the side if a negative/zero mean velocity has been determined
IF(vSF.LE.0.0) CYCLE
! Calculate the number density depending on the sampled velocity and target mass flow
ElemPartDensity = SF%AdaptiveMassflow / (vSF * SF%totalAreaSF * Species(iSpec)%MassIC)
! Skip the rest if the calculated number density is zero
IF(ElemPartDensity.EQ.0.) CYCLE
! Calculate the particle number per side
CALL RANDOM_NUMBER(RandVal1)
PartInsSubSidesAdapt(iSample,jSample,iSide) = INT(ElemPartDensity/Species(iSpec)%MacroParticleFactor*dt*RKdtFrac*nVFR+RandVal1)
IF(PartInsSubSidesAdapt(iSample,jSample,iSide).GE.0) THEN
PartInsSubSum = PartInsSubSum + PartInsSubSidesAdapt(iSample,jSample,iSide)
ELSE
IPWRITE(*,*) 'AdaptBCMacroVal', AdaptBCMacroVal(1:3,SampleElemID,iSpec), 'AdaptBCBackupVelocity', AdaptBCBackupVelocity(1:3,SampleElemID,iSpec)
IPWRITE(*,*) 'veloNormal', vSF, 'SF%totalAreaSF', SF%totalAreaSF
IPWRITE(*,*) 'ElemPartDensity', ElemPartDensity, 'MacroParticleFactor', Species(iSpec)%MacroParticleFactor
IPWRITE(*,*) 'dt', dt, 'nVFR', nVFR, 'PartInsSubSidesAdapt', PartInsSubSidesAdapt(iSample,jSample,iSide)
CALL abort(__STAMP__,'ERROR in AdaptiveBoundary_ConstMassflow_Weight: Negative number of particles to be inserted!')
END IF
! Skip the side if a negative/zero volume flow rate has been determined
IF(nVFR.LE.0.0) CYCLE
! Calculate the volume flow rate per side
SF%ConstMassflowWeight(iSample,jSample,iSide) = nVFR
! Calculate the total volume flow rate
nVFRTotal = nVFRTotal + nVFR
END DO; END DO
END DO

! Calculate the total volume flow rate over the whole BC across all processors
#if USE_MPI
CALL MPI_ALLREDUCE(MPI_IN_PLACE,PartInsSubSum,1,MPI_INTEGER,MPI_SUM,PartMPI%COMM,IERROR)
CALL MPI_ALLREDUCE(MPI_IN_PLACE,nVFRTotal,1,MPI_DOUBLE,MPI_SUM,PartMPI%COMM,IERROR)
#endif

IF(PartInsSubSum.GT.0) THEN
SF%ConstMassflowWeight(:,:,:) = REAL(PartInsSubSidesAdapt(:,:,:)) / REAL(PartInsSubSum)
! Determine the weight of each side compared to the total volume flow rate
IF(nVFRTotal.GT.0.) THEN
SF%ConstMassflowWeight(:,:,:) = SF%ConstMassflowWeight(:,:,:) / REAL(nVFRTotal)
ELSE
SF%ConstMassflowWeight(:,:,:) = 0.
END IF
Expand All @@ -1164,9 +1148,7 @@ SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight(iSpec,iSF)

END ASSOCIATE

SDEALLOCATE(PartInsSubSidesAdapt)

END SUBROUTINE AdaptiveBoundary_ConstMassflow_Weight
END SUBROUTINE CalcConstMassflowWeight


!===================================================================================================================================
Expand Down

0 comments on commit 6882692

Please sign in to comment.