Las pilas de lixiviación son comunes en la industria minera. En su diseño y gestión, comúnmente no se tienen en cuenta criterios basados en procesos hidrogeomecánicos-geoquímicos acoplados. Sin embargo, estos son de suma importancia por dos razones principales: (1) podrían desempeñar un papel importante en la estabilidad física de la almohadilla; y (2) son fundamentales para la eficiencia metalúrgica, especialmente a mediano y largo plazo. Como ejemplo de la relevancia de los fenómenos acoplados complejos, se puede mencionar que la estabilidad geotécnica de la pila está relacionada con la distribución de las presiones de los líquidos que, a su vez, está relacionada con la permeabilidad a través de la ley de Darcy. A su vez, la permeabilidad se verá afectada por los cambios de porosidad, que dependen del estado de consolidación y de la disolución/precipitación de los minerales. Obviamente, el rendimiento hidrometalúrgico del mineral de la pila también está fuertemente relacionado con los procesos acoplados mencionados anteriormente. En este trabajo, presentamos un modelo numérico para resolver los procesos acoplados hidro-geomecánico-geoquímicos, desarrollado con COMSOL-PHREEQC, una interfaz Java desarrollada para combinar COMSOL Multiphysics, una conocida plataforma FEM comercial, con PHREEQC, el simulador geoquímico de software abierto desarrollado en USGS y utilizado en todo el mundo. La herramienta numérica resultante es capaz de resolver simultáneamente el flujo de líquido insaturado, la tensión efectiva del suelo y la disolución geoquímica y la precipitación de SOLUCIONES DE LIXIVIACIÓN EN PILAS, 2014  LIMA, PERÚ 274. Los cambios geoquímicos que afectan la porosidad y la permeabilidad se acoplan explícitamente en los cálculos hidráulicos y geomecánicos.