The target of this work is to contribute to the enhancement of numerical methods for the simulation of complex thermal systems. Frequently, the factor that limits the accuracy of the simulations is the computing power: accurate simulations of complex devices require fine three-dimensional discretizations and the solution of large linear equation systems.Their efficient solution is one of the central aspects of this work. Low-cost parallel computers, for instance, PC clusters, are used to do so. The main bottle-neck of these computers is the notwork, that is too slow compared with their floating-point performance.Before considering linear solution algorithms, an overview of the mathematical models used and discretization techniques in staggered cartesian and cylindrical meshes is provided. The governing Navier-Stokes equations are solved using an implicit finite control volume method. Pressure-velocity coupling is solved with segregated approaches such as SIMPLEC.Different algorithms for the solution of the linear equation systems are reviewed: from incomplete factorizations such as MSIP, Krylov solvers such as BICGSTAB and GMRESR to acceleration techniques such as the Algebraic Multi Grid and the Multi Resolution Analysis with wavelts. Special attention is paid to preconditioned Krylov solvers for their application to parallel CFD problems.The fundamentals of parallel computing in distributed memory computers as well as implemetation details of these algorithms in combination with the domain decomposition method are given. Two different distributed memory computers, a Cray T3E and a PC cluster are used for several performance measures, including network throughput, performance of algebraic subroutines that affect to the overall efficiency of algorithms, and the solver performance. These measures are addressed to show the capabilities and drawbacks of parallel solvers for several processors and their partitioning configurations for a problem model.Finally, in order to illustrate the potential of the different techniques presented, a three-dimensional CFD problem is solved using a PC cluster. The numerical results obtained are validated by comparison with other authors. The speedup up to 12 processors is measured. An analysis of the computing time shows that, as expected, most of the computational effort is due to the pressure-correction equation,here solved with BiCGSTAB. The computing time algorithm , for different problem sizes, is compared with Schur-Complement and Multigrid. / El trabajo de tesis se centra en la solución numérica de las ecuaciones de navier-Stokes en regimen transitorio, tridimensional y laminar. Los algoritmos utilizados son del tipo segregado (SIMPLEC)y se basan en el uso de técnicas de volumenes finitos, con mallas estructurales del tipo staggered y discretizaciones temporales implícitas. En este contexto, el pricipal, problema son los elevados tiempos de cálculo de las simulaciones, que en buena parte se deben a la solución de los sistemas de ecuaciones lineales. Se hace una revisión de diferentes métodos utilizados típicamente en ordenadores secuenciales: GMRES, BICGSTAB, ACM, MSPIP. A fin de reducir los tiempos de cálculo se emplean ordenadores paralelos de memoria distribuida, basados en la agrupacion de ordenadores personales convencionales (PC clusters). Por lo que respecta a la potencia de cálculo por procesador, estos sistemas son comparables a los ordenadores paralelos de memoria distribuida convencionales (como el Cray T3E) siendo, su principal problema la baja capacidad de comunicación (elevada latencia, bajo ancho de banda). Este punto condiciona toda la estrategia computacional, obligando a reducir al máximo el número y el tamaño de los mensajes intercambiados. Este aspecto se cuantifica detalladamente en la tesis, realizando medidas de tiempos de cálculo en ambos ordenadores para diversas operaciones críticas para los algoritmos lineales. Tambien se miden y comparan los tiempos de cálculo y speed ups obtenidos en la solución de los sistemas lineales con diferentes algoritmos paralelos (Jacobi, MSIP, GMRES, BICGSTAB) y para diferentes tamaños de malla. Finalmente, se utilizan las técnicas anteriores para resolver el caso denominado driven cavity, en situacionies tridimensionales y con numeros de Reynolds de hasta 8000. Los resultados obtenidos se utilizan para validar los códigos desarrollados, en base a resultados de otros códigos y también se basa en la comparación con resultados experimentales procedentes de la bibliografía. Se utilizan hasta 12 procesadores, obteniendose spped ups de hasta 9.7 en el cluster de PCs. Se analizan los tiempos de cálculo de cada fase del código, señalandose areas para futuras mejoras. Se comparan los tiempos de cálculo con los algoritmos implementados en otros trabajos. La conclusión final es que los clusters de PCs son una plataforma de gran potencia en los cálculos de dinámica de fluidos computacional.

