[WUYing]

Hi, I have met a convergence problem in OF when simulating the FSI of cylinder. The CFL becomes larger until around 30,000 at the last timestep before crash. I also checked the pressure and velocity fields in Paraview and find they are obviously unreasonable at the very beginning.

The following is the error report:

#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::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/ying/OpenFOAM41/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
#8  ? in "/home/ying/OpenFOAM41/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/home/ying/OpenFOAM41/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
Floating point exception (core dumped)

I check the similar posts like:
Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOp
https://www.simscale.com/forum/t/cfd...-failing/82914
https://www.salome-platform.org/foru...m_10/549857363

but very limited information can be obtained....

I use a refine mesh to rerun and the case can last for more time steps (coarse mesh crash at time = 1.8 and refine mesh crash at time=18) but U and P are still incorrect at the very beginning.

The attachments include:
(1) case.zip:
0.zip (U, p, pointDisplacement)
system.zip (controlDict, fvScheme, fvSolution)
constant.zip (dynamicMeshDict, transportProperties, turbulenceProperties)
(2) log.zip (log.pimpleDyMFoam file)
(3) coarse.zip (coarse mesh, and corresponding (U, p) images at last timestep)
(4) refine.zip (refine mesh, and corresponding (U, p) images at last timestep)

Due to the size limit, I cannot upload my mesh and animation of flow fields, please send me email if you want to obtain the complete case file.

I am very appreciated for any suggestion!

[Fouch]

Hi,

Firstly, you should add maxDeltaT and maxCo in your controlDict to avoid any blow up of the time step.

I suggest to test it and revert to us.

Best Regards,

[WUYing]

Originally Posted by Fouch

Hi,

Firstly, you should add maxDeltaT and maxCo in your controlDict to avoid any blow up of the time step.

I suggest to test it and revert to us.

Best Regards,

Originally Posted by Fouch

Hi,

Firstly, you should add maxDeltaT and maxCo in your controlDict to avoid any blow up of the time step.

I suggest to test it and revert to us.

Best Regards,

Hi Fouch, thanks for your suggestion!

Actually, I originally set the adjustTimeStep as false and specify fixed deltaT, so the time step will not blow up. And I try to decrease deltaT but the case still crashes at very big CFL. For example, using the same mesh, the case crashes at time = 1.83 when deltaT = 0.01, crashes at time = 1.775 when deltaT = 0.005, crashes at time = 0.368 when deltaT = 0.001.

According to your suggestion, I also use the adjustable timestep by adding the following 3 lines before function in the controlDict:

adjustTimeStep  on; 

maxCo           2;  

maxDeltaT       0.1;

And the case would also crash at about time = 1.8. U and p are obviously incorrect.

I also attach part of log.pimpleDyMFoam (due to size limit) here.

Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: sixDoFRigidBodyMotion
Selecting sixDoFSolver Newmark
Translational constraint tensor (0 0 0 0 1 0 0 0 0)
Rotational constraint tensor (0 0 0 0 0 0 0 0 0)

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present

Reading/calculating face velocity Uf

No finite volume options present

Courant Number mean: 1.00047e-05 max: 0.000610799
fieldAverage fieldAverage1:
    Starting averaging at time 0

forces forces:
    Not including porosity effects
forces forceCoeffs:
    Not including porosity effects
forceCoeffs forceCoeffs:
    Not including porosity effects
Reading surface description:
    zNormal


Starting time loop

Courant Number mean: 1.00037e-05 max: 0.000610738
deltaT = 0.00119976
Time = 0.00119976

forces forces:
    Not including porosity effects
Restraint verticalSpring:  attachmentPt - anchor (0 0 0) spring length 0 force (-0 -0 -0)
6-DoF rigid body motion
    Centre of rotation: (0.280485 0 0)
    Centre of mass: (0.280485 0 0)
    Orientation: (1 0 0 0 1 0 0 0 1)
    Linear velocity: (0 0 0)
    Angular velocity: (0 0 0)
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.0196279, No Iterations 12
GAMG:  Solving for pcorr, Initial residual = 0.0199319, Final residual = 0.0199319, No Iterations 0
time step continuity errors : sum local = 6.71518e-05, global = -6.13618e-08, cumulative = -6.13618e-08
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 1.45212e-06, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0300958, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00738863, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.0194856, Final residual = 0.000119229, No Iterations 6
time step continuity errors : sum local = 4.03179e-07, global = -7.0946e-09, cumulative = -6.84564e-08
GAMG:  Solving for p, Initial residual = 0.0212682, Final residual = 0.00011576, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.0141619, Final residual = 8.47907e-05, No Iterations 6
time step continuity errors : sum local = 2.89751e-07, global = -5.20169e-09, cumulative = -7.36581e-08
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.555428, Final residual = 1.10649e-07, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.0382701, Final residual = 2.11528e-07, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0218124, Final residual = 0.000131796, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.0144855, Final residual = 8.75873e-05, No Iterations 6
time step continuity errors : sum local = 2.98957e-07, global = -5.5809e-09, cumulative = -7.9239e-08
GAMG:  Solving for p, Initial residual = 0.00406269, Final residual = 3.32292e-05, No Iterations 9
GAMG:  Solving for p, Initial residual = 0.00199958, Final residual = 4.88785e-08, No Iterations 17
time step continuity errors : sum local = 1.71296e-10, global = -1.77625e-12, cumulative = -7.92407e-08
ExecutionTime = 0.81 s  ClockTime = 1 s

    Reading/initialising field UMean
    Reading/initialising field pMean
    Reading/initialising field UPrime2Mean
    Reading/initialising field pPrime2Mean
