# Method for Predicting the Liquid-Sloshing Characteristics of a Microsatellite Propulsion System

Authors: L. Teng; T. Meng; Zh. H. Jin
Journal: International Journal of Aerospace Engineering
 (2019)

---

## Abstract

In this study, a method for establishing a prediction model for the liquid-sloshing characteristics of microsatellite propulsion systems is proposed. The ideal states of the liquid-sloshing characteristics are considered for the following three models: (1) storage tank, (2) coupled storage tank and microsatellite, and (3) coupled storage tank and microsatellite with a deployment mechanism. The smoothed-particle hydrodynamics method is implemented in ABAQUS to study the effect of sloshing on a storage tank and microsatellite disturbance for the above three cases. Relationship models between the sloshing time, sloshing angular velocity, amount of fluid filled, and satellite attitude angular velocity are established. The analysis results show that the disturbance angular velocity of the empty storage tank and sloshing angular velocity have a linear relationship. Furthermore, the disturbance angular velocity of the liquid-filled storage tank exhibits a surface relationship with the sloshing angular velocity and amount of fluid filled in a three-dimensional coordinate system. Additionally, the disturbance angular velocity in the liquid-filled state of the storage tank is higher than that of the empty storage tank, and the degree of disturbance decreases with the increase in the mass of the coupling. The disturbance of the storage tank and microsatellite with the deployment mechanism is 10−2°/s when the angular velocities of the satellite are ωx′=3, ωy′=3, and ωz′=2. The maximum and minimum deviations between the calculation and simulation results of the three models are 7.6% and 1.1%, respectively. The model is used to predict the disturbance angular velocity of the microsatellite. When the calculation results of the model are compared with the orbit satellite data, the maximum and minimum disturbance angular velocity deviations occur in the y and z directions with a deviation of 43.36% and 14.86%, respectively. This demonstrates the accuracy of the analysis and model. The results of this study can provide theoretical guidance for the engineering design and attitude and orbital control of a microsatellite propulsion system.


---

## Body

## 1. Introduction

The continuous improvement of modern technology has led to several important developments such as microelectronics, microcomputers, and new material development. Although large-satellite technology is advancing, the function and integration of microsatellites have gained increasing attention. Research in the field of space pertaining to modern microsatellites is rapidly expanding owing to several attractive features, including their low weight, small size, low cost, high performance, and short development period. Presently, the development of microsatellites is focused on network applications. Satellite formations or constellations consist of several small satellites, and thus, each satellite requires precise control of its orbit phase and attitude. In addition, because microsatellites are being used for increasingly difficult tasks, their pointing accuracy and precision have improved at a rapid rate. Furthermore, the demand for in-orbit propulsion technology of microsatellites has heightened along with that for miniaturized propulsion systems. The propulsion systems of microsatellites have the characteristics of high-volume utilization, high efficiency, adept control, and less environmental pollution. Specifically, the structure of a liquid micropropulsion system is relatively simple and has the desirable characteristics of low power consumption and high specific impulse. Thus, liquid micropropulsion systems have become the most common microsatellite propulsion systems.The liquid propulsion system of a microsatellite is an important component [1] having a direct effect on the operating life, reliability, orbit and attitude control, maneuverability, and position retention of the microsatellite. With the development of microsatellites, users are demanding improved satellite attitude and orbit control. However, instantaneous attitude and orbit adjustments lead to liquid sloshing in the storage tank of a satellite propulsion system. The dynamic characteristics of a storage tank will affect those of the satellite [2, 3], which in turn affect the stability and control precision of the motion system. Studying the sloshing characteristics of a liquid storage tank also requires understanding the effect of the liquid sloshing on the attitude and tracking control of the storage tank structure, which will assist in improving the control precision and capability.The occurrence of liquid sloshing in storage tanks is a concern across many fields such as aerospace, nautical science, and automotive engineering. Understanding the sloshing characteristics of storage tanks is essential for scientific research and engineering design. The smoothed-particle hydrodynamics (SPH) method is a Lagrange-type mesh-less method used to numerically simulate fluid dynamics. This method uses the kernel function to approximately describe the differential equation. There is no grid relation between the particles, and the value of a given point is locally approximated by its neighborhood node. This method is suitable for dealing with large deformations and fluid–solid coupling problems. Ni et al. [4] and Liu and Chang [5] studied the SPH method, and the findings showed that regular, uniform particle distributions yield more accurate results and increased the stability of the numerical calculations. Gingold and Monaghan [6] and Lucy [7] first proposed the SPH method and discussed its application in astrophysics research. The SPH method has been gradually applied to studying hydrodynamics and has shown significant advantages. Monaghan [8] was the first to apply the SPH method to solve a fluid dynamics problem, and since then, the method has been rapidly developed and applied. Additionally, the increase in the use of the SPH method has motivated the development of other hydrodynamic methods and theories.Jiang et al. [9] studied the liquid-sloshing characteristics of an elastic rectangular tank at different horizontal excitation frequencies and obtained the free-surface height and pressure change curves of the tank. Liu et al. [10] improved the SPH method to numerically simulate the liquid sloshing in a rhombus liquid tank at different filling depths, external frequencies, and amplitudes and analyzed strong nonlinear phenomena such as rolling and breaking. Huang et al. [11] studied the feasibility of implementing an SPH algorithm to analyze liquid sloshing and the effects of liquid volume, sloshing angle, liquid material property, and cycle. Additionally, the filling depth of the liquid tank under the maximum impact load was obtained. Hao et al. [12] used the SPH method to simulate the sloshing of the fuel in a missile tank under pitch and roll excitation. Furthermore, the impact force and impact torque produced by the sloshing liquid were determined, and the effects of the change law of the center of gravity of the liquid on the stability of the missile were examined. Cardoso-Ribeiro et al. [13] proposed port-Hamiltonian systems (PHs) for formulating the liquid-sloshing equations of moving containers and beam structural equations using piezoelectric actuators. The PH model proved useful for designing an active control law for the reducing the sloshing phenomenon. Hasheminejad and Soleimani [14] investigated using modified cylindrical Bessel functions the free three-dimensional (3D) sloshing in a finite-span circular cylindrical vessel that was partially filled, rigid-walled, and horizontally mounted. The numerical simulations illustrated the important effects of the container length and depth of the filled liquid on the calculated sloshing frequencies. Nicolsen et al. [15] developed a new total Lagrangian continuum-based liquid-sloshing model for a vehicle. The results demonstrated that liquid sloshing changed the contact forces between the tires and ground and thereby increased the forces on certain wheels while decreasing the forces on the other wheels. Farid and Gendelman [16] developed an equivalent mechanical model of the nonlinear liquid sloshing occurring in a partially liquid-filled vessel and explored its dynamic regimes. All the analytical predictions were in good agreement with the direct numerical simulations of the initial reduced-order model.The aforementioned research mainly introduced and solved the sloshing caused by a liquid-carrying storage tank on the surface of the earth. For aerospace field applications, the effect of sloshing on a vehicle and the attitude accuracy must be considered. If the amplitude of the liquid sloshing in a spacecraft increases, then its attitude angle will also increase under the sloshing inertia force, which will diverge the spacecraft attitude movement. Liquid sloshing will also cause a change in the liquid centroid, which will produce a change in the attitude angle. Concurrently, the change in the attitude angle will further aggravate the liquid sloshing. This coupled effect will destabilize the spacecraft centroid and affect the attitude control. Because of these challenges, the study of liquid sloshing is gaining attention of aerospace researchers.Considering the large amplitude of the liquid sloshing occurring in a spacecraft storage tank, Hang et al. [17] established a 3D surface-equivalent model of the centroid moving arbitrarily. The simulation results were compared with the numerical results based on the volume-of-fluid method to verify the effectiveness of the model for the analysis of liquid sloshing. Ahmad and Yue [18] developed a canonical liquid-sloshing Hamiltonian model for a container coupled with a spacecraft. The results demonstrated that the Hamiltonian–Casimir technique was effective for studying the spacecraft attitude dynamics, aerodynamics, and other fluid-carrying moving bodies that could be affected by liquid sloshing. Lazzarin et al. [19] studied the impact of propellant sloshing on the pointing stability of the EUCLID satellite. A new method was proposed and applied to a specific fluid system based on computational fluid dynamics instead of real experiments. Zhou and Huang [20] developed a 3D constraint surface mechanical model for the large-amplitude sloshing in propellant tanks. The model validation was further extended to include the in-orbit experimental data of the SloshSat FLEVO satellite, and the numerical results agreed with the experimental data. Chiba and Magata [21] studied the effect of liquid sloshing on the dynamics of flexible space structures. The results indicated that the vibration characteristics of the coupled system were dependent on the static contact angle of the liquid. Although these studies were aimed at the liquid-sloshing characteristics of spacecrafts and their effect on the structure and attitude of the orbiting spacecrafts, they did not provide a method for effectively establishing the relationship between the angular velocity of the satellite attitude and various disturbance factors.Herein, to determine the effect of liquid sloshing on the satellite attitude angular velocity and accuracy of the satellite attitude control during orbit, a method for predicting the liquid-sloshing characteristics of a microsatellite propulsion system is proposed. Three states are considered for the following models for determining the liquid-sloshing characteristics: (1) storage tank, (2) coupling of the storage tank and microsatellite, and (3) coupling of the storage tank and microsatellite with a deployment mechanism. The SPH method is used to study the effect of the sloshing on the storage tank and microsatellite disturbance for each of the three states. The relationship models between the sloshing time, sloshing angular velocity, amount of fluid filled, and satellite attitude angular velocity are established. The models are then used to predict the disturbance angular velocity of the microsatellite, which is compared with the orbital satellite data to illustrate the accuracy of the analysis and model. The results of this study can provide theoretical guidance for the engineering design and orbit control of a microsatellite propulsion system.


