# 2.1. Version 5.28 (latest - 2026/01/23)

New Features

  - Added a bulk node generalized natural convection BC for the porous enthalpy equation.

  - Added a volume weighted predictor corrector error for improved timestep sizes when using adaptivity and other mesh modification methods. To enable this, put `predictor-corrector weighting = volume_weighted` in the transient time block in your input deck.

  - Added new time-interpolation mode for MPMD coupling with SPARC when Aria is sub-stepping in time

  - Extended the `arialog_parse` script to work on log files containing solution control-driven nonlinear iteration loops

  - Added new configurations to the `teko_multiphysics` preset solver for automatically handling mesh displacement equations in ablation problems

  - Added new configurations to the `teko_multiphysics` preset solver for automatically splitting bulk nodes into a separate sub-block, improving performance and memory use with ILU preconditioners on the thermal solve

  - Memory usage improvements in `teko_multiphysics` preset solver: Reduce GMRES restart size from 300 to 30 (540 vector allocations saved) and avoid storing memory for adaptive subblock solvers that are no longer in use

Bug Fixes

  - Fixed bug in tftk_adapt where the log file would output the wrong element and node count before/after refinement

  - Fixed a bug where tftk_adapt would throw an unhelpful error when doing partial refinement on hex/quad meshes

  - Fixed a bug where the predictor corrector error was computed incorrectly for cylindrical coordinate systems

  - Fixed an issue when trying to apply contact enforcement to “sub-indexed” equations and disallowed using integers as species names

  - Fixed a contact enforcement bug when applying contact for multiple phases between a multiphase block and a single phase block

  - Fixed the mass-preserving density model to allow density-dependent properties to be used on sidesets