Numerical simulation of the evolution of viscous rotating drops

In this contribution we simulate numerically the evolution of a viscous fluid drop rotating about a fixed axis at constant angular velocity ? or constant angular momentum L, surrounded by another viscous fluid. The problem is considered in the limit of large Ekman number and small Reynolds number. In the lecture we will describe the numerical method we have used to solve the PDE system that describes the evolution of the drop (3D boundary element method). We will also present the results we have obtained, paying special attention to the stability/instability of the equilibrium shapes.

A three-dimensional self-adaptive hp finite element method for the characterization of waveguide discontinuities

We propose the use of a highly-accurate three-dimensional (3D) fully automatic hp-adaptive finite element method (FEM) for the characterization of rectangular waveguide discontinuities. These discontinuities are either the unavoidable result of mechanical/electrical transitions or deliberately introduced in order to perform certain electrical functions in modern communication systems. The proposed numerical method combines the geometrical flexibility of finite elements with an accuracy that is often superior to that provided by semi-analytical methods. It supports anisotropic refinements on irregular meshes with hanging nodes, and isoparametric elements. It makes use of hexahedral elements compatible with high-order H(curl)H(curl) discretizations. The 3D hp-adaptive FEM is applied for the first time to solve a wide range of 3D waveguide discontinuity problems of microwave communication systems in which exponential convergence of the error is observed.

Multiple solutions and numerical analysis to the dynamic and stationary models coupling a delayed energy balance model involving latent heat and discontinuous albedo with a deep ocean

We study a climatologically important interaction of two of the main components of the geophysical system by adding an energy balance model for the averaged atmospheric temperature as dynamic boundary condition to a diagnostic ocean model having an additional spatial dimension. In this work, we give deeper insight than previous papers in the literature, mainly with respect to the 1990 pioneering model by Watts and Morantine. We are taking into consideration the latent heat for the two phase ocean as well as a possible delayed term. Non-uniqueness for the initial boundary value problem, uniqueness under a non-degeneracy condition and the existence of multiple stationary solutions are proved here. These multiplicity results suggest that an S-shaped bifurcation diagram should be expected to occur in this class of models generalizing previous energy balance models. The numerical method applied to the model is based on a finite volume scheme with nonlinear weighted essentially non-oscillatory reconstruction and Runge–Kutta total variation diminishing for time integration.

On the design of spherical gradient index lenses

Classical spherical gradient index (GRIN) lenses (such as Maxwell Fish Eye lens, Eaton lens, Luneburg lens, etc.) design procedure using the Abel integral equation is reviewed and reorganized. Each lens is fully defined by a function called the angle of flight which describes the ray deflection through the lens. The radial refractive index distribution is obtained by applying a linear integral transformation to the angle of flight. The interest of this formulation is in the linearity of the integral transformation which allows us to derive new solutions from linear combinations of known lenses. Beside the review of the classical GRIN designs, we present a numerical method for GRIN lenses defined by the Abel integral equation with fixed limits, which is an ill-posed problem.

Design and kinematic analysis of 3PSS-1S wrist for needle insertion guidance

In this work it is presented a complete kinematic analysis of the 3PSS-1S parallel mechanism for its implementation as a spherical wrist for a needle insertion guidance robot. The spherical 3PSS-1S mechanism is a low weight and reduced dimension parallel mechanism that allows spherical movements providing the requirements needed for the serial–parallel robotic arm for needle insertion guidance. The solution of its direct kinematic is computed with a numerical method based on the Newton–Raphson formulation and a constraint function of the mechanism. The input–output velocity equation is obtained with the use of screw theory. Three types of singular postures are identified during simulations and verified in the real prototype. The 3PSS-1S can perform pure rotations of ±45°±45°, ±40°±40°, ±60°±60° along the View the MathML sourcex, View the MathML sourcey, View the MathML sourcez axes respectively.

Deformabilidad de escolleras empleadas en banquetas de cimentación de cajones portuarios

