Re: [deal.II] ABOUT ASSEMBLING RIGHT HAND SIDE

[孙翔]

Hi, Wolfgang,

I applied a Dirichlet boundary condition after I formed system matrixes and 
vectors, then use the method you taught in the lecture to change the 
entries in the system matrix and rhs vector. Could I upload the code here 
in order to make it clear? Thank you very much.

Best,

Xiang 


On Tuesday, 7 April 2020 16:03:10 UTC-7, Wolfgang Bangerth wrote:
>
> On 4/7/20 4:58 PM, 孙翔 wrote: 
> > 
> > Do you mean the Dirichlet boundary condition? I applied the boundary 
> condition 
> > after I form the right-hand side vector. I don't think the boundary 
> condition 
> > will lead to different solutions because the right-hand side vector has 
> been 
> > different before I add the boundary conditions. 
>
> "I don't think" is the wrong approach :-) You are seeing something you are 
> not 
> expecting. In other words, what you "expect" is not happening, and so you 
> need 
> to let go of what you think is happening -- you need to *verify* that what 
> you 
> think is happening is happening! That's the whole mind set of debugging 
> problems :-) 
>
> Best 
>   WB 
>
>
> -- 
>  
> Wolfgang Bangerth  email: bang...@colostate.edu 
>  
> www: http://www.math.colostate.edu/~bangerth/ 
>
>

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/13771a9b-4ee7-4a7e-91b7-e763ab0917a9%40googlegroups.com.

Re: [deal.II] ABOUT ASSEMBLING RIGHT HAND SIDE

[孙翔]

Hi, Wolfgang,

Thanks a lot. I didn't apply refinement in my code, so there is no hanging 
node problem. As you said, there is a big difference around the boundary. I 
set the boundary value always equal to the initial value which is equal to 
1. The M*old_solution method gives the correct answer which is 1, but the 
second method can not. The attached file is the solution. 

Best,

Xiang

On Tuesday, 7 April 2020 16:03:10 UTC-7, Wolfgang Bangerth wrote:
>
> On 4/7/20 4:58 PM, 孙翔 wrote: 
> > 
> > Do you mean the Dirichlet boundary condition? I applied the boundary 
> condition 
> > after I form the right-hand side vector. I don't think the boundary 
> condition 
> > will lead to different solutions because the right-hand side vector has 
> been 
> > different before I add the boundary conditions. 
>
> "I don't think" is the wrong approach :-) You are seeing something you are 
> not 
> expecting. In other words, what you "expect" is not happening, and so you 
> need 
> to let go of what you think is happening -- you need to *verify* that what 
> you 
> think is happening is happening! That's the whole mind set of debugging 
> problems :-) 
>
> Best 
>   WB 
>
>
> -- 
>  
> Wolfgang Bangerth  email: bang...@colostate.edu 
>  
> www: http://www.math.colostate.edu/~bangerth/ 
>
>

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/89db432f-1bbc-436a-a14c-c79b71d887a1%40googlegroups.com.

Re: [deal.II] ABOUT ASSEMBLING RIGHT HAND SIDE

[Wolfgang Bangerth]

On 4/7/20 4:58 PM, 孙翔 wrote:


Do you mean the Dirichlet boundary condition? I applied the boundary condition 
after I form the right-hand side vector. I don't think the boundary condition 
will lead to different solutions because the right-hand side vector has been 
different before I add the boundary conditions.


"I don't think" is the wrong approach :-) You are seeing something you are not 
expecting. In other words, what you "expect" is not happening, and so you need 
to let go of what you think is happening -- you need to *verify* that what you 
think is happening is happening! That's the whole mind set of debugging 
problems :-)


Best
 WB


--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/a5dac54b-9e2a-1f12-db7b-95c52babfd1a%40colostate.edu.

Re: [deal.II] ABOUT ASSEMBLING RIGHT HAND SIDE

[孙翔]

Hi, Jean-Paul,

Do you mean the Dirichlet boundary condition? I applied the boundary 
condition after I form the right-hand side vector. I don't think the 
boundary condition will lead to different solutions because the right-hand 
side vector has been different before I add the boundary conditions. 

Best,

Xiang

