#include <ERF_Fitch.H>

Inheritance diagram for Fitch:

Collaboration diagram for Fitch:

[legend]

Public Member Functions
	Fitch ()

virtual	~Fitch ()=default

void	advance (const amrex::Geometry &geom, const amrex::Real &dt_advance, amrex::MultiFab &cons_in, amrex::MultiFab &mf_vars_fitch, amrex::MultiFab &U_old, amrex::MultiFab &V_old, amrex::MultiFab &W_old, const amrex::MultiFab &mf_Nturb, const amrex::MultiFab &mf_SMark, const amrex::Real &time) override

void	source_terms_cellcentered (const amrex::Geometry &geom, const amrex::MultiFab &cons_in, amrex::MultiFab &mf_vars_ewp, const amrex::MultiFab &U_old, const amrex::MultiFab &V_old, const amrex::MultiFab &W_old, const amrex::MultiFab &mf_Nturb)

void	update (const amrex::Real &dt_advance, amrex::MultiFab &cons_in, amrex::MultiFab &U_old, amrex::MultiFab &V_old, const amrex::MultiFab &mf_vars_fitch)

Public Member Functions inherited from NullWindFarm
	NullWindFarm ()

virtual	~NullWindFarm ()=default

virtual void	set_turb_spec (const amrex::Real &rotor_rad, const amrex::Real &hub_height, const amrex::Real &thrust_coeff_standing, const amrex::Vector< amrex::Real > &wind_speed, const amrex::Vector< amrex::Real > &thrust_coeff, const amrex::Vector< amrex::Real > &power)

virtual void	set_turb_loc (const amrex::Vector< amrex::Real > &xloc, const amrex::Vector< amrex::Real > &yloc)

virtual void	set_turb_disk_angle (const amrex::Real &turb_disk_angle)

virtual void	set_blade_spec (const amrex::Vector< amrex::Real > &bld_rad_loc, const amrex::Vector< amrex::Real > &bld_twist, const amrex::Vector< amrex::Real > &bld_chord)

virtual void	set_blade_airfoil_spec (const amrex::Vector< amrex::Vector< amrex::Real >> &bld_airfoil_aoa, const amrex::Vector< amrex::Vector< amrex::Real >> &bld_airfoil_Cl, const amrex::Vector< amrex::Vector< amrex::Real >> &bld_airfoil_Cd)

virtual void	set_turb_spec_extra (const amrex::Vector< amrex::Real > &velocity, const amrex::Vector< amrex::Real > &C_P, const amrex::Vector< amrex::Real > &C_T, const amrex::Vector< amrex::Real > &rotor_RPM, const amrex::Vector< amrex::Real > &blade_pitch)

void	get_turb_spec (amrex::Real &rotor_rad, amrex::Real &hub_height, amrex::Real &thrust_coeff_standing, amrex::Vector< amrex::Real > &wind_speed, amrex::Vector< amrex::Real > &thrust_coeff, amrex::Vector< amrex::Real > &power)

void	get_turb_loc (amrex::Vector< amrex::Real > &xloc, amrex::Vector< amrex::Real > &yloc)

void	get_turb_disk_angle (amrex::Real &turb_disk_angle)

void	get_blade_spec (amrex::Vector< amrex::Real > &bld_rad_loc, amrex::Vector< amrex::Real > &bld_twist, amrex::Vector< amrex::Real > &bld_chord)

void	get_blade_airfoil_spec (amrex::Vector< amrex::Vector< amrex::Real >> &bld_airfoil_aoa, amrex::Vector< amrex::Vector< amrex::Real >> &bld_airfoil_Cl, amrex::Vector< amrex::Vector< amrex::Real >> &bld_airfoil_Cd)

void	get_turb_spec_extra (amrex::Vector< amrex::Real > &velocity, amrex::Vector< amrex::Real > &C_P, amrex::Vector< amrex::Real > &C_T, amrex::Vector< amrex::Real > &rotor_RPM, amrex::Vector< amrex::Real > &blade_pitch)

Protected Attributes
amrex::Vector< amrex::Real >	xloc

amrex::Vector< amrex::Real >	yloc

amrex::Real	hub_height

amrex::Real	rotor_rad

amrex::Real	thrust_coeff_standing

amrex::Real	nominal_power

amrex::Vector< amrex::Real >	wind_speed

amrex::Vector< amrex::Real >	thrust_coeff

amrex::Vector< amrex::Real >	power

Protected Attributes inherited from NullWindFarm
amrex::Vector< amrex::Real >	m_xloc

amrex::Vector< amrex::Real >	m_yloc

amrex::Real	m_turb_disk_angle

amrex::Real	m_hub_height

amrex::Real	m_rotor_rad

amrex::Real	m_thrust_coeff_standing

amrex::Real	m_nominal_power