fieldAverage fieldAverage1 write:
    Calculating averages

forces forces write:
    sum of forces:
        pressure : (2733.32 0.306721 2.04533e-15)
        viscous  : (12.2219 -0.0430355 -5.66497e-16)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-1.40488e-15 3.03509e-16 0.311642)
        viscous  : (4.34042e-16 -2.32209e-16 -0.00189707)
        porous   : (0 0 0)

....

Courant Number mean: 0.0190881 max: 3.12244
deltaT = 3.59939e-21
Time = 1.883575796235934

forces forces:
    Not including porosity effects
Restraint verticalSpring:  attachmentPt - anchor (0 -9.23421 0) spring length 9.23421 force (-0 98.8702 -0)
6-DoF rigid body motion
    Centre of rotation: (0.280485 -9.24123 0)
    Centre of mass: (0.280485 -9.24123 0)
    Orientation: (1 0 0 0 1 0 0 0 1)
    Linear velocity: (-0 -1.382e+18 -0)
    Angular velocity: (0 0 0)
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.017213, No Iterations 13
GAMG:  Solving for pcorr, Initial residual = 0.472833, Final residual = 0.0116254, No Iterations 5
time step continuity errors : sum local = 5.49285e-07, global = 1.49488e-08, cumulative = -1.04105e-05
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.0876311, Final residual = 0.000118407, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0125489, Final residual = 1.59195e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.389902, Final residual = 0.00248264, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.553004, Final residual = 0.00445544, No Iterations 7
time step continuity errors : sum local = 3.29021e-07, global = 1.29096e-08, cumulative = -1.03976e-05
GAMG:  Solving for p, Initial residual = 0.621976, Final residual = 0.00358929, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.446588, Final residual = 0.00358063, No Iterations 7
time step continuity errors : sum local = 2.32644e-07, global = 8.9411e-09, cumulative = -1.03886e-05
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.0111158, Final residual = 2.48451e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00147112, Final residual = 2.51958e-08, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.633415, Final residual = 0.00361362, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.465718, Final residual = 0.003714, No Iterations 7
time step continuity errors : sum local = 2.50486e-07, global = 9.62316e-09, cumulative = -1.0379e-05
GAMG:  Solving for p, Initial residual = 0.119757, Final residual = 0.00082459, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.0346993, Final residual = 7.58565e-08, No Iterations 26
time step continuity errors : sum local = 8.35769e-12, global = 4.16261e-14, cumulative = -1.0379e-05
ExecutionTime = 793.06 s  ClockTime = 1597 s

fieldAverage fieldAverage1 write:
    Calculating averages

forces forces write:
    sum of forces:
        pressure : (-1.50783e+37 -5.49819e+39 1.52322e+19)
        viscous  : (-8.90904e+18 2.87163e+19 -59.2447)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-1.36597e+20 -6.5933e+18 -6.88034e+39)
        viscous  : (548.733 7.40703 -1.02861e+20)
        porous   : (0 0 0)

forceCoeffs forceCoeffs write:
    Cm    = -1.12332e+40
    Cd    = -2.46177e+37
    Cl    = -8.97663e+39
    Cl(f) = -1.57215e+40
    Cl(r) = 6.74489e+39

Courant Number mean: 0.00837317 max: 1.43518
deltaT = 4.10099e-21
Time = 1.883575796235934

forces forces:
    Not including porosity effects
