Что За Ошибка В Cfx?

[ajax 0]

Расчет вылетел вот с такой ошибкой, кто нибудь может сказать чем может быть вызвана?

OUTER LOOP ITERATION = 33 CPU SECONDS = 1.271E+04

----------------------------------------------------------------------

| Equation | Rate | RMS Res | Max Res | Linear Solution |

+----------------------+------+---------+---------+------------------+

+--------------------------------------------------------------------+

| ERROR #001100279 has occurred in subroutine ErrAction. |

| Message: |

| Floating point exception: Overflow |

| |

| |

| |

| |

| |

+--------------------------------------------------------------------+

+--------------------------------------------------------------------+

| ERROR #001100279 has occurred in subroutine ErrAction. |

| Message: |

| Stopped in routine F_FPX_HANDLER |

| |

| |

| |

| |

| |

+--------------------------------------------------------------------+

+--------------------------------------------------------------------+

| An error has occurred in cfx5solve: |

| |

| The ANSYS CFX solver has terminated without writing a results |

| file. Command on host roman exited with return code 0. |

+--------------------------------------------------------------------+

End of solution stage.

+--------------------------------------------------------------------+

| The following user files have been saved in the directory |

| C:\cfx1\3_etap\evaporation\23_10\variant1\var1_001: |

| |

| PT1.trk, PT2.trk, mon |

+--------------------------------------------------------------------+

This run of the ANSYS CFX Solver has finished.

[Dick 0]

Из-за OVERFLOW вылетает. Причин может быть множество - от низкого качества сетки, до неправильных ГУ...

[digi 0]

А какое качество сетки считается приемлимым, чтобы из-за него эта ошибка не возникала? (я имею в виду минимальное допустимое положительное значение качества элементов сетки)

[Dick 0]

Когда солвер запускается на расчет в окне output file в самом начале пробегает статистика по сетке, например, вот такая:

+--------------------------------------------------------------------+

| Mesh Statistics |

+--------------------------------------------------------------------+

Domain Name : Default Domain

Total Number of Nodes = 47800

Total Number of Elements = 23562

Total Number of Hexahedrons = 23562

Total Number of Faces = 47798

Minimum Orthogonality Angle [degrees] = 85.3 OK

Maximum Aspect Ratio = 33.2 OK

Maximum Mesh Expansion Factor = 1.2 OK

Если написано "ОК" - значит сетка хорошая.

Если сетка "приемлемого" качества вместо "ОК" будут "ок". Если "низкого" качества - вместо "ОК" будут "!"

[dmiantosha 11]

А если качество на высоком уровне, но сетка редкая? Есть ли тут какой-то фактор. Имеется ввиду можно ли как-нибудь хотя бы косвено расчитать приемлимое число элементов, или в этом процессе можно полагаться только на опыт и интуицию?

[a_schelyaev 366]

Зависит от задачи. Многосеточная решалка любит монотонное решение.

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

Поэтому, только опыт и интуиция.

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

:)

[dmiantosha 11]

а на каких других то? Ведь решение можно получить толкьо в узле??? Между узлами решения то нет??? Поэтому пучть в узле А - 10 Па, в узле Б - тоже 10 Па, а между ними скачок в 20 Па. И как быть?

[digi 0]

Посмотрел статистику. Minimum Orthogonality Angle [degrees] = 16.3 ! и Maximum Mesh Expansion Factor = 3517.2 !. Как это можно исправить? При решении со средой "воздух" всё идёт нормально. Только меняю на "воду" - вылетает overflow. Может проблема не только в сетке?

[Parma 0]

Тут недавно наткнулся в хелпе на фразу: аспект ратио в пограничном слое может быть до величин 10в5 - 10в6 степени. Главное чтоб длинная сторона смотрела по потоку.

Сам ситаю сеточку с аспектом 15000 в хвосте погранслоя, ругается но считает. А шо делать, память от не резиновая...

С углом хуже - для каждого метода есть свой угол, меньше которого просто невозможно получить что-то похожее на физику. Сам стараюсь меньше 30 градусов не делать - полет нормальный, ну почти... ;)

[digi 0]

Люди! Как добиться нормального Minimum Orthogonality Angle и Maximum Mesh Expansion Factor? Это что, опять сетку перестраивать? Или refinement использовать? Подскажите пожалуйста!

