Jean Marie Gilliot

Serge Garlatti

An Adaptive and Context-Aware Architecture for Future Pervasive Learning Environments

Keywords: Pervasive learning, context-aware, adaptation, semantic web, linked data, ontology

I. Introduction

Nowadays, technology-enhanced learning systems can be build on learning management systems (LMS) or personal learning environments (PLE). On one hand, LMS provide administrative services (admissions, registration, etc.), trace tools of learners activities accessible for teachers, and serveral tools (discussion forums, file sharing, management of assignements, learning activities, chat, etc.). The main drawbacks are that those tools are pre defined and the correspond pedagogical approaches are implicitly frozen, thus choice is limited, and that they are bounded to a local area of participants.

On the other hand, PLE are distributed, personnalized and composed of separated tools, social media applications. From an educational perspective, social media applications including blogs, wikis, rich media sharing, etc. fit well with socio-constructivist learning approaches as they provide spaces for collaborative knowledge building and reflective practices. These social media tools are used in informal learning settings commonly found outside formal and institutional learning environments. In PLE, learning objects are distributed on the web in different tools that do not provide interoperability; activities are implicit and difficult to monitor; learners traces are almost impossible to get. Those environments are at present exclusive [START_REF] Dalsgaard | Social software: E-learning beyond learning management systems[END_REF].

Moreover, pervasive learning is recently becoming an important issue in technology-enhanced learning field. Pervasive computing has the ability to increase our capability to physically move computing tools and services with us and to inquire, detect and explore the surrounding environment in order to obtain information and to dynamically build context models that allow for supporting different aspects of the learning process [START_REF] Lyytinen | Issues and challenges in ubiquitous computing[END_REF]. For instance, computers can obtain information about the context of learning from the learning environment where small devices, sensors, pads, badges, learners, teachers or tutors, communities and so on, are embedded and communicate mutually. The physical environment is directly related to learning goals and activities. New situated learning activities may be achieved to enhance the learning process that were difficult to realized before (even impossible before). Context based activities, or pervasive learning, imply the ability to provide the right content, with the right tool within a well suited activity. Social interaction in such context make necessary the interoperablity of data across tools. Assessment of activities (either formative or summative) call for learning trace access. All of this also means that pedagogical activities should be adapted on the fly in an open, distributed environment, requires the definition of models, based on knowledge models and pedagogical theories, used by a distributed adaptation framework.

A major challenge for future technology-enhanced learning systems is to combine the benefits of learning management systems, personal learning environments, social media software and pervasive learning environments. According to recent studies, technology-enhanced learning (TEL) systems must have the capability to reuse learning resources and web services from large repositories, to take into account the context and to allow dynamic adaptation to different learners based on substantial advances in pedagogical theories and knowledge models [START_REF] Balacheff | 10 issues to think about the future of research on TEL[END_REF]. It is particularly true in pervasive learning. The reuse of learning resources and/or data requires interoperability at semantic level. In other words, it is necessary to have a semantic web approach to design TEL systems. Moreover, knowledge models and pedagogical theories can be fully represented by means of a semantic web approach.

Unfortunately, learning management systems and social media applications are data silos. In other words, data are unavailable on the web. Only people may have access to data, not computers. Reuse and exchange of data among LMS and social tools are only possible by means of APIthat is to say manually by mean of one API per tool. On the contrary, semantic web provides a common framework that allows data, information and knowledge to be shared and reused across applications, enterprises, and community boundaries. In such a framework, linked data describes a method of exposing, sharing, and connecting data, information and knowledge on the Web [START_REF] Bojaars | Using the Semantic Web for linking and reusing data across Web 2.0 communities[END_REF][START_REF] Gruber | Collective knowledge systems: Where the Social Web meets the Semantic Web[END_REF]. It provides a standardized, uniform and generic method for data discovery, distributed queries against several data repositories, integration or semantic mash-up, uniform access to metadata, data, information and knowledge. Some metadata can be generated automatically (sometimes on the fly) from the tool databases according to common vocabularies like Dublin Core, SKOS, SIOC, FOAF, etc. Most of these vocabularies are lightweight ontologies that can fit well database schemas. These vocabularies provide common semantic enabling computers to put queries on LMS and social media tools. Thus, the web can be viewed as a single global database. Users and/or computers can perform complex queries against this global database using the SPARQL language. Complex queries are queries over multiples pages / web sites / data repositories whatever the tool is. It only has to expose data on the common standard and vocabularies. Thus, future pervasive learning environments can be composed of lots of different LMS, social media tools exposing, sharing and connecting data, information and knowledge on the web.

Our vision is that future pervasive learning environments will be based on contextual adaptation of pedagogical activities and resources based on a semantic web approach to provide on the fly the distributed software environment, composed of appropriate standards tools, to enable the fullfillment of proposed or chosen assessable activities in a social environment.

II. Problems to solve

First of all, we highlight some functionalities that could be available by means of linked data. A tutor and/or a learner may gather the outcomes of learners' activities across multiple sites or relevant resources and tools may be recommended to a learner taking into account his social network and interests and/or his learning communities as found in social and learning web applications. Information available in a learner social networks and learning communities could be used to enrich his context features. Ontologies, like SIOC and FOAF, enable us to provide these type of functionalities. The related metadata are computed on the fly, for instance by specific plugins in social web applications. These ontologies fit well the corresponding attributes of social web databases.

Nevertheless, these ontologies are not sufficient to achieve all functionalities required to build future pervasive learning environments. Fulfilling this vision requires to solve a number of technical and pedagogical issues. Some of these issues are as follows:

1. Learning resources and web services retrieval 2. Pedagogical issues: activity modeling and management

Context modelling, management and adaptation policies

Learning resources and web services retrieval

At present, learning ressources are not indexed by means of semantic metatdata. To retrieve relevant resources, it is necessary to have at least a domain ontology modeling the topics of interest related to the learning goals and, of course, to setup corresponding metadata. In most of the cases, it is not possible to do it automatically. Specifc tools has to be developed to suggest relevant metadata to users.

Pedagogical issues: activity modeling and management

Activities have to be structured, coordinated and explictly represented in order to manage them and to define their relevant contexts.provide At least, two complementary approaches could be used: i) Recommendation of a priori defined and well structured activities that may be achieved in different ways according to the current situation; ii) Recommendation of activities and learning ressources according to the current situation. Consequently, ontologies must be designed to represent and manage learning activties. In future pervasive learning environments, an activity may be achieved in different tools and a tool can be used to accomplish different activties. It is necessary to link activities and tools to be able to monitor and to support learners according to context and learning goals.

Enrichment and management of context

Dourish [START_REF] Dourish | What we talk about when we talk about context[END_REF]) claims that context arises from the activities and cannot be separated from them. It is actively produced, maintained and enacted during activities. We claim that activities embedded in a particular physical world (or environment) are key issues to give us intention and meaning according to different situations and finally to determine the relevant features describing the different situations. Many definitions of context are given in literature. We mainly focus on some of them: "learning context is used to describe the current situation of a person related to a learning activity; in addition to attributes relying on the physical world model [START_REF] Derntl | Modeling Context-Aware E-learning Scenarios[END_REF]; " information and content in use to support a specific activity (being individual or collaborative) in a particular physical environment" [START_REF] Kurti | Designing and implementing ubuquitous learning activities supported by mobile and positioning technologies[END_REF]).

