Before starting using the modified version of CLOUDY and the code that computes the line profiles it is necessary to test them.
Bertout (1984) computed line profiles for spherically symmetric stellar envelopes expanding with a velocity dependence with radius given by the power law
where
is the wind velocity at the stellar surface. The results
obtained for a wind velocity of this form and for a constant turbulent
velocity, with the procedure described in the previous section are
displayed in Figure 6.1. The stellar
parameters used were those typical of a T Tauri star and referred to
above, i.e.
,
and Teff = 4000 K. The
various
profiles in Figure 6.1 are labeled by the ratio of the
maximum velocity in the flow to the turbulent velocity and by the
exponent
from equation 6.9.
Figure 6.1: Line profiles obtained by the procedure
outlined above for cases studied by Bertout (1984).
is the exponent in the velocity power law
(see text). The numbers labeling each line correspond to the ratio of the
maximum wind velocity to the (constant) turbulent velocity in the model.
Comparing Figure 6.1 with figures 1, 4 and 5 in
Bertout (1984), reproduced here in Figure 6.2, one
sees that there is a good agreement
between the results obtained with the implementation described above and
the results from Bertout's work. Note that an exact match between the
line profiles presented in Bertout (1984) and those presented
here would be difficult to achieve. This is mainly due to the fact that
the two models require a different set of input parameters. As an example, a
parameter in Bertout's work is
, which is defined as
, where h, c and
have
the usual meaning,
is the radius of the star and
and
are respectively the velocity at
and the lower-level number
density at
. Such a parameter is not directly controled in the
modeling presented here. The aim is therefore to obtain the
same general shape for the line profiles as that in Bertout
bertoutlp, as well as the same sort of behaviour when varying
and the ratio between the maximum velocity in the flow and the
turbulent velocity. This is achieved with the performed tests, indicating
that the implementation of the method outlined in Section 6.1
is producing the expected results and that it is therefore operational.
Figure 6.2: Line profiles results from
Bertout (1984). Top left panel -
and vwind(max) to vturb
ratios of 3, 5, 10, 30 and 50; Top right -
and vwind(max) to vturb ratios
of 3, 6, 9 and 15; Bottom panel -
and vwind(max) to vturb ratios
of 3, 10 and 30.