using System;

public class Program

{

    public static void Main()

    {

        //ForceConduction(0.2f);

        ForceConductionRevised(0.2f);

    }

​

    public static void ForceConduction(float dt)

    {

        float conductivity = System.Math.Min(0.6f, 20f) * dt;

​

​

        float totalEnergyA = 300.15f * 4.179f * 1000;

        Console.WriteLine("BEFORE Water energy: " + totalEnergyA);

        float startingTempA = 300.15f;

​

        float totalEnergyB = 673.15f * 1.0f * 100;

        Console.WriteLine("BEFORE Uranium energy: " + totalEnergyB);

        float startingTempB = 673.15f;

​

        float energyRatio = (totalEnergyB / (totalEnergyA + totalEnergyB)) * conductivity;

​

        //float newA = ((startingTempB - startingTempA) * energyRatio) + startingTempA;

        //float newB = ((startingTempA - startingTempB) * (1 - energyRatio)) + startingTempB;

        float newB = (totalEnergyA + totalEnergyB)/(4.179f * 1000 + 1.0f * 100f);

        float newA = newB;

        Console.WriteLine("Water: " + newA);

        Console.WriteLine("Uranium: " + newB);

​

        totalEnergyA = newA * 4.179f * 1000;

        Console.WriteLine("AFTER Water energy: " + totalEnergyA);

        Console.WriteLine("AFTER Water TEMP: " + ((totalEnergyA / 4.179f) / 1000));

​

        totalEnergyB = newB * 1.0f * 100;

        Console.WriteLine("AFTER Uranium energy: " + totalEnergyB);

        Console.WriteLine("AFTER Uranium TEMP: " + ((totalEnergyB / 1.0f) / 100));

    }

    public static void ForceConductionRevised(float dt)

    {

        float conductivity = System.Math.Min(0.6f, 20f) * dt;

​

        float aTemp =  300.15f;

        float bTemp =  673.15f;

        float aSHC =  4.179f;

        float bSHC =  1.0f;

        float aMass = 1000f;

        float bMass = 100f;

        float totalEnergyA = 300.15f * 4.179f * 1000;

        Console.WriteLine("BEFORE Water energy: " + totalEnergyA);

        float startingTempA = 300.15f;

​

        float totalEnergyB = 673.15f * 1.0f * 100;

        Console.WriteLine("BEFORE Uranium energy: " + totalEnergyB);

        float startingTempB = 673.15f;

​

        float massRatio = bMass / aMass;

        if(massRatio < 1)

            massRatio = 1 / massRatio;

        float addToA = (bTemp - aTemp) * (massRatio) * conductivity * dt;

        Console.WriteLine("Add to A: " + addToA);

        float newA = ((totalEnergyA + addToA) / aSHC) / aMass;

        float newB = ((totalEnergyB - addToA) / bSHC) / bMass;

        //float newB = (totalEnergyA + totalEnergyB)/(4.179f * 1000 + 1.0f * 100f);

        //float newA = newB;

        Console.WriteLine("Water: " + newA);

        Console.WriteLine("Uranium: " + newB);

​

        totalEnergyA = newA * 4.179f * 1000;

        Console.WriteLine("AFTER Water energy: " + totalEnergyA);

        Console.WriteLine("AFTER Water TEMP: " + ((totalEnergyA / 4.179f) / 1000));

​

        totalEnergyB = newB * 1.0f * 100;

        Console.WriteLine("AFTER Uranium energy: " + totalEnergyB);

        Console.WriteLine("AFTER Uranium TEMP: " + ((totalEnergyB / 1.0f) / 100));

    }

}