É necessária uma assinatura do JoVE para visualizar este conteúdo. Faça login ou inicie seu teste gratuito.

Artigo de método

Simulação CFD de geada em superfícies frias horizontais

574 vistas

DOI:

10.3791/68133

12 de setembro de 2025

Neste artigo

Resumo

Aqui, apresentamos um modelo numérico para simular a formação de geada em superfícies frias horizontais usando um modelo multifásico euleriano com a abordagem de mudança de fase de Lee. O modelo atualiza dinamicamente a fração de volume de gelo para capturar a variação de densidade e é validado em relação a dados experimentais sobre espessura, densidade e distribuição de gelo.

Resumo

A formação de geada é um fenômeno comum observado em muitos campos diferentes, incluindo refrigeração, construção e processamento de gás natural. No entanto, devido à sua natureza complexa, o desenvolvimento de um modelo numérico preciso e confiável continua sendo um desafio significativo. Apesar dos esforços anteriores para resolver esse problema, os modelos atuais ainda têm certas limitações. Este artigo apresenta um modelo numérico modificado para a formação de geada, desenvolvido com base nos mecanismos fundamentais subjacentes à formação de geada. O modelo utiliza uma abordagem de fluxo multifásico euleriana acoplada ao modelo de mudança de fase de Lee. Além disso, a abordagem para determinar a fração máxima do volume de gelo é atualizada, permitindo que o modelo considere as variações de densidade durante o processo de congelamento. O modelo é rigorosamente validado comparando-o com dados experimentais sobre espessura, densidade e distribuição de vários estudos. Os resultados indicam que o desvio relativo absoluto médio (MARD) para a espessura da geada é de 8,97%, enquanto o MARD para densidade é de 16,06%. Além disso, a morfologia da geada prevista pelo modelo corresponde às observações experimentais relatadas na referência.

Introdução

A formação de geada é um fenômeno comum observado em muitos campos diferentes. O acúmulo de gelo nas superfícies do trocador de calor prejudica significativamente a eficiência da transferência de calor1, obstrui o fluxo de fluido e interrompe o desempenho geral dos trocadores de calor2, dificultando sua operação normal3. Portanto, entender os mecanismos e o comportamento da formação de geada é fundamental para resolver esse problema em sistemas de refrigeração4. Nas últimas décadas, um corpo substancial de pesquisas tem se dedicado a investigar as causas e características da formação de geadas nesses sistemas.

Estudos experimentais demonstraram que a formação de geada é influenciada por vários fatores, incluindo temperatura do ar, umidade e temperatura da superfície fria 5,6,7,8,9,10. Numerosos achados experimentais indicaram que temperaturas mais baixas do ar de entrada tendem a resultar em camadas de gelo mais espessas6, enquanto níveis mais altos de umidade contribuem para a formação de camadas de gelo mais densas7. Song et al. estudaram a formação de geada em superfícies horizontais e descobriram que variações cíclicas de temperatura da superfície fria podem causar derretimento na interface da camada de geada, o que influencia significativamente a taxa de formação de geada, a espessura da camada de geada e a densidade dinâmica da geada8. Outra pesquisa examinou a morfologia e a distribuição da geada. Jeong et al. observaram em seus experimentos que a geada se forma inicialmente perto da enseada, levando à ocorrência de um fenômeno conhecido como colina de geada9. Noorshams et al. investigaram a formação de gelo em uma superfície de tubo circular horizontal e descobriram que as camadas de gelo nas superfícies frontal e traseira do cilindro eram mais espessas do que as das superfícies superiores10. Além disso, vários estudos 11,12,13,14,15,16 desenvolveram modelos para prever a espessura unidimensional da camada de gelo, utilizando padrões experimentais de formação de gelo junto com abordagens teóricas e empíricas. Jones e Parker desenvolveram um modelo preditivo para a espessura da geada com base na teoria da difusão molecular11. A discrepância entre o modelo e os dados experimentais permaneceu abaixo de 30% durante um período de 3 horas. Com a progressão da tecnologia computacional, um número crescente de pesquisadores recorreu à Dinâmica de Fluidos Computacional (CFD) para simular a formação de geadas. Ao contrário dos modelos unidimensionais tradicionais, as simulações CFD oferecem benefícios significativos, particularmente na visualização da distribuição da espessura do gelo e dos perfis de temperatura. Cui et al. realizaram simulações CFD de formação de geada com base na teoria da nucleação12. Suas previsões de espessura de geada mostraram um desvio de menos de 13% dos dados experimentais fornecidos por Lenic et al.13. Paralelamente, estudos de CFD sobre condensação em tubos estruturados mostraram que características geométricas como covinhas14 ou passos helicoidais15 aumentam a transferência local de calor e massa. Recentemente, You et al.16 desenvolveram um modelo CFD dinâmico baseado em malha que caracteriza a camada de gelo como um meio poroso em crescimento e incorpora a difusão de vapor diretamente, alcançando um desvio relativo de menos de 5%, mantendo baixo custo computacional. Essas descobertas ressaltam o potencial do CFD na resolução de fenômenos complexos de mudança de fase, oferecendo informações valiosas para a modelagem da formação de geada.

