Lattice Boltzmann modeling for conduction heat transfer based on the Qadyan method

Technische Samenvatting: Lattice Boltzmann-modellering voor geleidingswarmteoverdracht gebaseerd op de Qadyan-methode

Probleemstelling
Het artikel behandelt de simulatie van geleidingswarmteoverdracht met behulp van de Lattice Boltzmann Method (LBM). Hoewel LBM effectief is voor fluïdumdynamica, kent de toepassing ervan op warmtegeleiding beperkingen wat betreft roosteruniformiteit en numerieke stabiliteit. Standaard LBM-schema's, met name die gebaseerd op de Bhatnagar–Gross–Krook (BGK)-benadering, zijn vaak beperkt tot uniforme roosters met gelijke ruimtelijke stappen. Bovendien kunnen standaard ruimtelijke discretiemethoden, zoals tweede-orde eindverschil-schema's, numerieke oscillaties in de oplossing veroorzaken. De studie beoogt alternatieve ruimtelijke discretiestrategieën te evalueren om de nauwkeurigheid en stabiliteit te verbeteren zonder extra complexiteit toe te voegen aan het standaard LBM-raamwerk of terug te grijpen op niet-uniforme roosters.

Methodologie
De studie implementeert en vergelijkt drie verschillende benaderingen voor het oplossen van de lattice Boltzmann-vergelijking (LBE) voor warmtegeleiding:

1. Standaard LBM: Maakt gebruik van conventionele roostermodellen (D2Q9, D2Q8 en D2Q5) met standaard botsings-stroomprocessen. De studie merkt op dat voor zuiver diffusieproblemen modellen met minder snelheidssets (zoals D2Q8 en D2Q5) voldoende zijn, aangezien warmtediffusie geen directionele voorkeur heeft.
2. Finite Difference Lattice Boltzmann Method (FDLBM): Deze benadering discretiseert de convectieve term van de LBE met behodele eindverschil-schema's. Specifiek wordt een eerste-orde voorwaartse differentie gebruikt voor tijdsafgeleiden en een eerste-orde upwind-schema voor ruimtelijke afgeleiden om oscillaties te voorkomen, hoewel dit de ruimtelijke nauwkeurigheid nabij de grenzen beperkt tot de eerste orde.
3. De Qadyan-methode: Een nieuwe benadering voorgesteld door Ghadyani et al. die semi-discrete schema's combineert met de Differential Quadrature Method (DQM). In deze methode worden de ruimtelijke afgeleiden in de LBE gediscretiseerd met behulp van DQM-wegingscoëfficiënten. Om het upwind-effect te repliceren dat inherent is aan advectie-diffusieproblemen, modificeert de methode de wegingscoëfficiënten door meer gewicht toe te kennen aan upstream-knopen dan aan downstream-knopen. Dit maakt het mogelijk de LBE op te lossen als een systeem van gewone differentiaalvergelijkingen zonder gebonden te zijn aan een specif kind type discrete snelheidsmodellen.

De studie valideert deze methoden aan de hand van een eendimensionale transiënte warmtegeleiding in een plaat met zowel Dirichlet- (constante temperatuur) als Neumann- (geïsoleerd/nul flux) randvoorwaarden. Analytische oplossingen afgeleid van de warmtevergelijking dienen als ijkpunten. Foutnormen ($L_1$, $L_2$ en $L_\infty$) worden berekend om de nauwkeurigheid te beoordelen, en een 2D-plaatdiffusieprobleem wordt opgelost om de toepasbaarheid van de methode in meerdimensionale gevallen aan te tonen.

Belangrijkste Bijdragen en Resultaten

• Nauwkeurigheid versus Stabiliteit: De resultaten geven aan dat hoewel de traditionele LBM een hoge nauwkeurigheid biedt, de voorgestelde Qadyan-methode een robuust alternatief vormt dat de nauwkeurigheid verbetert ten opzichte van de FDLBM-benadering, met name voor Neumann-randvoorwaarden. Voor Dirichlet-condities vertoonde de standaard D2Q5 LBM de laagste foutnormen, maar de Qadyan-methode presteerde significant beter dan de FDLBM-formuleringen.
• Roosteronafhankelijkheid: De Qadyan-methode slaagt erin het nauwkeurigheids- en stabiliteitsbereik van oplossingen te verbreden, waardoor de nauwkeurigheid wordt verhoogd zonder extra complexiteit aan de standaard lattice Boltzmann-methode toe te voegen of terug te vallen op niet-uniforme roosters.
• Computationele Kosten: Een cruciale bevinding is de afruil tussen nauwkeurigheid en efficiëntie. De Qadyan-methode brengt hogere computationele kosten en een langere simulatietijd met zich mee vergeleken met zowel de standaard LBM als de FDLBM, vanwege het grotere aantal vereiste berekeningsoperaties voor de DQM-wegingscoëfficiënten.
• Parametersensitiviteit: De studie stelt vast dat voor standaard LBM en FDLBM de referentiewaarden voor de rooster-snelheid afhankelijk zijn van het type randvoorwaarde (temperatuur versus flux). In tegen tegenstelling hiertoe gebruikt de Qadyan-methode een consistente rooster-snelheidswaarde ($c=0.015$) ongeacht het type randvoorwaarde, wat de selectie van parameters vereenvoudigt.
• 2D-Validatie: In een tweedimensionale vierkante plaat-simulatie produceerde de Qadyan-methode isotherm-contouren en middellijn-temperatuurverdelingen die een goede overeenstemming vertoonden met de eindverschilmethode-resultaten, wat de capaciteit van de methode voor meerdimensionale warmteoverdrachtsproblemen aantoont.

Betekenis en Claims
Het artikel beweert dat de Qadyan-methode een geschikte benadering is voor een breed scala aan geleidingswarmteoverdrachtsproblemen. De primaire betekenis ligt in het vermogen om de numerieke nauwkeurigheid te verhogen zonder de complexiteit van niet-uniforme roosters te introduceren of het LBM-raamwerk te verlaten. De auteurs stellen dat de methode robuust en eenvoudig te implementeren is. Het artikel neemt echter een bescheiden standpunt in met betrekking tot efficiëntie en erkent dat de methode computationeel duurder is dan traditionele benaderingen. De studie concludeert dat hoewel de Qadyan-methode een langere rekentijd vereist, de verbeterde nauwkeurigheid en stabiliteit voor specifieke randvoorwaarden de toepassing ervan rechtvaardigen, waarmee het een competitief alternatief biedt voor standaard LBM-formuleringen voor geleidingswarmteoverdrachtsimulaties.