## 2. Liquid-Ammonia Propulsion System of a Microsatellite

The liquid-ammonia propulsion system of a microsatellite considered here as reference mainly consists of a storage tank, pipeline, filter, pressure sensor, injection port, isolation valve, battery valve, and thruster, as shown in Figure1. Because of the envelope-size limitation, the storage tank design is approximated as a column structure with an elliptical cross section. The storage tank has a height, width, length, and thickness of 66 mm, 130 mm, 112 mm, and 2.1 mm, respectively. The tank quality is 0.368 kg and the volume is 0.607 L. The ammonia propellant fills the inside of the storage tank under a pressure of 2.4 MPa, and it can be liquefied.Figure 1
Propulsion system composition.As previously mentioned, the following three models are established to analyze the effect of the propellant-filled storage tank on the attitude angular velocity of the storage tank and satellite: (1) storage tank model, (2) storage tank and microsatellite coupling model, and (3) storage tank and microsatellite with a deployment mechanism coupling model. The satellite has a 250 mm3 volume with a mass of 20 kg. The deployment mechanism is a hexagonal structure with a mass of 5 kg. Based on the installation of the microsatellite propulsion system, the storage tank propulsion system is connected to the satellite by a fixture tool, and the center of the storage tank coincides with the geometric center of the satellite. The geometry of the storage tank and a simplified illustration of its installation on the satellite are shown in Figure 2, and the moments of inertia of the three models when the storage tank is filled with a propellant are presented in Table 1.Figure 2
The shape and size of the storage tank and the position of the installation in the satellite: (a) the storage tank model; (b) the storage tank and microsatellite coupling model; (c) the storage tank and microsatellite with a deployment mechanism coupling model.Table 1
The moment of inertia of three cases.
Case
Weight (kg)
Ixx (kg m2)
Iyy (kg m2)
Izz (kg m2)
1
0.632
1.111 × 10−3
0.892 × 10−4
1.456 × 10−3
2
15.234
0.161
0.168
0.155
3
20
0.933
0.936
1.501


## 3. SPH Theory

The SPH method uses a series of randomly distributed disconnected particles to represent the problem domain. The integral method is used to approximate the field function. Furthermore, the particles are used to approximate the kernel equation. The sum of the superposition values of the adjacent particles in a local region is used to replace the integral of the field function and its derivative. The particle-approximation process is conducted at each time step, and therefore, the number of particles used depends on the current locally distributed particles. The particle approximation method is applied to all field functions, and a series of process quantities that are discretized with time are obtained. The process quantities are solved by the explicit integration method to obtain the fastest time integral.


### 3.1. Basic Principle of the SPH

The kernel approximation method is the gradual integration of any function with a smooth function. The integral form of any continuous function in variable fieldΩ can be expressed as
(1)fx=∫Ωfx′δx−x′dx′,where fx is a continuous function of coordinate vector x in variable field Ω, which is the integral volume including vectors x. δx−x′ is the Dirac δ function, which can be expressed as
(2)δx−x′=1,x=x′,0,x≠x′.Next,δx−x′ is replaced with Wx−x′,h, which has features similar to δ. In the SPH method, the angle brackets mark the kernel-approximation function, and the integral can be written in the following form:
(3)fx=∫Ωfx′Wx−x′,hdx′,where Wx−x′,h is the smoothing kernel function and h is the smoothing length, which is the distance between particles x and x′ and governs the size of the support domain of the kernel function.If the density of particlej is ρj, for j=1,2,…,N, the mass is mj, and N is the total number of particles present in the particle support domain, then the particle approximation function of particle i can be obtained by the discretization of (3).
(4)fxi=∑j=1NmjρjfxjWijx−x′,h.By using the differentiability of the kernel function, the particle approximation of the function-space derivative is derived as follows:(5)∇fxi=∑j=1Nmjρjfxj∇Wijx−x′,h.Equations (4) and (5) show that the particle-approximation equation converts the continuous integral expression of the function and derivative into a discrete sum of arbitrarily arranged particles, thereby realizing a mesh-less space.A commonly used kernel function is a piecewise third-order function that was proposed by Monaghan.(6)WR,h=ad23−R2+12R3,0≤R<1,162−R3,1≤R≤2,0,R>2,where R is the relative distance between particles x and x′ and ad is 1/h, 15/7πh2, and 3/2πh3 in the one-dimensional, two-dimensional, and 3D space, respectively.The mass, momentum, and energy conservation equations for continuum fluid dynamics are as follows:(7)dρdt=−ρ∇U,(8)dUdt=−ρ∇σ,(9)dEdt=−Pρ∇U,where ρ is the fluid density, U is the velocity vector, σ is the stress tensor, E is the energy, P is hydrostatic pressure, and t is the time.By the transformation of (7)–(9), the following equations are obtained:
(10)dρidt=−∑j=1Nmjui−uj∇Wij,(11)duidt=−∑j=1Nmjσiρi−σjρj∇Wij,(12)dEidt=∑j=1NmjPi+Pj2ρiρjUi−Uj∇Wij,where Wij=Wxi−xj,h and dρi/dt, dUi/dt, and dEi/dt represent the derivative of the corresponding physical quantity of particle i.


### 3.2. Calculation Process of the SPH Finite-Element Method

The SPH finite-element particle method expresses the spatial derivative of a function as an algebraic equation of particle summation, but a physical quantity is a continuous function of time. First, the continuum is discretized into a series of SPH particles having a mass. The velocity of the particles is obtained using the momentum conservation equation, and the space coordinates of the particles are calculated. Next, the strain rates are obtained using the geometric relations, and then the density of the particles is obtained from the mass conservation equation. Then, the specific internal energy of the particles is obtained using the energy equation. The hydrostatic pressure and stress tensor are calculated by applying the state equation, and the calculated data are stored. A loop calculation is performed until the end of the analysis step. The calculation process of the SPH finite-element method is displayed in Figure3.Figure 3
The calculation process of the SPH finite element method.


## 3.1. Basic Principle of the SPH

