For our suite of simulations exploring the effects of metallicity, we fix all of the simulation input parameters except for the metal content (represented by varying the line strength parameter \(\bar{Q}=\) 400, 700, 1000, 1300, 1700, 2000). For the initial condition of each model, we relax a smooth Sobolev solution with the corresponding \(\bar{Q}\) value, choosing an inflow density which matches the predicted CAK mass loss rate. In Table X we show the resulting mass loss rates under both the smooth Sobolev solution and under the SSF radiative force, across all models.

Figure XX shows simulation snapshots of the radial variation of velocity and density, corresponding to galactic (\(\bar{Q} = 2000\)), LMC (\(\bar{Q} = 1000\)), and SMC (\(\bar{Q} = 400\)) metal content. Each snapshot was taken after \(1.5 \times 10^6\) seconds of simulation, so transients from the initial conditions have fully propagated out of the flow. We see a decreasing mass loss rate as predicted by theory, which manifests as reductions in both the velocity and density.