Thermoacoustic functions are important parameters of one-dimensional codes used for the design of thermoacoustic engines. The thermal and viscous thermoacoustic functions allow the inclusion of three dimensional effects in one-dimensional codes. These functions are especially important in the regenerator of a thermoacoustic engine, where the thermoacoustic heat pumping occurs. Even though analytical solutions were derived for uniform pores, the thermoacoustic functions for complex geometries such as stacked screen or random fiber regenerators cannot be calculated analytically. In order to gain more insight into the geometry induced complex flow fields, the procedure of Udea et al. (2009) to estimate the thermoacoustic functions was applied in computational fluid-dynamic simulations. By using two measurement locations outside of the regenerator and modeling the regenerator as an array of uniform pores, it is possible to estimate the thermoacoustic functions for complex geometries. Furthermore, a correction method is proposed to quantify the entrance effects at the beginning and end of a regular pore. The simulations are first validated for a uniform cylindrical pore with the help of the analytical solution. Then the correction method is successfully applied to a cylindrical pore with the results closely matching the analytical solution. Finally, the method is applied to the model of a staggered, stacked screen regenerator.