شبیه سازی پدیده Vortex Shedding با توسعه حلگر viscoelasticFluidFlow

در این شبیه سازی پدیده ی vortex shedding اطراف استوانه ی ثابت و نوسانی برای سیال ویسکوالاستیک مدلسازی شد. برای این مدلسازی نیاز شد تا کتابخانه ی شبکه ی دینامیکی (dynamic mesh) به حلگر پایه یعنی viscoelasticFluidFlow از اپن فوم نسخه 3 اضافه شود.

 توضیحات بیشتر در ادامه مطلب:

ادامه نوشته

محاسبه ی نیروی درگ در جریان ویسکوالاستیک

برای محاسبه ی نیروی های وارد بر جسم در  openfoam-ext3 ابزار libforces.so به صورت زیر به فایل controlDict از پوشه ی سیستم افزوده می شود.

  body
    {
        type                forces;
        functionObjectLibs  ( "libforces.so" );
        outputControl       timeStep;
        outputInterval      1;
        patches             (wall_c);
        pName               p;
        UName               U;
        rhoName          rhoInf;
        log                 true;
        rhoInf         803.87097;   //set density for incompressible simulation
        CofR          (0.25 0.007 0.125);  // center of rotation, it is useless for calculation of drag
    }

 در اپن فوم با این ابزار نمی توان نیروهای وارد بر جسم در سیالهای ویسکوالاستیک را محاسبه کرد. یعنی با افزودن دیکشنری بالا حین اجرای حلگر viscoelasticFluidFlow نیروها و ممنتم وارد بر جسم گزارش نمی شود. اما می توان با افزودن چند خط کد به سورس برنامه در فایل forces.C این مشکل را برطرف کرد.