[deal.II] Re: Step-42

Dear All, Thanks very much for your help. I really appreciate your kindness if you let me know why this step cannot be run as other tutorials? Bests, Ben On Monday, May 15, 2017 at 4:12:17 PM UTC-5, Bruno Turcksin wrote: > > Benhour, > > do the following three steps: > cmake . > make > ./step-42

[deal.II] Re: Step-42

Dear Jean-Paul, If I did correct, first I call "cmake ." and then "make run ./step-42 p1_adaptive.prm" and did not use "make". With these changes, it did not work again. Thanks, Ben On Monday, May 15, 2017 at 3:27:29 PM UTC-5, Jean-Paul Pelteret wrote: > > Dear Benhour, > > I think that you mus

[deal.II] Re: Step-42

Dear Daniel, According to your noteworthy comments, I tried exactly what you said, however, I got exactly the same as before, as what you got. Thanks, Benhour On Monday, May 15, 2017 at 12:23:44 PM UTC-5, Daniel Arndt wrote: > > Benhour, > > `make run` gives me: > > Scanning dependencies of targ

[deal.II] Re: Step-42

Dear All, Unfortunately I could not run this step. In fact, none of my friends also could run this step. I was wondering if you let me know it works in your system. Thanks, Benhour On Monday, May 15, 2017 at 6:25:51 AM UTC-5, Jean-Paul Pelteret wrote: > > Dear Benhour, > > I would suggest remo

[deal.II] Re: Step-42

Hi Daniel, I tried with version 8.4.1 and also the new version, 8.5, however it does not run at all. The errors that I got are given as below: make[3]: *** [CMakeFiles/run] Error 1 make[2]: *** [CMakeFiles/run.dir/all] Error 2 make[1]: *** [CMakeFiles/run.dir/rule] Error 2 make: *** [run] Error 2

[deal.II] Step-42

Dear All, I really appreciate it if you let me know how I can run step 42 because whenever I ran it, I got two error without any information about that. This is exactly the same with "goal oriented mesh adaptivity in elastoplasticity proiblems" example which is written by Shahram. Looking forwar

[deal.II] distributed body force

Dear All, It would be very kind of you if you let me know how I can imply distributed load on the inner radius of a ring instead of using nodal load. Thanks, Benhour -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/foru

[deal.II] Re: Neumann boundary condition for the curved domain

Dear Wolfgang, It should not be zero. In addition, the local quadrature point is for the cell. Thanks, Benhour On Thursday, December 8, 2016 at 8:16:03 PM UTC-6, benhour@gmail.com wrote: > > Dear All, > I still have problem with implying boundary condition on the curved > domain. consider

[deal.II] Neumann boundary condition for the curved domain

Dear All, I still have problem with implying boundary condition on the curved domain. consider we produce quarter of a ring in 2D with the function "quarter_hyper_shell" and want to impose Neumann boundary condition on the outer surface. In fact we have a force normal to this surface. The other

[deal.II] general question

