github.com/archlabjp/eeslism-go@v0.0.0-20231109122333-4bb7bfcdf292/eeslism/blrmaceqcf.go

//This file is part of EESLISM.
//
//Foobar is free software : you can redistribute itand /or modify
//it under the terms of the GNU General Public License as published by
//the Free Software Foundation, either version 3 of the License, or
//(at your option) any later version.
//
//Foobar is distributed in the hope that it will be useful,
//but WITHOUT ANY WARRANTY; without even the implied warranty of
//MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.See the
//GNU General Public License for more details.
//
//You should have received a copy of the GNU General Public License
//along with Foobar.If not, see < https://www.gnu.org/licenses/>.

package eeslism

import "fmt"

/* -------------------------------------- */

var __Rmhtrcf_count int

func Rmhtrcf(exs *EXSFS, emrk []rune, rooms []*ROOM, sds []*RMSRF, wd *WDAT) {
	if rooms != nil {
		for i := range rooms {
			room := rooms[i]
			n := room.N
			brs := room.Brs
			sds := sds[brs : brs+n]

			if DEBUG {
				fmt.Printf("Room Name=%s\n", room.Name)
			}

			// 放射熱交換係数の計算
			if __Rmhtrcf_count == 0 || emrk[0] == '*' {
				radex(n, sds, room.F, room.Wradx)
			}
			if DEBUG {
				fmt.Printf("radex end\n")
			}

			// 放射熱伝達率の入れ替え
			Radcf0(room.Tsav, &room.alrbold, n, sds, room.Wradx, room.alr)
			if DEBUG {
				fmt.Printf("Radcf0 end\n")
			}

			Htrcf(room.alc, exs.Alosch, exs.Alotype, exs.Exs, room.Tr, n, room.alr, sds, &room.mrk, wd)
			if DEBUG {
				fmt.Printf("Htrcf end\n")
			}
		}

		__Rmhtrcf_count++
	}

	if sds != nil {
		for _, sd := range sds {
			if sd.mwtype == RMSRFMwType_C && sd.mwside == RMSRFMwSideType_i {
				// 内壁の場合は裏面室の熱伝達率を入れ替える
				nxsd := sd.nxsd
				sd.alo = nxsd.ali
				nxsd.alo = sd.ali
			}
		}
	}
}

/* ----------------------------------------------------------------- */

func Rmrdshfc(_Room []*ROOM, Sd []*RMSRF) {
	if len(_Room) == 0 {
		return
	}

	for i := range _Room {
		Room := _Room[i]
		brs := Room.Brs
		sd := Sd[brs:]

		Radshfc(Room.N, Room.FArea, Room.Area, sd, Room.tfsol, Room.eqcv, Room.Name, Room.fsolm)
	}
}

/* ----------------------------------------------------------------- */
func Rmhtrsmcf(_Sd []*RMSRF) {
	for n := range _Sd {
		Sd := _Sd[n]
		Sd.K = 1.0 / (Sd.Rwall + 1.0/Sd.ali + 1.0/Sd.alo)
	}
}