The kernel approximation method is the gradual integration of any function with a smooth function. The integral form of any continuous function in variable fieldΩ can be expressed as
(1)fx=∫Ωfx′δx−x′dx′,where fx is a continuous function of coordinate vector x in variable field Ω, which is the integral volume including vectors x. δx−x′ is the Dirac δ function, which can be expressed as
(2)δx−x′=1,x=x′,0,x≠x′.Next,δx−x′ is replaced with Wx−x′,h, which has features similar to δ. In the SPH method, the angle brackets mark the kernel-approximation function, and the integral can be written in the following form:
(3)fx=∫Ωfx′Wx−x′,hdx′,where Wx−x′,h is the smoothing kernel function and h is the smoothing length, which is the distance between particles x and x′ and governs the size of the support domain of the kernel function.If the density of particlej is ρj, for j=1,2,…,N, the mass is mj, and N is the total number of particles present in the particle support domain, then the particle approximation function of particle i can be obtained by the discretization of (3).
(4)fxi=∑j=1NmjρjfxjWijx−x′,h.By using the differentiability of the kernel function, the particle approximation of the function-space derivative is derived as follows:(5)∇fxi=∑j=1Nmjρjfxj∇Wijx−x′,h.Equations (4) and (5) show that the particle-approximation equation converts the continuous integral expression of the function and derivative into a discrete sum of arbitrarily arranged particles, thereby realizing a mesh-less space.A commonly used kernel function is a piecewise third-order function that was proposed by Monaghan.(6)WR,h=ad23−R2+12R3,0≤R<1,162−R3,1≤R≤2,0,R>2,where R is the relative distance between particles x and x′ and ad is 1/h, 15/7πh2, and 3/2πh3 in the one-dimensional, two-dimensional, and 3D space, respectively.The mass, momentum, and energy conservation equations for continuum fluid dynamics are as follows:(7)dρdt=−ρ∇U,(8)dUdt=−ρ∇σ,(9)dEdt=−Pρ∇U,where ρ is the fluid density, U is the velocity vector, σ is the stress tensor, E is the energy, P is hydrostatic pressure, and t is the time.By the transformation of (7)–(9), the following equations are obtained:
(10)dρidt=−∑j=1Nmjui−uj∇Wij,(11)duidt=−∑j=1Nmjσiρi−σjρj∇Wij,(12)dEidt=∑j=1NmjPi+Pj2ρiρjUi−Uj∇Wij,where Wij=Wxi−xj,h and dρi/dt, dUi/dt, and dEi/dt represent the derivative of the corresponding physical quantity of particle i.


## 3.2. Calculation Process of the SPH Finite-Element Method

The SPH finite-element particle method expresses the spatial derivative of a function as an algebraic equation of particle summation, but a physical quantity is a continuous function of time. First, the continuum is discretized into a series of SPH particles having a mass. The velocity of the particles is obtained using the momentum conservation equation, and the space coordinates of the particles are calculated. Next, the strain rates are obtained using the geometric relations, and then the density of the particles is obtained from the mass conservation equation. Then, the specific internal energy of the particles is obtained using the energy equation. The hydrostatic pressure and stress tensor are calculated by applying the state equation, and the calculated data are stored. A loop calculation is performed until the end of the analysis step. The calculation process of the SPH finite-element method is displayed in Figure3.Figure 3
The calculation process of the SPH finite element method.


## 4. Analysis of the Initial Liquid Sloshing in the Storage Tank

The attitude and orbit adjustment of a satellite will cause sloshing of a liquid-filled storage tank, leading to satellite disturbance, which will further affect the attitude and orbit control precision. The disturbance produced by excessive liquid sloshing will increase the difficulty of attitude control. To further clarify the effect of the storage tank sloshing on satellite attitude angular velocity, it is necessary to analyze the liquid sloshing of the storage tank under different working conditions. This analysis will assist in determining the effect of liquid sloshing on satellite attitude angular velocity and improve the accuracy of satellite attitude control during orbit.To analyze the attitude disturbance produced by the three different states of storage tank sloshing, the three models are rotated around the orthogonal-body-fixed reference-frame axes. Furthermore, when the storage tank-filled propellant is in orbit, the propellant distribution in the storage tank is different from that when it is on ground. To ignore the difference and for research convenience, we assume that the distribution of the propellant in the storage tank is same as that on the ground. This simplifies the boundary conditions of the contact between the storage tank and propellant. The inside of the storage tank is assumed to be pure liquid ammonia, and the state of the gas–liquid two-phase mixture is not considered. Additionally, the effects of gravity are not considered. A finite-element analysis model of the storage tank for different amounts of fluid fill is established, as presented in Figure4. The material of the storage tank is TC4, and the material properties are listed in Table 2. The tank is divided into 2384 elements and 3606 nodes by the C3D8R mesh type. For the SPH model, the liquid C3D8R solid element is converted to the PC3D particle element. The grid of the elements is 6 mm, the storage tank motion is coupled to the geometric center of storage tank, and a general contact is defined. Because the inner surface of the storage tank is rough, the friction formulation between the tank and liquid ammonia is of rough type. Double-precision analysis is applied in the ABAQUS finite-element software.Figure 4
The finite element analysis model of the storage tank.Table 2
The material properties for the storage tank.
Material
Density (kg/m3)
Modulus of elasticity (GPa)
Poisson ratio
Yield strength (MPa)
Tensile strength (MPa)
TC4
4500
120
0.3
850
1100The behavior of liquid materials during sloshing is very complex. It is necessary to use the equations of state to describe the pressure, volume, and energy characteristics of the materials. The isotropic pressure can be obtained from the Mie–Gruneisen equation of state,(13)Pρ,e=1−12Γ0ηPH+Γ0ρ0Em,(14)PH=ρ0c2η1−sη2,where H is the Hugoniot curve, Em is the internal energy of the unit mass, η=ρ/ρ0−1 is the nominal volume-compression rate, Γ0 is the Gruneisen constant, s is the limit volume-compression rate, and c is the sound velocity in liquid ammonia. When the nominal strain is small, the volume modulus is K=ρ0c2. c and s can be described by a linear relation between the shock velocity and particle velocity.
(15)Us=c+sUp,where Us-Up, the Hugoniot equation of state, provides a moving-liquid material model in which the volume strength of the material is determined. The storage tank propellant is liquid ammonia, whose performance parameters are listed in Table 3.Table 3
The performance parameters of the liquid ammonia.
TC4
Density (kg/m3)
The speed of sound (m/s)
Viscosity (kg/m s)
Coefficient of state equation
Liquid ammonia
729
1100
2.55 × 10−4
s
Γ0
0
0The rotation or disturbance of a storage tank during orbit is mainly caused by the angular velocity produced by the satellite attitude control and satellite separation from the rocket. Hence, the geometric center of the storage tank is assumed as the coordinate origin, the angular velocity is applied around the XYZ axes, and the angular velocity of the rotation process is as follows:(16)ωi=ωi0t,0≤t≤1,ωi00≤t≤1,where i=x,y,z and ωi0 is the initial angular velocity.To obtain an accurate disturbance angular velocity, the stable angular velocity is represented by the mean value of the disturbance angular velocity, which can be expressed as(17)ωi′=0,0≤t≤ti−1,∑j=1kωj′k,t>ti−1,where ωi′ is the mean value of the disturbance angular velocity, ωj′ is the output value of the disturbance angular velocity, k is the number of output values, and ti is the total time when i=x,y,z.The following two states are analyzed with a sloshing angular velocity of 10°/s and rotation time of 5 s: (1) empty storage tank and (2) liquid storage tank that is 50% full. Once the angular velocity is applied, the rotation constraint is removed. The angular velocity of the storage tank throughout the process is analyzed, and the analysis results are presented in Figures5 and 6 and Table 4.Figure 5
Moving speed of the storage tank and liquid: (a) the empty storage tank after 6 s; (b) the empty storage tank after 7 s; (c) the liquid storage tank after 6 s; (d) the liquid storage tank after 7 s.
(a)
(b)
(c)
(d)Figure 6
The disturbance angular velocity.Table 4
Simulation results of the empty storage tank and liquid storage tank.
Status
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx0′ (°/s)
ωy0′ (°/s)
ωz0′ (°/s)
Empty storage tank
5
10
10
10
0
3.495
3.69
3.428
Liquid storage tank
5
10
10
10
50
4.342
4.096
3.903Figures5 and 6 and Table 4 demonstrate that the liquid in the storage tank affects the disturbance angular velocity, which after the filling of the storage tank is higher than that of the empty storage tank. The fluid flow increases the disturbance angular velocity of the storage tank owing to the angular velocity of the storage tank rotation. Table 4 shows that the disturbance angular velocity is directly determined by the sloshing time, sloshing angular velocity, and amount of the fluid filled. The disturbance angular velocity directly governs the attitude-control precision of the satellite. Hence, further optimization and analysis are required to determine the effect of the satellite rotation angle, causing propellant sloshing, on the satellite attitude.


