Перейти к публикации

Расчет коэффициента лобового сопротивления (Cx) простых тел и сравнение полученного результата с экспериментом


dmiantosha

Рекомендованные сообщения

Всем Привет.

 

Проделал некоторый труд, продул и рассчитал Cx простых тел, сравнил с экспериментом.

 

Решил оформить все в этот топик, возможно некоторым будем полезно, как говорится я просто оставлю это здесь....

 

Все тела вращения продувал в цилиндре, брал сектор в 10 градусов. Создавал сначала плоскую блочную сетку на одной из стороне симметрии, затем либо выдавливал вращением блоки, либо сразу элементы в 6 слоев. 

Куб и куб повернутой ребром к потоку считал в половине параллелепипеда. Сетка сразу в 3D Blocking.

 

Все тела вращения имеют радиус r=0.1 метр.

Кубы со стороной a = 0.2 метра.

То есть примерно подобного объема тела.

 

post-11159-0-94888700-1423175988_thumb.jpg

post-11159-0-59280400-1423176001_thumb.jpg

 

Считал в Isotermal режиме, рабочее тело Air At 25C. Лучший результат дала модель тубрулентности k-epsilon, от интенсивности (low, medium, high) результат почти не зависел. 

Inltet->Subsonic->Normal Speed->5 м/с

Outlet->Average Static Pressure->0 Pa

Боковые поверхности цилиндра - Wall->Free Slip

Само тело - > Wall->No Slip

Симметрия на соответствующие грани.

 

Добавлял Mesh Adaptation по скорости в CFX-Pre, где было не много элементов 3-4 итерации, для кубов 1-2.

 

В итоге в расчетах получал подобную картину:

Сетка:

post-11159-0-02043300-1423176362_thumb.jpg

CFD-Post:

post-11159-0-97162800-1423176405_thumb.jpg:

 

Итоговая таблица с результатами:

post-11159-0-85996000-1423176437_thumb.jpg

 

Расчет с экспериментом сходится достаточно точно почти для всех тел, кроме сферы (я не знаю в чем тут дело, просто опустились руки) и конуса. У конуса Cx=0.5 при таком определение: "конус (2:1)", что это такое за соотношение 2 к 1, я не понял, взял 2h=D, возможно в этом причина. С каплевидном телом сложнее, само понятие слишком расплывчатое, нарисовал что-то вроде капли, но это с потолка, поэтому можно считать результат верным.

 

Также выкладываю сами файлы WB, возможно кому-то пригодятся.

 

Модель для тел вращения: http://rusfolder.com/42917782

Модель для Кубов: http://rusfolder.com/42917783

 

П.С: Если кому-то топик оказался полезным, приму благодарности на WMR: R275859695001  :smile: 

Ссылка на сообщение
Поделиться на других сайтах


Органолептически чую (а проверять влом): сопротивление сферы меньше сопротивления полусферы а у Вас по "экспериментальным" данным наоборот. Вот!

Ссылка на сообщение
Поделиться на других сайтах

Пожалуй для общего развития поинтересуюсь......... это коэффициент от набегающего потока, а как то учитывается статика тыльной стороны ? по аналогии с поршнем автомобиля - подвижная стенка с разностью давлений

Ссылка на сообщение
Поделиться на других сайтах

Пожалуй для общего развития поинтересуюсь......... это коэффициент от набегающего потока, а как то учитывается статика тыльной стороны ? по аналогии с поршнем автомобиля - подвижная стенка с разностью давлений

Строго говоря зависит от скорости  :biggrin:

 

На гиперзвуке почти не зависит, поэтому во всяких re-entry задачах (по русски - задача входа в атмосферу)  обычно считают только "переднюю" (forebody) часть тела ...

Ссылка на сообщение
Поделиться на других сайтах

Makarrr, сопротивление полусферы меньше, чем сопротивление сферы, если Вам проверять "влом", то поверьте на слово.

 

MotorManiac, Cx это отношение Силы лобового сопротивления к набегающему потоку. Лобовое сопротивление есть иксовая составляющая полной аэродинамической силы действующей на тело (составляющая вдоль оси Y - подъемная сила). Полная аэродинамическая сила есть интеграл распределения давления по всей поверхности тела, поэтому учитываются все площади.

