Yidali26

If the elasticity is removed, the result agree with the analytical solution. I wonder if the elasticity could be compatible with the Neumann BC. Thanks

Hi all, Is there any update or suggestion on this issue? I hope to set up some kind of open boundary with a visco-elastic problem, and this test is quite...

> There is something of a difficulty with non-linear stress boundary conditions / stress boundary conditions in which the element rheology plays a role because the usual integration by parts...

   

In the code appliedTractionField.data[mesh.specialSets["MaxK_VertexSet"].data,0] = time-eta_eff*previousStress.data[mesh.specialSets["MaxK_VertexSet"].data,4]/(miu*dt_e) should be replaced with appliedTractionField.data[mesh.specialSets["MaxK_VertexSet"].data,0] = time

> Hi @Yidali26, sorry for the late reply. > > > In the code > > appliedTractionField.data[mesh.specialSets["MaxK_VertexSet"].data,0] = time-eta_eff_previousStress.data[mesh.specialSets["MaxK_VertexSet"].data,4]/(miu_dt_e) > > should be replaced with > > appliedTractionField.data[mesh.specialSets["MaxK_VertexSet"].data,0] = time...

> Hi @Yidali26, > Thanks for the 2D version. Some thoughts on this: > > 1. The initial 1st solve is solving for 0 external force, so our stokes solvers...

Hi @naliboff, Thanks for your reply. Do I need to add trilinos back to the .cfg file? Previously I removed trilino in the local.cfg file since its installation will lead...

Hi @naliboff , Thanks, I tried so but still doesn't work. After removing once:trilinos candi doesn't install trilinos at all, and subsequent cmake .. for aspect fail with the same...

> @Yidali26 - were you able to find a solution for compiling Trilinos? If not, one possibility would be to manually specify `TRILINOS_DIR` in [dealii.package](https://github.com/dealii/candi/blob/master/deal.II-toolchain/packages/dealii.package#L171)? @naliboff Not yet. May I...