Introdução à Simulação de Gelo em Multifísica COMSOL

A formação e fusão de geada são fenômenos multifísicos complexos que desempenham um papel crítico em muitos sistemas industriais e naturais. Na engenharia aeroespacial, a geada em superfícies de aeronaves pode interromper o fluxo de ar e reduzir o elevador. Na refrigeração, o acúmulo de geada em bobinas evaporadoras degrada a eficiência de transferência de calor e aumenta o consumo de energia. Na ciência climática, o acúmulo de geada em vegetação e solo influencia os balanços de energia e umidade locais. A simulação precisa desses processos requer um tratamento acoplado de transferência de calor, transporte de massa e mudança de fase. COMSOL Multifísica, com seu módulo flexível de CFD e acoplamento multifísico integrado, proporciona um ambiente onde engenheiros e pesquisadores podem modelar a dinâmica de geada com alta fidelidade. Este artigo expande a configuração fundamental, métodos numéricos, abordagens de validação e aplicações práticas de formação de geada e simulações de fusão em COMSOL.

Princípios físicos da formação e fusão de geada

O gelo forma-se quando o vapor de água no ar se eleva diretamente sobre uma superfície sólida cuja temperatura está abaixo do ponto de geada. O ponto de geada depende da pressão e temperatura de vapor local; é tipicamente inferior ao ponto de orvalho porque a pressão de vapor de saturação sobre o gelo é menor do que sobre a água líquida. O derretimento ocorre quando a temperatura da superfície sobe acima de 0°C, fazendo com que o gelo se transfira para a água líquida, que pode então evaporar ou correr, dependendo das condições de superfície e de fluxo de ar. As equações que regem incluem a conservação da massa, do momento e da energia para a fase de fluido, uma equação de transporte para vapor de água, e um equilíbrio térmico que explica o calor latente da sublimação e fusão.

Equações Governantes para a Fase Fluída

No domínio do fluxo de ar, as equações de Navier- Stokes (RANS) médias de Reynolds são normalmente resolvidas para modelar o fluxo turbulento em torno das superfícies propensas a geadas. A equação da continuidade garante a conservação de massa, enquanto as equações de momento incorporam o fechamento de turbulência, usando frequentemente os modelos K- ε ou K- ω SST. A equação de energia inclui a transferência de calor convectiva e, quando associada ao transporte de humidade, o efeito da libertação ou absorção de calor latente. Para o campo de humidade, é resolvida uma equação de conversão- difusão para a fracção de massa de vapor de água, com um termo de fonte ou de dissipação representando sublimação ou deposição. A condição de fronteira na superfície de geada acopla o fluxo de vapor ao equilíbrio de calor através da relação Clausio- Clapeyron, que define a pressão de vapor de equilíbrio sobre o gelo. Estes acoplamentos não lineares exigem uma abordagem numérica robusta; a solução multifrontal da COMSOL ou a solução GMRES com pré-condicionamento é normalmente utilizada para lidar a rigidez.

Modelo de Mudança de Fase para Gelo