An history of learners and/or tutors activities across social web applications can be retrieved and/or discovered, as all user activities are time recorded. Such history could be used to manage seamless learning and learning across context issues. In the mid-time, the Semantic Sensor Network Incubator Group would provide ontologies that define the capabilities of sensors and sensor networks and develop semantic annotations of a key language used by services based sensor networks (http://www.w3.org/2005/Incubator/ssn/). Acquire data from network sensors at knowledge level will also enrich the learner context features.

III. Conclusion

Our vision is that future pervasive learning environments will be based on contextual adaptation of pedagogical activities and resources based on a semantic web approach to provide on the fly the distributed software environment, composed of appropriate standards tools, to enable the fullfillment of proposed or chosen assessable activities in a social environment.

In order to achieve this new generation of learning environments, appropriate ontologies for semantic web should be defined to enable development of new pervasive pedagogical activities and usage in an open environment shared between new generation of learners, teachers and learning communities.


Luca Pagliarini
email: luca.pagliarini@unipr.it

Luca Cattani
email: luca.cattani1@unipr.it

Vincent Ayel
email: vincent.ayel@ensma.fr

Maksym Slobodeniuk
email: maksym.slobodeniuk@ensma.fr

Cyril Romestant
email: cyril.romestant@ensma.fr

Fabio Bozzoli
email: fabio.bozzoli@unipr.it

Thermographic Investigation on Fluid Oscillations and Transverse Interactions in a Fully Metallic Flat-Plate Pulsating Heat Pipe

Keywords: flat-plate pulsating heat pipe, infrared thermographic acquisitions, wavelet method, fluid oscillations, transverse conduction

The present investigation deals with the quantification of fluid oscillation frequencies in a metallic pulsating heat pipe tested at varying heat loads and orientations. The aim is to design a robust technique for the study of the inner fluid dynamics without adopting typical experimental solutions, such as direct fluid visualizations through transparent inserts. The studied device is made of copper, and it is partially filled with a water-ethanol mixture (20 wt.% of ethanol). Heat fluxes locally exchanged between the working fluid and the device walls are first assessed through the inverse heat conduction problem resolution approach by processing outer wall temperature distributions acquired by thermography. The estimated local heat transfer quantities are therefore processed to quantify the fluid oscillatory behavior in every device branch during the intermittent flow and full activation regimes, thus providing a deeper insight into the heat transfer modes. After dealing with a further validation of the inverse approach in terms of oscillation frequency restoration capability, the wall-to-fluid heat fluxes referred to each channel are processed by means of the wavelet method. Scalograms and power spectra of the considered signals are presented for a time-based analysis of the working fluid oscillations, as well as for the identification of dominant oscillation frequencies. Fluid motion is then quantified in terms of the continuity of fluid oscillations and activity of channels by applying a scalogram denoising technique named K-means clustering method. Moreover, a statistical reduction of the channel-wise dominant oscillation frequencies is performed to provide useful references for the interpretation of the overall oscillatory behavior. The link between oscillations and transverse interactions is finally investigated. The vertical bottomheated mode exhibits stronger fluid oscillations with respect to the horizontal mode, with fluid oscillation frequencies ranging from 0.78 up to 1 Hz. Nonetheless, the fluid motion is more stable in terms of oscillation frequency between channels when the device operates in the horizontal orientation probably due to negligible buoyancy effects. Moreover, thermal interactions between adjacent channels are found to be stronger when the oscillatory behavior presents similar features from channel to channel in horizontal orientation. The proposed method for fluid oscillation analyses in fully metallic flat-plate pulsating heat pipes can be effectively adopted to other flat-plate layouts without any need for transparent windows, thus reducing the overall complexity of experimental set-ups and providing, at the same time, a good insight into the inner fluid dynamics.

Introduction

The self-sustained motion of the working fluid greatly affects the overall thermal behavior of passive two-phase heat transfer devices, such as heat pipes (HPs) and pulsating heat pipes (PHPs). In fact, in PHPs, fluid oscillations are responsible for efficient heat dissipation at the heated section (evaporator) [START_REF] Spinato | Understanding the self-sustained oscillating two-phase flow motion in a closed loop pulsating heat pipe[END_REF]. Specifically, in such passive systems, the fluid flow is promoted by temperature differences between the evaporator and the cooled section (condenser) rather than external auxiliary equipment (such as pumps). Since the fluid is mainly saturated, thermal instabilities result in pressure differences which act as driving forces for its self-motion [START_REF] Marengo | Pulsating Heat Pipes: Experimental Analysis[END_REF]. It is thus evident that such devices are governed by an extremely complex interplay between thermodynamic and fluid-dynamic behavior. The employment of HPs and PHPs has been found to be very attractive for many thermal management applications requiring a drastic reduction in electrical energy consumption [START_REF] Murshed | A critical review of traditional and emerging techniques and fluids for electronics cooling[END_REF]. Nevertheless, one of the main drawbacks of passive devices is represented by the need for a constant and self-induced backflow of liquid from the condenser to the evaporator, i.e., a continuously supplied evaporation process. While in HPs a persistent liquid flow to the evaporator is ensured by properly designed wick structures, such a phenomenon is only supported by non-equilibrium conditions in PHPs, thus largely increasing the complexity of their working behavior [START_REF] Nazari | A review on pulsating heat pipes: From solar to cryogenic applications[END_REF]. For this reason, the thermo-fluid dynamics underlying the oscillatory modes in PHPs are not yet completely understood despite the intensive experimental efforts of the last decades. In general, fluid oscillations are described in terms of amplitude, i.e., fluid displacement and frequency. However, a proper quantification of frequency in PHPs is hampered by the nonperiodicity of fluid oscillations [START_REF] Ishii | Temperature visualization and investigation inside evaporator of pulsating heat pipe using temperaturesensitive paint[END_REF]. Consequently, all the oscillation frequencies experienced by the flow pattern during the PHPs operation are usually represented by a single frequency, named dominant (or characteristic) oscillation frequency, defined as a peak in the power spectrum of a given signal [START_REF] Perna | Flow characterization of a pulsating heat pipe through the wavelet analysis of pressure signals[END_REF]. Generally, the dominant fluid oscillation frequency is evaluated by means of Fourier transform [START_REF] Xu | Start-up and steady thermal oscillation of a pulsating heat pipe[END_REF][START_REF] Karthikeyan | Understanding thermo-fluidic characteristics of a glass tube closed loop pulsating heat pipe: Flow patterns and fluid oscillations[END_REF][START_REF] Takawale | A comparative study of flow regimes and thermal performance between flat plate pulsating heat pipe and capillary tube pulsating heat pipe[END_REF][START_REF] Mameli | Local heat transfer measurement and thermo-fluid characterization of a pulsating heat pipe[END_REF], short-time Fourier transform [START_REF] Monroe | Analysis and comparison of internal and external temperature measurements of a tubular oscillating heat pipe[END_REF] or time-strip technique [START_REF] Fairley | Time-frequency analysis of flat-plate oscillating heat pipes[END_REF] on wall temperature measurements or optical acquisitions on the fluid flow. Perna et al. [START_REF] Perna | Flow characterization of a pulsating heat pipe through the wavelet analysis of pressure signals[END_REF] suggested the use of the wavelet method for a reliable time-frequency localization of dominant fluid oscillation frequencies, which is not always obtainable from the Fourier or short-time Fourier transforms. In their analysis, the fluid pressure signals, measured at two different points of a multi-turn PHP channel and sampled at 200 Hz, were proven to be suitable for the fluid oscillations assessment. However, the intrusive insertion of pressure transducers inside the fluid stream is generally hard to achieve from a manufacturing standpoint. To deal with this issue, Pagliarini et al. [START_REF] Pagliarini | Thermal characterization of a multiturn pulsating heat pipe in microgravity conditions: Statistical approach to the local wall-to-fluid heat flux[END_REF] replicated the analysis presented in [START_REF] Perna | Flow characterization of a pulsating heat pipe through the wavelet analysis of pressure signals[END_REF] on the same device under the same working conditions by adopting instead wall-to-fluid heat flux signals as inputs for the wavelet method, evaluated from thermographic acquisitions, i.e., non-intrusive measurements, on the outer wall surface. The outcomes denoted that the local wall-tofluid thermal interactions allow for a similar evaluation of dominant fluid oscillation frequencies to that provided by pressure signals. In addition, the experimental approach reliably quantified the oscillatory behavior in every device branch, thus representing an effective solution for the study of PHPs when fluid pressure measurements in multiple branches cannot be achieved. Iwata et al. [START_REF] Iwata | Characterization of thermal behavior of a micro pulsating heat pipe by local heat transfer investigation[END_REF] performed a similar analysis to characterize the oscillatory behavior in a tubular micro-PHP tested in a vertical position. The IR acquisitions were performed within the condenser section of the device, dissipating heat to the environment by radiation and free convection with air. Fluid oscillations became stronger at high heat loads both in terms of power and frequency, and every channel exhibited similar oscillatory behavior.

While most of the fluid oscillation frequency investigations have dealt with tubular PHPs, i.e., PHPs constituted by a pipe arranged in a serpentine manner, some researchers have focused on the thermofluidic behavior of other configurations, such as flat-plate PHPs, i.e., PHPs with generally rectangular/squared channels that are directly machined on metallic/polymeric plates [START_REF] Ayel | Flat plate pulsating heat pipes: A review on the thermohydraulic principles, thermal performances and open issues[END_REF]. In particular, FPPHPs exhibit a rather different thermal response with respect to tubular layouts. First, the conduction through their solid domain promotes heat dissipation at the evaporator, even though transverse interactions between adjacent channels may hamper the device's operation [START_REF] Khandekar | Thermofluid dynamic study of flat-plate closed-loop pulsating heat pipes[END_REF]. Secondly, the non-circular section of the channels may modify the overall thermofluidic behavior. Due to the manufacturing complexity of inserting pressure probes inside the fluid stream, the fluiddynamic behavior of FPPHPs is commonly studied by means of direct fluid visualizations. Slobodeniuk et al. [START_REF] Slobodeniuk | Experimental study of the flat plate pulsating heat pipe operation during dry-out and flow re-activation periods under microgravity conditions[END_REF] investigated stopovers of a seven-turn FPPHP under micro-gravity conditions through high-speed visualization. By means of a liquid-vapor interfaces tracking technique, previously introduced by Pietrasanta et al. [START_REF] Pietrasanta | Developing flow pattern maps for accelerated two-phase capillary flows[END_REF], the menisci displacement was assessed and both fluid velocity and fluid acceleration were evaluated. Miyazaki and Arikawa [START_REF] Miyazaki | Oscillatory flow in oscillating heat pipe[END_REF] highlighted that starting from the dried condition, the phase of liquid plugs oscillations is not independent in every channel, and the menisci displacement propagates in a waveform. Due to the spring effect of the vapor bubbles, the fluid motion is sequentially transmitted from one branch to another along the PHP. Moreover, they showed that the oscillation frequency is influenced by the heat flux and heat energy is accumulated in the solid walls. From a direct visualization of a horizontally oriented FPPHP, Borgmeyer and Ma [START_REF] Borgmeyer | Experimental Investigation of Oscillating Motions in a Flat Plate Pulsating Heat Pipe[END_REF] found that the nature of the working fluid significantly affects the oscillating motion of liquid slugs. With water, no oscillations were observed, whereas large oscillation amplitudes and velocities (up to 5 cm/s) were reached with ethanol under same operating conditions. On the other hand, Pagnoni et al. [START_REF] Pagnoni | Experimental and numerical analysis of a closed-loop flat-plate pulsating heat pipe: Effects of orientation with respect to gravity on the thermohydraulic behavior[END_REF] found maximum velocities around 0.6 and 1 m/s, and frequencies from 2.5 up to 4 Hz, for heat powers of 80 W and 200 W, respectively, in their vertically oriented FPPHP with horizontal channels and filled with ethanol. From a comparison between a multi-turn, tubular PHP and a FPPHP having transparent inserts, Takawale et al. [START_REF] Takawale | A comparative study of flow regimes and thermal performance between flat plate pulsating heat pipe and capillary tube pulsating heat pipe[END_REF] observed that oscillation amplitudes increased with the heat load to the evaporator in both configurations, although the tubular PHP exhibited much larger amplitudes. The absolute measured fluid velocities were about 8 to 10 times higher than those inside the FPPHP. The tubular PHP denoted more vigorous, large-amplitude oscillations with respect to the flat-plate layout. Jun and Kim [START_REF] Jun | Experimental study on a criterion for normal operation of pulsating heat pipes in a horizontal orientation[END_REF] proposed a mass-spring-damper model for the fluid displacement description in FPHPs during dry-out phases and until the large amplitude oscillations mode. The model was validated by coupling high-speed videos with fluid temperature/pressure measurements on a ten-turn micro FPPHP, thus highlighting optimal conditions for normal device operation in a horizontal orientation. Yoon and Kim [START_REF] Yoon | Experimental and theoretical studies on oscillation frequencies of liquid slugs in micro pulsating heat pipes[END_REF] designed and experimentally validated, with rather good accuracy, a mass-springdamper model that considered not only volume variations in slugs and plugs, but also mass variations due to evaporation/condensation phenomena. Such a model has also been proven to perform better than the usual ones where mass variations are generally neglected.

Although FPPHPs with transparent windows have been extensively studied in terms of flow characteristics, to the authors' knowledge, a significant amount of data regarding fluid oscillations in fully metallic devices is missing in the literature. This represents an extreme limitation for the deep understanding of FPPHPs heat transfer modes. On one hand, designing post-processing techniques for only direct visualizations precludes the possibility of extending the investigations to every possible device. On the other hand, the conductive phenomena largely occurring in FPPHPs might be influenced by the different thermal properties of transparent inserts, thus preventing robust comparisons with fully metallic devices. Moreover, such a lack of insight into fluid oscillation frequencies undermines the validation and improvement of existing numerical and theoretical models on FPPHPs. For instance, Dreiling et al. [START_REF] Dreiling | Thermal resistance modeling of flat plate pulsating heat pipes[END_REF] designed a three-dimensional thermal conduction-based model for pseudo-steady state operation to predict the FPPHP thermal resistance. The effect of inner fluid oscillations has been considered negligible due to the high uniformity of wall temperature usually observed in flat layouts. However, such an assumption may not be always true, especially when the device undergoes strong transient behavior in terms of intermittent operation or the occurrence of local hot spots [START_REF] Odagiri | Three-dimensional heat transfer analysis of flat-plate oscillating heat pipes[END_REF]. In [START_REF] Yoon | Experimental and theoretical studies on oscillation frequencies of liquid slugs in micro pulsating heat pipes[END_REF], since the investigation was conducted for micro-geometries, the resulting fluid oscillation frequencies were found to be much higher than those observed in regular, metallic layouts, thus hampering the capability of the adopted mass-spring-damper model to predict the oscillatory behavior in metallic FPPHPs. For the sake of completeness, it has to be stressed that some works have dealt with the direct observation of fluid motion inside fully metallic FPPHPs through neutron radiography [START_REF] Thompson | Investigation of a flat-plate oscillating heat pipe with Tesla-type check valves[END_REF][START_REF] Yasuda | Visualization of the working fluid in a flat-plate pulsating heat pipe by neutron radiography[END_REF][START_REF] Zhang | Effect of Hydrophilic Nanostructured Cupric Oxide Surfaces on the Heat Transport Capability of a Flat-Plate Oscillating Heat Pipe[END_REF]. However, such an experimental method is not easily implementable, especially due to the extremely high cost and complexity of instrumentations and peripheral facilities. Hence, the present work proposes a novel experimental approach for the evaluation of dominant fluid oscillation frequencies in fully metallic FPPHPs by means of infrared (IR) thermography. The aim is to provide a straightforward and effective technique for the estimation of fluid oscillation frequencies in fully opaque flat-plate layouts (i.e., where direct fluid visualizations cannot be carried out) to obtain interesting pieces of data regarding the inner fluid dynamics and to enlarge the overall literature dataset. The procedure is designed as an extension of the one presented in [START_REF] Pagliarini | Thermal characterization of a multiturn pulsating heat pipe in microgravity conditions: Statistical approach to the local wall-to-fluid heat flux[END_REF] for a metallic tubular PHP, as well as a natural evolution of the experimental work described in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] on the same device under the same working conditions. Here, the authors thoroughly explained a novel method for the evaluation of wall-to-fluid heat fluxes along the adiabatic section of an eight-turn FPPHP made of copper, starting from IR acquisitions on the outer wall surface at varying heat loads to the evaporator and tilt angles. In the current analysis, the evaluated wall-to-fluid heat flux distributions are used as inputs for the wavelet method to quantify dominant fluid oscillation frequencies in every branch of the device. The time evolution of oscillatory behavior is further quantified by such a non-intrusive technique, which does not rely on transparent windows nor on fluid pressure/temperature sensors. A statistical reduction of the results provides a complete description of the overall FPPHP oscillatory behavior, in accordance with the global and local heat transfer characteristics highlighted in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]. Finally, a link between transverse interactions and fluid oscillations is presented to provide an overview on the thermo-fluid dynamics of the studied device. To the authors' knowledge, this represents one of the first attempts of experimentally quantifying fluid oscillatory features in metallic FPPHPs through infrared observations. The newly proposed method, based on the reduction in local heat fluxes for the estimation of fluid oscillation frequencies in FPPHPs, provides a good degree of insight into the oscillatory behavior of the investigated devices without any direct fluid visualization. This strongly simplifies the experimental set-ups generally used for the study of fluid dynamic quantities in such devices. The adopted technique can be applied to a wide range of other flat-plate systems having fully metallic walls for a robust and repeatable study of the inner fluid dynamics, i.e., even when no transparent windows are present.

Experimental Apparatus

The studied device is an eight-turn FPPHP (overall dimensions: 80 × 200 × 4 mm) entirely made of copper and with 3 × 3 mm² squared channels (Figure 1a). One side of the device is coated with a black high-emissivity paint (𝜀 = 0.92) to allow thermographic measurements on the outer wall temperature within the adiabatic section. This first plate was covered and brazed with tin to a second copper plate having the same dimensions and 0.5 mm of thickness. The reader is referred to [START_REF] Ayel | Thermal performances of a flat-plate pulsating heat pipe tested with water, aqueous mixtures and surfactants[END_REF] for the description of the manufacturing process.

The opaque paint emissivity was evaluated by performing dedicated IR acquisitions on a coated isothermal surface (same opaque paint, at controlled temperature and monitored through thermocouples). The emissivity parameter set in the IR camera software was regulated until a match between the temperature measured by thermocouples and the one by the infrared camera was reached in the range of absolute uncertainty of the camera (± 1 °C). Since the emissive layer had a thickness lower than 100 μm, its thermal resistance was considered negligible, thus reasonably assuming that the temperature measured on top of the coating corresponded to the one of the device's surfaces. To guarantee full accessibility of the adiabatic section in terms of thermographic measurements, such a portion was not covered with insulating material.

Power is provided to the evaporator zone by a wire electrical heater (Thermocoax ® Type NcAc15, 1 mm external diameter), connected to an electrical power supply (ELC ® ALR3220, ±10 mV). The condenser section is cooled by a loop connected to a copper cold plate (width: 100 mm, length: 80 mm), directly brazed on the rear side of the FPPHP (Figure 1a). The loop temperature is controlled by a cryostat (Huber ® CC240 wl). Two pipes are brazed on the top and bottom sides of the device, allowing the filling and vacuuming of the system, respectively. The FPPHP was filled with a water/ethanol mixture (20 wt.% of ethanol) with a filling ratio of 50% (at ambient temperature). Specifically, such a working fluid was chosen since it led, as proven in [START_REF] Ayel | Thermal performances of a flat-plate pulsating heat pipe tested with water, aqueous mixtures and surfactants[END_REF]-for the same FPPHP and among a selection of various fluids (pure water, pure ethanol, water aqueous mixtures with 5% of butan-1-ol, 5% of butan-2-ol, 16% and 20% of ethanol, 28% of methanol and 0.5% of surfactants-Tween ® 20 and Tween ® 40)-to better thermal performances of the device at different tilt angles and working conditions, making the adopted water/ethanol aqueous mixture suitable for robust and repeatable fluid oscillation analyses. Part of the adiabatic section was monitored by means of a high-speed Medium Weave InfraRed camera (MWIR, FLIR ® SC5500, 320 × 256 pixels resolution). In Table 1, the uncertainties of the adopted experimental equipment are listed. The reference for the channels' numeration is reported in Figure 1b, together with the heat transfer sections, i.e., evaporator (red box, 40 mm long) and condenser (blue box, 80 mm long), and the area framed by the camera (green box, 40 mm long). The reference for the axial direction z is also shown. The analysis was performed for two device orientations, namely the vertical bottom heated mode (BHM) and the horizontal orientation (Figure 2), and different heat loads Q applied to the evaporator, spanning from 50 to 250 W. During the tests, IR acquisitions, each lasting 20 s, were performed with a sampling frequency of 50 Hz in the so-called pseudo-steady state operation of the device, i.e., after about 30 min from the delivery of each power step. Specifically, IR thermography was chosen, among other experimental methods, since it provided an accurate estimation of wall temperature variations for effective wall-to-fluid heat flux evaluations. In fact, using other experimental means, such as thermocouples attached on the outer wall of the device, would lead to a poorer assessment of thermal transients due to intrinsic limitations, including slower thermal response, lower accuracy and high contact resistance. When carrying out IR acquisitions, the effects of IR radiations coming from the environment and partially reflected by the surface of interest (belonging to the coated adiabatic section of the device) were considered as negligible, i.e., they did not substantially affect the IR measurements. This reasonable assumption is supported by the fact that during the tests, the device wall temperature was at least 13 °C higher than the environmental one, and the adopted, opaque coating was poorly reflective.

Methods

Wall-to-Fluid Heat Flux Estimation

The procedure for the estimation of the heat flux exchanged between the FPPHP thin walls and the working fluid is based on the post-processing of the thermographic data by the inverse heat conduction problem (IHCP) resolution approach, coupled with a threedimensional Gaussian filter for the regularization of the raw space-time temperature maps. Such a procedure is carried out within the adiabatic section of the device, and it is thoroughly explained in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF], together with its error estimation. In Figure 3, the areas of interest for the wall-to-fluid heat flux estimation, i.e., where the FPPHP channels present thin walls, are highlighted in red. The resulting wall-to-fluid heat flux distributions q are bi-dimensional maps referred to each FPPHP channel over 0.04 m (length of the section of interest) and 20 s (duration of each MWIR acquisition). In particular, by considering the energy balance equation in transient conditions at the thin wall element and by rearranging all the terms, the wall-tofluid heat flux finally assumes the following formulation [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]:

𝑞 = 𝑠 𝜌𝑐 𝜕𝑇 𝜕𝑡 -𝑘 𝜕 𝑇 𝜕𝑦 -𝑘 𝜕 𝑇 𝜕𝑧 + ℎ (𝑇 -𝑇 ) (1)
where s is the wall thickness, ρ, cp and k are density, specific heat and conductivity of copper, respectively, henv is the outer heat transfer coefficient by free convection and radiation to the environment, and Tenv is the environmental temperature. To obtain consistent values of heat flux q, the input temperature is filtered by following the steps below [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]:


The three-dimensional, raw space-time temperature map is converted to the frequency domain by the discrete fast Fourier transform (DFFT); 

The transfer function for the Gaussian filter is applied to the frequency image of the input temperature map. A first attempt cut-off frequency is initially guessed; 

The temperature map is converted back to the space-time domain by applying the inverse DFFT; 

The residual between the raw input temperature map and the filtered temperature is computed; 

Steps 2-4 are iteratively replicated by varying the cut-off frequency until the evaluated residual equals the standard deviation of the raw temperature data, σ = 0.04 K (i.e., fulfilment of the Morozov's discrepancy principle [START_REF] Morozov | Methods for Solving Incorrectly Posed Problems[END_REF]).

The goodness of the estimation procedure is therefore assessed by means of the following steps [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]:



The direct problem is implemented in a COMSOL Multiphysics ® environment. The solid domain of the device within the adiabatic section is simulated by using, as a boundary condition at the inner wall-fluid interface, a given space-time heat flux distribution that reflects the one exchanged in the real device. Equations for threedimensional transient conduction are solved by the in-built solver; 

The synthetic temperature distribution resulting from the simulation is then spoiled with a noise level equal to the standard deviation of the raw data, and processed by means of the filtering technique; 

The wall-to-fluid heat flux is estimated back by means of Equation (1), i.e., leading to the inverse problem solution. The estimated wall-to-fluid heat flux is therefore compared with the one given as a boundary condition in Step 1, allowing the evaluation of the estimation error for the procedure (quantified in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF], equal to about 15% of the heat flux given in Step 1).

