Trouble running adjointOptimisationFoam

lowdo · August 3, 2023, 11:28

Hi,

I'm fairly new to OpenFOAM so apologies in advance if I'm making basic errors.

I have a model which is simply an isolated wing, which I have managed to run successfully in simpleFoam (or at least, forces seemed roughly in line with what I was expecting). I now want to carry out an optimisation run with the adjoint solver, but am finding that after a few iterations of the primal solver it blows up. As far as I can tell, I've matched settings for the two cases so I'm a bit confused as to why this is happening.

If it helps, this is the output I'm getting from the adjoint solver:

Starting time loop

* * * * * * * * * * * * * * * * *
Optimisation cycle 1
* * * * * * * * * * * * * * * * *

Solving equations for solver op1

Time = 0.001

smoothSolver: Solving for Ux, Initial residual = 0.9999999999803243, Final residual = 0.07772845287810633, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.07684974481884306, No Iterations 5
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.07688831165285209, No Iterations 5
GAMG: Solving for p, Initial residual = 1, Final residual = 0.009605708128934338, No Iterations 35
GAMG: Solving for p, Initial residual = 3.535826797840461e-05, Final residual = 2.858242669232686e-07, No Iterations 2
GAMG: Solving for p, Initial residual = 7.126941704764867e-06, Final residual = 7.550606703926305e-08, No Iterations 2
GAMG: Solving for p, Initial residual = 2.102514646551255e-06, Final residual = 3.288155961099869e-08, No Iterations 2
time step continuity errors : sum local = 1.229229421408113e-07, global = -2.289088115992725e-09, cumulative = -2.289088115992725e-09
smoothSolver: Solving for omega, Initial residual = 0.0348456690990019, Final residual = 0.001951784080296146, No Iterations 4
bounding omega, min: -311.2361631697304 max: 14065451671.24115 average: 24208956.33683804
smoothSolver: Solving for k, Initial residual = 0.9999999961607311, Final residual = 0.06096080113120984, No Iterations 5

drag : 193.4351321638371
lift : -30.96791059040253
ExecutionTime = 258.68 s ClockTime = 258 s

Time = 0.002

smoothSolver: Solving for Ux, Initial residual = 0.03753072802341328, Final residual = 0.002853946768200444, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.07759706471823789, Final residual = 0.005018013040632728, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.1517784302589093, Final residual = 0.01254678018543783, No Iterations 4
GAMG: Solving for p, Initial residual = 0.01756602361905522, Final residual = 0.0001509626001795971, No Iterations 6
GAMG: Solving for p, Initial residual = 0.001284318366671466, Final residual = 1.284235130522751e-05, No Iterations 33
GAMG: Solving for p, Initial residual = 0.001739859030886709, Final residual = 1.697019549732851e-05, No Iterations 13
GAMG: Solving for p, Initial residual = 0.000308714611519725, Final residual = 3.041660057729567e-06, No Iterations 68
time step continuity errors : sum local = 5.271164838488112e-07, global = 8.257266560506197e-09, cumulative = 5.968178444513471e-09
smoothSolver: Solving for omega, Initial residual = 0.01206122025547818, Final residual = 0.0011679470510121, No Iterations 1
bounding omega, min: -23660024.04731201 max: 171181639029.5126 average: 1093346004.441741
smoothSolver: Solving for k, Initial residual = 0.9219684983364, Final residual = 0.003708304919700916, No Iterations 1

drag : 130.3339349542009
lift : -21.91366791043179
ExecutionTime = 320.58 s ClockTime = 320 s

Time = 0.003

