Articles de revista
http://hdl.handle.net/2117/3481
20180317T12:30:46Z

Geosynchronous inclined orbits for highlatitude communications
http://hdl.handle.net/2117/109270
Geosynchronous inclined orbits for highlatitude communications
Fantino, Elena; Flores Le Roux, Roberto Maurice; Di Carlo, Marilena; Di Salvo, Alessio; Cabot Soler, Elisenda
We present and discuss a solution to the growing demand for satellite telecommunication coverage in the highlatitude geographical regions (beyond 55°N), where the signal from geostationary satellites is limited or unavailable. We focus on the dynamical issues associated to the design, the coverage, the maintenance and the disposal of a set of orbits selected for the purpose. Specifically, we identify a group of highly inclined, moderately eccentric geosynchronous orbits derived from the Tundra orbit (geosynchronous, eccentric and critically inclined). Continuous coverage can be guaranteed by a constellation of three satellites in equally spaced planes and suitably phased. By means of a highprecision model of the terrestrial gravity field and the relevant environmental perturbations, we study the evolution of these orbits. The effects of the different perturbations on the ground track (which is more important for coverage than the orbital elements themselves) are isolated and analyzed. The physical model and the numerical setup are optimized with respect to computing time and accuracy. We show that, in order to maintain the ground track unchanged, the key parameters are the orbital period and the argument of perigee. Furthermore, corrections to the right ascension of the ascending node are needed in order to preserve the relative orientation of the orbital planes. A stationkeeping strategy that minimizes propellant consumption is then devised, and comparisons are made between the cost of a solution based on impulsive maneuvers and one with continuous thrust. Finally, the issue of endoflife disposal is discussed.
© 2017. This version is made available under the CCBYNCND 4.0 license http://creativecommons.org/licenses/byncnd/4.0/
20171026T14:15:28Z
Fantino, Elena
Flores Le Roux, Roberto Maurice
Di Carlo, Marilena
Di Salvo, Alessio
Cabot Soler, Elisenda
We present and discuss a solution to the growing demand for satellite telecommunication coverage in the highlatitude geographical regions (beyond 55°N), where the signal from geostationary satellites is limited or unavailable. We focus on the dynamical issues associated to the design, the coverage, the maintenance and the disposal of a set of orbits selected for the purpose. Specifically, we identify a group of highly inclined, moderately eccentric geosynchronous orbits derived from the Tundra orbit (geosynchronous, eccentric and critically inclined). Continuous coverage can be guaranteed by a constellation of three satellites in equally spaced planes and suitably phased. By means of a highprecision model of the terrestrial gravity field and the relevant environmental perturbations, we study the evolution of these orbits. The effects of the different perturbations on the ground track (which is more important for coverage than the orbital elements themselves) are isolated and analyzed. The physical model and the numerical setup are optimized with respect to computing time and accuracy. We show that, in order to maintain the ground track unchanged, the key parameters are the orbital period and the argument of perigee. Furthermore, corrections to the right ascension of the ascending node are needed in order to preserve the relative orientation of the orbital planes. A stationkeeping strategy that minimizes propellant consumption is then devised, and comparisons are made between the cost of a solution based on impulsive maneuvers and one with continuous thrust. Finally, the issue of endoflife disposal is discussed.

