Next Article in Journal
Damage Detection for Conveyor Belt Surface Based on Conditional Cycle Generative Adversarial Network
Next Article in Special Issue
Characterization of Biocomposites and Glass Fiber Epoxy Composites Based on Acoustic Emission Signals, Deep Feature Extraction, and Machine Learning
Previous Article in Journal
Demonstration of VOC Fenceline Sensors and Canister Grab Sampling near Chemical Facilities in Louisville, Kentucky
Previous Article in Special Issue
An Assessment of the Effect of Progressive Water Absorption on the Interlaminar Strength of Unidirectional Carbon/Epoxy Composites Using Acoustic Emission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Speed and High-Temperature Calorimetric Solid-State Thermal Mass Flow Sensor for Aerospace Application: A Sensitivity Analysis

by
Lucas Ribeiro
1,*,
Osamu Saotome
1,
Roberto d’Amore
1 and
Roana de Oliveira Hansen
2
1
Instituto Tecnológico de Aeronáutica, Pça. Mal. Eduardo Gomes 50, São José dos Campos 12228-900, Brazil
2
NanoSYD Centre, Mads Clausen Institute, University of Southern Denmark—SDU, Alsion 2, 6400 Sønderborg, Denmark
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(9), 3484; https://doi.org/10.3390/s22093484
Submission received: 16 March 2022 / Revised: 22 April 2022 / Accepted: 30 April 2022 / Published: 3 May 2022
(This article belongs to the Special Issue Damage Detection Systems for Aerospace Applications)

Abstract

:
A high-speed and high-temperature calorimetric solid-state thermal mass flow sensor (TMFS) design was proposed and its sensitivity to temperature and airflow speed were numerically assessed. The sensor operates at 573.15 Kelvin (300 °C), measuring speeds up to 265 m/s, and is customized to be a transducer for an aircraft Air Data System (ADS). The aim was to enhance the system reliability against ice accretion on pitot tubes’ pressure intakes, which causes the system to be inoperative and the aircraft to lose protections that ensure its safe operation. In this paper, the authors assess how the distance between heater and thermal sensors affects the overall TMFS sensitivity and how it can benefit from the inclusion of a thermal barrier between these elements. The results show that, by increasing the distance between the heater and temperature sensors from 0.1 to 0.6 mm, the sensitivity to temperature variation is improved by up to 80%, and that to airspeed variation is improved by up to 100%. In addition, adding a thermal barrier made of Parylene-N improves it even further, by nearly 6 times, for both temperature and air speed variations.

1. Introduction

In an aircraft, critical flight data are provided to the cockpit by the Air Data System (ADS). The consolidated technology used to transduce the external airspeed measurements into digital information for other aircraft systems is the pitot tube, which performs this task by reading the dynamic and static intake air pressures at flight altitude [1,2].
Although the ADS is rigorously designed and qualified to operate in harsh environments, incidents or even accidents may occur that are related to the freezing of the pitot tubes’ orifice, despite its heating capabilities [3]. Once ice crystals accumulate in the transducer, the data become noisy (or lost), leading the air data computer software to ignore the defective probe. In recent decades, several incidents associated with this failure event have been reported [4,5,6]. To address this scenario, the Federal Aviation Administration (FAA) issued Advisory Circular (AC) 91-74B to pilots, which provides instructions about how to identify the atmospheric conditions that lead to this particular risk and how to plan the flight route to avoid supercooled clouds [7]. These instructions impose operational restrictions on the airplane mission, such as changes to a scheduled route, thereby increasing the associated operational costs.
A calorimetric thermal mass flow sensor (TMFS) is designed to transduce the external aircraft airspeed information into electrical signals for an ADS system, and aims to add an extra measurement technique to improve the data voting and reliability when flying in atmospheric conditions that can lead to ice accretion on the pitot’s orifice. The technology related to TMFS has been explored for many aerodynamic applications. These mainly focus on the capability to detect and control the airflow pattern on airfoils with the intention to improve its performance by reducing the drag caused by turbulent flows [8,9,10,11,12,13]. Sturm et al. and Leu et al. present TMFS designs that aim to detect the airflow boundary layer separation from the airfoil while pitching the surface to better control its angle of attack. For the application proposed herein, boundary layer detachment must be avoided so the speed readings do not become noisy. Therefore, the sensor is ideally installed in an articulated airfoil that can ensure the airflow is parallel to the TMFS regardless of the aircraft’s maneuvers.
A review of the literature shows that the maximum measured airflow speed reported by similar TMFS designs is limited to up to 100 m/s, which is much lower than the level necessary for the proposed application, which may be as high as 265 m/s [14,15,16,17,18,19,20,21,22,23]. Additionally, most of these designs have a heater built into a suspended membrane over a cavity to better manage the conductive heat in the sensor body, which is not recommended for the intended application due to high speeds, vibrations, and other mechanical stress factors that may lead to the collapse of the membrane. TMFS designs exist that do not rely on this technique to build the heater element; examples include the design of Glatzl et al., which operates at low temperatures to avoid thermal conductivity issues that can collapse the sensor, or those of Kaltsas et al., Petropoulos et al., Kim et al., Mayer et al., and Bruschi and Piotto, which make use of flow channels to ensure that the fluid flow is laminar, consequently reducing the power demand from the heater once the flow volume is controlled by the tube diameter. However, these are not suitable for the proposed application due to the aforementioned ice condition [24,25,26,27,28,29,30,31].
The TMFS sensitivity is affected by many design constraints, such as the materials employed, the technology used to build the thermal sensors, the sensors’ intrinsic sensitivity to temperature variation, and how well the thermal sensor is isolated from the heater. Regarding TMFS geometry, Mehmood et al. presents a relatively recent review of over 60 different TMFS designs, assesses how the membrane shape affects the sensor sensitivity, and determines the best membrane to heater length ratio (MHR) to be 3.35 for square membranes and 3.30 for circular membranes. Kim et al. evaluates how the distance between heater and temperature sensors affects the sensitivity by increasing the distance in four steps of 4 mm for a sensor built inside a flow channel. The study conclusion is that the TMFS sensitivity worsens when exposed to nitrogen gas as the distance from the heater to the temperature sensor increases [23,27].
In this study, the proposed TMFS calorimetric sensitivity was evaluated by incrementally varying the distance between the heater and the temperature sensors, and by including a thermal barrier made of Parylene-N to better isolate the heater and temperature sensors. The insulator was chosen due to its low thermal conductivity and compatibility with the microfabrication process. In this process, filling of the insulation trench in the SiC substrate can be achieved by the Chemical Vapor Deposition (CVD) process, and the excess is further etched using the Reactive Ion Etching (RIE) process via oxygen gas [28,29]. No change in the heater geometry was made or assessed. The assessment was accomplished through COMSOL Multiphysics and Simulink high fidelity simulations. The former is a Computational Fluid Dynamic (CFD) simulation, which is static and represents a single point of operation, and the second is a dynamic simulation covering an assumed flight envelope, thus allowing the sensitivity to be calculated. The CFD results validate the dynamic simulation when the same boundary conditions apply. No tuning between the models was performed. The results presented herein define the TMFS configuration in which its sensitivity benefits most from the assessed configurations.
The presented calorimetric TMFS design is innovative regarding the manner in which the conductive heat is managed. No similar designs making use of solid materials to isolate sensor elements have been proposed, and no existing designs operate in open circuits that do not rely on flow channels or tubes to take the flow measurements. In this sense, the proposed high-speed and high-temperature solid-state calorimetric TMFS will overcome a necessary gap to become a viable solution for an ADS application.