En esta tesis doctoral se presenta una investigación sobre el comportamiento deformacional de las escolleras empleadas en banquetas de cimentación de obras portuarias de cajones. El trabajo aborda el estudio de la deformabilidad de escolleras portuarias combinando (i) investigación mediante ensayos de laboratorio; (ii) análisis del comportamiento in situ de las banquetas de escolleras y (iii) cálculos realizados con modelos numéricos. Se expone en primer lugar la investigación experimental realizada en el Laboratorio de Geotecnia del Centro de Estudios y Experimentación de Obras Públicas (CEDEX) para estudiar la deformabilidad de las escolleras mediante ensayos a gran escala, habida cuenta las grandes dimensiones de las partículas de escollera. Se ha tratado de establecer una metodología de ensayo que reproduzca las solicitaciones de las escolleras colocadas en banquetas de cimentación de obras portuarias. Asimismo, se ha hecho una interpretación exhaustiva de los resultados de los ensayos con el fin de establecer unos valores que caractericen la deformabilidad de las escolleras analizadas. Es posible dar un intervalo de valores de la compresibilidad de las escolleras portuarias que, dada la escasez de literatura existente, constituyen unos valores de referencia. Asimismo, se ha propuesto una metodología para para estimar la deformabilidad de escolleras portuarias in situ. La información disponible ha permitido realizar estudios de la deformabilidad in situ en dos muelles españoles con semejanzas estructurales y constructivas. La interpretación conjunta de los resultados ha sugerido unos valores de deformabilidad in situ. Conviene destacar que la práctica ausencia de rangos de valores de compresibilidad in situ para estos rellenos empleados en obras portuarias pone de manifiesto la importancia de los resultados obtenidos. Evidencias de diferencias de comportamiento de las escolleras empleadas en banquetas de cimentaci ón de cajones portuarios en laboratorio e in situ han sido documentadas. La evaluación conjunta del comportamiento tenso-deformacional de las escolleras en laboratorio e in situ ha estimulado la búsqueda de una correlación entre la compresibilidad de las escolleras en ambos escenarios. Finalmente, se ha elaborado un modelo numérico con la formulación matemática del método sincrético (Perucho (2004, 2008)) que supone una opción interesante para evaluar la deformabilidad de los rellenos granulares. En la práctica, el empleo del modelo sincrético requiere la determinación de unos microparámetros. La disponibilidad de numerosos resultados de laboratorio realizados en las escolleras portuarias ha permitido calibrar el modelo realizado. De esta manera, se dispone de una herramienta de cálculo para evaluar la deformabilidad de los relleno granulares con un método numérico. The focus of this Thesis is to explore the deformational behavior of large rock fill materials used as rock mattress foundations for gravity caissons structures. The determination of the compressibility of large granular media focuses on (i) laboratory testing, (ii) in situ performance analysis of rock mattress foundations for caissons, and (iii) numerical modelling. First, the results of the large-scale laboratory research program, conducted at the Geotechnical Laboratory for the Center for Studies and Experimentation for Public Works (CEDEX), to determine the deformability of large rock fill materials is presented. The testing procedure was specifically designed to reproduce the loading sequence of in situ rubble mound foundations. A thoughtful analysis of the laboratory testing results suggests a range of compressibility for large granular media. The lack of currently available information regarding large rock fill deformability places a certain emphasis on the results of the testing program. Second, the results of this research includes a procedure for evaluating in situ rock fill deformational behavior. Data, collected from monitoring two caisson-type quays in Spain, provides information to study in situ rock mattress foundations. Careful interpretation of in situ data reveals a range of deformability of rock mattress foundations in caisson-type quays. Based upon a review of available literature, assessments on the behavior of rock mattress foundations for caissons using in situ analysis are quite limited. The data from this research are likely to contribute to the knowledge of the in situ behavior of rock mattress foundations for caissons. Additionally, findings indicate an appreciable variation between the laboratory and the in situ behaviour of materials from rock mattress foundations for caissons. Dissimilarities between laboratory and in situ moduli of deformation are examined in detail. Correlations between laboratory and in situ values are made. Finally, numerical modeling, based upon the research of Perucho (2004, 2008), is presented to predict the deformation behavior of large granular media. The determination of microparameters that control macropropierties requires extensive calibration effort. The calibration process was carried out using the results of large-scale laboratory testing available from previous analysis. The presented numerical method is both versatile and attractive as it reasonably predicts the compressibility of large rock fill materials.

PIC computation of electron current collection to a moving bare tether in the mesothermal condition

