RE.PUBLIC@POLIMI

The sloshing of liquids in microgravity is a relevant problem of applied mechanics with important implications for spacecraft design. A magnetic settling force may be used to avoid the highly non-linear dynamics that characterize these systems. However, this approach is still largely unexplored. This paper presents a quasi-analytical low-gravity sloshing model for magnetic liquids under the action of external inhomogeneous magnetic ﬁelds. The problems of free and forced oscillations are solved for axisymmetric geometries and loads by assuming a linear behavior. The model may be of particular interest for the development of magnetic sloshing damping devices in space, whose behavior can be easily predicted and quantiﬁed with standard mechanical analogies.

The sloshing of liquids in microgravity is a relevant problem of applied mechanics with important implications for spacecraft design.A magnetic settling force may be used to avoid the highly non-linear dynamics that characterize these systems.However, this approach is still largely unexplored.This paper presents a quasi-analytical low-gravity sloshing model for magnetic liquids under the action of external inhomogeneous magnetic fields.The problems of free and forced oscillations are solved for axisymmetric geometries and loads by assuming a linear behavior.The model may be of particular interest for the development of magnetic sloshing damping devices in space, whose behavior can be easily predicted and quantified with standard mechanical analogies.Primitive functions of ζ i {r, θ, z} Cylindrical coordinates of the system {u r , u θ , u z } {u r , u θ , u z } Cylindrical reference system centered at the vertex of the meniscus

Introduction
The term sloshing refers to the forced movement of liquids in partially filled tanks [1].Propellant sloshing has been a major concern for space engineers since the beginning of the space era.During launch, it can result in the partial or total loss of control of the spacecraft [2].In a low-gravity environment, the liquid tends to adopt a random position inside the tank and mixes with pressurizing gas bubbles.This results in a complicated propellant management system design, often increasing the inert mass of the vehicle [1].
Low-gravity sloshing is characterized by the dominant role of surface tension, that produces a curved equilibrium free surface (or meniscus) and a complex interaction with the walls of the vessel that contains the liquid.The first solution of the low-gravity free surface oscillation problem was given in 1964 by Satterlee and Reynolds for cylindrical containers [3].In the context of the Space Race, a significant effort was made to study low-gravity sloshing in cylindrical [3][4][5][6][7][8][9], spheroidal [8,10,11] or axisymmetric [12][13][14][15] tanks.A nonextensive list of modern works includes numerical models for cryogens [16,17], coupled non-linear implementations [18] or Computational Fluid Dynamics (CFD) simulations [19].Analytical solutions of the free and forced oscillations problem were found by Utsumi [20][21][22][23].
Different active and passive strategies have been traditionally employed to mitigate liquid sloshing in microgravity.Active approaches settle the propellant by imposing an adequate inertial force with a set of thrusters.Passive techniques make use of surface tension or membranes to hold the liquid at a certain position and reduce the effect of random accelerations.The resulting technical implementations, named Propellant Management Devices, are currently used to grant adequate liquid propellant feeding in case of in-orbit ignition of chemical propulsion units [24].
Since the absence of a settling volume force is the main characteristic of low-gravity sloshing, the problem could be attacked by reproducing the force of gravity with electromagnetic fields if the liquid can answer to such stimulus.The use of dielectrophoresis, a phenomenon on which a force is exerted on dielectric particles in the presence of a nonuniform electric field, was explored by the US Air Force with dielectric propellants in 1963 [25].The study unveiled a high risk of arcing inside the tanks and highlighted the need of large, heavy and noisy power sources.Approaches exploring Magnetic Positive Positioning have also been suggested to exploit the inherent magnetic properties of paramagnetic (oxygen) and diamagnetic (hydrogen) liquids.Relatively recent studies employed numerical simulations and microgravity experiments to validate this concept [26,27].
Ferrofluids are colloidal suspensions of magnetic nanoparticles treated with a surfactant to prevent from agglomeration.As a result, they exhibit high magnetic susceptibility.Their invention is attributed to Steve Papell, who in 1963 proposed to "provide an artificially imposed gravity environment" with ferrofluid-based magnetic propellants [28].The basic equations governing the dynamics of ferrofluids were presented in 1964 by Neuringer and Rosensweig [29], giving rise to the field of Ferrohydrodynamics [30].Although since then ferrofluids have found numerous applications on Earth, works addressing their original purpose are scarce.A rare exception is the NASA Magnetically Actuated Propellant Orientation (MAPO) experiment, that studied the magnetic positioning of liquid oxygen and validated a custom CFD model with a series of parabolic flight experiments with ferrofluids [31].Subsequent publications presented numerical models to study and generalize the measurements for space applications [32][33][34].
One of the main drawbacks of magnetic sloshing control is the rapid decay of magnetic fields, that limits its applicability to small and compact tanks.In this context, the increasing number of propelled microsatellites may benefit from this technology since a direct control of liquid sloshing could be achieved with small low-cost magnets.Once the liquid is positioned, magnetic fields could also be used to tune the natural frequencies and damping ratios of the system.This approach has been adopted for terrestrial applications, such as Tuned Magnetic Liquid Dampers [35,36].Onground research exploring axisymmetric sloshing [37,38], the frequency shifts due to the magnetic interaction [39], two-layer sloshing [40] or the swirling phenomenon [41], among others, has been carried out in the past with notable results.The sloshing of ferrofluids in low-gravity was indirectly studied in 1972 with a focus on gravity compensation techniques [42].
This work addresses the free and forced oscillations of magnetic liquids in axisymmetric containers when subjected to an external inhomogeneous magnetic field in microgravity.A ferrohydrodynamic model is developed to predict the natural frequencies and modal shapes of the system and a Álvaro Romero-Calvo et al.
2 JAM-19-1431 Fig. 1: Geometry of the system under analysis, composed by a magnetic liquid that fills a container in microgravity while subjected to an external magnetic field.S' and C' refer to the meniscus surface and contour, while S and C are the dynamic fluid surface and contour, respectively.O denotes the vertex of the meniscus, W is the vessel wall and V denotes the fluid volume.
case of application is presented.Unlike non-magnetic lowgravity sloshing, the presence of a restoring force ensures that the hypothesis of small oscillations (linear sloshing) is satisfied in a wide range of operation.