2. Proposed Calorimetric TMFS Application

The envisioned calorimetric TMFS installation must be located external to the aircraft and on the surface of an articulated airfoil similar to that used for angle of attack sensors. Thus, the non-heated start length is known, and the air flow is maintained to be as stable as possible and is independent of airplane maneuvers. Figure 1 presents an illustration of the calorimetric TMFS installation.
The calorimetric TMFS design is composed of a heater element with two temperature sensors at its opposite ends, so that the air flow crosses these elements perpendicularly. The temperature sensors were designed as Schottky diodes working as differential proportional to absolute temperature (PTAT) sensors. This type of diode has the advantage of a linear voltage–temperature dependence, long-term stability, low turn-on voltage, and faster recovery time when compared with the junction diode [32,33,34]. Additionally, by making use of silicon carbide (SiC) as the substrate to fabricate the Schottky diodes, it is possible to enhance the temperature sensor operation to up to 1000 °C [33,34,35,36,37].
With the objectives of increasing the TMFS sensitivity to heat convection and minimizing the thermal noise due to heat conduction through the TMFS substrate, the idealized calorimetric TMFS was assessed by incrementally increasing the distance between the heater and thermal sensors, and further adding a thermal insulation barrier. Figure 2 presents an illustration of the proposed sensor geometry.
The temperature sensors and the heating element are formed by gold and tungsten, respectively, to take advantage of the metals’ electrical properties. The selected substrate for the TMFS is SiC, which enables the PTAT Schottky temperature sensors to be fabricated.

3. Physics of the Application

A few basic laws are fundamental to the subject of this study: the law of conservation of mass, the law of conservation of momentum, and the law of conservation of energy. The interaction of the fluid flow with the surface over which the fluid flows creates the boundary layer, which has three horizontal flow patterns: laminar, transition, and turbulent; and three vertical flow patterns: viscous sublayer, buffer layer, and turbulent. Above the streamline, the fluid flow is free from any dynamic interference from the surface, and is called the free stream region. Despite its simplicity, the parallel flow over a flat plate occurs in numerous engineering applications, including this work; Figure 3 presents an image of how it develops.
The relative longitudinal position of the calorimetric TMFS heater on a flat plate determines the relation between the heat transferred by convection and the flow speed providing the dimensioning requirements for the heating element, i.e., the amount of power the heater has to generate to allow the speed measurement (heat transfer rate). The conductive heat from the calorimetric TMFS heater to temperature sensors through the substrate affects the overall sensitivity of the sensor, which should ideally measure only the convective heat.

3.1. Heat Convection

A measure of the flow behavior is given by the Reynolds Number, which defines the distance at which the flow becomes turbulent; at this point, more power is demanded to calculate the velocity once the heat convection with the air intensifies. The flow transition point is then derived from the Reynolds Number equation, which represents the ratio of the inertia to viscous forces and is described as [38]:
R e x = ρ u x μ ,
where u is the free stream velocity, x is the point over the flat plate where it is desired to obtain the Reynolds Number, ρ is the fluid density, and μ is the fluid viscosity. The transition from laminar to turbulent flow is obtained from Equation (1) based on the transition distance x c , instead of x ; this defines the critical Reynolds Number, which can vary from approximately 10 5 to 3 × 10 6 , depending on the surface roughness and the turbulence level of the free stream.
To determine the amount of heat transfer needed to measure the desired speed, it is necessary to solve the boundary layer equations for the flat plate, starting with the development at the leading edge ( x ), and then obtaining the temperature profile, considering that there is a non-heated starting length (ξ) until the flow reaches the heating element ( q s ). Figure 4 provides an illustration of this, where δ is the fluid boundary layer and δ t is the thermal boundary layer.
With a known non-heated starting length (ξ), the Nusselt Number, a dimensionless parameter used to characterize the enhancement of the heat transfer due to convection, can be calculated. The Nusselt Number, N u L , is defined as the ratio of the heat transferred from a surface to the heat conducted away by the fluid being defined, for a flat plate with a constant average heat transfer, as [39]:
N u L ¯ = N u L ¯ ξ = 0 L L ξ 1 ξ L p + 1 p + 2 p p + 1 ,
where ξ is the distance of the non-heated part from the leading edge of the flat plate to the heating element; p is a constant that represents the flow pattern, which is equal to 2 for laminar flow and 8 for turbulent flow; and N u L ¯ ξ = 0 is the average Nusselt Number for a plate of length L when heating starts at the leading edge, which is expressed by [39]:
N u L ¯ = 0.680 · R e L 1 2 · P r 1 3 ,
where P r is known as the Prandtl Number, which is defined as the ratio of the kinematic viscosity to the thermal diffusivity and is, therefore, a fluid property.
Then, with Reynolds and Nusselt Numbers in place, it is possible to obtain the needed average heat transfer rate, in W/m2 K, to measure the airflow as [39]:
h L ¯ = N u L ¯ · κ L ,
where κ is the fluid thermal conductivity.
The electrical power needed to be generated by the heated element to measure the intended flow speed is equal to the total heat transfer rate times the heating element area times the temperature difference between the heater and the fluid, and is expressed by [39]:
q c o n v = h L ¯ · A h · T h T f ,
where A h is the heater surface area, T h is the heater temperature, and T f is the fluid temperature.

3.2. Heat Conduction