Em conclusão, um número significativo de estudos 5,6,7,8,9,10,11,12,13 explorou a formação de geada em superfícies frias, contribuindo para uma compreensão evolutiva dos padrões de formação de geada sob vários parâmetros. Embora vários modelos numéricos tenham sido desenvolvidos, incorporando várias dimensões e mecanismos, eles frequentemente carecem de validação abrangente. A maioria dos estudos 6,7,8,9,10,11,12,13 valida principalmente modelos usando espessura de geada, limitando sua aplicabilidade mais ampla 16. Para superar essas limitações, este artigo apresenta um modelo numérico que integra o modelo de mudança de fase de Lee e o modelo de fluxo multifásico euleriano, com ênfase nos mecanismos centrais subjacentes à formação de geada. Além disso, é introduzido um novo método para calcular o limite superior da fração de volume de geada, tendo em conta as alterações da densidade de geada dependentes do tempo, colmatando assim as deficiências dos modelos anteriores. A precisão e a confiabilidade do modelo proposto são avaliadas de vários ângulos, incluindo espessura de geada, flutuações de densidade e padrões de formação de geada em diferentes condições experimentais. Essa extensa validação oferece uma estrutura teórica robusta para prever com mais precisão o comportamento do gelo em aplicações do mundo real.

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Protocolo

1. Modelo físico e malha

  1. Abra o SpaceClaim, selecione a guia Esboço e escolha a opção Retângulo na função Criar.
  2. Crie um modelo geométrico 2D no plano XOY com um comprimento de 500 mm ao longo do eixo x e uma largura de 15 mm ao longo do eixo y.
  3. Abra o ICEM, vá para a guia Arquivo, selecione a guia Geometria , escolha a opção Abrir geometria e importe o modelo de geometria 2D.
  4. Selecione o limite do modelo 2D, abra a função Criar peça na guia Peças e atribua nomes aos diferentes limites.
  5. Selecione a guia Bloqueio e clique em Criar bloco. Na janela Criar bloco, marque a opção Herdar nome da peça . Em seguida, vá para a guia Criar bloco e selecione Inicializar bloco.
  6. Na guia Inicializar blocos, selecione Bloqueio de superfície 2D como o tipo e marque Inicializar com configurações.
  7. Na guia Escolher, selecione Superfície de formulário para a opção Método. Em seguida, em Superfícies, clique no ícone Selecionar superfície(s) e escolha o Plano 2D na janela de desenho.
  8. Na guia Bloqueio de superfície, selecione Principalmente mapeado para a opção Método, All Quad para o tipo de malha de face livre e ICEM CFD Quad para o método de malha de face livre.
  9. Na guia Mesclar blocos entre curvas, selecione Tudo para o Método e defina Ignorar tamanho como 0,0.
  10. Na guia Parâmetros de pré-malha, em Parâmetros de geração de malha, clique em Parâmetros de borda. Em Edge, clique no ícone Select Edge(s) e escolha o Edge na direção X na interface gráfica.
  11. Em Parâmetros de geração de malha, defina Nós como 500, Espaçamento 1 a 1e+10, Proporção 1 a 2, Espaçamento 2 a 1e+10, Proporção 2 a 2 e Espaço máximo como 1e+10. Em seguida, marque a opção Copiar parâmetros .
  12. Na opção Parâmetros de malha, defina a Lei da malha como Bigeométrica e clique em Aplicar. Na guia Parâmetros de pré-malha, na opção Parâmetros de malha, clique em Parâmetros de borda.
  13. Para a opção Aresta em Parâmetros de geração de malha, clique no ícone Selecionar aresta (s) e escolha a aresta na direção y na interface gráfica.
  14. Em Parâmetros de malha, defina Nós como 150, Espaçamento 1 a 1e+10, Proporção 1 a 2, Espaçamento 2 a 1e+10, Proporção 2 a 2 e Espaço Máximo como 1e+10. Em seguida, marque a opção Copiar parâmetros .
  15. Na opção Parâmetros de malha, defina a Lei da malha como Bigeométrica e clique em Aplicar. Na opção Bloqueio, ative a função Pré-malha e clique em Sim quando solicitado.
  16. Clique com o botão direito do mouse na pré-malha e selecione Converter em malha não estruturada no menu de contexto.
  17. Vá para a guia Malha de saída e clique em Configuração do solucionador. Na Configuração do Solver, selecione ANSYS Fluent como o solver e clique em Aplicar.
  18. Na guia Saída, clique em Gravar entrada. Na Nova Janela 1, clique em Salvar; na nova janela 2, clique em Sim e depois em Salvar; na nova janela 3, clique em Abrir; e na nova janela 4, clique em Concluído.

