Het Computerprogramma

Om in staat te zijn te verifiëren wat er eigenlijk gebeurd werd een simulatieprogramma ontwikkeld. Dit programma laat ons volgende simulaties toe:

• Gebruik van meerdere brandstofblokken, met twee opties:

  • geen verbranding aan de top en onderzijde

  • verbranding vanuit alle niet-geïsoleerde zijden (bates)

• Cylindrisch- of stervormidge vormen

• Toevoegiging van barsten, gedeeltelijk of totaal, op elk ogenblik tijdens het branden (in het programma gaan deze barsten geheel door de korrel)

• Erosie

• Ontstekingsvertraging

• Initiële onsteking

Met erosie bedoelen wij de situatie waar de regressie van het brandoppervlak sneller is dan normaal wegens de hoge uitstroomsnelheid van de reactieproducten op brandoppervlak. In ons programma hebben wij de volgende factor aan de brandsnelheid toegevoegd:

(1+ erosiecoëfficient x uistroomsnelheid)

Voor de berekening van de brandsnelheid (zonder erosie) werd er gebruik gemaakt van de zogenaamde wet van St. Roberts.

r=b+kPⁿ

Om de simulaties niet te ingewikkeld te maken werd "b"altijd als nul beschouwd.
Het gevolg van erosie is dat de regressie van het brandoppervlak niet dezelfde zal zijn over het gehele brandstofblok. Sommige delen zullen veel sneller dan anderen branden, dusdanig dat de lege ruimte spoedig van een cilindrische vorm zal afwijken. Om aan dit effect het hoofd te bieden werd de brandstofblok lengte verdeeld in 50 segmenten met constante snelheid.

Het programma houdt ook rekening met de drukveranderingen and het vrije volume. Er wordt echter geen rekening gehouden met de drukdaling over het brandstofblok van boven tot uitlaatpijp.

Voor de thermodynamische waarden van de reactieproducten (moleculegewicht, temperatuur, verhouding van de specifieke hitte...) werden de resultaten van GUIPEP gebruikt. In overeenstemming met de berekeningen van GUIPEP zal een groot deel van de reactieproducten waarschijnlijk uit condenseerbaar materiaal bestaan (hoofdzakelijk K₂C0₃) welke ook in acht genomen worden.