To examine the heat transfer by conduction, it is necessary to relate it to mechanical, thermal, or geometrical properties so the general form of the heat conduction equation can be expressed by the Fourier’s law [38]:
2 T x 2 + 2 T y 2 + 2 T z 2 + q k = 1 α T t
where k is the material thermal conductivity, q is the volumetric heat generation rate, and α = k ρ c p is known as the material thermal diffusivity composed of ρ , the material density and, c p the material specific heat capacity.
The solution for Fourier’s equation when the TMFS boundary conditions apply can be summarized as a unidimensional problem with the heat being conducted in the medium as a perpendicular vector flowing from the hottest to the cooler parts of the sensor. The implications of this approach are that Equation (6) is reduced to:
α 2 T x 2 + q ρ c p = 0 .
The temperature sensor measurements are then affected by the amount of heat being conducted through the TMFS body and the amount of heat being exchanged by convection. The upstream and downstream temperature sensors sense this phenomenon in different manners when the upstream sensor is measuring cooler air that has not been heated by the TMFS heater, which is not the case for the downstream sensor.

4. Materials and Methods

The results presented herein were obtained by numerical simulations. The COMSOL Multiphysics modeling enables CFD analysis by numerically resolving the thermal conductivity and convection equations based on the TMFS geometry when it interacts with air flowing on its top within fixed boundary conditions. The numerical solution provided by this software is highly reliable and represents, with confidence, the physical phenomena involved; however, it is computationally demanding. The high-fidelity Simulink model was created to dynamically observe the TMFS behavior when operating in the different configurations under evaluation, allowing the sensitivity to be calculated. In addition, the cost of the Simulink computational is much lower than that of the COMSOL Multiphysics simulation. The results obtained by both simulations were compared when operating under the same boundary conditions so to validate the Simulink model.

4.1. COMSOL Multiphysics

The designs presented in Figure 2a,b were modeled using the COMSOL Multiphysics tool, where the junction temperature of each thermal sensor metal was measured after the problem was numerically resolved by the software.
For comparison, 2D and 3D simulations were conducted and showed that the y-direction temperature dependence on the calorimetric TMFS body temperature distribution was limited to a small magnitude. This was because the low cross-section area available for the thermal conductive path formed in a way that could not significantly impact the results obtained from a 2D simulation. Therefore, each proposed design was modeled by creating bidimensional geometries, as presented in Figure 5, to save computational resources, and to allow a solution for heat transfer in a Solids and Fluids Multiphysics simulation to be computed. Additionally, this approach is highly used in CFD simulations because it reduces the size of the meshes generated without compromising the outcome results [40,41]. The TMFS element properties used during COMSOL Multiphysics simulations were as presented at Table 1 [42].
The CFD simulation represents a static solution for a specific moment in the flow pattern under analysis as a picture at an instant in time. This particular moment is represented by the air flow properties, which comprise the boundary conditions for the numeric solution to be computed, as shown in Table 2, with an air speed equal to 0.3 Mach (101.03 m/s) [43].
On top of the object under study, a two-dimensional volume was created that represents the mass of air traveling on the calorimetric TFMS, where the air intake and outtake were the left and right face of the volume, respectively. The non-heated starting length ξ was defined as 1   cm . The mesh was generated in the software by user-controlled options, where the maximum element size was 1.6   mm , the minimum element size was 0.1   μ m , the maximum element growth rate was 1.05, the curvature factor was 0.2, and the resolution of narrow regions was equal to 1. The TMFS geometric dimensions were the same as presented in [44].

4.2. Simulink

The Simulink model was created by capturing all physical and dynamic characteristics of the atmospheric variation while altitude and speed were varied from 1000 to 10,000 m and 101.03 to 300 m/s, respectively, using the “ISA Atmospheric Model” library block. The International Standard Atmosphere (ISA) model, provided by the International Civil Aviation Organization (ICAO), supplies the simulation with air properties such as external temperature, thermal conductivity, specific heat, density, viscosity, and speed of sound at a specific altitude [43].
The calorimetric TMFS interaction with the airflow (convection) and the sensor body (conduction) were modeled using elements from the Simscape Simulink thermal library. Figure 6 depicts the block diagram representing the Simulink model. The airflow parameters are obtained from the “Atmospheric Data” block, which then feed the “Thermal Sensor” blocks, upstream and downstream, where the convective and conductive thermic data are calculated as presented in Section 3. Later, the metal junction temperatures are used to determine the Schottky diode variations, which in turns feed the “PTAT readout circuitry”, where “Vout” is obtained. The “Heater” block implements the 132 Ω thermal resistance dimensioned using Equation (5), which generates the needed electrical power of 1.7 Watts to measure cruise speeds [44]. The dynamic behavior simulated by this model provides evidence about the TMFS sensor operation that supports the proposed air temperature and speed sensitivity analysis.
The speed and altitude inputs were modeled as being constant in the first 0.2 s so the atmospheric/boundary conditions were the same as those used in the COMSOL Multiphysics simulation; then, they were increased during the period from 0.2 to 2 s by varying the speed range up to 300 m/s and the altitude range up to 10,000 m, as previously described.
To better visualize how the Schottky diode junction temperature is simulated in the thermal sensor blocks in Figure 6, the conductive heat in the calorimetric TMFS body is represented by thermal resistances, as in Figure 7, for both assessed cases: without and with a thermal barrier. The model was divided into two layers: one layer representing the SiC wafer having a 0.3 mm thickness, which does not interact with the airflow (heat transfer on solids only), and another layer on top of the SiC wafer where the calorimetric TMFS elements are built in, having a 1 μm thickness, which interacts with the airflow (heat transfer on solids with convection to a fluid).
The Simscape block used to simulate the heat transfer on solids in the first layer of the SiC wafer is the conductive heat transfer, which represents the heat transfer by conduction between two layers of the same material, and is governed by the Fourier law. This remarkable law is used to derive a particular solution from the equations in Section 3.2. Figure 8 presents an image of the conductive heat transfer block and its associated parameters.
For the second layer of the calorimetric TMFS, the heat conduction problem has an amount of energy that is exchanged with the fluid by convection. Simulink Simscape does not have a block that models the conjugate heat transfer by conduction and convection with a fluid. Therefore, to determine the Schottky diode metal junction temperature, Equation (7) was resolved for the unidimensional heat conduction with heat transfer by convection to a fluid, and its solution is as follows:
T s = T q + h d / k T f h d / k + 1 ,
where T s is the Schottky diode metal junction temperature; T q is the conductive temperature from the SiC wafer; h is the heat transfer rate calculated using Equation (4), which varies with airflow speed (due to Reynolds Number influence) and the distance L from the flat plate leading edge; d is the Schottky metal thickness; k is the fluid thermal conductivity; and T f is the fluid temperature.
The differential PTAT circuit is depicted in Figure 9, where D1 and D2 are the upstream and downstream Schottky diodes, respectively, I1 and I2 are current references to polarize the diodes, Vin is equal to 12 Vdc, and A1 is an amplifier with gain equal to 12.
During the simulation, the heater was kept at a constant temperature of 573.15 Kelvin (300 °C). The Schottky diode Simulink model was created using the design proposed by Krasnov et al. [45], in which a 4H-SiC extreme temperature sensor is depicted. With the intention to validate the created Schottky diode model, the voltage–temperature dependence curves for electrical currents varying from 1   pA to 0.1   mA obtained from the Simulink model are presented in Figure 10, in addition to their respective sensitivities in mV / K .
All Simulink simulations, from the Schottky diode validation above to the TMFS dynamic behavior presented in Figure 6, were executed using a fixed-step of 1   μ s and automatic solver selection.

