Research Article  Open Access
Method for Predicting the LiquidSloshing Characteristics of a Microsatellite Propulsion System
Abstract
In this study, a method for establishing a prediction model for the liquidsloshing characteristics of microsatellite propulsion systems is proposed. The ideal states of the liquidsloshing 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 smoothedparticle 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 liquidfilled storage tank exhibits a surface relationship with the sloshing angular velocity and amount of fluid filled in a threedimensional coordinate system. Additionally, the disturbance angular velocity in the liquidfilled 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 , , and . 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 and 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.
1. Introduction
The continuous improvement of modern technology has led to several important developments such as microelectronics, microcomputers, and new material development. Although largesatellite 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 inorbit propulsion technology of microsatellites has heightened along with that for miniaturized propulsion systems. The propulsion systems of microsatellites have the characteristics of highvolume 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 smoothedparticle hydrodynamics (SPH) method is a Lagrangetype meshless 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 liquidsloshing characteristics of an elastic rectangular tank at different horizontal excitation frequencies and obtained the freesurface 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. CardosoRibeiro et al. [13] proposed portHamiltonian systems (PHs) for formulating the liquidsloshing 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 threedimensional (3D) sloshing in a finitespan circular cylindrical vessel that was partially filled, rigidwalled, 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 continuumbased liquidsloshing 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 liquidfilled vessel and explored its dynamic regimes. All the analytical predictions were in good agreement with the direct numerical simulations of the initial reducedorder model.
The aforementioned research mainly introduced and solved the sloshing caused by a liquidcarrying 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 surfaceequivalent model of the centroid moving arbitrarily. The simulation results were compared with the numerical results based on the volumeoffluid method to verify the effectiveness of the model for the analysis of liquid sloshing. Ahmad and Yue [18] developed a canonical liquidsloshing 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 fluidcarrying 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 largeamplitude sloshing in propellant tanks. The model validation was further extended to include the inorbit 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 liquidsloshing 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 liquidsloshing characteristics of a microsatellite propulsion system is proposed. Three states are considered for the following models for determining the liquidsloshing 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. LiquidAmmonia Propulsion System of a Microsatellite
The liquidammonia 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 Figure 1. Because of the envelopesize 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.
As previously mentioned, the following three models are established to analyze the effect of the propellantfilled 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 mm^{3} 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.

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 particleapproximation 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 where is a continuous function of coordinate vector in variable field , which is the integral volume including vectors . is the Dirac function, which can be expressed as
Next, is replaced with , which has features similar to . In the SPH method, the angle brackets mark the kernelapproximation function, and the integral can be written in the following form: where is the smoothing kernel function and is the smoothing length, which is the distance between particles and and governs the size of the support domain of the kernel function.
If the density of particle is , for , the mass is , and is the total number of particles present in the particle support domain, then the particle approximation function of particle can be obtained by the discretization of (3).
By using the differentiability of the kernel function, the particle approximation of the functionspace derivative is derived as follows:
Equations (4) and (5) show that the particleapproximation equation converts the continuous integral expression of the function and derivative into a discrete sum of arbitrarily arranged particles, thereby realizing a meshless space.
A commonly used kernel function is a piecewise thirdorder function that was proposed by Monaghan. where is the relative distance between particles and and is , , and in the onedimensional, twodimensional, and 3D space, respectively.
The mass, momentum, and energy conservation equations for continuum fluid dynamics are as follows: where is the fluid density, is the velocity vector, is the stress tensor, is the energy, is hydrostatic pressure, and is the time.
By the transformation of (7)–(9), the following equations are obtained: where and , , and represent the derivative of the corresponding physical quantity of particle .
3.2. Calculation Process of the SPH FiniteElement Method
The SPH finiteelement 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 finiteelement method is displayed in Figure 3.
4. Analysis of the Initial Liquid Sloshing in the Storage Tank
The attitude and orbit adjustment of a satellite will cause sloshing of a liquidfilled 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 orthogonalbodyfixed referenceframe axes. Furthermore, when the storage tankfilled 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 twophase mixture is not considered. Additionally, the effects of gravity are not considered. A finiteelement analysis model of the storage tank for different amounts of fluid fill is established, as presented in Figure 4. 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. Doubleprecision analysis is applied in the ABAQUS finiteelement software.