[a_schelyaev 366]

а на каких других то? Ведь решение можно получить толкьо в узле??? Между узлами решения то нет??? Поэтому пучть в узле А - 10 Па, в узле Б - тоже 10 Па, а между ними скачок в 20 Па. И как быть?

В других это в FlowVision. Глядя на какой сетке народ задачи считает, я мало понимаю почему это не сделать на FlowVision. Только разве что из-за работы на приатском софте.

Скачок 20 Па ерунда. Тут можно никак не быть.

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

У CFX на грубой тетра сетке вообще могут чудеса твориться, поэтому он очень любит подробную тетра сетку.

Люди! Как добиться нормального Minimum Orthogonality Angle и Maximum Mesh Expansion Factor? Это что, опять сетку перестраивать? Или refinement использовать? Подскажите пожалуйста!

Братан, нужно в генераторе сетки провести сглаживание по качеству или аспектному соотношению или углам.

Берешь ICEM CFD и спишь спокойно.

По идее он в процессе построения сетки уже проводит 5 итераций (по умолчанию) по сглаживанию сетки.

Только если у тебя не узкие места какие-нибудь, куда генератор как смог со скрипом сетку засунул.

[Dick 0]

Сам ситаю сеточку с аспектом 15000 в хвосте погранслоя, ругается но считает

Если "поставить галку" сжимаемый газ - считать скорее всего перестанет.

Посмотрел статистику. Minimum Orthogonality Angle [degrees] = 16.3 ! и Maximum Mesh Expansion Factor = 3517.2 !. Как это можно исправить? При решении со средой "воздух" всё идёт нормально. Только меняю на "воду" - вылетает overflow. Может проблема не только в сетке?

См. мой первый пост.

Надеюсь понятно, что вода и воздух немного различаются по свойствам. Соответственно, чувствительность к косячности сетки тоже немного другая.

Поднять качество по углу поможет сглаживание. Чтоб был нормальный Expansion factor - желательно генерить более или менее изотропную сетку.

[a_schelyaev 366]

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

[Dick 0]

На геометрии с непростой топологией накосячить очень легко.

Maximum Mesh Expansion Factor = 3517.2 !. - говорит о том, что максимальный элемент в 3517 раз больше минимального, т.е. если мимнимальный элемент порядка 0.1 мм, то максимальный будет порядка метра. В принципе это не настолько критично, в задачах внешней аэродинамики такая сетка вполне обычное дело. А вот, например, в случае многофазных течений или химических реакций сетка должна быть практически изотропной и такая сетка не годится.

[a_schelyaev 366]

Автор, картинку давай и можешь текстовку положить из out файла, с описанием варианта, которео в начале идет.

[digi 0]

Это моя сетка

Это файл вывода:

This run of the CFX-11.0 Solver started at 21:54:49 on 9 Dec 2008 by

user Борис on B9C2EA08AE1C463 (intel_p4.sse2_winnt5.1) using the

command:

"C:\Program Files\ANSYS Inc\v110\CFX\bin\perllib\cfx5solve.pl"

-stdout-comms -batch -ccl -

Setting up CFX Solver run ...

+--------------------------------------------------------------------+

| |

| CFX Command Language for Run |

| |

+--------------------------------------------------------------------+

LIBRARY:

MATERIAL: Water

Material Description = Water (liquid)

Material Group = Water Data, Constant Property Liquids

Option = Pure Substance

Thermodynamic State = Liquid

PROPERTIES:

Option = General Material

Thermal Expansivity = 2.57E-04 [K^-1]

ABSORPTION COEFFICIENT:

Absorption Coefficient = 1.0 [m^-1]

Option = Value

END

DYNAMIC VISCOSITY:

Dynamic Viscosity = 8.899E-4 [kg m^-1 s^-1]

Option = Value

END

EQUATION OF STATE:

Density = 997.0 [kg m^-3]

Molar Mass = 18.02 [kg kmol^-1]

Option = Value

END

REFERENCE STATE:

Option = Specified Point

Reference Pressure = 1 [atm]

Reference Specific Enthalpy = 0.0 [J/kg]

Reference Specific Entropy = 0.0 [J/kg/K]

Reference Temperature = 25 [C]

END

REFRACTIVE INDEX:

Option = Value