2. Operação do software de simulação de formação de geada

  1. Abra o Ansys Fluent. Vá para a guia Arquivo e selecione a opção Malha em Ler. Vá para a malha de escala e defina mm em A malha foi criada em. O procedimento aqui apresentado baseia-se nas condições experimentais definidas no Caso 1 da Tabela 1.
  2. Nas configurações do Solver, escolha Tipo baseado em pressão , Formulação de velocidade absoluta e Tempo transiente.
  3. Defina a aceleração gravitacional na direção Y para -9,81. Em Modelos, clique em Energia e ative a Equação de Energia.
  4. Em Modelos, ative Viscoso e selecione k-épsilon (2 eqn). Na seção Modelo k-épsilon, selecione Padrão.
  5. Para Tratamento próximo à parede, escolha Funções de parede padrão. Em Modelo multifásico de turbulência, escolha Mistura.
  6. Na opção Constantes do modelo, defina Cmu, C1-Epsilon, C2-Epsilon, Número de Prandtl TKE, Número de Prandtl TDR, Número de Prandtl de Dispersão, Número de Prandtl de Energia, Número de Prandtl de Parede e Número de Schmidt Turbulento como 0,09, 1,44, 1,92, 1, 1,3, 0,75, 0,85, 0,85 e 0,7, respectivamente.
  7. Em Funções definidas pelo usuário, defina a mistura de viscosidade turbulenta, fase 1, fase 2 como nenhuma. Em Modelos, clique em Espécies e ative Transporte de Espécies.
  8. Na janela Modelo de transporte de espécies, marque Fonte de energia de difusão e selecione a fase 1.
  9. Em Modelos, ative Multifásico e escolha o modelo Euleriano . Na janela Modelo multifásico, selecione Parâmetros de fração de volume Formulação implícita e defina o Número de fases eulerianas como 2.
  10. Na janela Modelo multifásico, selecione a guia Fases , selecione a fase 1 - Fase primária, defina Nome como fase 1 e defina Material de fase como modelo de mistura.
  11. Na janela Modelo multifásico, selecione a guia Fases , selecione a fase 2 - Fase secundária, defina Nome como fase 2, defina Material da fase como gelo e ative Granular.
  12. Na janela de configuração de fases da guia Fases, selecione Propriedade de fase como Modelo de temperatura granular, defina Diâmetro como 0,0001.
  13. Na janela Propriedades granulares da guia Fases, defina Viscosidade granular como 1e-05, defina Viscosidade granular em massa como 30, defina Pressão de sólidos como lun-et-al, defina Temperatura granular como algébrica, defina Viscosidade de atrito como nenhuma, defina Limite de empacotamento como definido pelo usuário, defina Distribuição radial como lun-et-al, defina Módulo de elasticidade como derivado.
  14. Na guia Interação de fases da janela Modelo multifásico, selecione o identificador Forças . Em seguida , selecione a fase 1, fase 2, defina o coeficiente para wen-yu.
  15. Na guia Interação de fases da janela Modelo multifásico, selecione o identificador Forças . Em seguida , selecione a fase 2, fase 2, defina o coeficiente de restituição para 0,9.
  16. Na guia Interação de fases da janela Modelo multifásico, selecione a guia Área interfacial . Em seguida, escolha a opção ia-simétrica .
  17. Em Condições de limite, clique em Entrada e selecione a guia fase 1, defina a opção Fase como fase 1.
  18. Na janela Entrada de velocidade da guia fase 1, selecione a guia Momentum , defina o Método de especificação de velocidade como Magnitude, Normal para Limite, defina o Quadro de referência como Absoluto e defina a Magnitude da velocidade como 0,6.
  19. Na janela Entrada de velocidade da fase 1, selecione Térmico tag, defina Temperatura para 292.8. Na janela Entrada de velocidade da tag phase-1, selecione Tag Species , defina h2o como 0,008202.
  20. Em Condições de limite, clique em Entrada e selecione a tag de fase 2 , defina a opção Fase como fase 2.
  21. Na janela Entrada de Velocidade da tag fase 2, selecione a tag Momentum , defina o Método de Especificação de Velocidade como 'Magnitude, Normal ao Limite', defina o Quadro de Referência como Absoluto, defina a Magnitude da Velocidade como 0, defina a Temperatura Granular como 0,0001.
  22. Na janela Entrada de velocidade da tag de fase 2, selecione a guia Térmica , defina Temperatura como 273. Na janela Entrada de velocidade da tag fase 2, selecione a guia Multifase , defina Fração de volume como 0.
  23. Em Condições de limite, clique em Saída e selecione a tag de fase 1 , defina a opção Fase como fase 1.
  24. Na janela Saída de pressão da guia fase 1, selecione a guia Térmica , defina a Temperatura total de refluxo como 300.
  25. Na janela Saída de pressão da guia fase 1, selecione a guia Espécies , defina h2o como 0. Em Condições de limite, clique em Saída e selecione a guia fase 2, defina a opção Fase como fase 2.
  26. Na janela Saída de pressão da guia fase 2, selecione a guia Térmica e defina a Temperatura total de refluxo como 300.
  27. Na janela Saída de pressão da guia fase 2, selecione a guia Multifásico , defina Temperatura granular de refluxo como 0,0001, defina Método de especificação de fração de volume como Fração de volume de refluxo e defina Fração de volume de refluxo como 0.
  28. 2.28 Em Condições de limite, clique em Parede e selecione a guia parede fria , defina a opção Fase como mistura.
  29. Na janela de parede da guia parede fria, selecione a guia Momentum , defina a opção Movimento da parede como Parede estacionária, defina a opção Modelos de rugosidade da parede como Padrão, defina a Altura da rugosidade como 0 e defina a Constante de rugosidade como 0,5.
  30. Na janela de parede da guia de parede fria, selecione a guia Térmica , selecione Temperatura como Condições térmicas, defina Temperatura como 252,65 e defina Material como aço.
  31. Na Solução, clique em Métodos e abra a janela Métodos da Solução. Na janela Métodos de solução, selecione Fase acoplada SIMPLE como Esquema de acoplamento pressão-velocidade, selecione Mínimos quadrados baseados em células como Discretização espacial de gradiente, selecione Segunda ordem como Discretização espacial de pressão, selecione Upwind de primeira ordem como Discretização espacial de densidade, selecione Upwind de primeira ordem como Discretização espacial de momento, selecione Upwind de primeira ordem como Discretização espacial de fração de volume, selecione Upwind de primeira ordem como Discretização Espacial de Energia Cinética Turbulenta, selecione Contra o Vento de Primeira Ordem como Discretização Espacial da Taxa de Dissipação Turbulenta, selecione Contra o Vento de Primeira Ordem como Discretização Espacial de Energia, selecione Contra o Vento de Primeira Ordem como Discretização Espacial h2o de fase 1, selecione Implícito de Primeira Ordem como Formulação Transitória.
  32. Na Solução, clique em Controles e abra a janela Controles da Solução. Defina o Fator de Sub-Relaxamento de Pressão para 0.4, defina o Fator de Sub-Relaxamento de Densidade para 1, defina o Fator de Sub-Relaxamento das Forças Corporais para 1, defina o Fator de Sub-Relaxamento de Momentum para 0.4, defina o Fator de Sub-Relaxamento da Fração de Volume para 0.4, defina o Fator de Sub-Relaxamento de Temperatura Granular para 0.3, defina o Fator de Sub-Relaxamento de Energia Cinética Turbulenta para 0.3, defina o Fator de Sub-Relaxamento da Taxa de Dissipação Turbulenta para 0.3, defina o Fator de Sub-Relaxamento de Viscosidade Turbulenta para 0.3, defina o Fator de Sub-Relaxamento de Energia para 0.4, defina o Fator de Sub-Relaxamento de fase 1 h2o para 0.4.
  33. Na Solução, clique em Inicialização e abra a janela Inicialização da Solução. Selecione Inicialização padrão como o método de inicialização, selecione Relativo à zona da célula como quadro de referência.
  34. Na janela Inicialização da solução, defina a pressão manométrica como 0, defina a energia cinética turbulenta como 0,00135, defina a taxa de dissipação turbulenta como 0,001143987, defina a velocidade X da fase 1 como 0, defina a velocidade Y da fase 1 como 0, defina a fase 1 h2o como 0,008202, defina a velocidade X da fase 2 como 0, defina a velocidade Y da fase 2 como 0, defina a fração de volume da fase 2 como 0, defina a temperatura granular da fase 2 para 0.0001, defina a temperatura da fase 2 para 273. Em seguida, clique em Inicializar.

