In order to model this Aeolian transport of sand, JenNins and Pasini ... flux and concentration within the transport layer are very thin and close to the bed.
5HYLVWDGHOD)DFXOWDGGH,QJHQLHUtD8&99RO1SS
NUMERICAL SIMULATION OF AEOLIAN SALTATION WITHIN THE SEDIMENT TRANSPORT LAYER USING GRANULAR KINETIC THEORY JUAN PEDRO MARVAL LUGO1, LUIS ROJAS-SOLORZANO2 AND JENNIFER CURTIS3 8QLYHUVLGDG1DFLRQDO([SHULPHQWDO)UDQFLVFRGH0LUDQGD'HSDUWPHQWRI(QHUJHWLFMSPDUYDO#JPDLOFRP 2 8QLYHUVLGDG6LPRQ%ROLYDU'HSDUWPHQWRI&RQYHUVLRQRI(QHUJ\UURMDV#XVEYH 3 8QLYHUVLGDGGH)ORULGD'HSDUWPHQWRI&KHPLFDO(QJLQHHULQJMFXUWLV#FKHXÀHGX
1
Recibido: marzo 2008
5HFLELGRHQIRUPD¿QDOUHYLVDGR RFWXEUH ABSTRACT
0RVWJUDQXODUÀRZVDWHQYLURQPHQWDOFRQGLWLRQVDUHXQVWHDG\DQGH[KLELWDFRPSOH[SK\VLFDOEHKDYLRU'XQHIRUPDWLRQ DQG PLJUDWLRQ LQ WKH GHVHUW DUH FRQWUROOHG QRW RQO\ E\ WKH ÀRZ RI VDOWDWLQJ SDUWLFOHV RYHU WKH VDQG EHG EXW DOVR E\ WXUEXOHQWDWPRVSKHULFDLUÀRZ,QIDFWVHGLPHQWVDUHWUDQVSRUWHGE\WKHDWPRVSKHULFDLUÀRZZLWKLQDWKLQOD\HURQO\DIHZ FHQWLPHWHUVDERYHWKHVDQG\VXUIDFH7KHVHMXPSLQJSDUWLFOHVUHDFKDPD[LPXPVHGLPHQWPDVVÀX[OHYHODWDFHUWDLQGHOD\ WLPH NQRZQDVWKH³VDWXUDWLRQWLPH´ DIWHUWKHLQLWLDOPRYHPHQWE\VOLGLQJDQGUROOLQJEHJLQV8QOLNHVHGLPHQWWUDQVSRUW LQZDWHUZKHUHWKHSDUWLFOHVDUHOLIWHGE\WKHWXUEXOHQWVXVSHQVLRQWKHVDOWDWLQJSDUWLFOHVDUHNHSWDOLYHLQWKHOD\HUPDLQO\ GXHWRSDUWLFOHSDUWLFOHDQGSDUWLFOHEHGFROOLVLRQV,QRUGHUWRPRGHOWKLV$HROLDQWUDQVSRUWRIVDQG-HQNLQVDQG3DVLQL SURSRVHGDWZRÀXLGPRGHO RQHGLPHQVLRQDODQGVWHDG\VWDWH XVLQJ*UDQXODU.LQHWLF7KHRU\ *.7 WRGHVFULEH WKHVROLGSKDVHVWUHVV7KHSUHVHQWZRUNH[WHQGVWKHRULJLQDOLGHDRI-HQNLQVDQG3DVLQLE\XVLQJDPRUHUREXVWPRGHORI *.7IRUWKHNLQHWLFFROOLVLRQDOFRQWULEXWLRQVWRWKHVROLGSKDVHVWUHVVWHQVRUWRJHWKHUZLWKDIULFWLRQPRGHOIRUVXVWDLQHG FRQWDFWVEHWZHHQSDUWLFOHV,QDGGLWLRQDVWDQGDUGNİWXUEXOHQFHPRGHOIRUWKHDLUDQGDGUDJPRGHOIRUWKHLQWHUDFWLRQ EHWZHHQWKHSKDVHVLVHPSOR\HG$UHFWDQJXODU'JHRPHWU\ZDVFKRVHQZLWKDW\SLFDOORJDULWKPLFSUR¿OHIRUWKHLQOHWDLU YHORFLW\DORQJZLWKDQLQLWLDODPRXQWRIVDQGDWUHVWLQWKHORZHUSDUWRIWKHVLPXODWLRQGRPDLQUHVHPEOLQJWKHSDUWLFOH VDOWDWLQJÀRZFRPPRQO\VHHQLQWKHYHUWLFDOPLGGOHSODQHZLWKLQ6DOWDWLRQZLQGWXQQHOV7KLVPRGHOLVYDOLGDWHGZLWK H[SHULPHQWDOGDWDIURP/LXDQG'RQJ DQGWKHUHVXOWVJLYHQE\3DVLQLDQG-HQNLQV DQG0DUYDOHWDO $JRRGHVWLPDWLRQIRUWKHSDUWLFOHHURVLRQDQGPDVVÀX[LQWKH6DOWDWLRQOD\HULVSUHGLFWHGHYHQWKRXJKWKHSUR¿OHVRIPDVV ÀX[DQGFRQFHQWUDWLRQZLWKLQWKHWUDQVSRUWOD\HUDUHYHU\WKLQDQGFORVHWRWKHEHG Keywords 6DOWDWLRQOD\HU*UDQXODUNLQHWLFWKHRU\&)'
SIMULACION NUMERICA DE LA SALTACION AEREA DENTRO DE LA CAPA DE TRANSPOPRTE DE SEDIMENTOS USANDO TEORIA CINETOCA GRANULAR RESUMEN /D PD\RUtD GH ORV ÀXMRV JUDQXODUHV TXH RFXUUHQ HQ FRQGLFLRQHV DPELHQWDOHV VRQ LQHVWDEOHV \ H[KLEHQ XQ FRPSOHMR FRPSRUWDPLHQWR ItVLFR /D IRUPDFLyQ \ PLJUDFLyQ GH GXQDV HQ ORV GHVLHUWRV VRQ FRQWURODGDV QR VyOR SRU HO ÀXMR GH SDUWtFXODVHQVDOWDFLyQVREUHHOOHFKRDUHQRVRVLQRWDPELpQSRUHOÀXMRWXUEXOHQWRGHDLUHDWPRVIpULFR(VWDVSDUWtFXODVHQ VDOWDFLyQDOFDQ]DQXQÀXMRPiVLFRGHVHGLPHQWDFLyQPi[LPROXHJRGHFLHUWRWLHPSROXHJRGHODPRYLOL]DFLyQLQLFLDOSRU GHVSOD]DPLHQWR\URGDPLHQWR$GLIHUHQFLDGHOWUDQVSRUWHGHVHGLPHQWRVHQDJXDHQHOFXDOODVSDUWtFXODVVRQHOHYDGDVSRU ODVXVSHQVLyQWXUEXOHQWDODVSDUWtFXODVHQVDOWDFLyQVRQDYLYDGDVHQODFDSDSULQFLSDOPHQWHSRUORVFKRTXHVH[LVWHQWHVGH SDUWtFXODSDUWtFXOD\SDUWtFXODOHFKRDUHQRVR3DUDHOPRGHODGRGHODWHRUtDGH7UDQVSRUWH(yOLFRGH$UHQD-HQNLQV\3DVLQL SURSXVLHURQXQPRGHORELIiVLFR XQHVWDGRHVWDEOH\XQLGLPHQVLRQDO XVDQGROD7HRUtD&LQpWLFD*UDQXODU 7&* SDUDGHVFULELUORVHVIXHU]RVHQODIDVHVyOLGD(OSUHVHQWHWUDEDMRH[WLHQGHODLGHDRULJLQDOGH-HQNLQV\3DVLQLXVDQGRXQ PRGHORPiVUREXVWRGH7&*SDUDODVFRQWULEXFLRQHVFLQpWLFDVFROLVLRQDOHVDORVHVIXHU]RVWHQVRUHVGHODIDVHVyOLGDMXQWR a un modelo de fricción para los contactos continuos entre partículas. Fue seleccionada una geometría 2D rectangular con XQSHU¿OORJDUtWPLFRWtSLFRSDUDODYHORFLGDGGHHQWUDGDGHODLUH(VWHPRGHORHVYDOLGDGRFRQGDWRVH[SHULPHQWDOHVGH/LX \'RQJ \ORVUHVXOWDGRVGDGRVSRU3DQVLQL\-HQNLQV \0DUYDO 8QDEXHQDHVWLPDFLyQSDUDODHURVLyQ GHODVSDUWtFXODV\HOÀXMRPiVLFRGHVDOWDFLyQHVSUHGLFKD Palabras Clave 6DOWDWLRQOD\HU*UDQXODUNLQHWLFWKHRU\&)' 80
INTRODUCTION &RPSXWDWLRQDO )OXLG '\QDPLFV &)' KDV UHFHQWO\ LPSDFWHG PDQ\ ¿HOGV RI UHVHDUFK HJ DWPRVSKHULF ÀRZ PLFURÀXLGLFV ZKHUHEHIRUHLWVDSSOLFDWLRQZDVXQWKLQNDEOH PDLQO\GXHWRWKHSUREOHPVDVVRFLDWHGWRWKHZLGHUDQJHLQ OHQJWKDQGWLPHVFDOHVDVZHOODVWKHJUHDWFRPSXWDWLRQDO HIIRUW QHFHVVDU\ WR REWDLQ UHOLDEOH UHVXOWV 6DOWDWLQJ ÀRZV DUH D W\SH RI PXOWLSKDVH ÀRZ ZKHUH &)' FDQ QRZ KDYH DJUHDWLPSDFW7KHÀRZRIVHGLPHQWLQWKHVDOWDWLQJOD\HU LVDPL[WXUHRIWZRGLIIHUHQWSKDVHV$LULVWKHFRQWLQXXP SKDVH RU SULPDU\ SKDVH DQG VDQG LV WKH GLVSHUVHG SKDVH :KHQ WKH VDQG SDUWLFOHV DUH FRQVLGHUHG DV DQ LPPLVFLEOH ³ÀXLG´ LW LV SRVVLEOH WR XVH WKHRULHV GHYHORSHG IRU WZR SKDVHÀRZ,QWKLVFDVHWKHGLVSHUVHGSKDVHRISDUWLFOHVDQG WKHÀXLGSKDVHDUHLQWHUSHQHWUDWLQJFRQWLQXD5HFHQWO\WKH VDOWDWLRQOD\HUKDVEHHQPRGHOHGZLWKPL[WXUHPRGHOV -LHW DO $OKDMUDI QHJOHFWLQJWKHFROOLVLRQDOHIIHFWV intrinsic to the motion of the particles and the coupling of WKH SDUWLFOH LQWHUDFWLRQV ZLWK WKH DLU VWUHDP +RZHYHU D complete description of sediment transport should consider WKUHH GLIIHUHQW ÀRZ UHJLRQV SUHVHQW LQ WKH SKHQRPHQRQ DQG WKHLU LQWHUDFWLRQV 7KH UHJLRQV DUH DV IROORZV D ERWWRPUHJLRQZKHUHWKHSDUWLFOHVDUHDWUHVW VDQGEHG DUHJLRQDERYHWKHVDQGEHGZKHUHWKHSDUWLFOHVDUHMXPSLQJ VDOWDWLQJOD\HU DQG DWRSUHJLRQZKHUHWKHUHLVRQO\DLU 7KHUHIRUHWKHPL[WXUHPRGHOLVLQDGHTXDWHIRUPRGHOLQJ this phenomenon because it does not permit a description of WKHSUHYDLOLQJVWUDWL¿FDWLRQ0XOWL(XOHULDQPRGHOVHQDEOH WKHVLPXODWLRQRIWKHWKUHHFRH[LVWLQJUHJLRQVDQGFRQVLGHU WKHPRYHPHQWRIWKHFORXGRISDUWLFOHVZLWKHTXDWLRQVWKDW DUHVHSDUDWHIURPWKHJDVHTXDWLRQV7KLVDSSURDFKLQYROYHV to incorporate a description for the solid-phase stress tensor XVLQJWKH*UDQXODU.LQHWLF7KHRU\ *.7 The movement of sand dunes is often a problem in desert DUHDV7KHVHGXQHVPRYHRQWRFRPPXQLFDWLRQURXWHVOLNH KLJKZD\VIUHHZD\VRURWKHUSODFHVRILQWHUHVW HJFLWLHV ODNHV DQG SDUNV 7KH GXQH ¿HOGV LQ VRXWKHUQ 0RURFFR 6DXHUPDQQ DQG 'XQHV RI &RUR LQ 9HQH]XHOD DUH H[DPSOHVRIWKHVHSUREOHPV,QWKHVHFDVHVWKHVDQGPXVWEH VKRYHOHGGDLO\IURPWKHXSZLQGVLGHWRWKHGRZQZLQGVLGH RIWKHURDGXVLQJODUJHEXOOGR]HUVREVWUXFWLQJWKHÀRZRI WUDI¿F7KHVLWXDWLRQDVVRFLDWHGZLWKWKH9HQH]XHODGXQHVLV PRUHFRPSOH[EHFDXVHWKHKLJKZD\UXQVRQWRDWKLQSRUWLRQ RI ODQG FDOOHG DQ LVWKPXV ZKLFK PDNHV LW LPSRVVLEOH WR FRQVWUXFW DQRWKHU IUHHZD\ ZLWKRXW SDVVLQJ WKURXJK WKH GXQHV 7KXV RQO\ D EHWWHU XQGHUVWDQGLQJ RI $HROLDQ transport on a micro- and meso-scale, and generating LPSURYHG PDWKHPDWLF PRGHOV ZLOO DLG LQ SUHGLFWLQJ ZLWK UHDVRQDEOH SUHFLVLRQ WKH ÀRZ EHKDYLRU RQ PDFURVFDOH These models enable the design and development of
control devices for dunes motion or in the construction of KLJKZD\VRYHUSDVVHVWXQQHOV 7KH SUHVHQW ZRUN LQWURGXFHV D QXPHULFDO PRGHO IRU SUHGLFWLQJWKHÀRZRIVDQGSDUWLFOHVLQDVDOWDWLQJVWDWHRYHU a sand bed using an Eulerian model for both phases. The gas SKDVHWXUEXOHQFHFORVXUHPRGHOLVWKHVWDQGDUGNİPRGHO DQG WKH VROLGSKDVH VWUHVV LV GHVFULEHG XVLQJ *.7 ZLWK D IULFWLRQDO PRGHO ZKLFK LV DSSOLHG IRU VXVWDLQHG SDUWLFOH contacts. The mathematical model is implemented in the commercial code ANSYS-Fluent® and computed results DUHFRPSDUHGZLWKH[SHULPHQWDOGDWDJLYHQE\/LX 'RQJ DQGWKHPRGHOSUHGLFWLRQVJLYHQE\3DVLQL -HQNLQV DQG0DUYDOHWDO OVERVIEW OF SALTATION :KHQ VDQG SDUWLFOHV DUH GUDJJHG E\ D ORZ YHORFLW\ ZLQG WZR GLIIHUHQW PHFKDQLVPV RI WUDQVSRUW DUH UHODWHG WR WKH VXVSHQVLRQ 7KHVH WZR PHFKDQLVPV DUH NQRZQ DV ³6DOWDWLRQ´ DQG ³5HSWDWLRQ´ $QGUHRWWL HW DO UHFRJQL]HGE\%DJQROG PRUHWKDQ\HDUVDJRLQ KLVFODVVLFDOERRN³7KHSK\VLFVRIEORZQVDQGDQGGHVHUW GXQHV´6DOWDWLRQDQG5HSWDWLRQWRJHWKHUFRQWULEXWHWRWKH SDUWLFOH PDVV ÀX[ RYHU WKH VDQG\ VXUIDFH DQG FRQWURO WKH G\QDPLFVRIDGXQHJLYHQLWVSDUWLFXODUVKDSHFKDUDFWHULVWLFV LQWKHGHVHUW%DVHGRQSUHYLRXVUHVHDUFKLWLVNQRZQWKDW WKHSDUWLFOHVLQWHUDFWZLWKHDFKRWKHUDQGZLWKWKHVDQGEHG $QGUHRWWLHWDO )URPWKHHDUO\ZRUNRI%DJQROG¶V XQWLO QRZ PXFK UHVHDUFK HIIRUW ERWK H[SHULPHQWDO DQG WKHRUHWLFDOQXPHULFDOPRGHOLQJKDVEHHQGHYRWHGWRWU\LQJ WRGHVFULEHWKHFRPSOH[SKHQRPHQRORJ\RIVDOWDWLRQ 0DQ\ H[SHULPHQWDO ZRUNV KDV JLYHQ LQVLJKW LQWR GHWDLOV DERXWWKH6DOWDWLRQOD\HU /LX 'RQJ 1DOSDQLVHW DO 6SLHVHWDO %XWWHU¿HOG 6WRXW ,YHUVHQ 5DVPXVVHQ 5DVPXVVHQ 6RUHQVHQ 2NROL 1LVKLPXUD +XQW 1L HW DO 'RQJHWDO %DXHUHWDO 5HVHDUFK UHVXOWV DUH SUHVHQWHG LQ WHUPV RI WZR OHQJWK VFDOHV 7KH small length scale is on the order of the particle diameter and consists of measurements of the ejection and impact velocities, ejection and impact angles, trajectories of the GLIIHUHQWW\SHVRISDUWLFOHV VQRZPXVWDUGVHHGTXDUW] DQG SDUWLFOHVL]H7KHPHDQHMHFWLRQYHORFLW\LVSURSRUWLRQDOWR WKHIULFWLRQYHORFLW\X ZKLFKFKDUDFWHUL]HVWKHÀRZDQG KDVGLPHQVLRQVRIVSHHGEXWLVDPHDVXUHRIVKHDU(TQ DQGGRHVYHU\PXFKZLWKWKHJUDLQVL]H
u* =
x
t
81
2Q DYHUDJH WKH VDQG SDUWLFOHV DUH OLIWHG LQ WKH DLUÀRZ DW DQ DQJOH EHWZHHQ Û DQG Û DQG GHVFHQG ZLWK DQ LPSDFW DQJOH EHWZHHQ Û DQG Û ,Q JHQHUDO WKHVH DQJOHV GHFUHDVH VOLJKWO\ ZKHQ WKH SDUWLFOHV DUH ODUJHU DQG LQFUHDVHVLJQL¿FDQWO\ZKHQWKHIULFWLRQYHORFLW\LVVOLJKWO\ UDLVHG 2Q WKH ODUJH VFDOH WKH IROORZLQJ PHDVXUHPHQWV DUH UHSRUWHG SDUWLFOH PDVV ÀX[ SDUWLFOH YHORFLW\ DQG FRQFHQWUDWLRQ SUR¿OHV VDOWDWLRQ OD\HU KHLJKW OHQJWK RI VDWXUDWLRQ DQG RWKHUV PDFURVFRSLF ÀRZ IHDWXUHV 7KH PDMRULW\ RI UHVHDUFKHUV DJUHH WKDW WKH SDUWLFOH PDVV ÀX[ and concentration of particles in saturation decrease H[SRQHQWLDOO\ZLWKWKHKHLJKWDERYHWKHVDQGEHGXQOLNHWKH SDUWLFOHYHORFLW\SUR¿OHZKLFKLQFUHDVHVH[SRQHQWLDOO\ZLWK WKHKHLJKW :KLWH )HZVWXGLHVKDYHIRFXVHGRQWKH VDWXUDWLRQOHQJWK KRZHYHUWKHJHQHUDODUJXPHQWLVWKDWWKH HTXLOLEULXPVWDWHLVUHDFKHGLQDVKRUWWLPHIRUDGHWHUPLQDWH OHQJWK 2QO\ UHFHQWO\ 'RQJ HW DO FRQFOXGHG WKDW D FRPSOHWHO\ VDWXUDWHG VWDWH ZLWKLQ WKH 6DOWDWLRQ OD\HU LV DFKLHYHG DIWHU D IHZ PHWHUV GHSHQGLQJ RQ WKH IULFWLRQ YHORFLW\
The main purpose of the researchers has been to obtain UHOLDEOH H[SUHVVLRQV WKDW DOORZ WR FDOFXODWH WKH WRWDO PDVV ÀX[ XQGHU VSHFL¿F FRQGLWLRQV %DJQROG ZDV WKH ¿UVW WR UHSRUW DQ H[SUHVVLRQ IRU WKH VDWXUDWHG PDVV ÀX[ RYHU WKH VDQG EHG ZKLFK VKRZV GHSHQGHQF\ RQ WKH WKLUG SRZHURIWKHIULFWLRQYHORFLW\7KLVHTXDWLRQGRHVQRWWDNH LQWR DFFRXQW WKDW WKH PDVV ÀX[ RQ WKH 6DOWDWLRQ OD\HU LV RQO\ SUHVHQW ZKHQ WKH IULFWLRQ YHORFLW\ LV JUHDWHU WKDQ D WKUHVKROG YHORFLW\ IRU WKH LQLWLDWLRQ RI SDUWLFOH PRYHPHQW RU ZKHQ WKH IULFWLRQ YHORFLW\ LV JUHDWHU WKDQ WKH LPSDFW YHORFLW\IRUJUDLQVDOUHDG\LQDVDOWDWLQJVWDWH([SHULPHQWDO UHVHDUFKKDVLPSURYHGRXUXQGHUVWDQGLQJRIVDOWDWLQJÀRZ EXW WKHVH VWXGLHV KDYH PDLQO\ HQDEOHG WKH GHYHORSPHQW RI PDWKHPDWLFDO H[SUHVVLRQV DQDO\WLFDO DQG H[SHULPHQWDO FRUUHODWLRQV IRU SUHGLFWLQJ WKH PDVV ÀX[ DW VDWXUDWLRQ FRQGLWLRQV 6RUHQVHQ KDV UHSRUWHG DQ LQWHUHVWLQJ VXPPDU\ RI HTXDWLRQV IRU SUHGLFWLQJ WKH VDWXUDWLRQ PDVV ÀX[WKLVFRPSLODWLRQVKRZVPDWKHPDWLFDOH[SUHVVLRQVZLWK DGHSHQGHQF\DFFRUGLQJWRWKHIROORZLQJUHODWLRQVKLS
Qg = f (V) tu 3*
Figure 1. Saltation layer skew. )LJXUH VKRZV D VXPPDU\ RI WKH PDLQ FKDUDFWHULVWLFV RI WKH6DOWDWLRQOD\HUEDVHGRQH[SHULPHQWDOREVHUYDWLRQV$OO WKHVNHWFKHGSUR¿OHVDUHSORWWHGRQDORJDULWKPLFVFDOHLQ WKH YHUWLFDO GLUHFWLRQ 7KH ZLQG EORZV IURP OHIW WR ULJKW ZLWK WKH ZHOONQRZQ ORJDULWKPLF SUR¿OH 6FKOLFKWLQJ DQGHURGHVWKHVDQGEHGVWDUWLQJDWSRLQW$%HWZHHQ SRLQWV $ DQG % WKH DLU YHORFLW\ SUR¿OH LV PRGL¿HG DV momentum is transferred from the air to jumping particles XQWLO HTXLOLEULXP LV UHDFKHG DIWHU SRLQW % %DVLFDOO\ WKH VDWXUDWLRQ OHQJWK GHSHQGV RQ WKH IULFWLRQ YHORFLW\ DQG particle diameter. After the saturation length, the saltating SDUWLFOH ÀRZ UHDFKHV D IXOO\ GHYHORSHG VWDWH DQG DOO ÀRZ YDULDEOHV RQO\ GHSHQG RQ WKH \GLUHFWLRQ 5HFHQW H[SHULPHQWDO ZRUN 'RQJ HW DO KDV FRQ¿UPHG WKDW LQVLGHWKH6DOWDWLRQOD\HUWKHDLUSUR¿OHLVPRGL¿HGE\WKH SUHVHQFHRIWKHSDUWLFOHVDQGWKDWWKHDLUSUR¿OHKDVDFRQYH[ VKDSHQHDUWRWKHVDQGEHG \ įPRP 7RZDUGVWKHWRSRI WKH6DOWDWLRQOD\HU \!įPRP įPRPLVWKHOD\HUZKHUHWKH DSSUHFLDEOH PRPHQWXP LQWHUFKDQJH RFFXUV EHWZHHQ ERWK SKDVHV WKHDLUSUR¿OHLVZHDNO\LQÀXHQFHGE\WKHJUDLQV,Q IDFWWKHSDUWLFOHFRQFHQWUDWLRQDQGPDVVÀX[QHDUWRWKHWRS RIWKH6DOWDWLRQOD\HUWHQGWR]HUR 82
:KHUH 9 LV WKH UDWLR EHWZHHQ WKH IULFWLRQ YHORFLW\ DQG WKUHVKROGYHORFLW\DQGI 9 LVGHWHUPLQHGH[SHULPHQWDOO\ )RU%DJQROG¶VHTXDWLRQWKLVIXQFWLRQVLVDFRQVWDQWZKLFK JLYHV DQ XQWUXH PDVV ÀX[ ZKHQ WKH IULFWLRQ YHORFLW\ LV VPDOOHU WKDQ WKH SLFNXS YHORFLW\ ZKLFK LV WKH YHORFLW\ UHTXLUHGWRVWDUWWKHPRYHPHQWRIWKHSDUWLFOHV(DFKRIWKRVH H[SUHVVLRQVIRUFDOFXODWLQJWKHPDVVÀX[VKDUHVLQFRPPRQ WKDW WKH VDOWDWLQJ SDUWLFOHV WUDYHO LQ WKH ZLQG GLUHFWLRQ WRWDOO\IUHHZLWKRXWREVWDFOHVRUZDOOVWKDWEORFNVLWDORQJ RI LWV WUDMHFWRU\7KLV IHDWXUH ZDV GHVFULEHG SUHYLRXVO\ WR LQGLFDWHWKDWWKHVHNLQGVRIPDFURH[SUHVVLRQVDUHLQYDOLG IRUSUHGLFWLQJDÀX[RIVDOWDWLQJSDUWLFOHVLQWHUDFWLQJZLWK DNLQGRIVWDWLFVROLGVXUIDFH,QDQ\FDVHWKHWRWDOPDVV ÀX[RYHUDYHUWLFDODUHDFDQEHFDOFXODWHGE\PHDQRIWKH IROORZLQJH[SUHVVLRQ
d sal
Q=
# 0
d sal
q sal (y) dy =
#
u s (y) z s (y) dy
0
,I ZH ZDQW WR DSSO\ WKH (TQ >@ LW LV QHFHVVDU\ WR NQRZ D SULRUL WKH VDQG YHORFLW\ DQG VDQG FRQFHQWUDWLRQ SUR¿OHV SUHYLRXVO\ PRGL¿HG E\ WKH SUHVHQFH RI REVWDFOHV DQG WKH IRUPRIVDQGEHGHURGHGDQGWKHUHIRUHRQO\DPXOWLSKDVH numerical simulation could be able to give a result close to WKHUHDOLW\
NOMENCLATURE ǻ[ ǻ\ ds a H Fr J go l / /[ /\ Q 3 p ps R 5H 4 T X 8 8J Y 9U 9 [ [ \
[PHVKVSDFLQJP \PHVKVSDFLQJP Particle diameter, m Particle radius, m 3DUWLFOHFRHI¿FLHQWRIUHVWLWXWLRQ Empirical material constant, Pa *UDYLWDWLRQDOFRQVWDQWPus-2 Radial distribution function Length for the averaging process, m 6KRUWHVWVLJQL¿FDQWPDFURVFRSLFOHQJWKP 'RPDLQKRUL]RQWDOGLPHQVLRQP 'RPDLQYHUWLFDOGLPHQVLRQP (PSLULFDOFRQVWDQWLQIULFWLRQDOSUHVVXUHHTXDWLRQ (PSLULFDOFRQVWDQWLQIULFWLRQDOSUHVVXUHHTXDWLRQ Pressure, Pa Solid pressure or granular pressure, Pa Characteristic length-scale, m 5H\QROGVQXPEHU 7RWDOPDVVÀX[NJāP-3 Time, s )ULFWLRQYHORFLW\ $LUIUHHVWUHDPYHORFLW\Pus-1 $LUYHORFLW\SUR¿OHPus-1 9HORFLW\YHFWRUPus-1 5DWLRRIWHUPLQDOYHORFLW\RIDJURXSRISDUWLFOHVWR that of an isolated particle 5DWLR EHWZHHQ WKH IULFWLRQ YHORFLW\ DQG WKUHVKROG YHORFLW\ 3RVLWLRQYHFWRUP +RUL]RQWDOFRRUGLQDWH 9HUWLFDOFRRUGLQDWH
*UHHNV\PERO Į ȕ ij ȖV țV ȜV ȜPIS ȝ ȝW ȡ įVDO į%/ įPRP į\
9ROXPHIUDFWLRQ ,QWHUSKDVHGUDJFRHI¿FLHQWNJum3us-1 $QJOHRILQWHUQDOIULFWLRQ 'LVVLSDWLRQRIJUDQXODUHQHUJ\NJum-3us-1 &RQGXFWLYLW\RIJUDQXODUHQHUJ\NJum-1us-1 6ROLGEXONYLVFRVLW\3Dus 0HDQIUHHSDWKP 9LVFRVLW\3Dus 7XUEXOHQWYLVFRVLW\3Dus 'HQVLW\NJum3 6DOWDWLRQOD\HUWKLFNQHVVP %RXQGDU\OD\HUWKLFNQHVVP 0RPHQWXPLQWHUFKDQJHOD\HUP ,QLWLDOWKLFNQHVVRIVDQGDWUHVWP
Stress tensor Num-2 ș 6KLHOGSDUDPHWHU Ĭ *UDQXODUWHPSHUDWXUHP2us-2 Subscripts fri J V p NLQ col sal
Frictional *DVSKDVH DLU 6ROLGSKDVH VDQG Particle .LQHWLF Collisional Saltation
AEOLIAN SALTATION CLOUD MODELS 7KH$HROLDQ SURFHVV EHJLQV ZKHQ WKH DHURG\QDPLF IRUFH LHGUDJ GLVORGJHVDIHZJUDLQVZKLFKUROODQGVOLGHRYHU WKHVDQGEHGXQWLOWKH\HQWHULQWR6DOWDWLRQ7KHHQWUDLQHG JUDLQV DUH DFFHOHUDWHG E\ WKH ZLQG DORQJ WKHLU WUDMHFWRU\ PDLQO\ E\ WKH GUDJ IRUFH EHIRUH WKH\ RQFH DJDLQ LPSDFW WKH EHG 7KH SDUWLFOHV ZKLFK LPSDFW WKH EHG SUREDEO\ rebound or dislodge other particles, creating a grain chain reaction. The visualization of particle collisions in H[SHULPHQWV 1DOSDQLV HW DO 0F(ZDQ :LOOHWV ZDVGHULYHGVWDWLVWLFDOO\WRJHWWKH³VSODVKSURFHVV´ WKH LQWHUDFWLRQ EHWZHHQ DQ LPSDFWLQJ JUDLQ DQG WKH EHG DQGWKLVLVPRGHOHGE\WKHVSODVKIXQFWLRQ 8QJDU +DII 7KH WUDYHOLQJ JUDLQV H[WUDFW PRPHQWXP IURP WKH DLU DQG FRQVHTXHQWO\ WKH DLU GHFHOHUDWHV DQG WKH SURFHVV RIGLUHFWDHURG\QDPLFHQWUDLQPHQWGHFUHDVHV7KLVQHJDWLYH IHHGEDFNPHFKDQLVPUHDFKHVHTXLOLEULXPLQWKHVDWXUDWLRQ WLPH WKH WLPH DW ZKLFK WKH DLU ZLOO KROG WKH PD[LPXP DPRXQWRIVDQGSDUWLFOHVLQVDOWDWLRQ DQGWKHQWKHJUDLQV DUHNHSWDOLYHLQWKH6DOWDWLRQVWDWHRQO\E\WKHLUFROOLVLRQV ZLWK WKH VDQG EHG7KH UDWH RI SDUWLFOH WUDQVSRUW EHFRPHV FRQVWDQW ZKHQ HTXLOLEULXP LV UHDFKHG 7KH VWDWLVWLFDO PRGHOV DUH EDVHG RQ WKH IRXU VXESURFHVVHV SUHYLRXVO\ GHVFULEHG DHURG\QDPLFVHQWUDLQPHQWWKHJUDLQWUDMHFWRU\ WKH JUDLQFROOLVLRQ DQG WKH ZLQG YHORFLW\ PRGL¿FDWLRQ $FRPPRQIHDWXUHRIDOOWKHPRGHOVLVWKDWWKH\FRQVLGHU VDOWDWLQJÀRZLQDIUHHFRQGLWLRQRIPRYHPHQWZLWKRXWDQ\ UHVWULFWLRQVRUREVWDFOHV7KLVSUHPLVHPDNHVWKHVWDWLVWLFDO model useless for the prediction of the patterns of erosion DQGRUGHSRVLWLRQDURXQGVROLGREMHFWVFDXVHGE\WKHJUDLQ ZDOOLQWHUDFWLRQ CFD MODELS AND GRANULAR KINETIC THEORY &XUUHQWO\D&)'PHWKRGRORJ\IRUWKHDSSURSULDWHWUHDWPHQW RI WKH VHGLPHQW WUDQVSRUW E\ $HROLDQ 6DOWDWLRQ GRHV QRW H[LVW5HFHQWO\-LHWDO DQG$OKDUDI SURSRVHG 83
interesting models to simulate the large-scale Aeolian VHGLPHQWWUDQVSRUWZLWKRXWFRQVLGHULQJLQGHWDLOWKHÀXLG SDUWLFOH LQWHUDFWLRQ DQG ZLWKRXW WDNLQJ LQWR DFFRXQW WKH FROOLVLRQDO LQWHUDFWLRQ EHWZHHQ SDUWLFOHV DV D IXQGDPHQWDO SDUDPHWHULQWKHLUFRQVLGHUDWLRQV,QERWKZRUNVWKH\XVHG DKRPRJHQRXVPRGHOZLWKDVLQJOHVHWRIHTXDWLRQVDQGWKH GLVSHUVHSKDVHZDVPRGHOHGZLWKDVOLSYHORFLW\EDVHGRQD GUDJIXQFWLRQ$OWKRXJKWKHHIIHFWRISDUWLFOHVRQWKHZLQG SUR¿OHLVZHDNWKHDLUYHORFLW\SUR¿OHZLWKLQWKH6DOWDWLRQ OD\HULVPRGL¿HGE\WKHPRPHQWXPWUDQVIHUEHWZHHQWKHWZR SKDVHV7KHUHIRUHLQRUGHUWRFRQVLGHUWKHPL[WXUHPRGHO LWZRXOGLQYROYHNQRZLQJDSULRULWKHFRQFHQWUDWLRQSUR¿OH YHORFLW\ SUR¿OH DQG RWKHU FKDUDFWHULVWLFV RI WKH PL[WXUH WKDWDUHFRPPRQO\XQDYDLODEOH7KHUHIRUHDOWKRXJKWKHVH SUHYLRXVZRUNVSUHVHQWVLJQL¿FDQWDGYDQWDJHVLQWKHXVHRI &)'WRVLPXODWHWKH6DOWDWLRQOD\HUWKH\GRQRWWDNHLQWR DFFRXQWVRPHIXQGDPHQWDODVSHFWVSUHVHQWLQWKHÀRZ-L HW DO XWLOL]HG WKH ORJDULWKPLF DQG WRWDO HQWUDLQHG SDUWLFOH ÀX[ SUR¿OHV REVHUYHG LQ H[SHULPHQWV LQ RUGHU WR LQFOXGHIXQFWLRQDOUHODWLRQVKLSVLQWRKLVPRGHO0HDQZKLOH $OKDMUDI IRFXVHG KLV ZRUN RQ DQ DOJRULWKP IRU WUHDWLQJDPRYLQJERXQGDU\ HURVLRQGHSRVLWLRQ EDVHGRQ KDQGOLQJDQDUWL¿FLDOFRPSXWDWLRQDOGRPDLQ 2Q WKH RWKHU KDQG WKHUH LV VRPH XQFHUWDLQW\ DERXW WKH SDUWLFOH VXVSHQVLRQ PHFKDQLVPV SUHVHQW LQ VKHHW ÀRZV ZDWHU DQGVDOWDWLQJÀRZV DLU ZKLOHH[SHULPHQWDOUHVXOWV RQO\RIIHUVRPHTXDOLWDWLYHFKDUDFWHULVWLFVRIWKHWZRW\SHV RIÀRZV-HQNLQV +DQHV LPSOHPHQWHGIRUWKH¿UVW WLPHDVLPSOHRQHGLPHQVLRQDODQGVWHDG\WZRÀXLGPRGHO XVLQJ*UDQXODU.LQHWLF7KHRU\IRUGHVFULELQJWKHLQWHUDFWLRQ DPRQJSDUWLFOHVLQWKHVHGLPHQWWUDQVSRUWOD\HUERWKLQZDWHU DQGDLU7KHLUUHVXOWVVKRZHGDTXDOLWDWLYHJRRGHVWLPDWLRQ RI WKH VROLGV HURVLRQ ZKLFK VXJJHVWV WKDW WKH *UDQXODU .LQHWLF7KHRU\FRXOGEHXVHGWRVLPXODWHWKLVSKHQRPHQRQ ZLWKVXFFHVV3DVLQL -HQNLQV H[WHQGHGWKHSUHYLRXV ZRUN E\ DGGLQJ DQ H[WUD WHUP LQ WKH JDV PRPHQWXP HTXDWLRQV RULJLQDWLQJ IURP D VHFRQG DYHUDJLQJ SURFHVV JLYHQ LQ +VX HW DO ZKLFK GHVFULEHG DQ DGGLWLRQDO mechanism of suspension due to turbulence effects together ZLWK WKH JUDQXODU SUHVVXUH JUDGLHQW )XUWKHUPRUH 3DVLQL -HQNLQV FRQVLGHUHG D ÀRZ UHJLPH EHWZHHQ WKH Saltation regime and turbulent suspension regime called WKHFROOLVLRQDOUHJLPHZKLOHDWVDPHWLPHWKH\UHFRJQL]HG WKDWWKLVUHJLPHKDGQRW\HWEHHQREVHUYHGH[SHULPHQWDOO\ 7KHLUUHVXOWVJDYHDQRYHUHVWLPDWLRQRIWKHVDQGPDVVÀX[ LQ FRPSDULVRQ ZLWK H[SHULPHQWDO UHVXOWV IURP SDUWLFOHV LQ the Saltation state.
ÀXFWXDWLRQYHORFLW\ 6DXHUPDQQ 3\H 7VRDU GHPRQVWUDWLQJ WKDW LQ WKH 6DOWDWLRQ OD\HU WKH WXUEXOHQW VXVSHQVLRQ LV QHJOLJLEOH ,QVWHDG WKH SDUWLFOHV DUH RQO\ NHSW LQ PRYHPHQW LQ WKH 6DOWDWLRQ OD\HU E\ EHGSDUWLFOH collisions and momentum transfer from the air stream to the SDUWLFOHV,QDQ\FDVHWKH6DOWDWLRQOD\HUUHSUHVHQWVDGLOXWH UHJLPH ZKHUH WKH PD[LPXP YROXPH IUDFWLRQ LV DURXQG /LX 'RQJ LQGLFDWLQJ WKDW WKH VDOWDWLQJ SDUWLFOHVFRQWULEXWHRQO\WRWKHNLQHWLFFRQWULEXWLRQWRWKH VROLGVYLVFRVLW\DQGQRWWKHFROOLVLRQDOUHJLPHOLNH-HQNLQV +DQHV DQG3DVLQL -HQNLQV VXJJHVW MATHEMATICAL MODEL $FFRUGLQJ WR 'RQJ WKUHH ZHOONQRZQ UHJLRQV DUH considered in our simulations, according to the solid volume IUDFWLRQDQGYHORFLW\¿HOGV2QHUHJLRQFRUUHVSRQGVWRWKH VDQG EHG ZKHUH WKH VROLGV SDFNLQJ OHYHO LV KLJK DQG WKH JDVDQGVROLGYHORFLWLHVDUHHTXDOWR]HUR7KLV]RQHFDQQRW PRYH OLNH D ÀXLG LW LV D FRPSOHWHO\ VROLG UHJLRQ ZKLFK LV VXSHU¿FLDOO\ HURGHG E\ WKH ZLQG DFWLRQ 7KH VHFRQG UHJLRQ LV IHZ FHQWLPHWHUV DERYH WKH VDQG EHG 6DOWDWLRQ OD\HU ZKHUHSDUWLFOHVWUDYHOTXLFNO\GULYHQE\DWXUEXOHQW VKHDULQJDLUÀRZ,QWKHWKLUGUHJLRQRXWVLGHRIWKH6DOWDWLRQ OD\HU WKHUH LV RQO\ DLU PRYLQJ SDUDOOHO WR WKH EHG ,Q WKLV ZRUN ZH SURSRVHG DQ (XOHULDQ PRGHO LQKRPRJHQHRXV PRGHO IRUERWKSKDVHVEDVHGRQWKHIDFWWKDWWKHUHDUHWKUHH ]RQHVZLWKGLIIHUHQWFRQFHQWUDWLRQDQGYHORFLWLHV VLPLODUWR DVWUDWL¿HGÀRZ )XUWKHUPRUHWKH6WRNHVQXPEHUIRUWKH 6DOWDWLRQOD\HUFRQGLWLRQLVHVWLPDWHGWREHJUHDWHUWKDQ DFFRUGLQJWRWKHXVHU¶VPDQXDORI)/8(17$16
