Scalar electrodynamics can be used to investigate the formation of cosmic strings in the early universe. We present the results of lattice Monte Carlo simulations of an effective three-dimensional U(1)+Higgs theory that describes the equilibrium properties of finite-temperature scalar electrodynamics near the transition. A gauge-invariant criterion for the existence of a vortex is used in measuring the properties of the vortex network in the equilibrium state both in the Coulomb and in the Higgs phase of the system. The naive definition of the vortex density becomes meaningless in the continuum limit and special care is needed in extracting physical quantities. Numerical evidence for a physical discontinuity in the vortex density is given.

We examine the first-order cosmological QCD phase transition for a large class of parameter values, previously considered unlikely. We find that the hadron bubbles can nucleate at very large distance scales, they can grow as detonations as well as deflagrations, and that the phase transition may be completed without reheating to the critical temperature. For a subset of the parameter values studied, the inhomogeneities generated at the QCD phase transition might have a noticeable effect on nucleosynthesis.

We study how bubbles grow after the initial nucleation event in generic first-order cosmological phase transitions characterised by the values of latent heat, interface tension and correlation length, and driven by a scalar order parameter $\phi$. Equations coupling $\phi$ and the fluid variables $v$ and $T$ and depending on a dissipative constant $\Gamma$ are derived and solved numerically in the 1+1 dimensional case starting from a slightly deformed critical bubble configuration. Parameters corresponding to QCD and electroweak phase transitions are chosen and the whole history of the bubble with formation of combustion and shock fronts is computed as a function of $\Gamma$. Both deflagrations and detonations can appear depending on the values of the parameters. Reheating due to collisions of bubbles is also computed.