Study of dissipative quantum phase transitions in the Ohmic spin-boson model is numerically challenging in a dense limit of environmental modes. In this work, large-scale numerical simulations are carried out based on the variational principle. The validity of variational calculations, spontaneous breakdown of symmetries, and quantum fluctuations and correlations in the Ohmic bath are carefully analyzed, and the critical coupling as well as exponents are accurately determined in the weak tunneling and continuum limits. In addition, quantum criticality of the Ohmic bath is uncovered both in the delocalized phase and at the transition point.