/* ----------------------------------------------------------------- */
// 透過日射、相当外気温度の計算
func Rmexct(Room []*ROOM, Sd []*RMSRF, Wd *WDAT, Exs []*EXSF, Snbk []*SNBK, Qrm []*QRM, nday, mt int) {
	if len(Room) == 0 {
		return
	}

	// 部位ごとの日射吸収比率のスケジュール対応（比率入力部位の日射入射比率初期化）
	for _, rm := range Room {

		// 室内部位の日射吸収比率の計算
		// 2017/12/25毎時計算へ変更
		// 家具の日射吸収割合
		if rm.fsolm != nil {
			rm.tfsol = *(rm.fsolm)
		} else {
			rm.tfsol = 0.0
		}

		for _, rsd := range rm.rsrf[:rm.N] {
			// 床の場合
			if rsd.ble == BLE_Floor || rsd.ble == BLE_InnerFloor {
				// どの部位も日射吸収比率が定義されていない場合
				if rm.Nfsolfix == 0 {
					// 床の日射吸収比率は固定
					rsd.ffix_flg = '*'
					rsd.fsol = CreateConstantValuePointer(*rm.flrsr * rsd.A / rm.FArea)
				}
			}

			// fsolが規定されている部位についてfsolを合計する
			if rsd.ffix_flg == '*' {
				rm.tfsol += *rsd.fsol // fsolの合計値計算
			}
		}
	}

	// 室内部位の日射吸収比率の計算（毎計算ステップへ変更）2017/12/25
	Rmrdshfc(Room, Sd)
	for i, rm := range Room {
		Q := Qrm[i]

		rm.Qgt = 0.0
		rm.Qsolm = 0.0
		rm.Qsab = 0.0
		rm.Qrnab = 0.0
		Q.Solo = 0.0
		Q.Solw = 0.0
		Q.Asl = 0.0

		RmSd := Sd[rm.Brs : rm.Brs+rm.N]

		for n, Sdn := range RmSd {

			if Sdn.ble == BLE_InnerWall || Sdn.ble == BLE_InnerFloor || Sdn.ble == BLE_Ceil {
				continue
			}

			e := Exs[Sdn.exs]

			Sdn.RSsol = 0.0
			Sdn.RSsold = 0.
			Fsdw := 0.0
			Qgtn := 0.0
			Qga := 0.0

			var Idre float64     // 直逹日射  [W/m2]
			var Idf float64      // 拡散日射  [W/m2]
			var RN float64       //  夜間輻射  [W/m2]
			if Sdn.Sname == "" { /*---higuchi 070918---start-*/
				Sdn.Fsdworg = 0.0

				sb := Sdn.sb
				if sb >= 0 && e.Cinc > 0.0 {
					S := Snbk[sb]

					// 日よけの影面積率 [-]
					Fsdw = FNFsdw(S.Type, S.Ksi, e.Tazm, e.Tprof, S.D, S.W, S.H, S.W1, S.H1, S.W2, S.H2)
					Sdn.Fsdworg = Fsdw

					if DEBUG {
						fmt.Printf(" xxx Rmexct xxx i=%d n=%d  sb=%d  type=%d  tazm=%f  tprof=%f Fsdw=%f\n",
							i, n, sb, S.Type, e.Tazm, e.Tprof, Fsdw)
					}
				} else {
					Fsdw = 0.0
				}

				Idre = e.Idre
				Idf = e.Idf
				RN = e.Rn

				if DEBUG {
					fmt.Printf("1:Fsdw=%f,Idre=%f,Idf=%f,rn=%f Sb=%d Cinc=%f\n", Fsdw, Idre, Idf, RN, sb, e.Cinc)
				}

			} else { /*---higuchi 070918 end--*/ /*--higuchi 070918 start--*/
				Fsdw = Sdn.Fsdw
				//                  Idre = Sdn.Idre ;  090131 higuchi Sdn.Idre が影をすでに考慮していたため、下に変更
				Idre = e.Idre /*--090131 higuchi  --*/
				Idf = Sdn.Idf
				RN = Sdn.rn
			} /*---higuchi 070918 end --*/

			switch Sdn.ble {
			case BLE_Window:
				// 通常窓の場合
				/*--higuchi add--*/
				Qgtn, Qga = Glasstga(Sdn.A, Sdn.tgtn, Sdn.Bn,
					e.Cinc, Fsdw, Idre, Idf, Sdn.window.Cidtype, e.Prof, e.Gamma)
				Rab := Sdn.Eo * RN / Sdn.alo // 夜間放射熱取得 [W]

				Sab := Qga / Sdn.A // [W/m2]
				Sdn.TeEsol = Sab
				Sdn.TeErn = -Rab
				Sdn.TeEsol = Sab / Sdn.K
				Sdn.Te = Sab/Sdn.K - Rab + Wd.T // 外表面の相当外気温
				Sdn.Qgt = Qgtn                  // 開口部の透過日射熱取得 [W]
				Sdn.Qga = Qga                   // 開口部の吸収日射熱取得 [W]
				Sdn.Qrn = -Rab                  // 開口部の夜間放射熱取得 [W/m2]

				// 部屋rm の日射
				rm.Qgt += Qgtn                  // 部屋rmの透過日射熱取得 [W]
				rm.Qsab += Sab * Sdn.A          // 部屋rmの吸収日射熱取得 [W]
				rm.Qrnab += Rab * Sdn.A * Sdn.K // 部屋rmの夜間放射による熱損失 [W]

				Q.Solw += Sdn.A * (Idre + Idf) /*--higuchi add  --*/
				break

			case BLE_ExternalWall, BLE_Floor, BLE_Roof: // このあたりを参考に修正（相当外気温度の計算）
				if Sdn.typ != RMSRFType_E && Sdn.typ != RMSRFType_e {
					/*---higuchi add---*/
					Sab := Sdn.as * (Idre*(1.0-Fsdw) + Idf) / Sdn.alo
					Rab := Sdn.Eo * RN / Sdn.alo // 長波長
					/*------------------*/

					Sdn.TeEsol = Sab
					Sdn.TeErn = -Rab

					// 建材一体型空気集熱器のための相当外気温度修正
					if Sdn.rpnl != nil && Sdn.rpnl.Type == 'C' {
						//wall := Sdn.mw.wall
						Sdn.Te = Sdn.Tcole
						Sdn.Iwall = Idre*(1.0-Fsdw) + Idf
					} else {
						Sdn.Te = Sab - Rab + Wd.T
					}

					rm.Qsab += Sab * Sdn.A * Sdn.K
					Sdn.Qga = Sab * Sdn.A * Sdn.K
					rm.Qrnab += Rab * Sdn.A * Sdn.K
					Q.Solo += Sdn.A * (Idre + Idf)
					Q.Asl += Sdn.as * Sdn.A * (Idre + Idf)
				} else {
					Sdn.Te = e.Tearth
					Sdn.TeEsol = 0.0
					Sdn.TeErn = 0.0
				}
				break

			case BLE_InnerWall, BLE_InnerFloor, BLE_Ceil, BLE_d:
				if Sdn.nxrm < 0 {
					// 隣室が無い場合
					Tr := Sdn.room.Trold
					Eo := Sdn.room.cmp.Elouts[0]
					if Eo.Control == LOAD_SW {
						Tr = Sdn.room.rmld.Tset
					}

					// 相当外気温を設定
					Sdn.Te = Sdn.c*Tr + (1.0-Sdn.c)*Wd.T
				} else {
					// 隣室がある場合
					Tr := Sdn.nextroom.Trold
					Eo := Sdn.nextroom.cmp.Elouts[0]
					if Eo.Control == LOAD_SW {
						Tr = Sdn.nextroom.rmld.Tset
					}

					// 相当外気温を設定
					Sdn.Te = Sdn.c*Tr + (1.0-Sdn.c)*Wd.T
				}
				Sdn.TeEsol = 0.0
				Sdn.TeErn = 0.0
				break
			}
		} // 表面ループ

		// 室内部位への入射日射の計算（吸収日射ではない）
		for _, Sdn := range RmSd {
			// 室内部位への入射日射量の計算
			Sdn.RSsold = rm.Qgt * Sdn.srg
		}
	} // 室ループ終了

	// 透過日射の室内部位の最終計算（隣接室への日射分配、透過日射のうちガラスから屋外に放熱される分も考慮）
	for _, rm := range Room {
		// 透過間仕切りなど、隣接空間への透過日射分配の計算
		for _, Sdn := range Sd[rm.Brs : rm.Brs+rm.N] {

			// 隣室への日射の分配率が指定されており、分配する入射がある
			if Sdn.tnxt > 0. && Sdn.RSsold > 0. {
				// 分配する入射量 RSsol
				Rmnxt := Room[Sdn.nxrm]
				RSsol := Sdn.RSsold * Sdn.tnxt

				// 入射日射×透過率が当該室の透過日射熱取得より減ずる
				rm.Qgt -= RSsol
				if Sdn.nextroom != nil {
					// 外皮でない場合は隣室の透過日射熱取得に透過分を加算
					Rmnxt.Qgt += RSsol
				}
			}

			// 透過日射が入射したときに屋外に放熱されるときには、表面吸収日射はゼロとする
			if Sdn.RStrans {
				rm.Qgt -= Sdn.RSsold
				Sdn.RSsol = 0.0
			}
		}
	}

	for i, rm := range Room {
		Q := Qrm[i]
		RmSd := Sd[rm.Brs : rm.Brs+rm.N]

		// 室内部位の短波長吸収量の計算
		for n, Sdn := range RmSd {
			Sdn.RSsol = rm.Qgt * Sdn.srg2 / Sdn.A

			Sdn.RS = (Sdn.RSsol*Sdn.A + rm.Hr*Sdn.srh + rm.Lr*Sdn.srl + rm.Ar*Sdn.sra + rm.Qeqp*Sdn.eqrd) / Sdn.A
			Sdn.RSin = (rm.Hr*Sdn.srh + rm.Lr*Sdn.srl + rm.Ar*Sdn.sra + rm.Qeqp*Sdn.eqrd) / Sdn.A
			Sdn.RSli = rm.Lr * Sdn.srl / Sdn.A

			if DEBUG {
				fmt.Printf("----- Rmexct  i=%d n=%d rm.Qgt=%f Fsdw=%f Fsdworg=%f Qgt=%f Qga=%f srg2=%f RSsol=%f RS=%f RSin=%f RSli=%f\n",
					i, n, rm.Qgt, Sdn.Fsdw, Sdn.Fsdworg, Sdn.Qgt, Sdn.Qga, Sdn.srg2, Sdn.RSsol, Sdn.RS, Sdn.RSin, Sdn.RSli)
			}
		}

		// 室の透過日射熱取得を再度積算（透明間仕切りによる隣接空間からの透過日射を考慮するため）
		if rm.rsrnx {
			for _, Sdn := range RmSd {
				if Sdn.ble == BLE_Ceil || Sdn.ble == BLE_InnerFloor {
					if Sdn.nxn >= 0 {
						Sdn.Te += Sd[Sdn.nxn].RS / Sdn.alo
					}
				}
			}
		}

		Q.Tsol = rm.Qgt
		Q.Asol = rm.Qsab
		Q.Arn = rm.Qrnab

		// 家具の日射吸収量の計算
		rm.Qsolm = 0.
		if rm.fsolm != nil {
			rm.Qsolm = rm.Qgt * rm.Srgm2
		}

	} // 室ループ

	for n, Sdn := range Sd {
		// 共用壁の場合の外表面の相当外気温 Te [C] の計算
		if Sdn.mwtype == RMSRFMwType_C {
			Sdnx := Sdn.nxsd
			Sdn.Te = (Sdnx.alir*Sdnx.Tmrt + Sdnx.RS) / Sdnx.ali

			if DEBUG {
				fmt.Printf("----- Rmexct  n=%d Te=%f nxalir=%f nxTmrt=%f nxSdnx->RS=%f nxali=%f\n",
					n, Sdn.Te, Sdnx.alir, Sdnx.Tmrt, Sdnx.RS, Sdnx.ali)
			}
		}
	}
}

