Расчет рабочего колеса центробежного насоса в ansys CFX 11

[pugach 8]

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

проблема:

Создал модель канала в solidworks, экспортировал icem cfd 11 (parasolid) .

Там создал тетра сетку с призматическими элементами на лопасти и дисках РК.

Экспортировал в ansys cfx11

запускал на счет, проходит какое-то время и пишет ошибку (см. ниже). Помогите пожалуйста. Заранее благодарю.

Пробовал менять модель турбулентности на sst, пробовал изменять физическое время, пробовал сетки различной размерности, пробовал брать не канал колеса а весь поток (без условия периодичности): во всех случаях при создании файла .def никакие ошибки не обнаруживаются, расчет прекращается на первой итерации.

Текст log:

+------------------------­-------------------------­-------------------+

| |

| 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: R1

Coord Frame = Coord 0

Domain Type = Fluid

Fluids List = Water

Location = CREATED_MATERIAL_15

Number of Passages in 360 = 6

Number of Passages in Component = 1

BOUNDARY: R1 Blade

Boundary Type = WALL

Frame Type = Rotating

Location = BLADE

BOUNDARY CONDITIONS:

WALL INFLUENCE ON FLOW:

Option = No Slip

END

WALL ROUGHNESS:

Option = Rough Wall

Roughness Height = 1.6 [micron]

END

END

END

BOUNDARY: R1 Hub

Boundary Type = WALL

Frame Type = Rotating

Location = HUB.1

BOUNDARY CONDITIONS:

WALL INFLUENCE ON FLOW:

Option = No Slip

END

WALL ROUGHNESS:

Option = Rough Wall

Roughness Height = 1.6 [micron]

END

END

END

BOUNDARY: R1 Inlet

Boundary Type = INLET

Frame Type = Stationary

Location = INLET

BOUNDARY CONDITIONS:

FLOW DIRECTION:

Option = Normal to Boundary Condition

END

FLOW REGIME:

Option = Subsonic

END

MASS AND MOMENTUM:

Mass Flow Rate = 0.0614 [kg s^-1]

Option = Mass Flow Rate

END

TURBULENCE:

Option = Medium Intensity and Eddy Viscosity Ratio

END

END

END

BOUNDARY: R1 Outlet

Boundary Type = OUTLET

Frame Type = Stationary

Location = OUTLET

BOUNDARY CONDITIONS:

FLOW REGIME:

Option = Subsonic

END

MASS AND MOMENTUM:

Option = Average Static Pressure

Relative Pressure = 6 [atm]

END

PRESSURE AVERAGING:

Option = Average Over Whole Outlet

END

END

END

BOUNDARY: R1 Shroud

Boundary Type = WALL

Frame Type = Rotating

Location = SHROUD

BOUNDARY CONDITIONS:

WALL INFLUENCE ON FLOW:

Option = No Slip

END

WALL ROUGHNESS:

Option = Rough Wall

Roughness Height = 1.6 [micron]

END

END

END

BOUNDARY: R1 to R1 Periodic 1 Side 1

Boundary Type = INTERFACE

Location = PERIODIC1

BOUNDARY CONDITIONS:

MASS AND MOMENTUM:

Option = Conservative Interface Flux

END

TURBULENCE:

Option = Conservative Interface Flux

END

END

END

BOUNDARY: R1 to R1 Periodic 1 Side 2

Boundary Type = INTERFACE

Location = PERIODIC2

BOUNDARY CONDITIONS:

MASS AND MOMENTUM:

Option = Conservative Interface Flux

END

TURBULENCE:

Option = Conservative Interface Flux

END

END

END

DOMAIN MODELS:

BUOYANCY MODEL:

Option = Non Buoyant

END

DOMAIN MOTION:

Alternate Rotation Model = true

Angular Velocity = -1450 [rev min^-1]

Option = Rotating

AXIS DEFINITION:

Option = Coordinate Axis

Rotation Axis = Coord 0.3

END

END

MESH DEFORMATION:

Option = None

END

REFERENCE PRESSURE:

Reference Pressure = 1 [atm]

END

END

FLUID MODELS:

COMBUSTION MODEL:

Option = None

END

HEAT TRANSFER MODEL:

Option = None

END

THERMAL RADIATION MODEL:

Option = None

END

TURBULENCE MODEL:

Option = k epsilon

END

TURBULENT WALL FUNCTIONS:

Option = Scalable

END

END

END

DOMAIN INTERFACE: R1 to R1 Periodic 1

Boundary List1 = R1 to R1 Periodic 1 Side 1

Boundary List2 = R1 to R1 Periodic 1 Side 2

Interface Type = Fluid Fluid

INTERFACE MODELS:

Option = Rotational Periodicity

AXIS DEFINITION:

Option = Coordinate Axis

Rotation Axis = Coord 0.3

END

END

MESH CONNECTION:

Option = Automatic

END

END

OUTPUT CONTROL:

RESULTS:

File Compression Level = Default

Option = Standard

END

END

SOLVER CONTROL:

ADVECTION SCHEME:

Option = High Resolution

END

CONVERGENCE CONTROL:

