[deal.II] Elasticity probelm: same matrices - different solution vector.

Dear all, I have been trying to implement orthotropic elasticity (similar as in step-18) and in order to verify my results I decided to set all E and G modulus values as in an isotropic case and to compare the outcome with the one for the same isotropic elastic problem but with another formulat

[deal.II] Interface finite elements

Dear All, I am investigating interface phenomena, taking place between porous absorbers and metal plates. I am trying to create a model that would account for adhesive, slip and roughness effects on such interfaces. In this respect I was wondering, whether some models with interface finite el

[deal.II] Coupled poroelastic-elastic problem

Dear All, I have been working with acoustic porous media and would like to extend my script to deal with two-layered acoustic absorbers, in this case I am interested in adding an elastic layer on top of my porous media. I took the step-46 as an example for matching two physical domains, everyt

Re: [deal.II] PETSc error 63

2022 at 5:47:08 PM UTC+1 Wolfgang Bangerth wrote: > On 2/4/22 04:42, Мария Бронзова wrote: > > > > Would you possibly have a hint on why that could be happening? > > This is how VTU files are supposed to look like. It's the format Visit and > Paraview understand. Wha

Re: [deal.II] PETSc error 63

e of the matrix. You need to find out why that is > so. > Run the program in a debugger and then inspect why the arguments you give > to > that function are out of range. > > Best > W. > > > On 2/3/22 07:03, Мария Бронзова wrote: > > *** Caution: EXTERNAL Sender

[deal.II] PETSc error 63

Deal all, I have been having troubles with running my script with PETSc. I have adjusted my code for parallel computing based on the example from step-40. But I am getting the following error: Exception on processing:

Re: [deal.II] PETScWrappers::MPI::BlockSparseMatrix

MakeFiles/Makefile2:271: CMakeFiles/run.dir/rule] Error 2 make: *** [Makefile:215: run] Error 2 среда, 19 января 2022 г. в 17:42:10 UTC+1, Wolfgang Bangerth: > On 1/19/22 03:41, Мария Бронзова wrote: > > I checked the member function documentation for the solve() function of &

Re: [deal.II] PETScWrappers::MPI::BlockSparseMatrix

Dear Wolfgang, Thank you a lot for the explanation, I got PETSc preinstalled already and didn't have a chance to see there were different compiling options. Now it's clear! Anyway, I ran into a further issue with my script that I wanted to parallelize. I am dealing with non spd matrices in ter

[deal.II] PETScWrappers::MPI::BlockSparseMatrix

