Re: [deal.II] Re: interpolate different solutions on different parts of the domain

2018-08-28 Thread Bruno Turcksin
Jane, Le mar. 28 août 2018 à 07:50, Jane Lee a écrit : > That is correct. How would i use fe_values.get_function_values to > interpolate though? Since I have to get the values off of domain 1 only (as > my domain 2 has been added to the domain 1, and I only have the solution f > - or as you

[deal.II] Design question about multiscale FEM with deal.ii

2018-08-28 Thread Konrad
Hi deal.ii community, I am having a little issue with a (smart) design of a (not so simple) multiscale FEM. The problem in its simple form is the following: I have a coarse mesh and on each coarse cell I replace the standard basis (for example Q1) with basis functions that reflect some

Re: [deal.II] Adding to / changing domain at every timestep. merging mesh.

2018-08-28 Thread Jane Lee
Thank you so much. I will have a look! On Wednesday, August 1, 2018 at 1:58:27 PM UTC+1, Jean-Paul Pelteret wrote: > > Dear Jane, > > The method described by Claire is exactly the approach used by my > colleagues to model additive manufacturing processes. See, for example, > this paper >

Re: [deal.II] Error in running Quasi_Static Compressible example on Ubuntu

2018-08-28 Thread Lucky Qin
Dear Dr. Bangerth, Thank you for the quick reply! After installing LAPACK, the code works! Thanks, Qin 在 2018年8月28日星期二 UTC+8上午4:54:52,Wolfgang Bangerth写道: > > On 08/27/2018 08:29 AM, Lucky Qin wrote: > > > > I am trying to run the example in code gallery named > >

Re: [deal.II] interpolate different solutions on different parts of the domain

2018-08-28 Thread Jane Lee
Thanks for this suggestion - I will have a go at this as well. On Tuesday, July 31, 2018 at 8:05:11 PM UTC+1, Jean-Paul Pelteret wrote: > > Hi Jane, > > I’m not quite sure that I understand exactly what you’re trying to do, but > as an additional consideration to Bruno’s suggestions is the >

Re: [deal.II] Re: interpolate different solutions on different parts of the domain

2018-08-28 Thread Jane Lee
HI Bruno, That is correct. How would i use fe_values.get_function_values to interpolate though? Since I have to get the values off of domain 1 only (as my domain 2 has been added to the domain 1, and I only have the solution f - or as you call alpha - in a domain that consists of only domain

Re: [deal.II] How to give dof per each node?

2018-08-28 Thread chandra sekhar yaswanth devarakonda
thanks a lot. I got it now. best, chandra On Mon, Aug 6, 2018 at 11:51 PM Jean-Paul Pelteret wrote: > Dear Chandra, > > Its a little unclear to me as to what you mean when you say > > to give dof per each element/node > > > Do you mean that you assign DoFs to the mesh? If so, then what

[deal.II] step-8

2018-08-28 Thread chandra sekhar yaswanth devarakonda
Where can i get theory behind step-8? Can anyone suggest any book? -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups

Re: [deal.II] step-8

2018-08-28 Thread Jean-Paul Pelteret
Dear Chandra, For me, one of the most instructive textbooks on linear elasticity (and that fundamentals of FEM) is this one: @Book{Hughes2000a, title = {The Finite Element Method: Linear Static and Dynamic Finite Element Analysis}, publisher = {Dover Publications Inc.}, year =

[deal.II] Inhomogeneous neumann bc for step-20

2018-08-28 Thread Jane Lee
Dear all, I am trying to solve the equations in step-20 with inhomogeneous neumann bcs on one of the boundaries and getting something very bizarre. I have a rectangular domain with the following: 1. Top boundary has homogeneous conditions: this is applied into the weak form like in step-20

Re: [deal.II] step-8

2018-08-28 Thread chandra sekhar yaswanth devarakonda
Thanks a lot. I will go through it. On Tue, Aug 28, 2018, 7:15 PM Jean-Paul Pelteret wrote: > Dear Chandra, > > For me, one of the most instructive textbooks on linear elasticity (and > that fundamentals of FEM) is this one: > > @Book{Hughes2000a, > title = {The Finite Element Method:

[deal.II] step-8

2018-08-28 Thread chandra sekhar yaswanth devarakonda
I am using step-8 to build 2d cantilever beam.Can anyone help me whether it is correct or not? actually i got error.i am new to finite element and dealii as well. for elemental stiffness matrix=k=integration[E.I B.transpose(B)] .. so i did like this cell_matrix(i,j) +=

Re: [deal.II] How to give dof per each node?

2018-08-28 Thread chandra sekhar yaswanth devarakonda
thanks Rakesh, thats really helped best chandra On Fri, Aug 3, 2018 at 6:04 PM Rakesh wrote: > Sry, I meant Stiffness matrix* not Stiffness element(in Step-3) > > > On Fri, Aug 3, 2018 at 8:26 AM, Rakesh > wrote: > >> Hello Yaswanth, >> >> To create a dof for a particular mesh, these are

[deal.II] Re: step-8

2018-08-28 Thread Loylick
Essentially There is no need in changing anything in assemble_system function from step-8. At least for stiffness matrix, you should change boundary conditions and triangulation to suit for your task, that is all. However you should not mix up Young's modulus with Lame parameters which are

[deal.II] Re: Inhomogeneous neumann bc for step-20

2018-08-28 Thread Loylick
ingomogeneous neumann BC are imposed through including them into a stiffness matrix. You should not use MatrixTools::apply_boundary_values for this purpose I guess. Look the Step-7 trutorial for the info on how to calculate boundary integral for neumann BC. -- The deal.II project is

[deal.II] Re: step-8

2018-08-28 Thread Loylick
Essentially, there is no need in changing anything in assemble_system function from step-8, at least for stiffness matrix. You should change boundary conditions and triangulation to suit for your task, that is all. However you should not mix up Young's modulus with Lame parameters which are

Re: [deal.II] Inhomogeneous neumann bc for step-20

2018-08-28 Thread Wolfgang Bangerth
On 08/28/2018 08:07 AM, Jane Lee wrote: I am trying to solve the equations in step-20 with inhomogeneous neumann bcs on one of the boundaries and getting something very bizarre. step-20 uses a mixed formulation in which both the pressure and the velocity (in essence, the gradient of the

Re: [deal.II] how to remove rows from system matrix

2018-08-28 Thread Wolfgang Bangerth
On 08/27/2018 06:58 PM, Loylick wrote: Thank you, Wolfgang, for your suggestion. I've read step-11 tutorial but I still can not figure out what kind of constraints I can impose in my case. What does the literature suggest how to handle this situation? Here is a picture of a void with

Re: [deal.II] Re: step-8

2018-08-28 Thread chandra sekhar yaswanth devarakonda
Thank you.Thats really helped On Wed, Aug 29, 2018, 1:18 AM Loylick wrote: > Essentially, there is no need in changing anything in assemble_system > function from step-8, > at least for stiffness matrix. You should change boundary conditions and > triangulation to suit for your task, > that is

Re: [deal.II] Design question about multiscale FEM with deal.ii

2018-08-28 Thread Wolfgang Bangerth
On 08/28/2018 04:47 AM, Konrad wrote: Is it possible to have system matrices, dofs and solution vectors etc on two different refinement levels of a mesh? Say my coarse FEM holds data on refinement level 5 and my basis function(s) on each coarse cell live on the leaves of each coarse cell on

[deal.II] Re: Publications based on deal.II

2018-08-28 Thread Alberto Salvadori
Dear Wolfgang, we just published the paper: A. Salvadori Et Al, "Modeling and Simulation of VEGF Receptors Recruitment in Angiogenesis," Mathematical Problems in Engineering, vol. 2018, Article ID 4705472, 10 pages, 2018. https://doi.org/10.1155/2018/4705472/. using and quoting deal.ii .

Re: [deal.II] projecting values of the solution on the boundary

2018-08-28 Thread Alberto Salvadori
Thank you JP, your hints are highly appreciated. I will work on creating my own data structure, to see how complex and effective it can be. In the meanwhile though I made some tests with "almost" zero volumes cohesive interfaces and they are indeed simple and effective. So, if there is a chance