Nonlinear loads model for harmonics flow prediction, using multivariate regression
http://hdl.handle.net/2117/108748
Nonlinear loads model for harmonics flow prediction, using multivariate regression
Lamich Arocas, Manuel; Balcells Sendra, Josep; Corbalán Fuertes, Montserrat; Griful Ponsati, Eulàlia
This paper describes a method for obtaining a model of a single or a set of nonlinear loads (NLL) connected to a certain point of an electrical network. The basic assumption is that the network supplying the NLL has significant series impedances and is disturbed by other parallel, random, and unknown neighbor loads, sharing part of the supply system with the NLL. The main interest for obtaining the model is its further use to predict the amount and flow of harmonic currents generated by the NLL, in the case of adding a filter to reduce the harmonics distortion. The modeling technique used in the paper is based on multivariate multiple outputs regression and leads to a set of equations giving the NLL behavior (one for each of the harmonic currents). The model is obtained from data taken at measuring point and is only valid to predict the NLL behavior when new loads are connected at this point. The modeling method was first tested with V, I data coming from simulations using a MATLABSimulink SimPowerSystems toolbox. Finally, the method has been validated using V, I data taken in a real installation with different neighbor loads and under different load conditions.
20171017T12:00:15Z
Lamich Arocas, Manuel
Balcells Sendra, Josep
Corbalán Fuertes, Montserrat
Griful Ponsati, Eulàlia
This paper describes a method for obtaining a model of a single or a set of nonlinear loads (NLL) connected to a certain point of an electrical network. The basic assumption is that the network supplying the NLL has significant series impedances and is disturbed by other parallel, random, and unknown neighbor loads, sharing part of the supply system with the NLL. The main interest for obtaining the model is its further use to predict the amount and flow of harmonic currents generated by the NLL, in the case of adding a filter to reduce the harmonics distortion. The modeling technique used in the paper is based on multivariate multiple outputs regression and leads to a set of equations giving the NLL behavior (one for each of the harmonic currents). The model is obtained from data taken at measuring point and is only valid to predict the NLL behavior when new loads are connected at this point. The modeling method was first tested with V, I data coming from simulations using a MATLABSimulink SimPowerSystems toolbox. Finally, the method has been validated using V, I data taken in a real installation with different neighbor loads and under different load conditions.

Influence of data resolution in nonlinear loads model for harmonics prediction
http://hdl.handle.net/2117/106333
Influence of data resolution in nonlinear loads model for harmonics prediction
Balcells Sendra, Josep; Lamich Arocas, Manuel; Griful Ponsati, Eulàlia; Corbalán Fuertes, Montserrat
This paper describes the influence of data resolution in the agreement of models to predict harmonics generated by nonlinear loads (NLL), basically formed by single phase and three phase rectifiers, eventually combined with linear loads. We assume that the network supplying the NLL has significant impedances and that it is disturbed by other parallel, random and unknown neighbor loads, sharing part of the supply system. The aim of building NLL models is to make predictions on the amount and flow paths of harmonic currents generated by such NLL in case of using parallel filters. In this paper, the models are obtained from sets of (V,I) data taken at a certain point, called measuring point (MP) and are valid to predict the NLL behavior when random known or unknown parallel loads are connected upstream of this point. The technique used to obtain the models studied here is based on Multivariate Multiple Outputs Regression (MMOR) and will not be described in detail in this paper. This method allows obtaining a set of equations giving the current harmonics as a function of voltage harmonics observed at the measuring point (MP). The accordance between model and the experimental results is very dependent on the resolution and accuracy of V and I measurements at the MP and is the core matter of this paper.
20170711T07:59:53Z
Balcells Sendra, Josep
Lamich Arocas, Manuel
Griful Ponsati, Eulàlia
Corbalán Fuertes, Montserrat
This paper describes the influence of data resolution in the agreement of models to predict harmonics generated by nonlinear loads (NLL), basically formed by single phase and three phase rectifiers, eventually combined with linear loads. We assume that the network supplying the NLL has significant impedances and that it is disturbed by other parallel, random and unknown neighbor loads, sharing part of the supply system. The aim of building NLL models is to make predictions on the amount and flow paths of harmonic currents generated by such NLL in case of using parallel filters. In this paper, the models are obtained from sets of (V,I) data taken at a certain point, called measuring point (MP) and are valid to predict the NLL behavior when random known or unknown parallel loads are connected upstream of this point. The technique used to obtain the models studied here is based on Multivariate Multiple Outputs Regression (MMOR) and will not be described in detail in this paper. This method allows obtaining a set of equations giving the current harmonics as a function of voltage harmonics observed at the measuring point (MP). The accordance between model and the experimental results is very dependent on the resolution and accuracy of V and I measurements at the MP and is the core matter of this paper.

