Albany is an implicit, unstructured grid, finite element code for the solution and analysis of partial differential equations. Albany integrates the FASTMath adaptive unstructured mesh capabilites with the linear and nonlinear solvers, preconditioners, continuation, and time integration methods available in Trilinos, directly supporting a host of physics capabilities through the use of template-based generic programming. The Albany test suite contains 230+ examples of problems that include thermal, fluids, and structural mechanics applications that demonstrate a wide variety of models and implementations that include coupled problems, plasticity, Cahn-Hilliard solutions, poro-elasticity, and Poisson-Schrodinger equations. Albany demonstrates how to combine these capabilities to solve problems on parallel machines and illustrates the use of automatic differentiation to support model derivatives, construction of the Jacobian operator for Newton’s method, use of embedded uncertainty quantification, and several other useful software construction techniques.