AW: AW: [deal.II] step-42 clarification

2018-01-17 Thread Frohne, Joerg
>>On 01/17/2018 02:45 PM, Frohne, Joerg wrote: >> obviously this is a mistake. I checked the source files which we have used >> for results in corresponding paper. >> There I have found the following coding: >> >>resid_old = resid; >> >>resid_vector = system_rhs_newton; >>

Re: AW: [deal.II] step-42 clarification

2018-01-17 Thread Wolfgang Bangerth
On 01/17/2018 02:45 PM, Frohne, Joerg wrote: obviously this is a mistake. I checked the source files which we have used for results in corresponding paper. There I have found the following coding: resid_old = resid; resid_vector = system_rhs_newton; resid_vecto

AW: [deal.II] step-42 clarification

2018-01-17 Thread Frohne, Joerg
Dear Wolfgang, dear community, obviously this is a mistake. I checked the source files which we have used for results in corresponding paper. There I have found the following coding: resid_old = resid; resid_vector = system_rhs_newton; resid_vector.compress(VectorO

Re: [deal.II] Re: the viscous term in SUPG stabilization terms

2018-01-17 Thread Feimi Yu
Thanks a lot, Daniel and Wolfgang! What I understood was that since the SUPG and PSPG terms are integrated over the cell interiors so the laplacians can be ignored while using Q1/Q1 element without any singularity problems on the cell edges. Thanks! Feimi On Wednesday, January 17, 2018 at 2:36

[deal.II] Error with Transfinite Interpolation

2018-01-17 Thread Juan Carlos Araujo Cabarcas
Dear all, Since the inclusion of Transfinite interpolation, I have been successful on working with this powerful technique in my research. I had coded a mesh implementing concentric circles, where the inner most is shifted a small distance s. All concentric circles are labeled 100+i, where i is

Re: [deal.II] Re: the viscous term in SUPG stabilization terms

2018-01-17 Thread Wolfgang Bangerth
On 01/17/2018 11:59 AM, Daniel Arndt wrote: These terms are simply zero if you consider Q1 elements for the velocity. I That's true only if you're on affine meshes. In general, the Laplacian of a Q1 function on an arbitrary mesh is not zero. That doesn't mean, though, that it may not be vali

[deal.II] Re: the viscous term in SUPG stabilization terms

2018-01-17 Thread Daniel Arndt
Feimi, u^h \cdot grad{w^h} \cdot laplacian{u^h} for SUPG and > grad{q^h} \cdot laplacian{u^h} for PSPG. > > Where u^h is the trial function and w^h, q^h refer to the velocity and > pressure test function, respectively. > It seems impossible to calculate these terms with linear elements (because

[deal.II] the viscous term in SUPG stabilization terms

2018-01-17 Thread Feimi Yu
Hi, I'm working on a solver to solve Navier-Stokes equations using the Streamline Upwind Petrov-Galerkin (SUPG) and Pressure Stabilization Petrov-Galerkin (PSPG) method, combined with Newton's iteration and FGMRES iterative method. However, I had some difficulties with implementing the formula