View Issue Details

IDProjectCategoryView StatusLast Update
0001217OpenFOAMBugpublic2014-03-28 20:33
Reporteruser885Assigned Touser2 
PriorityhighSeveritymajorReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSUbuntuOS Version12.x
Summary0001217: problem with cyclicAMI
Descriptionrotation of cyclicAMI does not work for complex geometries
Steps To Reproducetry a complex axial turbomachinery geometry
Additional Informationit can be fixed changing the code in:
./meshTools/lnInclude/cyclicAMIPolyPatch.C
where a subtraction is needed instead of a sum (see fix#1 and fix#2) in the following



 [blablabla.....]

 case ROTATIONAL:
        {
            tensor revT = tensor::zero;

            if (rotationAngleDefined_)
            {
                tensor T(rotationAxis_*rotationAxis_);

                tensor S
                (
                    0, -rotationAxis_.z(), rotationAxis_.y(),
                    rotationAxis_.z(), 0, -rotationAxis_.x(),
                    -rotationAxis_.y(), rotationAxis_.x(), 0
                );

                tensor RPos
                (
                    T
// originalcode#1 + cos(rotationAngle_)*(tensor::I + T)
fix#1 + cos(rotationAngle_)*(tensor::I - T)
                  + sin(rotationAngle_)*S
                );
                tensor RNeg
                (
                    T
//originalcode#2 + cos(-rotationAngle_)*(tensor::I + T)
fix #2 + cos(-rotationAngle_)*(tensor::I - T)
                  + sin(-rotationAngle_)*S
                );

                // check - assume correct angle when difference in face areas
                // is the smallest
                vector transformedAreaPos = sum(half0Areas & RPos);
                vector transformedAreaNeg = sum(half0Areas & RNeg);
                vector area1 = sum(half1Areas);
                reduce(transformedAreaPos, sumOp<vector>());
                reduce(transformedAreaNeg, sumOp<vector>());
                reduce(area1, sumOp<vector>());

[....blablabla....]
TagsNo tags attached.

Relationships

has duplicate 0001249 resolveduser2 Transformation Tensor in AMI is calculated wrong when using given rotation angle 

Activities

user2

2014-03-28 16:01

  ~0002986

Thanks for the report - fixed by commit b99b1d2

Issue History

Date Modified Username Field Change
2014-03-12 15:49 user885 New Issue
2014-03-28 15:43 user2 Relationship added has duplicate 0001249
2014-03-28 16:01 user2 Note Added: 0002986
2014-03-28 16:01 user2 Status new => resolved
2014-03-28 16:01 user2 Fixed in Version => 2.3.x
2014-03-28 16:01 user2 Resolution open => fixed
2014-03-28 16:01 user2 Assigned To => user2