amrex::Vector< amrex::Real >	m_wind_speed

amrex::Vector< amrex::Real >	m_thrust_coeff

amrex::Vector< amrex::Real >	m_power

amrex::Vector< amrex::Real >	m_bld_rad_loc

amrex::Vector< amrex::Real >	m_bld_twist

amrex::Vector< amrex::Real >	m_bld_chord

amrex::Vector< amrex::Vector< amrex::Real > >	m_bld_airfoil_aoa

amrex::Vector< amrex::Vector< amrex::Real > >	m_bld_airfoil_Cl

amrex::Vector< amrex::Vector< amrex::Real > >	m_bld_airfoil_Cd

amrex::Vector< amrex::Real >	m_velocity

amrex::Vector< amrex::Real >	m_C_P

amrex::Vector< amrex::Real >	m_C_T

amrex::Vector< amrex::Real >	m_rotor_RPM

amrex::Vector< amrex::Real >	m_blade_pitch

Additional Inherited Members
Static Public Member Functions inherited from NullWindFarm
static AMREX_GPU_DEVICE bool	find_if_marked (amrex::Real x1, amrex::Real x2, amrex::Real y1, amrex::Real y2, amrex::Real x0, amrex::Real y0, amrex::Real nx, amrex::Real ny, amrex::Real d_hub_height, amrex::Real d_rotor_rad, amrex::Real z)

Constructor & Destructor Documentation

◆ Fitch()

Fitch::Fitch ( )

inline

13 {}

◆ ~Fitch()

virtual Fitch::~Fitch ( )

virtualdefault

Member Function Documentation

◆ advance()

void Fitch::advance	(	const amrex::Geometry &	geom,
		const amrex::Real &	dt_advance,
		amrex::MultiFab &	cons_in,
		amrex::MultiFab &	mf_vars_fitch,
		amrex::MultiFab &	U_old,
		amrex::MultiFab &	V_old,
		amrex::MultiFab &	W_old,
		const amrex::MultiFab &	mf_Nturb,
		const amrex::MultiFab &	mf_SMark,
		const amrex::Real &	time
	)

overridevirtual

Implements NullWindFarm.

 {
     AMREX_ALWAYS_ASSERT(W_old.nComp() > 0);
     AMREX_ALWAYS_ASSERT(mf_SMark.nComp() > 0);
     AMREX_ALWAYS_ASSERT(time > -1.0);
     source_terms_cellcentered(geom, cons_in, mf_vars_fitch, U_old, V_old, W_old, mf_Nturb);
     update(dt_advance, cons_in, U_old, V_old, mf_vars_fitch);
 }

◆ source_terms_cellcentered()

void Fitch::source_terms_cellcentered	(	const amrex::Geometry &	geom,
		const amrex::MultiFab &	cons_in,
		amrex::MultiFab &	mf_vars_ewp,
		const amrex::MultiFab &	U_old,
		const amrex::MultiFab &	V_old,
		const amrex::MultiFab &	W_old,
		const amrex::MultiFab &	mf_Nturb
	)

 {
  
   get_turb_spec(rotor_rad, hub_height, thrust_coeff_standing,
                   wind_speed, thrust_coeff, power);
  
   auto dx = geom.CellSizeArray();
  
   // Domain valid box
   const amrex::Box& domain = geom.Domain();
   int domlo_z = domain.smallEnd(2);
   int domhi_z = domain.bigEnd(2) + 1;
  
  
   //Real sum = 0.0;
   //Real *sum_area = &sum;
  
   // The order of variables are - Vabs dVabsdt, dudt, dvdt, dTKEdt
   mf_vars_fitch.setVal(0.0);
   Real d_hub_height = hub_height;
   Real d_rotor_rad = rotor_rad;
   Gpu::DeviceVector<Real> d_wind_speed(wind_speed.size());
   Gpu::DeviceVector<Real> d_thrust_coeff(thrust_coeff.size());
  
   // Copy data from host vectors to device vectors
   Gpu::copy(Gpu::hostToDevice, wind_speed.begin(), wind_speed.end(), d_wind_speed.begin());
   Gpu::copy(Gpu::hostToDevice, thrust_coeff.begin(), thrust_coeff.end(), d_thrust_coeff.begin());
  
   for ( MFIter mfi(cons_in,TilingIfNotGPU()); mfi.isValid(); ++mfi) {
  
         const Box& gbx = mfi.growntilebox(1);
         auto fitch_array = mf_vars_fitch.array(mfi);
         auto Nturb_array = mf_Nturb.array(mfi);
         auto u_vel       = U_old.array(mfi);
         auto v_vel       = V_old.array(mfi);
         auto w_vel       = W_old.array(mfi);
  
         const Real* wind_speed_d     = d_wind_speed.dataPtr();
         const Real* thrust_coeff_d   = d_thrust_coeff.dataPtr();
         const int n_spec_table = d_wind_speed.size();
  
         ParallelFor(gbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
             int kk = amrex::min(amrex::max(k, domlo_z), domhi_z);
  
             Real z_k   = kk*dx[2];
             Real z_kp1 = (kk+1)*dx[2];
  
             Real A_ijk = compute_Aijk(z_k, z_kp1, d_hub_height, d_rotor_rad);
  
             // Compute Fitch source terms
  
             Real Vabs = std::pow(u_vel(i,j,k)*u_vel(i,j,k) +
                                  v_vel(i,j,k)*v_vel(i,j,k) +
                                  w_vel(i,j,kk)*w_vel(i,j,kk), 0.5);
  
             Real C_T = interpolate_1d(wind_speed_d, thrust_coeff_d, Vabs, n_spec_table);
             Real C_TKE = 0.0;
  
             fitch_array(i,j,k,0) = Vabs;
             fitch_array(i,j,k,1) =  -0.5*Nturb_array(i,j,k)/(dx[0]*dx[1])*C_T*Vabs*Vabs*A_ijk/(z_kp1 - z_k);
             fitch_array(i,j,k,2) = u_vel(i,j,k)/Vabs*fitch_array(i,j,k,1);
             fitch_array(i,j,k,3) = v_vel(i,j,k)/Vabs*fitch_array(i,j,k,1);
             fitch_array(i,j,k,4) = 0.5*Nturb_array(i,j,k)/(dx[0]*dx[1])*C_TKE*std::pow(Vabs,3)*A_ijk/(z_kp1 - z_k);
  
                  //amrex::Gpu::Atomic::Add(sum_area, A_ijk);
         });
     }
         //std::cout << "Checking sum here...." <<"\n";
         //printf("%0.15g, %0.15g\n", *sum_area , PI*R*R);
         //exit(0);
 }