smoothSolver: Solving for Ux, Initial residual = 9.513599076584994e-05, Final residual = 8.41811716628903e-06, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.0001238538077827657, Final residual = 9.601944630361708e-06, No Iterations 4
smoothSolver: Solving for Uz, Initial residual = 0.0002937497769088864, Final residual = 2.886902843679186e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.02271539311957872, Final residual = 0.0002132644793497413, No Iterations 7
GAMG: Solving for p, Initial residual = 0.0005770742137256225, Final residual = 5.745209906667744e-06, No Iterations 66
GAMG: Solving for p, Initial residual = 0.0003408951909632294, Final residual = 3.327806393769065e-06, No Iterations 17
GAMG: Solving for p, Initial residual = 4.422513155526133e-05, Final residual = 4.388078456171613e-07, No Iterations 92
time step continuity errors : sum local = 7.755764485770775e-07, global = 2.613534972884946e-08, cumulative = 3.210352817336293e-08
smoothSolver: Solving for omega, Initial residual = 0.03777015413861053, Final residual = 0.0002811468050693801, No Iterations 1
bounding omega, min: -9386755640.150707 max: 6052647289138.035 average: 417880268.3301896
smoothSolver: Solving for k, Initial residual = 0.03655561646644971, Final residual = 0.001358430484242084, No Iterations 1
bounding k, min: -0.124148961509257 max: 13.54980413364046 average: 0.5762411249125505

drag : -219.7520678407049
lift : -36.38409258194373
ExecutionTime = 405.41 s ClockTime = 405 s

Time = 0.004

smoothSolver: Solving for Ux, Initial residual = 0.3135682229449583, Final residual = 0.01065606675022723, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.2952628168401448, Final residual = 0.007776041656604311, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.2294808132606518, Final residual = 0.005939966102981783, No Iterations 2
GAMG: Solving for p, Initial residual = 0.01142465009285283, Final residual = 0.0001126659433648541, No Iterations 19
GAMG: Solving for p, Initial residual = 1.289434024409072e-13, Final residual = 1.289434024409072e-13, No Iterations 0
GAMG: Solving for p, Initial residual = 1.289434024409072e-13, Final residual = 1.289434024409072e-13, No Iterations 0
GAMG: Solving for p, Initial residual = 1.289434024409072e-13, Final residual = 1.289434024409072e-13, No Iterations 0
time step continuity errors : sum local = 0.001043798483872672, global = 4.972041248134074e-08, cumulative = 8.182394065470367e-08
smoothSolver: Solving for omega, Initial residual = 0.9964082542167786, Final residual = 0.068985186050248, No Iterations 6
bounding omega, min: -180036774609.3387 max: 6.329118058827089e+27 average: 1.089617791183109e+22
smoothSolver: Solving for k, Initial residual = 0.1882230893833461, Final residual = 5.367243120940624e-05, No Iterations 1
bounding k, min: -2.518099445090961 max: 95.21876859574435 average: 1.27445342716669

drag : -2869.298567197414
lift : -168.7879710025276
ExecutionTime = 429.36 s ClockTime = 429 s

Time = 0.005

smoothSolver: Solving for Ux, Initial residual = 0.1448256781172823, Final residual = 0.003454284159427664, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.3306804680012086, Final residual = 0.006427471134016752, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.3215369458304032, Final residual = 0.006512543729490321, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0002613290470658523, Final residual = 2.170007262345748e-06, No Iterations 5
GAMG: Solving for p, Initial residual = 8.867425281744662e-06, Final residual = 9.66013942785868e-08, No Iterations 12
GAMG: Solving for p, Initial residual = 4.241496950591247e-06, Final residual = 9.65159104422725e-08, No Iterations 5
GAMG: Solving for p, Initial residual = 2.057181654485832e-06, Final residual = 9.279290081773349e-08, No Iterations 2
time step continuity errors : sum local = 1024.804129831097, global = -0.07002320644511717, cumulative = -0.07002312462117652
smoothSolver: Solving for omega, Initial residual = 0.5386999546680851, Final residual = 1.456360013262606e-09, No Iterations 1
bounding omega, min: -1.669470785195285e+25 max: 8.034849210266171e+32 average: 7.548360186268818e+25
smoothSolver: Solving for k, Initial residual = 0.07480861504188824, Final residual = 2.47095546180343e-08, No Iterations 1
bounding k, min: -15.45926119247786 max: 603.0051979285942 average: 2.288263708045382

drag : 7221162522.83777
lift : -5361883058.663867
ExecutionTime = 454.64 s ClockTime = 454 s

Time = 0.006

smoothSolver: Solving for Ux, Initial residual = 0.3928042942650474, Final residual = 0.03693958221338808, No Iterations 17
smoothSolver: Solving for Uy, Initial residual = 0.3655420738309382, Final residual = 0.02841612700001563, No Iterations 18
smoothSolver: Solving for Uz, Initial residual = 0.3344446943495146, Final residual = 0.02705236061683774, No Iterations 17
GAMG: Solving for p, Initial residual = 0.9999876144347357, Final residual = 165.5982447656887, No Iterations 1000
GAMG: Solving for p, Initial residual = 1.744015098587018e-09, Final residual = 1.744015098587018e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.744015098587018e-09, Final residual = 1.744015098587018e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1.744015098587018e-09, Final residual = 1.744015098587018e-09, No Iterations 0
time step continuity errors : sum local = 3263977569297249, global = -1136250716544.243, cumulative = -1136250716544.313
smoothSolver: Solving for omega, Initial residual = 0.5384603217155158, Final residual = 8.440472411405353e-08, No Iterations 1
bounding omega, min: -5.634206076641107e+31 max: 2.410454853511547e+32 average: 3.159597015085345e+25
smoothSolver: Solving for k, Initial residual = 0.1272320936674846, Final residual = 0.0104019729237369, No Iterations 4
bounding k, min: -129.7591077230691 max: 637.4025949508549 average: 2.099637272606892

drag : 1.097508704882627e+27
lift : -1.381264669290234e+27
ExecutionTime = 847.5599999999999 s ClockTime = 848 s

Time = 0.007

--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

Process name: [[59045,1],1]
Exit code: 142
--------------------------------------------------------------------------

Thanks,

Chris

lowdo · August 21, 2023, 04:36

Is there anyone that can help with this?

