Colder Body to Hotter

In[]:=
SeedRandom[34535];ArrayPlot[#,ColorRules->{0->White,1->Lighter[Orange],2->Darker[Orange]},ImageSize->{Automatic,500}]&@ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},300],RandomChoice[{.2,.8}->{0,2},600],RandomChoice[{.8,.2}->{0,2},300]],1000]
Out[]=
In[]:=
Mean[{{a,b},{c,d}}]
Out[]=

a+c
2
,
b+d
2

In[]:=
SeedRandom[34535];ListLinePlot[MovingAverage[Mean[#],50]&/@Partition[ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},100],RandomChoice[{.2,.8}->{0,2},200],RandomChoice[{.8,.2}->{0,2},100]],1000],50],PlotTheme->"Minimal"]
Out[]=
In[]:=
SeedRandom[34535];ListDensityPlot[MovingAverage[Mean[#],50]&/@Partition[ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},100],RandomChoice[{.2,.8}->{0,2},200],RandomChoice[{.8,.2}->{0,2},100]],1000],50],MeshFunctions->(#3&),Mesh->True]
Out[]=
In[]:=
SeedRandom[34535];ListPlot3D[MovingAverage[Mean[#],150]&/@Partition[ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},300],RandomChoice[{.1,.9}->{0,2},600],RandomChoice[{.8,.2}->{0,2},300]],4000],150],ColorFunction->"Rainbow",MeshFunctions->(#3&),Mesh->10,Ticks->None]
Out[]=
SeedRandom[34535];ListPlot3D[MovingAverage[Mean[#],150]&/@Partition[ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},300],RandomChoice[{.1,.9}->{0,2},600],RandomChoice[{.8,.2}->{0,2},300]],4000],150],ColorFunction->"Rainbow",MeshFunctions->(#3&),Mesh->10,Ticks->None]
In[]:=
SeedRandom[34535];ListPlot3D[MovingAverage[Mean[#],250]&/@Partition[ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},300],RandomChoice[{.1,.9}->{0,2},600],RandomChoice[{.8,.2}->{0,2},300]],4000],150],ColorFunction->"Rainbow",MeshFunctions->(#3&),Mesh->10,Ticks->None,Filling->Bottom]
Out[]=
In[]:=
SeedRandom[34535];ListPlot3D[MovingAverage[Mean[#],250]&/@Partition[ResourceFunction["BlockCellularAutomaton"][{{2,2}{1,1},{1,1}{2,2},{1,2}{1,2},{2,1}{2,1},{2,0}{0,2},{1,0}{1,0},{0,2}{2,0},{0,1}{0,1},{0,0}{0,0}},Join[RandomChoice[{.8,.2}->{0,2},300],RandomChoice[{.1,.9}->{0,2},600],RandomChoice[{.8,.2}->{0,2},300]],4000],150],ColorFunction->"Rainbow",MeshFunctions->(#2&),Mesh->10,Ticks->None,Filling->Bottom]
Out[]=

Heat Engine

isothermal, adiabatic, isothermal, adiabatic