Added a very nice d3q27 viscoplastic model. Two example xml-scripts a… #102
Conversation
…re available for tests
| AddNodeType("CoutVelocity","BOUNDARY") | ||
| AddNodeType("PCoutVelocity","BOUNDARY") | ||
| AddNodeType("WVelocityEq","BOUNDARY") | ||
| AddNodeType("WVelocityBB","BOUNDARY") |
There was a problem hiding this comment.
Change BC names.
[N/S/W/E/I/O]Velocity_CleverVariant would be proper way
If you use In/Out (z direction) BC's, check naming with newest cumulant form Dmitro
There was a problem hiding this comment.
actually, there are no options for "WVelocityEq" and "WVelocityBB" -will be removed. Naming convention will be followed
| AddSetting(name="ForceX",default=0, comment='Force force X') | ||
| AddSetting(name="ForceY",default=0, comment='Force force Y') | ||
| AddSetting(name="ForceZ",default=0, comment='Force force Z') | ||
| AddSetting(name="SigmaY",default=0, comment='Yield stress') |
| AddDensity(name="nu_app", dx=0, dy=0, dz=0) | ||
| AddDensity(name="yield_stat", dx=0, dy=0, dz=0) | ||
|
|
||
| AddSetting(name="GalileanCorrection",default=1.,comment='Galilean correction term') |
There was a problem hiding this comment.
Unused. Will be removed.
| AddQuantity( name="P",unit="Pa") | ||
| AddQuantity( name="U",unit="m/s",vector=T) | ||
| AddQuantity( name="nu_app",unit="m2/s") | ||
| AddQuantity( name="yield_stat") |
There was a problem hiding this comment.
Is it possible to speficy "empty" or "-" unit? yield_stat shows if the region is yielded or not (value 0 or 1)
There was a problem hiding this comment.
Then don't specify any, its OK. And maybe add comment to future generations :]
| AddGlobal(name="YZrho2", comment='Volume flux', unit="kg/m") | ||
| AddGlobal(name="YZarea", comment='Volume flux', unit="m2") | ||
|
|
||
| #Adding terms for supporting time-correlation for synthetic turbulence |
| AddGlobal(name="YZvz", comment='Volume flux', unit="m3/s") | ||
| AddGlobal(name="YZrho1", comment='Volume flux', unit="kg/m") | ||
| AddGlobal(name="YZrho2", comment='Volume flux', unit="kg/m") | ||
| AddGlobal(name="YZarea", comment='Volume flux', unit="m2") |
There was a problem hiding this comment.
Do you use/intend to implement all those globals and nodetypes?
There was a problem hiding this comment.
Yes, I do. Very often when I run simulations I have to integrate fluxes along principal directions
1. naming convention as suggested by mdzik 2. removed SymmetryBottom and symmetryTop - thay act as SymmetryY - no need to distinguish that 3. units added at some quantities 4. two test ".xml" files modified adequately to follow the namig convention
| AddNodeType("NPressure_ZouHe", "BOUNDARY") | ||
| AddNodeType("SPressure_ZouHe", "BOUNDARY") | ||
| AddNodeType("EPressure_ZouHe", "BOUNDARY") | ||
| AddNodeType("WPressure_ZouHe", "BOUNDARY") |
There was a problem hiding this comment.
OK, my fault:
ZouHe is "default" way to handle BC's, so no need for suffix. Use suffix only if you have different "variants" of bc, like:
EVelocity (standard one)
EVelocity_ViscosityCorrection
EVelocity_ConstantFlux
EVelocity_Kowalski
or something like that
once again, sorry for double work. (I was reffering rather to PCoutVelocity)
|
I'm merging, as I think the model is more or less ready - and I'm planning a code cleanup soon. |
added d3q27_viscoplastic model. Please test it, two test scripts are available. Enjoy!