## 5. Sloshing Analysis

Owing to the limited energy available on a microsatellite, the attitude-control ability is poor. An attitude-control subsystem imposes corresponding requirements on the initial attitude after separation from the rocket to place the satellite in orbit. The detailed requirements are listed in Table5.Table 5
Technical indices of attitude angular velocity for microsatellite.
No.
Technical specification
Design value
1
Rolling angular velocity/(°/s)
Δωγ ≤ 2
2
Pitch velocity/(°/s)
Δωϕ ≤ 3
3
Yaw velocity/(°/s)
Δωφ ≤ 3To determine the law of the liquid sloshing in the storage tank, several simulation analyses of the storage tank under different states are discussed below.


### 5.1. Sloshing Analysis of the Storage Tank

To study the disturbance angular velocity of the liquid-sloshing process of the storage tank at different sloshing times, sloshing angular velocities, and liquid-fill volumes, the sloshing time range, sloshing angular velocity, and amount of the filled liquid are limited to 1–12 s, 2–20°/s, and 10–100%, respectively, based on the actual situation of such a satellite in orbit. Table6 presents the factor parameter levels without considering the mixing of the parameters. A uniform test table containing five factors and ten levels is established and summarized in Table 7. Ten storage tank models are established and input to the ABAQUS software for the simulation analysis, and the final results and computation time are summarized in Table 8. The fitting results of the linear regression analysis are shown in Figure 7, and the values of ωxf1′, ωyf1′, and ωzf1′ from the linear regression equations are presented in Table 9.Table 6
Parameter-level table.
Level
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
1
1
2
2
2
10
2
2
4
4
4
20
3
3
6
6
6
30
4
4
8
8
8
40
5
5
10
10
10
50
6
6
12
12
12
60
7
7
14
14
14
70
8
8
16
16
16
80
9
10
18
18
18
90
10
12
20
20
20
100Table 7
Five factors and ten levels uniform test table.
Level
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
1
5
16
4
16
90
2
6
2
12
14
10
3
10
10
14
12
100
4
3
12
6
4
20
5
4
4
20
6
80
6
1
18
16
10
40
7
2
8
8
20
60
8
8
14
18
18
30
9
9
6
2
8
50
10
7
20
10
2
70Table 8
Simulation results of empty storage tank.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ωx1′ (°/s)
ωy1′ (°/s)
ωz1′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
6.458
1.999
5.174
147
2
6
2
12
14
0.807
4.245
4.707
166
3
10
10
14
12
4.021
5.215
4.328
218
4
3
12
6
4
4.887
2.275
1.289
117
5
4
4
20
6
2.011
6.653
1.563
133
6
1
18
16
10
6.356
5.678
3.101
99
7
2
8
8
20
3.015
2.411
6.391
106
8
8
14
18
18
5.125
6.456
5.717
186
9
9
6
2
8
2.189
0.939
2.561
216
10
7
20
10
2
6.894
3.988
0.715
189Figure 7
Fitting polynomial and fitting results of empty storage tank forωxf1′, ωyf1′, and ωzf1′.Table 9
ωxf1′, ωyf1′, and ωzf1′ of the linear regression equations.
No.
Regression equation
MultipleR
1
ωxf1′=0.385+0.345ωx
0.988
2
ωyf1′=0.375+0.328ωy
0.99
3
ωzf1′=−0.039+0.327ωz
0.994To further analyze the disturbance angular velocity of the liquid sloshing of the storage tank under different sloshing times, sloshing angular velocities, and liquid-fill volumes, similar to the previous analyses, which did not consider the mixing among the various factors, the sloshing time range of the tank, sloshing angular velocity, and amount of the filled liquid in the storage tank are limited to 1–12 s, 2–20°/s, and 10–100%, respectively. Additionally, a uniform test table of five factors and ten levels is generated, and ten storage tank models are established and analyzed by ABAQUS. The final simulation analysis results and computation time are listed in Table10, fitting results are shown in Figure 8, and values of ωxf2′, ωyf2′, and ωzf2′ from the surface regression equations are presented in Table 11.Table 10
Simulation results of liquid storage tank.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx2′ (°/s)
ωy2′ (°/s)
ωz2′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
90
7.519
2.401
7.037
1346
2
6
2
12
14
10
0.938
4.774
6.038
432
3
10
10
14
12
100
4.765
6.058
5.876
2580
4
3
12
6
4
20
5.568
2.587
1.744
323
5
4
4
20
6
80
2.366
7.722
2.521
1172
6
1
18
16
10
40
7.367
6.477
4.299
377
7
2
8
8
20
60
3.582
2.926
8.476
723
8
8
14
18
18
30
5.936
7.291
7.482
838
9
9
6
2
8
50
2.644
1.162
3.549
1330
10
7
20
10
2
70
8.135
4.596
1.238
1570Figure 8
Fitting polynomial and surface fitting results of liquid storage tank for (a)ωxf2′, (b) ωyf2′, and (c) ωzf2′.
(a)
(b)
(c)Table 11
ωxf2′, ωyf2′, and ωzf2′ of the surface regression equations.
No.
Regression equation
MultipleR
1
ωxf2′=0.189+0.394ωx+6.434×10−3ν
0.993
2
ωyf2′=0.264+0.37ωy+4.721×10−3ν
0.992
3
ωzf2′=−0.045+0.415ωz+5.469×10−3ν
0.997


### 5.2. Sloshing Analysis of the Storage Tank in Microsatellite

The analysis results in Sloshing Analysis of Storage Tank indicate that it is necessary to analyze the sloshing by coupling the storage tank and microsatellite to evaluate the effect of the sloshing on the satellite attitude angular velocity. Accordingly, a finite-element analysis model that couples the storage tank and microsatellite with different amounts of the filled liquid is established, as shown in Figure9. The parameters of the microsatellite materials are listed in Table 12. The size of the storage tank, grid divisions, and performance of the liquid material are the same as those in Section 5.1. The satellite is divided into 1408 elements and 1014 nodes using the C3D8R mesh type.Figure 9
Finite element analysis model that couples the storage tank and microsatellite.Table 12
Material characteristics for the microsatellite.
Material
Density (kg/m3)
Modulus of elasticity (GPa)
Poisson ratio
Yield strength (MPa)
Tensile strength (MPa)
2A12
2780
72
0.3
270
440Using the sloshing times, sloshing angular velocity, and amounts of filled liquid in the storage tank that are presented in Section5.1, a uniform test table is established containing four factors and ten levels. The final simulation analysis results and computation time obtained using ABAQUS are presented in Table 13, fitting results are shown in Figure 10, and values of ωxf3′, ωyf3′, and ωzf3′ from the linear regression equations are listed in Table 14.Table 13
Simulation results of empty storage tank and microsatellite.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ωx3′ (°/s)
ωy3′ (°/s)
ωz3′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
7.97
1.963
8.108
238
2
6
2
12
14
0.799
6.163
7.462
248
3
10
10
14
12
5.478
7.036
6.264
347
4
3
12
6
4
6.063
3.176
2.11
181
5
4
4
20
6
1.934
10.699
3.305
196
6
1
18
16
10
8.976
8.481
5.196
131
7
2
8
8
20
3.801
4.206
10.137
156
8
8
14
18
18
7.7
9.287
9.232
303
9
9
6
2
8
2.951
1.002
4.069
317
10
7
20
10
2
10.436
4.868
1.068
259Figure 10
Fitting polynomial and fitting results of empty storage tank and microsatellite forωxf3′, ωyf3′, and ωzf3′.Table 14
ωxf3′, ωyf3′, and ωzf3′ of the linear regression equations.
No.
Regression equation
MultipleR
1
ωxf3′=−0.169+0.525ωx
0.996
2
ωyf3′=−0.14+0.529ωy
0.999
3
ωzf3′=0.137+0.505ωz
0.999The coupled sloshing of the propellant and microsatellite is concurrently analyzed. The final simulation results and computation time obtained from ABAQUS are listed in Table15, fitting results are shown in Figure 11, and values of ωxf4′, ωyf4′, and ωzf4′ of the surface regression equations are provided in Table 16.Table 15
Simulation results of liquid storage tank and microsatellite.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx4′ (°/s)
ωy4′ (°/s)
ωz4′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
90
8.203
2.088
8.273
1402
2
6
2
12
14
10
0.828
6.252
7.548
514
3
10
10
14
12
100
5.644
7.227
6.416
3000
4
3
12
6
4
20
6.216
3.231
2.149
354
5
4
4
20
6
80
2.009
10.897
3.406
1224
6
1
18
16
10
40
9.195
8.616
5.287
360
7
2
8
8
20
60
3.928
4.307
10.298
772
8
8
14
18
18
30
7.868
9.431
9.357
882
9
9
6
2
8
50
3.049
1.073
4.156
1422
10
7
20
10
2
70
10.729
5.004
1.139
1615Figure 11
Fitting polynomial and surface fitting results of liquid storage tank and microsatellite for (a)ωxf4′, (b) ωyf4′, (c) ωzf4′.
(a)
(b)
(c)Table 16
ωxf4′, ωyf4′, and ωzf4′ of the surface regression equations.
No.
Regression equation
MultipleR
1
ωxf4′=−0.215+0.536ωx+1.394×10−3ν
0.997
2
ωyf4′=−0.06+0.536ωy−3.449×10−4ν
0.999
3
ωzf4′=0.196+0.511ωz−1.897×10−4ν
0.999