On Tuesday, 7 April 2020 13:27:17 UTC-7, Jean-Paul Pelteret wrote:
>
> Hi Xiang, 
>
> You’ve also not mentioned whether you’re comparing vectors before or after 
> you’ve distributed the constraints. Perhaps you get the same result if you 
> call constraints.distribute() on both of your vectors? 
>
> Best, 
> Jean-Paul 
>
> > On 07 Apr 2020, at 22:22, Wolfgang Bangerth  > wrote: 
> > 
> > On 4/5/20 12:32 AM, Xiang Sun wrote: 
> >> For correction, the right code for the second method is 
> > > [...] 
> > 
> > Are you saying that you solved your own problem? Or are you still 
> unclear what is happening? 
> > 
> > Best 
> > W. 
> > 
> > 
> > -- 
> >  
> > Wolfgang Bangerth  email: bang...@colostate.edu 
>  
> >   www: http://www.math.colostate.edu/~bangerth/ 
> > 
> > -- 
> > 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 "deal.II User Group" group. 
> > To unsubscribe from this group and stop receiving emails from it, send 
> an email to dea...@googlegroups.com . 
> > To view this discussion on the web visit 
> https://groups.google.com/d/msgid/dealii/77b269f1-b91b-1dc3-2145-b023242ba43c%40colostate.edu.
>  
>
>
>

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/dbbae2fe-3a1d-4120-9922-c24463ff5e93%40googlegroups.com.

Re: [deal.II] Re: ABOUT ASSEMBLING RIGHT HAND SIDE

[Wolfgang Bangerth]

On 4/7/20 4:53 PM, 孙翔 wrote:


I am still unclear about what is happening. That zero is just a typo. Even if 
I remove the zero, these two methods can not get the same solution. What I did 
is just using two methods to calculate the right-hand side generated by time 
discretization. One way is to form a mass matrix then time the old solution, 
the other way is to form the local right-hand side in element level then 
assemble it into system right-hand side array.


Right, these two approaches should result in the same solution. So the 
question is: how do they differ? Do they result in different right hand side 
vectors, or in different solutions? If they differ in the solutions, once you 
visualize them, is it obvious that the problem is at the boundary? At hanging 
nodes? These are useful tools to figure out where the bug may be.


Best
 W.


--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/4f4e2b7f-11bc-7d89-69d4-938d672b20f6%40colostate.edu.

Re: [deal.II] Re: ABOUT ASSEMBLING RIGHT HAND SIDE

[孙翔]

Hi, Wolfgang,

I am still unclear about what is happening. That zero is just a typo. Even 
if I remove the zero, these two methods can not get the same solution. What 
I did is just using two methods to calculate the right-hand side generated 
by time discretization. One way is to form a mass matrix then time the old 
solution, the other way is to form the local right-hand side in element 
level then assemble it into system right-hand side array.  

Best,

Xiang

On Tuesday, 7 April 2020 13:22:55 UTC-7, Wolfgang Bangerth wrote:
>
> On 4/5/20 12:32 AM, Xiang Sun wrote: 
> > For correction, the right code for the second method is 
>  > [...] 
>
> Are you saying that you solved your own problem? Or are you still unclear 
> what 
> is happening? 
>
> Best 
>   W. 
>
>
> -- 
>  
> Wolfgang Bangerth  email: bang...@colostate.edu 
>  
> www: http://www.math.colostate.edu/~bangerth/ 
>
>

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/8c36b7a6-7fd0-4976-8aea-d6a52ab28670%40googlegroups.com.

Re: [deal.II] ABOUT ASSEMBLING RIGHT HAND SIDE

[Jean-Paul Pelteret]

Hi Xiang,

You’ve also not mentioned whether you’re comparing vectors before or after 
you’ve distributed the constraints. Perhaps you get the same result if you call 
constraints.distribute() on both of your vectors?

Best,
Jean-Paul

> On 07 Apr 2020, at 22:22, Wolfgang Bangerth  wrote:
> 
> On 4/5/20 12:32 AM, Xiang Sun wrote:
>> For correction, the right code for the second method is
> > [...]
> 
> Are you saying that you solved your own problem? Or are you still unclear 
> what is happening?
> 
> Best
> W.
> 
> 
> -- 
> 
> Wolfgang Bangerth  email: bange...@colostate.edu
>   www: http://www.math.colostate.edu/~bangerth/
> 
> -- 
> 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 
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit 
> https://groups.google.com/d/msgid/dealii/77b269f1-b91b-1dc3-2145-b023242ba43c%40colostate.edu.

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/E8F4F48E-9EE1-48B1-A301-F62441E0698B%40gmail.com.

Re: [deal.II] Re: ABOUT ASSEMBLING RIGHT HAND SIDE

[Wolfgang Bangerth]