/* ----------------------------------------------------------------- */

// 室の係数、定数項の計算
func Roomcf(mw []*MWALL, rooms []*ROOM, rdpnl []*RDPNL, wd *WDAT, exsf *EXSFS) {
	for _, rdpnl := range rdpnl {
		panelwp(rdpnl)
	}

	// 壁体係数行列の作成（壁体数RMSRF分だけループ）
	RMwlc(mw, exsf, wd)

	for i := range rooms {
		room := rooms[i]

		RMcf(room)
		RMrc(room) // 室の定数項の計算

		room.RMx = room.GRM / DTM
		room.RMXC = room.RMx*room.xrold + (room.HL+room.AL)/Ro

		room.RMt += Ca * room.Gvent
		room.RMC += Ca * room.Gvent * wd.T
		room.RMx += room.Gvent
		room.RMXC += room.Gvent * wd.X
	}

	for _, rdpnl := range rdpnl {
		Panelcf(rdpnl)
		rdpnl.EPC = Panelce(rdpnl)
	}
}

/* ----------------------------------------------------------------- */
// 前時刻の室温の入れ替え、OT、MRTの計算
func Rmsurft(rooms []*ROOM, sd []*RMSRF) {
	if rooms == nil {
		return
	}

	// 重み係数が未定義もしくは不適切な数値の場合の対処
	r := 0.5
	if rooms[0].OTsetCwgt != nil && *(rooms[0].OTsetCwgt) >= 0.0 && *(rooms[0].OTsetCwgt) <= 1.0 {
		r = *(rooms[0].OTsetCwgt)
	}

	if DEBUG {
		fmt.Printf("<Rmsurft> Start\n")
	}

	for i := range rooms {
		room := rooms[i]
		n := room.N
		brs := room.Brs
		sdr := sd[brs:]

		if DEBUG {
			fmt.Printf("Room[%d]=%s\tN=%d\tbrs=%d\n", i, room.Name, room.N, room.Brs)
		}

		// 前時刻の温度の入れ替え
		room.mrk = 'C'
		room.Trold = room.Tr
		room.xrold = room.xr

		if room.FunHcap > 0 {
			// 家具の温度の計算
			room.TM = room.FMT*room.Tr + room.FMC
			// 家具の吸放熱量の計算
			if room.CM != nil {
				room.QM = *room.CM * (room.TM - room.Tr)
			}

			room.oldTM = room.TM
		}

		if DEBUG {
			fmt.Printf("<Rmsurft>  RMsrt start\n")
		}

		// 室内表面温度の計算
		RMsrt(room)

		if DEBUG {
			fmt.Printf("<Rmsurft>  RMsrt end\n")
		}

		room.Tsav = RTsav(n, sdr)                // 平均表面温度 Ts,av
		room.Tot = r*room.Tr + (1.0-r)*room.Tsav // 作用温度 Tot
	}
}

