Multi-Granularity Modeling Method for Effectiveness Evaluation of Remote Sensing Satellites

: The effectiveness indicator system of remote sensing satellites includes various satellites capabilities. Effectiveness evaluation is the process of calculating these indicators in the digital world, involving many different physical parameters of multiple subsystems. Model-based simulation statistics method is the mainstream approach of effectiveness evaluation, and digital twin is currently the most advanced modeling method for simulation. The satellite digital twin model has the characteristics of multi-dynamic, multi-spatial scale and multi-physics ﬁeld coupling, which gives rise to challenges related to the stiff problem of ordinary differential equations and multi-scale problem of partial differential equations to the calculation process of indicators. It is difﬁcult to solve these problems by breakthroughs in numerical solution methods. This paper uses the sparsity of the satellite system to group each indicator of the effectiveness evaluation indicator system according to the change period. The satellite system model is decomposed into multiple modules according to the composition and structure, and a series of models with different simulation ﬁdelity are established for each module. The optimization schemes for selecting model granularity when calculating indicators by group is given. Simulation results show that this approach considers the coupling between systems, grasps the main contradiction of indicator calculation and overcomes the loss of indicator accuracy caused by the separate calculation of each subsystem under the neglect of coupling in the traditional method. Additionally, it avoids the difﬁculty in numerical calculation caused by coupling, while simultaneously balancing the accuracy and efﬁciency of the model simulations.


Introduction
Remote sensing satellites are a type of satellite that utilizes optical or microwave technology to collect information about the ground or air targets.They have numerous applications in areas such as agriculture, forestry, marine, land management, environmental protection, weather forecasting, and military operations [1][2][3][4][5][6][7][8].To create a comprehensive remote sensing satellite system, it is crucial to evaluate the effectiveness of these satellites [9,10].A rational and scientific effectiveness evaluation process can serve as a valuable reference for the requirements demonstration and development of remote sensing satellites.It is conducive to optimizing the design of remote sensing satellite and improving the capability of task execution [11,12].
There are several methods for effectiveness evaluation, and the simulation statistics methods have proven to be advantageous for the quantitative evaluation of complex problems [13].Simulation statistics methods are used to evaluate effectiveness through modeling and simulation.After many years of development, this field has mainly divided into two directions: functional direct modeling and principle modeling [14].Functional direct modeling involves building models based on the functions of the object being evaluated, resulting in a simple approach.However, this method often leads to high uncertainty and low accuracy in the evaluation.In contrast, principle modeling refers to build models through the mathematical analysis of the operating principles of the system [15,16].Computational resource requirements are high, but high-fidelity simulations of modeling can be accomplished [17].In 1994, Hughes and Loral Company in the United States(US) pioneered the use of dynamic simulation technology in satellite engineering by dividing subsystems into the structure, power, thermal control, Telemetry Tracking and Command(TTC), propulsion, attitude control, and payload.They also utilized model sharing to facilitate the exchange of work results among personnel in different positions [18].In 1997, the US Department of Defense proposed the concept of virtual prototyping (VP), which aimed to map digital systems to physical systems completely [19].In 2003, GRIEVES introduced the concept of Digital Twin (DT) and, in collaboration with the Air Force Research Laboratory, defined it as a simulation process that utilizes mathematical models, sensor updates, history operation data and integrates multidisciplinary, multi-physical quantities, and multiple scales and probabilities.This process aimed to achieve a complete mapping of physical systems in the digital world, reflecting the full life cycle of real physical systems [20].Considering multi-dynamic, multi-spatial scale and multi-physics coupling in the whole life cycle of the system is the core feature of DT modeling.The concept of DT has since been widely accepted and recognized as the cornerstone of the new industrial revolution [21].
Effectiveness evaluation reflects the overall system capability of remote sensing satellites.The indicator system is a set.Some indicators are system-wide, such as mass, cost-effectiveness ratio, etc. [22].Some indicators are decomposed by satellite components, focusing on satellite subsystem capabilities.Wang Xiji [23] established an evaluation indicator system for the functions and performance of commonly used satellite platforms and payloads, which reflects the system capabilities of satellites relatively completely.In this indicator system, the remote sensing payload subsystem indicators include parameters such as time coverage, payload data rate, and image signal-to-noise ratio (SNR).The attitude and orbit control indicators include attitude measurement accuracy, attitude pointing accuracy, and attitude stability accuracy.For the power subsystem's power consumption, the maximum output power of the solar array and other indicators are used to assess the performance.The thermal control subsystem indicators include the equilibrium temperature of the entire satellite and the temperature of typical devices.The propulsion subsystem indicators include specific impulse and propellant extrusion efficiency.The TTC subsystem indicators include TTC time coverage and TTC data transmission rate.It should be pointed out that although it can be classified by subsystem, each indicator shows the performance of this subsystem in a specific system and environment, and is the ability of this subsystem and the environment outside the system to interact with each other.For instance, when calculating attitude measurement accuracy, the error characteristics of the gyroscope are significantly influenced by the rotor shaft temperature.Any changes in external temperature will directly impact the measurement accuracy.Therefore, it is imperative to include temperature model to accurately reflect the state changes of the gyroscope [24].Momentum wheels also exhibit similar characteristics, where the friction force on their shaft is influenced by temperature and affects the control torque output [25].Simply establishing a model of this subsystem cannot reflect the relationship between this subsystem and other subsystems and the environment.Therefore, when evaluating these indicators, it is necessary to establish a model of a complete satellite system and environment to reflect the relationship between each subsystem and the environment.Without considering coupling, it cannot accurately reflect the core characteristics of digital satellites.
The challenge of digital satellite coupling modeling is that the state of the satellite has multiple dynamics and multiple spatial scales [26].Multi-dynamic refers to the inconsistency of the change cycle of satellite state parameters.When using the attitude and orbit control component-level model to calculate the attitude and orbit control indicators, the simulation step size is more appropriate in milliseconds [27].In an attempt to accurately calculate power indicators, a detailed model of the power supply system is necessary, which requires data on the sunlit and Earth shadow areas of the satellite, as well as data on several orbital periods.The charging and discharging processes of the power supply are relatively slow, requiring simulation step size measured in seconds [28].Conversely, when calculating thermal control indicators, a more detailed thermal control mesh must be established, the periods of thermal decay effects on components are long and data over extended periods of several months or even years is necessary for indicator calculation, with simulation step size measured in hours or days [29].It is possible to evaluate each indicator independently with each subsystem model, but to form a complete system model, it is necessary to couple all subsystem models such as the attitude and orbit control, power and thermal control together for calculation.At this time, it is difficult to choose the simulation step size.The multi-spatial scale problems is the same.Different problems involve different spatial scales.It is easy to choose the size of the finite element if it is solved independently, but if considering the coupling, the size of the finite element varies greatly, which causes difficulties to solve the partial differential equation.
A model is a simplified representation of real-world objects, and the depth of human understanding of real things is endless.For the specific problems to be studied, infinitely fine models are not necessary, as long as they meet the needs of the research problems.When evaluating the effectiveness of remote sensing satellites, it is not to adjust the model accuracy of each subsystem synchronously with the same accuracy, but to arrange the accuracy of the simulation model of each subsystem according to the indicators to be evaluated and the coupling degree of each subsystem.A fine-grained model with high precision is selected for key subsystems, and a coarse-grained model with low precision is selected for subsystems with relatively low correlation.When calculating attitude and orbit control indicators, instead of only building a fine-grained attitude and orbit control model and ignoring other subsystem models, it is appropriate to coarsen the granularity of the power, thermal control, and telemetry subsystems.When calculating the thermal control indicators, the division of the thermal meshes needs to be detailed enough, the model of the thermal control component should be refined, and the granularity of other subsystems should be appropriately coarsened.Additionally, for the power, thermal control, propulsion, TTC, and payload subsystems, it is necessary to refine their models based on their specific goals while coarsening the models of other systems.
To address these issues, this paper proposes a multi-granularity model effectiveness evaluation method for remote sensing satellites.This paper classifies indicators of the effectiveness evaluation indicator system according to the change period, establishes simulation models with different simulation fidelity for each subsystem.For each category of indicators, a satellite model with different granularities is selected, and a multi-granularity satellite effectiveness evaluation scheme is obtained.
The rest of this paper is organized as follows.Section 2 begins with a mathematical description of the evaluation indicators, followed by multi-granularity modeling.In Section 3, the model granularity corresponding to the three types of indicators is described.In Section 4, a discussion of the simulation case and results analysis is provided.Lastly, Section 5 concludes this paper and describes future work.