On 4/5/20 12:32 AM, Xiang Sun wrote:

For correction, the right code for the second method is

> [...]

Are you saying that you solved your own problem? Or are you still unclear what 
is happening?


Best
 W.


--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/77b269f1-b91b-1dc3-2145-b023242ba43c%40colostate.edu.

[deal.II] deal.II Newsletter #117

[Rene Gassmoeller]

Hello everyone!

This is deal.II newsletter #117.
It automatically reports recently merged features and discussions about the 
deal.II finite element library.


## Below you find a list of recently proposed or merged features:

#9851: do not build release on github actions OSX (proposed by tjhei) 
https://github.com/dealii/dealii/pull/9851

#9850: Document the mental model of DiscreteTime. (proposed by bangerth) 
https://github.com/dealii/dealii/pull/9850

#9849: Fix some warnings in tests with clang-10 (proposed by masterleinad) 
https://github.com/dealii/dealii/pull/9849

#9848: [CI]: add GitHub actions linux release build (proposed by tjhei) 
https://github.com/dealii/dealii/pull/9848

#9847: Rename NoncontiguousPartitioner::update_values() [WIP] (proposed by 
peterrum) https://github.com/dealii/dealii/pull/9847

#9846: Update documentation of NoncontiguousPartitioner (proposed by 
kronbichler; merged) https://github.com/dealii/dealii/pull/9846

#9845: Fix some Doxygen warnings (proposed by rezarastak; merged) 
https://github.com/dealii/dealii/pull/9845

#9844: Build debug and release with Actions (proposed by masterleinad; merged) 
https://github.com/dealii/dealii/pull/9844

#9843: Do not verify MPI_HAVE_MPI_SEEK_SET during cross compiling (proposed by 
dangars) https://github.com/dealii/dealii/pull/9843

#9842: Output step-62 results as a complex-valued vector (proposed by dangars; 
merged) https://github.com/dealii/dealii/pull/9842

#9841: Use more range-based for loops. (proposed by masterleinad; merged) 
https://github.com/dealii/dealii/pull/9841

#9840: Move implementation of DoFHandler and hp::DoFHandler (proposed by 
peterrum) https://github.com/dealii/dealii/pull/9840

#9838: step-47 - fix typos (proposed by jwitte08; merged) 
https://github.com/dealii/dealii/pull/9838

#9837: Update Epetra link (proposed by peterrum; merged) 
https://github.com/dealii/dealii/pull/9837

#9836: Make TpetraWrappers appear in Doxygen (proposed by peterrum; merged) 
https://github.com/dealii/dealii/pull/9836

#9835: Parallelize the main loops of MappingInfo::compute_mapping_q (proposed 
by kronbichler) https://github.com/dealii/dealii/pull/9835

#9834: Step-63: Fix two typos (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/9834

#9833: fixes issue 9722 by correctly writing time data into a aggregate pvd ... 
(proposed by krishnakumarg1984) https://github.com/dealii/dealii/pull/9833

#9832: This is the actual fix for issue #9722. Not the earlier PR #9754 
(proposed by krishnakumarg1984) https://github.com/dealii/dealii/pull/9832

#9831: Implement functionality from Hypercube on BoundingBox (proposed by 
simonsticko; merged) https://github.com/dealii/dealii/pull/9831

#9830: Doc of Tensor and SymmetricTensor improved. (proposed by rezarastak) 
https://github.com/dealii/dealii/pull/9830

#9829: FETools: Add additional lexicographic_to_hierarchic_numbering functions 
(proposed by kronbichler) https://github.com/dealii/dealii/pull/9829

#9828: Wrapper for SuperLU-MT linear solver (proposed by KParas) 
https://github.com/dealii/dealii/pull/9828

#9827: Fix bug in VectorTools::get_position_vector for high order mappings 
(proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/9827

#9826: Do not use FE_Q in MappingQGeneric (proposed by kronbichler) 
https://github.com/dealii/dealii/pull/9826

#9825: Fix warning in Utilities::MPI for serial release build (proposed by 
peterrum; merged) https://github.com/dealii/dealii/pull/9825

#9822: Place figures for 2d/3d quarter hyper ball in same line (proposed by 
kronbichler; merged) https://github.com/dealii/dealii/pull/9822

#9821: Remove deprecated class PArpackSolver::Shift (proposed by masterleinad; 
merged) https://github.com/dealii/dealii/pull/9821

#9819: Add Utilities::MPI::compute_set_union (proposed by peterrum; merged) 
https://github.com/dealii/dealii/pull/9819