Dear colleagues, I was trying to make adjustments in my script to enable parallel computing (as it takes unbearably long to calculate it even for small problems due to big sparsity matrix, I guess). In my script I am using block matrices and vectors with complex elements (i.e. BlockSparseMatr

Re: [deal.II] Dirichlet bc via apply_boundary_values function

Thank you a lot for the solution, good that this forum exists! Kind regards, Mariia четверг, 9 декабря 2021 г. в 04:44:01 UTC+1, Wolfgang Bangerth: > On 12/2/21 7:10 AM, Мария Бронзова wrote: > > > > Basically, now I only have one additional factor in my local_matrix >

Re: [deal.II] Dirichlet bc via apply_boundary_values function

u please give a hint, what is going on here? Thank you for your support! Here is the simplified code. Kind regards, Mariia понедельник, 25 октября 2021 г. в 22:21:31 UTC+2, Wolfgang Bangerth: > On 10/21/21 7:33 AM, Мария Бронзова wrote: > > > > I believe I answered you private

Re: [deal.II] Hinge boundary condition

12:45 PM, Мария Бронзова wrote: > > Oh, logic, thank you a lot! > > Would you like to suggest edits to the documentation of that function > that would have made that point clearer? Updating parts of the > documentation often make for excellent f

Re: [deal.II] Hinge boundary condition

Oh, logic, thank you a lot! среда, 3 ноября 2021 г. в 17:59:30 UTC+1, Wolfgang Bangerth: > On 11/3/21 10:00 AM, Мария Бронзова wrote: > > > > I was trying to prescribe displacement boundary conditions on a 3D model > > on certain faces but only for the x and y displ

[deal.II] Hinge boundary condition

Dear colleagues, I was trying to prescribe displacement boundary conditions on a 3D model on certain faces but only for the x and y displacement components. I used to apply the interpolate_boundary_values class before, but it demands dim+1 components inside my boundary_function (i have a mixe

[deal.II] Modelling of an infinite plate

Dear colleagues, I am trying to model an infinite plate in the context of the mixed formulation problem. To do so I need to set the x and y displacement component solution values to zero. I was looking for a function similar to apply_boundary_conditions, but for inner cells as well. Are there

Re: [deal.II] Dirichlet bc via apply_boundary_values function

<< "" << std::endl; std::cerr << "Exception on processing: " << std::endl << exc.what() << std::endl << "Aborting!"

Re: [deal.II] Dirichlet bc via apply_boundary_values function

;; SparseDirectUMFPACK A_direct; A_direct.initialize(system_matrix); A_direct.solve(system_matrix, system_rhs); solution = system_rhs; time.set_desired_next_step_size(1); } четверг, 14 октября 2021 г. в 11:15:44 UTC+2, Мария Бронзова: > Dear Wolfgang, > > Thank you for the valuabl

Re: [deal.II] Dirichlet bc via apply_boundary_values function

(0,0) (250,-3.92666e-06) (0,0) (0,0) (0,0) (0,0) (250,-5.66168) I am not quiet sure, what the problem might be. Would you possibly have an idea? Thank you a lot! Kind regards, Mariia среда, 13 октября 2021 г. в 17:35:12 UTC+2, Wolfgang Bangerth: > On 10/12/21 7:31 AM, Мария Бронзова wrote:

[deal.II] Dirichlet bc via apply_boundary_values function

Dear all, I"ve been trying to impose pressure boundary conditions in the context of a mixed displacement-pressure formulation problem, applying the "apply_boundary_values" function. I am dealing with complex numbers in my case, so I am imposing complex pressure values on certain faces with bou

Re: [deal.II] Dirichlet BC on block matrices

:49 AM, Мария Бронзова wrote: > > So, there are two boundary integrals in the formulation and I am trying > to > > implement them for the case of fixed displacement boundary condition. > The > > first boundary integral falls to zero for such a BC, as no displacement &

Re: [deal.II] Dirichlet BC on block matrices

*d_boundary_values* term? Or maybe even a smarter way to make it work? Please, let me know. Thank you a lot! четверг, 30 сентября 2021 г. в 14:28:26 UTC+2, Wolfgang Bangerth: > On 9/30/21 2:48 AM, Мария Бронзова wrote: > > > > Thank you a lot for the reply! I am trying to implement the Bi

Re: [deal.II] Dirichlet BC on block matrices

ay it"s done in step-20 example for the pressure boundary condition (so, in the local_rhs vector). Please, let me know, if it works so. Thank you for your support! With kind regards, Mariia Bronzova среда, 29 сентября 2021 г. в 19:23:16 UTC+2, Wolfgang Bangerth: > On 9/28/21 8

[deal.II] Dirichlet BC on block matrices

Deal all, I have an issue implementing Dirichlet boundary conditions in a mixed pressure-displacement problem. In step-20 there is an example for implementing pressure boundary conditions like that: *template * *class PressureBoundaryValues : public Function

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

lso tried fixing this within deal.II and discovered that some more > changes are necessary, see https://github.com/dealii/dealii/pull/12761. > Would you mind checking if the changes proposed there work for you? > > Best, > Daniel > > Am Fr., 17. Sept. 2021 um 06:25 Uhr sc

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

x_tools.cc#L240-L499 > > somewhere in your source code. > > Best, > Daniel > > Am Mi., 15. Sept. 2021 um 09:10 Uhr schrieb Мария Бронзова < > masia...@gmail.com>: > >> Dear Daniel, >> >> Thank you a lot for the hint, it makes p

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

std::map also uses values of type > std::complex. The compile error you shared shows that you use > doubles. > > Best, > Daniel > > Am Di., 14. Sept. 2021 um 07:30 Uhr schrieb Мария Бронзова < > masia...@gmail.com>: > >> Hi everyone, >> >> I ra

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

also uses values of type > std::complex. The compile error you shared shows that you use > doubles. > > Best, > Daniel > > Am Di., 14. Sept. 2021 um 07:30 Uhr schrieb Мария Бронзова < > masia...@gmail.com>: > >> Hi everyone, >> >> I ran into the

[deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

Hi everyone, I ran into the following difficulty: I have a vector-valued problem with complex numbers as entries for my BlockSparseMatrix system_matrix and BlockVector system_rhs and solution vectors and I am trying to apply Dirichlet boundary conditions, using the MatrixTools::apply_boundary_

[deal.II] Munmap_chunk(): invalid pointer

Hello everyone, I am trying to run the following script, but getting the error: 'munmap_chunk(): invalid pointer'. The eclipse debugger doesn't help me realize, where exactly the error occurs. In case you have an idea, how to correct it, please, let me know. I would be very gratefull! Kind re