4.3. Calorimetric TMFS Fabrication Process

The calorimetric TMFS thermal barrier can be fabricated by making use of standard cleanroom microfabrication techniques. The proposed process is only a coarse description of an actual fabrication process flow that is, in practice, much more complicated. This process can start with a commercial SiC wafer (300 µm thick). A photoresist is spin-coated onto the wafer, following a UV exposure using a photomask with the thermal barrier area design. The area for SiC etching is defined by the photolithography process, and a trench can be etched on SiC via Inductive Coupled Plasma RIE, as described in [45,46]. Then, a new photolithography process is conducted to define the areas for Parylene-N deposition, which can then be performed via CVD. The photoresist is then removed in an acetone bath, and the SiC wafer containing the Parylene-N trenches is cleaned via RCA cleaning.
The Schottky diode electrodes and heater can be defined by a subsequent photolithography process, followed by tungsten and gold deposition via e-beam metal deposition. Resist removal and lift-off are performed with an acetone bath. The final step is the dicing of the sensors with a diamond dicing saw and wiring-bonding to a PCB.

5. Results

Both COMSOL Multiphysics (static) and Simulink (dynamic) simulations were executed by incrementally increasing the distance between the heater and the sensing elements, i.e., distance “b” as defined in Figure 2a, by steps of 0.1   mm from the initial value of 0.1   mm to 0.6   mm . When reaching this value, the distance was held constant, but an insulation trench was added. This was filled with Parylene-N having an initial thickness, “ITh”, as defined in Figure 2b, of 0.2   mm ; this was incrementally increased by 0.1   mm to 0.4   mm . The results obtained from these simulations are presented hereafter.

5.1. COMSOL Multiphysics

By incrementing the distance between the heater and the temperature sensors, and later adding the insulation trench made with Parylene-N, the mesh was automatically generated in the tool as stated in Section 4.1; the quantity of the nodes and the minimum element quality of the domain and boundary elements are presented in Table 3.
The delta temperature Δ T between the upstream and downstream Schottky diodes’ metals was measured and is presented in Figure 11, in Kelvin.
Figure 12 presents an image of the COMSOL Multiphysics solution computed for each condition assessed, where, from letter (a) to (f), the distance between the heater and thermal sensors is increased from 0.1 to 0.6 mm, and a small variation in the shade of red can be seen; from letter (g) to (i), the thickness of the thermal barrier is varied from 0.2 to 0.4 mm, and a shade of yellow starts to develop, showing a lower temperature on the thermal sensor region and indicating the heat conduction was partially blocked through the TMFS substrate.

5.2. Simulink

As mentioned previously, in the first 0.2 s of the Simulink simulation, the airflow boundary conditions were kept at the same values as those used in the COMSOL Multiphysics CFD simulation for validation and comparison purposes. After this time, speed and altitude were linearly increased to 300 m/s and 10,000 m respectively. The obtained Δ T in Kelvin for each condition is presented in Figure 13.
The Δ T results obtained from COMSOL Multiphysics and Simulink simulations are then compared in Table 4. It is possible to notice differences due to the solver configurations of both software packages, which led to numerical discrepancies, and due to the thermal inertia captured in the Simulink model for the cases where the thermal barrier was present.
To calculate the overall TMFS sensitivity to temperature variation from the Simulink model, Vout at the PTAT circuit was linearized using the “polyfit” and “polyval” functions for t > 0.2, for each condition assessed. The calorimetric TMFS sensitivity values to temperature variation are presented in Figure 14, where the linearized voltage curves and their respective slopes, or sensitivity, are highlighted. It is possible to notice the improvement as the distance between the heater and thermal sensors increases (foreground graphic), and then with the addition of the thermal barrier (background graphic).
To calculate the overall TMFS sensitivity to air speed variation from the Simulink model, Vout at the PTAT circuit was linearized using the “polyfit” and “polyval” functions for t > 0.2, for each condition assessed. The calorimetric TMFS sensitivity to air speed variation is presented in Figure 15, where the linearized voltage curves and their respective slopes, or sensitivity, are highlighted. It is possible to notice the improvement as the distance between the heater and thermal sensors increases (foreground graphic), and then with the addition of the thermal barrier (background graphic).

6. Discussion

The simulation results presented in the previous section show a consistent improvement in the sensitivity to both temperature and speed variations of the proposed high-speed and high-temperature calorimetric solid-state TMFS when varying the distance between the heater and thermal sensors, and then including Parylene-N as a thermal barrier between them.
The results obtained via COMSOL Multiphysics CFD simulation presented in Figure 12 show how the heat propagates through the sensor body. As a result, it is possible to visualize the improvements in the Schottky diode metal junction temperature, which benefits the temperature difference between the upstream and downstream thermal sensors, as presented in Figure 11.
Due to the computational costs associated with COMSOL Multiphysics CFD simulation, which resolves the heat propagation problem by iteratively solving the heat transport from one mesh node to another until the convergence criteria are established by the software, the Simulink model was created. The solutions from both software packages were further compared, as presented in Table 4, for the same boundary conditions. The results were found to be similar, with an acceptable error margin, when the thermal barrier was included, validating the Simulink results for this condition. For the condition where there was no thermal barrier, the results from Simulink showed a significant error for b = 0.1 mm and b = 0.2 mm. This was caused by a numerical error due to the short distance “D” on the conductive heat blocks, which appears to be a limitation of the Simscape block. The sensor sensitivity to air temperature and air speed variation was then calculated from dynamic results obtained in Simulink, and showed significant improvements in the overall sensor sensitivity, as presented in Figure 14 and Figure 15. Furthermore, the addition of Parylene-N as a thermal barrier was proven to be highly efficient at reducing the heat conduction through the TMFS body, thus resolving the issue reported in previous work [44].
Here, it is important to comment that the results presented by Kim et al. [27] show the opposite behavior, where the sensitivity decreases as the distance between the temperature sensors and the heater increases. This is explained by the fact that the authors used a flow channel, which requires lower power from the heater (0.8 Watts) than the application proposed herein (1.7 Watts), and nitrogen gas, which is lighter than air, to perform the tests that support this conclusion. The heat convection with the fluid dissipated the energy from the heater operating at 160 °C; as a result, the more distant temperature sensors were unaffected, thus decreasing the overall sensor sensitivity. The application proposed herein does not rely on flow channels, operates at 573.15 Kelvin (300 °C), and has different boundary conditions. Therefore, the results presented by Kim et al. [27] do not form a basis for comparison with the presented results. Furthermore, previous work has shown that the thermal plume generated by the heater while operating at cruise speeds is prominent. This allows the temperature sensors placed at the distances used in the current study to measure the convective effect resulting from the airflow interaction with the proposed sensor [44].
Based on the results presented herein, the final TMFS dimensions, using a thermal barrier of 4 mm filled with Parylene-N and a distance between the heater and thermal sensors of 6 mm, can be defined, thus allowing it to be manufactured and tested in wind tunnels for the intended operation.

