Re: [deal.II] Re: How do I choose my spaces to enforce gradient continuity across an interface?

Dear Prof. Bangerth, Thank you very much for your reply! Apologies for the ambiguity. What I meant is for the interface condition to hold weakly on the boundary and that, as you correctly pointed out, the jump in the gradient should equal the pressure times the normal vector. As such, I spec

Re: [deal.II] Re: How do I choose my spaces to enforce gradient continuity across an interface?

and side to both sides (i.e. with zero pressure on both sides). This does not happen if I solve just stokes everywhere or just laplace (with the same right hand side :i.e. zero contribution from pressure). Kind regards, Georgios Sialounas On Wednesday, August 22, 2018 at 5:42:34 AM UTC+1, Wolfg

[deal.II] Re: How do I choose my spaces to enforce gradient continuity across an interface?

Dear Daniel, Thanks for your reply! Essentially what I was trying to do is patch together two product spaces: ([H1_0]^2, L^2_0) in one part of the domain and ([H1_0]^2, L^2_0) in another, with an interface condition to ensure solvability for my problem (Stokes in one part of the domain, laplac

Re: [deal.II] Step 46: pressure uniquely defined

Dear Time, Thank you very much! I've been using the preconditioned CG solver that came with the tutorial thus far. I'll have a crack with GMRES. Kind regards, Georgios On Sunday, August 19, 2018 at 1:23:29 AM UTC+1, Timo Heister wrote: > > Georgios, > > this is discussed very briefly in step-

Re: [deal.II] How do I calculate the (hdiv-semi) norm only on certain components of my (vector-valued) solution?

Dear Prof. Bangerth, Thank you very much! Kind regards, Georgios On Sun, Aug 12, 2018 at 10:05 PM, Wolfgang Bangerth wrote: > On 08/11/2018 09:37 AM, georgios.sialou...@gmail.com wrote: > >> >> I would be most grateful for your help on the following matter. I am >> solving a combined stokes-l

[deal.II] How do I choose my spaces to enforce gradient continuity across an interface?

_schur.cpp // stokes_cg // // Created by Georgios Sialounas on 19/08/2018. // Copyright © 2018 Georgios Sialounas. All rights reserved. // /* - * * Copyright (C) 2011 - 2015 by the deal.II authors * * This file is part of the deal.

[deal.II] Step 46: pressure uniquely defined

Dear all, In step 46 (and step 22 too if possible) for stokes, how to we ensure that the pressure is uniquely defined (i.e. that it is in L^2_0)? Kind regards, Georgios -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/

[deal.II] How do I calculate the (hdiv-semi) norm only on certain components of my (vector-valued) solution?

Dear all, I would be most grateful for your help on the following matter. I am solving a combined stokes-laplace problem (have a vector valued solution). I am trying to calculate the hdiv norm on the first two components of my solution (i.e. the two components of the velocity), cellwise.