Length Scale Option = Conservative

Maximum Number of Iterations = 100

Timescale Control = Auto Timescale

Timescale Factor = 1.0

END

CONVERGENCE CRITERIA:

Residual Target = 1.E-4

Residual Type = RMS

END

DYNAMIC MODEL CONTROL:

Global Dynamic Model Control = On

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: super

Host Architecture String = intel_p4.sse2_winnt5.1

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

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 = F:/program \

files/Job/Progs/splines/e­n1300-10-1/variant2/Turbo­PrePump.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: SUPER

Job started: Mon Feb 18 16:46:23 2008

+------------------------­-------------------------­-------------------+

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

+------------------------­-------------------------­-------------------+

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

Real 40882.4 686.81 255.83 159697.0 2747.24

Integer 14035.7 235.79 87.83 54826.9 943.18

Character 2837.8 47.67 17.76 2771.3 47.67

Logical 65.0 1.09 0.41 253.9 4.37

Double 908.0 15.25 5.68 7093.8 122.03

+------------------------­-------------------------­-------------------+

| Mesh Statistics |

+------------------------­-------------------------­-------------------+

Domain Name : R1

Total Number of Nodes = 59525

Total Number of Elements = 159800

Total Number of Tetrahedrons = 77420

Total Number of Prisms = 81000

Total Number of Pyramids = 1380

Total Number of Faces = 16652

Minimum Orthogonality Angle [degrees] = 8.0 !

Maximum Aspect Ratio = 765.7 ok

Maximum Mesh Expansion Factor = 5172.0 !

Domain Interface Name : R1 to R1 Periodic 1

Non-overlap area fraction on side 1 = 1.80E-02

Non-overlap area fraction on side 2 = 7.84E-03

+------------------------­-------------------------­-------------------+

| Average Scale Information |

+------------------------­-------------------------­-------------------+

Domain Name : R1

Global Length = 1.7534E+02

Minimum Extent = 2.3761E+02

Maximum Extent = 4.9243E+02

Density = 9.9700E+02

Dynamic Viscosity = 8.8990E-04

Velocity = 3.9570E+04

Advection Time = 4.4311E-03

Reynolds Number = 7.7730E+12

+------------------------­-------------------------­-------------------+

| Checking for Isolated Fluid Regions |

+------------------------­-------------------------­-------------------+

No isolated fluid regions were found.

+------------------------­-------------------------­-------------------+

| The Equations Solved in This Calculation |

+------------------------­-------------------------­-------------------+

Subsystem : Momentum and Mass

U-Mom

V-Mom

W-Mom

P-Mass

Subsystem : TurbKE and Diss.K

K-TurbKE

E-Diss.K

CFD Solver started: Mon Feb 18 16:46:46 2008

+------------------------­-------------------------­-------------------+

| Convergence History |

+------------------------­-------------------------­-------------------+

=========================­=========================­====================

| Timescale Information |

-------------------------­-------------------------­--------------------

| Equation | Type | Timescale |

+----------------------+-­-----------------------+-­-------------------+

| U-Mom | Auto Timescale | 6.58573E-04 |

| V-Mom | Auto Timescale | 6.58573E-04 |

| W-Mom | Auto Timescale | 6.58573E-04 |

+----------------------+-­-----------------------+-­-------------------+

| K-TurbKE | Auto Timescale | 6.58573E-04 |

| E-Diss.K | Auto Timescale | 6.58573E-04 |

+----------------------+-­-----------------------+-­-------------------+

=========================­=========================­====================

OUTER LOOP ITERATION = 1 CPU SECONDS = 1.461E+01

-------------------------­-------------------------­--------------------

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

+----------------------+-­-----+---------+---------­+------------------+

| U-Mom | 0.00 | 3.8E-02 | 5.8E-01 | 2.9E-03 OK|

| V-Mom | 0.00 | 2.5E-02 | 5.6E-01 | 3.4E-03 OK|

| W-Mom | 0.00 | 9.1E-03 | 1.8E-01 | 5.8E-02 OK|

| P-Mass | 0.00 | 2.2E-03 | 1.0E-01 | 15.0 8.5E-02 OK|

+----------------------+-­-----+---------+---------­+------------------+

+------------------------­-------------------------­-------------------+

| ****** Notice ****** |

| A wall has been placed at portion(s) of an OUTLET |

| boundary condition (at 18.3% of the faces, 5.6% of the area) |

| to prevent fluid from flowing into the domain. |

| The boundary condition name is: R1 Outlet. |

| The fluid name is: Water. |

| If this situation persists, consider switching |

| to an Opening type boundary condition instead. |

+------------------------­-------------------------­-------------------+

+------------------------­-------------------------­-------------------+

| ERROR #001100279 has occurred in subroutine ErrAction. |

| Message: |

| Floating point exception: Overflow |

| |

| |

| |

| |

| |

+------------------------­-------------------------­-------------------+

+------------------------­-------------------------­-------------------+

| Writing crash recovery file |

+------------------------­-------------------------­-------------------+

+------------------------­-------------------------­-------------------+

| An error has occurred in cfx5solve: |