In tethered satellite technology, it is important to estimate how many electrons a spacecraft can collect from its ambient plasma by a bare electrodynamic tether. The analysis is however very difficult because of the small but significant Geo-magnetic field and the spacecraft’s relative motion to both ions and electrons. The object of our work is the development of a numerical method, for this purpose. Particle-In-Cell (PIC) method, for the calculation of electron current to a positive bare tether moving at orbital velocity in the ionosphere, i.e. in a flowing magnetized plasma under Maxwellian collisionless conditions. In a PIC code, a number of particles are distributed in phase space and the computational domain has a grid on which Poisson equation is solved for field quantities. The code uses the quasi-neutrality condition to solve for the local potential at points in the plasma which coincide with the computational outside boundary. The quasi-neutrality condition imposes ne - ni on the boundary. The Poisson equation is solved in such a way that the presheath region can be captured in the computation. Results show that the collected current is higher than the Orbital Motion Limit (OML) theory. The OML current is the upper limit of current collection under steady collisionless unmagnetized conditions. In this work, we focus on the flowing effects of plasma as a possible cause of the current enhancement. A deficit electron density due to the flowing effects has been worked and removed by introducing adiabatic electron trapping into our model.

Simulación numérica computacional de fenómenos dinámicos

En este proyecto se trata la simulación numérica de un fenómeno dinámico, basado en el comportamiento de una onda transmitida a lo largo de una cuerda elástica de un instrumento musical, cuyos extremos se encuentran anclados. El fenómeno físico, se desarrolla utilizando una ecuación en derivadas parciales hiperbólicas con variables espacial y temporal, acompañada por unas condiciones de contorno tipo Dirichlet en los extremos y por más condiciones iniciales que dan comienzo al proceso. Posteriormente se han generado algoritmos para el método numérico empleado (Diferencias finitas centrales y progresivas) y la programación del problema aproximado con su consistencia, estabilidad y convergencia, obteniéndose unos resultados acordes con la solución analítica del problema matemático. La programación y salida de resultados se ha realizado con Visual Studio 8.0. y la programación de objetos con Visual Basic .Net In this project the topic is the numerical simulation of a dynamic phenomenon, based on the behavior of a transmitted wave along an elastic string of a musical instrument, whose ends are anchored. The physical phenomenon is developed using a hyperbolic partial differential equation with spatial and temporal variables, accompanied by a Dirichlet boundary conditions at the ends and more initial conditions that start the process. Subsequently generated algorithms for the numerical method used (central and forward finite differences) and the programming of the approximate problem with consistency, stability and convergence, yielding results in line with the analytical solution of the mathematical problem. Programming and output results has been made with Visual Studio 8.0. and object programming with Visual Basic. Net

Modelling of short runout propagation landslides and debris flows

This paper proposes an extension of methods used to predict the propagation of landslides having a long runout to smaller landslides with much shorter propagation distances. The method is based on: (1) a depth-integrated mathematical model including the coupling between the soil skeleton and the pore fluids, (2) suitable rheological models describing the relation between the stress and the rate of deformation tensors for fluidised soils and (3) a meshless numerical method, Smooth Particle Hydrodynamics, which separates the computational mesh (or set of computational nodes) from the mesh describing the terrain topography, which is of structured type – thus accelerating search operations. The proposed model is validated using two examples for which there are analytical solutions, and then it is applied to two short runout landslides which happened in Hong Kong in 1995, for which there is available information.

Analysis of Meshfree Methods for Lagrangian Fluid-Structure Interaction

In this dissertation a new numerical method for solving Fluid-Structure Interaction (FSI) problems in a Lagrangian framework is developed, where solids of different constitutive laws can suffer very large deformations and fluids are considered to be newtonian and incompressible. For that, we first introduce a meshless discretization based on local maximum-entropy interpolants. This allows to discretize a spatial domain with no need of tessellation, avoiding the mesh limitations. Later, the Stokes flow problem is studied. The Galerkin meshless method based on a max-ent scheme for this problem suffers from instabilities, and therefore stabilization techniques are discussed and analyzed. An unconditionally stable method is finally formulated based on a Douglas-Wang stabilization. Then, a Langrangian expression for fluid mechanics is derived. This allows us to establish a common framework for fluid and solid domains, such that interaction can be naturally accounted. The resulting equations are also in the need of stabilization, what is corrected with an analogous technique as for the Stokes problem. The fully Lagrangian framework for fluid/solid interaction is completed with simple point-to-point and point-to-surface contact algorithms. The method is finally validated, and some numerical examples show the potential scope of applications.

Dynamic analysis of payloads and structures with intermediate modal density

