Integration of new turbulence model into DAFOAM and encountering errors

[lnhang0618]

Professor He,

I have been working on integrating a complex turbulence model, kv2OmegaAr, into DAFOAM. The relevant files have been placed in the models and DAModel directories. My branch with the updates can be found at https://github.com/lnhang0618/dafoam.

I attempted to validate the integration by running the PH example, which I have also uploaded to the repository. However, I have encountered the following error during testing:

[[61757,1],0]: A high-performance Open MPI point-to-point messaging module
was unable to find any relevant network interfaces:

Module: OpenFabrics (openib)
  Host: 6f0950b74969

Another transport will be used instead, although this may result in
lower performance.

NOTE: You can disable this warning by setting the MCA parameter
btl_base_warn_component_unused to 0.
--------------------------------------------------------------------------

[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see https://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run 
[0]PETSC ERROR: to get more information on the crash.
[0]PETSC ERROR: Run with -malloc_debug to check if memory corruption is causing the crash.
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 50176059.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

I've searched through past discussions but have not found a lead that could help resolve these errors.

I would greatly appreciate your insights or any suggestions you might have regarding this issue.

Best regards,
Luo Yinhang

[lnhang0618]

kv2Omega_Ar.zip

[lnhang0618]

I have fixed the problem before，which is caused by 'v2s' definition.but I meet a new one：

--> FOAM FATAL ERROR: 

--> FOAM FATAL ERROR: 

    request for volScalarField::Internal kv2OmegaArFieldInversion:G from objectRegistry region0 failed
    available objects of type volScalarField::Internal are

50
(
alphaT
v2Res
USingleComponentData
v2
nut
yWall
((fSS*min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0))*v2)
v2Ref
v2ResRef
betaFieldInversion
((CRNAT*(1-exp((-phiNAT|ANAT))))*(1.41421*mag(skew(grad(U)))))
ArResRef
lambdaT
ReW
k
lambdaEff
((((fv*fINT)*Cmu)*sqrt(((fSS*min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0))*v2)))*lambdaEff)
surfaceFriction
Ar
omegaRef
((Cw2*sqr(min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0)))*omega)
ArRes
nu
kRef
(1-exp((((-CNC*sqrt(k))*yWall)|nu)))
(2*magSqr(symm(grad(U))))
omegaRes
kResRef
pData
omegaResPartDeriv
omegaResRef
(((((betaStar*2)*(1-(1-((1-F1)*fSS))))*sigmaW2)|max(omega,omegaMin))*(grad(k)&grad(omega)))
alphaPorosity
pRes
fv
kRes
kResPartDeriv
p
(((CR*(1-exp((-phiBP|ABP))))*omega)|max(min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0),ROTVSMALL))
v2ResPartDeriv
((sqr(min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0))*v2)|(nu*max(omega,omegaMin)))
v2l
(1.41421*mag(skew(grad(U))))
ArRef
omega
surfaceFrictionData
(1-((1-F1)*fSS))
min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0)
min((((((((C11*fTaul)*(1.41421*mag(skew(grad(U)))))*sqr(lambdaEff))*sqrt(v2l))*lambdaEff)|nu)+((((C12*betaTS)*pow4((lambdaEff|Clambda)))*sqr((1.41421*mag(skew(grad(U))))))|nu)),((0.5*(k-((fSS*min(pow((lambdaEff|max(lambdaT,SMALL)),0.666667),1.0))*v2)))|max(sqrt((2*magSqr(symm(grad(U))))),omegaMin)))
ArResPartDeriv
)

[lnhang0618]

In order to fix this problem，change the type of bottomWall in omega from omegawallfunction to zerogradient