Effect of Turbulence on the Drag and Lift of a Particle P. Bagchi and S

Total Page:16

File Type:pdf, Size:1020Kb

Effect of Turbulence on the Drag and Lift of a Particle P. Bagchi and S ¢¡¤£¦¥¨§ © §£¦¥¦£© §£! " $#% &§ © ' ()"*§+#¥,-£ .0/21436587%9;:36<>=?@/A1436BC367%9A36<;=A36D EGF$HA36DJILKMFN<OI,P6QR¦9;FNP8DSF$IL:T7U36B3V<>=XWYH>H>BT:CF$=ZFN7%9A36<;:C7N[ \]<;:T^8F$D_[S:TIa`bP6QdcaBCBT:C<>P6:C[43eI,\]D_f>36<A3egihj9A36KMHA36:C56<lk \]DSfA36<A3;k;cimon;prq8stp8kA\G?uW v!F$`xw¨P8D_=;[NyjZz>BTI_:CH>9>36[_FG{APrw|[$k2I_z>D_f;z>BCF$<>7NFVkAHA36DSI_:C7$z>B}3VI_FY{APrw|[Nk2=>D_365M36<>=BC:~QIUk HA36DJIL:C7$BCF%ILz>DSf>z>BTFN<>7$F!:C<OILF$DL367-IL:CP6< ] ;u2 W=>:CDSFN7$I,<z>KMFNDS:C7U3VB[_:CKz>B}3VI_:CP8<E!G?>¦:C[|z>[SFN=I_P[JILz>=t`I_9>F"F$@FN7$I]P6Q3MQD_F$FN[JILD_FN36K :T[_P6I_D_P8H>:T7]I_z>D_f>z;BCFN<OI4{>P&wP8<xI_9>FG=>D_36536<>=bBT:TQIjQP6D_7NF$[4P8<X3[SH>9>FNDS:C7N36BlHA36DJIL:T7NBCFV/R¦9>F!HA3VDJg I_:C7NBTFG=>:C36KF-ILF$D4:C[|36fP8z;IGp8/I_PpUsI_:CKMFN[jI_9>F"v!P8BCKMP858P8DSPr^[S7U36BTF6k;I_9>FHA36DJIL:C7$BCF!|F$`u<>P6BC=>[ <z>KfFNDY:C[!36fP8z;IGn8sILPxn8s8s;k*36<>=IL9;FQD_FNF$[SI_D_FU3VK ILz;D_f>z>BTFN<>7$F:C<OILF$<>[_:~Ia`:T[G36fP8z;IprsILP 6¢b/Xc£I:C["[_9>Prw|<¤I_9A3VII_9>FQDSFNF$[SILDSFU36K¥I_z>D_f>z;BCFN<;7NFx=>PuFN[ <>P6I 9A3U^8Fb3[_z>f;[SI%3V<I_:}36B¨36<>= ¡ [J`t[JILF$K3VI_:C7F$@FN7-IP8<)IL9;FMI_:CKMF-gi3U^6FNDL3V58FN=KFN36<=>DL365;/R¦9>FM[SIL36<>=A36DS=¤=>D_3657NP8DSD_FNBC3VIL:TP8< fA3V[_FN=¤P8<¤IL9;F:C<;[SI%3V<IL36<>F$P8z>[P6D KMFU36<¤[_BC:TH^6FNBTPt7$:TIa`D_F$[_z>B~IL[:T<¤3D_FU3V[_P8<A3Vf>BT`367$7Nz>D_3VILF H>DSFN=>:T7$I_:CP8<bP6QI_9>F!KFN36<b=>D_365P8f;I%3V:C<>F$=xQD_P8K¦IL9>FE!Y?@/>§,P&w¨F$^6FND$kI_9>F367N7$z>DL3V7$`bP6Q0H>DSF-g p =>:T7$I_:CP8<oPVQ,I_9>Fx:T<>[SIL36<OI%36<>F$P8z>[=>D_365=>F$7ND_FN36[_F$[Mw|:~IL9¨:T<>7NDSFU36[S:C<>5H>36DSI_:C7NBTF[S:C©NFV/¤ªAP8D IL9>F [SK36BTBCF$D,H>36DSI_:C7NBTFN[$kI_9>FBTP&w«QDSFN¬z>F$<>7$`P8[_7$:CBCBC3VIL:TP8<>[¦:T<IL9>FE!G?=;DL365b3VD_F w¨FNBTBd7N36H;I_z>D_F$= fO`I_9>F[SI%3V<>=A36DS=¤=>D_365>kf>z;IQP8DGI_9>FB}36DS58FNDYHA36DSI_:C7$BCFN[![S:C58<;:T­A7N36<OI=>:T@FNDSFN<>7$FN["F-®t:C[JIF$^8F$< QP8D|IL9;FBCPrw¯QDSFN¬z>F$<>7$`7NP8KMH@P6<>FN<OIL[$/|ca<>7NBTz>[_:TP8<P6QI_9>F3V=>=>FN=ugaK36[_[]36<>=9;:C[SI_P8DS`bQP6D_7NF$[Nk 7$P8KH;z;ILF$=°fA36[SFN= P6<°IL9>FX{Az>:T=±^6FNBCPu7N:~Ia`²3VIMIL9>F7NFN<OI_FNDPVQGIL9>FHA36DJIL:C7$BCFVk¨=>PuFN[<>P6I:CKg H>DSPr^8FI_9>FH>D_F$=>:C7-IL:TP8<l/xEG:~*F$D_F$<I"F$[SIL:TK3eILFN[P6Q¨IL9>F{Az>:T=^6FNBTPt7$:TIa`[SFNFN<¤f`I_9>FMHA36DSI_:C7$BCF 36DSFxF-®;36KM:C<;FN=l/)c£I:T[[_9>Prw|<¨IL9>3VIIL9;FbKFN36<o=>D_365:C[:T<>[_F$<>[_:~IL:~^8FI_PIL9>F{>z>:C=o^6FNBCPu7N:~Ia` KMFU36[Sz>D_F$=b3VI¨IL9>FGH>36DSI_:C7NBTF]7$FN<OILF$DNk>P6D¨P8ftI%36:T<>FN=xfO`3U^8FND_3658:T<>5Pr^8F$D¦3 {Az>:T=^6P8BCz>KMFYPVQIL9>F P8DS=>FNDP6Q,IL9;FH>36DSI_:C7NBTFx[_:T©NFV/²R¦9;F{Az;7$ILz>3VIL:TP8<>[=>:TK:T<>:C[S9¨36[IL9>FX[S:C©NFXP6Q,IL9;F3U^6FNDL3V58:C<>5 ^6P8BCz>KMF:T<>7NDSFU36[SFN[$/R¦9>FbF$@FN7-IP6Q,:C<>7$D_FN36[_:T<>5QDSFNFN[JILDSFU36K³ILz>DSf>z>BCF$<>7NFb:T<OILFN<;[_:TIa`QP8D"IL9>F [_36KF HA36DJIL:C7$BCF"[_:T©NF:T[Y[SI_z>=>:TFN=l/ªBTz>7$I_zA3VI_:CP8<>[]:C<IL9>F=;DL365X36<;=BC:~QI]QP8DS7NF$[!36D_F [_9>Prw|<ILP [S7U36BTFw|:~IL9IL9>F"KFN36<=>DL36536<>=XQD_F$FN[SI_D_FN36K :C<OILF$<>[_:~Ia`8/4R¦9>F [SIL36<>=A36DS==>D_365w|:TIL9;P8z;I,IL9>F 36=;=>FN=tg£K3V[_[,36<>=X9>:T[SILP6DS`xQP8DS7NF$[,H;D_Pr^u:C=>F$[|I_9>Ff@F$[SI]36H>H;D_PU®t:CK3VIL:TP8<xILPIL9;F E!Y?xD_FN[Sz>BTIN/ ¡ ´jµ¶´·x¸¹ºX»¼½¦¸¾NºM· EG:T[_HFND_[SFN=¨HA3VDSIL:T7Nz>BC3VILFM{APrw|[36fP8z><>=¨:T<o<A3VI_z>D_Fb3V<>=¨:C<¨F$<>58:C<;FNFNDS:C<>536H>H>BT:C7U3eIL:CP6<>[N/ca< KMP8[SI7U36[SFN[$kI_9>Fj{Az;:C=H>9>36[_F¦:T[dILz;D_f>z>BTFN<OIU/R¦9>F4:C<OI_FNDL3V7$IL:TP8<fF$Iaw¨FNF$<MI_9>Fj{Az>:T=H>9A3V[_F|36<>= I_9>FHA3VDSIL:T7Nz>BC3VILFH;9A36[_F:C[f>:~g£=>:CDSFN7-IL:CP6<A36B¿y!I_9>F7N36D_DS:CFNDigaH>9>36[_FILz>DSf>z>BTFN<>7$F:C<t{Az>FN<;7NFN[ IL9>F =>:T[_HFNDS[_:CP6<36<>=H>D_F-QFNDSFN<OIL:C36BA367$7Nz>Kz;B}3VI_:CP8<PVQI_9>F¦HA36DJIL:T7NBCF$[w|9>:T7%9M:T<ILz>DS<KPu=>z>BC3VILFÀIL9>F {Az;:C=I_z>D_f;z>BCF$<>7NFV/W4IjIL9;FGBCF-^8FNBP6Q3 [_:C<;58BCF|HA36DJIL:T7NBCFVkOIL9>FGF-*F$7$IjP6QlQD_F$FN[JILD_FN36KÁILz>DSf>z>BCF$<>7NF :T[4ILPKPu=>:~Q`I_9>F!=>DL3V5QP6D_7NF7$P8KH>36D_F$=XI_PIL9A3VI,:T<X3[JILFN36=;`Xz><>:~QP8D_K¦{APrwÁP6QILF$<7N36BCBTFN=  [JI%36<;=A36D_==>DL3V58Ã-/¦Ä!<IL9>F P6IL9;FND,9A36<>=lk@3H>36DSI_:C7NBTF!7U36<KPu=>:TQ`QD_F$FN[SI_D_FN36KÅILz>DSf>z>BCF$<>7NF fO`I_9>F4QP8D_K3VI_:CP8<PVQl3Gwj36Æ6FVktHFNDS:CPu=>:C7¨[_9>F$=>=>:C<;5 P6Q^8P6DSIL:T7NF$[Nku36<>=wj36Æ6F|ILz>DSf>z>BTFN<>7$F6/R¦9>F 7$P8BCBTFN7$I_:T^6FGF$@FN7-IYP6QI_9>FH>D_F$[_FN<;7NFPVQ3=>:T[SILDS:Cf>ztIL:CP6<xP6QdHA36DSI_:C7$BCFN[¦7N36<bQz>DSI_9>FND|KMPt=;:TQ`MIL9>F F-*F$7$IL:~^8FM=>DL3V5QP8DS7NFP8<¤3H>36DSI_:C7NBTF=>z>FILP[_7$D_F$FN<>:T<>5F-*F$7$I36<>=)IL9>F$D_FNfO`:C<t{Az>FN<;7NFIL9>F KMFU36<¤[_F-I_ILBT:C<>536<>=¤=>:C[SH@F$D_[S:CP8<¤7%9A36DL3V7$ILF$D_:T[SIL:T7N[$/?u:CKM:CB}3VD_BT`6k*IL9;F7NP6BCBCF$7$I_:T^8FMF$@FN7-IP6QjIL9>F =>:T[_HFNDS[_:CP6<P6QjH>36DSI_:C7NBTFN[Yw|:TBCBd=>F$I_FND_KM:C<;FIL9>FM3VISILF$<uz>3VIL:TP8<P8D36z>56KF$<IL3VIL:TP8<P6QÀIL9;F{Az>:C= I_z>D_f>z;BCFN<;7NF6/ ca<bIL9>:T[jHA3VH@F$D4wjFGw|:TBCB@36=>=>DSFN[_[4I_9>F!F$@FN7-I|P6QIL9>F]QD_FNF$[SI_D_FU3VKÇILz>DSf>z>BTFN<>7$F!P8<xIL9>FY=>DL365 36<;=¨BC:TQIQP8DS7NF$[P8<¨3[_:T<>58BTFHA3VDSIL:T7NBTF6/ȤFx7%9>PuP8[_FxI_P:C<O^6FN[SI_:C5O3eILFI_9>Fx[_:TKH;BCFN[JI7N36[_FbP6Q 36<):T[_P6I_D_P8H>:T7 QDSFNFN[JILDSFU36KÉI_z>D_f;z>BCF$<>7NFP6Q¨R3U`uBCP6DGK:T7NDSP8[_7N36BCF|F$`u<>P6BC=>[!<z>KfFNDF$¬uz>36BILP prnVÊ>/R¦9>F¦Iaw¨P:TKHP8DSIL36<OIHA3VDL36KMF$I_FND_[+P6QIL9;F,H;D_P8f>BTFNKË36D_F¦IL9;FN<IL9>F|DL3eIL:CP!PVQ@I_9>F|HA36DSI_:C7$BCF =>:C36KMF$ILF$DÌILPI_9>FvP6BCKMP858P8DSP&^[_7N36BCFÍ@kd36<>=¤I_9>FM:C<OILF$<>[_:~Ia`¤PVQ4QD_F$FN[SI_D_FN36K¥ILz>DSf>z>BCF$<>7NF =>F-­A<>F$=36[dIL9>F¦D_3VIL:TPGP6Q2IL9>F¦DSPtPVISgaKMFU3V<tga[S¬uz>36D_FD_KM[%0I_z>D_f;z>BCF$<I{Az>7$I_zA3VI_:CP8<I_PI_9>F4KMFU36< DSFNB}3eIL:T^6FM^6FNBTPt7$:TIa`)fF$Iaw¨FNFN<¨I_9>FHA3VDSIL:T7NBTF36<;=¤IL9>F[_z>DSD_P8z;<>=>:C<;5{Az>:T=l/XR¦9>Fx=>:C36KF-ILF$D P6Q I_9>FGHA36DJIL:T7NBCF]z><>=>F$D47NP6<>[_:T=>FND_3VIL:TP8<x:C[¨^63VD_:CF$=QD_P6KÇ36fP8z;IYpGILPbprs"IL:CKMFN[ÀIL9>FGvP6BCKMP858P8DSP&^ Î [S7U36BTF6k,36<>= I_9>FILz;D_f>z>BTFN<>7$F:C<OI_FN<>[S:TIa`±:C[^V36DS:CF$=°QD_P8K prs¤I_P 8¢x/hjP8<;[_FN¬z>F$<OILBT`°IL9>F ¡ HA3VDSIL:T7NBTF,F-`t<;P8BC=>[<z>KfFND$kdf>36[_F$=¤P8<)IL9>FMD_F$B}3VI_:T^6F^8FNBTPt7$:TIa`)3V<>=¤H>36DSI_:C7NBTF=>:}3VKF-ILFND$k:C[ :T<XI_9>FDL36<;58Fn8sILPMn8s8st/ hjP8<>[S:C=>F$DIL9;F7U3V[_FP6Q!3¤HA36DJIL:C7$BCF[_F-I_ILBT:C<>5)IL9;D_P8z>569Ï3¤ILz;D_f>z>BTFN<OI{>P&w/'R¦9>FKMFU36< [SF$I_I_BC:T<>5¦^8FNBTPt7$:TIa`P6QtIL9>FÀHA36DJIL:T7NBCFH;D_Pr^u:C=>F$[3|7NP8<O^8F$<>:CF$<OIKMFU36[Sz>D_FÀP6Q;I_9>FÀKFN36<=>D_365|QP8D_7$F6/ ca<bF®;HFNDS:CKMFN<OIL[¨IL9>FGKMFU3V<x=>DL3V57NPuF$ÐM7N:TFN<OI¦:T[j7NP6KH>ztILFN=bf>36[_F$=bP8<xIL9>FYKFN36[_z>DSFNKMFN<OI¦P6Q I_9>FGKMFU36<b[SF$I_I_BC:T<>5^8F$BCPu7N:~Ia`bÑ;Ò36<>=3QP8D_7$FGfA36BC36<>7$FGf@F-IawjF$FN<XIL9>F!56DL3U^u:TIa`36<>=bIL9>FY=>DL365 QP8DS7NF"36[$y p Ê Ó¦ÔÕ Jp& Ì@×ØÏp& ÑÙ Î@Ö Ò¤Ú w|9>F$D_F|×:T[IL9>F¦H>36DSI_:C7NBTF-g¿ILPVg¿{Az>:C==>FN<>[S:TIa`D_3VIL:TP>kO36<;= :C[I_9>F,367N7$FNBTFNDL3eIL:CP6<=>z>F¦ILP58D_3r^u:~Ia`8/ Ö c£QuIL9>FÀ36Kf>:CF$<I*{APrwo:C[[JI%3656<A36<OIUk6Ñ>Ò7NP6D_D_F$[_HP8<>=>[ILP4I_9>FILF$D_KM:C<A36Be^8F$BCPu7N:TIa`G36<>=IL9>F+3Vf@Pr^8F DSFNB}3eIL:CP6<`u:TFNBC=;[I_9>Fj[JI%36<;=A36D_= =>DL3V5Y7NPuF$ÐM7$:CFN<OId7NP8DSD_F$[_HP8<>=>:T<>5]I_PG3|z><>:TQP6D_K¯<>P6<I_z>D_f;z>BCF$<I {APrw/4ca<3I_z>D_f;z>BCF$<I,{APrwk29>PrwjF-^8F$DNk@I_9>FNDSF36DSF!IawjPMwjF$BCBz><>=>F$D_[JILPuPt=KF$7%9A36<>:T[_KM[|I_9A3VI :T<;{Az>F$<>7NF4I_9>F¦KFN36<[_F-I_ILBT:C<>5YDL3VI_F6/R¦9>F4­>D_[SIÀ:T[=>z>F4I_PGIL9>F|<>P8<ugaBC:T<>FU3VDd=>F$H@F$<>=>F$<>7NF,P6QIL9>F =>D_365P6<IL9>F]D_F$B}3VI_:T^8F|^8F$BCPu7N:~Ia`M3VIj­A<;:TILF]|F$`u<>P8BT=>[j<z>KfFNDS[N/+ªAP8D+IL9>F][L36KMFY=;FN<>[S:TIa`DL3VI_:CP ×3V<>=o=>:C36KF-ILF$D"Ì2k0IL9>FKFN36<¤[SF$I_I_BC:T<>5^6FNBTPt7$:TIa`):C<¤3I_z>D_f>z;BCFN<OI {APrwÁ:C["BCFN[S[ I_9A36<¤I_9A3VI :T<X3[JI%3656<A36<OI¦{APrwÁÛZFN:Ûk>WY=>DS:}36<XÜݧY36<;DL3VISIa`AÞ%/R¦9>F[_F$ISILBT:C<>5^6FNBCPu7N:~Ia`x=>F$7NDSFU36[SFN[|w|:TIL9 :T<>7NDSFU36[S:C<>5GI_z>D_f>z;BCFN<;7NF,:C<OILF$<>[_:~Ia`36<>=IL9;F,DSFN[Sz>BTI_:C<>5KMFU36<=>D_365"36[58:T^6FN<MfO`Jp&:T[9>:C589;FND I_9A36<IL9A3eI+fA36[_F$=P6<I_9>F4ILF$D_KM:C<A3VBt^8F$BCPu7N:~Ia`:C<3![JI%3658<>36<OI{APrw/R¦9>:C[+F-*F$7$I+w|:CBTBt=>F$7NDSFU36[SF w|:~IL9=;FN7NDSFU36[S:C<>5,F-`u<>P8BC=;[j<z>KfFNDj36<>=w|:TBCBAF$<I_:CDSFNBT`^V36<;:C[_9:C<IL9>F]BC:T<>FU36D¨?uI_P8Æ6F$[jBT:CKM:TIU/ R¦9>F![SFN7$P8<>=lkA3V<>=bKP6D_FY7$P8KMH>BCF®*ktKF$7%9A36<>:T[_K¦:C[4=;z>FYILPI_9>FGH>DSF$QF$D_FN<OI_:}36B*I_DL3&ßSF$7$I_P8D_:TFN[4P6Q QDSFNFNB~`)Q36BCBT:C<>5HA36DJIL:C7$BCF$[N/.d36DSI_:C7NBTFN[ =>P<>P6I[L36KMH>BCFMIL9>FILz>DSf>z>BCF$<OI{APrw¦z><>:TQP6D_KMBT`8kdf>z;I Ê 101 Uhlerr & Sinclair (1970) I = 0.4 Torobin & Gauvin (1961) X 0.3 0 0.2 Clamen & Gauvin (1969) D 10 XX C XX X Standard drag curve 10-1 0.2 I = 0.4 0.3 I = 0.1 Clift & Gauvin (1970) 100 101 102 103 104 105 106 Re ª:T58z>DSF¦p8yW²[Sz>KMK36DJ`!P6QtIL9>F+D_F$[_z>B~IL[0P8<IL9>FÀF$@FN7-IP6QtILz>DSf>z>BTFN<>7$F¨P8<"IL9>F+=>D_365,7NPuF$ÐM7N:TFN<OIU/ .+DSFN[SFN<OI,D_FN[Sz>BTI_[Náâ¢ã!P8D_FGÜËhjDSP&w¨F6áäÞ¿å!æç?t36<>Æe3658:TD_:@ÜÝ|z;¦áäÞ¿è]éMêl3VD_:C<ÜÝ,:C7%9>P6BC[Náìë]í à Ȥ36D_<;:C7U3îNï¨ð8ñ~òTáôóGõÅ|z>=>P8B~QÜË14367%9A36BTP>á÷ö¦ø1jDSz>7U3eILPîNïjð8ñ~òC/ å H>DSF$QF$DD_F$58:CP8<;[P6Q|=>Prw|<Ow43V[_9²7NP8KMHA36DSFN=oILPDSFN58:TP8<>[ P6Q|z>Htg£KPr^u:C<>5{Az>:C=¢[_F$FxÈo36<;5Ü Ù Z3&®;F-` /R¦9z>[¨I_9>FYKMFU36<{Az>:T=M^6FNBTPt7$:TIa`[_F$FN<xfO`M3"HA36DSI_:C7$BCF¦=>:~*F$D_[+QD_P8KËI_9>F|ILD_z;F]KMFU36< P8ftI%36:T<>FN=xfO`3U^8FND_3658:T<>5Pr^8F$D4IL9>FGF$<OIL:CDSFY^8P6BCz>KMFYP6Q{>z>:C=l/\]<>BT:CÆ6F,IL9>FYF$@FN7$I4P6Q<>P8<>BT:C<>FN36D =>D_365!=>F$H@F$<>=>FN<;7NF6kuIL9>F|F-*F$7$IÀP6QH>D_F-QFNDSFN<OIL:C36B>ILD_3&ßSFN7-ILP8DJ`:C[I_P:C<>7$D_FU3V[_F4IL9;F,KMFU3V<[_F$ISILBT:C<>5 ^6FNBCPu7N:~Ia`8/4R¦9z>[]IL9>F"=>DL3V57$PtF-Ð7$:CF$<IYF$^V36BTzA3VILF$=fA3V[_FN=P8<oip&4ILFN<;=>[]I_PfF
Recommended publications
  • Airfoil Drag by Wake Survey Using Ldv
    AIRFOIL DRAG BY WAKE SURVEY USING LDV 1 Purpose This experiment introduces the student to the use of Laser Doppler Velocimetry (LDV) as a means of measuring air flow velocities. The section drag of a NACA 0012 airfoil is determined from velocity measurements obtained in the airfoil wake. 2 Apparatus (1) .5 m x .7 m wind tunnel (max velocity 20 m/s) (2) NACA 0015 airfoil (0.2 m chord, 0.7 m span) (3) Betz manometer (4) Pitot tube (5) DISA LDV optics (6) Spectra Physics 124B 15 mW laser (632.8 nm) (7) DISA 55N20 LDV frequency tracker (8) TSI atomizer using 50cs silicone oil (9) XYZ LDV traversing system (10) Computer data reduction program 1 3 Notation A wing area (m2) b initial laser beam radius (m) bo minimum laser beam radius at lens focus (m) c airfoil chord length (m) Cd drag coefficient da elemental area in wake survey plane (m2) d drag force per unit span f focal length of primary lens (m) fo Doppler frequency (Hz) I detected signal amplitude (V) l distance across survey plane (m) r seed particle radius (m) s airfoil span (m) v seed particle velocity (instantaneous flow velocity) (m/s) U wake velocity (m/s) Uo upstream flow velocity (m/s) α Mie scatter size parameter/angle of attack (degs) δy fringe separation (m) θ intersection angle of laser beams λ laser wavelength (632.8 nm for He Ne laser) ρ air density at NTP (1.225 kg/m3) τ Doppler period (s) 4 Theory 4.1 Introduction The profile drag of a two-dimensional airfoil is the sum of the form drag due to boundary layer sep- aration (pressure drag), and the skin friction drag.
    [Show full text]
  • Drag Force Calculation
    DRAG FORCE CALCULATION “Drag is the component of force on a body acting parallel to the direction of relative motion.” [1] This can occur between two differing fluids or between a fluid and a solid. In this lab, the drag force will be explored between a fluid, air, and a solid shape. Drag force is a function of shape geometry, velocity of the moving fluid over a stationary shape, and the fluid properties density and viscosity. It can be calculated using the following equation, ퟏ 푭 = 흆푨푪 푽ퟐ 푫 ퟐ 푫 Equation 1: Drag force equation using total profile where ρ is density determined from Table A.9 or A.10 in your textbook A is the frontal area of the submerged object CD is the drag coefficient determined from Table 1 V is the free-stream velocity measured during the lab Table 1: Known drag coefficients for various shapes Body Status Shape CD Square Rod Sharp Corner 2.2 Circular Rod 0.3 Concave Face 1.2 Semicircular Rod Flat Face 1.7 The drag force of an object can also be calculated by applying the conservation of momentum equation for your stationary object. 휕 퐹⃗ = ∫ 푉⃗⃗ 휌푑∀ + ∫ 푉⃗⃗휌푉⃗⃗ ∙ 푑퐴⃗ 휕푡 퐶푉 퐶푆 Assuming steady flow, the equation reduces to 퐹⃗ = ∫ 푉⃗⃗휌푉⃗⃗ ∙ 푑퐴⃗ 퐶푆 The following frontal view of the duct is shown below. Integrating the velocity profile after the shape will allow calculation of drag force per unit span. Figure 1: Velocity profile after an inserted shape. Combining the previous equation with Figure 1, the following equation is obtained: 푊 퐷푓 = ∫ 휌푈푖(푈∞ − 푈푖)퐿푑푦 0 Simplifying the equation, you get: 20 퐷푓 = 휌퐿 ∑ 푈푖(푈∞ − 푈푖)훥푦 푖=1 Equation 2: Drag force equation using wake profile The pressure measurements can be converted into velocity using the Bernoulli’s equation as follows: 2Δ푃푖 푈푖 = √ 휌퐴푖푟 Be sure to remember that the manometers used are in W.C.
    [Show full text]
  • Chapter 4: Immersed Body Flow [Pp
    MECH 3492 Fluid Mechanics and Applications Univ. of Manitoba Fall Term, 2017 Chapter 4: Immersed Body Flow [pp. 445-459 (8e), or 374-386 (9e)] Dr. Bing-Chen Wang Dept. of Mechanical Engineering Univ. of Manitoba, Winnipeg, MB, R3T 5V6 When a viscous fluid flow passes a solid body (fully-immersed in the fluid), the body experiences a net force, F, which can be decomposed into two components: a drag force F , which is parallel to the flow direction, and • D a lift force F , which is perpendicular to the flow direction. • L The drag coefficient CD and lift coefficient CL are defined as follows: FD FL CD = 1 2 and CL = 1 2 , (112) 2 ρU A 2 ρU Ap respectively. Here, U is the free-stream velocity, A is the “wetted area” (total surface area in contact with fluid), and Ap is the “planform area” (maximum projected area of an object such as a wing). In the remainder of this section, we focus our attention on the drag forces. As discussed previously, there are two types of drag forces acting on a solid body immersed in a viscous flow: friction drag (also called “viscous drag”), due to the wall friction shear stress exerted on the • surface of a solid body; pressure drag (also called “form drag”), due to the difference in the pressure exerted on the front • and rear surfaces of a solid body. The friction drag and pressure drag on a finite immersed body are defined as FD,vis = τwdA and FD, pres = pdA , (113) ZA ZA Streamwise component respectively.
    [Show full text]
  • Drag Coefficients of Inclined Hollow Cylinders
    Drag Coefficients of Inclined Hollow Cylinders: RANS versus LES A Major Qualifying Project Report Submitted to the Faculty of the WORCESTER POLYTECHNIC INSTITUTE in partial fulfillment of the requirements for the Degree of Bachelor of Science In Chemical Engineering By Ben Franzluebbers _______________ Date: April 2013 Approved: ________________________ Dr. Anthony G. Dixon, Advisor Abstract The goal of this project was use LES and RANS (SST k-ω) CFD turbulence models to find the drag coefficient of a hollow cylinder at various inclinations and compare the results. The drag coefficients were evaluated for three angles relative to the flow (0⁰, 45⁰, and 90⁰) and three Reynolds numbers (1000, 5000, and 10000). The drag coefficients determined by LES and RANS agreed for the 0⁰ and 90⁰ inclined hollow cylinder. For the 45⁰ inclined hollow cylinder the RANS model predicted drag coefficients about 0.2 lower the drag coefficients predicted by LES. 2 Executive Summary The design of catalyst particles is an important topic in the chemical industry. Chemical products made using catalytic processes are worth $900 billion a year and about 75% of all chemical and petroleum products by value (U.S. Climate Change Technology Program, 2005). The catalyst particles used for fixed bed reactors are an important part of this. Knowing properties of the catalyst particle such as the drag coefficient is necessary to understand how fluid flow will be affected by the packed bed. The particles are inclined at different angles in the reactor and can come in a wide variety of shapes meaning that understanding the effect of particle shape on drag coefficient has practical importance.
    [Show full text]
  • On the Generation of a Reverse Von Kármán Street for the Controlled Cylinder Wake in the Laminar Regime
    On the generation of a reverse Von Kármán street for the controlled cylinder wake in the laminar regime Michel Bergmann,∗ Laurent Cordier, and Jean-Pierre Brancher LEMTA, UMR 7563 (CNRS - INPL - UHP) ENSEM - 2, avenue de la forêt de Haye BP 160 - 54504 Vandœuvre cedex, France (Dated: December 8, 2005) 1 Abstract In this Brief Communication we are interested in the maximum mean drag reduction that can be achieved under rotary sinusoidal control for the circular cylinder wake in the laminar regime. For a Reynolds number equal to 200, we give numerical evidence that partial control restricted to an upstream part of the cylinder surface may increase considerably the effectiveness of the control. Indeed, a maximum value of relative mean drag reduction equal to 30% is obtained when applying a specific sinusoidal control to the whole cylinder, where up to 75% of reduction can be obtained when the same control law is applied only to a well selected upstream part of the cylinder. This result suggests that a mean flow correction field with negative drag is observable for this controlled flow configuration. The significant thrust force that is locally generated in the near wake corresponds to a reverse Kármán vortex street as commonly observed in fish-like locomotion or flapping wing flight. Finally, the energetic efficiency of the control is quantified by examining the Power Saving Ratio: it is shown that our approach is energetically inefficient. However, it is also demonstrated that for this control scheme the improvement of the effectiveness goes generally with an improvement of the efficiency. Keywords: Partial rotary control ; Cylinder wake ; Drag minimization ; reverse Kármán vortex street.
    [Show full text]
  • Aerodynamics of High-Performance Wing Sails
    Aerodynamics of High-Performance Wing Sails J. otto Scherer^ Some of tfie primary requirements for tiie design of wing sails are discussed. In particular, ttie requirements for maximizing thrust when sailing to windward and tacking downwind are presented. The results of water channel tests on six sail section shapes are also presented. These test results Include the data for the double-slotted flapped wing sail designed by David Hubbard for A. F. Dl Mauro's lYRU "C" class catamaran Patient Lady II. Introduction The propulsion system is probably the single most neglect­ ed area of yacht design. The conventional triangular "soft" sails, while simple, practical, and traditional, are a long way from being aerodynamically desirable. The aerodynamic driving force of the sails is, of course, just as large and just as important as the hydrodynamic resistance of the hull. Yet, designers will go to great lengths to fair hull lines and tank test hull shapes, while simply drawing a triangle on the plans to define the sails. There is no question in my mind that the application of the wealth of available airfoil technology will yield enormous gains in yacht performance when applied to sail design. Re­ cent years have seen the application of some of this technolo­ gy in the form of wing sails on the lYRU "C" class catamar­ ans. In this paper, I will review some of the aerodynamic re­ quirements of yacht sails which have led to the development of the wing sails. For purposes of discussion, we can divide sail require­ ments into three points of sailing: • Upwind and close reaching.
    [Show full text]
  • Chapter 4: Immersed Body Flow [Pp
    MECH 3492 Fluid Mechanics and Applications Univ. of Manitoba Fall Term, 2017 Chapter 4: Immersed Body Flow [pp. 445-459 (8e), or 374-386 (9e)] Dr. Bing-Chen Wang Dept. of Mechanical Engineering Univ. of Manitoba, Winnipeg, MB, R3T 5V6 When a viscous fluid flow passes a solid body (fully-immersed in the fluid), the body experiences a net force, F, which can be decomposed into two components: a drag force F , which is parallel to the flow direction, and • D a lift force F , which is perpendicular to the flow direction. • L The drag coefficient CD and lift coefficient CL are defined as follows: FD FL CD = 1 2 and CL = 1 2 , (112) 2 ρU A 2 ρU Ap respectively. Here, U is the free-stream velocity, A is the “wetted area” (total surface area in contact with fluid), and Ap is the “planform area” (maximum projected area of an object such as a wing). In the remainder of this section, we focus our attention on the drag forces. As discussed previously, there are two types of drag forces acting on a solid body immersed in a viscous flow: friction drag (also called “viscous drag”), due to the wall friction shear stress exerted on the • surface of a solid body; pressure drag (also called “form drag”), due to the difference in the pressure exerted on the front • and rear surfaces of a solid body. The friction drag and pressure drag on a finite immersed body are defined as FD,vis = τwdA and FD, pres = pdA , (113) ZA ZA Streamwise component respectively.
    [Show full text]
  • Upwind Sail Aerodynamics : a RANS Numerical Investigation Validated with Wind Tunnel Pressure Measurements I.M Viola, Patrick Bot, M
    Upwind sail aerodynamics : A RANS numerical investigation validated with wind tunnel pressure measurements I.M Viola, Patrick Bot, M. Riotte To cite this version: I.M Viola, Patrick Bot, M. Riotte. Upwind sail aerodynamics : A RANS numerical investigation validated with wind tunnel pressure measurements. International Journal of Heat and Fluid Flow, Elsevier, 2012, 39, pp.90-101. 10.1016/j.ijheatfluidflow.2012.10.004. hal-01071323 HAL Id: hal-01071323 https://hal.archives-ouvertes.fr/hal-01071323 Submitted on 8 Oct 2014 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. I.M. Viola, P. Bot, M. Riotte Upwind Sail Aerodynamics: a RANS numerical investigation validated with wind tunnel pressure measurements International Journal of Heat and Fluid Flow 39 (2013) 90–101 http://dx.doi.org/10.1016/j.ijheatfluidflow.2012.10.004 Keywords: sail aerodynamics, CFD, RANS, yacht, laminar separation bubble, viscous drag. Abstract The aerodynamics of a sailing yacht with different sail trims are presented, derived from simulations performed using Computational Fluid Dynamics. A Reynolds-averaged Navier- Stokes approach was used to model sixteen sail trims first tested in a wind tunnel, where the pressure distributions on the sails were measured.
    [Show full text]
  • List of Symbols
    List of Symbols a atmosphere speed of sound a exponent in approximate thrust formula ac aerodynamic center a acceleration vector a0 airfoil angle of attack for zero lift A aspect ratio A system matrix A aerodynamic force vector b span b exponent in approximate SFC formula c chord cd airfoil drag coefficient cl airfoil lift coefficient clα airfoil lift curve slope cmac airfoil pitching moment about the aerodynamic center cr root chord ct tip chord c¯ mean aerodynamic chord C specfic fuel consumption Cc corrected specfic fuel consumption CD drag coefficient CDf friction drag coefficient CDi induced drag coefficient CDw wave drag coefficient CD0 zero-lift drag coefficient Cf skin friction coefficient CF compressibility factor CL lift coefficient CLα lift curve slope CLmax maximum lift coefficient Cmac pitching moment about the aerodynamic center CT nondimensional thrust T Cm nondimensional thrust moment CW nondimensional weight d diameter det determinant D drag e Oswald’s efficiency factor E origin of ground axes system E aerodynamic efficiency or lift to drag ratio EO position vector f flap f factor f equivalent parasite area F distance factor FS stick force F force vector F F form factor g acceleration of gravity g acceleration of gravity vector gs acceleration of gravity at sea level g1 function in Mach number for drag divergence g2 function in Mach number for drag divergence H elevator hinge moment G time factor G elevator gearing h altitude above sea level ht altitude of the tropopause hH height of HT ac above wingc ¯ h˙ rate of climb 2 i unit vector iH horizontal
    [Show full text]
  • Airfoils and Wings
    Airfoils and Wings Eugene M. Cliff 1 Introduction The primary purpose of these notes is to supplement the text material re- lated to aerodynamic forces. We are mainly interested in the forces on wings and complete aircraft, including an understanding of drag and related nomeclature. 2 Airfoil Properties 2.1 Equivalent Force Systems In some cases it’s convenient to decompose the forces acting on an airfoil into components along the chord (chordwise) and normal to it. These forces are related to lift and drag through the geometry shown in Figure 1. From l n d α c Figure 1: Force Systems the figure we have c(α)=cn(α)cosα − cc(α)sinα cd(α)=cn(α)sinα + cc(α)cosα Obviously, we can also express the normal and chordwise forces in terms of section lift and drag. 1 Figure 2: Flow Decomposition 2.2 Circulation Theory of Lift A typical flow about a lift-producing airfoil can be decomposed into a sum of two flows, as shown in Figure 2. The first flow (a) is ’symmetric’ flow and so produces no lift. The ’circulatory’ flow (b) is responsible for the net higher speed (and hence lower pressure) on the top of the airfoil (the suction side). This can be quantified by introducing the following line integral −→ −→ ΓC = u · d s C This is the circulation of the flow about the path C. It turns out that as long as C surrounds the airfoil (and doesn’t get too close to it), the the value of Γ is independent of C.
    [Show full text]
  • Calculation of Aerodynamic Drag of Human Being in Various Positions Calculation of Aerodynamic Drag of Human Being in Various Positions Mun Hon Koo*, Abdulkareem Sh
    EURECA 2013 - Calculation of Aerodynamic Drag of Human Being in Various Positions Calculation of Aerodynamic Drag of Human Being in Various Positions Mun Hon Koo*, Abdulkareem Sh. Mahdi Al-Obaidi Department of Mechanical Engineering, school of Engineering, Taylor’s University, Malaysia, *Corresponding author: [email protected] Abstract— This paper studies the aerodynamic drag of human 2. Methods being in different positions through numerical simulation using CFD with different turbulence models. The investigation 2.1. Theoretical Analysis considers 4 positions namely (standing, sitting, supine and According to Hoerner [4], the total aerodynamic drag of human squatting) which affect aerodynamic drag. Standing has the body can be classified into 2 components which are: highest drag value while supine has the lowest value. The = + numerical simulation was carried out using ANSYS FLUENT and compared with published experimental results. Where is pressure drag coefficient and is friction Aerodynamic drag studies can be applied into sports field related drag coefficient. Pressure drag is formed from the distribution of applications like cycling and running where positions optimising forces normal .to the human body surface [5]. The effects of are carried out to reduce drag and hence to perform better viscosity of the moving fluid (air) may contribute to the rising value during the competition. of pressure drag. Drag that is directly due to wall shear stress can be knows as friction drag as it is formed due to the frictional effect. The Keywords— Turbulence Models, Drag Coefficient, Human Being, friction drag is the component of the wall shear force in the direction Different Positions, CFD toward the flow, and it depends on the body surface area and the magnitude of the wall shear stress.
    [Show full text]
  • New Aproximate Equations to Estimate the Drag Coefficient of Different Particles of Regular Shape
    NEW APROXIMATE EQUATIONS TO ESTIMATE THE DRAG COEFFICIENT OF DIFFERENT PARTICLES OF REGULAR SHAPE A. D. SALMAN and A. VERBA Department of Mechanical Engineering for the Chemical Industry. Technical University, H-1521 Budapest Received October 20, 1986 Abstract Approximate equations for drag coefficients of different particle shapes are proposed in the form of Kaskas equation [I]. The constants of these equations have been determined by the least square method from experimental values [2, 3]. The constants of these equations are summarized in Table 1. These equations proved to be valid even at Reynolds numbers between 0,006 and 20000, the maximum relative error being 12~o' Introduction In many cases, to discuss the particle dynamics in a fluid, the relationship between the particle Reynolds number and the drag coefficient is indispensable for the calculation of the particle motion or its trajectories. For spherical particles extensive and adequate data have been collected and it has been found possible to present these data by empirical Equations [1,9, 10]. However, for nonspherical particles, although invaluable data have been published, they are presented either in Tables [2, 3J or plotted in Figures [4,6,7,8]. The sphericity employed by Waddel [4J is the most convenient among the several shape factors in use. This is defined as follows; 'l' = s/S, where s is the surface area of a sphere of the same volume as the particle; S is the actual surface area of the particle. The maximum value obtained by this formula is 1, which is the numerical expression for the degree of true sphericity of a sphere.
    [Show full text]