Here is the call graph for this function:

◆ update()

void Fitch::update	(	const amrex::Real &	dt_advance,
		amrex::MultiFab &	cons_in,
		amrex::MultiFab &	U_old,
		amrex::MultiFab &	V_old,
		const amrex::MultiFab &	mf_vars_fitch
	)

 {
  
     for ( MFIter mfi(cons_in,TilingIfNotGPU()); mfi.isValid(); ++mfi) {
  
         Box bx  = mfi.tilebox();
         Box tbx = mfi.nodaltilebox(0);
         Box tby = mfi.nodaltilebox(1);
  
         auto cons_array  = cons_in.array(mfi);
         auto fitch_array = mf_vars_fitch.array(mfi);
         auto u_vel       = U_old.array(mfi);
         auto v_vel       = V_old.array(mfi);
  
         ParallelFor(tbx, tby, bx,
         [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
         {
             u_vel(i,j,k) = u_vel(i,j,k) + (fitch_array(i-1,j,k,2) + fitch_array(i,j,k,2))/2.0*dt_advance;
         },
         [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
         {
             v_vel(i,j,k) = v_vel(i,j,k) + (fitch_array(i,j-1,k,3) + fitch_array(i,j,k,3))/2.0*dt_advance;
         },
         [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
         {
             cons_array(i,j,k,RhoKE_comp) = cons_array(i,j,k,RhoKE_comp) + fitch_array(i,j,k,4)*dt_advance;
         });
     }
 }

Member Data Documentation

◆ hub_height

amrex::Real Fitch::hub_height

protected

◆ nominal_power

amrex::Real Fitch::nominal_power

protected

◆ power

amrex::Vector<amrex::Real> Fitch::power

protected

◆ rotor_rad

amrex::Real Fitch::rotor_rad

protected

◆ thrust_coeff

amrex::Vector<amrex::Real> Fitch::thrust_coeff

protected

◆ thrust_coeff_standing

amrex::Real Fitch::thrust_coeff_standing

protected

◆ wind_speed

amrex::Vector<amrex::Real> Fitch::wind_speed

protected

◆ xloc

amrex::Vector<amrex::Real> Fitch::xloc

protected

◆ yloc

amrex::Vector<amrex::Real> Fitch::yloc

protected

The documentation for this class was generated from the following files:

Source/WindFarmParametrization/Fitch/ERF_Fitch.H
Source/WindFarmParametrization/Fitch/ERF_AdvanceFitch.cpp

Public Member Functions

Protected Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ Fitch()

◆ ~Fitch()

Member Function Documentation

◆ advance()

◆ source_terms_cellcentered()

◆ update()

Member Data Documentation

◆ hub_height

◆ nominal_power

◆ power

◆ rotor_rad

◆ thrust_coeff

◆ thrust_coeff_standing

◆ wind_speed

◆ xloc

◆ yloc