Aposteriori error estimation for the finite point method with applications to compressible flow
http://hdl.handle.net/2117/104305
Aposteriori error estimation for the finite point method with applications to compressible flow
Ortega, Enrique; Flores Le Roux, Roberto Maurice; Oñate Ibáñez de Navarra, Eugenio; Idelsohn Barg, Sergio Rodolfo
An aposteriori error estimate with application to inviscid compressible flow problems is presented. The estimate is a surrogate measure of the discretization error, obtained from an approximation to the truncation terms of the governing equations. This approximation is calculated from the discrete nodal differential residuals using a reconstructed solution field on a modified stencil of points. Both the error estimation methodology and the flow solution scheme are implemented using the Finite Point Method, a meshless technique enabling higherorder approximations and reconstruction procedures on general unstructured discretizations. The performance of the proposed error indicator is studied and applications to adaptive grid refinement are presented.
20170511T08:48:09Z
Ortega, Enrique
Flores Le Roux, Roberto Maurice
Oñate Ibáñez de Navarra, Eugenio
Idelsohn Barg, Sergio Rodolfo
An aposteriori error estimate with application to inviscid compressible flow problems is presented. The estimate is a surrogate measure of the discretization error, obtained from an approximation to the truncation terms of the governing equations. This approximation is calculated from the discrete nodal differential residuals using a reconstructed solution field on a modified stencil of points. Both the error estimation methodology and the flow solution scheme are implemented using the Finite Point Method, a meshless technique enabling higherorder approximations and reconstruction procedures on general unstructured discretizations. The performance of the proposed error indicator is studied and applications to adaptive grid refinement are presented.

In situ particle zeta potential evaluation in electroosmotic flows from timeresolved microPIV measurements
http://hdl.handle.net/2117/103563
In situ particle zeta potential evaluation in electroosmotic flows from timeresolved microPIV measurements
Sureda Anfres, Miquel; Miller, Andrew; Diez, F. J.
A timeresolved microPIV method is presented to measure in an EOF the particles zeta potential in situ during the transient startup of a microdevice. The method resolves the electrophoretic velocity of fluorospheres used as tracer particles in microPIV. This approach exploits the short transient regime of the EOF generated after a potential drop is imposed across a microchannel and before reaching quasisteady state. During the starting of the transient regime, the electrophoretic effect is dominant in the center of the channel and the EOF is negligible. By measuring the velocity of the tracer particles with a microPIV system during that starting period, their electrophoretic velocity is obtained. The technique also resolves the temporal evolution of the EOF with three regions identified. The first region occurs before the electroosmotic effect reaches the center of the channel, the second region extends until the EOF reaches steady state, and thereafter is the third region. The two time constants separating these regions are also obtained and compared to the theory. The zeta potential of 860 nm diameter polystyrene particles is calculated for different solutions including borate buffer, sodium chloride, and deionized water. Results show that the magnitudes of the electrophoretic and electroosmotic velocities are in the range of 300 to 700 µm/s for these measurements. The zeta potential values are compared to the wellestablished closed cell technique showing improved accuracy. The method also resolves the characteristic response time of the EOF, showing small but important deviations from current analytical predictions. Additionally, the measurements can be performed in situ in microfluidic devices under actual working EOF conditions and without the need for calibrations.
20170420T08:07:22Z
Sureda Anfres, Miquel
Miller, Andrew
Diez, F. J.
A timeresolved microPIV method is presented to measure in an EOF the particles zeta potential in situ during the transient startup of a microdevice. The method resolves the electrophoretic velocity of fluorospheres used as tracer particles in microPIV. This approach exploits the short transient regime of the EOF generated after a potential drop is imposed across a microchannel and before reaching quasisteady state. During the starting of the transient regime, the electrophoretic effect is dominant in the center of the channel and the EOF is negligible. By measuring the velocity of the tracer particles with a microPIV system during that starting period, their electrophoretic velocity is obtained. The technique also resolves the temporal evolution of the EOF with three regions identified. The first region occurs before the electroosmotic effect reaches the center of the channel, the second region extends until the EOF reaches steady state, and thereafter is the third region. The two time constants separating these regions are also obtained and compared to the theory. The zeta potential of 860 nm diameter polystyrene particles is calculated for different solutions including borate buffer, sodium chloride, and deionized water. Results show that the magnitudes of the electrophoretic and electroosmotic velocities are in the range of 300 to 700 µm/s for these measurements. The zeta potential values are compared to the wellestablished closed cell technique showing improved accuracy. The method also resolves the characteristic response time of the EOF, showing small but important deviations from current analytical predictions. Additionally, the measurements can be performed in situ in microfluidic devices under actual working EOF conditions and without the need for calibrations.

Ramair parachute simulation with panel methods and staggered coupling
http://hdl.handle.net/2117/91426
Ramair parachute simulation with panel methods and staggered coupling
Ortega, Enrique; Flores Le Roux, Roberto Maurice; Pons Prats, Jordi
20161103T13:43:42Z
Ortega, Enrique
Flores Le Roux, Roberto Maurice
Pons Prats, Jordi

