[rieuk]

Hey guys I'm simulating a (2D) transversely oscillating cylinder using a non-inertial reference frame using the following modification to icoFoam:

scalar Pi = mathematicalConstant::pi;
dimensionedScalar accCt("accCt",dimensionSet(0,1,-2,0,0,0,0),-3.0*pow(2.0*Pi/3.0,2)*Foam::cos(2.0*Pi*runTime.value()/3.0));
dimensionedVector acc("acc",dimensionSet(0,1,-2,0,0,0,0), vector::zero);

acc.component(vector::Y) = accCt;

        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(nu, U)
    ==
       acc
        );


        solve(UEqn ==  - fvc::grad(p));

However the run keeps blowing up (Floating Point Exception with ridiculously high velocity magnitude) at around the same time even for very small delta T... can anyone see what's wrong with my implementation?