Restraint verticalSpring:  attachmentPt - anchor (0 -9.24123 0) spring length 9.24123 force (-0 98.9453 -0)
6-DoF rigid body motion
    Centre of rotation: (0.280485 -9.24755 0)
    Centre of mass: (0.280485 -9.24755 0)
    Orientation: (1 0 0 0 1 0 0 0 1)
    Linear velocity: (0 -1.69837e+18 0)
    Angular velocity: (0 0 0)
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.01525, No Iterations 14
GAMG:  Solving for pcorr, Initial residual = 0.471027, Final residual = 0.0188302, No Iterations 4
time step continuity errors : sum local = 6.06274e-07, global = 1.3257e-08, cumulative = -1.03657e-05
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.220999, Final residual = 0.000216505, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0302858, Final residual = 2.76543e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.150527, Final residual = 0.00112312, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.226602, Final residual = 0.00217593, No Iterations 7
time step continuity errors : sum local = 1.41945e-07, global = -7.71576e-09, cumulative = -1.03735e-05
GAMG:  Solving for p, Initial residual = 0.128788, Final residual = 0.000958094, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.184062, Final residual = 0.00175177, No Iterations 7
time step continuity errors : sum local = 1.00278e-07, global = -5.32282e-09, cumulative = -1.03788e-05
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.0986745, Final residual = 2.10818e-08, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.021731, Final residual = 1.6546e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.162085, Final residual = 0.00142636, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.341383, Final residual = 0.00320235, No Iterations 7
time step continuity errors : sum local = 2.16111e-07, global = -1.16693e-08, cumulative = -1.03905e-05
GAMG:  Solving for p, Initial residual = 0.318819, Final residual = 0.00169682, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.0983, Final residual = 9.09694e-08, No Iterations 25
time step continuity errors : sum local = 2.31849e-12, global = 1.75734e-13, cumulative = -1.03905e-05
ExecutionTime = 793.68 s  ClockTime = 1599 s

fieldAverage fieldAverage1 write:
    Calculating averages

forces forces write:
    sum of forces:
        pressure : (-5.92266e+36 -8.01276e+38 3.23412e+18)
        viscous  : (-8.85916e+18 1.7958e+19 -23.3363)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-3.04516e+19 -8.96344e+17 -1.04569e+39)
        viscous  : (219.798 6.38085 -9.69291e+19)
        porous   : (0 0 0)

forceCoeffs forceCoeffs write:
    Cm    = -1.70725e+39
    Cd    = -9.66964e+36
    Cl    = -1.3082e+39
    Cl(f) = -2.36136e+39
    Cl(r) = 1.05315e+39

Courant Number mean: 0.00692353 max: 1.16656
deltaT = 4.80407e-21
Time = 1.883575796235934

forces forces:
    Not including porosity effects
Restraint verticalSpring:  attachmentPt - anchor (0 -9.24755 0) spring length 9.24755 force (-0 99.013 -0)
6-DoF rigid body motion
    Centre of rotation: (0.280485 -9.25888 0)
    Centre of mass: (0.280485 -9.25888 0)
    Orientation: (1 0 0 0 1 0 0 0 1)
    Linear velocity: (0 -3.01891e+18 0)
    Angular velocity: (0 0 0)
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.0194931, No Iterations 13
GAMG:  Solving for pcorr, Initial residual = 0.475357, Final residual = 0.0117041, No Iterations 5
time step continuity errors : sum local = 1.28399e-06, global = -3.8277e-08, cumulative = -1.04287e-05
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.0622515, Final residual = 0.000173185, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0108993, Final residual = 2.75533e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.693178, Final residual = 0.00518598, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.496134, Final residual = 0.00397948, No Iterations 7
time step continuity errors : sum local = 5.39988e-07, global = -1.98085e-08, cumulative = -1.04485e-05
GAMG:  Solving for p, Initial residual = 0.36376, Final residual = 0.00223436, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.304788, Final residual = 0.00242902, No Iterations 7
time step continuity errors : sum local = 3.75711e-07, global = -1.33077e-08, cumulative = -1.04618e-05
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.0296349, Final residual = 3.19176e-07, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.00491257, Final residual = 5.53792e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.29698, Final residual = 0.00258741, No Iterations 13
GAMG:  Solving for p, Initial residual = 0.23883, Final residual = 0.00196561, No Iterations 7
time step continuity errors : sum local = 3.33194e-07, global = -1.10358e-08, cumulative = -1.04729e-05
GAMG:  Solving for p, Initial residual = 0.0799564, Final residual = 0.000502528, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.0243955, Final residual = 6.41177e-08, No Iterations 25
time step continuity errors : sum local = 1.51374e-11, global = 4.87715e-14, cumulative = -1.04729e-05
ExecutionTime = 794.27 s  ClockTime = 1600 s

fieldAverage fieldAverage1 write:
    Calculating averages