Problem formulation
The system to be modeled is represented in Fig. 1.A volume V of a magnetic liquid fills an upright axisymmetric tank with radius a at the meniscus contour.The liquid is incompressible, Newtonian and is characterized by a density ρ, specific volume v = ρ −1 , kinematic viscosity ν, surface tension σ and magnetization curve M(H).H and M are respectively the modules of the magnetic field H and magnetization field M, which are assumed to be collinear.The liquid meets the container wall with a contact angle θ c .An applied inhomogeneous axisymmetric magnetic field H 0 is imposed by an external source (e.g. a coil located at the base of the container).The inertial acceleration g along the z axis is also considered.A non-reactive gas at pressure p g fills the free space.In the figure, s is a curvilinear coordinate along the meniscus with origin in the vertex O and the relative heights are given by w (fluid surface -vertex), f (meniscus -vertex) and h (fluid surface -meniscus).The container is subjected to a lateral displacement x(t).The meniscus is represented by a dashed line and the dynamic fluid surface is given by a solid line.The model here presented extends the works by Satterlee & Reynolds [3] and Yeh [12] by considering the magnetic interaction and the axisymmetric oscillations case.

Nonlinear formulation
A cylindrical reference system {u r , u θ , u z }, centered at the vertex of the meniscus, is subsequently considered.If an irrotational flow field is assumed, there exists a potential ϕ such that being v the flow velocity with the subindices denoting the derivative.The velocity potential satisfies Laplace's equation subjected to the non-penetration wall boundary condition An additional boundary condition at the free surface is given by the unsteady ferrohydrodynamic Bernoulli's equation, which for an isothermal system with collinear magnetization field M adopts the form [30,43] where g is the inertial acceleration, ψ is the magnetic force potential, β(t) is an arbitrary function of time and p * is the composite pressure, defined as [30] with the first, second and third terms being named thermodynamic, magnetostrictive and fluid-magnetic pressure, respectively.For magnetically diluted systems M ∼ ρ, where ρ is represented by the magnetic particles concentration for the case of ferrofluids.Under this additional assumption, both pressure-like components are approximately compensated, and hence p * ≈ p(ρ, T ) [30].
The canonical magnetic force per unit volume is given by µ 0 M∇H, with µ 0 = 4π • 10 −7 N/A 2 being the permeability of free space [30].It can be shown that, for an isothermal fluid, this force derives from the potential [29] Due to the discontinuity of the Maxwell stress tensor at the magnetic liquid interface, the ferrohydrodynamic boundary condition in the absence of viscous forces becomes Álvaro Romero-Calvo et al.