Oscillation Frequency Estimation

The fluid oscillation frequencies were evaluated starting from the wall-to-fluid heat flux distributions (Section 3.1), since the local thermofluidic interactions were believed to be strictly linked to the heat dissipation process underlying the FPPHP operation [START_REF] Pagliarini | Thermal characterization of a multiturn pulsating heat pipe in microgravity conditions: Statistical approach to the local wall-to-fluid heat flux[END_REF][START_REF] Cattani | An original look into pulsating heat pipes: Inverse heat conduction approach for assessing the thermal behaviour[END_REF]. As already stated, the oscillatory behavior quantification was based on the wavelet method; in particular, according to Perna et al. [START_REF] Perna | Flow characterization of a pulsating heat pipe through the wavelet analysis of pressure signals[END_REF], the Morlet wavelet, the shape of which is defined by Equation (2) for dimensionless time η and the characteristic parameter ω0, was adopted.

𝜓 (𝜂) = 1 √𝜋 𝑒 𝑒 (2)
The wavelet transform of a time-dependent signal g(t) reads as [START_REF] Lilly | Element analysis: A wavelet-based method for analysing time-localized events in noisy time series[END_REF] 𝑊 (𝑎, 𝜏) = 1

√𝑎 𝑔(𝑡)𝜓 * 𝑡 -𝜏 𝑎 𝑑𝑡 (3)
where a and τ are the wavelet scale (dilation of the wavelet function) and time shift, respectively, and the superscript * denotes the complex conjugate. The term 1/√𝑎 is due to the adopted L1 normalization of the wavelet, used for comparison purposes between different scales. The output of the wavelet analysis is the so-called magnitude scalogram, defined as the absolute value of the wavelet transform, plotted as a function of frequency and time. In particular, it represents the power related to each frequency over the sample duration. Such a power is strictly linked to the oscillation amplitude of the signal; in fact, high-amplitude oscillations result in high powers in the scalogram. On the other hand, the power spectrum referred to as the considered signal is defined as the time integral of the magnitude scalogram, for every frequency. While the scalogram describes the time evolution of the characteristic frequencies underlying the considered oscillatory phenomenon, the power spectrum allows for the identification of the dominant frequency related to the considered sample (identified by a peak in the power spectrum).

The wavelet parameters adopted in the present analysis are listed in Table 2. In particular, these parameters were chosen to guarantee a good trade-off between the frequency and time resolution of the method.

Table 2. Parameters used for the wavelet method.

A [1/Hz] 71 𝜔 [rad/s] 2π fmin [Hz] 0.1 fmax [Hz] 18

Validation of the Dominant Frequency Estimation Procedure

The method of Section 3.2 was validated following a similar validation procedure to that presented in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] and summarized in Section 3.1. Direct simulations were implemented in a COMSOL Multiphysics ® environment to solve conduction heat transfer equations within the FPPHP solid domain. In particular, the same analytical heat flux distributions 𝑞 𝑒𝑥𝑎𝑐𝑡 used for the validation procedure of [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] were given as a boundary condition at the wall-fluid interfaces. Specifically, the guessed distribution for the analytical heat flux 𝑞 𝑒𝑥𝑎𝑐𝑡 assumed the following sinusoidal form [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]:

𝑞 𝑒𝑥𝑎𝑐𝑡 = 𝐴𝑚𝑝 • 𝑐𝑜𝑠(2𝜋𝑡 • 𝑓 -𝜑) + 𝐴𝑚𝑝 2 - 𝑧 𝐻 • 𝐴𝑚𝑝 (4)
where Amp is the heat flux amplitude in W/m 2 , f is the oscillation frequency in Hz, φ is the phase shift (different from channel to channel to reflect the oscillatory behavior in the real device), z is the axial coordinate and H is the length of the portion of interest within the adiabatic section. The reader is referred to [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] for deeper details on the procedure. The resulting synthetic temperature maps were therefore spoiled with noise, filtered and used as inputs for the inverse procedure to evaluate back the heat flux, qrestored. In other words, noisy temperature distributions were used as inputs for Equation (1), solved through the finite element method. However, while in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] the validation procedure was focused on the percentage error estimation in terms of amplitude of the restored heat fluxes with respect to the exact ones (i.e., the guessed ones), the present validation aims at considering differences in oscillation frequency between qexact and qrestored. Hence, qexact and qrestored were used as inputs for the wavelet method.

In Figure 4, qexact (amplitude Amp = 2000 W/m 2 , frequency f = 1.2 Hz) is shown, together with the wavelet method outputs. Due to the fact that the analytical formulation of qexact, for which the time-dependent signal is presented in Figure 4a, is a sinusoid having a fixed frequency over time (please refer to [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]), the scalogram (Figure 4b) exhibits the same signal power (described by the color bar on the right) at each time instant of the observation window. It has to be stressed that the white dashed lines in scalograms delimit what is known as cone of influence (COI) [START_REF] Nobach | Review of some fundamentals of data processing[END_REF]. Specifically, the COI highlights areas (grey-shaded) potentially affected by edge-effects, i.e., where the wavelet extends over the boundaries of the considered time interval, and they might not be representative of an accurate time-frequency description of the considered samples. For this reason, the power spectrum evaluation does not consider such areas when integrating the power related to each frequency over time. In Figure 4c, the power spectrum is reported: here, the power related to low (<0.8 Hz) and high (>2 Hz) frequencies Pw is null, while it increases as long as the natural frequency of the treated signal (1.2 Hz) is approached. The absolute maximum of the power spectrum is found at 1.23 Hz. Hence, the error for the dominant frequency estimation procedure, i.e., the frequency resolution of the wavelet method for the present analysis, was assessed to be equal to 0.03 Hz. Such an error is intrinsic in the wavelet approach due to the shape of the considered wavelet function, having coefficients listed in Table 2. The lower limit of 0.5 Hz for the frequency scales has been chosen due to the limited observation window (20 s), which prevents the detection of a sufficient number of oscillations at lower frequencies, undermining the reliability of the method. In Figure 5, the same frequency analysis is presented for the restored heat flux referred to noise level σ = 0.04 K (noise level of the experimental set-up, as defined in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]). The scalogram of Figure 5b, when compared to that of Figure 4b, suggests that some lowand high-frequency components are introduced by the restoration procedure. In fact, the scalogram is no longer represented by just a clear, high-power horizontal band, but some other low/average-power spots are perceivable. This effect is clearly noticeable in the power spectrum (Figure 5c), where its maximum decreases of about 20% of the Pw peak in Figure 4c mean that the power of the exact signal has been scattered over a wider range

(a) (b) (c) [W/m 2 /Hz 1/2 ] P w [W/m 2 /Hz 3/2 ]
of frequencies. However, the maximum of the power spectrum is still at 1.23 Hz, suggesting that no significant errors are introduced by the heat flux restoration in terms of dominant frequency estimation. Such a remark was verified for all the considered qexact and qrestored, thus further proving that the presented procedure outcomes are not substantially dependent on the amplitude nor on the natural frequency of qexact.

Results and Discussion

Fluid Oscillatory Behavior