7. Conclusions

From the COMSOL Multiphysics and Simulink simulation results, it is possible to observe that the delta temperature Δ T measured between thermal sensors increases as the distance between the heater and the thermal sensors increases for airflow speeds and altitude variations representing a generic aircraft flight envelope. By subsequently adding a thermal barrier made of Parylene-N, the Δ T measured between thermal sensors increased even further, showing it to provide good thermal insulation and fit for its purpose.
The results obtained from COMSOL Multiphysics validate the Simulink simulation when compared under the same boundary conditions, with the exception of b = 0.1 mm and b = 0.2 mm, for the aforementioned reasons. The Δ T between the temperature sensors increases as the air temperature decreases, showing that the sensor operation at cruise speeds can yield a better flow measurement due to the high thermal convection caused by the air flow dynamics. It is important to also note that the thermal barrier has a thermic delay. Further investigation will be required about how this behaves when operating during long and continuous periods of time to ensure the benefit observed is permanent.
Furthermore, the improvements observed in the Δ T measured between the thermal sensors for the different simulated conditions show that the calorimetric TMFS sensitivity to air temperature and air speed variations was significantly improved.

Author Contributions

Data curation, L.R.; Formal analysis, L.R., O.S., R.d. and R.d.O.H.; Funding acquisition, R.d.O.H.; Investigation, L.R.; Methodology, L.R.; Resources, O.S., R.d. and R.d.O.H.; Software, L.R.; Supervision, O.S., R.d. and R.d.O.H.; Validation, L.R.; Visualization, L.R.; Writing–original draft, L.R.; Writing–review & editing, O.S., R.d. and R.d.O.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Independent Research Fund Denmark grant number 0217-00023B.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Models created for simulations herein presented can be shared upon request to authors as per the contact information provided above.

Acknowledgments