Reliability versus mass optimization of CO2 extraction technologies for long duration missions
http://hdl.handle.net/2117/89135
Reliability versus mass optimization of CO2 extraction technologies for long duration missions
Detrell Domingo, Gisela; Griful Ponsati, Eulàlia; Messerschmid, Ernst
The aim of this paper is to optimize reliability and mass of three CO2 extraction technologies/components: the 4Bed Molecular Sieve, the Electrochemical Depolarized Concentrator and the Solid Amine Water Desorption. The first one is currently used in the International Space Station and the last two are being developed, and could be used for future long duration missions. This work is part of a complex study of the Environmental Control and Life Support System (ECLSS) reliability. The result of this paper is a methodology to analyze the reliability and mass at a component level, which is used in this paper for the CO2 extraction technologies, but that can be applied to the ECLSS technologies that perform other tasks, such as oxygen generation or water recycling, which will be a required input for the analysis of an entire ECLSS. The key parameter to evaluate any system to be used in space is mass, as it is directly related to the launch cost. Moreover, for long duration missions, reliability will play an even more important role, as no resupply or rescue mission is taken into consideration. Each technology is studied as a reparable system, where the number of spare parts to be taken for a specific mission will need to be selected, to maximize the reliability and minimize the mass of the system. The problem faced is a MultiObjective Optimization Problem (MOOP), which does not have a single solution. Thus, optimum solutions of MOOP, the ones that cannot be improved in one of the two objectives, without degrading the other one, are found for each selected technology. The solutions of the MOOP for the three technologies are analyzed and compared, considering other parameters such as the type of mission, the maturity of the technology and potential interactions/synergies with other technologies of the ECLSS.
20160725T10:43:03Z
Detrell Domingo, Gisela
Griful Ponsati, Eulàlia
Messerschmid, Ernst
The aim of this paper is to optimize reliability and mass of three CO2 extraction technologies/components: the 4Bed Molecular Sieve, the Electrochemical Depolarized Concentrator and the Solid Amine Water Desorption. The first one is currently used in the International Space Station and the last two are being developed, and could be used for future long duration missions. This work is part of a complex study of the Environmental Control and Life Support System (ECLSS) reliability. The result of this paper is a methodology to analyze the reliability and mass at a component level, which is used in this paper for the CO2 extraction technologies, but that can be applied to the ECLSS technologies that perform other tasks, such as oxygen generation or water recycling, which will be a required input for the analysis of an entire ECLSS. The key parameter to evaluate any system to be used in space is mass, as it is directly related to the launch cost. Moreover, for long duration missions, reliability will play an even more important role, as no resupply or rescue mission is taken into consideration. Each technology is studied as a reparable system, where the number of spare parts to be taken for a specific mission will need to be selected, to maximize the reliability and minimize the mass of the system. The problem faced is a MultiObjective Optimization Problem (MOOP), which does not have a single solution. Thus, optimum solutions of MOOP, the ones that cannot be improved in one of the two objectives, without degrading the other one, are found for each selected technology. The solutions of the MOOP for the three technologies are analyzed and compared, considering other parameters such as the type of mission, the maturity of the technology and potential interactions/synergies with other technologies of the ECLSS.

A meshless finite point method for threedimensional analysis of compressible flow problems involving moving boundaries and adaptivity
http://hdl.handle.net/2117/86276
A meshless finite point method for threedimensional analysis of compressible flow problems involving moving boundaries and adaptivity
Ortega, Enrique; Oñate Ibáñez de Navarra, Eugenio; Idelsohn Barg, Sergio Rodolfo; Flores Le Roux, Roberto Maurice
A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwindbiased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dualtime steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and hadaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multicore performance of the proposed technique is also discussed through the examples provided.
This article may be used for noncommercial purposes in accordance with Wiley Terms and Conditions for SelfArchiving.
20160427T14:59:20Z
Ortega, Enrique
Oñate Ibáñez de Navarra, Eugenio
Idelsohn Barg, Sergio Rodolfo
Flores Le Roux, Roberto Maurice
A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwindbiased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dualtime steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and hadaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multicore performance of the proposed technique is also discussed through the examples provided.

