Development issueshttps://develop.openfoam.com/groups/Development/-/issues2023-11-21T14:27:42Zhttps://develop.openfoam.com/Development/openfoam/-/issues/3031BUG: multiple mapFields2023-11-21T14:27:42ZKutalmış BerçinBUG: multiple mapFields### Summary
The problem is that, apart from the first 'mapFields' function object, the following 'mapFields' function objects don't appear to be registering the fields they're working on to the correct region. For instance, in the probl...### Summary
The problem is that, apart from the first 'mapFields' function object, the following 'mapFields' function objects don't appear to be registering the fields they're working on to the correct region. For instance, in the problematic case below, the 'mapFields2' function object maps the 'UMean' from 'region0' to 'coarseMesh'. However, the object registry of 'coarseMesh' lacks any 'UMean'.
### Steps to reproduce
[2274-multiple-mapFields-21Nov23-GitLab.zip](/uploads/656e6c1e94ae2143f2373bad13ea4927/2274-multiple-mapFields-21Nov23-GitLab.zip)
### Environment information
api = 2308
HEAD = 0ae6141397
compiler = clang version 15.0.7
mpi = mpirun (Open MPI) 1.10.7.0.5e373bf1fd
OS = openSUSE Leap 15.5
opts = linux64ClangDPInt32Opt
@Mattijshttps://develop.openfoam.com/Development/openfoam/-/issues/3030Dead links on documentation page2024-01-10T16:12:56ZJohan RoenbyDead links on documentation pageAll pdf links are dead on https://www.openfoam.com/documentation/overviewAll pdf links are dead on https://www.openfoam.com/documentation/overviewhttps://develop.openfoam.com/Development/openfoam/-/issues/3029Adding the repositories and trying to update apt in WSL with Ubuntu 18.04 LTS...2023-11-16T15:27:35ZMatti ViljamaaAdding the repositories and trying to update apt in WSL with Ubuntu 18.04 LTS reports 404 Not found on some URLs even though the URLs work in a browser.<!--
*** Please read this first! ***
Before opening a new issue, make sure to search for keywords in the issues
filtered by the "bug" label and check to see if it has already been reported
You can see how your report will be re...<!--
*** Please read this first! ***
Before opening a new issue, make sure to search for keywords in the issues
filtered by the "bug" label and check to see if it has already been reported
You can see how your report will be rendered on the platform by using the
"preview" tab above
-->
<!--
All text between these markers are comments and will not be present in the
report
-->
### Summary
<!-- Summarize the bug encountered concisely -->
Adding the repositories and trying to update apt in WSL with Ubuntu 18.04 LTS reports 404 Not found on some URLs even though the URLs work in a browser.
### Steps to reproduce
<!-- How one can reproduce the issue - this is very important -->
I am using WSL with Ubuntu 18.04 LTS and the guide at:
https://develop.openfoam.com/Development/openfoam/-/wikis/precompiled/windows
and in particular:
https://develop.openfoam.com/Development/openfoam/-/wikis/precompiled/debian
I am doing:
`curl https://dl.openfoam.com/add-debian-repo.sh | sudo bash`
and adding trusted=yes to:
deb [trusted=yes arch=amd64] https://dl.openfoam.com/repos/deb bionic main
in
/etc/apt/sources.list.d/openfoam.list
(This fixes another problem with the repositories. You can see the problem by removing trusted=yes.)
After this, shell reports 404 on some urls when doing the `sudo apt-get update` even though the URL (https://downloads.sourceforge.net/project/openfoam/repos/deb/dists/bionic/main/binary-amd64/Packages?ts=gAAAAABlU0ddE33NYF2BYJCrUgKdzm3tdc4pT5s50s0bO4jq_uOT3e-imqmuDbyW9FG41PrmLX9L7rvax9Dfyt7RR5uaZ2ydbg==&use_mirror=master&r=) itself works in a browser. Full output:
```
sudo apt-get update
Hit:1 http://security.ubuntu.com/ubuntu bionic-security InRelease
Hit:2 http://archive.ubuntu.com/ubuntu bionic InRelease
Hit:3 http://archive.ubuntu.com/ubuntu bionic-updates InRelease
Hit:4 http://archive.ubuntu.com/ubuntu bionic-backports InRelease
Ign:5 https://sourceforge.net/projects/openfoam/files/repos/deb bionic InRelease
Ign:6 https://sourceforge.net/projects/openfoam/files/repos/deb bionic Release
Ign:7 https://dl.openfoam.com/repos/deb bionic/main amd64 Packages
Ign:8 https://dl.openfoam.com/repos/deb bionic/main all Packages
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:7 https://dl.openfoam.com/repos/deb bionic/main amd64 Packages
Ign:8 https://dl.openfoam.com/repos/deb bionic/main all Packages
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:7 https://dl.openfoam.com/repos/deb bionic/main amd64 Packages
Ign:8 https://dl.openfoam.com/repos/deb bionic/main all Packages
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:9 https://dl.openfoam.com/repos/deb bionic/main Translation-en
Ign:10 https://dl.openfoam.com/repos/deb bionic/main all c-n-f Metadata
Ign:11 https://dl.openfoam.com/repos/deb bionic/main amd64 c-n-f Metadata
Ign:7 https://sourceforge.net/projects/openfoam/files/repos/deb bionic/main amd64 Packages
Ign:7 https://sourceforge.net/projects/openfoam/files/repos/deb bionic/main amd64 Packages
Ign:7 https://dl.openfoam.com/repos/deb bionic/main amd64 Packages
Ign:7 https://dl.openfoam.com/repos/deb bionic/main amd64 Packages
Ign:8 https://sourceforge.net/projects/openfoam/files/repos/deb bionic/main all Packages
Ign:8 https://sourceforge.net/projects/openfoam/files/repos/deb bionic/main all Packages
Ign:8 https://dl.openfoam.com/repos/deb bionic/main all Packages
Ign:8 https://dl.openfoam.com/repos/deb bionic/main all Packages
Err:7 https://sourceforge.net/projects/openfoam/files/repos/deb bionic/main amd64 Packages
404 Not Found [IP: 82.71.205.33 443]
Reading package lists... Done
E: Failed to fetch https://downloads.sourceforge.net/project/openfoam/repos/deb/dists/bionic/main/binary-amd64/Packages?ts=gAAAAABlU0ddE33NYF2BYJCrUgKdzm3tdc4pT5s50s0bO4jq_uOT3e-imqmuDbyW9FG41PrmLX9L7rvax9Dfyt7RR5uaZ2ydbg==&use_mirror=master&r= 404 Not Found [IP: 82.71.205.33 443]
E: Failed to fetch https://downloads.sourceforge.net/project/openfoam/repos/deb/dists/bionic/main/binary-all/Packages?ts=gAAAAABlU0ddOvkkOuXCeg05riwpnOabnmAYq42Yl6AoBItkbuf95la7On9UFcH_JnZb9QeAH4IobNvAs2ekXvWON2UqNm0oVQ==&use_mirror=master&r= 404 Not Found [IP: 82.71.205.33 443]
E: Some index files failed to download. They have been ignored, or old ones used instead.
```
### Example case
<!--
If possible, please create a SMALL example and attach it to your report
If you are using an older version of OpenFOAM this will also determine
whether the bug has been fixed in a more recent version
-->
### What is the current *bug* behaviour?
<!-- What actually happens -->
### What is the expected *correct* behavior?
<!-- What you should see instead -->
No errors.
### Relevant logs and/or images
<!--
Paste any relevant logs - please use code blocks (```) to format console
output, logs, and code as it's very hard to read otherwise.
-->
### Environment information
<!--
Providing details of your set-up can help us identify any issues, e.g.
OpenFOAM version : v2306|v2212|v2206|v2112|v2106 etc
Operating system : ubuntu|openSUSE|centos etc
Hardware info : any info that may help?
Compiler : gcc|intel|clang etc
-->
- OpenFOAM version : v2306
- Operating system : WSL Ubuntu 18.04 LTS
- Hardware info :
- Compiler :
### Possible fixes
<!--
If you can, link to the line of code that might be responsible for the
problem
The "/label ~bug" text is a gitlab flag that will add the "bug" label to this
issue
-->
Based on some threads like:
https://unix.stackexchange.com/questions/148303/apt-get-install-gives-404-not-found-but-url-works
it could be a bug with apt.
Or it could be a problem with the installer script.https://develop.openfoam.com/Development/openfoam/-/issues/3028non-orthogonal faces to set nonOrthoFaces2023-11-16T15:22:11ZFarshad Ansarinon-orthogonal faces to set nonOrthoFacesHello,
When I import my data.unv with the ideasUnvToFoam command, everything is ok. But after meshCheck I get this error:
*Number of severely non-orthogonal (> 70 degrees) faces: 612.
Non-orthogonality check OK.
<<Writing 612 no...Hello,
When I import my data.unv with the ideasUnvToFoam command, everything is ok. But after meshCheck I get this error:
*Number of severely non-orthogonal (> 70 degrees) faces: 612.
Non-orthogonality check OK.
<<Writing 612 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.23757 OK.
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
what should I do?
thankshttps://develop.openfoam.com/Development/openfoam/-/issues/3027Evaporation of two liquid phases using icoReactingMultiphaseInterFoam2023-11-16T15:20:51ZMark YEvaporation of two liquid phases using icoReactingMultiphaseInterFoam### Summary
The Solver I used is the icoReactingMultiphaseInterFoam V2306, the Ubuntu version is 20.04.
I try to simulate two different liquid droplets, i.e., liquid1 and liquid2, evaporating to the gas phase. If I set two liquid drop...### Summary
The Solver I used is the icoReactingMultiphaseInterFoam V2306, the Ubuntu version is 20.04.
I try to simulate two different liquid droplets, i.e., liquid1 and liquid2, evaporating to the gas phase. If I set two liquid droplets, the vapor of the liquid, i.e., vapor1, occurs in both two locations, see the picture(the same as vapor2). Physically, vapor1 should occur in location 1 and vapour2 should occur at location 2.
![figure](/uploads/1b20be65014dcb0ad7784a44b6f28dee/figure.png)
I found the code to calculate the massSpeciesTransfer in multiphaseInter::MassTransferPhaseSystem.C. No matter what species is transferred from the liquid phase (liquid1 or liquid2), the value of Su is the same.
The attached is the setting case.[debug.tar.gz](/uploads/83421f905cc567ec6fd5249dac37fee9/debug.tar.gz)https://develop.openfoam.com/Development/openfoam/-/issues/3026rigidBodyMotion linearAxialAngularSpring theta angle appears to be incorrectl...2023-11-16T15:20:33ZChristian RohrrigidBodyMotion linearAxialAngularSpring theta angle appears to be incorrectly signed - spring acts in wrong direction### Summary
Theta angle used in (at least) `linearAxialAngularSpring.C` in the rigidBodyMotion solver appears to be signed incorrectly.
### Steps to reproduce
Apply the linearAxialAngularSpring restraint to a body, and observe it be p...### Summary
Theta angle used in (at least) `linearAxialAngularSpring.C` in the rigidBodyMotion solver appears to be signed incorrectly.
### Steps to reproduce
Apply the linearAxialAngularSpring restraint to a body, and observe it be pulled in the wrong direction. Have tested with pimpleDyMFoam. This has occured on another case which I cannot share, so I have created a new one from one of the tutorials.
Modifying the code for this restraint to ouptut the current angular position shows that the "theta" variable is signed incorrectly. For example a clockwise rotation around the z axis of 45° will be reported as 45° instead of -45°.
### Example case
Attached a case modified from the overset airfoil simpleFoam tutorial. It has been made transient, switched to rigidBody, and a spring constraint has been applied. Increasing the stiffness of this spring exacerbates the movement of the airfoil rather than assisting.
[testCase_axialSpring.zip](/uploads/676fc081bc3e63b66c9376639aa163e3/testCase_axialSpring.zip)
### What is the current _bug_ behaviour?
Spring moment acts in wrong direction.
### What is the expected _correct_ behavior?
A positive rotation should create a negative torque.
### Relevant logs and/or images
In the example case, I have attached a spring at the bottom of the domain to the airfoil with an Rz joint. The spring pulls the airfoil down and to the right, and this is accelerated for higher spring stiffnesses. Some movement in that direction is expected for the aerodynamic load, but it should reduce with increasing stiffness rather than worsen.
The case should be run with a stiffness of 1 and then with a stiffness of 100 to observe the issue.
Alternatively, add an `Info <<` message reporting the theta angle of the airfoil body at line 143 of `linearAxialAngularSpring.C` : it is signed incorrectly for the observed rotations.
![image.png](/uploads/71c2ddc6030e147d4b073b8df15267aa/image.png){width="314" height="299"}
### Environment information
- OpenFOAM version : v2112
- Operating system : WSL - OpenSUSE Leap 15.4
- Hardware info : x86_64
- Compiler : gcc
### Possible fixes
Manual flip to theta's sign on line `112` of `linearAxialAngularSpring.C`. But I have not been able to get to the bottom of why this is incorrectly signed to begin with. Could this just be due to case configuration?https://develop.openfoam.com/Development/openfoam/-/issues/3025stringOps::split with keep empty accidentally ignores non-empty trailing element2023-11-09T10:57:48ZMark OLESENstringOps::split with keep empty accidentally ignores non-empty trailing elementMark OLESENMark OLESENhttps://develop.openfoam.com/Development/openfoam/-/issues/3024compressibleInterFoam crahes when compiled with gcc 11.4.02024-01-19T20:23:54ZChris SesslercompressibleInterFoam crahes when compiled with gcc 11.4.0<!--
*** Please read this first! ***
Before opening a new issue, make sure to search for keywords in the issues
filtered by the "bug" label and check to see if it has already been reported
You can see how your report will be re...<!--
*** Please read this first! ***
Before opening a new issue, make sure to search for keywords in the issues
filtered by the "bug" label and check to see if it has already been reported
You can see how your report will be rendered on the platform by using the
"preview" tab above
-->
<!--
All text between these markers are comments and will not be present in the
report
-->
### Summary
<!-- Summarize the bug encountered concisely -->
The provided case crashes within the first 5 iterations due to "Negative initial temperature T0" when compiled with gcc 11.4.0. In contrast, when compiled with gcc 9.5.0 this problem doesn't occur.
### Steps to reproduce
<!-- How one can reproduce the issue - this is very important -->
Run the case after initialization with blockMesh and setFields
### Example case
<!--
If possible, please create a SMALL example and attach it to your report
If you are using an older version of OpenFOAM this will also determine
whether the bug has been fixed in a more recent version
-->
[compressibleInterFoam.zip](/uploads/f52c0f14bdf8eb473e1b5916e70af854/compressibleInterFoam.zip)
### What is the current *bug* behaviour?
<!-- What actually happens -->
This problem only arised with 3D cases, a different 2D case ran perfectly
### What is the expected *correct* behavior?
<!-- What you should see instead -->
Compilation with gcc 9.5.0 runs stable
### Relevant logs and/or images
<!--
Paste any relevant logs - please use code blocks (```) to format console
output, logs, and code as it's very hard to read otherwise.
-->
<details><summary>Error log</summary>
[1] [0]
[0]
[0] --> FOAM FATAL ERROR: (openfoam-2012)
[0] Negative initial temperature T0: -180.035335774
[0]
[0] From Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
[0] in file /home/itvcs/OpenFOAM/OpenFOAM-v2012/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 57.
</details>
### Environment information
<!--
Providing details of your set-up can help us identify any issues, e.g.
OpenFOAM version : v2306|v2212|v2206|v2112|v2106 etc
Operating system : ubuntu|openSUSE|centos etc
Hardware info : any info that may help?
Compiler : gcc|intel|clang etc
-->
- OpenFOAM version : v2012 commit 79e353b84e3e7cace6ea35c6b7570dfd198d0135
- Operating system : Ubuntu 22.04.3 LTS
- Hardware info : Intel 13th gen
- Compiler : gcc 9.5.0 and 11.4.0
### Possible fixes
<!--
If you can, link to the line of code that might be responsible for the
problem
The "/label ~bug" text is a gitlab flag that will add the "bug" label to this
issue
-->https://develop.openfoam.com/Development/openfoam/-/issues/3023object registry filtered objects uses derived name2023-11-07T21:01:55ZMark OLESENobject registry filtered objects uses derived nameAffects the implementation of csorted etc.
It checks the name from the cast pointer
```
matchName(ptr->name())
```
instead of from the object pointer
```
matchName(obj->name())
```
This only seems to affect finite-area directly, since t...Affects the implementation of csorted etc.
It checks the name from the cast pointer
```
matchName(ptr->name())
```
instead of from the object pointer
```
matchName(obj->name())
```
This only seems to affect finite-area directly, since the cast pointer will have "region0" as its name, whereas the object is actually registered as "faMesh". Cross-ref EP2297
Regression introduced with 95e2a2e887d but was not exposed until 129b738136f (develop branch)Mark OLESENMark OLESENhttps://develop.openfoam.com/Development/openfoam/-/issues/3022parallel snappyHexMesh errors2023-11-06T21:38:33ZMadis Listakparallel snappyHexMesh errors### Summary
Module: OpenFabrics (openib) problem, can't run snappyHexMesh in parallel
### Steps to reproduce
I am running openFoam2306 OpenSuse15.5 rpm packages and "tutorials/mesh/snappyHexMesh/gap_detection/" example
and I get an er...### Summary
Module: OpenFabrics (openib) problem, can't run snappyHexMesh in parallel
### Steps to reproduce
I am running openFoam2306 OpenSuse15.5 rpm packages and "tutorials/mesh/snappyHexMesh/gap_detection/" example
and I get an error when I try to execute Allrun and do snappyHexMesh in parallel
[[53808,1],7]: A high-performance Open MPI point-to-point messaging module
was unable to find any relevant network interfaces:
Module: OpenFabrics (openib)
Host: (my computer name)
Caught signal 4 (Illegal instruction: illegal operand)
==== backtrace (tid: 16643) ====
0 0x0000000000016910 __funlockfile() ???:0
1 0x000000000082eb52 fi_getinfo() ???:0
2 0x0000000000828d17 fi_getinfo() ???:0
3 0x0000000000839707 fi_getinfo() ???:0
4 0x0000000000817833 fi_getinfo() ???:0
5 0x000000000005369b fi_getinfo() ???:0
6 0x00000000000048bb ompi_mtl_ofi_progress_no_inline() ???:0
7 0x00000000000a5190 ompi_mtl_base_select() ???:0
8 0x000000000000541a mca_pml_cm_cancel() ???:0
9 0x00000000000ad645 mca_pml_base_select() ???:0
10 0x000000000004c797 ompi_mpi_init() ???:0
11 0x000000000006cbc8 PMPI_Init_thread() ???:0
12 0x000000000000dc42 Foam::UPstream::init() ???:0
13 0x0000000000413feb Foam::argList::argList() ???:0
14 0x000000000001955d main() ???:0
15 0x000000000003524d __libc_start_main() ???:0
16 0x000000000001f9ca _start() /home/abuild/rpmbuild/BUILD/glibc-2.31/csu/../sysdeps/x86_64/start.S:120
--------------------------------------------------------------------------
mpirun noticed that process rank 7 with PID 0 on node (my computer name) exited on signal 4 (Illegal instruction).
and snappyHexMesh is failing
Probably something is not installed or is not compiled correctly.https://develop.openfoam.com/Development/openfoam/-/issues/3021Error in application of parallel axis theorem in sixDoFRigidBodyMotion class2023-11-07T09:32:33ZJohan RoenbyError in application of parallel axis theorem in sixDoFRigidBodyMotion classThere is a bug in the application of the parallel axis theorem in the constructor here:
https://develop.openfoam.com/Development/openfoam/-/blob/master/src/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion.C?ref_type=hea...There is a bug in the application of the parallel axis theorem in the constructor here:
https://develop.openfoam.com/Development/openfoam/-/blob/master/src/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion.C?ref_type=heads#L148
where we do `momentOfInertia_ += mass_*diag(I*magSqr(R) - sqr(R));`
The matrix `sqr(R)` is generally not diagonal and by taking diag() of it, we omit its off-diagonal terms, which is an error.
Thus, simulations where the centre of rotation and centre of mass do not coincide must be expected to give wrong results.
The correct behaviour in such simulations should be tested and verified with a simple test case.Kutalmış BerçinKutalmış Berçinhttps://develop.openfoam.com/Development/openfoam/-/issues/3020momentOfInertia tensor should be symmTensor, not diagTensor in sixDofRigidBod...2023-11-06T20:44:24ZJohan RoenbymomentOfInertia tensor should be symmTensor, not diagTensor in sixDofRigidBodyMotion classThe moment of inertia is a symmetric tensor, not necessarily diagonal.
This should be reflected in its type in the sixDoFRigidBodyMotion class:
https://develop.openfoam.com/Development/openfoam/-/blob/master/src/sixDoFRigidBodyMotion/s...The moment of inertia is a symmetric tensor, not necessarily diagonal.
This should be reflected in its type in the sixDoFRigidBodyMotion class:
https://develop.openfoam.com/Development/openfoam/-/blob/master/src/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion.H#L111
and in how it is specified in the dynamicFvMeshDict in cases using it.
Indeed, for a symmetric matrix (with all eigenvalues non-zero) one can always find an orthogonal basis in which it is diagonal. This basis, however, often does not coincide with the natural coordinate system of a given simulation case.Kutalmış BerçinKutalmış Berçinhttps://develop.openfoam.com/Development/openfoam/-/issues/3019decompositionMethod calcCellCells with weighting looks odd2023-11-07T22:12:39ZMark OLESENdecompositionMethod calcCellCells with weighting looks oddIn calcCellCells, the weighting per cell is done by using the face area of some attached face. Eg,
// For internal faces is just offsetted owner and neighbour
for (label facei = 0; facei < mesh.nInternalFaces(); ++facei)
{
...In calcCellCells, the weighting per cell is done by using the face area of some attached face. Eg,
// For internal faces is just offsetted owner and neighbour
for (label facei = 0; facei < mesh.nInternalFaces(); ++facei)
{
const label own = agglom[faceOwner[facei]];
const label nei = agglom[faceNeighbour[facei]];
const label ownIndex = offsets[own] + nFacesPerCell[own]++;
const label neiIndex = offsets[nei] + nFacesPerCell[nei]++;
m[ownIndex] = globalAgglom.toGlobal(nei);
w[ownIndex] = mag(mesh.faceAreas()[facei]);
m[neiIndex] = globalAgglom.toGlobal(own);
w[ownIndex] = mag(mesh.faceAreas()[facei]);
}
The fact that `w[ownIndex]` is specified twice (not `w[neiIndex]`) looks weird, but seems to have been in there forever.
@MattijsMark OLESENMark OLESENhttps://develop.openfoam.com/Development/openfoam/-/issues/3018pTraits for complex looks incorrect2023-11-13T13:53:25ZMark OLESENpTraits for complex looks incorrectThe complex type is defined to have two components (ie, real, imag) but has cmptType as "complex" - doesn't really make sense, and has labelType as label anyhow.
It doesn't seem to be used anywhere but could be interesting with different...The complex type is defined to have two components (ie, real, imag) but has cmptType as "complex" - doesn't really make sense, and has labelType as label anyhow.
It doesn't seem to be used anywhere but could be interesting with different backends (eg ADIOS).
@andy @kutiv2406Kutalmış BerçinKutalmış Berçinhttps://develop.openfoam.com/Development/openfoam/-/issues/3017Bug liquidFilmStepWithSprinklers2023-11-06T21:29:20ZAlberto ceschinBug liquidFilmStepWithSprinklers### Summary
tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinklers breaks after about 1.9s . endTime in controlDict is 1, but particles are injected at 5 and it is suggested to increase endTime to 10. Same error in serial ...### Summary
tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinklers breaks after about 1.9s . endTime in controlDict is 1, but particles are injected at 5 and it is suggested to increase endTime to 10. Same error in serial and parallel. Tested on 2 different machines.
### Steps to reproduce
Increase endTime from 1 to 5s;`./Allrun`; At 1.9375s it will stop. Error below.
### Relevant logs and/or images
```plaintext
Solving3-D cloud reactingCloud1
Cloud: reactingCloud1
Current number of parcels = 0
Current mass in system = 0
Linear momentum = (0 0 0)
|Linear momentum| = 0
Linear kinetic energy = 0
Average particle per parcel = 0
Injector model1:
- parcels added = 0
- mass introduced = 0
Surface film:
- parcels absorbed = 0
- mass absorbed = 0
- parcels ejected = 0
- new splash parcels = 0
Parcel fate: system (number, mass)
- escape = 0, 0
Parcel fate: patch (number, mass) base
- escape = 0, 0
- stick = 0, 0
Parcel fate: patch (number, mass) outlet
- escape = 0, 0
- stick = 0, 0
Parcel fate: patch (number, mass) inlet
- escape = 0, 0
- stick = 0, 0
Parcel fate: patch (number, mass) sides
- escape = 0, 0
- stick = 0, 0
Parcel fate: patch (number, mass) top
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 0, 0
Mass transfer phase change = 0
Mass transfer devolatilisation = 0
Mass transfer surface reaction = 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
Evolving kinematicThinFilm for region region0
DILUPBiCGStab: Solving for Uf_filmx, Initial residual = 0.0118513, Final residual = 9.60929e-15, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmy, Initial residual = 0.0109228, Final residual = 6.0994e-11, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmz, Initial residual = 0.000543556, Final residual = 8.40472e-10, No Iterations 1
DILUPBiCGStab: Solving for hf_film, Initial residual = 0.00211817, Final residual = 4.50239e-12, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmx, Initial residual = 0.0122752, Final residual = 4.82225e-15, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmy, Initial residual = 0.00786743, Final residual = 4.88628e-11, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmz, Initial residual = 0.000531824, Final residual = 8.05986e-10, No Iterations 1
DILUPBiCGStab: Solving for hf_film, Initial residual = 0.000296377, Final residual = 5.39824e-13, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmx, Initial residual = 0.00201142, Final residual = 4.85637e-15, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmy, Initial residual = 0.000200432, Final residual = 3.70026e-13, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmz, Initial residual = 0.000302661, Final residual = 3.90512e-10, No Iterations 1
DILUPBiCGStab: Solving for hf_film, Initial residual = 5.02967e-05, Final residual = 1.05152e-13, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmx, Initial residual = 0.000488423, Final residual = 2.39743e-15, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmy, Initial residual = 9.81888e-05, Final residual = 2.3487e-13, No Iterations 1
DILUPBiCGStab: Solving for Uf_filmz, Initial residual = 0.000156307, Final residual = 1.99035e-10, No Iterations 1
DILUPBiCGStab: Solving for hf_film, Initial residual = 1.40484e-05, Final residual = 6.25537e-14, No Iterations 1
Film h min/max = (4.55048e-06 0.00165964)
Film mag(U) min/max = (0 2.51884)
smoothSolver: Solving for Ux, Initial residual = 0.00197335, Final residual = 2.55965e-07, No Iterations 5
smoothSolver: Solving for Uy, Initial residual = 0.00928135, Final residual = 8.40566e-07, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 0.00411528, Final residual = 5.61892e-07, No Iterations 5
DILUPBiCGStab: Solving for h, Initial residual = 0.00497396, Final residual = 1.94128e-07, No Iterations 1
T gas min/max = 299.53, 300.008
DICPBiCGStab: Solving for p_rgh, Initial residual = 0.0032935, Final residual = 0.000297896, No Iterations 11
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.43459e-07, global = -3.59919e-08, cumulative = -0.000505954
DICPBiCGStab: Solving for p_rgh, Initial residual = 0.00119447, Final residual = 0.000117951, No Iterations 10
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.54758e-07, global = -3.68368e-08, cumulative = -0.000505991
DICPBiCGStab: Solving for p_rgh, Initial residual = 0.000190021, Final residual = 8.92924e-07, No Iterations 69
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.92866e-09, global = 3.00806e-10, cumulative = -0.000505991
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.00229694, Final residual = 4.45827e-07, No Iterations 1
DILUPBiCGStab: Solving for k, Initial residual = 0.00441566, Final residual = 7.895e-07, No Iterations 1
ExecutionTime = 57.35 s ClockTime = 57 s
Courant Number mean: 0.466217 max: 0.868391
Max film Courant Number: 2.49983 Film velocity magnitude: 4.99966
deltaT = 0.0125
Time = 1.9375
Solving3-D cloud reactingCloud1
[stack trace]
=============
#1 Foam::sigSegv::sigHandler(int) in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::ThermoSurfaceFilm<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > > > > >::setParcelProperties(Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle> > > >&, int) const in ~/OpenFOAM/ceschia-v2306/platforms/linux64GccDPInt32Opt/lib/liblagrangianIntermediate.so
#4 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#5 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#6 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#7 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#8 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#9 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
#10 ? in /lib/x86_64-linux-gnu/libc.so.6
#11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12 ? in ~/OpenFOAM/OpenFOAM-v2306/platforms/linux64GccDPInt32Opt/bin/reactingParcelFoam
=============
Segmentation fault (core dumped)
```
### Environment information
- OpenFOAM version : 2306
- Operating system : ubuntu 22 and 20Kutalmış BerçinKutalmış Berçinhttps://develop.openfoam.com/Development/openfoam/-/issues/3016Paraview, OpenFoam Error.2023-11-09T10:59:00ZMathew DevaParaview, OpenFoam Error.I installed OpenFoam 2306 and paraview 5.11 from official sites. initaly i got the error as :
Invalid $PV_PLUGIN_PATH and -plugin-path= not defined
No supplementary ParaView/OpenFOAM reader modules
Using builtin reader: paraFoam -vtk
C...I installed OpenFoam 2306 and paraview 5.11 from official sites. initaly i got the error as :
Invalid $PV_PLUGIN_PATH and -plugin-path= not defined
No supplementary ParaView/OpenFOAM reader modules
Using builtin reader: paraFoam -vtk
Created temporary 'pitzDaily.foam'
qt.qpa.plugin: Could not load the Qt platform plugin "xcb" in "" even though it was found.
This application failed to start because no Qt platform plugin could be initialized. Reinstalling the application may fix this problem.
Available platform plugins are: xcb.
But after i executed sudo apt-get install libqt5x11extras5, now i am getting this error:
Invalid $PV_PLUGIN_PATH and -plugin-path= not defined
No supplementary ParaView/OpenFOAM reader modules
Using builtin reader: paraFoam -vtk
Created temporary 'pitzDaily.foam'
Help me with this. I'm a newbie to OpenFoam.https://develop.openfoam.com/Development/openfoam/-/issues/3015label 64 compilation failing2023-11-06T12:03:54ZPawan Ghildiyallabel 64 compilation failingHello,
While compiling the branch update-redistributePar with 64 label size, I encountered following error messagess
`M_DP -DWM_LABEL_SIZE=64 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof...Hello,
While compiling the branch update-redistributePar with 64 label size, I encountered following error messagess
`M_DP -DWM_LABEL_SIZE=64 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -Wno-unknown-pragmas -O3 -DNoRepository -ftemplate-depth-100 -I/home/pawan/OpenFOAM/OpenFOAM-v23.06-23W45-Redhat7.8/OpenFOAM-v23.06-23W45/build/linux64Gcc85DPInt64Opt/src/OpenFOAM -DHAVE_LIBZ -iquote. -IlnInclude -I/home/pawan/OpenFOAM/OpenFOAM-v23.06-23W45-Redhat7.8/OpenFOAM-v23.06-23W45/src/OpenFOAM/lnInclude -I/home/pawan/OpenFOAM/OpenFOAM-v23.06-23W45-Redhat7.8/OpenFOAM-v23.06-23W45/src/OSspecific/POSIX/lnInclude -fPIC -c primitives/chars/lists/charList.C -o /home/pawan/OpenFOAM/OpenFOAM-v23.06-23W45-Redhat7.8/OpenFOAM-v23.06-23W45/build/linux64Gcc85DPInt64Opt/src/OpenFOAM/primitives/chars/lists/charList.o
global/fileOperations/fileOperation/fileOperation.C:83:47: error: conflicting declaration 'Foam::label Foam::fileOperation::nProcsFilter_'
Foam::label Foam::fileOperation::nProcsFilter_(-1);
^
In file included from global/fileOperations/fileOperation/fileOperation.C:29:
global/fileOperations/fileOperation/fileOperation.H:158:20: note: previous declaration as 'int Foam::fileOperation::nProcsFilter_'
static int nProcsFilter_;
^~~~~~~~~~~~~`v2306Mattijs Janssens4-Mattijs@users.noreply.develop.openfoam.comMattijs Janssens4-Mattijs@users.noreply.develop.openfoam.comhttps://develop.openfoam.com/Development/openfoam/-/issues/3014ENH: patchSeed method to have more 'uniform' distribution2023-12-08T11:42:58ZPrashant SonakarENH: patchSeed method to have more 'uniform' distributioncross ref: EP#1875
- Maybe allow area based weighting to have much equal bias to generate a seed point ?cross ref: EP#1875
- Maybe allow area based weighting to have much equal bias to generate a seed point ?Andrew HeatherAndrew Heatherhttps://develop.openfoam.com/Development/openfoam/-/issues/3013ENH: mapFields FO2024-01-03T10:24:07ZPrashant SonakarENH: mapFields FOCross ref : EP#2230
- When using consistent = false, allow patchMap to be created automatically based on boundaries from target and source
- Probably filter out the patches listed in cuttingPatches entry
- any specific treatment - warn/...Cross ref : EP#2230
- When using consistent = false, allow patchMap to be created automatically based on boundaries from target and source
- Probably filter out the patches listed in cuttingPatches entry
- any specific treatment - warn/error on patches which are present in target and not in source and not in cuttingPatches ?
Maybe add-on tickets to consider #2984 #1939 #2333https://develop.openfoam.com/Development/openfoam/-/issues/3012variableHeightFlowRateInletVelocity & dynamic meshes2023-12-06T14:53:25ZShannon LeakeyvariableHeightFlowRateInletVelocity & dynamic meshes### Summary
A colleague (Bernardas Jankauskas) and I have been investigating a problem where `variableHeightFlowRateInletVelocity` doesn't seem to work with dynamic meshes. The problem has also been mentioned on [CFD Online](https://www....### Summary
A colleague (Bernardas Jankauskas) and I have been investigating a problem where `variableHeightFlowRateInletVelocity` doesn't seem to work with dynamic meshes. The problem has also been mentioned on [CFD Online](https://www.cfd-online.com/Forums/openfoam/220909-variableheightflowrateinletvelocity-whith-dynamicmesh.html) by another user.
### Steps to reproduce
1. Copy the `weirOverflow` tutorial for `interFoam`
2. Add a `constant/dynamicMeshDict` file
3. Add a `0.orig/pointDisplacement` file
4. Amend `system/fvSolution` to include a solver for `cellDisplacement`
5. Run without making any changes to the boundary condition (it crashes)
6. Run after changing `value` in the `inlet` boundary condition in `0.orig/alpha.water` to `0.000001` (it doesn't crash)
### Example case
[weirOverflow.zip](/uploads/cb578ffdd40af3da21be22650d40d716/weirOverflow.zip)
### What is the current *bug* behaviour?
Divide by zero [on this line](https://develop.openfoam.com/Development/openfoam/-/blob/master/src/finiteVolume/fields/fvPatchFields/derived/variableHeightFlowRateInletVelocity/variableHeightFlowRateInletVelocityFvPatchVectorField.C#L124)
### What is the expected *correct* behavior?
It shouldn't divide by zero
### Relevant logs
```
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib64/libpthread.so.0
#3 Foam::variableHeightFlowRateInletVelocityFvPatchVectorField::updateCoeffs() at ??:?
#4 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() at ??:?
#5 Foam::dynamicMotionSolverFvMesh::update() at ??:?
#6 ? at ??:?
#7 __libc_start_main in /lib64/libc.so.6
#8 ? at ??:?
```
### Environment information
- OpenFOAM version: v2112
### Possible fixes
`(alphap + VSMALL)` but that wouldn't explain why it thinks `alphap` is zero only when a dynamic mesh is used?