/* ----------------------------------------------------------------- */
// PCM収束計算過程における部位表面温度の計算
func Rmsurftd(_Room []*ROOM, Sd []*RMSRF) {
	var r float64

	if _Room == nil {
		return
	}

	Room := _Room[0]

	if Room.OTsetCwgt == nil || *(Room.OTsetCwgt) < 0.0 || *(Room.OTsetCwgt) > 1.0 {
		r = 0.5
	} else {
		r = *(Room.OTsetCwgt)
	}

	if DEBUG {
		fmt.Printf("<Rmsurft> Start\n")
	}

	for i := range _Room {
		Room := _Room[i]

		if DEBUG {
			fmt.Printf("Room[%d]=%s\tN=%d\tbrs=%d\n", i, Room.Name, Room.N, Room.Brs)
		}

		N := Room.N
		brs := Room.Brs
		sd := Sd[brs:]

		if DEBUG {
			fmt.Printf("<Rmsurft>  RMsrt start\n")
		}

		// 室内表面温度の計算
		RMsrt(Room)

		if DEBUG {
			fmt.Printf("<Rmsurft>  RMsrt end\n")
		}

		Room.Tsav = RTsav(N, sd)
		Room.Tot = r*Room.Tr + (1.0-r)*Room.Tsav
	}
}