Refractive Index = 1.0 [m m^-1]

END

SCATTERING COEFFICIENT:

Option = Value

Scattering Coefficient = 0.0 [m^-1]

END

SPECIFIC HEAT CAPACITY:

Option = Value

Specific Heat Capacity = 4181.7 [J kg^-1 K^-1]

Specific Heat Type = Constant Pressure

END

THERMAL CONDUCTIVITY:

Option = Value

Thermal Conductivity = 0.6069 [W m^-1 K^-1]

END

END

END

END

FLOW:

SOLUTION UNITS:

Angle Units = [rad]

Length Units = [m]

Mass Units = [kg]

Solid Angle Units = [sr]

Temperature Units = [K]

Time Units =

END

SIMULATION TYPE:

Option = Steady State

EXTERNAL SOLVER COUPLING:

Option = None

END

END

DOMAIN: ship

Coord Frame = Coord 0

Domain Type = Fluid

Fluids List = Water

Location = SOLID

BOUNDARY: FreeWalls

Boundary Type = WALL

Location = FreeWalls

BOUNDARY CONDITIONS:

WALL INFLUENCE ON FLOW:

Option = Free Slip

END

END

END

BOUNDARY: Inlet

Boundary Type = INLET

Location = INLET

BOUNDARY CONDITIONS:

FLOW REGIME:

Option = Subsonic

END

MASS AND MOMENTUM:

Normal Speed = 6.2 [m s^-1]

Option = Normal Speed

END

TURBULENCE:

Eddy Length Scale = 2.32e-4 [m]

Fractional Intensity = 0.05

Option = Intensity and Length Scale

END

END

END

BOUNDARY: Outlet

Boundary Type = OUTLET

Location = OUTLET

BOUNDARY CONDITIONS:

FLOW REGIME:

Option = Subsonic

END

MASS AND MOMENTUM:

Option = Static Pressure

Relative Pressure = 100 [Pa]

END

END

END

BOUNDARY: Ship

Boundary Type = WALL

Location = SHIP.1

BOUNDARY CONDITIONS:

WALL INFLUENCE ON FLOW:

Option = No Slip

END

END

END

BOUNDARY: SymP

Boundary Type = SYMMETRY

Location = SYMP

END

BOUNDARY: ship Default

Boundary Type = WALL

Location = TOP

BOUNDARY CONDITIONS:

WALL INFLUENCE ON FLOW:

Option = No Slip

END

END

END

DOMAIN MODELS:

BUOYANCY MODEL:

Option = Non Buoyant

END

DOMAIN MOTION:

Option = Stationary

END

MESH DEFORMATION:

Option = None

END

REFERENCE PRESSURE:

Reference Pressure = 1 [atm]

END

END

FLUID MODELS:

COMBUSTION MODEL:

Option = None

END

HEAT TRANSFER MODEL:

Fluid Temperature = 20 [C]

Option = Isothermal

END

THERMAL RADIATION MODEL:

Option = None

END

TURBULENCE MODEL:

Option = SST

TRANSITIONAL TURBULENCE:

Option = Fully Turbulent

END

END

TURBULENT WALL FUNCTIONS:

Option = Automatic

END

END

END

INITIALISATION:

Option = Automatic

INITIAL CONDITIONS:

Velocity Type = Cartesian

CARTESIAN VELOCITY COMPONENTS:

Option = Automatic with Value

U = 6.2 [m s^-1]

V = 0 [m s^-1]

W = 0 [m s^-1]

END

EPSILON:

Option = Automatic

END

K:

Option = Automatic

END

STATIC PRESSURE:

Option = Automatic

END

END

END

OUTPUT CONTROL:

RESULTS:

File Compression Level = Default

Option = Standard

END

END

SOLVER CONTROL:

ADVECTION SCHEME:

Option = High Resolution

END

BODY FORCES:

Body Force Averaging Type = Volume-Weighted

END

CONVERGENCE CONTROL:

Maximum Number of Iterations = 20

Physical Timescale = 2

Timescale Control = Physical Timescale

END

CONVERGENCE CRITERIA:

Residual Target = 0.001

Residual Type = RMS

END

DYNAMIC MODEL CONTROL:

Global Dynamic Model Control = On

HYDRO CONTROL:

Free Surface Harmonic Averaging = Yes