La necesidad de desarrollar técnicas para predecir la respuesta vibroacústica de estructuras espaciales lia ido ganando importancia en los últimos años. Las técnicas numéricas existentes en la actualidad son capaces de predecir de forma fiable el comportamiento vibroacústico de sistemas con altas o bajas densidades modales. Sin embargo, ambos rangos no siempre solapan lo que hace que sea necesario el desarrollo de métodos específicos para este rango, conocido como densidad modal media. Es en este rango, conocido también como media frecuencia, donde se centra la presente Tesis doctoral, debido a la carencia de métodos específicos para el cálculo de la respuesta vibroacústica. Para las estructuras estudiadas en este trabajo, los mencionados rangos de baja y alta densidad modal se corresponden, en general, con los rangos de baja y alta frecuencia, respectivamente. Los métodos numéricos que permiten obtener la respuesta vibroacústica para estos rangos de frecuencia están bien especificados. Para el rango de baja frecuencia se emplean técnicas deterministas, como el método de los Elementos Finitos, mientras que, para el rango de alta frecuencia las técnicas estadísticas son más utilizadas, como el Análisis Estadístico de la Energía. En el rango de medias frecuencias ninguno de estos métodos numéricos puede ser usado con suficiente precisión y, como consecuencia -a falta de propuestas más específicas- se han desarrollado métodos híbridos que combinan el uso de métodos de baja y alta frecuencia, intentando que cada uno supla las deficiencias del otro en este rango medio. Este trabajo propone dos soluciones diferentes para resolver el problema de la media frecuencia. El primero de ellos, denominado SHFL (del inglés Subsystem based High Frequency Limit procedure), propone un procedimiento multihíbrido en el cuál cada subestructura del sistema completo se modela empleando una técnica numérica diferente, dependiendo del rango de frecuencias de estudio. Con este propósito se introduce el concepto de límite de alta frecuencia de una subestructura, que marca el límite a partir del cual dicha subestructura tiene una densidad modal lo suficientemente alta como para ser modelada utilizando Análisis Estadístico de la Energía. Si la frecuencia de análisis es menor que el límite de alta frecuencia de la subestructura, ésta se modela utilizando Elementos Finitos. Mediante este método, el rango de media frecuencia se puede definir de una forma precisa, estando comprendido entre el menor y el mayor de los límites de alta frecuencia de las subestructuras que componen el sistema completo. Los resultados obtenidos mediante la aplicación de este método evidencian una mejora en la continuidad de la respuesta vibroacústica, mostrando una transición suave entre los rangos de baja y alta frecuencia. El segundo método propuesto se denomina HS-CMS (del inglés Hybrid Substructuring method based on Component Mode Synthesis). Este método se basa en la clasificación de la base modal de las subestructuras en conjuntos de modos globales (que afectan a todo o a varias partes del sistema) o locales (que afectan a una única subestructura), utilizando un método de Síntesis Modal de Componentes. De este modo es posible situar espacialmente los modos del sistema completo y estudiar el comportamiento del mismo desde el punto de vista de las subestructuras. De nuevo se emplea el concepto de límite de alta frecuencia de una subestructura para realizar la clasificación global/local de los modos en la misma. Mediante dicha clasificación se derivan las ecuaciones globales del movimiento, gobernadas por los modos globales, y en las que la influencia del conjunto de modos locales se introduce mediante modificaciones en las mismas (en su matriz dinámica de rigidez y en el vector de fuerzas). Las ecuaciones locales se resuelven empleando Análisis Estadístico de Energías. Sin embargo, este último será un modelo híbrido, en el cual se introduce la potencia adicional aportada por la presencia de los modos globales. El método ha sido probado para el cálculo de la respuesta de estructuras sometidas tanto a cargas estructurales como acústicas. Ambos métodos han sido probados inicialmente en estructuras sencillas para establecer las bases e hipótesis de aplicación. Posteriormente, se han aplicado a estructuras espaciales, como satélites y reflectores de antenas, mostrando buenos resultados, como se concluye de la comparación de las simulaciones y los datos experimentales medidos en ensayos, tanto estructurales como acústicos. Este trabajo abre un amplio campo de investigación a partir del cual es posible obtener metodologías precisas y eficientes para reproducir el comportamiento vibroacústico de sistemas en el rango de la media frecuencia. ABSTRACT Over the last years an increasing need of novel prediction techniques for vibroacoustic analysis of space structures has arisen. Current numerical techniques arc able to predict with enough accuracy the vibro-acoustic behaviour of systems with low and high modal densities. However, space structures are, in general, very complex and they present a range of frequencies in which a mixed behaviour exist. In such cases, the full system is composed of some sub-structures which has low modal density, while others present high modal density. This frequency range is known as the mid-frequency range and to develop methods for accurately describe the vibro-acoustic response in this frequency range is the scope of this dissertation. For the structures under study, the aforementioned low and high modal densities correspond with the low and high frequency ranges, respectively. For the low frequency range, deterministic techniques as the Finite Element Method (FEM) are used while, for the high frequency range statistical techniques, as the Statistical Energy Analysis (SEA), arc considered as more appropriate. In the mid-frequency range, where a mixed vibro-acoustic behaviour is expected, any of these numerical method can not be used with enough confidence level. As a consequence, it is usual to obtain an undetermined gap between low and high frequencies in the vibro-acoustic response function. This dissertation proposes two different solutions to the mid-frequency range problem. The first one, named as The Subsystem based High Frequency Limit (SHFL) procedure, proposes a multi-hybrid procedure in which each sub-structure of the full system is modelled with the appropriate modelling technique, depending on the frequency of study. With this purpose, the concept of high frequency limit of a sub-structure is introduced, marking out the limit above which a substructure has enough modal density to be modelled by SEA. For a certain analysis frequency, if it is lower than the high frequency limit of the sub-structure, the sub-structure is modelled through FEM and, if the frequency of analysis is higher than the high frequency limit, the sub-structure is modelled by SEA. The procedure leads to a number of hybrid models required to cover the medium frequency range, which is defined as the frequency range between the lowest substructure high frequency limit and the highest one. Using this procedure, the mid-frequency range can be define specifically so that, as a consequence, an improvement in the continuity of the vibro-acoustic response function is achieved, closing the undetermined gap between the low and high frequency ranges. The second proposed mid-frequency solution is the Hybrid Sub-structuring method based on Component Mode Synthesis (HS-CMS). The method adopts a partition scheme based on classifying the system modal basis into global and local sets of modes. This classification is performed by using a Component Mode Synthesis, in particular a Craig-Bampton transformation, in order to express the system modal base into the modal bases associated with each sub-structure. Then, each sub-structure modal base is classified into global and local set, fist ones associated with the long wavelength motion and second ones with the short wavelength motion. The high frequency limit of each sub-structure is used as frequency frontier between both sets of modes. From this classification, the equations of motion associated with global modes are derived, which include the interaction of local modes by means of corrections in the dynamic stiffness matrix and the force vector of the global problem. The local equations of motion are solved through SEA, where again interactions with global modes arc included through the inclusion of an additional input power into the SEA model. The method has been tested for the calculation of the response function of structures subjected to structural and acoustic loads. Both methods have been firstly tested in simple structures to establish their basis and main characteristics. Methods are also verified in space structures, as satellites and antenna reflectors, providing good results as it is concluded from the comparison with experimental results obtained in both, acoustic and structural load tests. This dissertation opens a wide field of research through which further studies could be performed to obtain efficient and accurate methodologies to appropriately reproduce the vibro-acoustic behaviour of complex systems in the mid-frequency range.

