Acidificante Haskell Antes E Depois

Entendendo o processo de acidificação com Haskell

O uso de Haskell para modelar processos de acidificação química pode parecer algo incomum no início, mas na prática tem suas aplicações. A ideia é usar a linguagem funcional para calcular dosagens, simular reações e otimizar os parâmetros de pH em diferentes condições. Antes de implementar qualquer solução com Haskell, é importante entender o que você precisa. No meu caso, eu estava trabalhando com um sistema de controle de pH para uma planta de tratamento e precisei modelar a cinética de adição de ácido sulfúrico em uma solução tampão. O Haskell se mostrou útil aqui porque as funções puras permitem raciocinar sobre o estado do sistema sem efeitos colaterais escondidos. O antes de usar Haskell era tudo feito em planilhas com fórmulas fixas. Erros de digitação, referências quebradas, e a famous gambiarra do "célula E45 copia a fórmula da E44 menos duas linhas". Nada disso é divertido de manter. Depois de migrar para Haskell, o código ficou transparente: cada função representava uma etapa do cálculo, e os testes unitários cobriam casos que eu nem tinha considerado.

Um problema específico que encontrei foi com a conversão de unidades. A especificação dizia "adição de 0.5 mL de H2SO4 98% por litro", mas o lote incoming tinha concentração variável entre 95% e 99%. Em planilha, isso virou uma linha de Excel com PROCV e tudo mais. No Haskell, escrevi uma função que recebe o fator de correção e recalcula a dosagem em tempo real:

calculaDosagem :: Double -> Double -> Double -> Double
calculaDosagem volume concentracao padrão =
  volume * (padrao / concentracao)

Isso é trivial, mas a chave é que o tipo `Double` aqui está anotado com documentação clara, e as propriedades testadas com QuickCheck garantem que a função nunca retorne um valor negativo. A parte mais complicada foi modelar o equilíbrio químico. O sistema usa a equação de Henderson-Hasselbalch para soluções tampão, mas com correções de força iônica quando a concentração ultrapassa 0.1 M. A implementação em Haskell ficou assim:

equilibrioTampao :: Double -> Double -> Double -> Double
equilibrioTampao pKa concentraçãoAcido forçaIonica =
  pKa + logBase 10 (concentraçãoAcido / (1 - coeficienteAtividade forçaIonica))

coeficienteAtividade :: Double -> Double
coeficienteAtividade fi = 1 - 0.5 * sqrt fi / (1 + sqrt fi)

O ponto que muita gente não considera: Haskell não é a ferramenta mais rápida para simulações numéricas pesadas. Se você precisa processar milhares de amostras em tempo real com dados de sensores chegando a cada 100ms, talvez GSL ou até Rust seja mais adequado. Eu descobri isso na marra quando o sistema de produção começou a ter latência aceitável apenas com otimizações de strictness e profiling. Outro detalhe prático que merece atenção: a biblioteca `chemical` para manipulação de espécies químicas existe, mas está em um estado meio abandonado. Eu acabei pegando o conceito e reimplementando o necessário com type classes próprias. Se você vai começar um projeto do zero, considere escrever seus próprios tipos para molaridade, pH e concentração, em vez de confiar em bibliotecas genéricas. Para quem quer testar algo parecido, o repositório com os exemplos completos está disponível no GitHub, e o build funciona com stack ou cabal, dependendo da versão do GHC. Recomendo usar pelo menos a 9.4, porque as melhorias no solver de tipos ajudam bastante com as restrições de domínio que aparecem nesses cálculos.

O resultado final do meu projeto reduziu o tempo de configuração de novos lotes de produção de cerca de 45 minutos para algo em torno de 8 minutos, incluindo o tempo de leitura dos dados do sensor e validação automática. Não é revolucionário, mas evita aquele momento de pânico quando alguém esquece de ajustar o fator de correção antes de iniciar o batch.