La tesis analiza el comportamiento anómalo de una servoválvula instalada en una máquina de electroerosión operando bajo ciertas condiciones de trabajo. En el primer capítulo de la tesis se describen una serie de experimentos llevados a cabo para: en primer lugar, llegar a obtener bajo que condiciones de funcionamiento se producen las inestabilices (vibraciones audibles) y en que parte de la máquina se originan. Como resultado de estas primeras investigaciones se obtiene que el sistema de transmisión de energía oleohidráulica es el generador de las perturbaciones. Una segunda fase de experimentos descritos asimismo en el primer capítulo, llevan a poder afirmar que las perturbaciones se originan en la servoválvula.El siguiente paso consistió en profundizar sobre las prestaciones de la servoválvula, analizando inicialmente los elementos de que consta y sus características teóricas y experimentales. Paralelamente se llevo a cabo un estudio crítico de la bibliografía sobre el tema, en particular sobre posibles anomalías de funcionamiento en servoválvulas. Algunos investigadores insinuaron que las fuerzas de flujo sobre las toberas eran una fuente de inestabilidades. Nuestras investigaciones confirman esta hipótesis y ponen a la luz la importancia del modelo de flujo tobera/paleta en este sentido. Se comprueba que los coeficientes de descarga y la distancia tobera paleta son definitorios en el modelo de flujo.En el tercer capítulo se estudian en profundidad los diversos modelos de flujo, se hace hincapié que en la servoválvula bajo estudio hay dos tipos de flujo: el denominado tobera paleta y el denominado paleta tobera. Se observa experimentalmente diferencias importantes entre los coeficientes de descarga para cada uno de estos dos tipos de flujo. En este mismo capítulo se realiza una modelización bidimensional por volúmenes finitos de la evolución del fluido dentro de una servoválvula, en donde se pueden diferenciar los diversos modelos de flujo en función de la distancia tobera paleta y la diferencia de presión entre la entrada y salida. Al final de este capítulo se comparan los resultados experimentales con los hallados vía modelización, observando un buen comportamiento del modelo. El capitulo cuarto se centra en la obtención de un modelo que simule el comportamiento de la servoválvula. La modelización se realiza utilizando el método de diagramas de enlaces "Bond Graph".Como innovación en este apartado cabe destacar la modelización de las fuerzas de flujo y la variabilidad de los coeficientes de descarga en función del número de Reynolds. Como complemento a la modelización y en lo que a las fuerzas de flujo hace referencia, se realiza un estudio comparativo de los dos modelos de servoválvula de cuatro toberas existentes: el Williams y el Moog 92, comprobando que el modelo actual, el Moog 92 es mas tolerante a las inestabilidades.De los resultados de la modelización se puede estimar con gran precisión el comportamiento en régimen permanente y en régimen transitorio de la servoválvula. Se demostró asimismo que la frecuencia propia fundamental de la servoválvula estaba alrededor de los 580 Hz. Esta frecuencia se corresponde con la obtenida en forma de espectro audible en los experimentos iniciales.En el capítulo cinco se analiza con detalle el origen de la perturbación audible de alrededor de 2200 Hz. Después de una serie de experimentos y de analizar los resultados que otros investigadores obtuvieron en sistemas análogos, se llega a la conclusión de que la perturbación de alrededor de 2200 Hz es de origen fluidodinámico, y capaz de excitar uno de los modos de vibración de la armadura. Esta inestabilidad se ve reforzada por un fenómeno de resonancia que es consecuencia de la morfología de la tobera / paleta.Como innovaciones de esta tesis cabe citar:1.- El estudio del flujo entrante paleta/tobera. En este sentido se han definido los modelos de flujo que tienen lugar y los parámetros que los definen: coeficientes de descarga en función de la distancia tobera paleta. Se ha observado que el origen de todas la inestabilidades se localizan cuando el flujo va de paleta a tobera, (flujo entrante). Estos modelos de flujo fueron modelizados mediante volúmenes finitos siendo los resultados obtenidos satisfactorios.2.- La modelización del comportamiento en régimen permanente y régimen transitorio de la servoválvula mediante técnicas de diagramas de enlaces (Bond Graph).3.- La metodología experimental seguida para poner en evidencia que fenómenos son el origen de las perturbaciones audibles, y que factores deben considerarse para minimizar sus efectos. / The PhD focuses on the analysis of the erratic behaviour of the oil hydraulic system installed in a sparkling machine and when operating under some particular conditions.The first chapter describes a set of experiments to find out under which conditions and were the instabilities (audible vibrations) are created. As a result, it is demonstrated that the servovalve is the origin of the instabilities. Some researchers pointed out that flow forces acting onto the servovalve flappers could originate vibrations. It is then decided to measure the discharge coefficients for flow going from nozzle to flapper and from flapper to nozzle, being the aim to better understand the flow models inside the four nozzle two flapper servovalve under study.The discharge coefficients demonstrate that two completely different models of flow appear inside the servovalve. To better understand such flow models a CFD computer modelisation of the nozzle-flapper and flapper-nozzle flows was undertaken, finding out that flow in direction flapper nozzle is prone to be instable. Good agreement was found when comparing discharge coefficients found experimentally and via CFD modelisation. The fourth chapter is dedicated to simulate the complete servovalve behaviour, the simulation is performed using the Bond Graph technique. It is interesting to point out that when the armature of the servovalve was excited with a frequency near the servovalve natural frequency the model gave an over excitation of the servovalve, as it occurs in reality. Such frequency was one of the perturbation peaks measured at the beginning of the study. On the last chapter is analysed experimentally the origin of a second perturbation peak with a frequency of 2200 Hz. Finding out that such instability has a fluid dynamic origin and is created by the flow going from flapper to nozzle. Such perturbation is amplified by the nozzle chamber and it is able to excite some vibration modes of the servovalve armature. The way to overcome such problem would be via preventing the amplification from being created. This main mean whether avoid flow from flapper to nozzle or modifying the nozzle chamber in the nozzles having flapper nozzle flow.