The authors would like to acknowledge the University of Southern Denmark (SDU) Sønderborg for the financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pratt, R.W. Flight Control Systems: Practical Issues in Design and Implementation; The Institution of Electrical Engineers and The American Institute of Aeronautics and Astronautics: Herts, UK, 2000. [Google Scholar]
  2. Nelson, R.C. Flight Stability and Automatic Control; McGraw-Hill Book Company: New York, NY, USA, 1989; ISBN 0-07-046218-6. [Google Scholar]
  3. Ali, A.; Souza, J.R.B.; Lisboa, K.M.; Cotta, R.M. Thermal analysis of heated pitot probes in atmospheric conditions of ice accretion. In Proceedings of the 23rd ABCM International Congress of Mechanical Engineering, COBEM, Rio de Janeiro, Brazil, 6–11 December 2015. [Google Scholar] [CrossRef]
  4. Bureau d’Enquêtes Et d’Analyses Pour La Sécurité De I’aviationcivile, Interim Report No 2 on the Accident on 1st June 2009 to the Airbus A330-203 Registered F-gzcp Operated by Air France Flight AF447 Rio de Janeiro-Paris. 2009. Available online: https://bea.aero/docspa/2009/f-cp090601.en/pdf/f-cp090601.en.pdf (accessed on 12 July 2018).
  5. Australian Transport Safety Bureau, In-Flight upset 154 Km West of Learmonth, WA 7 October 2008, VH-QPA, Airbus A330-303. 2008. Available online: https://www.atsb.gov.au/publications/investigation_reports/2008/aair/ao-2008-070.aspx (accessed on 12 July 2018).
  6. United Arab Emirates General Civil Aviation Authority, Serious Incident, Unreliable Airspeed Indication, Etihad Airways, A6-EHF, Airbus A340-600. 2013. Available online: https://www.gcaa.gov.ae/en/epublication/admin/iradmin/lists/incidents%20investigation%20reports/attachments/67/2013-2013%20-%20final%20report%20-%20aais%20case%20aifn-0005-2013%20-%20a6-ehf%20serious%20incident.pdf (accessed on 12 July 2018).
  7. Barbagallo, J. Pilot guide: Flight in ice condition. In FAA Advisory Circular Repository; [AC 91-74B]; Federal Aviation Administration: Washington, DC, USA, 2015. Available online: https://www.faa.gov/documentlibrary/media/advisory_circular/ac_91-74b.pdf (accessed on 23 August 2019).
  8. Tsao, T.; Jiang, F.; Miller, R.; Tai, Y.; Gupta, B.; Goodman, R.; Tung, S.; Ho, C. An integrated MEMS system for turbulent boundary layer control. In Proceedings of the International Solid State Sensors and Actuators Conference, Chicago, IL, USA, 19 June 1997; Volume 1, pp. 315–318. [Google Scholar] [CrossRef] [Green Version]
  9. Tung, S.; Maines, B.; Jiang, F.; Tsao, T. Development of a MEMS-Based Control System for Compressible Flow Separation. J. Microelectromechanical Syst. 2004, 13, 91–99. [Google Scholar] [CrossRef]
  10. Kumar, S.; Reynolds, W.; Kenny, T. MEMS based transducers for boundary layer control, Technical Digest. In Proceedings of the IEEE International MEMS 99 Conference, Orlando, FL, USA, 21 January 1999. [Google Scholar] [CrossRef]
  11. Sturm, H.; Dumstorff, G.; Busche, P.; Westermann, D.; Lang, W. Boundary layer separation and reattachment detection on airfoils by thermal flow sensors. Sensors 2012, 12, 14292–14306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Mansoor, M.; Haneef, I.; Akhtar, S.; Rafiq, M.A.; Ali, S.Z.; Udrea, F. SOI CMOS multi-sensors MEMS chip for aerospace applications. In Proceedings of the IEEE Conference on Sensors, Valencia, Spain, 2–5 November 2014. [Google Scholar] [CrossRef]
  13. Leu, T.S.; Yu, J.M.; Miau, J.J.; Chen, S.J. MEMS flexible thermal flow sensors for measurement of unsteady flow above a pitching wind turbine blade. Exp. Therm. Fluid Sci. 2016, 77, 167–178. [Google Scholar] [CrossRef]
  14. Balakrishnan, V.; Phan, H.P.; Dinh, T.; Dao, D.V.; Nguyen, N.T. Thermal flow sensors for harsh environments. Sensors 2017, 17, 2061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Que, R.; Zhu, R. Aircraft aerodynamic parameter detection using micro hot-film flow sensor array and BP neural network identification. Sensors 2012, 12, 10920–10929. [Google Scholar] [CrossRef] [Green Version]
  16. Lyons, C.; Friedberger, A.; Welser, W.; Muller, G.; Krotz, G.; Kassing, R. A high-speed mass flow sensor with heated silicon carbide bridges. In Proceedings of the MEMS 98, IEEE, Eleventh Annual International Workshop on Micro Electro Mechanical Systems. An Investigation of Micro Structures, Sensors, Actuators, Machines and Systems (Cat. No.98CH36176), Heidelberg, Germany, 25–29 January 1998; pp. 356–360. [Google Scholar] [CrossRef]
  17. Sun, J.; Qin, M.; Huang, Q. Flip-chip packaging for a two-dimensional thermal flow sensor using a copper pillar bump technology. IEEE Sens. J. 2007, 7, 990–995. [Google Scholar] [CrossRef]
  18. Kuo, J.T.W.; Yu, L.; Meng, E. Micromachined thermal flow sensors—A Review. Micromachines 2012, 3, 550–573. [Google Scholar] [CrossRef] [Green Version]
  19. Billat, S.; Storz, M.; Ashauer, H.; Hedrich, F.; Kattinger, G.; Lust, L.; Ashauer, M.; Zengerle, R. Thermal flow sensors for harsh environment applications. Procedia Chem. 2009, 1, 1459–1462. [Google Scholar] [CrossRef] [Green Version]
  20. Lange, P.; Weiss, M.; Warnat, S. Characteristics of a micro-mechanical thermal flow sensor based on a two hot wires principle with constant temperature operation in a small channel. J. Micromech. Microeng. 2014, 24, 125005. [Google Scholar] [CrossRef]
  21. Jiang, F.; Lee, G.B.; Tai, Y.C.; Ho, C.M. A flexible micromachine-based shear-stress sensor array and its application to separation-point detection. Sens. Actuators A Phys. 2000, 79, 194–203. [Google Scholar] [CrossRef]
  22. Cerimovic, S.; Talic, A.; Beigelbeck, R.; Sauter, T.; Kohl, F.; Schalko, J.; Keplinger, F. A novel calorimetric flow sensor implementation based on thermal sigma-delta modulation. In Proceedings of the SENSORS, 2009 IEEE, Christchurch, New Zealand, 25–28 October 2009; pp. 1923–1926. [Google Scholar] [CrossRef]
  23. Mehmood, Z.; Haneef, I.; Ali, S.Z.; Udrea, F. Sensitivity enhancement of silicon-on-insulator CMOS MEMS thermal hot-film flow sensors by minimizing membrane conductive heat losses. Sensors 2019, 19, 1860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Glatzl, T.; Beigelbeck, R.; Cerimovic, S.; Steiner, H.; Wenig, F.; Sauter, T.; Treytl, A.; Keplinger, F. A thermal flow sensor based on printed circuit technology in constant temperature mode for various fluids. Sensors 2019, 19, 1065. [Google Scholar] [CrossRef] [Green Version]
  25. Kaltsas, G.; Nassiopoulos, A.A.; Nassiopoulou, A.G. Characterization of a silicon thermal gas-flow sensor with porous silicon thermal isolation. IEEE Sens. J. 2002, 2, 463–475. [Google Scholar] [CrossRef]
  26. Petropoulos, A.; Kaltsas, G. Study and evaluation of a PCB-MEMS liquid microflow sensor. Sensors 2010, 10, 8981–9001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kim, T.H.; Dong-Kwon, K.; Sung, J.K. Study of the sensitivity of a thermal flow sensor. Int. J. Heat Mass Transf. 2009, 52, 7–8. [Google Scholar] [CrossRef]
  28. Wang, Y.; Dai, W.; Wu, Y.; Guo, Y.; Xu, R.; Yan, B.; Xu, Y. A bendable microwave GaN HEMT on CVD Parylene-C substrate. In Proceedings of the IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization (NEMO), Hangzhou, China, 7–9 December 2020. [Google Scholar]
  29. Li, L.; Mason, A.J. Post-CMOS parylene packaging for on-chip biosensor arrays. In Proceedings of the IEEE Sensors Conference, Waikoloa, HI, USA, 1–4 November 2010. [Google Scholar]
  30. Mayer, F.; Paul, O.; Baltes, H. Influence of Design Geometry and Packaging on the response of thermal CMOS flow sensors. In Proceedings of the International Solid-State Sensors and Actuators Conference—TRANSDUCERS, Stockholm, Sweden, 25–29 June 1995. [Google Scholar]
  31. Bruschi, P.; Piotto, M. Design issues for low power integrated thermal flow sensors with ultra-wide dynamic range and low insertion loss. Micromachines 2012, 3, 295–314. [Google Scholar] [CrossRef] [Green Version]
  32. Hoeneisen, B.; Mead, C.A. Power Schottky diode design and comparison with the junction diode. Solid-State Electron. 1971, 14, 1225–1236. [Google Scholar] [CrossRef]
  33. Pascu, R.; Pristavu, G.; Brezeanu, G.; Draghici, F.; Godignon, P.; Romanitan, C.; Matei, S.; Adrian, T. 60–700 K CTAT and PTAT temperature sensors with 4H-SiC Schottky diodes. Sensors 2021, 21, 942. [Google Scholar] [CrossRef]
  34. Mitra, J.; Feng, L.; Penate, L.; Dawson, P. An alternative methodology in Schottky diode physics. J. Appl. Phys. 2015, 117, 244501. [Google Scholar] [CrossRef]
  35. Josan, I.; Boianceanu, C.; Brezeanu, G.; Obreja, V.; Avram, M.; Puscasu, D.; Ioncea, A. Extreme environment temperature sensor based on silicon carbide Schottky Diode. In Proceedings of the International Semiconductor Conference, Sinaia, Romania, 12–14 October 2009. [Google Scholar] [CrossRef]
  36. Zhang, N.; Pisano, A. Harsh environment temperature sensor based on 4H-Silicon carbide PN diode. In Proceedings of the International Conference on Solid State Sensors and Actuators (TRANSDUCERS), Barcelona, Spain, 16–20 June 2013. [Google Scholar] [CrossRef]
  37. Rao, S.; Pangallo, G.; Pezzimenti, F.; Della Corte, F. High-performance temperature sensor based on 4H-SiC Schottky diodes. IEEE Electron Device Lett. 2015, 36, 720–722. [Google Scholar] [CrossRef]
  38. White, F.M. Viscous Fluid Flow, 3rd ed.; McGraw-Hill Book Company: New York, NY, USA, 2006; ISBN 007-124493-X. [Google Scholar]
  39. Bergman, T.L.; Lavine, A.S.; Incropera, F.P.; Dewitt, D.P. Fundamentals of Heat and Mass Transfer, 7th ed.; John Willey & Sons: Hoboken, NJ, USA, 2011; ISBN 13 978-0470-50197-9. [Google Scholar]
  40. Alam, J.M.; Walsh, R.P.; Hossain, M.A.; Rose, A.M. A computational methodology for two-dimensional fluid flows. Int. J. Numer. Meth. Fluids 2014, 75, 835–859. [Google Scholar] [CrossRef]
  41. Alizard, F.; Robinet, J.C. Modeling of optimal perturbations in flat plate boundary layer using global modes: Benefits and limits. Theory Comput. Fluid Dyn. 2011, 25, 147–165. [Google Scholar] [CrossRef]
  42. Cardarelli, F. Materials Handbook: A Concise Desktop Reference; Springer: Berlin/Heidelberg, Germany, 2005; ISBN 978-1-84628-668-1. [Google Scholar]
  43. ISA. Manual of the ICAO Standard Atmosphere. International Civil Aviation Organization. 1993. Available online: http://code7700.com/pdfs/icao/icao_doc_7488_standard_atmosphere.pdf (accessed on 9 May 2017).
  44. Ribeiro, L.C.; Korsa, M.T.; Saotome, O.; Adam, J.; Hansen, R.M.O. Design of a micro-machined flow sensor for aircraft air data systems application: Mechanical considerations. In Proceedings of the International Conference on Microelectronics and Microsystems Technology ICMMT, New York, NY, USA, 10–11 December 2020. [Google Scholar]
  45. Krasnov, V.A.; Shutov, S.V.; Yerochin, S.Y.; Demenskyi, O.M. High Temperature operation limit assessment for 4H-SiC schottky diode-based extreme temperature sensors. IEEE Sens. J. 2019, 19, 1640–1644. [Google Scholar] [CrossRef]
  46. Racka-Szmidt, K.; Stonio, B.; Zelazko, J.; Filipiak, M.; Sochacki, M. A review: Inductively coupled plasma reactive ion etching of silicon carbide. Materials 2022, 15, 123. [Google Scholar] [CrossRef]