The local thermal interactions between the working fluid and the FPPHP wall are representative of the fluid motion, i.e., fluid oscillations, inside the device, as demonstrated in [START_REF] Pagliarini | Thermal characterization of a multiturn pulsating heat pipe in microgravity conditions: Statistical approach to the local wall-to-fluid heat flux[END_REF]. Hence, the wavelet approach was used for the evaluation of dominant frequencies fd in the overall studied device, i.e., every FPPHP channel, by adopting wall-to-fluid heat flux distributions for single axial coordinates over time. In Figure 6, the post-processed heat flux signal referred to z = 0.02 m for channel 8 (reference of Figure 1b) and Q = 100 W in the horizontal orientation is shown as a sample for the wavelet method outcomes understanding of experimental data. The experimental-based signal (Figure 6a) highlights non-fixed amplitude and oscillation frequency over time. In fact, the scalogram of Figure 6b exhibits scattered spots having different power in the range from 0.8 to 2 Hz. From 0 to about 10 s, the high-power spots are mainly present around 1.5 Hz, denoting that the heat flux, i.e., the working fluid, oscillates with great amplitude at that frequency. From 10 to 14 s, the scalogram does not present significant power spots, thus suggesting the occurrence of weak heat flux oscillations. From 14 to 18 s, a high-power spot is perceivable around 1 Hz, denoting a substantial change in terms of oscillation frequency of the heat flux. Such a time change in oscillation frequency is also clear in the power spectrum (Figure 6c), where two peaks are present: one at about 1 Hz and one at about 1.5 Hz. Here, although there is a shorttime change in oscillatory behavior, the heat flux mainly oscillates with variable amplitudes at about 1.5 Hz over the entire observation window, and the absolute maximum of the power spectrum is located around such frequency. In particular, the dominant oscillation frequency fd of the signal is found at about 1.45 Hz.

It has to be pointed out that since the heat flux distributions evaluated in the present device significantly vary with space [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF], the wavelet method may provide different results depending on the considered axial coordinate. To investigate such dependency, the wall-to-fluid heat flux signals related to a single channel at different z coordinates were used as inputs for the wavelet approach, and the resulting power spectra were therefore compared. In Figure 7, such a comparison was carried out for channel 9 at z = 0.01 m, z = 0.02 m and z = 0.03 m, and both device orientations at Q = 250 W. Here, despite different amplitudes of the wall-to-fluid heat fluxes of Figures 7a-c, their oscillatory trends over time do not present any dissimilarities from axial point to axial point. In fact, all the power spectra of Figures 7b-d exhibit absolute peaks in correspondence with the same frequency, thus highlighting same the dominant frequency regardless of the considered axial position. For this reason, every processed heat flux was chosen, for the present analysis, to be at z = 0.02 m for every channel, representative of all other z coordinates.

Channel-Wise Frequency Analysis

The wavelet method was extended to the heat flux signals referred to all the FPPHP branches, thus assessing channel-wise fluid oscillation frequencies at varying heat loads and orientations.

In Figure 8, the scalograms referred to three meaningful channels, namely channels 4, 8 and 12 (Ch4, Ch8 and Ch12, respectively), which are shown for the vertical BHM at power inputs equal to 50 W, 150 W and 250 W. Specifically, the three channels only were chosen to achieve a good trade-off between a clear spatial overview of the global oscillatory behavior in the studied device and a better readability of the discussion. The power spectra related to every FPPHP channel are further presented for a comprehensive description of the oscillatory behavior in Figure 9 for the three considered heat loads. Note that magnitude scalograms in Figure 8 and power spectra in Figure 9 present, respectively, the same scale for better comparison. In Figure 8, Q = 50 W, the scalogram of Ch4 exhibits a medium-power area during the first seconds of acquisitions, suggesting the presence of weak oscillations. On the contrary, Ch8 and Ch12 do not present any significant power spots except for some low-power areas, thus highlighting a near absence of fluid motion. Such evidence is confirmed by the power spectra of Figure 9 (Q = 50 W), which are flat for most of the channels, i.e., no fluid oscillations can be perceived in almost the whole device. As the power input to the evaporator increases to 150 W, the scalograms of Figure 8, referred to the vertical BHM, present more persistent medium-power areas over time; additionally, all the power spectra of Figure 9 (Q = 150 W) exhibit perceivable peaks, suggesting that the fluid oscillates in every branch with dominant frequencies mainly below 1 Hz. For the maximum provided power input (Q = 250 W), the scalograms of Figure 8 highlight many high-power spots, denoting strong fluid oscillations over time. Nevertheless, the power spots are referred to different frequencies, suggesting that the oscillatory behavior is not stable over the observation window. The presence of strong oscillations described by high-power spots is confirmed by the power spectra of Figure 9 (Q = 250 W), which assume high peaks of Pw. Their absolute maxima are furthermore scattered within a wide range of frequencies between 0.6 and 1.6 Hz, thus denoting greatly different dominant fluid oscillation frequencies from channel to channel. The observed progressive increase in fluid motion in every channel confirms the typical flow regimes in PHPs [START_REF] Cui | Combination study of operation characteristics and heat transfer mechanism for pulsating heat pipe[END_REF][START_REF] Karthikeyan | Infrared thermography of a pulsating heat pipe: Flow regimes and multiple steady states[END_REF][START_REF] Mameli | Thermal instability of a Closed Loop Pulsating Heat Pipe: Combined effect of orientation and filling ratio[END_REF]: at low heat loads, the device operation is mainly intermittent, presenting fluid oscillations in few channels; meanwhile, at high heat loads, every branch presents strong oscillations, although they are not persistent over the observation window. Again, the scalograms and power spectra scale are adjusted to match the one at the highest heat load to the evaporator (Q = 250 W) to allow for a better and more straightforward comparison between different power inputs.

In Figure 10 at Q = 50 W, the scalograms do not present any significant power spots, denoting either the absence of fluid motion or very weak and intermittent oscillations. The device inactivity over the observation window is additionally highlighted by a complete absence of peaks in the power spectra of Figure 11 (Q = 50 W). It has to be stressed that since the color and Pw scales of the provided scalograms and power spectra are adjusted for comparison to match the ones at the highest power input, the values seem to be generally null. However, by observing, in Figure 12, a representative scalogram of Ch 4, Figure 10, 50 W, having an adjusted color scale around zero, it can be noted that the graph presents non-null power areas. Nonetheless, such non-null, low-amplitude regions account for residual noise in the input heat flux data rather than actual fluid oscillations. The same conclusion can be drawn for the other scalograms at 50 W, horizontal orientation, thus confirming the absence of fluid motion in the considered channels for the analyzed heat load. Such a remark agrees with what was observed in the previous literature [START_REF] Ayel | Thermal performances of a flat-plate pulsating heat pipe tested with water, aqueous mixtures and surfactants[END_REF][START_REF] Betancur-Arboleda | Experimental study of channel roughness effect in diffusion bonded pulsating heat pipes[END_REF][START_REF] Jun | Comparison of the thermal performances and flow characteristics between closed-loop and closed-end micro pulsating heat pipes[END_REF], where operating FPPHPs in a horizontal orientation led to delayed start-ups at higher power inputs due to a lack of pumping forces assisted by gravity effects. For Q = 150 W, the scalograms of Figure 10 exhibit medium-power stripes which persist with continuity during almost the whole acquisition in the considered channels, except for some higher-frequency spots that are especially evident for Ch4. High Pw are clearly noticeable from the power spectra of Figure 11 (Q = 150 W), confirming the presence of regular fluid oscillations with dominant frequencies around 0.9 Hz. At 250 W, the scalograms of Figure 10 show high-power spots which replicate at almost the same frequency over time, denoting the time regularity of the strong fluid oscillations occurring in such operating conditions. Such vigorous fluid oscillations in every FPPHP channel are also clear from the power spectra high peaks of Figure 11 (Q = 250 W), with dominant frequencies around 1.2 Hz. It is worth noting that, contrarily to the vertical BHM orientation, the peaks of Pw span a much narrower range of frequencies for every channel in the horizontal orientation for medium/high power inputs. This might be due to the fact that the buoyancy forces acting on the fluid in BHM increase the complexity of the device's operation and, consequently, of its oscillatory behavior. In fact, in the vertical BHM, a transition between slug-plug and semi-annular flow may occur, especially due to the squared section of the present FPPHP channels [START_REF] Khandekar | Thermofluid dynamic study of flat-plate closed-loop pulsating heat pipes[END_REF][START_REF] Kearney | Thermal performance of a PCB embedded pulsating heat pipe for power electronics applications[END_REF][START_REF] Qu | Flow visualization of silicon-based micro pulsating heat pipes[END_REF], thus increasing the chaoticity of the thermofluidic wall-to-fluid interactions. On the contrary, in the horizontal orientation, the flow pattern could probably be only slug-plug even for high power inputs to the evaporator, i.e., high inertial forces on the fluid driven by damped oscillations of the vapor bubbles, thus resulting in a similar oscillatory behavior between every FPPHP branch. In fact, a slug-plug flow is more likely to occur in a horizontal orientation due to absence of buoyancy forces which tend to deform the liquid-vapor interfaces. For instance, Ayel et al. [START_REF] Ayel | Visualization of Flow Patterns in Closed Loop Flat Plate Pulsating Heat Pipe Acting as Hybrid Thermosyphons under Various Gravity Levels[END_REF] observed a slug-plug flow during microgravity experiments on a closed-loop FPPHP with channel dimensions of 2.5 × 2.5 m² and FC-72 as a working fluid. Of note, weightlessness acts similarly to the horizontal orientation in terms of the negligible effects of buoyancy forces on the fluid stream, as proven in [START_REF] Mameli | Thermal response of a closed loop pulsating heat pipe under a varying gravity force[END_REF]. In [START_REF] Ayel | Visualization of Flow Patterns in Closed Loop Flat Plate Pulsating Heat Pipe Acting as Hybrid Thermosyphons under Various Gravity Levels[END_REF], when a transition from hyper-gravity to microgravity occurred, the flow pattern varied from annular to slugplug, confirming the higher stability of fluid menisci and overall fluid flow observed during horizontal operation in the present investigation.