#9818: A few doxygen-related fixes. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/9818

#9817: Require gcc 4.9 and clang 4.0. (proposed by bangerth) 
https://github.com/dealii/dealii/pull/9817

#9816: more boost warning guards (proposed by tjhei; merged) 
https://github.com/dealii/dealii/pull/9816

#9815: Allow outputting complex-valued vector and tensor fields (proposed by 
bangerth; merged) https://github.com/dealii/dealii/pull/9815

#9814: Add a missing assertion. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/9814

#9813: Simplify some code by using a range-based for loop. (proposed by 
bangerth; merged) https://github.com/dealii/dealii/pull/9813

#9812: Fix get_boundary_ids() and get_manifold_ids() for parallel 
triangulations (proposed by peterrum) https://github.com/dealii/dealii/pull/9812

#9811: Simplify instantiation of DoFHandler and hp::DoFHandler (proposed by 
peterrum; merged) https://github.com/dealii/dealii/pull/9811

#9810: Add functions accepting hp::DoFHandler to BlockInfo (proposed by 
peterrum; merged) https:/

Re: [deal.II] Polycrystal on deal ii

[Wolfgang Bangerth]

On 4/7/20 11:59 AM, A.Z Ihsan wrote:


Yes, thats what i want, just simply put the elastic anisotropy at the Gauss 
point level.


step-6 already does that: It uses different material parameters depending on 
the location of the quadrature point.


Best
 W.

--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/976f3e9d-a795-de54-c2b7-0b1e8ad32cc9%40colostate.edu.

Re: [deal.II] Polycrystal on deal ii

[A.Z Ihsan]

Hi Andrew, 

Yes, thats what i want, just simply put the elastic anisotropy at the Gauss 
point level. 
Thank you for the answer. 

Best, 
Ihsan

On Tuesday, April 7, 2020 at 7:41:42 PM UTC+2, mac wrote:
>
> In the classical continuum approach to plasticity there is no grain 
> boundary. The elastic anisotropy of the crystal is captured at the Gauss 
> point level. So you would start by associating different materials 
> properties (elastic / plastic) to Gauss points depending on their spatial 
> location. 
>
> Have a look at the CellDataStorage class and how it is used in step-44.
>
> Best
> Andrew
>
> On 7 Apr 2020, at 18:36, A.Z Ihsan > 
> wrote:
>
> Hi Wolfgang, 
>
> alright, suppose i have a box that divided into two regions, and each 
> region has different elastic constant (anisotropy material). 
> Perhaps the idea is when calculating the cell_matrix, inside the for-loop 
> there is a function for choosing in which elastic constant to be used 
> depending on the region. 
>
> I already have the elastic solver written specifically for one elastic 
> constant.  
> Do you have any hint where should i start? 
>
> BR, 
> Ihsan
> On Tuesday, April 7, 2020 at 6:17:59 PM UTC+2, Wolfgang Bangerth wrote:
>>
>> On 4/7/20 2:30 AM, A.Z Ihsan wrote: 
>> > 
>> > Would anyone give me a hint on how to implement the polycrystal 
>> problem, i.e. 
>> > more than one grain boundary? 
>> > I already have the tensor mechanics solver written in deal.ii, but now 
>> it 
>> > needs to be extended solveing also polycrystal problem. 
>>
>> Ihsan, 
>> that's a rather unspecific question. What have you tried already, and what
>>  
>> step specifically is it that you have trouble with? 
>>
>> Best 
>>   WB 
>>
>>
>> -- 
>>  
>> Wolfgang Bangerth  email: bang...@colostate.edu 
>> www: http://www.math.colostate.edu/~bangerth/
>>  
>>
>>
> -- 
> 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 
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to dea...@googlegroups.com .
> To view this discussion on the web visit 
> https://groups.google.com/d/msgid/dealii/aa866c2c-6d5b-4798-a3df-abcb6a811628%40googlegroups.com
>  
> 
> .
>
>
>

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/83805872-e30d-4a4d-9c38-1c91e85b855b%40googlegroups.com.

Re: [deal.II] Polycrystal on deal ii

[Andrew McBride]

In the classical continuum approach to plasticity there is no grain boundary. 
The elastic anisotropy of the crystal is captured at the Gauss point level. So 
you would start by associating different materials properties (elastic / 
plastic) to Gauss points depending on their spatial location. 

Have a look at the CellDataStorage class and how it is used in step-44.

