Numerical modeling of thin‐film bulk acoustic wave resonators using a Crank‐Nicolson finite‐difference time‐domain method