El desarrollo de esta tesis se centra en el estudio del sincronizado de actuadores lineales oleohidráulicos mediante válvulas divisoras de caudal.La mayoría de los estudios dedicados al análisis del comportamiento dinámico de una válvula divisora plantean las ecuaciones de estado (sistema de ecuaciones diferenciales no lineales) y las integran aplicando un método clásico de integración, tal como, Runge-Kutta. Otros utilizan programas comerciales de simulación, tales como el Bathfp, Simulink, etc. En cualquier caso, el objetivo de todos ellos es demostrar la viabilidad de determinados métodos y códigos para la modelización. En ninguno de ellos se profundiza en la evaluación de parámetros o variables de diseño.Acorde con todo lo anterior, los objetivos básicos de esta tesis se han centrado en la simulación mediante bond graph del comportamiento dinámico de las válvulas divisoras de caudal de doble funcionalidad con sistema antibloqueo al final de carrera de los actuadores, atendiendo al estudio de la influencia de las variables de diseño más importantes.Con esta metodología, se dispuso de una valiosa herramienta para simular, evaluar, predecir y plantear nuevos diseños de válvulas divisoras de caudal. Para la explotación sistemática de los modelos propuestos ha sido necesario desarrollar toda una serie de actividades experimentales con la finalidad de validarlos.Las conclusiones más destacadas son:La disposición de orificios en serie es la mejor configuración posible para el diseño de las válvulas divisoras/integradoras de caudal. El modelo de simulación desarrollado para la válvula divisora de caudal permite contemplar en un solo bond graph los dos modos de funcionamiento (unidad divisora/integradora), acoplable como submodelo a cualquier modelo jerárquicamente superior que represente un circuito real. Para valorar la precisión de las válvulas divisoras/integradoras de caudal se requiere considerar dos aspectos: el comportamiento dinámico del compensador de carga y el papel de las tolerancias de fabricación de las piezas con movimiento relativo, no solo en valor absoluto, sino también en valor relativo, es decir, nivel de simetría real del prototipo una vez fabricado.El dimensionado de los orificios principales fijos (chicles) definen el caudal nominal de la válvula, o rango útil de caudales. El nivel de precisión es directamente proporcional a la sección efectiva de estos chicles. Sin embargo reducir estas secciones impide mantener una aceptable caída de presión dentro del intervalo de caudales para el cual ha sido diseñada la válvula.El nivel de precisión es inversamente proporcional al cuadrado del diámetro de la corredera del compensador de carga, y se ve afectado ligeramente por el tamaño de los orificios de amortiguación.Para cargas simétricas, si el nivel de fugas es idéntico en ambas ramas, se produce un error de división de caudal nulo. Y el error absoluto entre desplazamientos de los cilindros, también, es nulo.Si los niveles de fugas son significativos y diferentes en ambas ramas, el error de división de caudal es apreciable y en consecuencia, el error absoluto entre desplazamientos de los actuadores aumenta con el tiempo.Para cargas asimétricas, si el nivel de fugas es idéntico en ambas ramas el error de división de caudal se compensa al cabo de un tiempo, y el error absoluto entre desplazamientos de los cilindros se mantiene constante (error caudal), una vez compensado el error de división de caudal. Sin embargo, la magnitud del error absoluto entre desplazamientos es proporcional al nivel de las fugas.Si el nivel de fugas es diferente en ambas ramas, el error de división de caudal tiende a compensarse, pero la influencia de la magnitud de las fugas en combinación con el estado de las cargas puede favorecer o desfavorecer el error absoluto de desplazamientos de los cilindros. / The development of this thesis is centered on the study of synchronization of hydraulic linear actuators by means of proportional flow divider-combiner valves.Most of the studies dedicated to the analysis of the dynamic behaviour of a flow divider valve raise the equations of state (system of nonlinear differential equations) and integrate them by applying a classic method of integration, such as, Runge-Kutta. Others use commercial simulation programs, such as the Bathfp, Simulink, etc. Even though the objective of all of them is to demonstrate the viability of certain methods and codes for their modelling, none of them deepen the evaluation of parameters or design variables.In agreement with the above, the basic objectives of this thesis have been centered on the simulation by means of the Bond graph method, of the dynamic behaviour of the flow divider-combiner valves of double functionality with system antiblockade at the end of the stroke of the actuators.With this methodology, it was a valuable tool to simulate, evaluate, predict and raise new designs for flow divider valves.For the systematic operation of the proposed models it has been necessary to develop a series of experimental activities with the purpose of validating them. The most outstanding conclusions are:The disposition of orifices in series is the best possible configuration for the design of the flow divider-combiner valves.The simulation model developed for the flow divider-combiner valve allows the analysis in a single Bond graph of both ways of operation (divider-combiner unit). It can be connected as a submodel to any hierarchically superior model that represents a real circuit.In order to value the precision of the flow divider-combiner valve it is necessary to consider two aspects: the dynamic behaviour of the load compensator and the influence of the manufacturing tolerances of the pieces with relative movement. However, both the single movements in absolute and relative valves must be considered.The sizing of the fixed primary holes defines the nominal flow rate of the valve, or the useful range of flows. The level of precision is directly proportional to the effective section of these holes. Nevertheless, reducing these sections prevents an acceptable drop in pressure within the interval of flow rate for which the valve has been designed.The level of precision is inversely proportional to the square of the diameter of the track of the load compensator, and it is affected slightly by the size of the damping orifices.For symmetrical loads, if the level of leakage is identical in both directions, a division error of zero flow takes place. The absolute error between displacements of the cylinders is also zero.If the leakage levels are significant and different in both directions, the flow error of the flow division is considerable and consequently, the absolute error between displacements of the actuators increases with time.For asymmetric loads, if the level of leakage is identical in both directions, the flow division error is compensated after a time and the absolute error between displacements of the cylinders stays constant (error of great flow), once flow division error is compensated. Nevertheless, the magnitude of the absolute error between displacements is proportional to the level of the leakage.If the level of leakage is different in both directions, the flow division error tends to be compensated, but the influence of the magnitude of the leakage in combination with the state of the loads can favour or work against the absolute error of displacements of the cylinders.

