Cantera: cálculos de ignição homogênea

Quando usar

Use para um reator de gás ideal, fechado, adiabático e homogêneo com um mecanismo cinético conhecido, temperatura inicial, pressão e composição molar. O utilitário incluído executa casos de volume constante e pressão constante e relata um atraso precisamente definido baseado na temperatura. Não é um resolvedor de chamas ou um construtor geral de rede de reatores.

Um cálculo que é concluído com sucesso estabelece a execução numérica, não a validade do mecanismo para o combustível, pressão, temperatura, diluente ou observável de ignição medido. Leia references/interpretation.md ao escolher um mecanismo, comparar experimentos ou interpretar ignição não resolvida/em dois estágios.

Fluxo de trabalho

1. Identifique o mecanismo e seu intervalo de condições validadas. Preserve sua origem, versão, citação e quaisquer modificações. Verifique se sua fase é ideal-gas e se cada reagente, diluente e espécie rastreada existe. Para YAML personalizado com importações, mantenha os arquivos de dependência originais, bem como o instantâneo de fase gerado. Extensões de taxa Python personalizadas também precisam de seu código original e ambiente para repetição.
2. Escolha volume constante ou pressão constante do experimento físico. Forneça K, Pa, segundos e quantidades molares explicitamente. mole_amounts é normalizado para frações molares; não é um mapeamento de fração de massa. O relatório inclui a composição inicial normalizada.
3. Copie assets/hydrogen-ignition.json e altere suas condições. O caso H2/O2/Ar fornecido usa o h2o2.yaml incluído do Cantera para um exemplo numérico executável; não é uma recomendação para todos os experimentos de hidrogênio.
4. Defina um horizonte de tempo longo o suficiente para observar o aumento da temperatura e o declínio do pico da taxa de aquecimento. Escolha um espaçamento de saída fino o suficiente para localizar esse pico. Defina um aumento mínimo de temperatura para distinguir a ignição de aquecimento insignificante ou ruído numérico.
5. Execute o utilitário e inspecione todos os quatro históricos e o relatório. Refine novamente se o atraso mudar materialmente, se o máximo se aproximar de um limite de tempo, ou se a conservação falhar. Compare os históricos de temperatura e espécies rastreadas com a definição real de ignição.
6. Relate o conjunto de condições, hash do mecanismo, restrição do reator, definição de atraso, espaçamento de saída, alterações numéricas e limites científicos juntamente com o atraso.

Execute o exemplo testado

Da raiz da coleção:

uv run --no-project --python 3.12 --with cantera==3.2.0 --with numpy==2.5.3 \
  python skills/cantera/scripts/ignition_delay.py \
  skills/cantera/assets/hydrogen-ignition.json hydrogen-result

Testado em Python 3.12, Cantera 3.2.0 e NumPy 2.5.3. Nenhum executável de solver externo ou credenciais são necessários. Caminhos de mecanismo relativos locais são resolvidos contra o diretório do arquivo de configuração antes da pesquisa de dados integrada do Cantera. Use um novo diretório de saída a cada execução.

O exemplo de volume constante de 1000 K, 101325 Pa, H2:O2:Ar = 2:1:7 dá cerca de 0,313 ms usando a definição max(dT/dt) declarada. Em 3 ms sua temperatura é de aproximadamente 2920,67 K e concorda com um cálculo de equilíbrio UV separado. Estes são valores de regressão do pacote, não dados de validação experimental.

Contrato de atraso exato e refinamento

O atraso é o tempo do máximo global de numpy.gradient(T, time, edge_order=2) em uma grade de saída uniforme. Ele é relatado apenas se o aumento máximo de temperatura atingir minimum_temperature_rise_k e o máximo estiver a pelo menos dois índices de amostra de cada borda. Caso contrário, delay_s é nulo e um status explica o motivo. Nenhum atraso além do horizonte da simulação é extrapolado.

O utilitário usa explicitamente o clone=True do Cantera 3.2 e lê propriedades em evolução de reactor.phase. O Solution original mantém o estado inicial; não o leia como o estado final do reator. ReactorNet.advance(t) solicita um tempo absoluto, e nenhum limite de avanço é configurado, então a grade de saída permanece uniforme.

Ele executa quatro reatores frescos independentes:

numerically_resolved exige que todas as execuções produzam atrasos, com mudanças relativas de atraso dentro de

…