| |

| The ANSYS CFX solver could not be started, or exited with return |

| code 255: . 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 F:\program |

| files\Job\Progs\splines\e­n1300-10-1\variant2\Turbo­PrePump_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. Note that a lock file was left for |

| the crash file when the solver exited, and so it is probably |

| incomplete. |

+------------------------­-------------------------­-------------------+

+------------------------­-------------------------­-------------------+

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

| F:\program |

| files\Job\Progs\splines\e­n1300-10-1\variant2\Turbo­PrePump_001: |

| |

| mon |

+------------------------­-------------------------­-------------------+

This run of the ANSYS CFX Solver has finished.

[dmiantosha 11]

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

[Lektor 0]

Я конечно не шибко дольшой специалист по ANSYS, поэтому ничего не могу сказать по поводу установок модели и граничных условий, но вот эти данные меня очень сильно смущают и вызывают подозрения

+------------------------­-------------------------­-------------------+

| Average Scale Information |

+------------------------­-------------------------­-------------------+

Domain Name : R1

Global Length = 1.7534E+02

Minimum Extent = 2.3761E+02

Maximum Extent = 4.9243E+02

Density = 9.9700E+02

Dynamic Viscosity = 8.8990E-04

Velocity = 3.9570E+04

Advection Time = 4.4311E-03

Reynolds Number = 7.7730E+12

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

[a_schelyaev 366]

Да воде с ГУ все в норме. РАзве что расставил не там.

BOUNDARY: R1 Inlet

Mass Flow Rate = 0.0614 [kg s^-1]

BOUNDARY: R1 Outlet

Boundary Type = OUTLET

Option = Average Static Pressure

Relative Pressure = 6 [atm]

BOUNDARY: R1 to R1 Periodic 1 Side 1

Boundary Type = INTERFACE

Location = PERIODIC1

BOUNDARY: R1 to R1 Periodic 1 Side 2

Boundary Type = INTERFACE

Location = PERIODIC2

DOMAIN MODELS:

Option = Non Buoyant

DOMAIN MOTION:

Alternate Rotation Model = true

Angular Velocity = -1450 [rev min^-1]

Rotation Axis = Coord 0.3

Option = k epsilon

DOMAIN INTERFACE: R1 to R1 Periodic 1

Boundary List1 = R1 to R1 Periodic 1 Side 1

Boundary List2 = R1 to R1 Periodic 1 Side 2

Interface Type = Fluid Fluid

INTERFACE MODELS:

Option = Rotational Periodicity

Option = Coordinate Axis

Rotation Axis = Coord 0.3

MESH CONNECTION:

Option = Automatic

Maximum Number of Iterations = 100

Timescale Control = Auto Timescale

Timescale Factor = 1.0

Minimum Orthogonality Angle [degrees] = 8.0 !

Maximum Aspect Ratio = 765.7 ok

Maximum Mesh Expansion Factor = 5172.0 !

Domain Interface Name : R1 to R1 Periodic 1

Non-overlap area fraction on side 1 = 1.80E-02

Non-overlap area fraction on side 2 = 7.84E-03

Последние несколько строк - у тебя сопрягамеые поверхности сильно отличаются - сетка не грубовата?

И еще, попробуй начальными условями поиграть - задай скорость осевую в области примерно соответствующей некой средней.

[Dick 0]

В принципе уже все возможные проблемы описали. Еще можно проверить не попутались ли размерностями при экспорте/импорте (Global Length = 1.7534E+02). А так, как и написал Щеляев - задать вероятную среднюю скорость, либо шаг по времени уменьшить раз в 100 (хотя он и так слишком маленький и вероятно дело не в нем, а все-таки в ГУ) . Еще вместо OUTLET можно попробовать задать ГУ OPENING если физика задачи позволяет.

[pugach 8]

Спасибо всем за советы.

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

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

[pugach 8]

Проблема возникла: пытаюсь поменять вид стенок (на шероховатую) или например задать начальные условия

Domain\ initialisation\ initial conditions

после нажатия кнопеи apply вылетает сообщение:

ERROR

CCL validation failed with message:

error: essential parametr 'number of passages in 360' is missing from /turbo post data/ turbo post component: R1

R1 - вся моя расчетная область

Нашел еще ошибку. Как мне правильно подсказали, оказывается импортировал сетку не в тех единицах (опять таки ошибка в 1000)

Так, вроде поперло :) перешло на 2ую итерацию, жаль пора уходить, вечером отпишусь.

Изменено 20 февраля 2008 пользователем pugach

[pugach 8]

Еще раз добрый день.

Посоветуйте пожалуйста: в моем расчете для увеличения расчетной области были продлены основной и покрывающий диски, причем продленные участки были заданы как отдельные части (parts в icem), при постановке граничных условий они так же выставляются отдельно от Hub и Shroud, но не удается убрать их вращение, т.е. если поставить им условие wall то они будут по умолчанию rotating и это не поменять .

И еще вопросик, после 320 итераций на выходе есть завихренности 2,2% элементов и 5,6% площади, хотя дополнительная область очень большая. Может это связано с вращением дисков?