The numerical simulation model presented is based on the integration of the fluid conservation equations (continuity, momentum and energy) in the whole compressor domain (compression chamber, valves, manifolds, mufflers, connecting tubes, parallel paths, etc.) using instantaneous local mean values for the different variables. It is interesting to remark how momentum equation has been taken into account in all compressor parts and the possibility to solve parallel paths, resonators, etc. Effective flow areas are evaluated considering multidimensional models based on modal analysis of fluid interaction in the valve. Then, second and third order vibration models of valve are also considered. The possibility to use compound bound has been also implemented.The force balances in the crankshaft connecting rod mechanical system are simultaneously solved at each time-step considered in the thermal and fluid dynamic compressor model. It allows to evaluate the instantaneous compression chamber volume and the different forces in the crankshaft connecting rod mechanical system. Mechanical system forces allows to know important information to predict possibleover-stresses in piston, piston pin, crankshaft, connecting rod, etc.The thermal analysis of the solid elements is based on global energy balances at each macro volume considered (shell, muffler, tubes, cylinder head, crankcase, motor, etc.). Some improvements can be implemented (shell conduction, heat transfer coefficient evaluation, etc.).The resulting governing equations (fluid flow, valve dynamics, conduction heat transfer in solids, etc.) are discretized by means of a fully implicit control volume formulation. The complete set of algebraic equations are coupled using the segregated he complete set of algebraic equations are coupled using the segregated pressure based algorithm Semi-Implicit Method for Pressure-Linked Equations(SIMPLEC) extended to compressible flow. Second and third time order schemes have been implemented for the transient terms.An extensive hermetic reciprocating compressor experimental validation has been presented and the experimental know-how acquired has been highlighted. Furthermore, two commercial hermetic reciprocating compressor have been instrumented in detail to obtain the thermal temperatures map and the pressure fluid evolutions along compressor for different working conditions. It is interesting to remark as a novelty, the use of very small absolute pressure transducers, instead of the standard relative transducers. They allow to know instantaneous absolute pressure inside compressor chamber, without the necessity of measurement an absolute pressure outside the compression chamber (as is usual in this kind of experimental works). The global comparative results have allowed to check the possibilities of the numerical simulation presented above and its accuracy compared with experimental data. After that, this work show the capabilities offered by the simulation presented and its final objective, a better understanding of the thermal and fluid dynamic compressor behaviour to improve the design of these equipments.Then, the objective has been to review and present different physically meaningful parameters that characterize the reciprocating compressor behaviour (volumetric efficiency, isentropic efficiency, heat transfer efficiency, mechanical, electrical and heat losses, Coefficient of Performance, etc.), their influence detachment and evolution under different working conditions, with the idea to predict the performance of hermetic reciprocating compressors under different working conditions using the above mentioned non-dimensional parameters.Finally, a parametric study of hermetic reciprocating compressors behaviour has been carried out. Results presented show the influence of different aspects (geometry, valves, motor, working conditions, etc.) in the compressor behaviour. The parametric studies and compressor characterization detachment allows also a better implementation of simplest models of the compressors in the thermal and fluid dynamic numerical simulation of vapour compressor cycles together with the rest of elements.

