100
6 Ratchet Dimer Brownian Motor with Hydrodynamic Interactions
C EFFICIENCY: STOKESEFF, ENERGCONVEFF, KINETICEFF
c ========================================================
c Spiechowicz, Hänggi, Luczka, 2015-Int. Conf. on Noise
c DT = eq. =kT/m = gamma*DEFF/mass
C ========================================================
c VCX2EQ = 19.56142
c VCX2EQ = gamma*DEFFX/mass
STOKESEFF = mvcx**2/(mvcx2 - VCX2EQ)
ENERGCONVEFF = abs(Fload*mvcx)/(gamma*(mvcx2 - VCX2EQ))
KINETICEFF = w*mvcx**2/(4*pi*gamma*(mvcx2 - VCX2EQ))
CUBEROEFF = abs(Fload*mvcx)/(abs(Fload) + gamma*(mvcx2 - VCX2EQ))
SUZUKIEFF = (abs(Fload)*mvcx + gamma*mvcx2)/(-sumFrat)
c ========================================================
WRITE(*,’(” Fload =”,f10.7)’) Fload
WRITE(*,’(” gamma =”,f10.7)’) gamma
WRITE(*,’(” mass =”,f10.7)’) mass
WRITE(*,’(” Fo =”,f10.7)’) Fo
WRITE(*,’(” w =”,f10.7)’) w
WRITE(*,’(” =”,f10.7)’) MVCX
c WRITE(*,’(” VC =”,f10.7)’) VC
c WRITE(*,’(” VCX/VC =”,f10.7)’) VCX/VC
WRITE(*,’(” =”,f10.2)’) MVCX2
WRITE(*,’(”eq=”,E12.5)’) VCX2EQ
WRITE(*,’(” - eq =”,f10.7)’) (mvcx2 - vcx2eq)
WRITE(*,’(” sumFrat =”,f10.2)’) sumFrat
WRITE(*,’(” !Fload!* =”,f10.7)’) abs(Fload)*mvcx
WRITE(*,’(” gamma* =”,f10.7)’) gamma*mvcx2
WRITE(*,’(”RMEAN =”,f10.5)’) RMEAN
WRITE(*,’(”RMEANSQ =”,f10.5)’) RMEANSQ
WRITE(*,’(”DELTAR =”,f10.5)’) DELTAR
WRITE(*,’(” DEFF =”,f10.5)’) DEFF
WRITE(*,’(” DEFFX =”,f10.5)’) DEFFX
WRITE(*,’(” PECLET =”,f10.5)’) PECLET
WRITE(*,’(” STOKESEFF =”,f10.5)’) STOKESEFF
WRITE(*,’(” ENERGCONVEFF =”,f10.5)’) ENERGCONVEFF
WRITE(*,’(” KINETICEFF =”,f10.5)’) KINETICEFF
WRITE(*,’(” SUZUKIEFF =”,f10.5)’) SUZUKIEFF
WRITE(*,’(” CUBEROEFF =”,f10.5)’) CUBEROEFF
WRITE(4,’(/” ** brownian dynamics ends ** ”///)’)
c ** calculate and write out running averages **
nstep = nstep - 7e5
avv = acv / real ( nstep )
acvsq = ( acvsq / real ( nstep ) ) - avv ** 2
avr = acr / real ( nstep )
avrsq = acrsq / real (nstep )