### 5.3. Analysis of the Storage Tank and Microsatellite with the Deployment Mechanism

Based on the analysis results of Sections5.1 and 5.2, the disturbance angular velocity of the storage tank varies in different states. The satellite has an unfolding mechanism for use during orbit. Therefore, it is necessary to analyze the effect of the sloshing in the storage tank and microsatellite with the deployment mechanism on the satellite attitude angular velocity. A finite-element analysis model with different filled liquid levels in the storage tank and microsatellite with the deployment mechanism is established, as shown in Figure 12. The material parameters of the deployment mechanism are listed in Table 17. The deployment mechanism is divided into 344 elements and 1218 nodes using the C3D8R mesh type. A test table of four factors and ten levels is established to obtain the simulation results and computation time in ABAQUS, as presented in Table 18. The fitting results are displayed in Figure 13, and the values of ωxf5′, ωyf5′, and ωzf5′ from the linear regression equations are listed in Table 19.Figure 12
Finite element model of storage tank and microsatellite with the deployment mechanism.Table 17
Material parameters of the deployment mechanism.
Material
Density (kg/m3)
Modulus of elasticity (GPa)
Poisson ratio
Tensile strength (MPa)
Carbon fiber T700
1800
126
0.31
3500Table 18
Simulation results of empty storage tank and microsatellite with the deployment mechanism.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ωx5′ (°/s)
ωy5′ (°/s)
ωz5′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
2.596
0.677
1.826
1065
2
6
2
12
14
0.368
1.813
1.518
1247
3
10
10
14
12
1.489
2.212
1.231
1904
4
3
12
6
4
1.762
0.962
0.38
781
5
4
4
20
6
0.512
3.081
0.558
969
6
1
18
16
10
2.858
2.541
0.972
502
7
2
8
8
20
1.112
1.232
2.234
650
8
8
14
18
18
2.166
2.868
2.076
1548
9
9
6
2
8
0.809
0.345
0.774
1690
10
7
20
10
2
3.133
1.521
0.271
1429Figure 13
Fitting polynomial and fitting results of empty storage tank and microsatellite with the deployment mechanism forωxf5′, ωyf5′, and ωzf5′.Table 19
ωxf5′, ωyf5′, and ωzf5′ of the linear regression equations.
No.
Regression equation
MultipleR
1
ωxf5′=−0.109+0.163ωx
0.997
2
ωyf5′=0.022+0.155ωy
0.999
3
ωzf5′=−0.095+0.116ωz
0.995The sloshing in the liquid storage tank and microsatellite with the deployment mechanism is analyzed similarly. The results of the final simulation analysis and computation time obtained from ABAQUS are given in Table20, fitting results are shown in Figure 14, and values of ωxf6′, ωyf6′, and ωzf6′ from the surface regression equations are presented in Table 21.Table 20
Simulation results of liquid storage tank and microsatellite with the deployment mechanism.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx6′ (°/s)
ωy6′ (°/s)
ωz6′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
90
2.649
0.716
1.883
6993
2
6
2
12
14
10
0.375
1.839
1.547
2696
3
10
10
14
12
100
1.536
2.273
1.278
15,109
4
3
12
6
4
20
1.793
0.98
0.391
1725
5
4
4
20
6
80
0.541
3.145
0.593
5664
6
1
18
16
10
40
2.906
2.589
1.001
1708
7
2
8
8
20
60
1.143
1.265
2.289
3565
8
8
14
18
18
30
2.203
2.917
2.119
4708
9
9
6
2
8
50
0.835
0.366
0.802
7356
10
7
20
10
2
70
3.191
1.563
0.298
7333Figure 14
Fitting polynomial and surface fitting results of liquid storage tank and microsatellite with the deployment mechanism for (a)ωxf6′, (b) ωyf6′, and (c) ωzf6′.
(a)
(b)
(c)Table 21
ωxf6′, ωyf6′, and ωzf6′ of the surface regression equations.
No.
Regression equation
MultipleR
1
ωxf6′=−0.089+0.165ωx+1.678×10−4ν
0.997
2
ωyf6′=0.165+0.157ωy+4.274×10−4ν
0.999
3
ωzf6′=−0.095+0.118ωz+2.908×10−4ν
0.995


### 5.4. Analysis of the Sloshing Characteristics and Comparison of the Calculation and Simulation