Ссылка на сообщение
Поделиться на других сайтах
  • 4 недели спустя...

Всем Привет.

 

Проделал некоторый труд, продул и рассчитал Cx простых тел, сравнил с экспериментом.

 

Решил оформить все в этот топик, возможно некоторым будем полезно, как говорится я просто оставлю это здесь....

 

Все тела вращения продувал в цилиндре, брал сектор в 10 градусов. Создавал сначала плоскую блочную сетку на одной из стороне симметрии, затем либо выдавливал вращением блоки, либо сразу элементы в 6 слоев. 

Куб и куб повернутой ребром к потоку считал в половине параллелепипеда. Сетка сразу в 3D Blocking.

 

Все тела вращения имеют радиус r=0.1 метр.

Кубы со стороной a = 0.2 метра.

То есть примерно подобного объема тела.

 

attachicon.gifsetka.jpg

attachicon.gifsetkaover.jpg

 

Считал в Isotermal режиме, рабочее тело Air At 25C. Лучший результат дала модель тубрулентности k-epsilon, от интенсивности (low, medium, high) результат почти не зависел. 

Inltet->Subsonic->Normal Speed->5 м/с

Outlet->Average Static Pressure->0 Pa

Боковые поверхности цилиндра - Wall->Free Slip

Само тело - > Wall->No Slip

Симметрия на соответствующие грани.

 

Добавлял Mesh Adaptation по скорости в CFX-Pre, где было не много элементов 3-4 итерации, для кубов 1-2.

 

В итоге в расчетах получал подобную картину:

Сетка:

attachicon.gifmesh.jpg

CFD-Post:

attachicon.gifpost.jpg:

 

Итоговая таблица с результатами:

attachicon.giftable.jpg

 

Расчет с экспериментом сходится достаточно точно почти для всех тел, кроме сферы (я не знаю в чем тут дело, просто опустились руки) и конуса. У конуса Cx=0.5 при таком определение: "конус (2:1)", что это такое за соотношение 2 к 1, я не понял, взял 2h=D, возможно в этом причина. С каплевидном телом сложнее, само понятие слишком расплывчатое, нарисовал что-то вроде капли, но это с потолка, поэтому можно считать результат верным.

 

Также выкладываю сами файлы WB, возможно кому-то пригодятся.

 

Модель для тел вращения: http://rusfolder.com/42917782

Модель для Кубов: http://rusfolder.com/42917783

 

П.С: Если кому-то топик оказался полезным, приму благодарности на WMR: R275859695001  :smile: 

Находили данные экспериментов (обтекания простых тел; на низких скоростях и возможно на высоких) с подробным описанием? - не книжек где просто цифры коэффициентов указанны, а данные продувок в А. трубах с указанием всех параметров потока.

Ссылка на сообщение
Поделиться на других сайтах
  • 2 года спустя...
В 6.2.2015 в 02:59, dmiantosha сказал:

Всем Привет.

 

Проделал некоторый труд, продул и рассчитал Cx простых тел, сравнил с экспериментом.

 

Решил оформить все в этот топик, возможно некоторым будем полезно, как говорится я просто оставлю это здесь....

 

Все тела вращения продувал в цилиндре, брал сектор в 10 градусов. Создавал сначала плоскую блочную сетку на одной из стороне симметрии, затем либо выдавливал вращением блоки, либо сразу элементы в 6 слоев. 

Куб и куб повернутой ребром к потоку считал в половине параллелепипеда. Сетка сразу в 3D Blocking.

 

Все тела вращения имеют радиус r=0.1 метр.

Кубы со стороной a = 0.2 метра.

То есть примерно подобного объема тела.

 

post-11159-0-94888700-1423175988_thumb.jpg

post-11159-0-59280400-1423176001_thumb.jpg

 

Считал в Isotermal режиме, рабочее тело Air At 25C. Лучший результат дала модель тубрулентности k-epsilon, от интенсивности (low, medium, high) результат почти не зависел. 

Inltet->Subsonic->Normal Speed->5 м/с

Outlet->Average Static Pressure->0 Pa