Transition Iteration = 10

END

END

EQUATION CLASS: continuity

CONVERGENCE CONTROL:

Physical Timescale = 2

Timescale Control = Physical Timescale

END

END

END

END

COMMAND FILE:

Version = 11.0

Results Version = 11.0

END

EXECUTION CONTROL:

INTERPOLATOR STEP CONTROL:

Runtime Priority = Standard

EXECUTABLE SELECTION:

Double Precision = Off

END

MEMORY CONTROL:

Memory Allocation Factor = 1.0

END

END

PARALLEL HOST LIBRARY:

HOST DEFINITION: b9c2ea08ae1c463

Installation Root = C:\Program Files\ANSYS Inc\v%v\CFX

Host Architecture String = intel_p4.sse2_winnt5.1

END

END

PARTITIONER STEP CONTROL:

Multidomain Option = Independent Partitioning

Runtime Priority = Standard

EXECUTABLE SELECTION:

Use Large Problem Partitioner = Off

END

MEMORY CONTROL:

Memory Allocation Factor = 1.0

END

PARTITIONING TYPE:

MeTiS Type = k-way

Option = MeTiS

Partition Size Rule = Automatic

END

END

RUN DEFINITION:

Definition File = C:/ship.def

Interpolate Initial Values = Off

Run Mode = Full

END

SOLVER STEP CONTROL:

Runtime Priority = Standard

EXECUTABLE SELECTION:

Double Precision = Off

END

MEMORY CONTROL:

Memory Allocation Factor = 1.0

END

PARALLEL ENVIRONMENT:

Number of Processes = 1

Start Method = Serial

END

END

END

+--------------------------------------------------------------------+

| |

| Solver |

| |

+--------------------------------------------------------------------+

+--------------------------------------------------------------------+

| |

| ANSYS CFX Solver 11.0 |

| |

| Version 2007.01.15-19.20 Mon Jan 15 19:24:21 GMTST 2007 |

| |

| Executable Attributes |

| |

| single-int32-32bit-novc6-optimised-supfort-noprof-nospag-lcomp |

| |

| Copyright 1996-2007 ANSYS Europe Ltd. |

+--------------------------------------------------------------------+

+--------------------------------------------------------------------+

| Job Information |

+--------------------------------------------------------------------+

Run mode: serial run

Host computer: B9C2EA08AE1C463

Job started: Tue Dec 09 21:55:16 2008

+--------------------------------------------------------------------+

| Memory Allocated for Run (Actual usage may be less) |

+--------------------------------------------------------------------+

Data Type Kwords Words/Node Words/Elem Kbytes Bytes/Node

Real 28605.6 267.56 288.82 111740.7 1070.24

Integer 8813.1 82.43 88.98 34426.3 329.73

Character 2820.4 26.38 28.48 2754.3 26.38

Logical 65.0 0.61 0.66 253.9 2.43

Double 1196.5 11.19 12.08 9347.9 89.53

+--------------------------------------------------------------------+

| Mesh Statistics |

+--------------------------------------------------------------------+

Domain Name : ship

Total Number of Nodes = 106913

Total Number of Elements = 99044

Total Number of Hexahedrons = 99044

Total Number of Faces = 15422

Minimum Orthogonality Angle [degrees] = 16.3 !

Maximum Aspect Ratio = 490.8 ok

Maximum Mesh Expansion Factor = 3517.2 !

+--------------------------------------------------------------------+

| Average Scale Information |

+--------------------------------------------------------------------+

Domain Name : ship

Global Length = 1.0936E+05

Minimum Extent = 4.0199E+04

Maximum Extent = 6.7420E+05

Density = 9.9700E+02

Dynamic Viscosity = 8.8990E-04

Velocity = 6.2000E+00

Advection Time = 1.7638E+04

Reynolds Number = 7.5962E+11

+--------------------------------------------------------------------+

| Checking for Isolated Fluid Regions |

+--------------------------------------------------------------------+

No isolated fluid regions were found.

+--------------------------------------------------------------------+

| The Equations Solved in This Calculation |

+--------------------------------------------------------------------+

Subsystem : Wall Scale

Wallscale

Subsystem : Momentum and Mass

U-Mom

V-Mom

W-Mom

P-Mass

Subsystem : TurbKE and TurbFreq

