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

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

Народ привет!

Есть проблема с инсталляцией.... кто может помочь пишите в личку или в аську.

Спасибо!

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


хочу сделать расчет движения воздуха по трубе длиной 1,5 км и диаметром=1м. Один вход, один выход, на входе воздух загоняется в трубу со статич. Давлением =1000 Ра

Однако лихо (см. размерность давления) ... на входе ~0.01 кг/см2...

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

по меньшей мере это странно:

при таких начальных давлениях P1=1000 Pa и P1*=1550.5 Pa их отношение равно 0.645, что для воздуха (k=1.4) соответствует числу Маха M=0.817. С другой стороны, имея объемный расход Q=23.7 м^3/сек через диаметер d=1м, получаем скорость u=30.2 m/s. Теперь, чтобы получить такое число Маха при такой скорости, необходимо иметь температуру T=3.4 K, как довести воздух до такого состояния я не представляю. :-)

Во-вторых, при движении воздуха по трубе постоянного сечения, при наличии только трения, скорость газа будет увеличиваться, а давление соответственно, уменьшяться.

Например, если принять начальные условия M1=0.09 (u1=31.2 m/s), T1*=300 K, P1*=1550.5 Pa (P1=1541.7), d=1 m, то из одномерных уравнений можно получить следующее: на длине L=1500 m u2=40.77 m/s, M2=0.117, P2*=1191.57 Pa, P2=1180.13 Pa

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

На всякий случай напомню важные моменты относительно Fluent.

1. Абсолютное давление равно Operating Pressure + Gauge Pressure. Gauge Pressure - это то, что задается в граничных условиях, выводится на картинках и т.д. (оно часто называется просто Pressure). Operating Pressure задается в Define->Operating Conditions и по-умолчанию равно 101325 Па. (Судя по Вашим параметрам, наверное у Вас Operating Pressure близко 101325 Па, а значения давления Вы приводите для Gauge Pressure).

2. Если использовать для плотности закон идеального газа, то, во первых, скорость к выходу будет расти, поскольку плотность из-за падения давления будет уменьшаться, а во вторых, в процессе расчета, пока решение еще не установилось, может быть дисбаланс массового расхода между входом и выходом - для установления нужно просто подольше посчитать.

А какие граничные условия стоят у Вас на входе и выходе? И какая задана температура на входе, выходе и в начальных условиях?

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

по меньшей мере это странно:

при таких начальных давлениях P1=1000 Pa и P1*=1550.5 Pa их отношение равно 0.645, что для воздуха (k=1.4) соответствует числу Маха M=0.817. С другой стороны, имея объемный расход Q=23.7 м^3/сек через диаметер d=1м, получаем скорость u=30.2 m/s. Теперь, чтобы получить такое число Маха при такой скорости, необходимо иметь температуру T=3.4 K, как довести воздух до такого состояния я не представляю. :-)

Во-вторых, при движении воздуха по трубе постоянного сечения, при наличии только трения, скорость газа будет увеличиваться, а давление соответственно, уменьшяться.

Например, если принять начальные условия M1=0.09 (u1=31.2 m/s), T1*=300 K, P1*=1550.5 Pa (P1=1541.7), d=1 m, то из одномерных уравнений можно получить следующее: на длине L=1500 m u2=40.77 m/s, M2=0.117, P2*=1191.57 Pa, P2=1180.13 Pa

Привет, ну что-то совсем не то вы написали. :smile:

как так? скорость на входе значит = u1=31.2 m/s, P1*=1550.5 Pa согласна, а вот на счет P1=1541.7 это не верно. эт значит что динамич. давление =8,8 Ра. вроде динамическое давление на входе равно 598,66 при постоянной плотности=1,23 и при такой скорости.

На всякий случай напомню важные моменты относительно Fluent.

1. Абсолютное давление равно Operating Pressure + Gauge Pressure. Gauge Pressure - это то, что задается в граничных условиях, выводится на картинках и т.д. (оно часто называется просто Pressure). Operating Pressure задается в Define->Operating Conditions и по-умолчанию равно 101325 Па. (Судя по Вашим параметрам, наверное у Вас Operating Pressure близко 101325 Па, а значения давления Вы приводите для Gauge Pressure).

2. Если использовать для плотности закон идеального газа, то, во первых, скорость к выходу будет расти, поскольку плотность из-за падения давления будет уменьшаться, а во вторых, в процессе расчета, пока решение еще не установилось, может быть дисбаланс массового расхода между входом и выходом - для установления нужно просто подольше посчитать.

А какие граничные условия стоят у Вас на входе и выходе? И какая задана температура на входе, выходе и в начальных условиях?

Здраствуйте, спасибо за советы. :smile:

посоветуйте пожалуйта если сталкивались вот на какую тему:

задача ротор-статор. для взаимодействия выбрана модель miming plane. на статоре часть проточной зоны закрыта заслонками. и получается много обратного течения при решении задачи на границе состыковки ротор -статор. и как следствие не выполняется условие сохранения массы вещества.

как можно поступить в таком случае. как сохранить массовый расход?

гр. услович такие:

на вход ротора- pressure inlet=атм=101325

выход ротор -pressure outlet (mix plane)

на вход статор pressure inlet (mix plane)

выход статор -pressure outlet =102325 Ра

в панели Operating Conditions Operating Pressure =101325 Па

может поставить на входе в область статора не pressure inlet а mass flow inlet как условие сохранения массы вещества?

Изменено пользователем Fl-lena
Ссылка на сообщение
Поделиться на других сайтах

Добрый день!

Диспут приобретает немного бестолковый характер, но умолчать не могу :-)

Привет, ну что-то совсем не то вы написали.

как так? скорость на входе значит = u1=31.2 m/s, P1*=1550.5 Pa согласна, а вот на счет P1=1541.7 это не верно. эт значит что динамич. давление =8,8 Ра. вроде динамическое давление на входе равно 598,66 при постоянной плотности=1,23 и при такой скорости.

Значит так. Если указанные Вами давления рассматривать как избыточные к атмосферному, то при M1=0.09, T1.tot=300 K,

P1.tot=(1550.5+101325) Pa, для воздуха k=1.4, получаем,

что P1=(1+(k-1)/2*M1**2)**(-k/(k-1))*P1.tot=102294.316 Pa, T1=(1+(k-1)/2*M1**2)**(-1)*T1.tot=299.515 K.

Более того, из уравнения состояния rho1=P1/(R*T1)=1.188 kg/m**3, где R=287.4 J/(kg*K). В этом случае динамическое давление rho1*u1**2/2=580.009 Pa.

Надеюсь, с этим Вы согласитесь? ;-)

Что касаемо второго вопроса.

Не совсем понятна постановка задачи. Может быть иллюстрация что-то прояснила?

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

Здравствуйте.

Модель Mixing Plane у Fluenta очень капризна относительно сходимости, что связано с явной схемой сопряжения решений в разных подобластях, поэтому для начала убедитесь в том, что решение действительно сошлось. Если сходимость есть, а дисбаланс большой, то можно попробовать Ваш вариант с массовым расходом (правда в этом случае может получится большой разрыв полного давления). Если сходимости нет, то попробуйте уменьшить коэффициент релаксации. Еще сходимость может ухудшаться, если взято слишком много точек на профилях Mixing Plane (их должно быть не больше, чем количество ячеек вдоль выбранного направления на границе, а лучше немного меньше).

Вообще, я лично стараюсь всячески избегать использования mixing plane model - слишком уж много недостатков у ее технической реализации. Про один из них я уже писал выше - это явная схема сопряжения, другой - требование заранее определиться с направлением потока. То есть если в течении появляются значительные обратные токи на границе, то считать нельзя. В этом случае нужно попробовать MRF или Sliding Mesh. А лучше всего сделать так: для начала считаем MRF, и быстро получаем стационарное решение, а затем переключаемся на нестационар, включаем вращение сетки и считаем Sliding Mesh.

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

Здравствуйте.

Модель Mixing Plane у Fluenta очень капризна относительно сходимости, что связано с явной схемой сопряжения решений в разных подобластях, поэтому для начала убедитесь в том, что решение действительно сошлось. Если сходимость есть, а дисбаланс большой, то можно попробовать Ваш вариант с массовым расходом (правда в этом случае может получится большой разрыв полного давления). Если сходимости нет, то попробуйте уменьшить коэффициент релаксации. Еще сходимость может ухудшаться, если взято слишком много точек на профилях Mixing Plane (их должно быть не больше, чем количество ячеек вдоль выбранного направления на границе, а лучше немного меньше).

Вообще, я лично стараюсь всячески избегать использования mixing plane model - слишком уж много недостатков у ее технической реализации. Про один из них я уже писал выше - это явная схема сопряжения, другой - требование заранее определиться с направлением потока. То есть если в течении появляются значительные обратные токи на границе, то считать нельзя. В этом случае нужно попробовать MRF или Sliding Mesh. А лучше всего сделать так: для начала считаем MRF, и быстро получаем стационарное решение, а затем переключаемся на нестационар, включаем вращение сетки и считаем Sliding Mesh.

Спасибо большое за совет, уже считаю с MRF :smile:

можно при решении одной задачи использовать и mixing plane и MRF? :blink:

Изменено пользователем Fl-lena
Ссылка на сообщение
Поделиться на других сайтах

Если Mixing Plane и MRF-интерфейсы заданы на разных границах, то можно.

если вам не трудно подскажите еще вот что:

что то не так, массовый расход через интерфейсы для MRF = 0. почему так может быть?

я связала выходную поверхность от секрора (периодический кусок области от ротора) и полную поверхность на входе в статор. может так нельзя?

(правда сетка на одной поверхности в два раза мельче чем на другой, но вроде это допустимо)

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

если вам не трудно подскажите еще вот что:

что то не так, массовый расход через интерфейсы для MRF = 0. почему так может быть?

я связала выходную поверхность от секрора (периодический кусок области от ротора) и полную поверхность на входе в статор. может так нельзя?

(правда сетка на одной поверхности в два раза мельче чем на другой, но вроде это допустимо)