Optimal laser intensity profiles for a uniform target illumination in direct-drive inertial confinement fusion

A numerical method providing the optimal laser intensity profiles for a direct-drive inertial confinement fusion scheme has been developed. The method provides an alternative approach to phase-space optimization studies, which can prove computationally expensive. The method applies to a generic irradiation configuration characterized by an arbitrary number NB of laser beams provided that they irradiate the whole target surface, and thus goes beyond previous analyses limited to symmetric configurations. The calculated laser intensity profiles optimize the illumination of a spherical target. This paper focuses on description of the method, which uses two steps: first, the target irradiation is calculated for initial trial laser intensities, and then in a second step the optimal laser intensities are obtained by correcting the trial intensities using the calculated illumination. A limited number of example applications to direct drive on the Laser MegaJoule (LMJ) are described.

Flight dynamics and stability of kites in steady and unsteady wind conditions

The flight dynamics and stability of a kite with a single main line flying in steady and unsteady wind conditions are discussed. A simple dynamic model with five degrees of freedom is derived with the aid of Lagrangian formulation, which explicitly avoids any constraint force in the equations of motion. The longitudinal and lateral–directional modes and stability of the steady flight under constant wind conditions are analyzed by using both numerical and analytical methods. Taking advantage of the appearance of small dimensionless parameters in the model, useful analytical formulas for stable-designed kites are found. Under nonsteady wind-velocity conditions, the equilibrium state disappears and periodic orbits occur. The kite stability and an interesting resonance phenomenon are explored with the aid of a numerical method based on Floquet theory.