K-TurbKE

O-TurbFreq

CFD Solver started: Tue Dec 09 21:55:35 2008

+--------------------------------------------------------------------+

| Convergence History |

+--------------------------------------------------------------------+

======================================================================

| Timescale Information |

----------------------------------------------------------------------

| Equation | Type | Timescale |

+----------------------+------------------------+--------------------+

| U-Mom | Physical Timescale | 2.00000E+00 |

| V-Mom | Physical Timescale | 2.00000E+00 |

| W-Mom | Physical Timescale | 2.00000E+00 |

+----------------------+------------------------+--------------------+

| K-TurbKE | Physical Timescale | 2.00000E+00 |

| O-TurbFreq | Physical Timescale | 2.00000E+00 |

+----------------------+------------------------+--------------------+

======================================================================

OUTER LOOP ITERATION = 1 CPU SECONDS = 4.922E+00

----------------------------------------------------------------------

| Equation | Rate | RMS Res | Max Res | Linear Solution |

+----------------------+------+---------+---------+------------------+

| Wallscale | 0.00 | 5.8E-04 | 2.6E-03 | 23.6 2.4E-01 ok|

+----------------------+------+---------+---------+------------------+

| U-Mom | 0.00 | 3.2E+00 | 3.8E+01 | 4.7E-02 OK|

| V-Mom | 0.00 | 1.2E-05 | 1.4E-04 | 1.7E+04 ok|

| W-Mom | 0.00 | 7.1E-06 | 1.4E-04 | 5.1E+04 ok|

| P-Mass | 0.00 | 2.9E-04 | 3.4E-03 | 16.8 7.2E-02 OK|

+----------------------+------+---------+---------+------------------+

+--------------------------------------------------------------------+

| ERROR #001100279 has occurred in subroutine ErrAction. |

| Message: |

| Floating point exception: Overflow |

| |

| |

| |

| |

| |

+--------------------------------------------------------------------+

+--------------------------------------------------------------------+

| Writing crash recovery file |

+--------------------------------------------------------------------+

Details of error:-

----------------

Error detected by routine POPDIR

CRESLT = ILEG

Current Directory : /FLOW/NAMEMAP

+--------------------------------------------------------------------+

| An error has occurred in cfx5solve: |

| |

| The ANSYS CFX solver exited with return code 1. No results file |

| has been created. |

+--------------------------------------------------------------------+

End of solution stage.

+--------------------------------------------------------------------+

| Warning! |

| |

| The ANSYS CFX Solver has written a crash recovery file. This file |

| has been saved as C:\ship_001.res.err and may be an aid to |

| diagnosing the problem or restarting the run. More details should |

| be available in the solver output section of the output file. |

+--------------------------------------------------------------------+

+--------------------------------------------------------------------+

| The following user files have been saved in the directory |

| C:\ship_001: |

| |

| mon |

+--------------------------------------------------------------------+

This run of the ANSYS CFX Solver has finished.

[Dick 0]

Minimum Extent = 4.0199E+04

Maximum Extent = 6.7420E+05

длина домена 674.2 км?

PS. Сгущение сетки должно быть плавнее и равномерней.

[a_schelyaev 366]

1. С ед.изм что ли попутал?

2. Вот это не понял Free Surface Harmonic Averaging - это всегда по умолчанию?

3. Minimum Orthogonality Angle [degrees] = 16.3 ! - если мне склероз не изменяет, то для гексы минимальный угол скоса раен 20.

4. Никто не запрещал уменьшить шаг по времени и посмотреть что будет.

5. Ну и Дик прав - такие перепады по сетке просто недопустимы.

[digi 0]

1. Как бороться с перепадами сетки? Refinement? Делаю сглаживание по углу. Задаю мин. угол 20 градусов, но он не дотягивает до этого значения (поднимает до 17 градусов)

2. Ед. изм. - мм. Сам не понимаю, откуда берутся такие величины. Длина домена, по идее, 674.2 м.

3. Сглаживание по expansion factor ничего не даёт.

[Dick 0]

2. Ед. изм. - мм. Сам не понимаю, откуда берутся такие величины. Длина домена, по идее, 674.2 м.

Может все-таки стоит в препроцессоре или постпроцессоре по "линейке" со шкалой посмотреть реальный размер домена?

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