Commit 634a8fdf authored by Lluis Jofre Cruanyes's avatar Lluis Jofre Cruanyes
Browse files

Pressure accuracy improved

parent 00df5bf6

src/FlowSolverRHEA.cpp

+2 −0
Original line number Diff line number Diff line
@@ -922,6 +922,7 @@ void FlowSolverRHEA::calculateThermodynamicsFromPrimitiveVariables() { 
	            T = T_field[I1D(i,j,k)]; 

                    thermodynamics->calculateTemperatureFromPressureDensityWithInitialGuess( T, P_field[I1D(i,j,k)], rho_field[I1D(i,j,k)] );

                    T_field[I1D(i,j,k)]    = T;

		    P_field[I1D(i,j,k)]    = thermodynamics->calculatePressureFromTemperatureDensity( T_field[I1D(i,j,k)], rho_field[I1D(i,j,k)] );	/// Added to improve pressure accuracy

                    e  = thermodynamics->calculateInternalEnergyFromPressureTemperatureDensity( P_field[I1D(i,j,k)], T_field[I1D(i,j,k)], rho_field[I1D(i,j,k)] ); 

                    ke = 0.5*( pow( u_field[I1D(i,j,k)], 2.0 ) + pow( v_field[I1D(i,j,k)], 2.0 ) + pow( w_field[I1D(i,j,k)], 2.0 ) ); 

                    E_field[I1D(i,j,k)]    = e + ke;
@@ -962,6 +963,7 @@ void FlowSolverRHEA::calculateThermodynamicsFromPrimitiveVariables() { 
                    thermodynamics->calculatePressureTemperatureFromDensityInternalEnergy( P, T, rho_field[I1D(i,j,k)], e );

                    P_field[I1D(i,j,k)]   = P; 

                    T_field[I1D(i,j,k)]   = T; 

		    P_field[I1D(i,j,k)]   = thermodynamics->calculatePressureFromTemperatureDensity( T_field[I1D(i,j,k)], rho_field[I1D(i,j,k)] );	/// Added to improve pressure accuracy

                    sos_field[I1D(i,j,k)] = thermodynamics->calculateSoundSpeed( P_field[I1D(i,j,k)], T_field[I1D(i,j,k)], rho_field[I1D(i,j,k)] );

                    thermodynamics->calculateSpecificHeatCapacities( c_v, c_p, P_field[I1D(i,j,k)], T_field[I1D(i,j,k)], rho_field[I1D(i,j,k)] );

                    c_v_field[I1D(i,j,k)] = c_v;