Боковые поверхности цилиндра - Wall->Free Slip

Само тело - > Wall->No Slip

Симметрия на соответствующие грани.

 

Добавлял Mesh Adaptation по скорости в CFX-Pre, где было не много элементов 3-4 итерации, для кубов 1-2.

 

В итоге в расчетах получал подобную картину:

Сетка:

post-11159-0-02043300-1423176362_thumb.jpg

CFD-Post:

post-11159-0-97162800-1423176405_thumb.jpg:

 

Итоговая таблица с результатами:

post-11159-0-85996000-1423176437_thumb.jpg

 

Расчет с экспериментом сходится достаточно точно почти для всех тел, кроме сферы (я не знаю в чем тут дело, просто опустились руки) и конуса. У конуса Cx=0.5 при таком определение: "конус (2:1)", что это такое за соотношение 2 к 1, я не понял, взял 2h=D, возможно в этом причина. С каплевидном телом сложнее, само понятие слишком расплывчатое, нарисовал что-то вроде капли, но это с потолка, поэтому можно считать результат верным.

 

Также выкладываю сами файлы WB, возможно кому-то пригодятся.

 

Модель для тел вращения: http://rusfolder.com/42917782

Модель для Кубов: http://rusfolder.com/42917783

 

П.С: Если кому-то топик оказался полезным, приму благодарности на WMR: R275859695001  :smile: 

интересный пост,жаль нет продувки конусов,пирамид и прочих тел - поперек потока, например при 100км в час,было бы любопытно посмотреть результаты..

что касается каплевидной формы, она повторяет форму капли дождя, летящей вниз к земле, вода принимает форму наиболее обтекаемую!

Ссылка на сообщение
Поделиться на других сайтах
  • 1 месяц спустя...
В 06.02.2015 в 01:59, dmiantosha сказал:

Добавлял Mesh Adaptation по скорости в CFX-Pre, где было не много элементов 3-4 итерации, для кубов 1-2.

А как адаптацию подключили? Можете показать в каком это разделе пре-процессора? 

Ссылка на сообщение
Поделиться на других сайтах

http://docs.cntd.ru/document/1200084848   В СП немного другие цифры особенно с учетом наветренного и подветренного давлений... :) 

Ссылка на сообщение
Поделиться на других сайтах
  • 3 года спустя...
Богдан Романович

Добрый день. Проконсультируйте пожалуйста как посчитать коэффициент лобового сопротивления просечно-витяжной сетки которую устанавливают в бампер автомобиля для защиты радиатора. Можно ее ставить, нельзя, будет она мешать воздушному потоку и ухудшать охлаждение радиатора или не будет, будет ли сетка на скорости приводить к завихрениям воздушного потока ухудшая тем самым доступ воздуха к радиатору ? С завода сетки не ставят, ставят уже сами автолюбители для защиты радиатора от мусора и крупных насекомых. Спасибо за ответ.

Ссылка на сообщение
Поделиться на других сайтах
28.06.2021 в 01:19, Богдан Романович сказал:

коэффициент лобового сопротивления просечно-витяжной сетки

 

Идельчик И. Е. Справочник по гидравлическим сопротивлениям там все есть.

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