The 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, where is the Hugoniot curve, is the internal energy of the unit mass, is the nominal volumecompression rate, is the Gruneisen constant, is the limit volumecompression rate, and is the sound velocity in liquid ammonia. When the nominal strain is small, the volume modulus is . and can be described by a linear relation between the shock velocity and particle velocity. where , the Hugoniot equation of state, provides a movingliquid 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.

The 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: where and 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 where is the mean value of the disturbance angular velocity, is the output value of the disturbance angular velocity, is the number of output values, and is the total time when .
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 Figures 5 and 6 and Table 4.
(a)
(b)
(c)
(d)

Figures 5 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 attitudecontrol 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 attitudecontrol ability is poor. An attitudecontrol 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 Table 5.

To 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 liquidsloshing process of the storage tank at different sloshing times, sloshing angular velocities, and liquidfill 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. Table 6 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 , , and from the linear regression equations are presented in Table 9.




To further analyze the disturbance angular velocity of the liquid sloshing of the storage tank under different sloshing times, sloshing angular velocities, and liquidfill 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 Table 10, fitting results are shown in Figure 8, and values of , , and from the surface regression equations are presented in Table 11.

(a)
(b)
(c)

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 finiteelement analysis model that couples the storage tank and microsatellite with different amounts of the filled liquid is established, as shown in Figure 9. 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.

Using the sloshing times, sloshing angular velocity, and amounts of filled liquid in the storage tank that are presented in Section 5.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 , , and from the linear regression equations are listed in Table 14.


The coupled sloshing of the propellant and microsatellite is concurrently analyzed. The final simulation results and computation time obtained from ABAQUS are listed in Table 15, fitting results are shown in Figure 11, and values of , , and of the surface regression equations are provided in Table 16.

(a)
(b)
(c)

5.3. Analysis of the Storage Tank and Microsatellite with the Deployment Mechanism
Based on the analysis results of Sections 5.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 finiteelement 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 , , and from the linear regression equations are listed in Table 19.



The 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 Table 20, fitting results are shown in Figure 14, and values of , , and from the surface regression equations are presented in Table 21.

(a)
(b)
(c)

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 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 Figures 15–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.



