[deal.II] Re: Finite Element Method for Two Phase Flow

I would like to add an idea I intend to apply but right now am clueless of how to apply. In the earlier post I mentioned issue of 'spurious currents'. In the past researchers have tried quite a few approaches to handle this. One of them involves reconstruction of computed pressure field. In ess

Re: [deal.II] Re: run deal.ii on stampede

Zhenlin, The petsc configuration I posted was for complex numbers, which is enabled with --with-scalar-type=complex. If you have not removed this option, then it makes sense that you get the error - the data_out classes do not yet support complex numbers (https://github.com/dealii/dealii/issues

[deal.II] Re: Finite Element Method for Two Phase Flow

Daniel, I had went through step-21 and it appears step-43 extends step-21 with the intentions of improving preconditioning et cetera. In essence both the tutorials solve the same set of equations, namely Darcy's Law, conservation of mass and saturation equation. But for this set of approximate

Re: [deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

> Just post a question/request and hope for the best ;-) I'm pretty sure that > Timo will read this at some point (won't you Timo?). Wait, what? :-) I am skimming most threads of course. Sorry I didn't reply yet. Are any questions left after JP's reply? > Yes, I agree that the fact that there's n

Re: [deal.II] Distribute constraints when using Workstream class

I just understand what you meant ! Thanks ! 2017년 4월 6일 목요일 오후 7시 24분 1초 UTC-5, Wolfgang Bangerth 님의 말: > > On 04/06/2017 02:55 PM, Jaekwang Kim wrote: > > template > > > > void > > > > Step4::copy_local_to_global (constAssemblyCopyData ©_data) > > > > { > > > > for(unsignedinti=0;

[deal.II] Re: Finite Element Method for Two Phase Flow

Sumedh, In my project I am modeling two-phase flow using phase-field model (an > interface capturing method). Now I am done with first draft of my code and > I am facing the well-known issue of 'spurious currents in two-phase flow > using FEM' (paper 1 >

Re: [deal.II] Re: Arbitrary mapping or Manifold?

I will see places where it may fit, I will come back to you. Hej, Toby! nice to hear that it may help, looking forward to see the result =) El viernes, 7 de abril de 2017, 15:42:13 (UTC+2), Wolfgang Bangerth escribió: > > On 04/07/2017 07:24 AM, Juan Carlos Araujo Cabarcas wrote: > > In the

[deal.II] Finite Element Method for Two Phase Flow

Hello, In my project I am modeling two-phase flow using phase-field model (an interface capturing method). Now I am done with first draft of my code and I am facing the well-known issue of 'spurious currents in two-phase flow using FEM' (paper 1

[deal.II] Re: Inhomogeneous Dirichlet bounday conditions with Raviart-Thomas (RT) elements

Sumedh, I intend to use RT elements for my fluid flow problem and the test case > (2D) is shear flow. Consequently the I need to apply constant shear/sliding > velocity (u_x = shear_velocity, u_y = 0) at top and bottom boundaries and > velocities of form u_y = 0, u_x = y*shear_velocity for the

[deal.II] Inhomogeneous Dirichlet bounday conditions with Raviart-Thomas (RT) elements

Hello, I intend to use RT elements for my fluid flow problem and the test case (2D) is shear flow. Consequently the I need to apply constant shear/sliding velocity (u_x = shear_velocity, u_y = 0) at top and bottom boundaries and velocities of form u_y = 0, u_x = y*shear_velocity for the right an

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

Hi Alex, > Yes I began this project by reviewing step-8 > > and > the Handling vector valued problems >

[deal.II] Re: Indexing a Tensor<2,dim> representing the gradient of a Tensor<1,dim>

Yes I began this project by reviewing step-8 and the Handling vector valued problems

Re: [deal.II] Re: Arbitrary mapping or Manifold?

On 04/07/2017 07:24 AM, Juan Carlos Araujo Cabarcas wrote: In the same spirit I would like to share a code that uses ChartManifold, for generating a mapping of an ellipse applied to HyperBall. For both of these cases, do you think that there is a good place somewhere in the documentation where

Re: [deal.II] Re: Arbitrary mapping or Manifold?

Thank you Juan, for sharing your code. I am you have solved your problem. You've also partially solved (and reminded me of) a similar problem I was working on many years (where I failed). I will return to that problem in the near future and I will be implementing your solution to it. :-) Thanks.

Re: [deal.II] Re: Arbitrary mapping or Manifold?

Thank you Juan, for sharing your code. I am you have solved your problem. You've also partially solved (and reminded me of) a similar problem I was working on many years (where I failed). I will return to that problem in the near future and I will be implementing your solution to it. :-) Thanks.

[deal.II] Re: Arbitrary mapping or Manifold?

In the same spirit I would like to share a code that uses ChartManifold, for generating a mapping of an ellipse applied to HyperBall. El miércoles, 5 de abril de 2017, 14:47:17 (UTC+2), Juan Carlos Araujo Cabarcas escribió: > > Dear all, > I recently found the thread: Something wrong with Chart

Re: [deal.II] Re: Crack propagation

Dear Seyed Ali, very good! However, when you write that you solve "exactly" my/our (Timo and I) code then I would assume it is really "exactly" that version. When you add other things, of course this might accelerate or slow down the code :) Next time, please write first if you solve "exactly"

[deal.II] Re: Crack propagation

Dear Thomas, I found the problem. Your version works exactly like discussed in your paper, even faster with my Intel i7, namely around 1400 s. The timeconsuming part is the computation of my B-operator. The B-operator computation is implemented inside* assemble_system (bool residual_only):*

[deal.II] Re: Crack propagation

Hi, > Each Core i7 CPU has maximum 4 physical cores. So, in MPI mode, solving > the simulation using 4 core is faster than 8 cores in your computer. > This is not true. I checked it, it is slower. Best, Seyed Ali -- The deal.II project is located at http://www.dealii.org/ For mailing list/

Re: [deal.II] Re: fe_enriched and step-47

> On 7 Apr 2017, at 03:09, Wolfgang Bangerth wrote: > > On 04/06/2017 06:14 AM, Daniel Arndt wrote: >>I don't know if Wolfgang, the author of the tutorial, plan to finish it. >>Maybe we can team up here. >> >> step-47 seems to be completed in the paper On the implementation of the >> eX