[generalKEE]

When I used my own written boundary conditions for calculations, I had a strange error. Could someone help me?

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
                                                                                                 
                        
Create mesh for time = 0
                                                                                                 
                        
                                                                         
                        
PISO: Operating solver in PISO mode
                                                                                                 
                        
Reading transportProperties
                                                                                                 
                        
Reading field p
                                                                                                 
                         

Reading field U
                                                                                                  
                        
Reading field C
                                                                                                 
                        
Reading/calculating face flux field phi
                                                                                                 
                        
                                                                         
                        
Starting time loop
                                                                                                 
                        
Time = 5e-05
                                                                                                 
                        
Courant Number mean: 7.5e-05 max: 0.194817
                                                                                                smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 2.38974e-06, No Iterations 8
                                                                                                smoothSolver:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::operator/(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#5  Foam::newFvPatchScalarField::updateCoeffs() at ??:?
#6  Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#7  Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#9  ? at ??:?
 #10  ? at ??:?
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

#12  ? at ??:?

[mtgoncalves]

This line:

#1  Foam::sigFpe::sigHandler(int) at ??:?

indicates that an error occured related to floating-point arithmetics.

This line:

#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

indicates that the error occured during the call to the function Foam::divide. Whether you called this function yourself or not you are facing a division error. You are probably dividing by 0 or by NaN. Make sure that you are not accidently dividing by zero. Check all the values of the variables by which you are dividing.

[generalKEE]

thanks for your help. Today, i meet another problem

myicoFoam: symbol lookup error: /home/a/OpenFOAM/a-5.x/platforms/linux64GccDPInt32Opt/lib/libnewFvPatchScalarField.so: undefined symbol: _ZN4Foam21newFvPatchScalarFieldC1ERKNS_7fvPatchERKNS_16DimensionedFieldIdNS_7volMeshEEERKNS_10dictionaryE

[mtgoncalves]

Undefined symbol means that there is a type which is not defined. I see you have developed your own solver myicoFoam and your own library for a BC called newFvPatchScalarField. As the lib is being picked up during compiling but some symbol is unknown I suspect that maybe the header file is not being picked up. It is hard to tell without seeing your Make/files and Make/options files.