CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Calculating fvm::laplacian(tensor,scalar) on cylinder

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 3, 2014, 15:05
Default Calculating fvm::laplacian(tensor,scalar) on cylinder
  #1
New Member
 
James Montague
Join Date: Jan 2014
Location: Vermont, United States
Posts: 7
Rep Power: 12
Montague is on a distinguished road
I am using OpenFoam 2.2.2. I'm getting the following error when trying to solve laplacianFoam for a cylinder using a slight modification in which laplacianFoam takes a tensor as the coefficient D.

#0 Foam::error::printStack(Foam::Ostream&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::fvc::surfaceIntegrate<double>(Foam::Field<do uble>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<double>(Foam::Geometri cField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gaussLaplacianScheme<double, Foam::Tensor<double> >::fvmLaplacian(Foam::GeometricField<Foam::Tensor< double>, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::laplacianScheme<double, Foam::Tensor<double> >::fvmLaplacian(Foam::GeometricField<Foam::Tensor< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/bin/laplacianFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/bin/laplacianFoam"
Floating point exception (core dumped)

The problem appears to be with the mesh that I am using because the same solver will run just fine for the provided flange example. The solver also works on the cylinder mesh if a scalar is used for D.

I have attached both the model files in cylinder.zip and the slightly modified solver so that it uses a tensor coefficient in laplacianFoamT.zip.

Any insight into how I can get this to work would be much appreciated.

Thanks,
Jim
Attached Files
File Type: zip cylinder.zip (9.8 KB, 4 views)
File Type: zip laplacianFoamT.zip (3.3 KB, 23 views)
Montague is offline   Reply With Quote

Old   June 20, 2014, 11:12
Default
  #2
New Member
 
James Montague
Join Date: Jan 2014
Location: Vermont, United States
Posts: 7
Rep Power: 12
Montague is on a distinguished road
It appears that the original mesh that I used with blockMesh was not suitable. I have remeshed the cylinder using snappyHexMesh and the swiftSnap tool developed for blender (it works really well!)

The revised case is attached below.

-James
Attached Files
File Type: zip cylinder.zip (20.4 KB, 10 views)
Montague is offline   Reply With Quote

Reply

Tags
cylinder mesh, laplacianfoam, surfaceintegrate, tensor

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Calculating Sloshing mass in a cylinder satellite_control FLOW-3D 0 June 20, 2011 05:34
[blockMesh] Specifying boundary faces failes in blockMesh blaise OpenFOAM Meshing & Mesh Conversion 0 May 10, 2010 03:56
calculating the pressure on cylinder for immersed boundary technique chaithanya Main CFD Forum 0 May 19, 2009 09:28
How to update polyPatchbs localPoints liu OpenFOAM Running, Solving & CFD 6 December 30, 2005 17:27
Calculating shedding frequency of square cylinder zonexo Main CFD Forum 5 November 17, 2005 18:12


All times are GMT -4. The time now is 06:46.