A mudança de fase de vapor para gelo (deposição) e de gelo para líquido (derretimento) é modelada utilizando um método de porosidade entalpia ou capacidade de calor equivalente. Em COMSOL, a mudança de fase pode ser implementada através de uma função de passo suavizada que define a fração de massa de gelo em função da temperatura e concentração de vapor. O calor latente é incluído como um termo fonte volumétrica na equação de energia, proporcional ao tempo derivado da fração de massa de gelo. Para o crescimento da geada, um submodelo rastreia a porosidade e condutividade térmica da camada de geada, ambos evoluindo com densidade e espessura. As correlações empíricas para a condutividade térmica de geada (por exemplo, do trabalho de Hayashi ou O'Neal e Árvore) são frequentemente introduzidas como funções de densidade, que dependem da taxa de deposição e do tempo. Estes submodelos são críticos porque uma camada de geada densa e baixa porosidade conduz calor diferentemente a uma camada de alta porosidade.

Configurar o modelo COMSOL passo a passo

Criação de geometria e considerações sobre malha

Comece definindo a geometria no ambiente CAD integrado da COMSOL ou importando de uma ferramenta externa. Um modelo típico inclui uma placa plana que representa uma extremidade ou superfície da asa do trocador de calor, colocada dentro de um canal de fluxo que define o domínio do ar. Para geometrias axissimétricas (por exemplo, um tubo de arrefecimento), uma configuração axissimétrica 2D pode reduzir o custo computacional. A malha deve ser refinada perto da superfície de cobertura para resolver os gradientes agudos de temperatura, concentração de vapor e velocidade. Uma malha de camada de contorno com 5-10 camadas de prisma e uma primeira espessura de célula de modo que y+ □1 é recomendada para modelos de turbulência que resolvam a subcamada viscosa. Os estudos de independência da grelha devem ser realizados comparando os resultados de espessura de geada e de transferência de calor em malhas mais grossas, médias e finas.

Definindo Interfaces Físicas

COMSOL fornece interfaces físicas dedicadas que podem ser combinadas:

  • Transferência de calor em fluidos (ht) – para distribuição de temperatura no ar e no sólido.
  • Transferência de calor em sólidos (ht) – para a parede ou substrato.
  • Transportes de Espécies Diluídas (tds) – para fração de massa de vapor de água no ar.
  • Fluxo turbulento, k-ε ou k-ω (spf) – para transporte de momentum.
  • Acoplamentos multifísicos – acoplamento de fluxo não isotérmico (pseudo-acoplamento) e um acoplamento personalizado “Deposição de Frutos” que liga o fluxo de vapor ao equilíbrio térmico através da condição limite.

Para incorporar o crescimento da geada como um limite móvel, pode-se usar o recurso Movendo a Mesh] ou Deformado . Alternativamente, para simplicidade, pode-se usar um método de grade fixa com porosidade efetiva, onde a camada de geada é tratada como um meio poroso com espessura e propriedades térmicas dependentes do tempo.

Propriedades do Material

Os dados precisos do material são essenciais. Para o ar, a densidade segue a lei do gás ideal, a condutividade térmica e a viscosidade são dependentes da temperatura (usar expressões incorporadas). A difusividade do vapor de água no ar pode ser obtida a partir de fórmulas empíricas (por exemplo, de Cussler ou do Manual de Perry). Para a camada de geada, a densidade ρ frost normalmente varia de 50 a 300 kg/m3 dependendo das condições de deposição. Uma expressão como ρ frost = 650 − 2.0×106/(T frost+273.15) ou correlação de Hayashi (1977) pode ser usada. A condutividade térmica do geada k frost (W/(m·K)) é frequentemente dada por k frost = 0,024 + 7,52×10−4 ρ frost + 4,1×10−9 ρ frost3. O calor latente da sublimação é aproximadamente 2,83×106 J/kg, e para derretiá- lo é 3,34×105 J/kg.

Condições de Limite e Inicial

Para a simulação de cobertura normal:

  • Início: velocidade constante (1–5 m/s típico para refrigeração), temperatura (por exemplo, -5 a 10°C), e humidade relativa (50–90%).
  • Saída: estado de pressão ou de saída de gás zero.
  • Superfície de refrigeração: parede com temperatura prescrita (por exemplo, -15°C para uma placa fria) ou acoplada a um balanço energético que inclui convecção e calor latente. A concentração de vapor na parede é definida para o valor de saturação sobre o gelo a essa temperatura de parede (através de uma função ou tabela de pesquisa).
  • Condições iniciais: todo o domínio de fluxo definido para entrada de temperatura e umidade, fração de massa zero de gelo. A simulação é dependente do tempo; os tempos de execução típicos vão de minutos a horas em tempo físico, com passos de tempo na ordem de 1-10 segundos.

Simulando a formação de geada: fluxo de trabalho numérico

Uma vez que o modelo é construído, execute um estudo dependente do tempo. A sequência do solucionador da COMSOL segrega automaticamente a física — primeiro fluxo de resolução (estacionário ou transitório) em seguida, as equações escalares — mas, para casos fortemente acoplados, uma abordagem totalmente acoplada pode ser mais estável. Monitore as seguintes saídas:

  • Espessura de geada em função do tempo (pós-processar a interface em movimento ou o campo de fração de massa de gelo).
  • Coeficiente de transferência de calor superficial (computado do fluxo de calor total através da parede dividido pela diferença de temperatura).
  • Massa total de geada por unidade de área.
  • Perfis de temperatura e velocidade perto da parede.

