Cosmos::Physics::ThermalPropagator Class Reference

#include <physicsclass.h>

Inheritance diagram for Cosmos::Physics::ThermalPropagator:

Collaboration diagram for Cosmos::Physics::ThermalPropagator:

Public Member Functions
	ThermalPropagator (locstruc *locp, physicsstruc *physp, double idt, float itemperature)
int32_t	Init (float temp)
int32_t	Propagate ()
Public Member Functions inherited from Cosmos::Physics::Propagator
	Propagator (locstruc *locp, physicsstruc *physp, double idt)
int32_t	Increment (double nextutc)

Public Attributes
float	temperature
Public Attributes inherited from Cosmos::Physics::Propagator
double	dt
double	dtj
locstruc	oldloc
locstruc *	newloc
physicsstruc	oldphys
physicsstruc *	newphys
Type	type

Additional Inherited Members
Public Types inherited from Cosmos::Physics::Propagator
enum	Type {   None = 0, PositionInertial = 10, PositionIterative = 11, PositionGaussJackson = 12,   AttitudeInertial = 20, AttitudeIterative = 21, AttitudeLVLH = 22, Thermal = 30,   Electrical = 40 }

Constructor & Destructor Documentation

Cosmos::Physics::ThermalPropagator::ThermalPropagator ( locstruc *  locp, physicsstruc *  physp, double  idt, float  itemperature  )

inline

                : Propagator{ locp, physp, idt }, temperature{itemperature}
            {
                type = AttitudeLVLH;
                newphys->temp = temperature;
                newphys->heat = newphys->temp * (newphys->mass * newphys->hcap);
                newphys->radiation = SIGMA * pow(temperature, 4.);
                for (trianglestruc& triangle : newphys->triangles)
                {
                    triangle.temp = temperature;
                    triangle.heat = triangle.temp * (triangle.mass * triangle.hcap);
                }
            }

Member Function Documentation

int32_t Cosmos::Physics::ThermalPropagator::Init

(

float

temp

)

int32_t Cosmos::Physics::ThermalPropagator::Propagate

(

)

        {
            Vector units = Quaternion(newloc->att.icrf.s).irotate(Vector(newloc->pos.icrf.s).normalize());
            Vector unite = Vector(newloc->pos.eci.s).normalize(-1.);
            unite = Vector(newloc->pos.eci.s).normalize(-1.);
            unite = Quaternion(newloc->att.icrf.s).irotate(Vector(newloc->pos.eci.s).normalize(-1.));
            double energyd;
            double sdot;
            double edot;
            newphys->temp = newphys->heat / (newphys->mass * newphys->hcap);
            newphys->heat = 0.;
            double ienergyd = newphys->radiation / newphys->area;
            for (trianglestruc& triangle : newphys->triangles)
            {
                triangle.temp = triangle.heat / (triangle.mass * triangle.hcap);

                // First do all inputs

                if (triangle.external)
                {
                    // External Normal face
                    // Solar effects
                    sdot = units.dot(triangle.normal);

                    // Earth effects
                    edot = acos(unite.dot(triangle.normal) / triangle.normal.norm()) - RADOF(5.);
                    if (edot < 0.)
                    {
                        edot = 1.;
                    }
                    else
                    {
                        edot = cos(edot);
                    }

                    energyd = 0.;
                    if (newloc->pos.sunradiance && sdot > 0)
                    {
                        energyd +=  sdot * triangle.irradiation * dt;
                    }

                    if (edot > 0)
                    {
                        energyd += edot * dt * SIGMA * pow(290.,4);
                    }

                    // Area not covered with cells
                    triangle.heat += (1. - triangle.pcell) * triangle.area * triangle.abs * energyd;
                    if (triangle.pcell > 0.)
                    {
                        // Area covered with cells
                        if (triangle.ecellbase > 0.)
                        {
                            double efficiency = triangle.ecellbase + triangle.ecellslope * triangle.temp;
                            triangle.heat += (triangle.pcell) * triangle.area * (triangle.abs - efficiency) * energyd;
                        }
                        else {
                            triangle.heat += (triangle.pcell) * triangle.area * triangle.abs * energyd;
                        }
                    }

                    if (triangle.external == 2)
                    {
                        // External Anti-Normal face
                        // Solar effects
                        sdot = units.dot(-triangle.normal);

                        // Earth effects
                        edot = acos(unite.dot(-triangle.normal) / -triangle.normal.norm()) - RADOF(5.);
                        if (edot < 0.)
                        {
                            edot = 1.;
                        }
                        else
                        {
                            edot = cos(edot);
                        }

                        energyd = 0.;
                        if (newloc->pos.sunradiance && sdot > 0)
                        {
                            energyd +=  sdot * triangle.irradiation * dt;
                        }

                        if (edot > 0)
                        {
                            energyd += edot * dt * SIGMA * pow(290.,4);
                        }

                        triangle.heat += triangle.area * triangle.abs * energyd;
                    }
                    else {
                        // Internal Anti-Normal face
                        triangle.heat += triangle.iabs * triangle.area * ienergyd;
                        newphys->radiation -= triangle.iabs * triangle.area * ienergyd;
                    }
                }
                else {
                    // Internal Normal face
                    triangle.heat += triangle.iabs * triangle.area * ienergyd;
                    newphys->radiation -= triangle.iabs * triangle.area * ienergyd;

                    // Internal Anti-Normal face
                    triangle.heat += triangle.iabs * triangle.area * ienergyd;
                    newphys->radiation -= triangle.iabs * triangle.area * ienergyd;
                }

                // Self emission
                //                energyd = dt * SIGMA * pow(triangle.temp ,4);
                //                triangle.heat -= triangle.iemi * triangle.area * energyd;
                //                newphys->radiation += triangle.iemi * triangle.area * energyd;

                // Then do outputs
                energyd = dt * SIGMA * pow(triangle.temp ,4);
                if (triangle.external)
                {
                    // External Normal face
                    triangle.heat -= triangle.emi * triangle.area * energyd;
                    if (triangle.external == 2)
                    {
                        // External Anti-Normal face
                        triangle.heat -= triangle.emi * triangle.area * energyd;
                    }
                    else {
                        // Internal Anti-Normal face
                        triangle.heat -= triangle.iemi * triangle.area * energyd;
                        newphys->radiation += triangle.iemi * triangle.area * energyd;
                    }
                }
                else {
                    {
                        // Internal Normal face
                        triangle.heat -= triangle.iemi * triangle.area * energyd;
                        newphys->radiation += triangle.iemi * triangle.area * energyd;

                        // Internal Anti-Normal face
                        triangle.heat -= triangle.iemi * triangle.area * energyd;
                        newphys->radiation += triangle.iemi * triangle.area * energyd;

                    }
                }

                newphys->heat += triangle.heat;
                triangle.temp = triangle.heat / (triangle.mass * triangle.hcap);
            }
            newphys->temp = newphys->heat / (newphys->mass * newphys->hcap);
            temperature = newphys->temp;

            return 0;
        }

Member Data Documentation

float Cosmos::Physics::ThermalPropagator::temperature

---

The documentation for this class was generated from the following files:

• physicsclass.h
• physicsclass.cpp