[deal.II] How to solve a 2D problem on the boundary in 3D simulations

Hi, I am using DGP to solve a 3D problem which requires me to first solve a 2D problem (boundary condition) on the boundary at each time step. My question is how to assemble a 2D problem on the boundary in 3D simulations. I cannot find similar tutorials to start. Do I just use shape functions

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

0021-9991", doi = "https://doi.org/10.1016/j.jcp.2018.11.029";, url = "http://www.sciencedirect.com/science/article/pii/S0021999118307733";, author = "Jiaqi Zhang and Pengtao Yue", keywords = "Hamilton–Jacobi equation, Numerical flux, Second-derivative lim

Re: [deal.II] A bug of visualization with VISIT

Hey Prof. Bangerth, Thanks for your help again, I understand now :) Best, Jiaqi 在 2017年2月22日星期三 UTC-5下午11:58:06，Wolfgang Bangerth写道： > > On 02/22/2017 06:41 AM, Jiaqi ZHANG wrote: > > Soluiton2 is fine. For Solution1, I checked the result with MATLAB, the > > contour plo

Re: [deal.II] A bug of visualization with VISIT

e triangles are aligned with the solution, > while in case 1 there are not (the domain is split in two triangle with a > common edge going from top left to bottom right). > > Best, > > Bruno > > On Wednesday, February 22, 2017 at 8:41:54 AM UTC-5, Jiaqi ZHANG wrote: >&g

Re: [deal.II] obtaining solution values at any point using Quadrature< dim >::Quadrature(const Point< dim > & point)

Got it. Thanks so much, Dr. Bangerth. Jiaqi 2016-12-16 0:06 GMT-05:00 Wolfgang Bangerth : > On 12/15/2016 06:11 PM, Jiaqi ZHANG wrote: > >> std::cout<<"testing quadrature points:\n"; >> const Point v0 = cell->vertex(0); >> Quadrature q_point(v0); &

[deal.II] obtaining solution values at any point using Quadrature< dim >::Quadrature(const Point< dim > & point)

Dear all, I was trying to write a roots finding method along cell faces in 2D, so I need to know the solution values at any point of cell faces. Is this the right way to use Quadrature < dim >:: Quadrature

Re: [deal.II] Interpolating values of the neighbor cells in DG

Thank you so much, Praveen. Your codes are really helpful. Regards, Jiaqi 2016-09-20 0:11 GMT-04:00 Praveen C : > > On Tue, Sep 20, 2016 at 9:34 AM, Jiaqi ZHANG > wrote: > >> Hello Praveen, >> >> Thank you so much. I am wondering if you could tell me which >&

Re: [deal.II] Interpolating values of the neighbor cells in DG

s, 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+unsubs

Re: [deal.II] Interpolating values of the neighbor cells in DG

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. > For more options, visit https://groups.google.com/d/optout. > -- Jiaqi Zhang Graduate

Re: [deal.II] Interpolating values of the neighbor cells in DG

to extend the solutions from neighbors through interpolation? I am not sure how to do that in deal.II. Thanks, Jiaqi 2016-09-19 17:36 GMT-04:00 Wolfgang Bangerth : > On 09/19/2016 03:15 PM, Jiaqi ZHANG wrote: > >> >> I was trying to get the values of the solution polynomial

[deal.II] Interpolating values of the neighbor cells in DG

Dear all, I was trying to get the values of the solution polynomial in the neighbor cells in DG, for example, I want to compute [image: View the MathML source] where p_i are the solution polynomials of the neighbors, \Delta_0 is the current cell. I cannot find any function to do this in Fini

Re: [deal.II] Re: How to find all the partial derivatives of the solution polynomial in deal.II

rum/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. > For more options, v

[deal.II] How to find all the partial derivatives of the solution polynomial in deal.II

Dear all, Is there any function in deal.II that can compute all the partial derivatives of the solution polynomial? For example, I want to compute the following term: [image: View the MathML source] where \tilda{p(x,y)} is the DG solution in one cell in terms of basis functions. In order to