JAM-19-1431
where p n = µ 0 M 2 n /2 is the magnetic normal traction, M n is the magnetization component normal to the fluid surface, and p c is the capillary pressure.The last is defined by the Laplace-Young equation p c = σK, where (8) is the curvature of the surface [1].Since at Eq. ( 4) only the spatial derivatives of the velocity potential have a physical meaning (e.g.Eq. ( 1)), any function of time can be added to ϕ whenever it is convenient.From a physical viewpoint, the absolute value of p remains undetermined under the incompressible flow assumption [43].The integration constant β(t) can be then absorbed into the definition of ϕ.By arbitrarily selecting β(t) = p g /ρ, the dynamic interface condition is obtained In an inertial reference system, the vertical displacement w of a surface point lying at (r, θ) in the interface z = w(r, θ,t) is given by If the velocity components relative to the tank dw/dt, dr/dt and rdθ/dt are expressed as a function of the potential given by Eq. ( 1), the kinematic interface condition that relates the last with the shape of the free surface is The continuity equation given by Eq. ( 2), together with the kinematic relation in Eq. ( 11) and the boundary conditions at Eq. (3) and Eq. ( 9), conform the problem to be solved after imposing the contact angle at the wall (θ c ) and a contact hysteresis parameter that will be described later in the text.

Equilibrium free surface shape
Due to the axisymmetry of geometry and loads, the static equilibrium surface of the fluid (S ) is also axisymmetric.Its shape can be easily determined from the balance of vertical forces in a circular crown of inner radius r and infinitesimal width dr along the surface [24].This results in the following set of dimensionless differential equations dF dS and boundary conditions where R = r/a, F = f /a, S = s/a, Bo = ρga 2 /σ is the Bond number, λ = a(p g − p 0 )/σ, being p 0 the liquid pressure at the free surface vertex, ψ includes the magnetic potential and magnetic normal traction through and the static contact angle with respect to the vertical θ c is given by A numerical solution can be easily computed by (1) setting an initial vertex position, (2) calculating the value of λ iteratively in order to satisfy the contact angle condition given by Eq. (12d), (3) solving the system with an ODE solver, and (4) obtaining the new height of the vertex through volume conservation.The procedure is repeated until the vertex height converges with a prescribed relative variation.When non-trivial magnetic setups are involved, a FEM simulation must be included in the loop.

Linear equations
The dynamic and kinematic conditions at Eq. ( 9) and Eq. ( 11) are highly nonlinear.The standard analytical procedure overcomes this difficulty by linearizing the problem and restricting the analysis to small oscillations.If the wave position is expressed as the sum of the static equilibrium solution and a small perturbation Álvaro Romero-Calvo et al.

JAM-19-1431
it will be possible to express the system of equations and boundary conditions as a Taylor's series expansion around the equilibrium surface S .If second-order terms are neglected, the boundary-value problem becomes Equation 16e assumes that the slope of the perturbation field at the wall is related to the magnitude of the perturbation at the same point through the parameter γ.The free-edge condition is characterized by γ = 0, while the stuck-edge condition is associated with γ → ∞ [12].This assumption is far from being rock-solid, and has indeed motivated a strong debate in the past.It has been suggested that the contact angle hysteresis condition depends not only on the position of the wave but also on its velocity [4] or the state of the wall [7].In the absence of a clear criteria, some studies assume the freeedge condition or intermediate approaches, generally obtaining a reasonable agreement with experimental data [10].
The only difference between the previous formulation and the classical problem without magnetic interactions is given by the magnetic term in Eq. (16c).The effective gravity acceleration includes both inertial and magnetic components and is given by where it can be observed that the magnetic contribution is a function of the radius.The magnitude and relative importance of the magnetic terms depend on the magnetic configuration and gravity level of the system under analysis.In particular, the magnetic component will be more significant in the absence of gravity.The magnetic field modifies the effective gravity acceleration of the system and shifts its natural frequencies, as reported in normal-gravity works [35,44].If the magnetic term was approximately constant in R, like in the case of a linear magnetic field and a flat surface, the problem would be equivalent to the non-magnetic system [42].In this analysis, however, an inhomogeneous magnetic field is being considered.

