13.4 LiguidCoolingofPhotovoltaicpanel光伏板液冷仿真Focus: Compared with previous examples, the focusof this example is about thermal management inrenewableenergy systems.6/33
13.4 Liquid Cooling of Photovoltaic panel 光伏板液冷仿真 Focus: Compared with previous examples, the focus of this example is about thermal management in renewable energy systems. 6/33
13.4 Single-phase fluid flow and heat transfer forthermal management of photovoltaic panelKnown: Cold water at T=20C flows into thechannel inlet with velocity u=0.0lm/s. The PV panelgenerates a volumetric heat source of g = 27333W/m3Assumption: (1) steady state, (2) laminar flow, (3)incompressible fluid, (4) constant fluid properties(5)negligibleradiativeheattransfer7/33
13.4 Single-phase fluid flow and heat transfer for thermal management of photovoltaic panel Known ： Cold water at Tf =20℃ flows into the channel inlet with velocity u=0.01m/s. The PV panel generates a volumetric heat source of q = 27333W/m3 . Assumption: (1) steady state, (2) laminar flow, (3) incompressible fluid, (4) constant fluid properties, (5) negligible radiative heat transfer 7/33
The domain of liguid cooling channelThe"Liquid cooling channel"Size:1000*900(mm*mm)Outerwallthickness:25 (mm)Innerwallthickness:50 (mm)Shell number:6Total height:100 (mm)Photovoltaicpanelthickness:30(mm)OutletInletThe liguid cooling channel is adjacent to the back plate of the photovoltaic panelThesolarirradiationis1000W/m28/33
The liquid cooling channel is adjacent to the back plate of the photovoltaic panel. The solar irradiation is 1000W / m ². The “Liquid cooling channel” Size:1000*900 (mm*mm) Outer wall thickness: 25 (mm) Inner wall thickness: 50 (mm) Shell number: 6 Total height: 100 (mm) Photovoltaic panel thickness: 30 (mm) Inlet Outlet The domain of liquid cooling channel 8/33
Find: Temperature, velocity and pressure distributions in the system.Governingequations:duavawContinuityequation:0axayazduduaudu10p(aua2uHMomentum equations:uVXW+ayazaxax2dy20z2praxPraauaua2v1210p1020Hf十u1Wdyozax2dy2axaz2PrdyPfawawa2waw10pa2wa2wH21Waayax0x2PrdzQy20z2pr9/33
Find: Temperature, velocity and pressure distributions in the system. Governing equations: Continuity equation: Momentum equations: 𝜕𝑢 𝜕𝑥 + 𝜕𝑣 𝜕y + 𝜕𝑤 𝜕𝑧 = 0 𝑢 𝜕𝑢 𝜕𝑥 + 𝑣 𝜕𝑢 𝜕𝑦 + 𝑤 𝜕𝑢 𝜕𝑧 = − 1 𝜌𝑓 𝜕𝑝 𝜕𝑥 + 𝜇𝑓 𝜌𝑓 𝜕 2𝑢 𝜕𝑥 2 + 𝜕 2𝑢 𝜕𝑦 2 + 𝜕 2𝑢 𝜕𝑧 2 𝑢 𝜕𝑣 𝜕𝑥 + 𝑣 𝜕𝑣 𝜕𝑦 + 𝑤 𝜕𝑣 𝜕𝑧 = − 1 𝜌𝑓 𝜕𝑝 𝜕𝑦 + 𝜇𝑓 𝜌𝑓 𝜕 2𝑣 𝜕𝑥 2 + 𝜕 2𝑣 𝜕𝑦 2 + 𝜕 2𝑣 𝜕𝑧 2 𝑢 𝜕𝑤 𝜕𝑥 + 𝑣 𝜕𝑤 𝜕𝑦 + 𝑤 𝜕𝑤 𝜕𝑧 = − 1 𝜌𝑓 𝜕𝑝 𝜕𝑧 + 𝜇𝑓 𝜌𝑓 𝜕 2𝑤 𝜕𝑥 2 + 𝜕 2𝑤 𝜕𝑦 2 + 𝜕 2𝑤 𝜕𝑧 2 9/33
Energy equation:(a2Ta(prCpuT)a(prCpruT)apfCprwT)a2Ta2Tax2axdyaz0z2a1where T,is the fluid temperature, Cpris fluid specific heat andk,isfluidthermal conductivityEnergy equation for the PV panel and liquid channel wall:a2Tpv(a2TpuaTp(a2T+a2Ts+a2Ts0 =ks0=kpvaz2ax2az2ax2y2ay2where T,, and T, is solid temperature, while kp and k, is solidthermal conductivity10/33
Energy equation: where Tf is the fluid temperature, cPf is fluid specific heat and kf is fluid thermal conductivity. Energy equation for the PV panel and liquid channel wall: where Tpv and Ts is solid temperature, while 𝒌𝒑𝒗 and ks is solid thermal conductivity 𝜕(𝜌𝑓𝐶𝑝𝑓𝑢𝑓𝑇𝒇 ) 𝜕𝑥 + 𝜕 𝜌𝑓𝐶𝑝𝑓𝑣𝑓𝑇𝒇 𝜕𝑦 + 𝜕 𝜌𝑓𝐶𝑝𝑓𝑤𝑓𝑇𝒇 𝜕𝑧 =𝑘𝑓 𝜕 2𝑇𝑓 𝜕𝑥 2 + 𝜕 2𝑇𝒇 𝜕𝑦 2 + 𝜕 2𝑇𝒇 𝜕𝑧 2 0 = 𝑘𝑝𝑣 𝜕 2𝑇𝒑𝒗 𝜕𝑥 2 + 𝜕 2𝑇𝒑𝒗 𝜕𝑦2 + 𝜕 2𝑇𝒑𝒗 𝜕𝑧 2 +s 0 = 𝑘𝑠 𝜕 2𝑇𝑠 𝜕𝑥 2 + 𝜕 2𝑇𝑠 𝜕𝑦2 + 𝜕 2𝑇𝑠 𝜕𝑧 2 10/33