مدل‌ سازی عددی اندرکنش موج و موج ‌شکن متخلخل مستغرق با روش احجام محدود

به منظور شبیه سازی مسائل مربوط به اندر کنش موج و سازه های متخلخل، یک مدل عددی غیرهیدرواستاتیک در فضای دوبعدی قائم (2dv) توسعه داده شده است. معادلات حاکم بر مدل، معادلات ناویر-استوکس اصلاح‌شده برای محیط متخلخل می‌باشند که در آن مقاومت محیط متخلخل با در نظر گرفتن نیروهای دراگ و اینرسی در معادلات ناویر-استوکس متداول لحاظ شده است. گسسته‌سازی روابط حاکم بر مدل، با روش احجام محدود (fvm) در سیستم اختیاری لاگرانژی-اویلری (ale) صورت گرفته است. آلگوریتم حل عمومی بر اساس روش تفکیک زمانی (پروجکشن) می‌باشد. معادلات حاکم بر مدل در دو گام اصلی حل شده‌اند. ابتدا جمله‌ی فشار از معادلات مقدار حرکت حذف شده و معادلات انتقال، پخشیدگی و نیروی دراگ به منظور محاسبه‌ی سرعت‌های میانی حل می‌شوند. در گام دوم، با اعمال معادله‌ی پیوستگی در معادلات مقدار حرکت با در نظر گرفتن جمله‌ی فشار، معادله‌ی پواسون فشار به دست می‌آید. گام دوم با اصلاح مقادیر سرعت و محاسبه‌ی تراز سطح آزاد آب تکمیل می‌گردد. به‌منظور ارزیابی عملکرد مدل در غیاب محیط متخلخل، نتایج مدل با نتایج تحلیلی موجود برای انتشار موج یکتا در آب با عمق ثابت مقایسه و سپس مدل عددی حاضر برای شبیه‌سازی اندرکنش موج و موج‌شکن متخلخل مستغرق به کار گرفته شده است. مقایسه‌ نتایج عددی با نتایج آزمایشگاهی برای تاریخچه‌ی زمانی تراز سطح آزاد آب، توزیع مکانی سطح آب، میدان‌های سرعت و پروفیل‌های قائم سرعت، نشان‌دهنده‌ توانایی مدل عددی حاضر در شبیه‌سازی اندرکنش موج و سازه‌ متخلخل می‌باشد.

numerical modeling of wave interaction with submerged pearmeable breakwater by finite volume method

this paper presents a non-hydrostatic two-dimensional vertical (2dv) numerical model for the simulation of wave-porous structure problems. the flow in both porous and pure fluid regions is described by the extended navier-stokes equations, in which the resistance to flow through a porous medium is considered by including the additional terms of drag and inertia forces. the finite volume method (fvm) in an arbitrary lagrangian-eulerian (ale) description is employed to discretize the flow and transport equations. a two-step fractional method has been deployed to solve the governing equations. in the first step, the momentum equations in the absence of pressure field were solved to compute an intermediate velocity. the second fractional step consisted of bringing the pressure terms back into the equations, and calculating the pressure field by solving the extended continuity equation and the momentum equations excluding advective and diffusive terms and drag force components. by substitution of the approximations of the pressure derivatives into momentum equations, and subject to the continuity constraint, the pressure poisson equation was obtained. the solution of the pressure poisson equation led to a linear system of equations in the form of a block tri-diagonal matrix with the pressures as unknowns. the second step was completed by computing the updated velocity values. in the present numerical model, two types of boundary conditions, namely dirichlet and neumann boundary conditions were adapted to solve the governing equations. the dirichlet boundary condition was set to zero for normal velocities at impermeable bottom and the neumann boundary condition was considered to be equal to zero for normal gradient of the tangential velocities at impermeable bed and also the left side of the computational domain. at open boundaries, where required, by setting the dynamic pressure equal to zero at the end of the numerical domain, a free exit for water was considered. the newly developed model in the absence of porous medium was verified by comparing the numerical simulations with the analytical solutions of a solitary wave propagation in a constant water depth. the newly developed model was then employed to simulate the solitary wave interaction with a permeable submerged breakwater. based on the numerical results, when the solitary wave front reaches the offshore side of the submerged breakwater, due to the hydraulic jump formation, the flow is separated from the top of the obstacle and small clockwise vortices are generated at the leading edge of the breakwater. as the wave passes over the breakwater, the primary vortex grows in size and penetrates into the deeper layers of water. it was also seen that, due to the drag and inertia resistance forces of the porous medium, the velocity inside the permeable breakwater was noticeably smaller than that on the top of the breakwater. the comparisons between the numerical results and experimental measurements for time histories of water displacements, spatial distributions of free surface elevation, velocity fields and velocity profiles in both horizontal and vertical components, showed the capability of the newly developed model in predicting wave interaction with permeable submerged breakwater.