As principais percepções de um modelo de formação de geada incluem o período de crescimento rápido inicial (dominado por transferência de massa), seguido de uma fase de crescimento mais lenta, à medida que a camada de geada isola a superfície. O modelo também pode prever a distribuição espacial da geada – muitas vezes mais espessa perto da borda superior de uma placa plana devido a taxas de difusão de vapor mais elevadas. As visualizações da camada de geada usando parcelas de corte ou gráficos de linha ao longo da superfície são inestimáveis para comparação de design.

Modelação de Frost Melting e degelo Ciclos

As simulações de fusão requerem uma mudança de fase reversa: a camada de gelo passa para água líquida, que então flui ou evapora. Em COMSOL, uma abordagem é incluir um modelo de fusão de camada fina [] que ativa quando a temperatura da superfície excede 0°C. A frente de fusão pode ser rastreada com um método de degelo ou de interface, mas para muitas aplicações de engenharia, um método integral de equilíbrio térmico simplificado é suficiente. Imponha um fluxo de calor ou uma condição de limite de temperatura elevada na parede (por exemplo, 30°C para um aquecedor de descongelamento elétrico) e resolva pelo tempo necessário para derreter completamente uma determinada espessura de geada. A simulação deve ser considerada para a energia consumida pelo calor latente de fusão, bem como aquecimento sensível da geada porosa e do substrato. A física adicional pode incluir o fluxo de filme de água impulsionado pela gravidade e cisalhamento do fluxo de ar, que pode ser modelado com o Thin-Film Flow[FT]: interface vapor de transferência de calor e vapor.

A otimização do ciclo de descongelamento pode ser realizada parametrizando a potência de aquecimento, a duração do aquecimento e o período de drenagem, e usando o módulo de otimização da COMSOL para minimizar o uso de energia, garantindo a remoção completa da geada. Por exemplo, pode-se estudar a penalidade energética de deixar geada residual versus a energia de superaquecimento.

Validação e Boas Práticas

A validação de modelos de geada contra dados experimentais é crucial. Os benchmarks publicados incluem o trabalho de Liu e Jacobi (2013) sobre o crescimento de geada em placas planas em convecção forçada e os estudos de Padhmanaban e Sherif (2000) sobre propriedades de geada. Tipicamente, comparações são feitas para:

  • Espessura de geada vs. tempo
  • Temperatura da superfície de gelo
  • Redução global do coeficiente de transferência de calor devido à geada

Os usuários devem calibrar as correlações empíricas para a condutividade térmica e densidade de geada para suas faixas específicas de temperatura e umidade. Uma análise de sensibilidade sobre esses parâmetros ajuda a quantificar a incerteza. Estudos de refinamento de malha devem direcionar uma incerteza de espessura de geada de <5%. Additionally, for turbulent flows, ensure that the y+ value at the first cell is checked—the built‑in wall functions in the k-ε model may require the first cell to lie in the log‑law region (y+ > 30), enquanto um modelo de turbulência de baixa-Re (como k-ω) requer y+ ~1. Mismatch pode levar a previsões erradas de transferência de calor e massa.

Considerações Avançadas: Turbulência, Acoplamento e Demanda Computacional

A interação entre geada e turbulência é dupla: a rugosidade do geada altera as estruturas turbulentas de parede próxima e os turbulentos turbulentos aumentam o transporte de vapor. Alguns modelos avançados incorporam uma função de rugosidade no tratamento da parede de turbulência, aumentando a rugosidade da superfície à medida que a geada se acumula. O COMSOL permite funções definidas pelo utilizador para a altura de rugosidade em função da densidade ou do tempo de geada, que podem ser adicionadas às condições de contorno da parede da interface de fluxo.

Solução totalmente acoplada de fluxo, calor e transporte de umidade com frente de geada móvel pode ser computacionalmente pesada. Durante longos tempos físicos (por exemplo, horas de crescimento de geada), recomenda-se o passo do tempo adaptativo com uma tolerância de 0,001. Os usuários podem aproveitar as varreduras paramétricas eficientes da COMSOL para avaliar múltiplas condições de entrada. Para grandes modelos 3D, uma malha grossa nas regiões não-frosting e uma malha fina apenas perto da superfície pode reduzir a contagem de células. O tempo total de simulação para um modelo 2D (vários mil elementos) com 30 minutos de tempo físico pode levar 1-4 horas em uma estação de trabalho, enquanto um modelo 3D de trocador de calor de ponta e tubo pode levar 24 horas ou mais.