August 3, 2023, 11:28	Trouble running adjointOptimisationFoam	#1
lowdo New Member Chris Join Date: Feb 2013 Posts: 17 Rep Power: 13	Hi, I'm fairly new to OpenFOAM so apologies in advance if I'm making basic errors. I have a model which is simply an isolated wing, which I have managed to run successfully in simpleFoam (or at least, forces seemed roughly in line with what I was expecting). I now want to carry out an optimisation run with the adjoint solver, but am finding that after a few iterations of the primal solver it blows up. As far as I can tell, I've matched settings for the two cases so I'm a bit confused as to why this is happening. If it helps, this is the output I'm getting from the adjoint solver: Starting time loop * * * * * * * * * * * * * * * * * Optimisation cycle 1 * * * * * * * * * * * * * * * * * Solving equations for solver op1 Time = 0.001 smoothSolver: Solving for Ux, Initial residual = 0.9999999999803243, Final residual = 0.07772845287810633, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.07684974481884306, No Iterations 5 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.07688831165285209, No Iterations 5 GAMG: Solving for p, Initial residual = 1, Final residual = 0.009605708128934338, No Iterations 35 GAMG: Solving for p, Initial residual = 3.535826797840461e-05, Final residual = 2.858242669232686e-07, No Iterations 2 GAMG: Solving for p, Initial residual = 7.126941704764867e-06, Final residual = 7.550606703926305e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 2.102514646551255e-06, Final residual = 3.288155961099869e-08, No Iterations 2 time step continuity errors : sum local = 1.229229421408113e-07, global = -2.289088115992725e-09, cumulative = -2.289088115992725e-09 smoothSolver: Solving for omega, Initial residual = 0.0348456690990019, Final residual = 0.001951784080296146, No Iterations 4 bounding omega, min: -311.2361631697304 max: 14065451671.24115 average: 24208956.33683804 smoothSolver: Solving for k, Initial residual = 0.9999999961607311, Final residual = 0.06096080113120984, No Iterations 5 drag : 193.4351321638371 lift : -30.96791059040253 ExecutionTime = 258.68 s ClockTime = 258 s Time = 0.002 smoothSolver: Solving for Ux, Initial residual = 0.03753072802341328, Final residual = 0.002853946768200444, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.07759706471823789, Final residual = 0.005018013040632728, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.1517784302589093, Final residual = 0.01254678018543783, No Iterations 4 GAMG: Solving for p, Initial residual = 0.01756602361905522, Final residual = 0.0001509626001795971, No Iterations 6 GAMG: Solving for p, Initial residual = 0.001284318366671466, Final residual = 1.284235130522751e-05, No Iterations 33 GAMG: Solving for p, Initial residual = 0.001739859030886709, Final residual = 1.697019549732851e-05, No Iterations 13 GAMG: Solving for p, Initial residual = 0.000308714611519725, Final residual = 3.041660057729567e-06, No Iterations 68 time step continuity errors : sum local = 5.271164838488112e-07, global = 8.257266560506197e-09, cumulative = 5.968178444513471e-09 smoothSolver: Solving for omega, Initial residual = 0.01206122025547818, Final residual = 0.0011679470510121, No Iterations 1 bounding omega, min: -23660024.04731201 max: 171181639029.5126 average: 1093346004.441741 smoothSolver: Solving for k, Initial residual = 0.9219684983364, Final residual = 0.003708304919700916, No Iterations 1 drag : 130.3339349542009 lift : -21.91366791043179 ExecutionTime = 320.58 s ClockTime = 320 s Time = 0.003 smoothSolver: Solving for Ux, Initial residual = 9.513599076584994e-05, Final residual = 8.41811716628903e-06, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.0001238538077827657, Final residual = 9.601944630361708e-06, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.0002937497769088864, Final residual = 2.886902843679186e-05, No Iterations 4 GAMG: Solving for p, Initial residual = 0.02271539311957872, Final residual = 0.0002132644793497413, No Iterations 7 GAMG: Solving for p, Initial residual = 0.0005770742137256225, Final residual = 5.745209906667744e-06, No Iterations 66 GAMG: Solving for p, Initial residual = 0.0003408951909632294, Final residual = 3.327806393769065e-06, No Iterations 17 GAMG: Solving for p, Initial residual = 4.422513155526133e-05, Final residual = 4.388078456171613e-07, No Iterations 92 time step continuity errors : sum local = 7.755764485770775e-07, global = 2.613534972884946e-08, cumulative = 3.210352817336293e-08 smoothSolver: Solving for omega, Initial residual = 0.03777015413861053, Final residual = 0.0002811468050693801, No Iterations 1 bounding omega, min: -9386755640.150707 max: 6052647289138.035 average: 417880268.3301896 smoothSolver: Solving for k, Initial residual = 0.03655561646644971, Final residual = 0.001358430484242084, No Iterations 1 bounding k, min: -0.124148961509257 max: 13.54980413364046 average: 0.5762411249125505 drag : -219.7520678407049 lift : -36.38409258194373 ExecutionTime = 405.41 s ClockTime = 405 s Time = 0.004 smoothSolver: Solving for Ux, Initial residual = 0.3135682229449583, Final residual = 0.01065606675022723, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.2952628168401448, Final residual = 0.007776041656604311, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.2294808132606518, Final residual = 0.005939966102981783, No Iterations 2 GAMG: Solving for p, Initial residual = 0.01142465009285283, Final residual = 0.0001126659433648541, No Iterations 19 GAMG: Solving for p, Initial residual = 1.289434024409072e-13, Final residual = 1.289434024409072e-13, No Iterations 0 GAMG: Solving for p, Initial residual = 1.289434024409072e-13, Final residual = 1.289434024409072e-13, No Iterations 0 GAMG: Solving for p, Initial residual = 1.289434024409072e-13, Final residual = 1.289434024409072e-13, No Iterations 0 time step continuity errors : sum local = 0.001043798483872672, global = 4.972041248134074e-08, cumulative = 8.182394065470367e-08 smoothSolver: Solving for omega, Initial residual = 0.9964082542167786, Final residual = 0.068985186050248, No Iterations 6 bounding omega, min: -180036774609.3387 max: 6.329118058827089e+27 average: 1.089617791183109e+22 smoothSolver: Solving for k, Initial residual = 0.1882230893833461, Final residual = 5.367243120940624e-05, No Iterations 1 bounding k, min: -2.518099445090961 max: 95.21876859574435 average: 1.27445342716669 drag : -2869.298567197414 lift : -168.7879710025276 ExecutionTime = 429.36 s ClockTime = 429 s Time = 0.005 smoothSolver: Solving for Ux, Initial residual = 0.1448256781172823, Final residual = 0.003454284159427664, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.3306804680012086, Final residual = 0.006427471134016752, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.3215369458304032, Final residual = 0.006512543729490321, No Iterations 2 GAMG: Solving for p, Initial residual = 0.0002613290470658523, Final residual = 2.170007262345748e-06, No Iterations 5 GAMG: Solving for p, Initial residual = 8.867425281744662e-06, Final residual = 9.66013942785868e-08, No Iterations 12 GAMG: Solving for p, Initial residual = 4.241496950591247e-06, Final residual = 9.65159104422725e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 2.057181654485832e-06, Final residual = 9.279290081773349e-08, No Iterations 2 time step continuity errors : sum local = 1024.804129831097, global = -0.07002320644511717, cumulative = -0.07002312462117652 smoothSolver: Solving for omega, Initial residual = 0.5386999546680851, Final residual = 1.456360013262606e-09, No Iterations 1 bounding omega, min: -1.669470785195285e+25 max: 8.034849210266171e+32 average: 7.548360186268818e+25 smoothSolver: Solving for k, Initial residual = 0.07480861504188824, Final residual = 2.47095546180343e-08, No Iterations 1 bounding k, min: -15.45926119247786 max: 603.0051979285942 average: 2.288263708045382 drag : 7221162522.83777 lift : -5361883058.663867 ExecutionTime = 454.64 s ClockTime = 454 s Time = 0.006 smoothSolver: Solving for Ux, Initial residual = 0.3928042942650474, Final residual = 0.03693958221338808, No Iterations 17 smoothSolver: Solving for Uy, Initial residual = 0.3655420738309382, Final residual = 0.02841612700001563, No Iterations 18 smoothSolver: Solving for Uz, Initial residual = 0.3344446943495146, Final residual = 0.02705236061683774, No Iterations 17 GAMG: Solving for p, Initial residual = 0.9999876144347357, Final residual = 165.5982447656887, No Iterations 1000 GAMG: Solving for p, Initial residual = 1.744015098587018e-09, Final residual = 1.744015098587018e-09, No Iterations 0 GAMG: Solving for p, Initial residual = 1.744015098587018e-09, Final residual = 1.744015098587018e-09, No Iterations 0 GAMG: Solving for p, Initial residual = 1.744015098587018e-09, Final residual = 1.744015098587018e-09, No Iterations 0 time step continuity errors : sum local = 3263977569297249, global = -1136250716544.243, cumulative = -1136250716544.313 smoothSolver: Solving for omega, Initial residual = 0.5384603217155158, Final residual = 8.440472411405353e-08, No Iterations 1 bounding omega, min: -5.634206076641107e+31 max: 2.410454853511547e+32 average: 3.159597015085345e+25 smoothSolver: Solving for k, Initial residual = 0.1272320936674846, Final residual = 0.0104019729237369, No Iterations 4 bounding k, min: -129.7591077230691 max: 637.4025949508549 average: 2.099637272606892 drag : 1.097508704882627e+27 lift : -1.381264669290234e+27 ExecutionTime = 847.5599999999999 s ClockTime = 848 s Time = 0.007 -------------------------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpirun detected that one or more processes exited with non-zero status, thus causing the job to be terminated. The first process to do so was: Process name: [[59045,1],1] Exit code: 142 -------------------------------------------------------------------------- Thanks, Chris

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM.org] Error running #codeStream/functions wmake lnInclude fails to open *.H Files	tarik	OpenFOAM Installation	8	June 14, 2020 04:12
[foam-extend.org] Foam Extend 3.2 in ubuntu 16.04 LTS	JMosiezny	OpenFOAM Installation	15	December 9, 2017 09:20
OpenFOAM parallel on VM? pimpleDyMFoam propeller tutorial has been running for hours.	MBttR	OpenFOAM Running, Solving & CFD	0	June 13, 2016 11:58
Hypersonic Sim (running at Mach 7): floating point exception has occurred	bungusbeefcake	STAR-CCM+	10	March 31, 2015 07:59
Fluent 14.0 file not running in parallel mode in cluster	tejakalva	FLUENT	0	February 4, 2015 08:02

August 21, 2023, 04:36		#2
lowdo New Member Chris Join Date: Feb 2013 Posts: 17 Rep Power: 13	Is there anyone that can help with this?