forces forces write:
    sum of forces:
        pressure : (-2.26653e+37 7.39699e+39 -5.44652e+19)
        viscous  : (-8.76194e+18 4.3249e+19 130.447)
        porous   : (0 0 0)
    sum of moments:
        pressure : (5.08955e+20 2.40672e+19 8.91389e+39)
        viscous  : (-1195.8 -87.193 -1.10644e+20)
        porous   : (0 0 0)

forceCoeffs forceCoeffs write:
    Cm    = 1.45533e+40
    Cd    = -3.70045e+37
    Cl    = 1.20767e+40
    Cl(f) = 2.05916e+40
    Cl(r) = -8.51492e+39

Courant Number mean: 0.0170929 max: 2.75666
deltaT = 3.48543e-21
Time = 1.883575796235934

forces forces:
    Not including porosity effects
Restraint verticalSpring:  attachmentPt - anchor (0 -9.25888 0) spring length 9.25888 force (-0 99.1343 -0)
6-DoF rigid body motion
    Centre of rotation: (0.280485 -9.26765 0)
    Centre of mass: (0.280485 -9.26765 0)
    Orientation: (1 0 0 0 1 0 0 0 1)
    Linear velocity: (-0 -2.01578e+18 -0)
    Angular velocity: (0 0 0)
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.0151311, No Iterations 14
GAMG:  Solving for pcorr, Initial residual = 0.476142, Final residual = 0.0191458, No Iterations 4
time step continuity errors : sum local = 2.22023e-07, global = 4.91909e-09, cumulative = -1.0468e-05
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.140065, Final residual = 0.000267872, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0211073, Final residual = 3.56625e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.267493, Final residual = 0.002458, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.514274, Final residual = 0.00461277, No Iterations 7
time step continuity errors : sum local = 3.23868e-07, global = 2.21413e-08, cumulative = -1.04458e-05
GAMG:  Solving for p, Initial residual = 0.318334, Final residual = 0.00275394, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.507438, Final residual = 0.00468354, No Iterations 7
time step continuity errors : sum local = 2.35224e-07, global = 1.26793e-08, cumulative = -1.04331e-05
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.037576, Final residual = 9.85391e-08, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.00479482, Final residual = 1.87923e-07, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.283688, Final residual = 0.00266262, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.550086, Final residual = 0.00506526, No Iterations 7
time step continuity errors : sum local = 3.24593e-07, global = 1.73568e-08, cumulative = -1.04158e-05
GAMG:  Solving for p, Initial residual = 0.238318, Final residual = 0.00237296, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.0725705, Final residual = 5.86767e-08, No Iterations 26
time step continuity errors : sum local = 3.54871e-12, global = -2.73579e-13, cumulative = -1.04158e-05
ExecutionTime = 794.89 s  ClockTime = 1601 s

fieldAverage fieldAverage1 write:
    Calculating averages

forces forces write:
    sum of forces:
        pressure : (-2.60458e+37 -3.00623e+39 1.41423e+19)
        viscous  : (-9.85558e+18 4.19069e+19 -13.4112)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-1.37712e+20 1.01025e+17 -3.87417e+39)
        viscous  : (140.312 -3.8486 -1.19684e+20)
        porous   : (0 0 0)

forceCoeffs forceCoeffs write:
    Cm    = -6.32518e+39
    Cd    = -4.25238e+37
    Cl    = -4.90814e+39
    Cl(f) = -8.77925e+39
    Cl(r) = 3.87111e+39

[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  void Foam::vtkSurfaceWriter::writeData<double>(Foam::Ostream&, Foam::Field<double> const&) at ??:?
[0] #4  void Foam::vtkSurfaceWriter::writeTemplate<double>(Foam::fileName const&, Foam::fileName const&, Foam::Field<Foam::Vector<double> > const&, Foam::List<Foam::face> const&, Foam::word const&, Foam::Field<double> const&, bool, bool) const at ??:?
[0] #5  void Foam::sampledSurfaces::writeSurface<double>(Foam::Field<double> const&, int, Foam::word const&, Foam::fileName const&) at ??:?
[0] #6  void Foam::sampledSurfaces::sampleAndWrite<double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[0] #7  void Foam::sampledSurfaces::sampleAndWrite<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >(Foam::IOobjectList const&) at ??:?
[0] #8  Foam::sampledSurfaces::write() at ??:?
[0] #9  Foam::functionObjects::timeControl::write() at ??:?
[0] #10  Foam::functionObjectList::execute() at ??:?
[0] #11  Foam::Time::run() const at Time.C:?
[0] #12  ? in "/home/ying/OpenFOAM41/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/pimpleDyMFoam"
[0] #13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #14  ?[ubuntu:62701] *** Process received signal ***

U and p fields at last time step are also attached.

Thanks!