3. Configuração de pós-processamento e exportação de dados

  1. Em Resultados, clique em Contornos para abrir a janela Contornos. Na janela Contornos, ative as opções Preenchido, Valores de nó, Valores de limite, Intervalo global e Intervalo automático . Selecione Fases como Contornos do Tipo e selecione Fração de volume. Em seguida, escolha Fase 2 como Fase e clique em Salvar/Exibir.
  2. Em Atividades de Cálculo, clique em Animações de Solução para abrir a janela Definição de Animações.
  3. Na janela Definição de animações, defina Gravar após cada como 1 e escolha o intervalo de tempo. Selecione Arquivo HSF como o Tipo de armazenamento. Selecione o contorno-1 nas opções de Objeto de animação e clique em OK.
  4. Em Resultados, escolha a opção Superfícies e, em seguida, clique na opção Nova superfície de linha/inclinação para abrir a janela Superfície de linha/inclinação.
  5. Na janela Linha/Superfície de inclinação, ative Linha, defina x0 [m] como 0,21, defina x1 [m] como 0,21, defina y0 [m] como 0, defina y1 [m] como 0,015 e clique em Criar.
  6. Na guia Arquivo, escolha a opção Durante o cálculo em Exportar. Em seguida, clique em Dados da Solução para abrir a janela Exportação Automática.
  7. Na janela Exportação automática, escolha ASCII como o Tipo de arquivo. Selecione Centro da célula como local, selecione Espaço como delimitador. Defina Exportar dados a cada como 1 e escolha Intervalo de tempo. Na opção Superfícies, escolha linha-1. Na opção Quantidades, selecione Densidade (fase 2) e Fração de volume (fase 2). Em seguida, clique em Procurar, abra a janela Selecionar arquivo. Na janela Selecionar arquivo, clique em OK. Na janela Exportação automática, clique em OK.
  8. Em Atividades de Cálculo, clique em Salvamento Automático (Cada Tempo de Fluxo) para abrir a janela Salvamento Automático. Em seguida, defina Salvar arquivo de dados a cada [s] como 100, selecione Tempo de fluxo, selecione Somente se modificado como o tipo de arquivo de caso associado a Salvar e clique em OK.
  9. Na Solução, clique em Executar cálculo para abrir a janela Executar cálculo. Em seguida, selecione Fixo como o Tipo de Avanço de Tempo e Especificado pelo Usuário como o Método de Avanço de Tempo.
  10. Na janela Executar cálculo, defina Número de etapas de tempo como 7200, defina Tamanho da etapa de tempo como 1, defina Máximo de iterações/Etapa de tempo como 20, defina Intervalo de relatório como 1 e defina Intervalo de atualização de perfil como 2. Em seguida, clique em Calcular.

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Resultados