Threedimensional numerical simulation of an external gear pump with decompression slot and meshing contact point
http://hdl.handle.net/2117/85440
Threedimensional numerical simulation of an external gear pump with decompression slot and meshing contact point
Castilla López, Roberto; Gámez Montero, Pedro Javier; Campo Sud, David del; Raush Alviach, Gustavo Adolfo; García Vilchez, Mercedes; Codina Macià, Esteban
Recently several works have been published on numerical simulation of an external gear pump (EGP). Such kinds of pumps are simple and relatively inexpensive, and are frequently used in fluid power applications, such as fluid power in aeronautical, mechanical, and civil engineering. Nevertheless, considerable effort is being undertaken to improve efficiency and reduce noise and vibration produced by the flow and pressure pulsations. Numerical simulation of an EGP is not straightforward principally for two main reasons. First, the gearing mechanism between gears makes it difficult to handle a dynamic mesh without a considerable deterioration of mesh quality. Second, the dynamic metalmetal contact simulation is important when high pressure outflow has to be reproduced. The numerical studies published so far are based on a twodimensional (2D) approximation. The aim of the present work is to contribute to the understanding of the fluid flow inside an EGP by means of a complete threedimensional (3D) parallel simulation on a cluster. The 3D flow is simulated in a LINUX cluster with a solver developed with the OPENFOAM Toolbox. The hexahedral mesh quality is maintained by periodically replacing the mesh and interpolating the physical magnitudes fields. The meshing contact point is simulated with the viscous wall approach, using a viscosity model based on wall proximity. The results for the flow rate ripples show a similar behavior to that obtained with 2D simulations. However, the flow presents important differences inside the suction and the discharge chambers, principally in the regions of the pipes' connection. Moreover, the decompression slot below the gearing zone, which can not be simulated with a 2D approximation, enables a more realistic simulation of a contact ratio greater than 1. The results are compared with experimental measurements recently published.
20160408T14:13:47Z
Castilla López, Roberto
Gámez Montero, Pedro Javier
Campo Sud, David del
Raush Alviach, Gustavo Adolfo
García Vilchez, Mercedes
Codina Macià, Esteban
Recently several works have been published on numerical simulation of an external gear pump (EGP). Such kinds of pumps are simple and relatively inexpensive, and are frequently used in fluid power applications, such as fluid power in aeronautical, mechanical, and civil engineering. Nevertheless, considerable effort is being undertaken to improve efficiency and reduce noise and vibration produced by the flow and pressure pulsations. Numerical simulation of an EGP is not straightforward principally for two main reasons. First, the gearing mechanism between gears makes it difficult to handle a dynamic mesh without a considerable deterioration of mesh quality. Second, the dynamic metalmetal contact simulation is important when high pressure outflow has to be reproduced. The numerical studies published so far are based on a twodimensional (2D) approximation. The aim of the present work is to contribute to the understanding of the fluid flow inside an EGP by means of a complete threedimensional (3D) parallel simulation on a cluster. The 3D flow is simulated in a LINUX cluster with a solver developed with the OPENFOAM Toolbox. The hexahedral mesh quality is maintained by periodically replacing the mesh and interpolating the physical magnitudes fields. The meshing contact point is simulated with the viscous wall approach, using a viscosity model based on wall proximity. The results for the flow rate ripples show a similar behavior to that obtained with 2D simulations. However, the flow presents important differences inside the suction and the discharge chambers, principally in the regions of the pipes' connection. Moreover, the decompression slot below the gearing zone, which can not be simulated with a 2D approximation, enables a more realistic simulation of a contact ratio greater than 1. The results are compared with experimental measurements recently published.

A meshless finite point method for three dimensional analysis of compressible flow problems involving moving boundaries and adaptivity
http://hdl.handle.net/2117/28569
A meshless finite point method for three dimensional analysis of compressible flow problems involving moving boundaries and adaptivity
Ortega, Enrique; Oñate Ibáñez de Navarra, Eugenio; Idelsohn Barg, Sergio Rodolfo; Flores Le Roux, Roberto Maurice
A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwindbiased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dualtime steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and hadaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multicore performance of the proposed technique is also discussed through the examples provided.
20150713T10:44:47Z
Ortega, Enrique
Oñate Ibáñez de Navarra, Eugenio
Idelsohn Barg, Sergio Rodolfo
Flores Le Roux, Roberto Maurice
A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwindbiased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dualtime steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and hadaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multicore performance of the proposed technique is also discussed through the examples provided.