Dear All, I hope you have had a great time. I have a general question. I define total strain and deviatoric part of the elastic strain tensor in two different ways as below: 1- Total strain tensor: const SymmetricTensor<2, dim> total_strain = 0.5*(Grad_u + transpose(Grad_u)); (Tensor<2, dim>

Re: [deal.II] Re: Apply Neumann Boundary condition

Dear Jean, Thanks very much for your help and support. Most of my problem is related to applying axisymmetry boundary condition for a curved domain. I attached my code for your consideration. According to the fact that this is my first project with dealii, Please accept my lack of knowledge in t

[deal.II] Re: Apply Neumann Boundary condition

Dear Jean, Thanks very much for your help and support. I do really appreciate your time and know you are really busy. The fact is that my code ran just for the first time step and can not solve for the other time step. when I refine the global mesh twice, this problem happened, however, with in

[deal.II] Re: Apply Neumann Boundary condition

Dear Jean, I know that this website is not the code reviewer. I want to apply Neumann boundary condition on the curved side of my domain. It should be note that the domain is symmetry. I have considered half of a circle. It would be very kind of you if you let me know whether it is applied corre

[deal.II] Apply Neumann Boundary condition

Dear All, I want to apply non-zero neumann boundary condition on the curved domain of half of the circle. The way that I define my geometry, the boundary id and apply the Neumann B.C is written as follow: *Defining boundary id for different sections of domain (half_hyper_ball):* const d

[deal.II] Re: Error with the content "dealii::numbers::is_finite(sum)"

Dear All, I am running in one process. What should i do for solving this error? Does it mean that my vector has complex component? On Tuesday, October 25, 2016 at 2:29:19 PM UTC-5, Daniel Arndt wrote: > > Benhour, > > are you running in serial or with more than one process? In the first > case,

[deal.II] Error with the content "dealii::numbers::is_finite(sum)"

Dear All, I am solving a thermoelastic problem in cylindrical coordinate system. I used CG solver for thermal problem and Newton-Raphson method for elasticity section. For solving linear system, I used direct solver SparseDirectUMFPACK with SSOR as a preconditioner. Unfortunately, I got some er

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Dear All, I want to solve my thermoelasticity problem in cylindrical coordinate systems. my domain includes two quarter of a circle in 2D and two quarter of two spheres in 3D. They are bulk. I attached my domain that was created by quarter_hyper_shell for 2D. But, I really do not kno how I can c

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Dear All, I still have some problem in solving my equations on a curved domain (sphere in 3D or circle in 2D). I got error with this content: invalid initialization of reference of type ‘dealii::Triangulation<2, 3>&’ from expression of type ‘dealii::Triangulation<3, 3>’ Does it mean that I should

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Yes exactly. The only difference is in weak form of my equation which should be multiplied by 2*pi*r which I define r_value = fe_values_eta.quadrature_point(q_point)[0]; In fact, I should solve my problem on a sphere (circle). According to your noteworthy comments, the only thing that I have to

Re: [deal.II] Re: Defining tensor in different dim and spacedim

Yes You are right. I mean curved domain. My question is should i use codim instead of mapping because of the fact that there is only two tutorials in dealii which solved a differential equation on c curved domain according to the fact that my derived weak expressions are in cylindrical coordinat

[deal.II] Re: Defining tensor in different dim and spacedim

Dear Jean-Paul, Let me put it in a such way. I want to solve my equation on a curved edge. According to dealii library, I should use MappingQ whenever I have nonlinear domain. Am I correct? In addition, Should I replace all dim in defining my tensors to spacedim? It should be noted that dim = sp

[deal.II] Re: Defining tensor in different dim and spacedim

Dear All, I want to solve a thermoelastic problem in cylindrical coordinate system. My domain is a circle in 2D or Sphere in 3D. I am following step 44 with considering the fact that I am solving for small strain. Because I have curved domain and I should Impose a boundary condition that change

[deal.II] Defining tensor in different dim and spacedim

Dear All, I am solving an elastic problem coupled with heat equation transfer (thermal-elasticity). I should solve the problem in a curved domain (circle in 2d and sphere in 3d) . I am developing step 44 in this case. I have to define some symmetric tensors. I have attached some parts of my code

[deal.II] Tensor with variable components

Dear All, I have a question. I want to define a tensor with components like the below expression, where eta is a scalar function in cylindrical coordinate system. I defined this component in dealii like this: FE_Qfe const QGauss quadrature_formula(3)

[deal.II] Derivative of a scalar function in the code

Dear All, I want to extract the first derivative and the gradient of a scalar function in dealii. Maybe I should use the void FEValuesViews::Scalar < dim, spacedim >::get_function_gradients. However, I really do not

[deal.II] Re: geometry

Daniel, Thanks very much for your help. For improving my knowledge and information, It would be very kind of you if you let me know why my code does not produce right geometry and which part is wrong. Sincerely, Benhour On Sunday, July 17, 2016 at 7:10:41 AM UTC-5, Daniel Arndt wrote: > > Benho

[deal.II] Re: geometry

Daniel, In that code, I had a half circle in the right part of the coordinate axis. By writing v(1) <=0, I wanted to choose the part with negative y components from the geometry, however, The attached file has been obtained. I really do not know how I should modify the source you sent to me beca

[deal.II] Re: geometry

Dear J-P, Thanks for your response. I used the code for creating a quarter of a circle that comes as follow: Triangulation<2> triangulation; const Point<2> center; const double radius = 1.; GridGenerator::half_hyper_ball(triangulation, center, radius); Triangulation<2>::active_cell_iterator cell

[deal.II] axis-symmetric problem

Dear All, Is it possible to define a problem in spherical or cylindrical coordinate system instead of Cartesian coordinate system?It would be very kind of you if you refer me to a tutorial in which there is a problem that has been solved in spherical or cylindrical coordinate system. Best, Benh

[deal.II] geometry

Dear All, I want to create a quarter of a circle in deal ii. In this case, I could create half of a circle with half_hyper_ball. I have tried GridTools::delete_unused_vertices or GridTools::transformation to remove or shift the same vertices to eliminate half of the circle. It would be very ki

[deal.II] Re: Geometry and boundary conditions

Dear Daniel, Thanks very much for your response. I should model a whole circle, however for simplicity I want to model one quarter of it. In fact I have 2 quarter of circles with different radius and with the same center. I have a boundary load on the perimeter(curve) side of the circle, axial s

[deal.II] Geometry and boundary conditions

Dear All, I want to define a new geometry (a quarter of a circle) and apply symmetry boundary condition on two sides of the quarter. It would be very kind of you if you help me for that. Best, Benhour -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, s

[deal.II] Re: Elasticity

Dear all, Thanks very much for your response and help. If I can make it simple, in my equation, the lame constants in stress-strain tensor are related to another function that changes with solution vector. For instance, consider a situation in which you have a constant rank-2 tensor for initial

[deal.II] Elasticity

Dear all, In step 18, the elasticity problem has been solved with this assumption that there was not initial strain in the equation. I want to add the initial strain with consideration of the fact that this tensor is a function of the solution. In other words, this tensor should be multiplied