View Issue Details

IDProjectCategoryView StatusLast Update
0002957OpenFOAMBugpublic2018-05-27 11:27
Reporterpumr Assigned Tohenry  
PrioritynormalSeverityfeatureReproducibilityalways
Status closedResolutionsuspended 
PlatformGNU/LinuxOSUbuntuOS Version14.04
Product Versiondev 
Summary0002957: Coupled solver for U together with symmetry patch gives unphysical results
DescriptionWhen Coupled solver for U is enabled, then symmetry (or symmetryPlane) patches with normals having non-zero X component give unphysical values of U.

Additionally, there is an inconsistency in capital letters convention:
Type coupled requires capital "S" for smoothSolver
        type coupled;
        solver Smo
Steps To ReproduceAny case having symmetry patches not parallel with X-axis.
Additional InformationYou can test it on the attached example - windAroundBuildings tutorial which is
modified to have inlet direction parallel to Y-axis (and with type coupled for solving U).
TagsNo tags attached.

Activities

pumr

2018-05-25 07:47

reporter  

henry

2018-05-25 08:05

manager   ~0009647

> Additionally, there is an inconsistency in capital letters convention:
> Type coupled requires capital "S" for smoothSolver
> type coupled;
> solver Smo

In what way is this inconsistent? This is named consistently with the name of the class which is capitalized because it is templated on the type it block-solves. How would you rather it were named to distinguish it from the component solver?

pumr

2018-05-25 09:38

reporter   ~0009648

Now I see, sorry for this comment. Thx for the explanation...

henry

2018-05-25 14:33

manager   ~0009649

> patches with normals having non-zero X component give unphysical values of U

Could you provide a more complete description of the issue? I have run the case and while the results on the symmetry boundary are different it is not clear that the results for the coupled solver are unphysical, the velocity is parallel to the symmetry surface as expected.

henry

2018-05-25 15:02

manager   ~0009650

I have studied the case and current implementation of the vector coupled solver and it is optimized for cases with a scalar diagonal in the momentum equation, i.e. does not support directional transformations, in particular symmetry planes. The underlying templated matrix and solver implementation supports vector and tensor diagonal forms so it would be possible to add support for symmetry planes etc. but this vector coupled solver option has not proved particularly beneficial so I don't think there is a lot of point, at least not for this particular purpose.

henry

2018-05-27 11:27

manager   ~0009659

Feature request, requires funding:

https://openfoam.org/news/funding-2018/
https://openfoam.org/maintenance/

Issue History

Date Modified Username Field Change
2018-05-25 07:47 pumr New Issue
2018-05-25 07:47 pumr File Added: windAroundBuildings-Yinlet.tar.gz
2018-05-25 08:05 henry Note Added: 0009647
2018-05-25 09:38 pumr Note Added: 0009648
2018-05-25 12:26 henry Severity major => minor
2018-05-25 14:33 henry Note Added: 0009649
2018-05-25 14:58 henry Severity minor => feature
2018-05-25 15:02 henry Note Added: 0009650
2018-05-27 11:27 henry Assigned To => henry
2018-05-27 11:27 henry Status new => closed
2018-05-27 11:27 henry Resolution open => suspended
2018-05-27 11:27 henry Note Added: 0009659