Extraction of tank motion
The potential ϕ has been referred to an inertial reference system.In order to analyze the movement of the free surface in the tank reference frame, it would be convenient to split this potential into rigid body (φ 0 ) and perturbed (φ) components, so that The boundary-value problem can be then expressed as a function of the perturbed potential and becomes 3 Free oscillations problem 3.1 Dimensionless linear equations In [3,12] it is proposed to split the potentials φ and h into spatial and temporary components, the second being a cyclic function of time with a circular frequency ω.The resulting dimensionless boundary-value problem is Álvaro Romero-Calvo et al.

JAM-19-1431 (20c
) where , Ω 2 = ρa 3 ω 2 /σ, Γ = aγ and g 0 is the acceleration of gravity at ground level [1].The Magnetic Bond Number has been defined as and accounts for the effects of the external magnetic field on the liquid.

Variational formulation
By following the procedure described in [3,12], it is possible to develop a variational principle equivalent to Eq. (20b) and Eq.(20c) as subjected to where G and its associated terms arise from the application of the wall boundary condition given by Eq. (20b) as detailed in [12].The obtention of this variational formulation follows the procedure described in [12,45].

Ritz method
The previous set of equations can only be analytically solved for simplified configurations in the absence of magnetic fields, like the case of a cylindrical container with flat bottom and flat fluid surface (θ c = 90 • ) [3].For other physical systems, Ritz approximations [1,12] or finite differences approaches [5,6] have been proposed to approximate the eigenfunctions of the problem.The basic formulation of the first approach is subsequently developed based on [3,12].
By following Ritz's method, the eigenfunctions Φ (n) can be approximated as the linear combination of admissible functions Φ i (R, θ, Z) that satisfy the boundary conditions of the problem described by Eq. (22b) to Eq. ( 22e).This results in where N is the size of the set of admissible functions.In the same way, the eigenfunctions H (n) i and G (n) are approximated by ζ i (R, θ) and ξ i (R, θ) through The sets of admissible functions are linked through Eq. ( 22c) and Eq.(22d).If i , the extremum condition represented by Eq. (22a) requires that ∂I ∂C (n) i = 0, (i=1,2,...,N), (26) which results in the system of equations (28d) The system has a nontrivial solution only when its determinant is zero.The eigenvalues Ω 2 n , and therefore the corresponding modal circular frequencies ω n , are then computed by means of the characteristic equation Once solved, the eigenfunctions of the problem are obtained from Eq. ( 23) to Eq. ( 25).

Forced lateral oscillations and mechanical analogies
In order to solve the forced lateral oscillations case, a modal solution for the linearized boundary-value problem given by Eq. ( 19) is built from the eigenmodes Φ (n) and H (n) .This solution is expressed as where the modal coordinates A n (t) and B n (t) are computed from their corresponding modal equations.The reader is referred to [12] for a full description of this procedure.
Since the modal equations are linear, the forced sloshing problem can be assimilated as the superposition of several linear oscillators.These are usually assumed to be a series of spring-mass systems on which a linear damper is included a posteriori.The employment of a mechanical analogy simplifies the integration of the sloshing problem into the equations of motion of the vehicle.In the case here analyzed, an additional magnetic pressure should be considered.However, in virtue of Newton's action-reaction principle, if the magnetic source is rigidly coupled to the tank then the distribution of magnetic pressures cannot produce torque in the assembly.That is, non-magnetic mechanical analogies can be extended for the magnetic case by simply employing the new magnetic eigenmodes Φ (n) and H (n) with their corresponding eigenfrequencies.Some possibilities are the models developed by Dodge [46] or Utsumi [23].