По-моему, для MRF и Sliding периоды сетки должны быть одинаковые с обеих сторон от интерфейса.

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Aiche
      Да, сменили, ошибка ушла, референтные получается не правильно сделали, не обнуляли 4 и 5 бит, все остальное сделали по инструкции и ошибки референсов все равно ушли, просто в 0 выставили и перезагрузили, теперь другая проблема, программа не запускается и инструмент не сменяется по командам М16, М19.
    • mmaag
      Сварные базовые детали станков и машин. Л.П.Черкассая, Л.М.Финкель, 1981
    • lux59
      Доброго всем дня! Подскажите, а можно в солиде поделить линию на равные отрезки, сделать я так могу 2 способами, на линии расставить точки и выставить равенство или сразу рисовать линию из нескольких отрезков и сделать им равенство, есть ли проще путь? В двух вариантах страдает редактируемость. В Sketch AP делается это просто, можно так же в солиде сделать?
    • ДОБРЯК
      Если до решения сделать такой сдвиг, для вырожденной матрицы [К] [К]+w*[М], то эта матрица уже не будет вырождена. То после решения вместо нулевой первой собственной частоты, какую частоту мы получим?   Например w=100, какая будет первая собственная частота вместо нуля?
    • Fedor
      https://ru.wikipedia.org/wiki/LU-разложение  Вот что делается... :)   https://en.wikipedia.org/wiki/LU_decomposition
    • ДОБРЯК
    • Горыныч
      Только Вы решаете. Я Вам из собственного опыта привел рекомендации. Если страшно, то не берите, только и всего. Сломать можно все что угодно. 
    • Orchestra2603
      Ладно. Я примерно понял, что вы имеете в виду. В теории, даже чтобы выделить n-ный базис нужно зафиксировать какие-то значения n раз. Только надо так их задавать, чтобы каждый раз получались линейно-независимые вектора. Я бы не называл это "закреплением", потому что возникает ложная кажущаяся параллель со статическим решением, где мы ищем единственное решение и для этого фиксируем какие-то компоненты, хотя здесь цель расчета принципиально иная. Это, как по мне, приводит к путанице, но если вам нравится это так называть, то - пожалуйста. Кто я такой, чтобы вам запрещать)   Я не могу сказать, что решать нужно всегда. В моем понимании в определенных итерационных схемах собственные вектора получаются полседовательными приближениями сами по себе. Применяя всякие манипуляции к матрицам (сдвиги и т.д.) можно контролировать сходимость к наибольшему, наименьшему, i-му и т.д. собственному вектору. Я не прав?   Что касается факторизаций, то для любых квадратных матриц существует разложение Шура, где ортогональными преобразованиями матрица приводится к треугольной форме. В частности, для обобщенной задачи, вот выдержка из статьи. Алгоритм получаения такого разложения для обобщенной задачи называется QZ алгоритм. Он вовсю используется в Матлабею Я не вижу никаких препятствий принципиальных, чтобы получить такую факторизацию матриц.     Т.е. похоже, что, да, все же приходится решать СЛАУ, в общем случае для треугольной матрицы A-lambda*B. И покуда lambda будет собственным значением, вся матрица A-lambda*B будет вырождена. В частности, для нулевых СЗ A-lambda*B = A, и А, сама собой, вырождена.   Но.. не могу пока сказать точно, но вроде как если A и B - симметричны, то QAZ и QBZ из треугольных превращаются в диагональные, и тогда столбцы Q становятся автоматически собственными векторами (разложение Шура переходит в спектральное), т.е. и решать ничего не приходится. Но где-то пишут, что нужна положительная определенность, где-то пишут, что только симметрии достаточно... не могу пока точно сказать.  
    • ДОБРЯК
      Я вижу, что вы не поняли. А говорили, что это знали еще в детском саду... https://ru.wikipedia.org/wiki/%D0%9C%D0%B5%D1%82%D0%BE%D0%B4_%D0%93%D0%B0%D1%83%D1%81%D1%81%D0%B0 Читайте внимательно. :=)   Ме́тод Га́усса — классический метод решения системы линейных алгебраических уравнений (СЛАУ). Назван в честь немецкого математика Карла Фридриха Гаусса. Это метод последовательного исключения переменных, когда с помощью элементарных преобразований система уравнений приводится к равносильной системе треугольного вида, из которой последовательно, начиная с последних (по номеру), находятся все переменные системы.  
    • Shura762
      (дисклеймер ) не взлетит... особенно учитывая его отношение, как человека колбасит, и почему в ТФ "узел" а не точка... перепутано местами "ты" , "на" и пропустили слово "пошёл" ... можно сказать я работаю в таком коллективе, и ЧЁ, плевали на ТФ с высокой колокольни, если молодежь можно еще как-то замотивировать, показать как надо работать, то старшее поколение делает на "отвали моя черешня" чертежи выдали тяп-ляп и что тебе надо то еще .... (дальше идиоматические выражения в рамках законодательства) p.s. я догадываюсь почему "узел" и считаю что правильное название дали.
×
×
  • Создать...