Re: [deal.II] Re: How to apply a point source condition (point boundary condition) in Discontinuous Galerkin methods

orTools.html# > ad03b858b1a3b59003a76f6224e67efc7) > > Best, > > Bruno > > > On Monday, August 29, 2016 at 10:39:27 AM UTC-4, Jiaqi Zhang wrote: >> >> Dear All, >> >> I was trying to use DG to solve a Hamilton-Jacobi equation, and the >> problem has a poi

[deal.II] How to apply a point source condition (point boundary condition) in Discontinuous Galerkin methods

Dear All, I was trying to use DG to solve a Hamilton-Jacobi equation, and the problem has a point source condition for example, the center of a unit square pointCenter = 0. First I thought of using ZeroSourcePoint.add_line(cell->vertex_dof_index(v,0)); but later found that DG elements aren't d

Re: [deal.II] Re: Step-35 Poiseuille flow in a 3D pipe

and their > implications (go read the relevant papers) but rather more technical > aspects regarding the implementation of finite element discretizations in > deal.II. See also here: > https://groups.google.com/forum/#!topic/dealii/GRZMUTLIm2I > > Best, > Martin > > On 03.08.

Re: [deal.II] Re: Step-35 Poiseuille flow in a 3D pipe

> On Friday, July 29, 2016 at 9:19:29 PM UTC-4, Jiaqi ZHANG wrote: >> >> Hello David Wells, >> >> Sorry to bother you, but I have a problem about Step-35. I have posted >> the question >> >> a few days ago, and no one answered me, so I have to sea

Re: [deal.II] Re: A question about velocity correction in Step-35

> email to dealii+unsubscr...@googlegroups.com. > For more options, visit https://groups.google.com/d/optout. > -- Jiaqi Zhang Graduate student Department of Mathematics Virginia Tech -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see h

[deal.II] Re: Step-35 Poiseuille flow in a 3D pipe

Hello David Wells, Sorry to bother you, but I have a problem about Step-35. I have posted the question a few days ago, and no one answered me, so I have to search in the mailing list and found that you have been using it a lot. I was wondering if you could help me, and the following is my q

[deal.II] A question about velocity correction in Step-35

Dear all, I cannot find anything about velocity correction in the code of Step-35. I read the references listed there and some other papers about projection methods, and know that after obtaining the intermediate values for velocity, say v_n, we need to correct the velocity by u_n = v_n + g

Re: [deal.II] Re: How to add values of solutions from different FE spaces (separate DoFHandlers)

t-the-mapped-position-of-support-points-of-my-element > > > Am Montag, 25. Juli 2016 19:33:48 UTC+2 schrieb Jiaqi ZHANG: >> >> Dear all, >> >> I am trying to solve the NS equation using project methods, and follow >> the idea of Step-35 to treat veloci

[deal.II] How to add values of solutions from different FE spaces (separate DoFHandlers)

Dear all, I am trying to solve the NS equation using project methods, and follow the idea of Step-35 to treat velocity and pressure separately, but when I am not sure how to implement the correction of the velocity "u_n_plus_1 = u_star - dt * dp", where "u_star" and "dp" are from the diffusio

[deal.II] Re: Output problem of DG with degree =0

tml#ace4b76e565ba0701c4d32c26075ed3b9> > (solution, "u", > DataOut::type_dof_data > <https://dealii.org/8.4.0/doxygen/deal.II/classDataOut.html>); > > I hope that this solves your problem! > J-P > > On Saturday, July 16, 2016 at 6:04:39 PM UTC+2, Jiaqi Zhang w

[deal.II] Output problem of DG with degree =0

Hey all, I am trying to solve a level set equation with DG. When I use degree=1, everything is ok, but when I change it to degree=0, it cannot output the results, and says: An error occurred in line <957> of file in function void dealii::DataOut_DoFData::add_data_vector(const VectorType&