The above analysis demonstrated a linear relationship between the disturbance angular velocity of the empty storage tank and sloshing angular velocity. Furthermore, the disturbance angular velocity generated by the liquid storage tank has a surface relationship with the sloshing angular velocity and amount of the filled liquid. The two sloshing models are independent of time, and the multiple-R value of the model is larger than 0.98, indicating the validity of the model. For the three considered cases, the disturbance angular velocity of the liquid storage tank is higher than that of the empty storage tank, and the disturbance decreases with an increase in the mass of the coupling.The simulation results of the three states of the storage tanks obtained when the amount of the filled liquid propellant in the storage tank and microsatellite during orbit is 100%, and the separation attitude angular velocity is satisfied according to the separation requirements are shown in Figures15–20. Moreover, the results of the simulation analysis and errors are compared in Tables 22–24. The tables show that the maximum and minimum deviation of the three states is 7.6% and is 1.1%, respectively, which indicates that the regression model of the tank with sloshing is effective.Figure 15
Simulation results of empty storage tank forωx1′, ωy1′, and ωz1′ when ωx′=3, ωy′=3, and ωz′=2.Figure 16
Simulation results of liquid storage tank forωx2′, ωy2′, and ωz2′ when ωx′=3, ωy′=3, ωz′=2, and ν=100.Figure 17
Simulation results of empty storage tank and microsatellite forωx3′, ωy3′, and ωz3′ when ωx′=3, ωy′=3, and ωz′=2.Figure 18
Simulation results of liquid storage tank and microsatellite forωx4′, ωy4′, and ωz4′ when ωx′=3, ωy′=3, ωz′=2, and ν=100.Figure 19
Simulation results of empty storage tank and microsatellite with the deployment mechanism forωx5′, ωy5′, and ωz5′ when ωx′=3, ωy′=3, and ωz′=2.Figure 20
Simulation results of liquid storage tank and microsatellite with the deployment mechanism forωxf6′, ωyf6′, and ωzf6′ when ωx′=3, ωy′=3, ωz′=2, and ν=100.Table 22
Relative errors of calculation and simulation results for storage tank.
Calculation results
Simulation results
Relative errors
ωxf1′=1.419
ωx1′=1.386
2.4%
ωyf1′=1.36
ωy1′=1.345
1.1%
ωzf1′=0.615
ωz1′=0.599
2.7%
ωxf2′=2.016
ωx2′=1.911
5.5%
ωyf2′=1.848
ωy2′=1.789
3.3%
ωzf2′=1.333
ωz2′=1.305
2.1%Table 23
Relative errors of calculation and simulation results for storage tank and microsatellite.
Calculation results
Simulation results
Relative errors
ωxf3′=1.407
ωx3′=1.38
1.9%
ωyf3′=1.449
ωy3′=1.418
2.2%
ωzf3′=1.148
ωz3′=1.106
3.8%
ωxf4′=1.535
ωx4′=1.553
1.2%
ωyf4′=1.512
ωy4′=1.633
7.4%
ωzf4′=1.198
ωz4′=1.155
3.7%Table 24
Relative errors of calculation and simulation results for storage tank and microsatellite with the deployment mechanism.
Calculation results (°/s)
Simulation results (°/s)
Relative errors
ωxf5′=0.379
ωx5′=0.374
1.3%
ωyf5′=0.486
ωy5′=0.457
6.3%
ωzf5′=0.138
ωz5′=0.144
4.2%
ωxf6′=0.389
ωx6′=0.383
1.6%
ωyf6′=0.529
ωy6′=0.507
4.3%
ωzf6′=0.17
ωz6′=0.158
7.6%When the sloshing angular velocity inputs for the calculation and simulation analysis areωx=3°/s, ωy=3°/s, and ωz=2°/s to generate positive results, the results of the calculation are ωxf6′=0.389°/s, ωyf6′=0.529°/s, and ωzf6′=0.17°/s, respectively, and the results of the simulation are ωx6′=0.383°/s, ωy6′=0.507°/s, and ωz6′=0.158°/s, respectively.


## 5.1. Sloshing Analysis of the Storage Tank

To study the disturbance angular velocity of the liquid-sloshing process of the storage tank at different sloshing times, sloshing angular velocities, and liquid-fill volumes, the sloshing time range, sloshing angular velocity, and amount of the filled liquid are limited to 1–12 s, 2–20°/s, and 10–100%, respectively, based on the actual situation of such a satellite in orbit. Table6 presents the factor parameter levels without considering the mixing of the parameters. A uniform test table containing five factors and ten levels is established and summarized in Table 7. Ten storage tank models are established and input to the ABAQUS software for the simulation analysis, and the final results and computation time are summarized in Table 8. The fitting results of the linear regression analysis are shown in Figure 7, and the values of ωxf1′, ωyf1′, and ωzf1′ from the linear regression equations are presented in Table 9.Table 6
Parameter-level table.
Level
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
1
1
2
2
2
10
2
2
4
4
4
20
3
3
6
6
6
30
4
4
8
8
8
40
5
5
10
10
10
50
6
6
12
12
12
60
7
7
14
14
14
70
8
8
16
16
16
80
9
10
18
18
18
90
10
12
20
20
20
100Table 7
Five factors and ten levels uniform test table.
Level
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
1
5
16
4
16
90
2
6
2
12
14
10
3
10
10
14
12
100
4
3
12
6
4
20
5
4
4
20
6
80
6
1
18
16
10
40
7
2
8
8
20
60
8
8
14
18
18
30
9
9
6
2
8
50
10
7
20
10
2
70Table 8
Simulation results of empty storage tank.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ωx1′ (°/s)
ωy1′ (°/s)
ωz1′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
6.458
1.999
5.174
147
2
6
2
12
14
0.807
4.245
4.707
166
3
10
10
14
12
4.021
5.215
4.328
218
4
3
12
6
4
4.887
2.275
1.289
117
5
4
4
20
6
2.011
6.653
1.563
133
6
1
18
16
10
6.356
5.678
3.101
99
7
2
8
8
20
3.015
2.411
6.391
106
8
8
14
18
18
5.125
6.456
5.717
186
9
9
6
2
8
2.189
0.939
2.561
216
10
7
20
10
2
6.894
3.988
0.715
189Figure 7
Fitting polynomial and fitting results of empty storage tank forωxf1′, ωyf1′, and ωzf1′.Table 9
ωxf1′, ωyf1′, and ωzf1′ of the linear regression equations.
No.
Regression equation
MultipleR
1
ωxf1′=0.385+0.345ωx
0.988
2
ωyf1′=0.375+0.328ωy
0.99
3
ωzf1′=−0.039+0.327ωz
0.994To further analyze the disturbance angular velocity of the liquid sloshing of the storage tank under different sloshing times, sloshing angular velocities, and liquid-fill volumes, similar to the previous analyses, which did not consider the mixing among the various factors, the sloshing time range of the tank, sloshing angular velocity, and amount of the filled liquid in the storage tank are limited to 1–12 s, 2–20°/s, and 10–100%, respectively. Additionally, a uniform test table of five factors and ten levels is generated, and ten storage tank models are established and analyzed by ABAQUS. The final simulation analysis results and computation time are listed in Table10, fitting results are shown in Figure 8, and values of ωxf2′, ωyf2′, and ωzf2′ from the surface regression equations are presented in Table 11.Table 10
Simulation results of liquid storage tank.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx2′ (°/s)
ωy2′ (°/s)
ωz2′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
90
7.519
2.401
7.037
1346
2
6
2
12
14
10
0.938
4.774
6.038
432
3
10
10
14
12
100
4.765
6.058
5.876
2580
4
3
12
6
4
20
5.568
2.587
1.744
323
5
4
4
20
6
80
2.366
7.722
2.521
1172
6
1
18
16
10
40
7.367
6.477
4.299
377
7
2
8
8
20
60
3.582
2.926
8.476
723
8
8
14
18
18
30
5.936
7.291
7.482
838
9
9
6
2
8
50
2.644
1.162
3.549
1330
10
7
20
10
2
70
8.135
4.596
1.238
1570Figure 8
Fitting polynomial and surface fitting results of liquid storage tank for (a)ωxf2′, (b) ωyf2′, and (c) ωzf2′.
(a)
(b)
(c)Table 11
ωxf2′, ωyf2′, and ωzf2′ of the surface regression equations.
No.
Regression equation
MultipleR
1
ωxf2′=0.189+0.394ωx+6.434×10−3ν
0.993
2
ωyf2′=0.264+0.37ωy+4.721×10−3ν
0.992
3
ωzf2′=−0.045+0.415ωz+5.469×10−3ν
0.997


## 5.2. Sloshing Analysis of the Storage Tank in Microsatellite