Best
Andrew

> On 7 Apr 2020, at 18:36, A.Z Ihsan  wrote:
> 
> Hi Wolfgang, 
> 
> alright, suppose i have a box that divided into two regions, and each region 
> has different elastic constant (anisotropy material). 
> Perhaps the idea is when calculating the cell_matrix, inside the for-loop 
> there is a function for choosing in which elastic constant to be used 
> depending on the region. 
> 
> I already have the elastic solver written specifically for one elastic 
> constant.  
> Do you have any hint where should i start? 
> 
> BR, 
> Ihsan
> On Tuesday, April 7, 2020 at 6:17:59 PM UTC+2, Wolfgang Bangerth wrote:
> On 4/7/20 2:30 AM, A.Z Ihsan wrote: 
> > 
> > Would anyone give me a hint on how to implement the polycrystal problem, 
> > i.e. 
> > more than one grain boundary? 
> > I already have the tensor mechanics solver written in deal.ii, but now it 
> > needs to be extended solveing also polycrystal problem. 
> 
> Ihsan, 
> that's a rather unspecific question. What have you tried already, and what 
> step specifically is it that you have trouble with? 
> 
> Best 
>   WB 
> 
> 
> -- 
>  
> Wolfgang Bangerth  email: bang...@ <>colostate.edu 
>  
> www: http://www.math.colostate.edu/~bangerth/ 
>  
> 
> 
> -- 
> 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 
> "deal.II User Group" group.
> To unsubscribe from this group and stop receiving emails from it, send an 
> email to dealii+unsubscr...@googlegroups.com 
> .
> To view this discussion on the web visit 
> https://groups.google.com/d/msgid/dealii/aa866c2c-6d5b-4798-a3df-abcb6a811628%40googlegroups.com
>  
> .

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/FF965E85-F69E-4F8C-AE25-6B88AE330BCF%40gmail.com.

Re: [deal.II] Polycrystal on deal ii

[A.Z Ihsan]

Hi Wolfgang, 

alright, suppose i have a box that divided into two regions, and each 
region has different elastic constant (anisotropy material). 
Perhaps the idea is when calculating the cell_matrix, inside the for-loop 
there is a function for choosing in which elastic constant to be used 
depending on the region. 

I already have the elastic solver written specifically for one elastic 
constant.  
Do you have any hint where should i start? 

BR, 
Ihsan
On Tuesday, April 7, 2020 at 6:17:59 PM UTC+2, Wolfgang Bangerth wrote:
>
> On 4/7/20 2:30 AM, A.Z Ihsan wrote: 
> > 
> > Would anyone give me a hint on how to implement the polycrystal problem, 
> i.e. 
> > more than one grain boundary? 
> > I already have the tensor mechanics solver written in deal.ii, but now 
> it 
> > needs to be extended solveing also polycrystal problem. 
>
> Ihsan, 
> that's a rather unspecific question. What have you tried already, and what 
> step specifically is it that you have trouble with? 
>
> Best 
>   WB 
>
>
> -- 
>  
> Wolfgang Bangerth  email: bang...@colostate.edu 
>  
> www: http://www.math.colostate.edu/~bangerth/ 
>
>

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/aa866c2c-6d5b-4798-a3df-abcb6a811628%40googlegroups.com.

Re: [deal.II] Polycrystal on deal ii

[Wolfgang Bangerth]

On 4/7/20 2:30 AM, A.Z Ihsan wrote:


Would anyone give me a hint on how to implement the polycrystal problem, i.e. 
more than one grain boundary?
I already have the tensor mechanics solver written in deal.ii, but now it 
needs to be extended solveing also polycrystal problem.


Ihsan,
that's a rather unspecific question. What have you tried already, and what 
step specifically is it that you have trouble with?


Best
 WB


--

Wolfgang Bangerth  email: bange...@colostate.edu
   www: http://www.math.colostate.edu/~bangerth/

--
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 "deal.II User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/1d131e00-a447-b925-1a2c-cc5e8adb1264%40colostate.edu.

[deal.II] Polycrystal on deal ii

[A.Z Ihsan]

Hi, 

Would anyone give me a hint on how to implement the polycrystal problem, 
i.e. more than one grain boundary?
I already have the tensor mechanics solver written in deal.ii, but now it 
needs to be extended solveing also polycrystal problem. 

THank you. 

Best, 
ihsan

-- 
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 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/eced32b3-fcec-4adf-b24f-3bf868d23743%40googlegroups.com.