Aplicações e Relevância Industrial

Sistemas anti-gelo Aeroespacial

O gelo difere do esmalte ou do gelo de rima, mas a sua formação em asas durante as operações no solo pode afectar o desempenho da descolagem. Os modelos COMSOL ajudam a projetar sistemas anti-gelo electrotérmicos prevendo a potência necessária para manter superfícies acima do ponto de geada. Estudos paramétricos sobre o layout do aquecedor e o ciclo de trabalho podem ser realizados no mesmo modelo, incluindo o efeito da fusão de geada e do escoamento de água.

Desempenho da bomba de calor e refrigeração

Em bobinas evaporadoras de bombas de calor de fonte de ar, a acumulação de geada reduz o fluxo de ar e a transferência de calor, levando à degradação do desempenho. Os engenheiros usam COMSOL para determinar os critérios de iniciação de descongelamento ideais com base na espessura ou na queda de pressão de geada medida. O modelo também pode explorar novos revestimentos de superfície (hidrofílicos ou hidrofóbicos) que alteram a nucleação e o crescimento de geada. Ao modificar a condição de contorno do ângulo de contato no transporte de umidade, os usuários podem avaliar como a energia superficial afeta a deposição.

Sistemas de construção eficientes em termos energéticos

A geada em ventiladores de recuperação de calor (VFCs) em climas frios impede a operação. Simulações ajudam na concepção de estratégias de bypass ou pré-aquecimento que minimizem o acúmulo de geada, mantendo o conforto. Modelos de construção-energia associados podem ser simplificados usando a co-simulação COMSOL com ferramentas como EnergyPlus (via interface externa).

Proteção contra o gelo da turbina do vento

O gelo de geada e de rima nas lâminas de turbinas reduzem a potência e causam desequilíbrio. Modelos multifísicos que incluem rotação de lâminas (usando interfaces rotativas de máquinas) e cargas de vento estocásticas podem simular a acumulação de geadas em condições atmosféricas variadas. Estes modelos apoiam o desenvolvimento de sistemas de desfrimento ativos, como aquecimento resistivo ou vibração ultrasssônica, validados no ambiente COMSOL.

Desafios e Instruções de Pesquisa Emergentes

Apesar do poder da COMSOL, a modelagem de geada continua a ser um desafio devido à complexidade da nucleação microescala e à natureza estocástica do crescimento do cristal de geada. Muitos modelos atuais dependem de correlações empíricas que são validadas apenas em intervalos de parâmetros estreitos. Pesquisas futuras estão se movendo para a simulação numérica direta (DNS) de interações de gotas de congelamento e formação de grades de cristal, embora estes sejam computacionalmente caros para problemas de macroescala. Outra tendência emergente é o uso de aprendizado de máquina para substituir partes do modelo de geada – por exemplo, substituindo a correlação de condutividade térmica porosa com uma rede neural treinada em dados experimentais.

Conclusão

Simulando a formação e fusão de geadas em modelos CFD COMSOL fornece aos engenheiros uma estrutura robusta para estudar interações multifísicas que são difíceis de reproduzir em experimentos. Ao configurar cuidadosamente a geometria, interfaces físicas e propriedades materiais, os usuários podem obter informações temporais e espaciais detalhadas sobre o crescimento de geadas, degradação de transferência de calor e comportamento de fusão. A validação contra dados experimentais continua sendo um passo crítico para garantir credibilidade de simulação. Com os avanços contínuos nos recursos computacionais e submodelos físicos, o COMSOL continuará a ser uma ferramenta chave para projetar sistemas resistentes à geada e otimizar estratégias de descongelamento em setores aeroespacial, refrigerador e energético. Para leitura adicional, consulte as Bibliotecas de Aplicações COMSOL para exemplos de mudanças de fase e de transporte de umidade, e reveja os artigos de COMSOL sobre modelagem de geadas em trocadores de calor e e [FLT].