Figure 1. Idealized TMFS installation.
Figure 1. Idealized TMFS installation.
Sensors 22 03484 g001
Figure 2. Top view of the idealized TMFS layout (a) without and (b) with thermal barrier, where the red square “A” represents the heated area, the yellow rectangles are the temperature sensor metal, “b” is the distance between the heater and temperature sensors, and “ITh” represents the insulation barrier thickness.
Figure 2. Top view of the idealized TMFS layout (a) without and (b) with thermal barrier, where the red square “A” represents the heated area, the yellow rectangles are the temperature sensor metal, “b” is the distance between the heater and temperature sensors, and “ITh” represents the insulation barrier thickness.
Sensors 22 03484 g002
Figure 3. Velocity boundary layer development on a flat plate.
Figure 3. Velocity boundary layer development on a flat plate.
Sensors 22 03484 g003
Figure 4. Thermal boundary condition for on a flat plate with unheated starting length.
Figure 4. Thermal boundary condition for on a flat plate with unheated starting length.
Sensors 22 03484 g004
Figure 5. TMFS geometry with thermal barriers for COMSOL Multiphysics simulation, showing the temperature sensors and heater. Figure is not to scale.
Figure 5. TMFS geometry with thermal barriers for COMSOL Multiphysics simulation, showing the temperature sensors and heater. Figure is not to scale.
Sensors 22 03484 g005
Figure 6. Simulink model block diagram.
Figure 6. Simulink model block diagram.
Sensors 22 03484 g006
Figure 7. Side view of idealized TMFS layout with thermal resistances and heat flux for (a) without and (b) with thermal barrier.
Figure 7. Side view of idealized TMFS layout with thermal resistances and heat flux for (a) without and (b) with thermal barrier.
Sensors 22 03484 g007aSensors 22 03484 g007b
Figure 8. Conductive heat Simulink Simscape block implementation by Mathworks.
Figure 8. Conductive heat Simulink Simscape block implementation by Mathworks.
Sensors 22 03484 g008
Figure 9. TMFS differential PTAT circuit.
Figure 9. TMFS differential PTAT circuit.
Sensors 22 03484 g009
Figure 10. Schottky diode Simulink model validation: IRef in amps and the correspondent Schottky diode sensitivity in mV/K.
Figure 10. Schottky diode Simulink model validation: IRef in amps and the correspondent Schottky diode sensitivity in mV/K.
Sensors 22 03484 g010
Figure 11. Calorimetric TMFS upstream and downstream temperature difference for each assessed condition, where the dashed line separates the increase in “b” and the introduction and variation in “ITh” thickness for a fixed “b” of 0.6 mm—COMSOL Multiphysics.
Figure 11. Calorimetric TMFS upstream and downstream temperature difference for each assessed condition, where the dashed line separates the increase in “b” and the introduction and variation in “ITh” thickness for a fixed “b” of 0.6 mm—COMSOL Multiphysics.
Sensors 22 03484 g011
Figure 12. CFD results: (a) b = 0.1   mm , (b) b = 0.2   mm , (c) b = 0.3   mm , (d) b = 0.4   mm , (e) b = 0.5   mm , (f) 0.6   mm , (g) b = 0.6   mm   |   I T h = 0.2   mm , (h) b = 0.6   mm   |   I T h = 0.3   mm , and (i) b = 0.6   mm   |   I T h = 0.4   mm —COMSOL Multiphysics.
Figure 12. CFD results: (a) b = 0.1   mm , (b) b = 0.2   mm , (c) b = 0.3   mm , (d) b = 0.4   mm , (e) b = 0.5   mm , (f) 0.6   mm , (g) b = 0.6   mm   |   I T h = 0.2   mm , (h) b = 0.6   mm   |   I T h = 0.3   mm , and (i) b = 0.6   mm   |   I T h = 0.4   mm —COMSOL Multiphysics.
Sensors 22 03484 g012aSensors 22 03484 g012b
Figure 13. Measured Δ T between upstream and downstream sensors. The plot in the background presents the Δ T results with the thermal barrier and the plot in the foreground presents the Δ T results without the thermal barrier. The arrows point to the temperature difference measured at 0.2 s of the simulation—Simulink.
Figure 13. Measured Δ T between upstream and downstream sensors. The plot in the background presents the Δ T results with the thermal barrier and the plot in the foreground presents the Δ T results without the thermal barrier. The arrows point to the temperature difference measured at 0.2 s of the simulation—Simulink.
Sensors 22 03484 g013
Figure 14. Linearized Vout from PTAT circuit as a function of air temperature and respective sensitivities to temperature variation. The plot in the background presents the TMFS Vout and sensitivity results with the thermal barrier and the plot in the foreground presents the TMFS Vout and sensitivity results without the thermal barrier. The arrows point to the sensitivity to air temperature variation of each TMFS design assessed during the simulation—Simulink.
Figure 14. Linearized Vout from PTAT circuit as a function of air temperature and respective sensitivities to temperature variation. The plot in the background presents the TMFS Vout and sensitivity results with the thermal barrier and the plot in the foreground presents the TMFS Vout and sensitivity results without the thermal barrier. The arrows point to the sensitivity to air temperature variation of each TMFS design assessed during the simulation—Simulink.
Sensors 22 03484 g014
Figure 15. Linearized Vout from the PTAT circuit as a function of air speed and respective sensitivities to speed variation. The plot in the background presents the TMFS Vout and sensitivity results with the thermal barrier and the plot in the foreground presents the TMFS Vout and sensitivity results without the thermal barrier. The arrows point to the sensitivity to air speed variation of each TMFS design assessed during simulation—Simulink.
Figure 15. Linearized Vout from the PTAT circuit as a function of air speed and respective sensitivities to speed variation. The plot in the background presents the TMFS Vout and sensitivity results with the thermal barrier and the plot in the foreground presents the TMFS Vout and sensitivity results without the thermal barrier. The arrows point to the sensitivity to air speed variation of each TMFS design assessed during simulation—Simulink.
Sensors 22 03484 g015
Table 1. Material properties.
Table 1. Material properties.
PropertyUnitSubstrate (SiC)Insulation (Parylene N)Heater (W)Thermal Sensor (Au)
Electrical ResistivityΩm @ 20 °C10-5.65 × 10−82.44 × 10−8
Thermal ConductivityW/m—K1140.126174300
Young’s ModuleNmm24.15 × 105-4.11 × 10579 × 103
Poisson’s RatioN/A0.16-0.2800.42
Densitygcm−3 @ 25 °C3.161.1119.319.32
Melting Point°C2830-36871064
Specific Heat CapacityJ/g-°C0.6701.30.1340.133
Table 2. Air properties as per ISA standard.
Table 2. Air properties as per ISA standard.
PropertyUnitValue
Dynamic ViscosityμPa.s17.606
Sound Speedm/s336.766
Kinematic Viscosityμm2/s15.705
Thermal ConductivityW/m·K0.024874
Altitudem1000
TemperatureK282.207
PressurehPa90,813
Table 3. Mesh node numbers and quality.
Table 3. Mesh node numbers and quality.
b / I T h   in   mm 0.10.20.30.40.50.60.6/0.20.6/0.30.6/0.4
Domain1,427,3291,464,0801,500,4541,535,4961,571,8201,607,1781,602,0221,602,8761,601,974
Boundary29,37829,96730,54931,12631,69832,26632,52032,52632,540
Minimum Quality0.48750.44850.44910.43860.47480.44840.44840.44840.4484
Table 4. Comparison of COMSOL Multiphysics and Simulink Δ T results.
Table 4. Comparison of COMSOL Multiphysics and Simulink Δ T results.
B   in   mm IT h in   mm Δ T from   COMSOL   Multiphysics   in   Kelvin Δ T from   Simulink   in   Kelvin Error in Kelvin (%)
0.101.12.61.5 (58%)
0.202.323.51.18 (34%)
0.303.694.40.71 (16%)
0.405.185.40.22 (4%)
0.506.796.40.39 (6%)
0.608.487.50.98 (13%)
0.60.185.1287.52.38 (3%)
0.60.2107.89114.97.01 (6%)
0.60.3126.541358.46 (6%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ribeiro, L.; Saotome, O.; d’Amore, R.; de Oliveira Hansen, R. High-Speed and High-Temperature Calorimetric Solid-State Thermal Mass Flow Sensor for Aerospace Application: A Sensitivity Analysis. Sensors 2022, 22, 3484. https://doi.org/10.3390/s22093484

AMA Style

Ribeiro L, Saotome O, d’Amore R, de Oliveira Hansen R. High-Speed and High-Temperature Calorimetric Solid-State Thermal Mass Flow Sensor for Aerospace Application: A Sensitivity Analysis. Sensors. 2022; 22(9):3484. https://doi.org/10.3390/s22093484

Chicago/Turabian Style

Ribeiro, Lucas, Osamu Saotome, Roberto d’Amore, and Roana de Oliveira Hansen. 2022. "High-Speed and High-Temperature Calorimetric Solid-State Thermal Mass Flow Sensor for Aerospace Application: A Sensitivity Analysis" Sensors 22, no. 9: 3484. https://doi.org/10.3390/s22093484

APA Style

Ribeiro, L., Saotome, O., d’Amore, R., & de Oliveira Hansen, R. (2022). High-Speed and High-Temperature Calorimetric Solid-State Thermal Mass Flow Sensor for Aerospace Application: A Sensitivity Analysis. Sensors, 22(9), 3484. https://doi.org/10.3390/s22093484

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop