The investigation of super-cooled droplet impingement characteristics is the most important step for aircraft icing and anti-icing/de-icing analyses. The Lagrangian method and the Eulerian method are widely used to compute the droplet motion and collection efficiency, and the two methods are considered to obtain almost the same results for surface impingement characteristics under icing conditions. The models and implementation approaches of the two methods were established in this work, and the simulations of droplet motion were carried out for a NACA 0012 airfoil, a 2D section of an A320 head, a multi-element airfoil, and an icing wind tunnel. The collection efficiencies of the NACA 0012 airfoil obtained by the present Lagrangian and Eulerian methods show good agreement with the results in the literature, validating the established methods. The droplet impingement characteristics of the two methods are consistent for the aircraft surfaces without upstream trajectory deflections. However, when the droplet motion is deflected by the frontal body before hitting the rear surfaces, the results obtained by the two methods are different whether the droplet trajectories intersect or not, which subverts the traditional opinion that the Lagrangian and Eulerian methods would obtain the same result of the droplet impingement characteristics. The reason is studied in detail according to the droplet motion results in the icing wind tunnel. The findings of this work are helpful for the accuracy of aircraft icing and anti-icing/de-icing simulations, and useful for the development of airworthiness certification.

The impact of supercooled water droplets is the cause of aircraft icing, and the acquisition of impingement characteristics is the key to icing prediction and the design of ice protection systems. The introduction of water droplet collection efficiency is required to obtain the characteristics for the Lagrangian method. In this work, a traditional flow tube method, a local flow tube method, and a statistical method are established to calculate the local collection efficiency, based on Lagrangian droplet trajectories. Through the numerical simulations of the air–droplet flow field around an NACA 0012 airfoil, the accuracies of the three methods in regard to collection efficiency are verified. Then, these three methods are applied to obtain the results for water droplet trajectories and the collection efficiency of an S-shaped duct, a 2D engine cone section and an icing wind tunnel. The results show that the distributions of water droplet collection efficiency obtained by the three methods are consistent and the three methods are all feasible when the water droplets do not overlap or cross before hitting the aircraft surfaces. When the water droplets are shadowed by upstream surfaces or blown by air injection, the droplet trajectories might overlap or even cross, and the local collection efficiencies obtained by the traditional flow tube method, local flow tube method, and statistical method might differ. The statistical method is relatively accurate. However, not all the droplet impingement characteristics obtained by the three methods are different due to these effects, and the non-crossing of the droplet trajectories is not a necessary condition for the use of the flow tube method. The effects of trajectory crossings are analyzed and discussed in detail in different situations for the three methods. This work is helpful for understanding and accurately calculating the droplet impingement characteristics and is of great significance for simulations of the aircraft icing process and anti/de-icing range.

The occurrence of aircraft icing can significantly affect flight performance. One of the most important aspects in the study of anti-icing technology for aircraft is the distribution of overflow water. Owing to the external airflow pressure, shear stress, and surface tension, the water film breaks up to form steady rivulets. Experiments on NACA0012 airfoil surfaces were conducted based on an open straight-flow and low-speed wind tunnel. Simultaneously, an engineered three-dimensional rivulet model considering the surface roughness was established based on the energy-minimum principle. A comparison between the simulation and experimental results shows that the errors in the water film breakup location and the flow velocity of rivulets are less than 20%, and the errors in the spacing and width of rivulets are less than 40%. In addition, the effects of surface temperature and uniform roughness on water film breakup were investigated. Furthermore, the rivulet model was applied to the numerical calculation of the thermal performance of hot-air anti-icing systems. The simulations reveal that the uniform roughness of the wing surface causes the water film to break earlier. As the surface roughness increases, the thickness, spacing, and width of the rivulets increase, and the rivulet flow velocity decreases.

The electronic oxygen regulator (EOR) is a new type of aviation oxygen equipment which uses electronic servo control technology to control breathing gas pressure. In this paper, the control method of EOR was studied, and the dynamic model of the aviation oxygen system was established. A disturbance-observer-based controller (DOBC) was designed by the backstepping method to achieve the goal of stable and fast breath pressure control. The sensitivity function was proposed to describe the effect of inspiratory flow on breath pressure. Combined with the frequency domain analysis of the input sensitivity function, the parameters of the DOBC were analyzed and designed. Simulation and experiment studies were carried out to examine the control performance of DOBC in respiratory resistance and positive pressurization process under the influence of noise and time delay in the discrete electronic control system, which could meet the aviation physiology requirements. The research results not only verified the rationality of the application of DOBC in the breath control of EOR, but also proved the effectiveness of the control parameters design method according to the frequency domain analysis, which provided an important design basis for the subsequent study of EOR.

Developing renewable energy, especially solar energy related, is of great importance for securing our future energy society. Steam generation in nanofluids based on solar radiation has been increasingly studied. It has been determined that the efficiency of steam generation is significantly enhanced when nanoparticles are seeded into the fluid owing to their unique radiative heat transfer performance. The nanoparticles trap solar energy inside the fluid and convert it into thermal form, which dramatically accelerates the steam generation process. In this study, we experimentally investigated different nanofluids that directly absorb solar energy to generate steam. Ag nanofluid, Au nanofluid and MWCNT nanofluid with different concentration have been carefully investigated. We analyzed the temperature increase and steam generation combined with the calculation of the efficiency factor from radiative heat transfer. The heating power and steam generation power of different nanofluids and the same nanofluid with different concentrations were compared. For Au nanofluid with concentration of 0.5 wt‰, the absorbed solar energy for heating the volume and generating steam is 6 and 40 times higher than those of pure water, respectively. We concluded that localized boiling generates steam rapidly in nanofluids based on the observation of three types of nanofluids. Furthermore, the heating power and steam generation power of different nanofluids increase with concentration. Moreover, the difference between the efficiency factors results in varied volume heating and steam generation efficiencies for different nanofluids despite identical concentrations.

To simulate aircraft thermal anti-icing systems and solve the conjugate heat transfer of air-droplet flow and solid skin, the heat and mass transfer model of the runback water on the anti-icing surface was combined with the heat conduction equation of the skin by loosely coupled methods. According to the boundary conditions used for the runback water conservation equations, two loose-coupling methods for the heat exchange between the runback water and the solid skin were developed based on surface heat flow and surface temperature, respectively. The anti-icing and ice accretion results of a NACA 0012 electro-thermal anti-icing system were obtained by the two loose-coupling methods. The heat flow-based method directly solves the thermodynamic model of the runback water without any extra assumptions, but the convergence rate is relatively slow. On the other hand, the temperature-based method achieves higher calculation speed, but the freezing point is extended to an artificial temperature range between water and ice phases. When the value of the artificial temperature range is small, the results obtained by the temperature-based method are consistent with those of the heat flow-based method, indicating that the effect of freezing point extension can be ignored for thermal anti-icing simulation. Furthermore, the solutions of the two methods are in acceptable and comparable agreement with the experimental and simulative results in the literature, confirming their feasibility and effectiveness. In addition, it is found that the ice thicknesses and ice shapes rise obviously near the runback water limits as a result of the transverse heat conduction of the solid skin.