La presente tesis doctoral trata sobre el accionamiento y control del posicionado de aperos agrícolas que configuran aspectos funcionales de máquinas tales como palas cargadoras frontales, máquinas cosechadoras de caña de azúcar (alzadoras y combinadas), etc., siendo el objetivo fundamental desarrollar un sistema que accione y controle dos o más actuadores óleo-hidráulicos de manera simultánea, coordinada y sincronizada que mejore la funcionalidad de elevación de mecanismos tipo péndulo doble.El estudio de la bibliografía consultada nos permite hacer un análisis detallado de los trabajos más importantes sobre los sistemas mecánicos de elevación, inclinación y volteo de tales máquinas, así como, las diferentes alternativas de sincronización.Acorde con lo anterior, los objetivos básicos de la tesis se han centrado en: estudiar el comportamiento de un bloque de sincronización oleohidráulica constituido por dos unidades divisoras de caudal, de relación de partición variable, con compensación de carga así como dos unidades de regulación de presión (frenado y limitación de presión) pilotadas externamente, simular su interacción con las prestaciones del mecanismo de elevación utilizando el método de bond graph, y por último, validar y optimizar el diseño y la metodología empleada a través de ensayos experimentales. La elección del bond graph como método de simulación obedece a dos objetivos básicos: evaluación del comportamiento cinemático y dinámico del mecanismo, así como permitir su interacción con el modelo del bloque de sincronización oleohidráulico. Es decir, disponer de un modelado físico en un sistema multidominio. La explotación sistemática de esta metodología teórico experimental ha permitido colegir las siguientes conclusiones:Desde el punto de vista mecánico/oleohidráulico, a pesar de que el comportamiento cinemático-dinámico del sistema durante los ciclos de subida y bajada son diferentes, es posible utilizar unidades divisoras de caudal de características inherentes idénticas (curva partición de caudal en función de la posición de la aguja). Los valores obtenidos en laboratorio y los obtenidos mediante la simulación numérica difieren en el peor de los casos, en un valor menor del 10 % (caso de bajada).Se ha constatado que, a excepción de que el mecanismo del péndulo doble incluidos los actuadores sea un mecanismo homotético, cualquier otra configuración geométrica del mismo, presenta una relación no lineal entre el cociente de velocidades de los actuadores en función del ángulo de elevación. Relación que limita drásticamente la sincronización. Las configuraciones no homotéticas son una constante en el accionamiento de aperos agrícolas y vienen impuestas, en la mayoría de las veces, por conveniencias estratégicas, tales como: ventajas mecánicas, visibilidad del operario, espacio de trabajo, estabilidad, transportabilidad, etc.Para el prototipo construido y ensayado, los valores de partición de caudal óptimos se corresponden a una posición de aguja de 1.2-1.6, para el ciclo de bajada y de 7-7.1, para el ciclo de subida. Sin embargo, no es posible compensar la no linealidad con válvulas de aguja controladas manualmente.También se ha constatado, experimentalmente y numéricamente, que la bondad del sincronizado es muy sensible a la posición de la aguja que define la partición de caudal.En general, se puede considerar que la solución oleohidráulica ofrece unas buenas prestaciones aunque se ven penalizadas ya que la unidad hidráulica de sincronizado acusa una pérdida de presión excesiva.Habida cuenta la gran sensibilidad de la posición de la aguja frente a la división de caudal mencionada anteriormente, se ha constatado, experimentalmente, una acusada diferencia de comportamiento entre los ciclos de subida y bajada. En algunos casos, estas diferencias han sido acumulativas y, en consecuencia, transcurrido un elevado número de ciclos, se pierde la capacidad de sincronización. / This doctoral thesis deals with the drive and control of the positioning of agricultural implements that form functional aspects of machines such as front-end loaders, sugar cane machines harvesters (gatherers and combined), etc. The main aim is to develop a system that drives and controls two or more hydraulic actuators in a simultaneous, coordinated and synchronised way that improves the elevation functionality of double pendulum type mechanisms.The study of the consulted bibliography allows us to make a detailed analysis of the most important works of the elevation, inclination and turn around mechanical systems of such machines, as well as, the different synchronization alternatives.In agreement with the above, the basic objectives of this thesis have been centered on: the study of the behavior of a hydraulic synchronization block consisting of two variable partition flow divider units with load compensation as well as two externally piloted pressure regulation units (restrained and limitation of pressure); the simulation of their interaction with the elevating mechanism characteristics using the bond graph method, and finally, the validation and optimization the design and the methodology used through experimental tests. The election of the bond graph as a simulation method obeys the two basic objectives: evaluation of the kinematic and dynamic behavior of the mechanism, as well as its interaction with the model of the hydraulic synchronization block. That is to say, to have a physical model in a multidominion system. The systematic operation of this experimental theoretical methodology has lead to the following conclusions:From the mechanical / hydraulic point of view, although the kinematic-dynamic behavior of the system during the cycles of ascent and descent is different, it is possible to use flow divider units of identical inherent characteristics (the flow partition curve as a function of the position of the needle). The values obtained in the laboratory and these obtained by means of the numerical simulation differ in the worst case by a value smaller than 10 % (case of descent).It has been stated that all geometric configurations of the double pendulum mechanisms including actuators (with the exception of a homothetic type) present a nonlinear relationship between the quotients of speeds of the actuators based on the elevation angle. A relationship that limits the synchronization drastically. The nonhomothetics configurations are constant in agricultural drive implements and are imposed, most of the time, by strategic advantages, such as: mechanical advantages, visibility of the worker, space of work, stability, transportability, etc.For the manufactured and tested prototype, the optimal partition of flow values correspond to a position of the needle of 1.2-1.6, for the cycle of descent and 7-7.1, for the cycle of ascent. Nevertheless, it is not possible to manually compensate the nonlinearity with needle valves.It has also been stated, experimentally and numerically, that the synchronous is very sensitive to the position of the needle that defines the flow partition.In general, it is possible to consider that the hydraulic solution offers good performance characteristics although they are penalized as the synchronised hydraulic unit causes an excessive loss of pressure.Taking into account the great sensitivity of the position of the needle as opposed to the flow division previously mentioned, it has been stated experimentally that there is a significant behavior difference between the cycles of ascent and descent. In some cases, these differences have been cumulative and consequently, after a high number of cycles, the synchronization capacity is lost.