/*--------------------------------------------------------------------------------------------*/

// 室の熱取得要素の計算
func Qrmsim(Room []*ROOM, Wd *WDAT, Qrm []*QRM) {
	for i := range Room {
		Q := Qrm[i]
		rm := Room[i]

		// 人体・照明・機器の顕熱 [W]
		Q.Hums = rm.Hc + rm.Hr
		Q.Light = rm.Lc + rm.Lr
		Q.Apls = rm.Ac + rm.Ar
		Q.Hgins = Q.Hums + Q.Light + Q.Apls

		// 人体・機器の潜熱 [W]
		Q.Huml = rm.HL
		Q.Apll = rm.AL

		// 熱負荷 [W]
		Q.Qinfs = Ca * rm.Gvent * (Wd.T - rm.Tr)
		Q.Qinfl = Ro * rm.Gvent * (Wd.X - rm.xr)
		Q.Qeqp = rm.Qeqp
		Q.Qsto = rm.MRM * (rm.Trold - rm.Tr) / DTM
		Q.Qstol = rm.GRM * Ro * (rm.xrold - rm.xr) / DTM

		// 電力の消費量 [W]
		if rm.AEsch != nil {
			Q.AE = rm.AE * *rm.AEsch
		} else {
			Q.AE = 0.0
		}

		// ガスの消費量 [W]
		if rm.AGsch != nil {
			Q.AG = rm.AG * *rm.AGsch
		} else {
			Q.AG = 0.0
		}
	}
}