Selection of admissible and primitive functions
The set of admissible functions for Φ, H and G, related through Eq. (22c) and Eq.(22d), satisfy by definition Eq. (22b) to Eq. (22e) and form truncated series that approximate the eigenfunctions of the problem.A set of primitives should be previously defined as [3] with k p being the roots of the equation where J n is the Bessel function of first kind and order n.This index is used to study the axisymmetric (n = 0) or lateral (n = 1) case, while m defines the circumferential symmetry of the problem.Axisymmetric primitive functions will be characterized by n = m = 0, while lateral sloshing functions can be obtained with n = m = 1.However, the previous set of primitives does not satisfy Eq. (22e).The set of admissible functions is then created as a linear combination of the previous The N + 1 − p coefficients a ip for each i value are determined by imposing (i) a normalization condition, (ii) a contact angle value, and (iii) a Lagrange minimization problem designed to produce Bessel-like functions.These condition are respectively expressed as [3] Álvaro Romero-Calvo et al.
= 0 (j=i, i+1, ..., N, N+1), where λ 1i and λ 2i are the Lagrange multipliers of the minimization problem and Once the system is solved, the admissible set can be used to solve the eigenvalue problem.The success of this method depends on finding an adequate set of admissible functions Φ i such that the eigenfunctions Φ (n) can be represented with a reduced number of elements.The Z term in the primitives ϑ p , evaluated at the equilibrium surface, grows exponentially when F(R) departs significantly from Z = 0.This is the case of low Bond numbers and small contact angles.In [13] it is stated that for contact angles lower than 15 • in the case of freeedge condition (Γ = 0) or lower than 60 • for the stuckedge condition (Γ → ∞) the system may become numerically ill-conditioned.Furthermore, the comparison between this method and a finite differences approach showed significant divergences in the shape of the eigenfunctions Φ (n) for particular cases.
A possible solution would be finding a set of primitive functions without exponential terms.To the best knowledge of the authors, and considering the attempts made in [3], an alternative has not yet been proposed.It should also be noted that the magnetic force generally flattens the equilibrium surface, hence mitigating the effect of the exponential term in Eq. (31a).

Case of Application 4.1 System description
CubeSats are a particular class of nanosatellites composed of standardized 10×10 cm cubic units (U) [47].Their subsystems are designed to fit this standard, lowering the manufacturing costs and enhancing adaptability (e.g. the Aerojet MPS-130 propulsion module is offered in 1U or 2U formats).Due to the rapid decay of magnetic fields with distance, small spacecrafts with small propellant tanks may be particularly well suited for magnetic sloshing control implementations.Their development requires a dedicated feasibility study that is beyond the scope of this paper; however, the performance of an hypothetical system is here addressed.A 1U cylindrical container with 10 cm height, 5 cm radius and filled with a ferrofluid solution up to half of its height is subsequently considered.In order to generate a downward restoring force, a high-end cylindrical neodymium magnet magnetized at 1500 kA/m in the vertical direction is placed under the tank.The magnet has a hole of r i = 2.5 mm radius at the center (liquid outlet), a height h m , an external radius r e and a density ρ m = 7010 kg/m 3 .A contact angle of 67 • is assumed in microgravity conditions.The sketch of this setup is given in Fig. 2.
The liquid is a 1:10 volume solution of the Ferrotec water-based EMG-700 ferrofluid with a 0.58% volume concentration of magnetic nanoparticles.Its magnetic properties were measured with a MicroSense EZ-9 Vibrating Sample Magnetometer at the Physics Department of Politecnico di Milano.The corresponding magnetization curve is depicted in Fig. 3 and shows an initial susceptibility χ = 0.181 and saturation magnetization M s = 3160 A/m.Viscosity is assumed to have a negligible effect on the free sloshing problem.

Magnetic modeling
The magnetic system is modeled in Comsol Multiphysics, which is interfaced with the model developed in section 2. The FEM simulation is employed to estimate the fields H and M for a given magnet and equilibrium surface shape (meniscus).The last is computed iteratively by means of Eq. ( 12) with a FEM-in-the-loop implementation.Equation 21 is then employed to calculate the magnetic Bond number at the surface, that determines the solution of the free oscillation problem.The eigenfrequencies and eigenmodes of the system are finally obtained by solving Eq. (22).
To simulate the magnetic field, the model solves the stationary non-electric Maxwell equations where A is the magnetic vector potential produced by the magnetized materials.The constitutive relation is applied to the magnet with M = [0, 0, 1500] kA/m and to the surrounding air with M = 0.The magnetization curve M = f (H) given in Fig. 3 is applied to the ferrofluid volume.
The simulation domain is a rectangular 1×2 m region enclosing the container.An axisymmetric boundary condition is applied to the symmetry axis, while the tangential magnetic potential is imposed at the external faces through n × A = n × A 0 .A 0 is the dipole term of the magnetic vector potential generated by the magnetization fields of the magnet and ferrofluid.Consequently, A 0 is computed as the potential vector generated by two point dipoles applied at the Fig. 4: Mesh and overall dimensions of the magnetic FEM simulation domain.Units in mm.centers of the magnetization distributions and whose dipole moments are those of said distributions.While the dipole associated to the magnet can be calculated beforehand, the ferrofluid dipole needs to be approximated iteratively by integrating M in the ferrofluid volume.The relative error in the magnetic vector potential due the dipole approximation is estimated to be below 0.03% at the boundary of the domain with respect to the exact value generated by the equivalent circular loop.
The mesh is composed by 53000 irregular triangular elements and is refined at the meniscus, as shown in Fig. 4. Mean and minimum condition numbers of 0.983 and 0.766 are measured.
Figure 5 shows a particular configuration of analysis.Positive and negative curvatures are observed at the meniscus due to the high intensity of the magnetic field.Weaker magnets would result in convex equilibrium surfaces, as in the non-magnetic case.It should be noted that the magnetic Bond number rapidly decreases with distance to the source, spanning between 0 and 20 at the free surface.

