[deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-17 Thread Sclah
Hello everyone, I am using step-46 as starting point to implement a time dependent fsi problem using block vectors and block matrices. I successfully written the time iteration and changed everything to block structures but I am struggling to find a method to impose the fluid velocity equal to

Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-18 Thread Sclah
le 04:27:19 UTC+1 Wolfgang Bangerth ha scritto: > On 12/16/24 16:12, Sclah wrote: > > > > I am using step-46 as starting point to implement a time dependent fsi > problem > > using block vectors and block matrices. > > I successfully written the time iteration and

Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-18 Thread Sclah
fgang Bangerth ha scritto: > On 12/16/24 16:12, Sclah wrote: > > > > I am using step-46 as starting point to implement a time dependent fsi > problem > > using block vectors and block matrices. > > I successfully written the time iteration and changed everything t

Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-19 Thread Sclah
fgang Bangerth ha scritto: > On 12/18/24 16:00, Sclah wrote: > > > > for (unsigned int i = 0; i < local_face_dof_indices_fluid.size(); ++i) > > { > > if (fe_fluid->face_system_to_component_index(i).first < dim) // select > just > > fluid v

Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2025-01-10 Thread Sclah
UTC+1 Wolfgang Bangerth ha scritto: > On 12/19/24 11:52, Sclah wrote: > > Thank you for your help.Yes, I am using the same finite element for the > two > > variables. > > I just checked and face_component_to_system_index() doesn't seem to > exist > > inside

[deal.II] Compute the divergence of a tensor

2025-04-04 Thread Sclah
Hello everyone, I am trying to implement a continuous mechanics simulation and I am encountering troubles trying to compute the divergence of the deformation gradient inverse. I need it to assemble the variational form terms. My code reads something like the following: - take the gradient of the

Re: [deal.II] Compute the divergence of a tensor

2025-04-04 Thread Sclah
(u)^T)] Finv* Then I guess that using get_function_hessians() on *u* is what I need to assemble correctly my weak form. Thank you again! Il giorno venerdì 4 aprile 2025 alle 17:49:05 UTC+2 Wolfgang Bangerth ha scritto: > On 4/4/25 09:24, Sclah wrote: > > > > - take the

Re: [deal.II] Compute the divergence of a tensor

2025-04-04 Thread Sclah
Yes, thank you for the warning, that could be an issue indeed. In my case I am using P2 elements, I'll check the best way to approach this, thank you! Il giorno venerdì 4 aprile 2025 alle 18:29:26 UTC+2 Wolfgang Bangerth ha scritto: > On 4/4/25 10:25, Sclah wrote: > > ** >