[hebbo]

Hello I try to use adjointShapOptimzationFoam with simple geometry in 3D (cube with intel and outlet), and when I launch the solver it stopped after some iterations
someone has any idea why time step continuity errors are very big

Time = 130

smoothSolver: Solving for Ux, Initial residual = 0.781005072918, Final residual = 0.0233840030481, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.712067168093, Final residual = 0.036631069425, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.515832116941, Final residual = 0.00395781721415, No Iterations 2
GAMG: Solving for p, Initial residual = 0.953967377828, Final residual = 0.000799590557173, No Iterations 2
time step continuity errors : sum local = 2.73866399586e+022, global = 46589779899.5, cumulative = 320771204800
smoothSolver: Solving for Uax, Initial residual = 0.727543426035, Final residual = 0.0374367869541, No Iterations 2
smoothSolver: Solving for Uay, Initial residual = 0.514450892256, Final residual = 0.0226805370686, No Iterations 2
smoothSolver: Solving for Uaz, Initial residual = 0.739721444744, Final residual = 0.0380601695909, No Iterations 2
GAMG: Solving for pa, Initial residual = 0.897579380387, Final residual = 0.0102502059598, No Iterations 1000
Adjoint continuity errors : sum local = 1.07202261853e+013, global = -39486445.3566, cumulative = -39486445.3826
smoothSolver: Solving for epsilon, Initial residual = 0.182253772738, Final residual = 6.35942170634e-013, No Iterations 2
bounding epsilon, min: -6.90189475985e+038 max: 1.90004723707e+047 average: 5.82196874261e+042
smoothSolver: Solving for k, Initial residual = 0.513925080089, Final residual = 0.00744794499962, No Iterations 2
bounding k, min: -1.81554379651e+024 max: 1.66279311428e+035 average: 1.55756260899e+031
ExecutionTime = 7.798 s

Time = 140

smoothSolver: Solving for Ux, Initial residual = 3.81727510694e-009, Final residual = 3.11489470734e-010, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1.25511052017e-008, Final residual = 1.16892199155e-010, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 3.48885137386e-009, Final residual = 2.45941972485e-011, No Iterations 2
GAMG: Solving for p, Initial residual = 0.570942089878, Final residual = 0.00260936431887, No Iterations 3
time step continuity errors : sum local = 2.18696253546e+031, global = 4.72054741864e+027, cumulative = 4.72054741864e+027
smoothSolver: Solving for Uax, Initial residual = 0.000279435712908, Final residual = 1.41112899367e-005, No Iterations 2
smoothSolver: Solving for Uay, Initial residual = 0.0155509364012, Final residual = 0.000178602548244, No Iterations 2
smoothSolver: Solving for Uaz, Initial residual = 0.000603922931721, Final residual = 4.56602634073e-006, No Iterations 2
GAMG: Solving for pa, Initial residual = 0.998558574882, Final residual = 0.134155188336, No Iterations 1000
Adjoint continuity errors : sum local = 4.44802225035e+020, global = 1.02952664845e+019, cumulative = 1.02952664844e+019
smoothSolver: Solving for epsilon, Initial residual = 3.04170188515e-013, Final residual = 3.04170188515e-013, No Iterations 0
bounding epsilon, min: 1.04375481111e-016 max: 3.14166387103e+052 average: 1.48808714507e+048
smoothSolver: Solving for k, Initial residual = 0.000238301599633, Final residual = 8.07487490198e-006, No Iterations 4
ExecutionTime = 12.457 s


Time = 150

smoothSolver: Solving for Ux, Initial residual = 0.0470877727959, Final residual = 0.00130757778594, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.505192425985, Final residual = 0.0251640250965, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0231806609558, Final residual = 2.79628829666e-005, No Iterations 2
GAMG: Solving for p, Initial residual = 0.79276045501, Final residual = 0.00589673078511, No Iterations 2
time step continuity errors : sum local = 3.13972580956e+036, global = 3.35553345637e+026, cumulative = 5.05610076428e+027
smoothSolver: Solving for Uax, Initial residual = 0.0264664497917, Final residual = 0.000434151868154, No Iterations 2
smoothSolver: Solving for Uay, Initial residual = 0.0105814970443, Final residual = 0.000156570829142, No Iterations 2
smoothSolver: Solving for Uaz, Initial residual = 0.0125595018958, Final residual = 0.000399637645081, No Iterations 2
Generating stack trace...

[MaryBau]

Hi hebbo,

I am not familiar with adjointShapOptimzationFoam, but here are some ideas based on the fact that your simulation has big errors right from the start:

- checkMesh -allTopology -allGeometry. Does your mesh have any errors?

- should you reduce your time step?

- Are your initial field and/or boundary conditions are well defined?

I hope it helps, good luck,

Mary