Parametric analysis
Figure 6 depicts the fundamental sloshing frequency ω 1 , corresponding to the lowest root of Eq. ( 29), as a function of the external radius r e and height h m of the magnet.N = 7 admissible functions were employed in the computation.The mass of the magnet is given in a second scale, reflecting the technical trade-off between mass and sloshing frequency.Larger magnets result in stronger restoring forces and higher sloshing frequencies.For example, a downward 7 N force and a 100% increase in the fundamental sloshing frequency would be achieved with a 3 mm height, 30 mm external ra-Álvaro Romero-Calvo et al. dius and 60 grams magnet.Unlike the non-magnetic case, the presence of a significant restoring force ensures that the assumption of small oscillations (linear sloshing) is not violated for moderate displacements of the container.A small fluctuation in the frequency plot appears due to numerical errors.The procedures for solving Eq. ( 12) are highly dependent on the initial estimation of λ.This behavior has been extensively reported in the bibliography [3,24] Fig. 7: Magnetic and non-magnetic sloshing modes shape for the case under analysis and is further complicated by the magnetic interaction.In addition, the eigenvalue problem that provides the natural frequencies of the system becomes ill-conditioned if the actual modal shapes diverge significantly from the primitive functions [13].For strong magnetic fields, this may certainly be the case.The problem would be solved if a finite differences approach is employed instead of Ritz's method.
The modal shapes for the 3 mm height and 30 mm external radius magnet are represented in Fig. 7 together with the non-magnetic modes (i.e. the ones obtained when the magnet is removed).Although the profiles are essentially the same, it is interesting to observe how the fundamental mode slightly reduces and increases its displacement where the magnetic Bond number is greater and smaller, respectively.This is consistent with the aforementioned stabilizing role of the magnetic force.

Conclusions
A quasi-analytical model has been developed to study the sloshing of magnetic liquids in low-gravity conditions.The magnetic interaction modifies the shape of the meniscus and the effective inertial acceleration of the system as shown in Eq. ( 13) and Eq. ( 17), respectively.As a consequence, a shift of the eigenfrequencies and a modification of the eigenmodes of the free oscillations problem is produced.The framework here presented extends the models developed in previous works [3,12] by adding the magnetic and axisymmetric cases.
The small oscillations assumption is generally not valid for non-magnetic sloshing in microgravity, which is usually characterized by complex non-linear deformations driven by surface tension.In the magnetic case, however, the linear treatment of the problem is endorsed by a significant magnetic restoring force.It would be reasonable to ask how strong magnetic fields affect the shape of the eigenmodes and hence the reliability of Ritz's approximation.Since this procedure has been historically discussed [13], the computation of a finite-differences-based solution becomes particularly convenient for future implementations.Ritz's method represents, however, the simplest tool to solve the variational formulation given by Eq. ( 22) and has been presented for illustrative purposes.
From the engineering viewpoint, the magnetic sloshing concept represents an opportunity to develop new sloshing control devices for microsatellites.Unlike the non-magnetic case, the response of the system can be easily predicted, quantified and simulated by means of mechanical analogies.These simplified models can be easily embedded in a controller (e.g. a linear observer) used to predict and compensate the sloshing disturbances of a spacecraft in orbit.The spacecraft would then benefit from a significantly improved pointing performance.

Fig. 2 :
Fig. 2: Sketch of the case of application.Units in mm.

Fig. 5 :
Fig. 5: Equilibrium configuration of the magnetic sloshing damping system for h m = 5 mm and r e = 2.8 cm.The magnetic Bond number is represented in the color scale for the range of interest.

Fig. 6 :
Fig. 6: Fundamental sloshing frequency ω 1 (top) and mass of the magnet (bottom) as a function of the height h m and external radius r e of the magnet Modal coefficients used by Φ (n) , H (n) and G(n) i Internal radius of the cylindrical magnet Álvaro Romero-Calvo et al. 1 JAM-19-1431