An implicit time marching Galerkin method for the simulation of icing phenomena with a triple layer model

Chauvin, Rémi; Bennani, Lokman; Trontin, Pierre; Villedieu, Philippe

doi:10.1016/j.finel.2018.07.003

Cited by 22 publications

(14 citation statements)

References 28 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…During this power cycle ice may build up or melt and liquid water may run downstream under the effect of the aerodynamic forces. This may lead to several possible states as illustrated in Figure 7b [7].…”

Section: Application To Electrothermal De-icingmentioning

confidence: 99%

“…In this work, a finite volume solver called ETIPS2D is used to simulate the electro-thermal system (see [2] for a description). Concerning the unsteady ice accretion and melting problem, an unsteady mixed finite volume-Galerkin solver called MiLeS2D is used (see [7] for a detailled description). A global illustration is shown in Figure 7b.…”

Section: Application To Electrothermal De-icingmentioning

confidence: 99%

“…(b) Global illustration of the case. Several states are possible for the icing process leading to different mass and energy transfer terms (see [7] for details) Before moving to the test case, note that MiLeS2D has a complex algorithmic structure, the explicit tracking of phase change fronts requiring to take into account several possible states and possible shifts between these states. The six possible states, illustrated in Figure 7b, are defined as:…”

Section: Application To Electrothermal De-icingmentioning

confidence: 99%

See 2 more Smart Citations

A non-overlapping optimized Schwarz method for the heat equation with non linear boundary conditions and with applications to de-icing

Bennani

Trontin

Chauvin

et al. 2020

Computers & Mathematics with Applications

View full text Add to dashboard Cite

When simulating complex physical phenomena such as aircraft icing or de-icing, several dedicated solvers often need to be strongly coupled. In this work, a non-overlapping Schwarz method is constructed with the unsteady simulation of de-icing as the targetted application. To do so, optimized coupling coefficients are first derived for the one dimensional unsteady heat equation with linear boundary conditions and for the steady heat equation with non-linear boundary conditions. The choice of these coefficients is shown to guarantee the convergence of the method. Using a linearization of the boundary conditions, the method is then extended to the case of a general unsteady heat conduction problem. The method is tested on simple cases and the convergence properties are assessed theoretically and numerically. Finally the method is applied to the simulation of an aircraft electrothermal de-icing problem in two dimensions.

show abstract

Section: Application To Electrothermal De-icingmentioning

confidence: 99%

Section: Application To Electrothermal De-icingmentioning

confidence: 99%

Section: Application To Electrothermal De-icingmentioning

confidence: 99%

See 1 more Smart Citation

A non-overlapping optimized Schwarz method for the heat equation with non linear boundary conditions and with applications to de-icing

Bennani

Trontin

Chauvin

et al. 2020

Computers & Mathematics with Applications

View full text Add to dashboard Cite

show abstract

“…The formation of wall film during the impingement process has been modelled and investigated with different methods, although the classic Messinger model for icing simulation does not consider the water film flow on the ice layer. Some researches (Myers and Charpin, 2004) (Nakakita et al, 2010) (Cao et al, 2016) (Chauvin et al, 2018) adopted the lubrication equation based on the assumption of continuous and thin liquid film to model the film height, which is validated for the phase where the film has formed well and is unable to consider the initial icing roughness, which has great influence on calculating the shape of ice accurately. The global "equivalent sand grain roughness" whose height does not change with time and space was used to simulate effects of the surface roughness on the flow and heat transfer process.…”

Section: Introductionmentioning

confidence: 99%

Modeling and Analysis Thin Wall Film Formation with Integrated Discrete and Continuous Phase Simulation Approach

Peng¹,

Liu²

2019

Proceedings of Global Power &Amp; Propulsion Society

View full text Add to dashboard Cite

Water film will be initiated and formed under the continuous impingement of incoming supercooled droplets, and directly affect the ice formation, which is always vital for the performance and safety of aero-engine. A hybrid algorithm integrating Discrete Phase Model (DPM) and Volume of Fluid (VOF) was established to simulate all the stages from droplet initiation to the film well-formed which can make the numerical simulation of icing more realistic and accurate. The transformation criterion from the particle to liquid volume fraction and the distribution of source terms has been improved compared with similar previous approaches. Three verification cases, namely single droplet train impact, staggered droplets impact and two glycerine-solution droplet impact, were conducted to verify the mass conservation and the applicability to both the sparse and dense two-phase flow. The results showed that the relative error of mass transformation is no more than 0.5%. The comparison between simulation and experiment result demonstrated that the hybrid algorithm has the ability to simulate the flow behaviour of droplets accurately. The thin water film formation on a typical wing of NACA0012 airfoil was then analysed, the results of which illustrated that the integrated algorithm could successfully capture the characteristics of liquid film formation and motion.

show abstract

“…The Messinger model has been implemented in the first generation of ice crystal icing codes, for both non-porous ( [3], [7]) and porous ( [8]) accretions. It has also been used in cases where unsteady heating from an anti-icing or deicing system is modelled; in this case a three-layer (water-ice-water) model may occur [9].…”

Section: Introductionmentioning

confidence: 99%

A Three-Layer Thermodynamic Model for Ice Crystal Accretion on Warm Surfaces: EMM-C

Bucknell¹,

McGilvray²,

Gillespie³

et al. 2019

SAE Technical Paper Series

View full text Add to dashboard Cite

Ingestion of high altitude atmospheric ice particles can be hazardous to gas turbine engines in flight. Ice accretion may occur in the core compression system, leading to blockage of the core gas path, blade damage and/or flameout. Numerous engine powerloss events since 1990 have been attributed to this mechanism. An expansion in engine certification requirements to incorporate ice crystal conditions has spurred efforts to develop analytical models for phenomenon, as a method of demonstrating safe operation. A necessary component of a complete analytical icing model is a thermodynamic accretion model. Continuity and energy balances are performed using the local flow conditions and the mass fluxes of ice and water that are incident on a surface to predict the accretion growth rate. In this paper, a new thermodynamic model for ice crystal accretion is developed through adaptation of the Extended Messinger Model (EMM) from supercooled water conditions to mixed phase conditions (ice crystal and supercooled water). A novel three-layer accretion structure is proposed and the underlying equations described. The EMM improves upon the original model for airframe icing, the Messinger Model, by permitting a linear temperature gradient through the ice and water layers. This in turn allows prediction of the time over which water exists in isolation on an initially warm surface, before an ice layer forms. This is of particular interest to engine icing, as surfaces may initially be significantly above freezing temperature, before cooling on exposure to ice particles. The method is solved in a multi-step approach, where the overall exposure time is divided into discrete windows, and the calculation performed over each window. This allows the local flow conditions to be updated between windows, permitting the incorporation of a reducing flow enthalpy due to particle evaporation, as well as transient engine operation. Model results are then compared to experimental results. Comparisons are made to solutions generated using the standard Messinger Model.

show abstract

An implicit time marching Galerkin method for the simulation of icing phenomena with a triple layer model

Cited by 22 publications

References 28 publications

A non-overlapping optimized Schwarz method for the heat equation with non linear boundary conditions and with applications to de-icing

A non-overlapping optimized Schwarz method for the heat equation with non linear boundary conditions and with applications to de-icing

Modeling and Analysis Thin Wall Film Formation with Integrated Discrete and Continuous Phase Simulation Approach

A Three-Layer Thermodynamic Model for Ice Crystal Accretion on Warm Surfaces: EMM-C

Contact Info

Product

Resources

About