atenamsv

Thanks for your tips, but unfortunately the error does not vanish. here are the required files. [createFields.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911102/createFields.txt) [pa.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911145/pa.txt) [T.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911141/T.txt) [U.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911142/U.txt) [Ua.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911143/Ua.txt) [p.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911144/p.txt) [PadjointShapeOptimizationFoam.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4911149/PadjointShapeOptimizationFoam.txt)

Here i could see in each loop pressure is calculating 2 time, that's where density goes minus and causes crash. [errorLog2.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4912471/errorLog2.txt)

I want to use adjointShapeOptimization for a compressible flow. [adjointShapeOptimizationFoam.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4919314/adjointShapeOptimizationFoam.txt) [createFields.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4919316/createFields.txt) The default solver in openFoam is for incompressible flow, and uses SIMPLE algorithm. I tried to remove simple solver...

[PimpleAdjointShapeOptimizationFoam.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4923937/PimpleAdjointShapeOptimizationFoam.txt) [hEqn.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4923949/hEqn.txt) [pEqn.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4923950/pEqn.txt) [UEqn.txt](https://github.com/unicfdlab/hybridCentralSolvers/files/4923951/UEqn.txt) I tried to alter the solver, pEqn, UEqn and hEqn. now it is working, and goes some steps forward, but now the problem is that the...