When the sloshing angular velocity inputs for the calculation and simulation analysis are °/s, °/s, and °/s to generate positive results, the results of the calculation are °/s, °/s, and °/s, respectively, and the results of the simulation are °/s, °/s, and °/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 Table 5. 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 microsatelliteB are (−0.215, 0.238), (−0.248, 0.269), and (−0.138, 0.148) for each axis.
(a)
(b)
The deviation between the calculation and experimental results is expressed by the following equation: where is the deviation, , denotes the calculation results, and denotes the experimental results.
The calculation and experimental results of microsatellites A and B are summarized in Table 25, where test results A and B are calculated according to (17). Table 25 demonstrates that that the maximum disturbance deviation is in the direction and is 43.36% whereas the minimum deviation is in the direction and is 14.86%.

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 , , and .
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 highprecision 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 liquidsloshing model and generate a mixture model that considers the gas–liquid twophase mixture of the propellant. Furthermore, this will help predict the disturbance angular velocity and provide the highprecision and realtime 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 liquidsloshing characteristics of a microsatellite propulsion system is proposed. The ideal states of the liquidsloshing 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 , , and . 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 direction, whereas the minimum deviation of 14.86% occurs in the direction, which demonstrates the accuracy of the analysis and model.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was funded by the China National Funds for Distinguished Young Scientists (grant number 61525403) and Fundamental Research Funds for the Central Universities (grant number 2017QNA4035).
References
 X. L. Ling, “Micropropulsion system for modern small satellites,” Spacecraft Engineering, vol. 19, no. 6, pp. 13–20, 2010. View at: Google Scholar
 M. Ihmsen, J. Cornelis, B. Solenthaler, C. Horvath, and M. Teschner, “Implicit incompressible SPH,” IEEE Transactions on Visualization and Computer Graphics, vol. 20, no. 3, pp. 426–435, 2014. View at: Publisher Site  Google Scholar
 H. Zhang, Research on Fluid Movement and SolidFluid Coupling Simulation Based on SPH, Master’s Thesis, University Of Science and Technology of China, Hefei, China, 2015.
 G. X. Ni, R. L. Wang, Z. Lin, and Q. W. Ge, “Equidistribution of particles in arbitrary domain,” Chinese Journal of Computational Mechanics, vol. 24, no. 4, pp. 408–413, 2007. View at: Google Scholar
 M. B. Liu and J. Z. Chang, “Particle distribution and numerical stability in smoothed particle hydrodynamics method,” Acta Physica Sinica, vol. 59, no. 6, pp. 3654–3662, 2010. View at: Google Scholar
 R. A. Gingold and J. J. Monaghan, “Smoothed particle hydrodynamics: theory and application to nonspherical stars,” Monthly Notices of the Royal Astronomical Society, vol. 181, no. 3, pp. 375–389, 1977. View at: Publisher Site  Google Scholar
 L. B. Lucy, “A numerical approach to the testing of the fission hypothesis,” The Astronomical Journal, vol. 82, pp. 1013–1024, 1977. View at: Publisher Site  Google Scholar
 J. J. Monaghan, “Simulating free surface flows with SPH,” Journal of Computational Physics, vol. 110, no. 2, pp. 399–406, 1994. View at: Publisher Site  Google Scholar
 M. R. Jiang, B. Ren, H. J. Wen, and Y. X. Wang, “Experimental study on sloshing in elastic tanks,” Ocean Engineering, vol. 31, no. 5, pp. 1–10, 2013. View at: Google Scholar
 M. Liu, R. Zhou, and J. Shao, “Numerical simulation of liquid sloshing in aprismatic tank with SPH method,” Journal of Hohai University, vol. 42, no. 3, pp. 257–261, 2014. View at: Google Scholar
 Y. T. Hang, F. X. Du, H. Liu, J. Wang, Z. C. Ou, and X. H. Yao, “Analysis on sloshing characteristics of liquid in elastomer tank based on SPH method,” Computer Aided Engineering, vol. 24, no. 2, pp. 36–41, 2015. View at: Google Scholar
 T. Hao, X. Jian, and Z. Jianhui, “Study of liquid sloshing in missile fuel tank based on SPH method,” Wireless Internet Technology, vol. 3, pp. 76–79, 2015. View at: Google Scholar
 F. L. CardosoRibeiro, D. Matignon, and V. PommierBudinger, “A portHamiltonian model of liquid sloshing in moving containers and application to a fluidstructure system,” Journal of Fluids and Structures, vol. 69, pp. 402–427, 2017. View at: Publisher Site  Google Scholar
 S. M. Hasheminejad and H. Soleimani, “An analytical solution for free liquid sloshing in a finitelength horizontal cylindrical container filled to an arbitrary depth,” Applied Mathematical Modelling, vol. 48, pp. 338–352, 2017. View at: Publisher Site  Google Scholar
 B. Nicolsen, L. Wang, and A. Shabana, “Nonlinear finite element analysis of liquid sloshing in complex vehicle motion scenarios,” Journal of Sound and Vibration, vol. 405, pp. 208–233, 2017. View at: Publisher Site  Google Scholar
 M. Farid and O. V. Gendelman, “Response regimes in equivalent mechanical model of strongly nonlinear liquid sloshing,” International Journal of NonLinear Mechanics, vol. 94, no. 4, pp. 146–159, 2017. View at: Publisher Site  Google Scholar
 H. Hang, L. Yang, H. Zhang, and G. J. Qu, “Three dimensional centroid surface equivalent model for large liquid sloshing in spacecraft tank,” Journal Astronautics, vol. 31, no. 1, pp. 55–59, 2010. View at: Google Scholar
 S. Ahmad and B. Yue, “Motion stability dynamics for spacecraft coupled with partially filled liquid container,” Theoretical and Applied Mechanics Letters, vol. 2, no. 1, article 013002, 2012. View at: Publisher Site  Google Scholar
 M. Lazzarin, M. Biolo, A. Bettella, M. Manente, R. da Forno, and D. Pavarin, “EUCLID satellite: development of a lumped parameter model for sloshing description,” Aerospace Science and Technology, vol. 39, pp. 491–500, 2014. View at: Publisher Site  Google Scholar
 Z. Zhou and H. Huang, “Constraint surface model for large amplitude sloshing of the spacecraft with multiple tanks,” Acta Astronautica, vol. 111, pp. 222–229, 2015. View at: Publisher Site  Google Scholar
 M. Chiba and H. Magata, “Influence of liquid sloshing on dynamics of flexible space structures,” Journal of Sound and Vibration, vol. 401, pp. 1–22, 2017. View at: Publisher Site  Google Scholar
 L. Teng and Z. Jin, “Structural performance optimization and verification of an improved thinwalled storage tank for a Picosatellite,” Applied Sciences, vol. 7, no. 11, p. 1168, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 L. Teng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.