To clearly provide a quantitative description of the channel-wise oscillatory behavior, all the evaluated dominant oscillation frequencies from the power spectra curves of Figures 9-11 are listed in Tables 3 and4 for the vertical BHM and horizontal orientations, respectively. In Table 3, at low heat loads applied to the evaporator, the vertical BHM orientation denotes several channels in which the fluid does not present any dominant oscillation within the observation window, thus suggesting that the high intermittency of the overall device operation prevents any identification of the dominant frequencies. At average-high heat loads, a dominant oscillation frequency can be perceived in every channel, confirming the transition from intermittent flow to full activation of the device. In the horizontal orientation (Table 4), a lack of fluid motion hampers a clear identification of a dominant frequency in the overall device at Q = 50 W, while every FPPHP branch denotes fluid oscillations from 100 W up to the maximum provided power input, similarly confirming the transition between the two mentioned working regimes. The listed fluid oscillation frequencies are in agreement with those typically found with direct fluid visualizations [START_REF] Spinato | Understanding the self-sustained oscillating two-phase flow motion in a closed loop pulsating heat pipe[END_REF][START_REF] Xu | Start-up and steady thermal oscillation of a pulsating heat pipe[END_REF][START_REF] Ayel | Flat plate pulsating heat pipes: A review on the thermohydraulic principles, thermal performances and open issues[END_REF][START_REF] Kim | The study on pressure oscillation and heat transfer characteristics of oscillating capillary tube heat pipe[END_REF][START_REF] Khandekar | Multiple quasi-steady states in a closed loop pulsating heat pipe[END_REF]. In addition, they confirm the working regimes identified in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] by means of standard techniques, i.e., evaporator temperature fluctuations, and statistical reduction in wall-to-fluid heat fluxes, for which the overview is reported in Table 5.

Table 5. Device working regimes identified in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]. "I": intermittent flow, "F": full activation.

Orientation Q [W]

Vertical

Time Evolution of Fluid Oscillations

The FPPHP oscillatory behavior over time, which was qualitatively observed in scalograms of Section 5.2, was assessed by investigating some of its features, such as continuity and channel activity. Specifically, the continuity tc is here defined as the persistency of fluid oscillations over time around the dominant frequency, i.e., in the range Bf = fd ± 15% fd. It can be conceptually expressed, for a given device channel, as

𝑡 = 𝑇𝑖𝑚𝑒 𝑑𝑢𝑟𝑖𝑛𝑔 𝑤ℎ𝑖𝑐ℎ 𝑓𝑙𝑢𝑖𝑑 𝑜𝑠𝑐𝑖𝑙𝑙𝑎𝑡𝑒𝑠 𝑎𝑡 𝑎 𝑓𝑟𝑒𝑞𝑢𝑒𝑛𝑐𝑦 𝑐𝑙𝑜𝑠𝑒 𝑡𝑜 𝑡ℎ𝑒 𝑑𝑜𝑚𝑖𝑛𝑎𝑛𝑡 𝑜𝑛𝑒 𝑇𝑜𝑡𝑎𝑙 𝑎𝑐𝑞𝑢𝑖𝑠𝑖𝑡𝑖𝑜𝑛 𝑡𝑖𝑚𝑒 (5)
Channels activity ta is instead referred to presence of fluid oscillations:

𝑡 = 𝑇𝑖𝑚𝑒 𝑑𝑢𝑟𝑖𝑛𝑔 𝑤ℎ𝑖𝑐ℎ 𝑓𝑙𝑢𝑖𝑑 𝑜𝑠𝑐𝑖𝑙𝑙𝑎𝑡𝑖𝑜𝑛𝑠 𝑎𝑟𝑒 𝑝𝑒𝑟𝑒𝑐𝑒𝑣𝑖𝑏𝑎𝑙𝑒 𝑇𝑜𝑡𝑎𝑙 𝑎𝑐𝑞𝑢𝑖𝑠𝑖𝑡𝑖𝑜𝑛 𝑡𝑖𝑚𝑒 (6)
Both properties are believed to complete the fluid oscillations description in terms of time evolution, providing a better quantification of the device heat transfer capabilities. The evaluation of continuity and activity is based on the time-wise determination of the main oscillation frequencies, i.e., referred to as the power peaks in scalograms. However, a reliable and unambiguous evaluation of scalogram peaks is greatly hampered by the presence of low-power peaks, which are not always representative of significant fluid oscillations. Hence, a procedure for scalogram denoising, aimed at discerning meaningful from neglectable power peaks, is required. In the present analysis, the K-means clustering was employed to highlight high-power regions in scalograms. In particular, such a method relies on an iterative data-partitioning algorithm that assigns N observations to K separate clusters defined by centroids [START_REF] Frank | The Data Analysis Handbook[END_REF]. K is preliminarily chosen, depending on the specific application. The observations assignment relies on within-cluster variance minimizations through the evaluation of squared Euclidean distances from centroids, progressively moved until the minimization condition is satisfied. Given a set of observations O = {O1, O2, … On}, the partition into k ≤ K subsets S = {S1, S2, … Sk} is performed according to the following minimization problem:

𝑎𝑟𝑔𝑚𝑖𝑛 𝑺 ‖𝑠 -𝜇 ‖ ∈ 𝑺 𝒌 (7)
where μk is the mean of the k-th subset. When applied to scalograms, clusters containing residual noise, wavelet artifacts or non-relevant oscillations will present lower k coefficients [START_REF] Zhang | Effect of Hydrophilic Nanostructured Cupric Oxide Surfaces on the Heat Transport Capability of a Flat-Plate Oscillating Heat Pipe[END_REF]. Areas including significant oscillatory phenomena will instead present high k coefficients. Consequently, the only partitions defined by k = K were considered as meaningful regions, while other partitions were discarded. For the present analysis, the value K = 4 was chosen using, as a reference, the one provided by the Calinski-Harabasz criterion [START_REF] Caliński | A dendrite method for cluster analysis[END_REF]. Since the convergence of the K-means algorithm is affected by the initial position of centroids, their initialization was replicated until an optimal final partition of scalograms, i.e., a partition capable of identifying with good accuracy the high-power regions, was obtained. In Figure 13a,b, the scalograms of Figure 8 (250 W, Ch 12, BHM) and 10 (250 W, Ch 12, horizontal) are proposed to clearly present, from a visual standpoint, the steps of the present peak identification procedure. In Figure 13c,d, the corresponding denoised scalograms are shown: the application of the K-means method highlights high-power peaks, while the scalograms low-power background is set to 0. Once the denoising procedure was applied, the scalograms' maximum peaks, i.e., main oscillatory phenomena, could be identified, time instant by time instant. In Figure 13e,f, the time location of frequencies referred to scalogram peaks is shown over time, together with the dominant frequencies (solid lines) and Bf boundaries (dotted lines). As is noticeable, time instants that do not present high-power peaks in the threated scalograms are not assigned to any frequencies (gaps in Figure 13e,f). In addition, time periods exceeding the COIs limits are not considered. When compared to the vertical BHM (Figure 13e), the signal referring to the horizontal orientation (Figure 13f) denotes longer periods of oscillation within Bf, suggesting greater continuity over time. In particular, the continuity time tt, defined as the total time during which oscillations occur around the dominant frequency, was assessed to be equal to 2.14 s and 8.22 s for the vertical and horizontal orientations, respectively. This confirms the qualitative information provided in Section 5.2. On the other hand, the vertical BHM exhibits an overall longer activity, since oscillations at given frequencies are present for more time within the observation window. The channel activity time was evaluated to be equal to 17.32 s and 12.76 s for the vertical and horizontal orientations, respectively. The investigation on the continuity of fluid oscillations and channel activity was extended to all FPPHP channels for average and high-power inputs to the evaporator, namely Q = 150 W and Q = 250 W. Lower heat loads were not considered since either weak or absent fluid oscillations prevented a clear discernment of high-power peaks from lowpower backgrounds in scalograms. In Tables 6 and7, tc and ta are listed for the mentioned study cases to present a channel-wise description of continuity and activity. It has to be stressed that the uncertainty, referring to the evaluated times, is strictly linked to the time resolution of the obtained magnitude scalograms. In fact, the time resolution determines the error made in the estimation of the time instants at which maximum powers are identified in scalograms after the denoising procedure. For the given wavelet parameters (Table 2), the time resolution equals 0.01 s, corresponding to less than 0.1 % of the total observation windows. For both heat loads in the vertical BHM (Table 6), the continuity generally assumes different values from channel to channel with respect to the horizontal orientation (Table 7), confirming that the fluid oscillatory motion is more scattered over a wider range of frequencies. For what concerns activity, the vertical BHM denotes more similar ta values, especially at high heat load, underlining that the fluid oscillations are promoted in every channel by such orientation.

Table 6. Continuity and activity evaluated for the vertical BHM.

Channel Both continuity and activity values shown in Tables 6 and7 were therefore averaged over the entire device, i.e., for all channels, resulting in tc,av and ta,av, respectively. In Figure 14, tc,av is shown for both considered heat loads and device orientations. The vertical BHM exhibits lower values of tc,av (around 23 and 27% of the total observation window), confirming the greater continuity of oscillations in the horizontal orientation, where tc,av is equal to about 37 and 34% for average and high heat loads, respectively. However, a certain decrease in average continuity time with the increase in heat load is perceivable, denoting a reduction in fluid oscillation around the dominant frequency probably due to the onset of higher inertial forces, i.e., less stable oscillatory behavior [START_REF] Ayel | Flat plate pulsating heat pipes: A review on the thermohydraulic principles, thermal performances and open issues[END_REF][START_REF] Khandekar | Thermofluid dynamic study of flat-plate closed-loop pulsating heat pipes[END_REF][START_REF] Qu | Flow visualization of silicon-based micro pulsating heat pipes[END_REF]. ta,av is instead shown in Figure 15 for the test conditions. Here, the vertical BHM presents generally greater channel activity with respect to the horizontal orientation, since the gravity-assisted mode promotes the fluid motion in the overall device without significant stopovers, confirming what was observed in other works on FPPHPs, where the horizontal orientation usually results in a less stable operation of the device over time [START_REF] Jun | Experimental study on a criterion for normal operation of pulsating heat pipes in a horizontal orientation[END_REF][START_REF] Jun | Comparison of the thermal performances and flow characteristics between closed-loop and closed-end micro pulsating heat pipes[END_REF][START_REF] Ayel | Visualization of Flow Patterns in Closed Loop Flat Plate Pulsating Heat Pipe Acting as Hybrid Thermosyphons under Various Gravity Levels[END_REF]. Additionally, for such an orientation, the activity time increases with the power input, while the horizontal mode presents similar activity times. This remark suggests that above a certain heat load, channel activity saturates in the horizontal orientation probably due to the absence of significant buoyancy forces. Higher activity times of the vertical BHM furthermore confirm the better thermal performance assessed in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] for such an orientation.

Q [W] 150 250 tc [%] ta [%] tc [%] ta [%] 1

Statistical Description of the Oscillatory Behavior

After providing a quantification of the main properties defining the oscillatory behavior in the studied device, the evaluated channel-wise dominant frequencies were statistically reduced for a complete description of fluid oscillations. To this aim, the arithmetic mean (fd,av) and the standard deviation (fd,std) of fd were computed by considering every FPPHP channel. In Figure 16, fd,av and fd,std are shown for the only power inputs to the evaporator in which fd could be evaluated in every channel of the device, i.e., during its full activation. Here, the average dominant oscillation frequency (Figure 16a) exhibits a general increasing trend, in accordance with other studies also concerning tubular layouts [START_REF] Borgmeyer | Experimental Investigation of Oscillating Motions in a Flat Plate Pulsating Heat Pipe[END_REF][START_REF] Yasuda | Visualization of the working fluid in a flat-plate pulsating heat pipe by neutron radiography[END_REF][START_REF] Zhang | Advances and unsolved issues in pulsating heat pipes[END_REF]. For the vertical BHM orientation, fd,av ranges from about 0.78 Hz to 1 Hz, while the horizontal orientation highlights values of fd,av that increase from a minimum of about 0.84 Hz to a maximum of 1.25 Hz, confirming the typical fluid oscillation frequencies found in the literature [START_REF] Marengo | Pulsating Heat Pipes: Experimental Analysis[END_REF]. It has to be observed that while fd,av denotes a linear increment with the power input to the evaporator in the horizontal orientation, it seems to approach an asymptotic value with the increase in Q in the vertical BHM. This suggests that the semi-annular or annular flow, promoted by the gravity effects on the device at high heat loads, results in different oscillatory behavior with respect to that in the horizontal orientation (where a slug-plug flow is more likely to occur), probably due to significantly different phase-change phenomena. In fact, at high heat loads, significant fluid motion towards the condenser may occur, while liquid backflow by gravity may be, to some extent, less predominant during operation, affecting the overall oscillatory behavior, which is in agreement with [START_REF] Yasuda | Visualization of the working fluid in a flat-plate pulsating heat pipe by neutron radiography[END_REF]. Additionally, the standard deviation of fd (Figure 16b) highlights a different behavior in terms of oscillations from channel to channel in the vertical BHM orientation, which persists even during the full activation of the device. For the horizontal orientation, fd,std seems to decrease instead with the increasing of the power input, suggesting that during the full activation, every channel approaches the same oscillatory behavior. This is probably due to the fact that when buoyancy forces are negligible, as in the horizontal orientation, flow reactivations may occur in a waveform over the entire device, leading to similar fluid oscillation frequencies [START_REF] Miyazaki | Oscillatory flow in oscillating heat pipe[END_REF][START_REF] Ayel | Visualization of Flow Patterns in Closed Loop Flat Plate Pulsating Heat Pipe Acting as Hybrid Thermosyphons under Various Gravity Levels[END_REF][START_REF] Kim | Flow visualization of oscillation characteristics of liquid and vapor flow in the oscillating capillary tube heat pipe[END_REF].

Link between Fluid Oscillations and Transverse Interactions

During its operation, the thermal performance of a FPPHP might be significantly influenced by thermal interactions occurring between adjacent channels, which lead to a homogenization of temperature and, consequently, of the driven force for the working fluid motion in the overall two-phase device [START_REF] Khandekar | Thermofluid dynamic study of flat-plate closed-loop pulsating heat pipes[END_REF][START_REF] Noh | Investigation into the effect of transverse condudction on the thermal performance of a flat-plate pulsating heat pipe[END_REF]. A local heat transfer analysis could thus be required to attempt a description of conduction between adjacent channels in FPPHPs. To this aim, the temperature distributions acquired by means of thermography, and properly regularized as thoroughly explained in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF], were adopted to provide an estimation of the heat flux exchanged between the n-th and (n+1)-th channels by conduction qFou,n→n+1, according to Fourier's law (along y-axis as defined in Figure 3):

𝑞 , → (𝑧, 𝑡) = -𝑘 𝜕𝑇(𝑦, 𝑧, 𝑡) 𝜕𝑦 (8)
In particular, the partial derivative of T is evaluated in the median point between adjacent channels for every axial coordinate z as a function of time. The values averagely assumed by qFou over time were therefore adopted to estimate the power transferred by transverse conduction 𝑄 between every adjacent channel pair in the overall device along the adiabatic section (green box in Figure 1b) as follows:

𝑄 = ∑ {| ∑ 𝑚𝑒𝑎𝑛[𝑞 , → (𝑧, 𝑡 = 1 … 𝑁)] • ∆𝑧 • 𝑇 |} 𝑛𝑢𝑚 -1 (9)
where num is the total number of FPPHP channels, mean is the arithmetic mean operator, N is the total number of time samples, 𝐿 is the length of the adiabatic section framed by the IR camera, ∆z is the distance between two consecutive pixels in the z direction of the thermographic data, and T_FPPHP is the device thickness. Finally, 𝑄 was normalized with respect to the net power input to the evaporator divided by the number of FPPHP channels, resulting in 𝑄 = 𝑄 /( ). Specifically, 𝑄 defines the percentage of transverse conduction with respect to the power potentially transferrable by a single FPPHP branch from the evaporator to the condenser. In Figure 17, 𝑄 is shown as a function of Q for both considered device orientations. For the vertical BHM, the effects of transverse conduction become more and more negligible as the power input to the evaporator increases. For the horizontal orientation, 𝑄 is almost equal to 15% at high power inputs, denoting instead strong conductive effects between adjacent channels in agreement with other literature findings [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF][START_REF] Ayel | Thermal performances of a flat-plate pulsating heat pipe tested with water, aqueous mixtures and surfactants[END_REF][START_REF] Jun | Comparison of the thermal performances and flow characteristics between closed-loop and closed-end micro pulsating heat pipes[END_REF][START_REF] Yang | Performance characteristics of pulsating heat pipes as integral thermal spreaders[END_REF]. In fact, operating a FPPHP in horizontal orientation generally leads to a higher evaporator temperature due to lower thermal performances, eventually increasing the temperature of the fluid flowing towards the condenser. This may result in greater temperature gradients between adjacent branches, leading to higher thermal interactions between the working fluid and the side walls of the channels and, consequently, to stronger transverse interactions. The highlighted variation in terms of transverse conduction from vertical BHM to horizontal orientation may be due to different oscillatory and local heat transfer behaviors of the FPPHP, which are strictly linked to wall-to-fluid heat transfer interactions. To investigate the interplay between the local wall-to-fluid heat flux and transverse conduction qFou,7→8, the heat flux q exchanged in two branches, i.e., channels 7 and 8, was compared with the corresponding heat flux qFou,7→8 (z = 0.02 m) in terms of oscillation frequency. In Figure 18a, all the mentioned heat fluxes related to the vertical BHM are shown for Q = 200 W. In Figure 18b, the power spectra resulting from the wavelet method of Section 3.2 are instead provided, highlighting no dominant frequency indicative of the qFou,7→8 signal. The scalograms that referred to the analyzed signals were not considered, since only the power spectra were believed to be more accurate in underlying differences or similarities between the compared, overall oscillatory behaviors. A similar comparison was carried out in Figure 19 for the horizontal orientation, where all considered signals exhibited the same dominant frequency. In fact, in such an orientation, the fluid similarly oscillates in both FPPHP branches and the regular slug-plug motion results in stronger transverse conduction between adjacent channels. This negative effect is furthermore enhanced by greater amplitudes of q during the horizontal operation, which lead to higher wall temperature variations within the adiabatic section. On the contrary, in the vertical BHM, the greater chaoticity of the fluid motion and different flow pattern, coupled with low-amplitude values of q, seem to prevent the channels from thermally interacting with each other, leading to weaker transverse conduction.

Conclusions

The presented study concluded the analysis on the global and local thermal behavior of a fully metallic flat-plate pulsating heat pipe, first discussed in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF]. The aim of the present work was to design a novel approach for the study of the inner fluid dynamics in FPPHPs without taking advantage of any direct fluid visualizations, which represent most of the experimental studies on flow behavior in such heat transfer systems. The wall-tofluid heat fluxes, evaluated within the adiabatic section of the device by means of thermo- graphic acquisitions on its outer surface, were processed by the wavelet method to estimate dominant fluid oscillation frequencies at varying heat loads and orientations. A validation procedure on synthetic data was carried out to assess the reliability of the estimation procedure. Scalograms (i.e., time-frequency representation of the wavelet transform) were introduced to give useful references for the interpretation of fluid motion in terms of time evolution of the oscillatory phenomena. The characteristic frequencies of the system were furthermore assessed in every device branch by means of power spectra of the considered signals, thus providing a channel-wise description of the fluid motion. The Kmeans clustering method was used for peak identification in scalograms; two quantitative parameters, namely the average continuity time of oscillations tt,av and the average activity time of channels, were consequently proposed for the description of the device oscillatory behavior over time. The average fd,av and standard fd,std channel-wise dominant frequencies were therefore evaluated for the whole device to give an overview on the global oscillatory modes. Finally, thermal interactions between adjacent channels were estimated through the Qratio parameter and explained in terms of fluid oscillation frequencies. The following remarks can be drawn from the present study:



The higher average activity of channels found in the vertical bottom heated mode confirms its better heat transfer capability with respect to the horizontal orientation, even though values of tc,av equal to 23 and 27% of the total observation window at average and high heat loads, respectively, highlight a lower continuity of fluid oscillations in the vertical operation; 

The generally greater continuity of oscillations in the horizontal orientation is also established by the fd,std values, which suggest more similar dominant fluid oscillation frequencies among the device branches especially at high heat loads, where fd,std is set below 0.05 Hz;  For the vertical BHM orientation, fd,av ranges from 0.78 Hz up to 1 Hz during the device's full activation, while the horizontal orientation exhibits values of fd,av that increase from a minimum of about 0.84 Hz to a maximum of 1.25 Hz; 

The heat transferred by conduction between adjacent channels is higher in the horizontal orientation, where Qratio ranges from about 60% at low heat loads to 15% at high heat loads;  Similar oscillatory behaviors in adjacent channels increase the magnitude of transverse thermal interactions.

In conclusion, the presented results quantified the oscillatory behavior of the investigated device. The working regimes assessed in [START_REF] Pagliarini | Novel Infrared Approach for the Evaluation of Thermofluidic Interactions in a Metallic Flat-Plate Pulsating Heat Pipe[END_REF] of the analysis were confirmed by the data, thus giving a deeper insight into the heat transfer modes of the device without adopting any transparent sections for direct fluid visualizations. In view of the results, the adopted technique provided a good insight into the inner fluid dynamics and the overall thermal behavior of the studied device, in accordance with what was observed in previous experimental works relying on different experimental approaches. The here-proposed method for fluid flow assessment can thus be extended to other flat-plate layouts, having walls made entirely of metallic material to give interesting pieces of data for both the augmentation of the literature dataset in terms of fluid oscillation frequencies, and the improvement of existing numerical and theoretical models. To provide a comprehensive summarization of all the results collected in both parts of the experimental work ([29] and present investigation), the main evaluated quantities are listed in Tables 8 and9 for the vertical bottom heated mode and the horizontal orientation, respectively. Here, Req is the evaluated equivalent thermal resistance of the device, cvt and cvs are the coefficients of variation of the wall-to-fluid heat flux over time and along space, respectively, and q80 is the 80th percentile of the heat flux.

Figure 1 .

 1 Figure 1. Rear side of the studied FPPHP (a); sketch of the machined path for the working fluid and reference for the evaporator/condenser/adiabatic sections' locations, as well as for the channels' numeration and area framed by the MWIR camera (green box) (b).

Figure 2 :

 2 Figure 2: Investigated device orientations with respect to the gravity field.

Figure 3 .

 3 Figure 3. Three-dimensional sketch of a FPPHP portion within the adiabatic section; in red, areas of interest for the wall-to-fluid heat flux evaluation.

Figure 4 .

 4 Figure 4. Exact heat flux over time (Amp = 2000 W/m 2 , f = 1.2 Hz), used as input for the wavelet method (a); corresponding scalogram (b) and power spectrum (c).

Figure 5 .

 5 Figure 5. Restored heat flux, used as input for the wavelet method (a); corresponding scalogram (b) and power spectrum (c).

Figure 6 .

 6 Figure 6. Heat flux for a single axial coordinate referred to channel 8 (z = 0.02 m), Q = 100 W, horizontal orientation (a), corresponding scalogram (b) and power spectrum (c).

Figure 7 .

 7 Figure 7. Wall-to-fluid heat flux signals over time related to channel 9 and Q = 250 W for three different axial coordinates and corresponding power spectra; vertical BHM (a,b) and horizontal orientation (c,d).

Figure 8 .

 8 Figure 8. Wavelet scalograms for channels 4 (Ch4), 8 (Ch8) and 12 (Ch12) of the wall-to-fluid heat fluxes for the fixed axial coordinate z = 0.02 m at different heat loads to the evaporator (vertical BHM).

Figure 9 .

 9 Figure 9. Power spectra of the wall-to-fluid heat fluxes for the fixed axial coordinate z = 0.02 m related to every FPPHP channel at different heat loads to the evaporator (vertical BHM).

Figure 10 .

 10 Figure 10. Wavelet scalograms for channels 4 (Ch4), 8 (Ch8) and 12 (Ch12) of the wall-to-fluid heat fluxes for the fixed axial coordinate z = 0.02 m at different heat loads to the evaporator (horizontal orientation).

Figure 11 .

 11 Figure 11. Power spectra of the wall-to-fluid heat fluxes for the fixed axial coordinate z = 0.02 m related to every FPPHP channel at different heat loads to the evaporator (horizontal orientation).

Figure 12 .

 12 Figure 12. Scalogram of Figure 10, Ch 4, having adjusted color scale around the zero value; lowpower areas are present due to residual noise in the input heat flux rather than actual oscillations of the working fluid.

Figure 13 .

 13 Figure 13. Scalograms referred to channel 12 for the vertical BHM (a) and horizontal orientation (b), Q = 250 W; denoised scalograms (c,d) and frequency identification over time (e,f).

Figure 14 .Figure 15 .

 1415 Figure 14. Mean continuity time for average and high heat loads and the two considered orientations.

Figure 16 .

 16 Figure 16. Average dominant fluid oscillation frequency in the overall device (a) and channel-wise standard deviation (b) for the two considered FPPHP orientations.

Figure 17 .

 17 Figure 17. 𝑄 as a function of the power input to the evaporator for both considered device orientations.

Figure 18 .

 18 Figure 18. Wall-to-fluid heat fluxes corresponding to channel 7 and 8, heat flux by transverse conduction qFou,7→8 (a), and corresponding power spectra (b) for the vertical BHM at Q = 200 W and z = 0.02 m.

Figure 19 .

 19 Figure 19. Wall-to-fluid heat fluxes corresponding to channel 7 and 8, heat flux by transverse conduction qFou,7→8 (a), and corresponding power spectra (b) for the horizontal orientation at Q = 200 W and z = 0.02 m.

Table 1 .

 1 Experimental equipment and corresponding uncertainty and operation range.

	Experimental Equipment	Uncertainty	Operating Range
	Power supply	±2 W	50-250 W
	Cryostat	±1 °C	20 °C
	MWIR camera	Absolute: ±1 °C Sensitivity: <25 mK	3-5 μm 5°C-300 °C

Table 3 .

 3 Dominant oscillation frequencies evaluated in the vertical BHM orientation for every channel at varying heat load. Dashes stand for absence of dominant fluid oscillation.

	fd [Hz]				Q [W]
	Channel	50	70	90	120	150	200	250
	1	-	-	0.81	1.15	1.07	1.32	1
	2	1.15	-	0.81	0.58	1.23	1.07	0.66

Table 4 .

 4 Dominant oscillation frequencies evaluated in the horizontal orientation for every channel at varying heat load. Dashes stand for absence of dominant fluid oscillation.

	fd [Hz]			Q [W]
	Channel	50	100	150	200	250
	1	-	0.66	0.87	1.07	1.32
	2	-	0.66	0.93	1.07	1.15
	3	-	0.62	1.87	1.07	1.32
	4	-	0.71	0.93	1	1.23
	5	-	0.71	0.93	1.07	1.32
	6	-	1.42	0.87	1	1.23
	7	-	1.42	0.87	1.07	1.23
	8	-	1.42	0.93	1	1.23
	9	-	0.87	0.81	1.07	1.23
	10	-	0.76	0.87	1.07	1.32
	11	-	0.66	0.87	1.07	1.23
	12	-	0.76	0.87	1	1.23
	13	-	0.71	1.74	1.07	1.23
	14	-	0.71	0.87	1.07	1.23
	15	-	0.66	0.87	1.07	1.32
	16	-	0.66	0.93	1.07	1.15

Table 7 .

 7 Continuity and activity evaluated for the horizontal orientation.

		17.7	52.1	28.6	89.9
	2	14	88.9	30.5	91.65
	3	27.1	88.1	35.6	73.7
	4	24.2	80.1	43.4	89.9
	5	5.7	47.1	31.9	79.9
	6	34.3	73.9	13.3	90.1
	7	9.6	65.1	52.9	85.7
	8	11.6	87.2	25.1	93.8
	9	13.5	74.4	28.8	92.7
	10	5	49.7	32.4	99.7
	11	69.5	100	15.5	94.8
	12	30.2	81.8	10.7	86.6
	13	36	58.4	8.7	80.5
	14	46.1	80.3	21	99
	15	12.5	54	31.4	94.1
	16	13.2	46.8	27.8	84.7
			Q [W]
	Channel	150		250
		tc [%]	ta [%]	tc [%]	ta [%]
	1	46	46

Table 8 .

 8 Evaluated quantities for the vertical bottom heated mode orientation. Green and red colors refer to the intermittent flow and the full activation regime, respectively.

	f	Frequency

Acknowledgments: The work has been carried out through a fruitful and effective collaboration between the University of Parma and PPRIME Institute-ENSMA (Poitiers). The authors would like to acknowledge the European Space Agency (ESA)'s support.

Data Availability Statement:

The data presented in this study are available on request from the corresponding authors. The data are not publicly available due to the large dataset.