Typical Satellite Effectiveness Evaluation Indicator
Time resolution refers to the minimum time step of the simulation process.The changes in physical quantities and the simulation duration required to calculate indicators will impact the selection of time resolution.In this section, the typical effectiveness indicators with classical characteristics are filtered from commonly used indicators, and are classified according to time resolution.

Millisecond-Level Time Resolution Indicator (1) Attitude measurement accuracy
The measurement accuracy of the three-axis attitude is equal to the maximum absolute value of the difference between the attitude angle and the reference attitude angle, and the real dynamic attitude angle error at all times.
where δ ϕ mesure , δ θ mesure , and δ ψ mesure represent the measurement accuracy of roll angle, pitch angle, and yaw angle, ϕ sensor , θ sensor and ψ sensor represent the roll angle, pitch angle, and yaw angle measured by the sensor, ϕ, θ and ψ are the real dynamic roll angle, pitch angle, and yaw angle.
(2) Attitude pointing accuracy The three-axis pointing accuracy is equal to the absolute value of the difference between the expected target attitude angle and the three-axis stable attitude angle.
where δ ϕ point , δ θ point , and δ ψ point represent pointing accuracy of roll angle, pitch angle, and yaw angle.
(3) Attitude stable accuracy The three-axis stability accuracy is equal to the difference between the expected target attitude angular velocity and the three-axis stable attitude angular velocity.
where ω ϕ stable , ω θ stable , and ω ψ stable represent stable accuracy of roll angular velocity, pitch angular velocity, and yaw angular velocity.