Вы можете опубликовать сообщение сейчас, а зарегистрироваться позже. Если у вас есть аккаунт, войдите в него для написания от своего имени.
Примечание: вашему сообщению потребуется утверждение модератора, прежде чем оно станет доступным.

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • Zergus
      Выкладывайте сюда свои "чудеса", может по модели получится понять, что в ней не так. А если есть возможность писать в техподдержку - кидайте им свой файл, думаю разберутся.
    • Богоманшин Игорь
      А есть подтверждение в каких либо нормативных документах? Или только из соображений "максимально минимальная себестоимость"?
    • ANT0N1DZE
      Это не поможет ответить на этот вопрос конкретными числами. Судя по всему, там металлические частички внутри материала. Подобные эластомеры существуют и применяются. 
    • gudstartup
      как перевести шпиндель снова первоначальный режим управления скоростью когда ему пофиг число импульсов на оборот? резьбу пока резать не будут а в подпрограмму смены можно что нить добавить или в реализацию m19  чтобы чпу из режима контроля позиционного кодера вышло?
    • gudstartup
    • mactepnew
      Вот я тоже так думал, пока не столкнулся с обратным) Как еще столкнусь, пришлю скриншот. Хаос еще в том, что порядок в дереве построения весьма условный. операция, которая стоит перед групой, имеет бОльший timestamp(пардон, не знаю как это в переводе), чем эта группа, и выполняется позже. При этом все операции в группе имеют timestamp на пару сотен меньше и выполняются раньше, чем две группы, стоящие перед ней. Логика?  Если в эскизе все элементы определены, то можно поставить еще один размер (хочешь переопределить эскиз - да пожалста!), изменить его - и NX самостоятельно выберет, какой размер или ограничение отпустить и весь эскиз летит к чертям. Это ли не хаос?   Назовете хоть один CAD, которым можно пользоваться так? Любой из встреченных мной имеет много неочевидных нюансов, которые конечно не проявляются при выполнении университетского или базового курса. В Creo я на протяжении 5 лет находил всё новые фишки, которые не прописаны в документации. Огромная благодарность Ruslan за пояснения тонкостей. Даже Компас и Solidworks, которые считаются самыми дружественными, оказалось непросто обуздать, когда проходил собеседования и делал тестовые. А ведь университетский курс знал на отлично!   Возвращаясь к NX - вот сегодняшний пример. Утром деталь была оранжевая, единая. Изменил пару уклонов и толщин, стала желтая с элементами оранжевого и частично твердотельной, частично листовой. Я даже не знаю, как сформулировать вопрос в поддержку. 
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • kkk
      Короче, как сделать автоматом не нашел. Раскрасил нужные грани нужными цветами. Затем на вкладке с цветом модели можно выбирать нужный цвет, тогда в модели выбираются грани, покрашенные этим цветом. Можно еще сохранить выбор этих граней, но смысла не вижу, т.к. все легко выбирается по цвету. Затем либо в инструменте Измерить смотрим площадь и как-то ее используем, либо можно макросом (что выкладывал Kir95 ) внести площадь в свойства (на вкладку Конфигурации). Можно еще внести выбранные грани в датчики, но как использовать их в свойствах не нашел.    Макрос прикрепляю.     Площадь выбранных граней в свойства.swp
    • IgnatvKudrin
      Для нашего производства требуются постоянно инструментальные пружины для штампов. Сейчас заказы обрабатываются довольно долго - до 2-х месяцев. А у нас сроки на штампы 1-2 месяца с испытанием, доработкой и т.п. Хотелось бы найти поставщика, у которого уже есть склад с данными пружинами или кто может доставить за 1 месяц.   Пружины стандартные по ISO, но на всякий случай приложил каталог. Вот список пружин, которые нужны очень срочно:   LR 32 x 44 - 1 шт. LR 25 x 51 - 2 шт. LB 63 x 89 - 7 шт. LB 32 x 102 - 10 шт. LB 32 x 89 - 20 шт. LB 32 x 76 - 12 шт. LB 32 x 51 - 6 шт. LB 25 x 127 - 9 шт.   Готов выслушать все предложения. Возможно даже есть те, кто сможет изготовить подобные пружины с такими же параметрами. 20231101.pdf
    • bubblyk
      Знатоки форума, подскажите, как смоделировать простую (на первый взгляд) задачу термомеханики в упругой постановке.   Есть брусок правильный формы ака кирпич, в котором генерируется НЕРАВНОМЕРНОЕ удельное тепловыделение (HGEN). Задняя и передняя поверхности конвективно охлаждаются с известными (но различными) коэффициентами теплоотдачи и температурами жидкости/газа. Свойства кирпича (теплофизика + механика) принимаются зависящими от температуры. Постановка - стационарная Собственно вопрос : как грамотно смоделировать температурное поле в кирпиче с учетом температурного расширения?   Исходный код приложил.     Надеюсь на коллективный разум и опыт! EPTH2D.mac p.s.   Если есть опыт решения этой задачи в Ansys WB, тоже можно сюда приложить командные вставки.
×
×
  • Создать...