We present a methodology to simulate the quantum thermodynamics of thermal machines which are built from an interacting working medium in contact with fermionic reservoirs at a fixed temperature and chemical potential. Our method works at a finite temperature, beyond linear response and weak system-reservoir coupling, and allows for nonquadratic interactions in the working medium. The method uses mesoscopic reservoirs, continuously damped toward thermal equilibrium, in order to represent continuum baths and a novel tensor-network algorithm to simulate the steady-state thermodynamics. Using the example of a quantum-dot heat engine, we demonstrate that our technique replicates the well-known Landauer-Büttiker theory for efficiency and power. We then go beyond the quadratic limit to demonstrate the capability of our method by simulating a three-site machine with nonquadratic interactions. Remarkably, we find that such interactions lead to power enhancement, without being detrimental to the efficiency. Furthermore, we demonstrate the capability of our method to tackle complex many-body systems by extracting the superdiffusive exponent for high-temperature transport in the isotropic Heisenberg model. Finally, we discuss transport in the gapless phase of the anisotropic Heisenberg model at a finite temperature and its connection to charge conjugation parity, going beyond the predictions of single-site boundary driving configurations.14 MoreReceived 7 January 2020Revised 1 June 2020Accepted 7 July 2020DOI:https://doi.org/10.1103/PhysRevX.10.031040Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasNonequilibrium statistical mechanicsQuantum thermodynamicsTransport phenomenaQuantum InformationCondensed Matter & Materials PhysicsAtomic, Molecular & OpticalStatistical Physics