Second-Level Time Resolution Indicator
(1) Bus power supply performance where I Loop is bus current and U Loop is bus voltage.
(2) Output power of solar array where I Wing is the output current of the solar array, and U Wing is the output voltage of the solar array.
(3) Maximum output power of solar array The maximum output power of a solar array over some time (4) Maximum discharge depth of the battery where C ba is the total designed capacity of the battery, C bat_min is the minimum remaining battery power during all periods, and C bat is the total remaining battery power.
(5) Efficiency of propellant where m EOL represents the total mass of remaining propellant at the end of its lifespan, and m p is the total mass of propellant at the beginning of its lifespan.Extrusion efficiency of oxidant and fuel storage tanks are η o and η f , the m o is the total mass of initial oxidant, and m f is the total mass of fuel.
(6) Specific impulse I SP is the specific impulse of propellant, the total impulse is as follows 2.1.3.Day-Level Time Resolution Indicator (1) Subsystem thermal power The thermal power of the thermal control system is equal to the sum of the thermal power of all devices of the thermal control subsystem where P thermal is thermal power of the thermal control subsystem, D j Thermal represents the devices of the thermal control subsystem, P(D j Thermal ) is device thermal power.For any device (the number is j the temperature T Dj at a certain time point can be obtained from the thermal control subsystem model, calculate the temperature at all time points t 0 , . . ., t i , . . ., t N T to obtain the minimum temperature T min Dj and maximum temperature T max Dj . T min Dj = min T Dj (t 0 ), . . ., T Dj (t i ), . . ., T Dj (t N T ) T max Dj = max T Dj (t 0 ), . . ., T Dj (t i ), . . ., T Dj (t N T ) (2) Node temperature By calculating the temperature at all time points of all hot nodes on the satellite, the equilibrium temperature T bal , minimum temperature T min sat , and maximum temperature T max sat of the entire satellite can be obtained.The formula is as follows, T min sat = min T 0 , . . ., T j , . . ., T M , T j = min T j (t 0 ), . . ., T j (t i ), . . ., T j (t N T ) T max sat = max T 0 , . . ., T j , . . ., T M , T j = max T j (t 0 ), . . ., T j (t i ), . . ., T j (t N T ) where N node represents the total number of thermal nodes, T j is the temperature of the node j, and t i is the time point.
(3) Payload time coverage The number of coverage times is equal to the number of coverage periods in the entire simulation time.
where t 0 is the simulation start time, t i = t 0 + i∆t is the time point of the ith step of the simulation start, t N is the time point of simulation end time point, c T (t i ) is switching identification from noncoverage time period to coverage period, t i−1 is the time point of noncoverage time, and t i is the time point of observing time.
The total coverage time is equal to the sum of each coverage time The coverage rate on target time is The average coverage time can be calculated by dividing the total coverage time by the number of coverage times.
(4) Discovery probability [30] n Discovery is the total number of tasks found in an orbital period after all targets are lost where t Lost i,j is each lost time of the target and t Discovery i,j is the last recent discovery time.(5) Discovery response time [30] T Discovery = 1 where V N is the satellite ground trajectory rate, θ x is payload scanning angle width, b it is the number of bits per sampling, s am is the average number of samples per pixel, and q sam is frame efficiency.(7) Coverage of telemetry and command The number of coverage times is equal to the number of coverage time intervals in the overall simulation duration.
where t 0 is the simulation start time, t i = t 0 + i∆t is the time point of the ith step of the simulation start, t N is the time point of simulation end time point, c TTC (t i ) is switching identification from noncoverage period to the coverage period, t i−1 is the time point of noncoverage time, and t i is the time point of observing time.
The total coverage time is equal to the sum of each coverage time The coverage rate on target time is The average coverage time can be calculated by dividing the total coverage time by the number of coverage times.

Multi-Granularity Modeling
Satellite multi-granularity modeling can be divided into the payload subsystem, attitude and orbit control subsystem, power subsystem, thermal control subsystem, propulsion subsystem, telemetry tracking and command subsystem.The specific structure of multigranularity modeling is shown in Figure 1.

Payload Subsystem
The payload subsystem model can be built with two granularities: ideal optical visible granularity and resolution-constrained granularity.
The former is characterized by the angle of visibility, which is determined by the

Payload Subsystem
The payload subsystem model can be built with two granularities: ideal optical visible granularity and resolution-constrained granularity.
The former is characterized by the angle of visibility, which is determined by the target and the light shaft and can be calculated based on the location of the camera and the angle of connection between the two.If the angle is smaller than the corner of the field, it is considered to be optical visibility, otherwise, it is not visible.The latter granularity is constrained by the resolution of the payload and is determined by the smallest resolvable element in the image.
where β is the angle of the target and the light shaft and β view is the visual angle.
Resolution-constrained granularity needs to consider data generation, data compression, and data storage of the payload camera.The rectangular field of view generated by the optical camera is shown in Figure 2. S is a satellite camera, and C1C2C3C4 is the rectangular field of view.Based on the camera's field of view, the positions of the four points in a rectangle can be calculated.Images can be converted to pixel information and data after the fit resolution is selected.If the pixel data is empty, the amount of data generated is constant, otherwise, the amount of data generated is as follows, where f is the resolution of the rectangle image's long side, d is the resolution of the rectangle image's short side, h is image depth, q is the data volume of the picture, the unit is the byte.
If the sampling frequency is  , the total data volume generated in the period t is Data are stored by compression; compressibility is c and the compressed data volume is When data are transmitted to the Earth station, link transmission must be calculated.This calculation process is managed by the track telemetry and control subsystem and is shown in the telemetry and control subsystem.

Attitude and Orbit Control Subsystem
The attitude and orbit control subsystem model can be described as the payload and Based on the camera's field of view, the positions of the four points in a rectangle can be calculated.Images can be converted to pixel information and data after the fit resolution is selected.If the pixel data is empty, the amount of data generated is constant, otherwise, the amount of data generated is as follows, where f is the resolution of the rectangle image's long side, d is the resolution of the rectangle image's short side, h is image depth, q is the data volume of the picture, the unit is the byte.
If the sampling frequency is η, the total data volume generated in the period t is Data are stored by compression; compressibility is c and the compressed data volume is When data are transmitted to the Earth station, link transmission must be calculated.This calculation process is managed by the track telemetry and control subsystem and is shown in the telemetry and control subsystem.

Attitude and Orbit Control Subsystem
The attitude and orbit control subsystem model can be described as the payload and orbital kinematics constraints granularity, capability constraints of attitude and orbit control subsystem granularity, and capability constraints of attitude and orbit control devices granularity.
The payload and orbital kinematics constraints granularity is a mathematical representation that accounts for the changes in orbit and attitude kinematics parameters over time, without taking into consideration the control process.Orbit is calculated as an elliptical orbit, and the kinematics formula is The first line in the formula represents the elliptical track motion equation, E is the eccentric anomaly, M is the mean anomaly, e is eccentricity, n is the average orbital angular velocity, and τ is perigee time.
Setting attitude mode as gaze mode, the expected attitude angle of the satellite can be calculated based on the projection R T of the target.
where ϕ, θ, and ψ are roll angle, pitch angle, and yaw angle, where ψ expect is the expected yaw angle.Capability constraints of the attitude and orbit control subsystem granularity require the closed-loop control process of sensors, controllers, and actuators, as well as the consideration of the influence of force and torque on the attitude and orbit control state parameters.
Sensor model [32] is where X Sensor is measured value, X Theory Sensor is real value, X Error Sensor is measured error.The measured error can be calculated using a normal distribution.f Theory Sensor is the theoretical model of the sensor, and f Error Sensor is the error model of the sensor.X Install Sensor is install information of sensor, X Cele is a celestial state, i.e., solar vector, Earth vector.
The output command of the satellite controller [32] X Controller is where f Controller is a controller algorithm including an exclusion algorithm, filtering algorithm, attitude determination algorithm, and control algorithm.X Mode Controller is the state parameters of the satellite controller, X Expected is the expected satellite state, and X Mode Controller is related to the measurement parameters of the sensor and telecommand, where X TC is telecommand.
The expected satellite state includes default parameters for the controller and parameters set through ground telecommand, f Expected is a relation function, where X 0 Expected is the controller default parameters.The real output value of the actuator [32] is as follows where X

Theory
Actuator is the ideal output value of the actuator, and X Error Actuator is the output error of the actuator.The output value of the actuator X Actuator includes force, torque, and angular momentum.f Theory Actuator is the theoretical model of the actuator, and f Error Actuator is the error model of the actuator.X Install Actuator is install information of actuator.For fixed actuators, the install parameter X Install Actuator is invariable, only related to the initial install information X Install Fixed .
For moving actuators, the install parameter is related to the initial state and motion parameter.
where X Install Movable is the install parameter of the moving actuator, X Movable is the motion parameter of the moving actuator, f Install Movable is the kinematics rule of the moving device, and X 0 Movable_Install is the initial install information of the moving device.The changes in motion parameters are related to the initial state, control commands, or telecommand Moving devices of satellites include solar arrays, antenna, etc.For solar arrays, motion parameters X Movable mainly include the speed of the solar array .θ wing and the rotation angle of the solar array θ wing .The install parameter X Install Movable mainly includes the normal vector of the solar array W. For the antenna, the motion parameter X Movable mainly includes the rotation rate of the antenna .α Ant , .β Ant , and the rotation angle of the antenna α Ant , β Ant .The install parameter X Install Movable is the normal vector of the antenna A. The dynamic equation is as follows where F o is environmental force, F t is control force.M o is environmental torque,M t is control torque, H is the angular momentum of the satellite.The capability constraints of attitude and orbit control devices granularity further considers the coupling effect between each device and the surrounding environment based on the previous granularity.
The error model of the sensor and actuator is where Taking a laser gyro as an example, considering the coupling with thermal field and radiation, the angular velocity measurement value ω m is where λ is the wavelength, n r is the normalized refractive index of the optical path, A laser , L laser are area and perimeter enclosed by a closed optical path, ∆ f m is a measurement of the frequency difference between the front and back beams of light, where ∆ f real is the real value of frequency difference, B 0 is zero bias error, B 1 is random walk error, and B 2 is an error caused by the accumulation of radiation fluence.Temperature couples between B 0 and thermal field calculation, B 2 coupled with particle radiation.The error model is as follows where a 0 , a 1 , a 2 are zero bias compensation coefficient obtained by fitting measurement data, T laser is the temperature of a laser gyro, Ω L is lock zone threshold, Ω Dm is peak jitter rate, S K is laser gyro scale factor, J(≥ E 0 ) is radiation flux along satellite orbit, f Rad_RLG is the relationship function between error and radiation flux.
Taking the wheel as an example, considering the coupling between friction coefficient and thermal field, the torque output by the wheel in the body coordinates is where M M is output torque of the motor, M f 0 is bearing static friction torque, K f is the frictional coefficient, ω w is the speed of the wheel, and r w is the unit installation vector of the wheel.The output torque is as follows, where K M is torque voltage ratio coefficient, V C is the input voltage.Temperature coupling between K f and thermal field calculation.The model is as follows where ρ w is the density of lube, G w is the amount of lube, C w is the specific heat of lube, ∆T w is the difference between temperature and nominal temperature of the wheel, and D w is the bearing diameter of the wheel.Taking into account the friction factor of bearings, the angular acceleration of the wheel is where I w is the inertia of the wheel.Device-level performance indicators must not only account for the impact of intersubsystem interactions on device performance but also consider the varying degrees of impact that these interactions have on different parts of the device.For instance, the measurement accuracy of gyroscopes is affected by shaft temperature, with higher temperatures leading to greater errors.Therefore, it is important to consider the calculation of the bearing temperature field.However, the temperature of other parts of the satellite has a smaller impact, and therefore a multi-scale approach with local refinement and overall coarseness should be adopted when dividing the mesh.

Power Subsystem
The power subsystem is characterized by a set of parameters including the output of the solar array, load power consumption, battery charging and discharging, and remaining capacity.The calculation model of each parameter is as follows (1) Output model of the solar array The output model of a solar array can be modeled with constant output power granularity and considering solar array state granularity.
Constant power granularity refers to the fact that the solar array generates zero electricity in the shadow area and the output power is calculated at a constant value in the illumination area.
The considering solar array state granularity refers to the relationship between the output power of the solar array and the satellite's state, taking into account the material characteristics of the solar array and the incidence angle of solar optics.
where C shadow is shadow area identification, S 0 is the power generated by sunlight, A sa is an area of the solar array, β P is the power temperature coefficient of the solar array, X Wing is other parameter, θ SW is the normal vector of the solar array and included angle of solar vector, ∆T wing is difference between working temperature and standard temperature.θ SW is not only related to the position and attitude of the satellite but also to the relative angle θ Wing of the sail relative to the body.
(2) Load power consumption The load power can be built with three granularities: constant load granularity, considering device switch granularity, and considering device state granularity.
The power consumption in the constant load granularity is calculated using fixed values, which can be either the rated power or the average power consumption over a specified period.
The considering device switch granularity involves the real-time statistical analysis of the switch state of each device, as well as the power consumption of each device, which is constant and remains the same over time.
where N PD is the total number of electrical devices, S i Device is the switch state of the device i, power on is 1, power off is 0 and P D i is the power of the device i.
The considered device state granularity refers to the relationship between the actual power consumption and the device's working state.For example, this includes the power consumption and speed of wheels.The formula is as follows where I wheel is the inertia of the wheel, and ω wheel is the angular velocity of the wheel, α wheel is the angular acceleration of the wheel.
(3) Battery charging and discharging The battery charging and discharging model includes ideal battery granularity and nonideal battery granularity.
In the ideal battery granularity model, when the power generated by the solar array exceeds the electrical power of the equipment and the battery is not fully charged, the battery is charged at constant power, with any excess electricity dissipated through the charging regulator.Conversely, when the Earth shadow or power supply of the solar array is insufficient, the battery is discharged at constant power.
The consumption power of the charging regulator is as follows where P MaxCharge is the maximum charging power of the battery.
where P Load is load power, Q MinBattery is minimum battery capacity, and Q MaxBattery is maximum battery capacity.Battery remaining capacity is as follows where ∆t is the simulation step.The power relationship in nonideal battery granularity is unaltered, but varying charging and discharging coefficients must be taken into account when calculating power consumption.
where k Charge is charging coefficients, k Discharge is discharging coefficients.

Thermal Control Subsystem
The thermal balance expression of satellites in space is where q 1 is the direct solar radiation heat absorbed by satellites, q 2 is the Earth infrared radiation heat absorbed by satellites, q 3 is the infrared radiation heat absorbed by the satellite, q 4 is space background heating amount, q 5 is the heat generated by the satellite, q 6 is a change in internal energy of satellites, q 7 is the heat emitted by satellites into space.Due to the low temperature of the spatial background and the small heating heat, the heating amount of the spatial background q 4 can be ignored.
The temperature field of the satellite is solved using the finite element method, and the thermal balance equation of each finite element node is [33] T where T i is node temperature, m i is node mass, c i is node specific heat capacity, q i S is solar radiation heat, q i E is Earth radiation heat, q i ER is Earth reflection radiation heat,q i Rad is heat radiation, q i Cond is heat conduction, P i is internal thermal power, q i OS is external surface heat dissipation, q i IS is internal surface heat dissipation, T 0 i is node initial temperature.The thermal control subsystem model includes three parts: finite element partitioning calculation, external heat flow calculation, and internal heat flow calculation.
(1) Finite element partition calculation model The heat balance equation for each node remains constant regardless of the specific research object, but the number of nodes distributed throughout the satellite may vary depending on the research object.The thermal balance equations of each node remain unchanged, but the division of nodes in different parts of satellites can differ based on the research object.Taking Figure 3 as an example, the satellite structure mainly includes the main body, solar array, and various devices.When temperature is not the primary factor affecting certain indicators, the temperature of each part can be considered at a coarser level of detail.Typically, all structures and devices can be divided into meshes with the same level of accuracy.However, when calculating indicators that involve the temperature characteristics of a specific part, a finer mesh should be used for that part, while other parts can still use a coarser granularity to increase computational efficiency.All meshing is performed through Ansys software.q is internal surface heat dissipation, The thermal control subsystem model includes three parts: finite element partitioning calculation, external heat flow calculation, and internal heat flow calculation.
(1) Finite element partition calculation model The heat balance equation for each node remains constant regardless of the specific research object, but the number of nodes distributed throughout the satellite may vary depending on the research object.The thermal balance equations of each node remain unchanged, but the division of nodes in different parts of satellites can differ based on the research object.Taking Figure 3 as an example, the satellite structure mainly includes the main body, solar array, and various devices.When temperature is not the primary factor affecting certain indicators, the temperature of each part can be considered at a coarser level of detail.Typically, all structures and devices can be divided into meshes with the same level of accuracy.However, when calculating indicators that involve the temperature characteristics of a specific part, a finer mesh should be used for that part, while other parts can still use a coarser granularity to increase computational efficiency.All meshing is performed through Ansys software.(2) External heat flow calculation [33] The external heat flow calculation model can be described as constant external heat flow granularity and node state external heat flow granularity.
The constant external heat flow granularity does not account for local temperature differences, and the radiation parameters are calculated using average values in the local film area and the light area, respectively.(2) External heat flow calculation [33] The external heat flow calculation model can be described as constant external heat flow granularity and node state external heat flow granularity.
The constant external heat flow granularity does not account for local temperature differences, and the radiation parameters are calculated using average values in the local film area and the light area, respectively.
The node state external heat flow granularity is a model that considers the orbit, attitude, and other states that need to be accounted for in the calculation of external heat flow, as well as material attributes.
Node i solar radiation heat is where α Si is solar absorption rate on the outer surface of satellites, F s is solar constant, A i is the surface area of node i, and φ 1i is the direct solar radiation angle coefficient of the node i.
where dot(r s , n i ) is the dot product of the vector r s and vector n i , C shadow is the shadow zone identification, r s is a solar vector, and n i is an illuminated surface normal vector.Earth radiation heat of node i is as follows where ε ei is the outer surface emission rate of the node i, E e is the average of surface infrared radiation density, and φ 3i is the Earth's infrared radiation angle coefficient of the node i. Define k = R e /(R e + h), where h is satellite orbital altitude, R e is the average radius of the Earth.When 0 Earth reflection radiation heat of node i is where E r is the average reflection density of the Earth's surface facing the solar radiation, and φ 2i is the Earth's infrared radiation factor of the node i.
where r e is the Earth vector.
(3) Internal heat flow calculation [33] The internal heat flow calculation model can be described as constant internal heat flow granularity and node state internal heat flow granularity.
where m 0 f uel is initial propellant mass in the tank, ∆m f uel is mass consumption of propellant, η L is propellant consumption rate, N Fire is ignition command.
Constant granularity refers to the calculation of specific impulse and fuel consumption rates based on constant values.
The constant pressure model granularity is concerned with the relationship between specific impulse, fuel quality, and the working status of temperature, pressure, and valve switching duration.
where µ f is flow coefficient, A F is the total area of the thruster injection hole, g is the gravitational acceleration, and γ f is the relative density of the propellant liquid.
The storage tank maintains constant gas pressure, that is P gas = P 0 gas .The gas mass in the tank can be obtained from the ideal gas equation.m gas = P gas V gas M gas /R gas T gas (87 where M gas is the average molar mass of gas, and R gas is the universal gas constant. In the depressurization model granularity, the working status of the specific impulse, fuel quality consumption and temperature, pressure, and valve switch duration is related.The mass of gas in the storage tank is constant, that is m gas = m 0 gas .The gas pressure in the storage tank can be calculated through the ideal gas equation P gas = m gas R gas T gas /V gas M gas (88) 2.2.6.TTC Subsystem TTC subsystem includes communication visible constraint granularity, Earth occlusion granularity, and SNR transmission granularity.
The communication visible constraint granularity focuses on the connection state of the satellite and earth radar station and does not consider the specific transmission process.
The Earth occlusion granularity focuses on geometric constraints.When the satellite and the Earth station are not obstructed by the Earth, their transmission ability is considered connecting, otherwise disconnected.The elevation between the satellite and the Earth station can be calculated.The elevation angle is the angle formed by the horizon horizontal line where the ground station is located at the center line of the antenna.When the actual elevation angle ε is greater than the minimum elevation angle ε min , it is determined that the two can communicate.Assuming that the position vector of satellites and Earth stations is R s and R e , the included angle R s and R e is α.
Transmission capacity is In the SNR transmission granularity, the transmission capacity is influenced by factors such as the relative position of the receiver and transmitter, performance, data compression, and transmission loss, the basic transfer equation is as follows where P Receiver is the receiving antenna power, P Sender is the transmitter power, G Sender is the transmitter antenna gain, G Receiver is receive antenna gain, L Space is the path loss, L x is other various losses.
The unit of operating frequency f is MHz, and The distance d between two interfaces is measured in kilometers.L x mainly includes the influence of antenna pointing.Loss of transmitter antenna pointing is L s , the loss of receiver antenna pointing L r .The loss is related to the directional pattern of the antenna, A computational model for a point beam antenna is as follows where θ s is the angle between the Z-axis of the transmitting antenna instrument coordinate system and the vector between the transmitting node and the receiving node.
The current carrier-to-noise ratio from sender to receiver ] min , the sender from the receiver is normal, otherwise it will not be connected.

Selection of Effectiveness Evaluation Indicator Models for Remote Sensing Satellites
Capitalizing on the sparse nature of the satellite system, the effectiveness evaluation metrics are categorized into three groups based on time intervals: milliseconds, seconds, and days.As depicted in Figure 1, the satellite architecture consists of multiple subsystems, including payload, attitude and orbit control, power, thermal control, propulsion, and TTC.Each subsystem has various levels of granularity, as illustrated in Figure 4.When assessing each group of indicators, it is recommended to utilize the corresponding granularity level shown in Figure 4.
The millisecond-level indicator is primarily used to evaluate the performance of attitude control, requiring the use of finer simulation granularity for attitude sensors, controllers, and actuators.In contrast, other subsystems use coarser granularity.For instance, the payload subsystem selects the ideal optical visible granularity, while the TTC subsystem selects the Earth occlusion granularity.The attitude and orbit control subsystem selects the capacity constraints of attitude and orbit control devices granularity.The power subsystem selects the constant granularity.As for the thermal control subsystem, the device mesh should be finely divided, while the structures such as the body and solar array should be coarser.The external and internal heat flows should be calculated using node state granularity.The propulsion subsystem selects constant granularity.When establishing simulation condition, it is recognized that the overall temperature change cycle is influenced by the annual cycle, which is much greater than the time required for millisecond-level simulations.Consequently, the two stages of local shadow and light are based on the two seasons of winter and summer, with simulation time spans of several thousands of seconds and step intervals of 0.1 s. metrics are categorized into three groups based on time intervals: milliseconds, seconds, and days.As depicted in Figure 1, the satellite architecture consists of multiple subsystems, including payload, attitude and orbit control, power, thermal control, propulsion, and TTC.Each subsystem has various levels of granularity, as illustrated in Figure 4.When assessing each group of indicators, it is recommended to utilize the corresponding granularity level shown in Figure 4.The millisecond-level indicator is primarily used to evaluate the performance of attitude control, requiring the use of finer simulation granularity for attitude sensors, controllers, and actuators.In contrast, other subsystems use coarser granularity.For instance, the payload subsystem selects the ideal optical visible granularity, while the TTC subsystem selects the Earth occlusion granularity.The attitude and orbit control subsystem selects the capacity constraints of attitude and orbit control devices granularity.The power subsystem selects the constant granularity.As for the thermal control subsystem, the device mesh should be finely divided, while the structures such as the body and solar array should be coarser.The external and internal heat flows should be calculated using node state granularity.The propulsion subsystem selects constant granularity.When establishing simulation condition, it is recognized that the overall temperature change cycle is influenced by the annual cycle, which is much greater than the time required for millisecond-level simulations.Consequently, the two stages of local shadow and light are based on the two seasons of winter and summer, with simulation time spans of several thousands of seconds and step intervals of 0.1 s.
The second-level indicator is primarily used to evaluate the performance of power and propulsion.The power subsystem selects the considering state granularity, while the propulsion subsystem selects the depressurization mode granularity.Other subsystems use coarser granularity, such as the payload subsystem selecting the ideal optical visible granularity, the TTC subsystem selecting the Earth occlusion granularity, and the attitude and orbit control subsystem selecting the capability constraints of attitude and orbit control subsystem granularity.The mesh of the thermal control subsystem can be broadly classified into two categories, and internal and external heat flows are computed using constant values.To collect simulation statistics, the light area and the shadow area are selected.The integration duration for each stage is set at several thousands of seconds, with a step size of 1 s.The second-level indicator is primarily used to evaluate the performance of power and propulsion.The power subsystem selects the considering state granularity, while the propulsion subsystem selects the depressurization mode granularity.Other subsystems use coarser granularity, such as the payload subsystem selecting the ideal optical visible granularity, the TTC subsystem selecting the Earth occlusion granularity, and the attitude and orbit control subsystem selecting the capability constraints of attitude and orbit control subsystem granularity.The mesh of the thermal control subsystem can be broadly classified into two categories, and internal and external heat flows are computed using constant values.To collect simulation statistics, the light area and the shadow area are selected.The integration duration for each stage is set at several thousands of seconds, with a step size of 1 s.
The day-level indicator is mainly used to evaluate load observation, TTC performance, and thermal control performance.The payload subsystem selects the resolution-constrained granularity, the TTC subsystem selects the SNR transmission granularity.For the thermal control subsystem, the devices and solar array are finely divided, while the body is coarser.The calculation of internal and external heat flow selects node states granularity.The granularity of other subsystems is relatively coarse.The attitude and orbit control subsystem selects the payload and orbital kinematics constraints granularity, the power subsystem and the propulsion subsystem select the constant granularity.When setting up the working conditions, the integration duration for each stage is about one year, and the integration step is 100 s.
When calculating indicators, the granularity selected by each model is not fixed, but selected according to a certain probability.Figure 4 shows the maximum possible recommended particle size.Given the unique characteristics of individual satellites, a tailored approach is necessary to identify the most suitable combination of indicators and model granularities.By implementing a multi-granularity model, iterative evaluation processes can be employed to refine simulation granularity at each iteration.In practical engineering applications, fine-tuned adjustments can be made based on evaluation results to achieve an optimal balance between computing accuracy and resource utilization.This approach enables the optimization of indicator grouping and model granularities for efficient performance assessment.

Simulation Case 4.1. Millisecond-Level Time Resolution Indicator Calculation
The orbital elements are as follows: the semi-major axis is 6,994,596.133m, the eccentricity is 0.0001537, the inclination is 97.8416 • , the right ascension of ascending node is 155.985• , the argument of perigee is 92.481 • , and the true anomaly is 0 • .
The millisecond-level indicators mainly include the attitude and orbit control subsystem.A comparison between the coupled model and the uncoupled model of attitude orbit control subsystem should be studied.The attitude angles of the attitude and orbit control subsystem without considering thermal coupling, considering the coarse mesh and considering the multi-scale mesh, are compared, as shown in Figure 5.When the coupling effect is not considered, the noise model of the device is a normal distribution model and the attitude angle changes randomly around a constant value.Considering the influence of temperature on attitude and orbit control devices, when the mesh is divided roughly, the device temperature changes with the temperature of the whole satellite over a long period.When the mesh is divided into different scales and the local temperature is considered, the components are more affected by the local temperature, and the temperature changes quickly there.The impact on the attitude control effect is also faster.
enables the optimization of indicator grouping and model granularities for efficie formance assessment.

Millisecond-Level Time Resolution Indicator Calculation
The orbital elements are as follows: the semi-major axis is 6,994,596.133m, the tricity is 0.0001537, the inclination is 97.8416°, the right ascension of ascending n 155.985°, the argument of perigee is 92.481°, and the true anomaly is 0°.
The millisecond-level indicators mainly include the attitude and orbit control s tem.A comparison between the coupled model and the uncoupled model of attitud control subsystem should be studied.The attitude angles of the attitude and orbit subsystem without considering thermal coupling, considering the coarse mesh an sidering the multi-scale mesh, are compared, as shown in Figure 5.When the co effect is not considered, the noise model of the device is a normal distribution mod the attitude angle changes randomly around a constant value.Considering the in of temperature on attitude and orbit control devices, when the mesh is divided ro the device temperature changes with the temperature of the whole satellite over period.When the mesh is divided into different scales and the local temperature sidered, the components are more affected by the local temperature, and the temp changes quickly there.The impact on the attitude control effect is also faster.The results are shown in Table 1.According to the indicator calculation results, it can be seen that the attitude measurement accuracy, attitude pointing accuracy, and attitude stable accuracy of the shadow area are generally lower than those of the light area.This is because the working temperature of the lighting area devices is close to the nominal temperature of 20 • C, while the ground shadow area temperature is around −20~−40 • C, far deviating from the nominal temperature.Therefore, the measurement error of the sensor is greater than that of the lighting area.The calculation results of the finest granularity model and the multi-granularity model are close.Comparing the simulation time of the three types of models, the coarse granularity model takes 93 s, the multi-granularity model takes 262 s, and the fine granularity model takes 1132 s.The calculation time of the multi-granularity model is much shorter than that of the finest granularity model, but the calculation accuracy is the same.

Second-Level Time Resolution Indicator Calculation
The orbit parameter selections are the same as above.The second-level indicators mainly include the power subsystem and propulsion subsystem.A comparison between the coupled model and uncoupled model should be studied.The parameters of the power subsystem and propulsion subsystem are compared in the following figures .As shown in the figures, comparing the state variables of the power subsystem under the granularity of the constant granularity, considering device switch granularity, and considering device state granularity, when only the average power generation, power consumption, and power storage are considered, the state is a straight line, and the switching granularity is considered to be stepped.Neither of these two models can fully reflect the real state change, but considering device state granularity can better reflect the coupling change between the working state of the device and the electric quantity.The propulsion subsystem models of constant value granularity, constant pressure mode granularity, and depressurization mode granularity also reflect the above rules.
state is a straight line, and the switching granularity is considered to be stepped.Neither of these two models can fully reflect the real state change, but considering device state granularity can better reflect the coupling change between the working state of the device and the electric quantity.The propulsion subsystem models of constant value granularity, constant pressure mode granularity, and depressurization mode granularity also reflect the above rules.of these two models can fully reflect the real state change, but considering device state granularity can better reflect the coupling change between the working state of the device and the electric quantity.The propulsion subsystem models of constant value granularity, constant pressure mode granularity, and depressurization mode granularity also reflect the above rules.Comparing the indicator calculations of the three models for the power subsystem and the propulsion subsystem, they are the coarse granularity model, the multi-granularity model, and the finest granularity model.The attitude and orbit control subsystem of the coarse granularity model is the subsystem capability constraint granularity, the power and propulsion are constant granularity, and the mesh division is coarse.The attitude and orbit control subsystem of the multi-granularity model is the capability constraints of attitude and orbit control subsystem granularity, the power subsystem is the granularity of considering the state changes of the devices, and the propulsion subsystem is the depressurization mode, and the mesh is divided into multiple granularities according to different parts.The finest granularity model: the attitude and orbit control subsystem is the capability constraint of the attitude and orbit control devices granularity, the power subsystem is the granularity of considering device state changes and the propulsion subsystem is depressurization mode granularity.The mesh division is relatively fine.
The working state was configured at UTC 12:40:50 on 21 June 2022.The simulation step was 1 s and the total simulation time was 2000 s.The results are shown in Table 2.
The simulation time of the three types of models is 154 s for the coarse granularity model, 469 s for the multi-granularity model, and 2018 s for the finest granularity model.The calculation time of the multi-granularity model is much shorter than that of the finest granularity model, but the calculation accuracy is the same.The orbit parameter selections are the same as above.
The day-level indicators mainly include the thermal control subsystem and TTC subsystem.A comparison between coupled model and uncoupled model should be studied.The parameters of the thermal control and TTC subsystem are compared in the following figures (Figures 11-15).Comparing the state variables of the thermal control subsystem under the constant granularity, considering device switching granularity and considering the device state granularity, it can be seen that the change law is similar to that of the power subsystem in Figures 11 and 12   Comparing the indicator calculation results of the three models, the coarse granularity model, multi-granularity model, and the finest granularity model, the difference between the models is that the mesh division of the coarse granularity model is rough, the  Comparing the indicator calculation results of the three models, the coarse granularity model, multi-granularity model, and the finest granularity model, the difference between the models is that the mesh division of the coarse granularity model is rough, the attitude and orbit control subsystem uses payload and orbital kinematics constraints gran- stored on the satellite continues to increase.Considering the visible granularity of the target is to take pictures when the target is within the field of view of the optical axis, and not to take pictures for the rest of the time.The finest granularity not only considers whether the target is visible, but also pictures can be compressed and downloaded to the ground, so the amount of data stored in it is much smaller than the previous two and can be cleared.As shown in Figures 14 and 15, the finer the granularity of the TTC models, the more constraints they are subject to, and the more stringent the conditions for communication with ground stations.When there is no constraint, the satellite, and the ground station can communicate all the time, and the SNR is constant.When considering the Earth occlusion granularity, only when the communication link between the two is not occluded by the Earth can they have communication capability, and the SNR is also constant.In the full constrain state, the SNR varies with the path conditions such as the distance between the satellite and the ground station.When it is less than a certain threshold, even if the two are on the same side of the Earth, they may not be connected.Comparing the indicator calculation results of the three models, the coarse granularity model, multi-granularity model, and the finest granularity model, the difference between the models is that the mesh division of the coarse granularity model is rough, the attitude and orbit control subsystem uses payload and orbital kinematics constraints granularity.
The mesh of the multi-granularity model is divided according to the size requirements of devices and structures, the attitude and orbit control subsystem in the multi-granularity model uses the payload and orbital kinematics constraints model.All meshes in the finest granularity model is very fine, and other systems such as attitude and orbit control subsystem are the finest granularity.
The working state was configured at UTC 12:40:50 on 21 June 2022.The simulation step was 100 s and the total simulation time was 1 year.The results are shown in Table 3.As shown in the table, as the life of the satellite decreases, the heat generated by the devices increases, increasing the temperature indicators of the satellite.The simulation time of the three types of models is 2986 s for the coarse granularity model and 5823 s for the multi-granularity model.The simulation time of the finest granularity model is slower than the real-time, and it is unrealistic to use this model for calculation.
In Figure 13, optical payload and the TTC models of different granularity are compared.The uncoupled optical model can always take pictures of the target, and the data stored on the satellite continues to increase.Considering the visible granularity of the target is to take pictures when the target is within the field of view of the optical axis, and not to take pictures for the rest of the time.The finest granularity not only considers whether the target is visible, but also pictures can be compressed and downloaded to the ground, so the amount of data stored in it is much smaller than the previous two and can be cleared.
As shown in Figures 14 and 15, the finer the granularity of the TTC models, the more constraints they are subject to, and the more stringent the conditions for communication with ground stations.When there is no constraint, the satellite, and the ground station can communicate all the time, and the SNR is constant.When considering the Earth occlusion granularity, only when the communication link between the two is not occluded by the Earth can they have communication capability, and the SNR is also constant.In the full constrain state, the SNR varies with the path conditions such as the distance between the satellite and the ground station.When it is less than a certain threshold, even if the two are on the same side of the Earth, they may not be connected.
Different models of payload and TTC indicators are calculated as follows.The ideal visibility model is that the optical payload is always visible to the target, and the satellite communicates with the ground station in real-time.Both the optical load and the communication occlusion of the Earth occlusion model consider the occlusion effect of the Earth.The finest model is the optical load considering compressibility and transmission, and the TTC subsystem considering the SNR ratio transmission granularity.The results are shown in Table 4.It can be seen from the table that as the granularity of the model becomes finer and more constraints are considered, the probability of target discovery will decrease, and the longer the target response time is, the lower the measurement and control coverage rate will be.

Conclusions
Reasonable and scientific performance evaluation provides a powerful reference for the demonstration and development of remote sensing satellite systems.Aiming at the problem of coupling multi-dynamic and multi-spatial scales in a satellite digital twin model in the effectiveness evaluation of remote sensing satellites, this paper proposes a method to calculate the effectiveness indicators by using multi-granularity modeling.The multigranularity model in this paper takes into account the coupling effects of device states and attitude and orbit control, power, thermal control, and other subsystems.From the simulation results, it can be seen that the change in parameters of the uncoupled model has nothing to do with its own working state over time, and the component running states have no effect on the change in model parameters, while the trend of the parameters of the coupled model over time is closer to reality.For the millisecond-level, second-level, and daylevel indicators, the simulation results of the multi-granularity model evaluation indicators compared with the single-granularity evaluation indicators show that the calculation accuracy of the multi-granularity model indicators is much higher than that of the coarsegrained model, and is close to the accuracy of the most fine-grained model.Multi-grained models run much faster than fine-grained models, but slower than coarse-grained models.In particular, for the simulation of the thermal control model, even the simulation time of the fine-grained model is slower than the real time, which loses the meaning of simulation.Compared with the single-granularity model, the accuracy of the evaluation results of the multi-granularity model meets the requirements and the calculation efficiency is higher, which verifies the feasibility and effectiveness of the multi-granularity model.
In order to strengthen the evaluation of multi-granularity modeling and its applicability, looking toward to the future, research can be carried out from two aspects.According to the missions of different satellites, multi-granularity models of different types of satellites will be constructed to complete the performance evaluation.For high real-time evaluation requirements, the machine learning evaluation model will be trained using the simulation data of the multi-granularity model to achieve rapid evaluation.

Figure 1 .
Figure 1.The tree diagram of multi-granularity modeling.

Figure 1 .
Figure 1.The tree diagram of multi-granularity modeling.

Figure 2 .
Figure 2. Optical area array camera rectangular viewing diagram.

Figure 2 .
Figure 2. Optical area array camera rectangular viewing diagram.

Figure 3 .
Figure 3. Multi-scale division of finite element mesh.

Figure 3 .
Figure 3. Multi-scale division of finite element mesh.

Figure 4 .
Figure 4. Granularity combination recommended for effectiveness evaluation indicators calculation.

Figure 4 .
Figure 4. Granularity combination recommended for effectiveness evaluation indicators calculation.

Figure 5 .
Figure 5.Comparison diagram of attitude angle measurement for different models.The calculation of indicators is verified under different granularity models and different working conditions.The models are the coarse granularity model, the multi-granularity model, and the finest granularity model.The coarsest granularity model of attitude and orbit control subsystem considers the device level, but the mesh is divided coarsely, which only reflects the temperature change of the whole satellite and cannot reflect the temperature of local devices.The multi-granularity model considers the characteristics of the devices, including the influence of local temperature on the noise.The mesh of the temperature of the device is finer, and the body and solar array are rough.The finest granularity model considers the characteristics of devices, and the mesh division of devices, body, and sailboard is very fine.The working states are statistically analyzed in four time periods: UTC 12:40:50 on 21 June 2022 (Summer Shadow Zone), 13:50:50 on 21 June 2022 (Summer Light Zone), 13:29:40 on 21 June 2022 (Winter Shadow Zone), and 14:25:50 on 21 December 2022 (Summer Light Zone).The simulation step is 0.1 s and the total simulation time is 1600 s.The results are shown in Table1.

Figure 8 .
Figure 8.Comparison diagram of remaining battery power.

Figure 8 .
Figure 8.Comparison diagram of remaining battery power.

Figure 8 .
Figure 8.Comparison diagram of remaining battery power.

Figure 9 .
Figure 9.Comparison diagram of fuel mass.

Figure 10 .
Figure 10.Comparison diagram of Helium pressure in the fuel tank.

Figure 9 .
Figure 9.Comparison diagram of fuel mass.

Figure 8 .
Figure 8.Comparison diagram of remaining battery power.

Figure 9 .
Figure 9.Comparison diagram of fuel mass.

Figure 10 .
Figure 10.Comparison diagram of Helium pressure in the fuel tank.

Figure 10 .
Figure 10.Comparison diagram of Helium pressure in the fuel tank. .

32 Figure 11 .
Figure 11.Comparison diagram of thermal power of devices.

Figure 12 .
Figure 12.Comparison diagram of momentum wheel temperature.

Figure 11 . 32 Figure 11 .
Figure 11.Comparison diagram of thermal power of devices.

Figure 12 .
Figure 12.Comparison diagram of momentum wheel temperature.

Figure 12 .
Figure 12.Comparison diagram of momentum wheel temperature.

Figure 13 .
Figure 13.Changes in data stored in different models.(a) Uncoupled model and optical visible model.(b) Considering data compression and transmission models.

Figure 15 .
Figure 15.Comparison diagram of SNR changes of different models.Different models of payload and TTC indicators are calculated as follows.The ideal visibility model is that the optical payload is always visible to the target, and the satellite communicates with the ground station in real-time.Both the optical load and the communication occlusion of the Earth occlusion model consider the occlusion effect of the Earth.The finest model is the optical load considering compressibility and transmission, and the TTC subsystem considering the SNR ratio transmission granularity.The results are shown in Table4.

Figure 14 . 32 Figure 14 .
Figure 14.Comparison diagram of link state changes in different models.

Figure 15 .
Figure 15.Comparison diagram of SNR changes of different models.Different models of payload and TTC indicators are calculated as follows.The ideal visibility model is that the optical payload is always visible to the target, and the satellite communicates with the ground station in real-time.Both the optical load and the communication occlusion of the Earth occlusion model consider the occlusion effect of the Earth.The finest model is the optical load considering compressibility and transmission, and the TTC subsystem considering the SNR ratio transmission granularity.The results are shown in Table4.

Figure 15 .
Figure 15.Comparison diagram of SNR changes of different models.
Actuator are coupling state variables of the sensor and actuator.There are mainly two types of state variables: one is transient state variables, which can be updated immediately after time, such as voltage and switch state.The other one is the steady-state state variable, which changes relatively slowly over time.
couplingSensor and X coupling

Table 1 .
Comparison of millisecond-level time resolution indicators.

Table 4 .
Payload and communication link indicator statistics.