The analysis results in Sloshing Analysis of Storage Tank indicate that it is necessary to analyze the sloshing by coupling the storage tank and microsatellite to evaluate the effect of the sloshing on the satellite attitude angular velocity. Accordingly, a finite-element analysis model that couples the storage tank and microsatellite with different amounts of the filled liquid is established, as shown in Figure9. The parameters of the microsatellite materials are listed in Table 12. The size of the storage tank, grid divisions, and performance of the liquid material are the same as those in Section 5.1. The satellite is divided into 1408 elements and 1014 nodes using the C3D8R mesh type.Figure 9
Finite element analysis model that couples the storage tank and microsatellite.Table 12
Material characteristics for the microsatellite.
Material
Density (kg/m3)
Modulus of elasticity (GPa)
Poisson ratio
Yield strength (MPa)
Tensile strength (MPa)
2A12
2780
72
0.3
270
440Using the sloshing times, sloshing angular velocity, and amounts of filled liquid in the storage tank that are presented in Section5.1, a uniform test table is established containing four factors and ten levels. The final simulation analysis results and computation time obtained using ABAQUS are presented in Table 13, fitting results are shown in Figure 10, and values of ωxf3′, ωyf3′, and ωzf3′ from the linear regression equations are listed in Table 14.Table 13
Simulation results of empty storage tank and microsatellite.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ωx3′ (°/s)
ωy3′ (°/s)
ωz3′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
7.97
1.963
8.108
238
2
6
2
12
14
0.799
6.163
7.462
248
3
10
10
14
12
5.478
7.036
6.264
347
4
3
12
6
4
6.063
3.176
2.11
181
5
4
4
20
6
1.934
10.699
3.305
196
6
1
18
16
10
8.976
8.481
5.196
131
7
2
8
8
20
3.801
4.206
10.137
156
8
8
14
18
18
7.7
9.287
9.232
303
9
9
6
2
8
2.951
1.002
4.069
317
10
7
20
10
2
10.436
4.868
1.068
259Figure 10
Fitting polynomial and fitting results of empty storage tank and microsatellite forωxf3′, ωyf3′, and ωzf3′.Table 14
ωxf3′, ωyf3′, and ωzf3′ of the linear regression equations.
No.
Regression equation
MultipleR
1
ωxf3′=−0.169+0.525ωx
0.996
2
ωyf3′=−0.14+0.529ωy
0.999
3
ωzf3′=0.137+0.505ωz
0.999The coupled sloshing of the propellant and microsatellite is concurrently analyzed. The final simulation results and computation time obtained from ABAQUS are listed in Table15, fitting results are shown in Figure 11, and values of ωxf4′, ωyf4′, and ωzf4′ of the surface regression equations are provided in Table 16.Table 15
Simulation results of liquid storage tank and microsatellite.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx4′ (°/s)
ωy4′ (°/s)
ωz4′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
90
8.203
2.088
8.273
1402
2
6
2
12
14
10
0.828
6.252
7.548
514
3
10
10
14
12
100
5.644
7.227
6.416
3000
4
3
12
6
4
20
6.216
3.231
2.149
354
5
4
4
20
6
80
2.009
10.897
3.406
1224
6
1
18
16
10
40
9.195
8.616
5.287
360
7
2
8
8
20
60
3.928
4.307
10.298
772
8
8
14
18
18
30
7.868
9.431
9.357
882
9
9
6
2
8
50
3.049
1.073
4.156
1422
10
7
20
10
2
70
10.729
5.004
1.139
1615Figure 11
Fitting polynomial and surface fitting results of liquid storage tank and microsatellite for (a)ωxf4′, (b) ωyf4′, (c) ωzf4′.
(a)
(b)
(c)Table 16
ωxf4′, ωyf4′, and ωzf4′ of the surface regression equations.
No.
Regression equation
MultipleR
1
ωxf4′=−0.215+0.536ωx+1.394×10−3ν
0.997
2
ωyf4′=−0.06+0.536ωy−3.449×10−4ν
0.999
3
ωzf4′=0.196+0.511ωz−1.897×10−4ν
0.999


## 5.3. Analysis of the Storage Tank and Microsatellite with the Deployment Mechanism

Based on the analysis results of Sections5.1 and 5.2, the disturbance angular velocity of the storage tank varies in different states. The satellite has an unfolding mechanism for use during orbit. Therefore, it is necessary to analyze the effect of the sloshing in the storage tank and microsatellite with the deployment mechanism on the satellite attitude angular velocity. A finite-element analysis model with different filled liquid levels in the storage tank and microsatellite with the deployment mechanism is established, as shown in Figure 12. The material parameters of the deployment mechanism are listed in Table 17. The deployment mechanism is divided into 344 elements and 1218 nodes using the C3D8R mesh type. A test table of four factors and ten levels is established to obtain the simulation results and computation time in ABAQUS, as presented in Table 18. The fitting results are displayed in Figure 13, and the values of ωxf5′, ωyf5′, and ωzf5′ from the linear regression equations are listed in Table 19.Figure 12
Finite element model of storage tank and microsatellite with the deployment mechanism.Table 17
Material parameters of the deployment mechanism.
Material
Density (kg/m3)
Modulus of elasticity (GPa)
Poisson ratio
Tensile strength (MPa)
Carbon fiber T700
1800
126
0.31
3500Table 18
Simulation results of empty storage tank and microsatellite with the deployment mechanism.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ωx5′ (°/s)
ωy5′ (°/s)
ωz5′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
2.596
0.677
1.826
1065
2
6
2
12
14
0.368
1.813
1.518
1247
3
10
10
14
12
1.489
2.212
1.231
1904
4
3
12
6
4
1.762
0.962
0.38
781
5
4
4
20
6
0.512
3.081
0.558
969
6
1
18
16
10
2.858
2.541
0.972
502
7
2
8
8
20
1.112
1.232
2.234
650
8
8
14
18
18
2.166
2.868
2.076
1548
9
9
6
2
8
0.809
0.345
0.774
1690
10
7
20
10
2
3.133
1.521
0.271
1429Figure 13
Fitting polynomial and fitting results of empty storage tank and microsatellite with the deployment mechanism forωxf5′, ωyf5′, and ωzf5′.Table 19
ωxf5′, ωyf5′, and ωzf5′ of the linear regression equations.
No.
Regression equation
MultipleR
1
ωxf5′=−0.109+0.163ωx
0.997
2
ωyf5′=0.022+0.155ωy
0.999
3
ωzf5′=−0.095+0.116ωz
0.995The sloshing in the liquid storage tank and microsatellite with the deployment mechanism is analyzed similarly. The results of the final simulation analysis and computation time obtained from ABAQUS are given in Table20, fitting results are shown in Figure 14, and values of ωxf6′, ωyf6′, and ωzf6′ from the surface regression equations are presented in Table 21.Table 20
Simulation results of liquid storage tank and microsatellite with the deployment mechanism.
No.
t (s)
ωx (°/s)
ωy (°/s)
ωz (°/s)
ν (%)
ωx6′ (°/s)
ωy6′ (°/s)
ωz6′ (°/s)
Computation time in ABAQUS (s)
1
5
16
4
16
90
2.649
0.716
1.883
6993
2
6
2
12
14
10
0.375
1.839
1.547
2696
3
10
10
14
12
100
1.536
2.273
1.278
15,109
4
3
12
6
4
20
1.793
0.98
0.391
1725
5
4
4
20
6
80
0.541
3.145
0.593
5664
6
1
18
16
10
40
2.906
2.589
1.001
1708
7
2
8
8
20
60
1.143
1.265
2.289
3565
8
8
14
18
18
30
2.203
2.917
2.119
4708
9
9
6
2
8
50
0.835
0.366
0.802
7356
10
7
20
10
2
70
3.191
1.563
0.298
7333Figure 14
Fitting polynomial and surface fitting results of liquid storage tank and microsatellite with the deployment mechanism for (a)ωxf6′, (b) ωyf6′, and (c) ωzf6′.
(a)
(b)
(c)Table 21
ωxf6′, ωyf6′, and ωzf6′ of the surface regression equations.
No.
Regression equation
MultipleR
1
ωxf6′=−0.089+0.165ωx+1.678×10−4ν
0.997
2
ωyf6′=0.165+0.157ωy+4.274×10−4ν
0.999
3
ωzf6′=−0.095+0.118ωz+2.908×10−4ν
0.995


## 5.4. Analysis of the Sloshing Characteristics and Comparison of the Calculation and Simulation