O modelo numérico modificado proposto captura efetivamente as principais características da formação de geada. Baseia-se nos mecanismos fundamentais subjacentes ao crescimento da geada e emprega uma abordagem de fluxo multifásico euleriana, juntamente com o modelo de mudança de fase de Lee. Essa abordagem permite que o modelo atualize a fração máxima do volume de gelo, levando em conta as variações de densidade ao longo do processo de congelamento. As simulações mostram que o modelo é rigorosamente validado comparando su...

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Discussão

Este estudo desenvolve um modelo numérico capaz de simular a formação de geada em uma superfície fria horizontal de baixa temperatura, ajustando dinamicamente o limite superior da fração volumétrica de acordo com o tempo e as condições de operação, reproduzindo variações na densidade da geada. Embora a validação apresentada aqui esteja confinada às temperaturas convencionais da superfície fria, o modelo garante que o desvio máximo entre a espessura da geada simulada e experimental não ex...

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Divulgações

Os autores declaram que não têm interesses financeiros concorrentes conhecidos ou relacionamentos pessoais que possam ter influenciado o trabalho relatado neste artigo.

Agradecimentos

Esta pesquisa é apoiada por (XLYC2203184), (U23A20657) e (LJ222410153082).

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Materiais

Lista de materiais utilizados neste artigo
NomeEmpresaNúmero de catálogoComentários
FluenteANSYS
ICEMSistema Dassault
Reivindicação espacialANSYS

