Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Wasim, Why do you have two different triangulations? Can you have single Triangulation and use FE_Nothing ? For example Step-70 solves two different equations on different parts of the domain using a single Triangulation. Best, Bruno Le lun. 21 nov. 2022 à 15:27, Wasim Niyaz Munshi ce21d400 <

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Thank you, Prof. Bruno. I used the above line of code and realized where I was going wrong. When solving the damage equation, I also call a function (let's call it H). This function is needed to find the element stiffness matrix of the damage equation. This function returns some value for a

Re: [deal.II] Solve block Trilinos matrix directly

On 11/21/22 07:45, 'yy.wayne' via deal.II User Group wrote: However, by copy to and from between LinearAlgebra::distributed::BlockVector and LinearAlgebra::distributed::Vector, it produces wrong result. Though I cannot confirm the wrong result comes from transfer between Vector and

Re: [deal.II] Implementation of Neumann boundary condition in 3d

On 11/21/22 03:04, HIMAL MAGAR wrote: 0-02-03-ecf50731d20a4fed58e332ff9a1702a1e061032d7400f34e6ee11a2b55bd6aea_2073f8a6123d5b.png Scale was the same. I am extremely happy about the result :-). Very nice! For now, I want to output norm_stress in the results. In step-18, it is said that the

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Wasim, Before the line that errors out can you do std::cout I am working with 2 DoFHandlers. I checked my code. I am using the correct > DoFHandlers for the 2 equations. > > On Mon, Nov 21, 2022 at 6:58 PM Bruno Turcksin > wrote: > >> Wasim, >> >> It's hard to

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

I am working with 2 DoFHandlers. I checked my code. I am using the correct DoFHandlers for the 2 equations. On Mon, Nov 21, 2022 at 6:58 PM Bruno Turcksin wrote: > Wasim, > > It's hard to say without seeing any code. Are you working with one or two > DoFHandler? If you are working with two

[deal.II] Solve block Trilinos matrix directly

I've post a similar question before(Solve distributed BlockVector with direct solver ), where I had a Trilinos SparseMatrix(the system has 2 components) and

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Wasim, In case, what you do is close to the phase-field fracture model (which is essentially a gradient damage model, and in this case a one-way coupled one), you might find the following master thesis and the associated code helpful: https://github.com/Jatandeep/Thesis Best regards, Paras On

[deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Wasim, It's hard to say without seeing any code. Are you working with one or two DoFHandler? If you are working with two DoFHandler, are you sure that you are using the correct one? Best, Bruno On Monday, November 21, 2022 at 5:54:22 AM UTC-5 ce21...@smail.iitm.ac.in wrote: > Hello

[deal.II] Solving elasticity and laplace(damage) equation, one after the other.

Hello everyone. I am trying to solve elasticity and *laplace (damage)* equations in the same program. The rhs of laplace equation depends on the elasticity solution vector. Both equations have already been solved in step3 and 8. I am following the same approach. I am creating my objects (like