The above analysis demonstrated a linear relationship between the disturbance angular velocity of the empty storage tank and sloshing angular velocity. Furthermore, the disturbance angular velocity generated by the liquid storage tank has a surface relationship with the sloshing angular velocity and amount of the filled liquid. The two sloshing models are independent of time, and the multiple-R value of the model is larger than 0.98, indicating the validity of the model. For the three considered cases, the disturbance angular velocity of the liquid storage tank is higher than that of the empty storage tank, and the disturbance decreases with an increase in the mass of the coupling.The simulation results of the three states of the storage tanks obtained when the amount of the filled liquid propellant in the storage tank and microsatellite during orbit is 100%, and the separation attitude angular velocity is satisfied according to the separation requirements are shown in Figures15–20. Moreover, the results of the simulation analysis and errors are compared in Tables 22–24. The tables show that the maximum and minimum deviation of the three states is 7.6% and is 1.1%, respectively, which indicates that the regression model of the tank with sloshing is effective.Figure 15
Simulation results of empty storage tank forωx1′, ωy1′, and ωz1′ when ωx′=3, ωy′=3, and ωz′=2.Figure 16
Simulation results of liquid storage tank forωx2′, ωy2′, and ωz2′ when ωx′=3, ωy′=3, ωz′=2, and ν=100.Figure 17
Simulation results of empty storage tank and microsatellite forωx3′, ωy3′, and ωz3′ when ωx′=3, ωy′=3, and ωz′=2.Figure 18
Simulation results of liquid storage tank and microsatellite forωx4′, ωy4′, and ωz4′ when ωx′=3, ωy′=3, ωz′=2, and ν=100.Figure 19
Simulation results of empty storage tank and microsatellite with the deployment mechanism forωx5′, ωy5′, and ωz5′ when ωx′=3, ωy′=3, and ωz′=2.Figure 20
Simulation results of liquid storage tank and microsatellite with the deployment mechanism forωxf6′, ωyf6′, and ωzf6′ when ωx′=3, ωy′=3, ωz′=2, and ν=100.Table 22
Relative errors of calculation and simulation results for storage tank.
Calculation results
Simulation results
Relative errors
ωxf1′=1.419
ωx1′=1.386
2.4%
ωyf1′=1.36
ωy1′=1.345
1.1%
ωzf1′=0.615
ωz1′=0.599
2.7%
ωxf2′=2.016
ωx2′=1.911
5.5%
ωyf2′=1.848
ωy2′=1.789
3.3%
ωzf2′=1.333
ωz2′=1.305
2.1%Table 23
Relative errors of calculation and simulation results for storage tank and microsatellite.
Calculation results
Simulation results
Relative errors
ωxf3′=1.407
ωx3′=1.38
1.9%
ωyf3′=1.449
ωy3′=1.418
2.2%
ωzf3′=1.148
ωz3′=1.106
3.8%
ωxf4′=1.535
ωx4′=1.553
1.2%
ωyf4′=1.512
ωy4′=1.633
7.4%
ωzf4′=1.198
ωz4′=1.155
3.7%Table 24
Relative errors of calculation and simulation results for storage tank and microsatellite with the deployment mechanism.
Calculation results (°/s)
Simulation results (°/s)
Relative errors
ωxf5′=0.379
ωx5′=0.374
1.3%
ωyf5′=0.486
ωy5′=0.457
6.3%
ωzf5′=0.138
ωz5′=0.144
4.2%
ωxf6′=0.389
ωx6′=0.383
1.6%
ωyf6′=0.529
ωy6′=0.507
4.3%
ωzf6′=0.17
ωz6′=0.158
7.6%When the sloshing angular velocity inputs for the calculation and simulation analysis areωx=3°/s, ωy=3°/s, and ωz=2°/s to generate positive results, the results of the calculation are ωxf6′=0.389°/s, ωyf6′=0.529°/s, and ωzf6′=0.17°/s, respectively, and the results of the simulation are ωx6′=0.383°/s, ωy6′=0.507°/s, and ωz6′=0.158°/s, respectively.


## 6. Results and Discussion

The Zhejiang University microsatellites are double satellites; microsatellite A is filled with a propellant, and microsatellite B is not filled with a propellant. Following the separation of the satellite from the rocket during orbit, the angular velocities of the three axes are satisfied with the separation precision listed in Table5. Then, the deployment mechanism on the satellites is executed, and the orbit data of the entire process is relayed to the ground. The changing process of the attitude angular velocity is displayed in Figure 21. Figure 21 illustrates that the attitude angular velocity of the initial satellite separation is high and stable following the attitude adjustment. The angular velocity ranges for the sloshing in microsatellite A are (−0.256, 0.278), (−0.319, 0.349), and (−0.157, 0.164) and in microsatellite-B are (−0.215, 0.238), (−0.248, 0.269), and (−0.138, 0.148) for each axis.Figure 21
Angular velocity of on orbit attitude of Zhejiang University microsatellites: (a) A satellite and (b) B satellite.
(a)
(b)The deviation between the calculation and experimental results is expressed by the following equation:(18)δi=ωif′−ωiωi×100%,where δi is the deviation, i=x,y,z, ωif′ denotes the calculation results, and ωi denotes the experimental results.The calculation and experimental results of microsatellites A and B are summarized in Table25, where test results A and B are calculated according to (17). Table 25 demonstrates that that the maximum disturbance deviation is in the y direction and is 43.36% whereas the minimum deviation is in the z direction and is 14.86%.Table 25
The comparison of calculation and test results.
Technical specification
Calculation results-A
Calculation results-B
Test results-A
Test results-B
Relative errors-A
Relative errors-B
Rolling angular velocity (°/s)
0.138
0.17
0.164
0.148
15.85%
14.86%
Pitch velocity (°/s)
0.379
0.389
0.278
0.338
36.33%
15.09%
Yaw velocity (°/s)
0.486
0.529
0.349
0.369
39.25%
43.36%Because the disturbance angular velocity is low, there is an error between the orbital data and calculation results. This deviation is mainly caused by the simplification of the model, precision of the model, precision of the attitude sensor, and unfolding of the deployment mechanism. However, the calculation results can reflect the trends of the satellite disturbance state and attitude angular velocity of the satellite during orbit. The effect of the initial angular velocity and liquid sloshing on the satellite attitude angular velocity is consistent with the calculation results. This effect is small and of the order of 10−2°/s when the angular velocities of the satellite are ωx′=3, ωy′=3, and ωz′=2.The satellite propellant weight accounts for 0.03% of the entire satellite weight, and the magnitude of the attitude disturbance angular velocity is 10−2°/s. This indicates that the disturbance of the attitude angular velocity of a microsatellite is directly determined by the liquid propellant volume. Thus, the results can be extended to larger satellites, and a quantitative evaluation of the angular velocity of the disturbance derived from sloshing can be obtained. For microsatellite formations and satellites that require high-precision attitude positioning, such as those acquiring photographs and videos, the effect of the disturbance angular velocity of the propellant cannot be ignored. Therefore, it is necessary to improve the accuracy of the liquid-sloshing model and generate a mixture model that considers the gas–liquid two-phase mixture of the propellant. Furthermore, this will help predict the disturbance angular velocity and provide the high-precision and real-time adjustment of the momentum wheel for restraining the satellite attitude disturbance. The method of increasing the structure of the storage tank to prevent liquid sloshing without affecting the tank performance of the propulsion system can also restrain the attitude disturbance. This complex network structure, generated by 3D printing, has been mentioned in the literature [22].


## 7. Conclusions

In this study, a method for establishing the prediction model of the liquid-sloshing characteristics of a microsatellite propulsion system is proposed. The ideal states of the liquid-sloshing characteristics are considered for three cases. The relationship models between the sloshing time, sloshing angular velocity, amount of fluid filled, and satellite attitude angular velocity are established. The analysis results show that the disturbance angular velocity of the empty storage tank and sloshing angular velocity have a linear relationship. Furthermore, the disturbance angular velocity of the liquid storage tank demonstrates a surface relationship with the sloshing angular velocity and amount of fluid filled. Additionally, the disturbance angular velocity for the liquid state of the storage tank is higher than that for the empty state, and the degree of disturbance decreases with increase in the mass of the coupling. The disturbance of the storage tank and satellite with the deployment mechanism is 10−2°/s when the angular velocities of the satellite are ωx′=3, ωy′=3, and ωz′=2. The maximum deviation between the calculation and simulation results of the three models is 7.6%, and the minimum deviation is 1.1%. The model is used to predict the disturbance angular velocity of the microsatellite. When compared with the orbital satellite data, the maximum disturbance angular velocity deviation of 43.36% occurs in the y direction, whereas the minimum deviation of 14.86% occurs in the z direction, which demonstrates the accuracy of the analysis and model.



---

*Source: 2702057-2019-02-03.xml*