Study of Polarized Wave with a Hydrodynamic Model and Fourier Spectral Method

<table>Time evolution for the hydrodynamics equation together with initial wave <svg height="19.924999" id="M73" style="vertical-align:-3.276pt" version="1.1" viewbox="0 0 324.48749 19.924999" width="324.48749" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,15.775)"><path d="M515 96q-37 -44 -84.5 -76t-71.5 -32q-12 0 -18 7t-6 34t10 75l19 89h-2q-87 -113 -181 -176q-43 -29 -83 -29q-46 0 -46 63q0 31 9 67l52 222q7 36 -1 36q-10 0 -76 -50l-13 24q47 45 92 71.5t67 26.5q18 0 21 -16.5t-8 -65.5l-56 -242q-16 -67 13 -67q41 0 114 74
t114 146l32 145l74 26h11q-52 -199 -80 -347q-8 -39 6 -39q7 0 32 17.5t47 39.5z" id="x1D462"></path></g><g transform="matrix(.017,-0,0,-.017,13.916,15.775)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,28.62,15.775)"><path d="M298 36l-68 -35q-25 -13 -42 -13q-98 0 -98 119v290h-69l-4 11l28 29h45v72l65 67l14 -2v-137h113q8 -8 6 -21t-12 -19h-107v-256q0 -96 59 -96q29 0 61 17z" id="x74"></path></g><g transform="matrix(.017,-0,0,-.017,33.805,15.775)"><path d="M433 39l-85 -51q-27 0 -47 20q-16 16 -24 44q-5 -3 -24 -16.5t-27 -18.5t-23 -13.5t-25 -12t-19 -3.5q-52 0 -86 36.5t-34 85.5q0 71 81 99q128 43 155 65v17q0 54 -23 83.5t-62 29.5q-28 0 -45 -19t-29 -66q-7 -23 -29 -23q-15 0 -29 13t-14 30t31 40q25 19 68.5 40.5
t81.5 29.5q49 0 82 -27q45 -39 45 -123v-185q0 -60 35 -60q18 0 36 11zM275 84v156q-12 -6 -49 -23t-41 -19q-31 -14 -46 -31t-15 -43q0 -35 22.5 -55.5t48.5 -20.5q19 0 42 10.5t38 25.5z" id="x61"></path></g><g transform="matrix(.017,-0,0,-.017,41.267,15.775)"><path d="M524 0h-209v26q44 5 55 18t11 65v167q0 114 -91 114q-50 0 -109 -51v-235q0 -50 10 -61.5t55 -16.5v-26h-217v26q51 5 62 17t11 61v206q0 47 -9.5 58.5t-50.5 19.5v23q82 15 139 40v-79l67 47q47 30 83 30q60 0 94.5 -40t34.5 -112v-193q0 -50 10 -61.5t54 -16.5v-26z
" id="x6E"></path></g><g transform="matrix(.017,-0,0,-.017,50.566,15.775)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,56.447,15.775)"><path d="M475 507q0 -83 -20 -172t-56 -167.5t-93.5 -129t-125.5 -50.5q-157 0 -157 227q0 78 21.5 164t59 161t96.5 123.5t126 48.5q79 0 114 -58t35 -147zM391 522q0 155 -81 155q-62 0 -111 -82.5t-73 -200.5h253q12 81 12 128zM373 346h-255q-12 -91 -12 -150q0 -72 20 -123
t63 -51q34 0 64 28.5t52.5 77t39 103t28.5 115.5z" id="x1D703"></path></g><g transform="matrix(.017,-0,0,-.017,64.913,15.775)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g><g transform="matrix(.017,-0,0,-.017,70.794,15.775)"><path d="M319 325l-25 -7q-33 99 -103 99q-29 0 -47 -19.5t-18 -49.5t22 -49.5t62 -36.5q63 -26 95 -57t32 -79q0 -64 -50 -101t-115 -37q-35 0 -67.5 10.5t-46.5 23.5q-5 11 -11 51t-6 67l27 5q14 -53 46.5 -88.5t75.5 -35.5q29 0 50 19.5t21 50.5t-19.5 51.5t-59.5 39.5
q-28 12 -46 22.5t-38.5 27t-30.5 38.5t-10 49q0 54 42.5 92t109.5 38q48 0 88 -18q6 -15 13 -50.5t9 -55.5z" id="x73"></path></g><g transform="matrix(.017,-0,0,-.017,77.169,15.775)"><path d="M380 106l15 -22q-42 -52 -87 -74t-76 -22q-89 0 -142 62t-53 151q0 103 66 175q67 73 146 73h1q61 0 105.5 -46.5t44.5 -105.5q0 -16 -4.5 -22.5t-18.5 -9.5q-54 -10 -262 -17q1 -95 46 -143.5t104 -48.5q61 0 115 50zM225 412q-35 0 -64.5 -32.5t-40.5 -90.5q83 0 170 3
q14 0 19 4.5t5 15.5q1 41 -23.5 70.5t-65.5 29.5z" id="x65"></path></g><g transform="matrix(.017,-0,0,-.017,84.478,15.775)"><path d="M390 111l17 -21q-34 -45 -80 -73.5t-89 -28.5q-91 0 -146 62t-55 147q0 118 101 195q74 57 149 57h1q59 0 90 -27q16 -14 16 -30q0 -15 -12 -29t-21 -14q-8 0 -19 11q-44 41 -101 41q-52 0 -87.5 -42.5t-35.5 -117.5q0 -49 15 -87t39 -58t49 -30t48 -10q33 0 60.5 12
t60.5 43z" id="x63"></path></g><g transform="matrix(.017,-0,0,-.017,91.584,15.775)"><path d="M513 0h-211v26q44 5 55 17t11 60v166q0 123 -98 123q-59 0 -104 -51v-238q0 -48 10 -60.5t55 -16.5v-26h-213l1 26q46 4 57 16t11 61v477q0 48 -9.5 58.5t-58.5 14.5v24q63 9 147 35v-336q73 73 143 73q64 0 101 -42t37 -129v-175q0 -49 11 -61t55 -16v-26z" id="x68"></path></g>
<g transform="matrix(.012,-0,0,-.012,100.65,7.613)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
<g transform="matrix(.017,-0,0,-.017,106.988,15.775)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,112.869,15.775)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.017,-0,0,-.017,126.145,15.775)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,139.897,15.775)"><path d="M383 397q0 -32 -35 -49q-12 -6 -23 8q-37 45 -84 45t-90 -71q-40 -65 -40 -167q0 -57 22 -86t59 -29q38 0 81.5 24.5t69.5 51.5l16 -21q-44 -53 -104 -84t-109 -31q-56 0 -89.5 41t-33.5 117q0 61 30 124t79 105q33 28 81 50.5t86 22.5q34 0 59 -15.5t25 -35.5z" id="x1D450"></path></g><g transform="matrix(.017,-0,0,-.017,146.799,15.775)"><path d="M324 430l-26 -36l-112 -4l-55 -265q-13 -66 7 -66q13 0 44.5 20t50.5 40l17 -24q-38 -40 -85.5 -73.5t-87.5 -33.5q-50 0 -21 138l55 262h-80l-2 8l25 34h66l25 99l78 63l10 -9l-37 -153h128z" id="x1D461"></path></g><g transform="matrix(.017,-0,0,-.017,152.595,15.775)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g><g transform="matrix(.017,-0,0,-.017,162.25,15.775)"><path d="M535 230h-212v-233h-58v233h-213v50h213v210h58v-210h212v-50z" id="x2B"></path></g><g transform="matrix(.017,-0,0,-.017,176.002,15.775)"><path d="M298 36l-68 -35q-25 -13 -42 -13q-98 0 -98 119v290h-69l-4 11l28 29h45v72l65 67l14 -2v-137h113q8 -8 6 -21t-12 -19h-107v-256q0 -96 59 -96q29 0 61 17z" id="x74"></path></g><g transform="matrix(.017,-0,0,-.017,181.187,15.775)"><path d="M433 39l-85 -51q-27 0 -47 20q-16 16 -24 44q-5 -3 -24 -16.5t-27 -18.5t-23 -13.5t-25 -12t-19 -3.5q-52 0 -86 36.5t-34 85.5q0 71 81 99q128 43 155 65v17q0 54 -23 83.5t-62 29.5q-28 0 -45 -19t-29 -66q-7 -23 -29 -23q-15 0 -29 13t-14 30t31 40q25 19 68.5 40.5
t81.5 29.5q49 0 82 -27q45 -39 45 -123v-185q0 -60 35 -60q18 0 36 11zM275 84v156q-12 -6 -49 -23t-41 -19q-31 -14 -46 -31t-15 -43q0 -35 22.5 -55.5t48.5 -20.5q19 0 42 10.5t38 25.5z" id="x61"></path></g><g transform="matrix(.017,-0,0,-.017,188.649,15.775)"><path d="M524 0h-209v26q44 5 55 18t11 65v167q0 114 -91 114q-50 0 -109 -51v-235q0 -50 10 -61.5t55 -16.5v-26h-217v26q51 5 62 17t11 61v206q0 47 -9.5 58.5t-50.5 19.5v23q82 15 139 40v-79l67 47q47 30 83 30q60 0 94.5 -40t34.5 -112v-193q0 -50 10 -61.5t54 -16.5v-26z
" id="x6E"></path></g><g transform="matrix(.017,-0,0,-.017,197.948,15.775)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,203.829,15.775)"><path d="M475 507q0 -83 -20 -172t-56 -167.5t-93.5 -129t-125.5 -50.5q-157 0 -157 227q0 78 21.5 164t59 161t96.5 123.5t126 48.5q79 0 114 -58t35 -147zM391 522q0 155 -81 155q-62 0 -111 -82.5t-73 -200.5h253q12 81 12 128zM373 346h-255q-12 -91 -12 -150q0 -72 20 -123
t63 -51q34 0 64 28.5t52.5 77t39 103t28.5 115.5z" id="x1D703"></path></g><g transform="matrix(.017,-0,0,-.017,212.295,15.775)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g><g transform="matrix(.017,-0,0,-.017,221.015,15.775)"><path d="M380 106l15 -22q-42 -52 -87 -74t-76 -22q-89 0 -142 62t-53 151q0 103 66 175q67 73 146 73h1q61 0 105.5 -46.5t44.5 -105.5q0 -16 -4.5 -22.5t-18.5 -9.5q-54 -10 -262 -17q1 -95 46 -143.5t104 -48.5q61 0 115 50zM225 412q-35 0 -64.5 -32.5t-40.5 -90.5q83 0 170 3
q14 0 19 4.5t5 15.5q1 41 -23.5 70.5t-65.5 29.5z" id="x65"></path></g><g transform="matrix(.017,-0,0,-.017,228.24,15.775)"><path d="M474 0h-214v26q32 4 39 12t-5 28l-74 117q-39 -58 -64 -106q-16 -25 -10 -35t44 -16v-26h-186v26q46 6 64 17t48 49t85 120l-87 136q-23 36 -40 47t-58 16v26h199v-26q-31 -4 -36.5 -12.5t9.5 -29.5l59 -93q37 56 58 94q13 23 7 30t-40 11v26h183v-26q-47 -5 -65 -13.5
t-45 -39.5q-24 -29 -78 -110l97 -148q28 -43 46.5 -55.5t63.5 -18.5v-26z" id="x78"></path></g><g transform="matrix(.017,-0,0,-.017,236.263,15.775)"><path d="M169 380l92 53q26 16 49 16q80 0 128.5 -56t48.5 -133q0 -108 -72.5 -180.5t-176.5 -91.5q-38 1 -69 18v-131q0 -57 13 -70t71 -18v-27h-236v27q49 5 61 17t12 61v452q0 48 -9.5 58.5t-56.5 16.5v24q66 11 145 43v-79zM169 346v-270q37 -39 94 -39q63 0 101.5 50
t38.5 134q0 78 -38.5 123t-95.5 45q-45 0 -100 -43z" id="x70"></path></g><g transform="matrix(.017,-0,0,-.017,248.009,15.775)"><path d="M290 -163h-170v866h170v-28q-79 -7 -94 -19.5t-15 -72.5v-627q0 -59 14.5 -71.5t94.5 -19.5v-28z" id="x5B"></path></g><g transform="matrix(.017,-0,0,-.017,253.874,15.775)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g><g transform="matrix(.017,-0,0,-.017,258.514,15.775)"><path d="M480 416q0 -21 -18 -41q-9 -11 -17 -7q-20 9 -42 9q-62 0 -140 -78q23 -69 88 -192q17 -31 27 -42t20 -11q16 0 62 46l17 -20q-64 -92 -119 -92q-35 0 -70 66q-41 73 -84 187q-36 -30 -62 -61q-27 -115 -35 -172q-41 -8 -78 -20l-6 6l140 612q7 28 0.5 34t-37.5 7l-34 1
l5 26q38 4 74 13.5t57 17t25 7.5q12 0 4 -32l-104 -443h2q35 38 97 93q39 35 65.5 56t62 41.5t58.5 20.5q19 0 30.5 -10t11.5 -22z" id="x1D458"></path></g><g transform="matrix(.017,-0,0,-.017,267.065,15.775)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,272.946,15.775)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.017,-0,0,-.017,286.222,15.775)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,299.974,15.775)"><path d="M383 397q0 -32 -35 -49q-12 -6 -23 8q-37 45 -84 45t-90 -71q-40 -65 -40 -167q0 -57 22 -86t59 -29q38 0 81.5 24.5t69.5 51.5l16 -21q-44 -53 -104 -84t-109 -31q-56 0 -89.5 41t-33.5 117q0 61 30 124t79 105q33 28 81 50.5t86 22.5q34 0 59 -15.5t25 -35.5z" id="x1D450"></path></g><g transform="matrix(.017,-0,0,-.017,306.876,15.775)"><path d="M324 430l-26 -36l-112 -4l-55 -265q-13 -66 7 -66q13 0 44.5 20t50.5 40l17 -24q-38 -40 -85.5 -73.5t-87.5 -33.5q-50 0 -21 138l55 262h-80l-2 8l25 34h66l25 99l78 63l10 -9l-37 -153h128z" id="x1D461"></path></g><g transform="matrix(.017,-0,0,-.017,312.672,15.775)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g><g transform="matrix(.017,-0,0,-.017,318.554,15.775)"><path d="M226 -163h-170v27q79 7 94 20t15 73v627q0 59 -15 72t-94 20v27h170v-866z" id="x5D"></path></g>
</svg>, <svg height="12.4375" id="M74" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 48.212502 12.4375" width="48.212502" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path d="M475 507q0 -83 -20 -172t-56 -167.5t-93.5 -129t-125.5 -50.5q-157 0 -157 227q0 78 21.5 164t59 161t96.5 123.5t126 48.5q79 0 114 -58t35 -147zM391 522q0 155 -81 155q-62 0 -111 -82.5t-73 -200.5h253q12 81 12 128zM373 346h-255q-12 -91 -12 -150q0 -72 20 -123
t63 -51q34 0 64 28.5t52.5 77t39 103t28.5 115.5z" id="x1D703"></path></g><g transform="matrix(.017,-0,0,-.017,13.237,12.162)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,27.94,12.162)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g><g transform="matrix(.017,-0,0,-.017,36.1,12.162)"><path d="M113 -12q-24 0 -39.5 16t-15.5 42q0 24 16 40.5t40 16.5t40 -16.5t16 -40.5q0 -26 -16 -42t-41 -16z" id="x2E"></path></g><g transform="matrix(.017,-0,0,-.017,39.976,12.162)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g>
</svg> and <svg height="11.125" id="M75" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 36.75 11.125" width="36.75" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,10.862)"><path d="M509 373q0 -52 -42.5 -118.5t-124.5 -159.5l-94 -107h-27q-22 209 -78 342q-21 48 -41 48q-18 0 -63 -62l-22 21q77 111 125 111q30 0 52 -21.5t38 -75.5q37 -130 52 -285q48 62 74.5 99.5t53.5 87t27 76.5q0 43 -17 60q-5 5 -5 15q0 17 12.5 30.5t28.5 13.5
q21 0 36 -22.5t15 -52.5z" id="x1D708"></path></g><g transform="matrix(.017,-0,0,-.017,13.814,10.862)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,28.518,10.862)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g>
</svg>.</table>

Modelling and Simulation in Engineering

fig3

Figure 3

Figure 3: Study of Polarized Wave with a Hydrodynamic Model and Fourier Spectral Method 