[Bish]

Hi Foamers!

I am currently using OpenFOAM-dev, build dev-b9b8a8ac124e to model an explosion (specifically pressure piling) of two interconnected vessels (essentially 2 cylinders connected via a small tube or orifice). The domain is homogeneously filled with a gas such as hydrogen and ignited from one end of the chambers.

I am looking to use adaptive mesh refinement (AMR) to refine the mesh in areas of high turbulence as the simulation is falling over as the velocity in certain areas is going very high and therefore the Courant number crashes and the timestep goes very low (using adaptive time step). The mesh is approximately 8 million cells and is made up of tetrahedral cells.

The simulation runs but as soon as the criteria is met to refine the mesh, the following error:
[13] --> FOAM FATAL ERROR:
[13]
request for volVectorField rhoU_0 from objectRegistry region0 failed
available objects of type volVectorField are

3
(
U_0
U
rhoU
)

I am not sure where rhoU_0 is being called from but none of the dictionaries are.

Within the dynamicMeshDict, correctFluxes has to be defined and I believe it could have to do with this?
// Fluxes to adapt. For newly created faces or split faces the flux
// gets estimated from an interpolated volVectorField ('velocity')
// First is name of the flux to adapt, second is velocity that will
// be interpolated and inner-producted with the face area vector.
correctFluxes ((phi none) (phi none));


Looking forward to your help,
Bish

[Bish]

Has anyone got any thoughts on this error?

Thanks,
Bish