pressure.h

Go to the documentation of this file.

/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
 * others.  For conditions of distribution and use see copyright notice in license.txt */
 
#ifndef PRESSURE_H_
#define PRESSURE_H_
 
#include "rt.h"
#include "rfield.h"
#include "doppvel.h"
#include "physconst.h"
#include "transition.h"
 
void PresTotCurrent(void);
 
inline double PressureRadiationLine( const TransitionProxy &t, realnum DopplerWidth )
{
        DEBUG_ENTRY( "PressureRadiationLine()" );
 
        /* return zero if below plasma frequency */     
        if( t.EnergyErg() / EN1RYD <= rfield.plsfrq )
                return 0.;
 
        /* RT_LineWidth gets line width in terms of RT effects */
        double width = RT_LineWidth(t, DopplerWidth);
 
        double PopOpc = t.Emis().PopOpc()/(*t.Lo()).g();
        /* return zero line radiation PressureReturned if line mases or has
         * zero opacity */
        /* \todo 2 1e-22 is arbtrary but roughly 1/kpc.  Replace with a cloud width if available? */
        if( PopOpc*t.Emis().opacity()/ DopplerWidth <= 1.e-22 || width<=0. )
                return 0.;
 
        double PressureReturned = PI8 * HPLANCK / 3. * POW4(t.EnergyWN()) *
                ((*t.Hi()).Pop()/(*t.Hi()).g())/PopOpc * width;
 
        /* this prevents line radiation PressureReturned from being very large when line 
         * is not optically thick but total opacity at that energy is large 
         * due to overlapping transitions */
        long int ipLineCenter = t.Emis().ipFine() + rfield.ipFineConVelShift;
        if( ipLineCenter > 0 && ipLineCenter < rfield.nfine && rfield.lgOpacityFine &&
                rfield.fine_opac_zone[ipLineCenter] > SMALLFLOAT )
        {
                double FractionThisLine = t.Emis().PopOpc() * t.Emis().opacity() / DopplerWidth/
                        rfield.fine_opac_zone[ipLineCenter];
                if( FractionThisLine<1e-5 )
                        FractionThisLine = 0.;
                /* fine opacities are only reevaluated one time per zone due
                 * to the expense - PopOpc is for the current solution - but the two
                 * may be out of step by a few percent, due to the variation in
                 * abundance from zone to zone.  This prevents the change
                 * in solution from increasing the radiation pressure.
                 * This correction is mainly an order of magnitude scaler to prevent
                 * optically thin lines from appearing to be optically thick due to
                 * overlapping lines */
                FractionThisLine = MIN2(1., FractionThisLine);
                ASSERT( FractionThisLine >= 0. && FractionThisLine <= 1.0 );
                PressureReturned *= FractionThisLine;
        }
 
        return PressureReturned;
}
 
struct t_pressure {
 
        realnum PresLow, 
          PresHigh;
 
        realnum PresPowerlaw;
 
        double PresRamCurr;
 
        double PresTurbCurr;
 
        double PresTotlCurr, 
          PresTotlError,
          PresGasCurr;
 
        double PresTotlInit;
 
        bool lgPressureInitialSpecified;
        double PressureInitialSpecified;
 
        double pres_radiation_lines_curr;
 
        bool lgContRadPresOn;
 
        realnum PresInteg, 
          pinzon;
 
        realnum PresIntegElecThin,
                pinzon_PresIntegElecThin;
 
        double RhoGravity_dark;
        double RhoGravity_self;
        double RhoGravity_external;
        double RhoGravity;
        double IntegRhoGravity;
        int gravity_symmetry;
        double self_mass_factor;
        vector<double> external_mass[3];
 
        bool lgPres_radiation_ON;
        bool lgPres_magnetic_ON;
        bool lgPres_ram_ON;
 
        realnum 
          RadBetaMax, 
          pbeta, 
          PresMax;
 
        long int ipPradMax_line;
 
        long int ipPradMax_nzone;
 
        char chLineRadPres[101];
 
        bool lgPradCap, 
          lgPradDen;
 
        bool lgLineRadPresOn;
 
        bool lgRadPresAbortOK;
 
        bool lgSonicPointAbortOK; 
 
        bool lgSonicPoint;
 
        bool lgStrongDLimbo;
 
        };
extern t_pressure pressure;
 
 
#endif /* PRESSURE_H_ */