Referências

  1. Saygin, A., Basol, A. M., Arik, M. An experimental study on the frost formation over a flat plate: Effect of frosting on heat transfer. Exp Therm Fluid Sci. 144, 110862(2023).
  2. Fang, X., et al. A new frictional pressure drop correlation based on flow patterns for hydrocarbon refrigerants condensation flow. Int J Refrig. 170, 214-223 (2025).
  3. Rong, X., et al. Experimental study on a multi-evaporator mutual defrosting system for air source heat pumps. Appl Energy. 332, 120528(2023).
  4. Jia, Y., Xu, X., Li, Y., Liang, X., Yao, M. Experimental studies on frost and defrost of fine tube bundles under coolant temperature between −20 and −5 °C. Int J Heat Mass Transf. 116, 617-620 (2018).
  5. Song, M., Dang, C. Review on the measurement and calculation of frost characteristics. Int. J. Heat Mass Transf. 124, 586-614 (2018).
  6. Lee, J., Lee, K. -S. The behavior of frost layer growth under conditions favorable for desublimation. Int J Heat Mass Transf. 120, 259-266 (2018).
  7. Lee, Y. B., Ro, S. T. Frost formation on a vertical plate in simultaneously developing flow. Exp Therm Fluid Sci. 26 (8), 939-945 (2002).
  8. Mengjie, S., Shangwen, L., Hosseini, S. H., Xiaoyan, L., Zhihua, W. An experimental study on the effect of horizontal cold plate surface temperature on frosting characteristics under natural convection. Appl Therm Eng. 211, 118416(2022).
  9. Jeong, H., Byun, S., Kim, D. R., Lee, K. S. Frost growth mechanism and its behavior under ultra-low temperature conditions. Int J Heat Mass Transf. 169, 120941(2021).
  10. Barzanoni, Y., Noorshams, O., Basirat Tabrizi, H., Damangir, E. Experimental investigation of frost formation on a horizontal cold cylinder under cross flow. Int J Refrig. 34 (4), 1174-1180 (2011).
  11. Jones, B. W., Parker, J. D. Frost formation with varying environmental parameters. J Heat Transf. 97 (2), 255-259 (1975).
  12. Cui, J., Li, W. Z., Liu, Y., Jiang, Z. Y. A new time- and space-dependent model for predicting frost formation. Appl Therm Eng. 31 (4), 447-457 (2011).
  13. Lenic, K., Trp, A., Frankovic, B. Transient two-dimensional model of frost formation on a fin-and-tube heat exchanger. Int J Heat Mass Transf. 52 (1-2), 22-32 (2009).
  14. Yu, J., Huo, R., Shen, H., Li, X., Zhu, Z. A simulation study on the condensation flow and thermal control characteristics of mixed refrigerant in a dimpled tube. Appl Therm Eng. 231, 120889(2023).
  15. Yu, J., Jiang, Y., Cai, W., Li, X., Zhu, Z. Condensation flow patterns and heat transfer correction for zeotropic hydrocarbon mixtures in a helically coiled tube. Int J Heat Mass Transf. 143, 118500(2019).
  16. You, Y., Wang, S., Lv, W., Chen, Y., Gross, U. A CFD model of frost formation based on dynamic meshes technique via secondary development of ANSYS Fluent. Int J Heat Fluid Flow. 89, 108807(2021).
  17. Cai, W., Fang, X., Li, S., Qiu, G. A modified CFD model for frosting on a horizontal plate. Int J Heat Mass Transf. 229, 125726(2024).
  18. Boyina, K. S., et al. Condensation frosting on meter-scale superhydrophobic and superhydrophilic heat exchangers. Int J Heat Mass Transf. 145, 118694(2019).

Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.

Reimpressões e permissões

Etiquetas

Forma o de GeadaModelo Num ricoEscoamento Multif sicoMudan a de Fase de LeeMorfologia da GeadaEspessura da GeadaDensidade da GeadaVaria o de DensidadeValida o Experimental