Stability of liquid bridges between an elliptical and a circular supporting disk

A numerical method has been developed to determine the stability limits for liquid bridges held between noncircular supporting disks and the application to a configuration with a circular and an elliptical disk subjected to axial acceleration has been made. The numerical method led to results very different from the available analytical solution which has been revisited and a better approximation has been obtained. It has been found that just retaining one more term in the asymptotic analysis the solution reproduces the real behavior of the configuration and the numerical results.

Estabilidad de puentes líquidos bajo el efecto de una rotación no axilsimétrica

El estudio de la influencia de perturbaciones de distinta naturaleza en configuraciones de puentes líquidos apoyados en dos discos coaxiales en rotación encuentra una importante motivación en el uso de dicha configuración en la fabricación de cristales semiconductores ultra-puros por la denominada técnica de zona flotante, en la que la rotación de los discos se utiliza para alcanzar temperaturas uniformes. El presente estudio muestra los resultados obtenidos mediante la aplicación de un método numérico en el análisis de la estabilidad de puentes líquidos en isorrotación sometidos al efecto de una fuerza axial uniforme (gravedad axial) y una excentricidad entre el eje de giro y el eje de los discos. Se analiza el efecto de la aplicación de estos factores tanto de forma conjunta como por separado. Aunque existen numerosos estudios previos sobre puentes líquidos sometidos a diversos efectos, el análisis del efecto combinado de la rotación con excentricidad y gravedad axial no ha sido realizado con anterioridad. Este estudio permite además entender los resultados del experimento a bordo de la misión TEXUS-23, en el que un puente líquido sujeto entre dos discos circulares y coaxiales es sometido al efecto de una rotación creciente en torno a un eje desplazado respecto al eje de los discos. Aunque en el experimento no se impone una fuerza axial controlada, la desestabilización y rotura del puente se produce de forma notablemente asimétrica, lo que no puede ser explicado con los estudios precedentes y sugiere una posible presencia de una aceleración axial residual. Se ha desarrollado por tanto un método de análisis de imágenes que permite comparar las formas obtenidas en el experimento con las calculadas numéricamente. En este estudio se muestran los detalles del procesado realizado en las imágenes de la misión TEXUS-23, y los resultados de su comparación con el análisis numérico, que permiten determinar el valor de la gravedad axial que mejor reproduce los resultados del experimento. Estos resultados ponen de manifiesto la importancia del conocimiento y la modelización de efectos cuya presencia (intencionada o no) afectan de forma visible a la estabilidad y la morfología de los puentes líquidos. ABSTRACT The study of the influence of various disturbances in configurations consisting of a liquid bridge supported by two co-axial disks in rotation has an important motivation in the use of this configuration in the fabrication of ultrapure semiconductor crystals via the so-called floating zone technique, in which the rotation of the disks is used to achieve a uniform temperature field. The present study shows the results obtained through the application of a numerical method in the analysis of the stability of liquid bridges in isorotation under the effect of a uniform axial force field (axial gravity) and an offset between the rotation axis and the axis of the supporting disks (eccentricity). The analysis studies the effect of both the combined and separate application of these factors. Although there are numerous studies on liquid bridges subject to various effects, the analysis of the combined effect of rotation with eccentricity and axial gravity has not been done before. Furthermore, this study allows us to understand the results from the experiment aboard the TEXUS-23 mission, in which a liquid bridge supported between two circular-shaped, co-axial disks is subject to the effect of an increasing rotation around an axis with an offset with respect to the axis of the disks. Although the experiment conditions do not include a controlled axial force field, the instability and breakage of the bridge occurs with a marked asymmetry, which cannot be explained by previous studies and suggests the possible presence of a residual axial gravity. Therefore, an image analysis method has been developed which allows to compare the shapes obtained in the experiment with those calculated with the numerical method. This study shows the details of the processing performed on the images from the TEXUS-23 mission and the results from their comparison with the numerical analysis, which allow to determine the axial gravity value which best recovers the experimental results. These results highlight the importance of the understanding and modelling of effects which, when present (intentionally or not), noticeably affect the stability and shape of the liquid bridges.