<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jan 18, 2018 at 10:46 PM, 서승진 (기계항공및원자력공학부) <span dir="ltr"><<a href="mailto:mistloin@unist.ac.kr" target="_blank">mistloin@unist.ac.kr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_-3039112721082268532divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<div id="m_-3039112721082268532divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Dear petsc-user.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Thanks for your respond.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"></p>
<div style="margin-top:0px;margin-bottom:0px"><b>2) About application boundary condition</b></div>
<div style="margin-top:0px;margin-bottom:0px"><strong>2-1) About my geometry</strong><br>
</div>
<p></p>
<p style="margin-top:0px;margin-bottom:0px">The reason why I use <span>DMPlexMarkBoundaryFaces()
</span>is to mark the boundary faces where dirichlet condition is applied. </p>
<p style="margin-top:0px;margin-bottom:0px">I accepted your advise, so there is no DMPlexMarkBoundaryFaces() <wbr>in my code.</p>
<p style="margin-top:0px;margin-bottom:0px"> </p>
<p style="margin-top:0px;margin-bottom:0px">I want to explain my problem with my geometry. I want to solve the heat transfer equation (only conduction with internal heat source Q) with four boundary condition.</p>
<p style="margin-top:0px;margin-bottom:0px"> Governing equation is </p>
<p style="margin-top:0px;margin-bottom:0px"><img width="156" height="26" id="m_-3039112721082268532423350_rs|3" style="width:156px;height:26.42px" size="1135" src="cid:8a4a5d6c-1ce7-41da-8085-a846856847f4"></p>
<p style="margin-top:0px;margin-bottom:0px"> Left and Top : Dirichlet condition (T=T1, T=T2) (T1 < T, T2 < T)</p>
<p style="margin-top:0px;margin-bottom:0px"> Right and Bottom : Neumann condition (grad(T)=0 or -k*grad(T)=q) </p>
<p style="margin-top:0px;margin-bottom:0px">Heat is supplied from bottom boundary and generated in internal region, and it is transferred to left and top boundary.
</p>
<p style="margin-top:0px;margin-bottom:0px">Internal heat source is function of temperature, so the problem is non-linear problem.
<br></p></div></div></div></blockquote><div><br></div><div>I recommend using a constant source to start, so you can use MMS and look at the actual error. This is invaluable for debugging, and</div><div>you will see that all my examples do it.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr"><p style="margin-top:0px;margin-bottom:0px">
</p>
<p style="margin-top:0px;margin-bottom:0px"><img width="336" height="434" id="m_-3039112721082268532158976_rs|3" style="width:336px;height:434.24px" size="76848" src="cid:5ba5a6d9-49e2-4e6b-9f27-5453037d2c1e"></p>
<p style="margin-top:0px;margin-bottom:0px"><strong>2-2) About error message</strong></p>
<p style="margin-top:0px;margin-bottom:0px">The error message which I mentioned in previous mail is that</p></div></div></div></blockquote><div><br></div><div>You are giving the name of a label that does not exist I think. You can check what labels are defined (and how many points for each value)</div><div>using -dm_view.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px"><img width="776" height="371" id="m_-3039112721082268532545388_rs|3" style="width:776px;height:371.48px" size="486428" src="cid:de6e76da-0661-4e8a-92f2-aec7512ea42b"></p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><strong>2-3) About inhomogeneous neumann boundary condition
</strong></p>
<p style="margin-top:0px;margin-bottom:0px"></p>
<div> You kindly answer my question; "If they are inhomogeneous, you must specify a part of the weak form, as is done in SNES ex12 and ex77."</div>
<p></p>
<p style="margin-top:0px;margin-bottom:0px"> Is that means I should specify terms "f0_bd_u" and "f1_bd_u" in weak form? (they are integrated along boundary line/surface in weak form)</p></div></div></div></blockquote><div><br></div><div>Yes.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px"><strong>2-4) About dirichlet boundary condition
</strong></p>
<p style="margin-top:0px;margin-bottom:0px"> Simplicial finite element examples (ex12, ex77, ex56) set the exact function ( = exact solution, for example, u=x^2+y^2 in ex12) as dirichlet condition.</p>
<p style="margin-top:0px;margin-bottom:0px"> Is the condition menas that we set the exact "u" values on the boundary (labeled as "marker")?</p></div></div></div></blockquote><div><br></div><div>The Dirichlet conditions is always the exact value on the boundary. That is what it means.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr"><div id="m_-3039112721082268532divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px"> Thanks for reading my email.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Best regard,</p>
<p style="margin-top:0px;margin-bottom:0px">Seungjin Seo<br>
</p>
</div>
<hr style="width:98%;display:inline-block">
<div id="m_-3039112721082268532divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size:11pt"><b>보낸 사람:</b> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
<b>보낸 날짜:</b> 2018년 1월 19일 금요일 오전 1:09:02<br>
<b>받는 사람:</b> 서승진 (기계항공및원자력공학부)<br>
<b>참조:</b> Smith, Barry F.; <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
<b>제목:</b> Re: [petsc-users] I have a question for PETSc example</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div class="m_-3039112721082268532x_gmail_extra">
<div class="m_-3039112721082268532x_gmail_quote">On Thu, Jan 18, 2018 at 6:10 AM, 서승진 (기계항공및원자력공학부) <span dir="ltr">
<<a href="mailto:mistloin@unist.ac.kr" target="_blank">mistloin@unist.ac.kr</a>></span> wrote:<br>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Dear petsc-user.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Hi. I am trying several parts of example code, (SNES ex12.c : Poisson Problem in 2d and 3d with simplicial finite element)</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">and I have several problems to understand.
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I appologize to sent a long and boring email to you...</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><strong>1) About setting of geometry</strong><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><img width="320" height="116" id="m_-3039112721082268532170597_rs|3" style="width:320px;height:115.8px" size="74067" src="cid:009d718e-c50b-4bd9-8d2c-714321b9fc46"></p>
<p style="margin-top:0px;margin-bottom:0px"><img width="470" height="73" id="m_-3039112721082268532156910_rs|3" style="width:470px;height:73.01px" size="43510" src="cid:7e0f732c-ecad-4ea7-a3f3-8d6a09bc5cc9"></p>
<p style="margin-top:0px;margin-bottom:0px"><img width="385" height="24" id="m_-3039112721082268532448976_rs|3" style="width:385px;height:24.41px" size="24165" src="cid:40e3a7a8-4191-4693-bc97-84a550d5152f"></p>
<p style="margin-top:0px;margin-bottom:0px">I referred the "DMPlexCreateBoxMesh" function to make simplicial box mesh.
</p>
<p style="margin-top:0px;margin-bottom:0px">I found that lower[] and upper[] are factors to set the size of the box mesh. Is it right?</p>
<p style="margin-top:0px;margin-bottom:0px">Actually, I want to make a rectangular geometry with size of 3m x 1m, and four boundary (two dirichlet and two neumann condition). </p>
<p style="margin-top:0px;margin-bottom:0px">If I wrong, please give me some tip to make a good geometry.</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Yes, you can use <span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">DMPlexCreateBoxMesh() </span>to <wbr>make
that. </div>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px"></p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><strong>2) About application boundary condition</strong><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><img width="532" height="221" id="m_-303911272108226853295096_rs|3" style="width:532px;height:220.79px" size="109380" src="cid:c6eed790-070c-4e66-8dbb-8b461fb06f10"></p>
<p style="margin-top:0px;margin-bottom:0px"><img width="602" height="148" id="m_-3039112721082268532843541_rs|3" style="width:602px;height:148.29px" size="134256" src="cid:37a6abda-13a5-45bd-b751-2d193f784071"></p>
<p style="margin-top:0px;margin-bottom:0px">I reffered the example "SNES ex69.c" to set a label and condition for the four boundarys.
</p>
<p style="margin-top:0px;margin-bottom:0px">First part of code is about label of boundary condition,</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>This is wrong. You should not have DMPlexMarkBoundaryFaces() in there. What did you want it to do?</div>
<div> </div>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">and second part is about setting of boundary condition.
</p>
<p style="margin-top:0px;margin-bottom:0px">However, due to second part, error occurs.</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>Always, always, always send the full error message. Otherwise, I am just guessing.</div>
<div> </div>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px"></p>
<p style="margin-top:0px;margin-bottom:0px">Is there any wrong argument?</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>If they are homogeneous Neumann conditions, there is no need to mark them. If they are inhomogeneous, you</div>
<div>must specify a part of the weak form, as is done in SNES ex12 and ex77.</div>
<div><br>
</div>
<div>Second, for the Dirichlet conditions, are you sure you only want to constrain the components that I am in ex69?</div>
<div>What would help is to fully specify your problem, so that we know what you want to do.</div>
<div> </div>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Please give me some tip to applied four boundary conditions on the mesh. </p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><strong>3) About input parameters of DSSetJacobian()</strong></p>
<p style="margin-top:0px;margin-bottom:0px">In DSSetJacobian (prob, f, g, g0, g1, g2, g3), I have a problem to understand g0~g3
</p>
<p style="margin-top:0px;margin-bottom:0px"> g0 : integrand for test function and basis function</p>
<p style="margin-top:0px;margin-bottom:0px"> g1 : integrand for test function gradient and basis function</p>
<p style="margin-top:0px;margin-bottom:0px"> g2 : integrand for test function and basis function gradient
</p>
<p style="margin-top:0px;margin-bottom:0px"> g3 : integrand for test function gradient and basis function gradient</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I have studied the FEM before one month, but I am very confused by above four terms.</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>There is a short explanation here:</div>
<div><br>
</div>
<div> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DT/PetscDSSetJacobian.html" target="_blank">http://www.mcs.anl.gov/<wbr>petsc/petsc-current/docs/<wbr>manualpages/DT/<wbr>PetscDSSetJacobian.html</a></div>
<div><br>
</div>
<div>and here you can see the four terms which are being integrated to arrive at the element matrix:</div>
<div><br>
</div>
<div> <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DT/PetscFEIntegrateJacobian.html" target="_blank">http://www.mcs.anl.gov/<wbr>petsc/petsc-current/docs/<wbr>manualpages/DT/<wbr>PetscFEIntegrateJacobian.html</a></div>
<div> </div>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">I should study this in my self, but I need to some tip or internet address which include a helpful description about above terms.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">(I think this issue is not good at this email, but I have been really suffered by above term...) <br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><strong>4) About useful tool to visualize the result and geometry<br>
</strong></p>
<p style="margin-top:0px;margin-bottom:0px"><img id="m_-3039112721082268532x_gmail-m_-9079360143004781812img348174" style="max-width:99.9%" size="8003" src="cid:80900b60-6b68-4c72-8f3e-df0d475795d7"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I also have searched the tool which can express the visualized result and geometry.</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>I use Paraview from KitWare, which is free. I output HDF5, using the generic PETSc Viewer options detailed in the manual,</div>
<div><br>
</div>
<div> -dm_view hdf5:lap.h5 -vec_view hdf5:lap.h5::append</div>
<div><br>
</div>
<div>and then I automatically create an XDMF file lap.xmf, describing the HDF5, which can be read by Paraview,</div>
<div><br>
</div>
<div> $PETSC_DIR/bin/petsc_gen_xdmf.<wbr>py lap.h5</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="m_-3039112721082268532x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Could you give me a tip to select the visualizing tool for better analysis?</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Thanks you for reading my email.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Best regards,</p>
<p style="margin-top:0px;margin-bottom:0px">Seungjin Seo<br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
</div>
<hr style="width:98%;display:inline-block">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size:11pt"><b>보낸 사람:</b> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
<b>보낸 날짜:</b> 2018년 1월 16일 화요일 오후 8:48:16<br>
<b>받는 사람:</b> 서승진 (기계항공및원자력공학부)<br>
<b>참조:</b> Smith, Barry F.; <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">
petsc-users@mcs.anl.gov</a><br>
<b>제목:</b> Re: [petsc-users] I have a question for PETSc example</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail_extra">
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail_quote">On Mon, Jan 15, 2018 at 8:28 PM, 서승진 (기계항공및원자력공학부)
<span dir="ltr"><<a href="mailto:mistloin@unist.ac.kr" target="_blank">mistloin@unist.ac.kr</a>></span> wrote:<br>
<blockquote class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-m_7279604020525382251divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">sorry, I mean that</p>
<p style="margin-top:0px;margin-bottom:0px"><span><strong>My geometry is a rectangle</strong>.</span></p>
<p style="margin-top:0px;margin-bottom:0px">That means, geometry is just rectangle box shape with four boundary.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">My question is that, how can I set the four boundary condition individually in example 12 of SNES?</p>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>You put a different marker on each side. The easiest way to do it is to give the</div>
<div><br>
</div>
<div> -dm_plex_separate_marker</div>
<div><br>
</div>
<div>option. This is done in SNES ex69.</div>
<div><br>
</div>
<div> Thanks,</div>
<div><br>
</div>
<div> Matt</div>
<div> </div>
<blockquote class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail_quote" style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div dir="ltr">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-m_7279604020525382251divtagdefaultwrapper" style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:12pt" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">(I think that the boundary condition should be added in right side of the weak form, where load vector f and test function term is included. However, in example of petsc, boundary condition is applied as general
form.)</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Seungjin Seo<br>
</p>
</div>
<hr style="width:98%;display:inline-block">
<div id="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-m_7279604020525382251divRplyFwdMsg" dir="ltr">
<font color="#000000" face="Calibri, sans-serif" style="font-size:11pt"><span class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-"><b>보낸 사람:</b> Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>><br>
</span><b>보낸 날짜:</b> 2018년 1월 15일 월요일 오후 11:06:18<br>
<b>받는 사람:</b> 서승진 (기계항공및원자력공학부)<br>
<b>참조:</b> Smith, Barry F.; <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">
petsc-users@mcs.anl.gov</a>
<div>
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-h5"><br>
<b>제목:</b> Re: [petsc-users] I have a question for PETSc example</div>
</div>
</font>
<div> </div>
</div>
<div>
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-h5">
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-m_7279604020525382251BodyFragment">
<font size="2"><span style="font-size:11pt">
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail-m_7279604020525382251PlainText">
<br>
The term "rectangular mesh" usually refers to a "structured grid mesh" so it is not clear if you truly need an unstructured grid. Take a look at src/snes/examples/tutorials/ex<wbr>19.c
<br>
<br>
Barry<br>
<br>
<br>
> On Jan 15, 2018, at 1:03 AM, 서승진 (기계항공및원자력공학부) <<a href="mailto:mistloin@unist.ac.kr" target="_blank">mistloin@unist.ac.kr</a>> wrote:<br>
> <br>
> Dear PETSc-user<br>
> <br>
> I am Seungjin Seo, researcher of Korean Advanced Institute of Science and Technology, South Korea.<br>
> <br>
> I am tring to understand the ex12.c in SNES examples, and I have a problem...<br>
> Could you give me some advise? Is there any good description or materials for PetscFE and DMPlex?<br>
> <br>
> I will make a code with unstructured grid with PetscFE and SNES. My geometry is a rectangle mesh (the size is not fixed so I think this is unstructured grid).<br>
> The equation is <br>
> I will set boundry conditions as follow;<br>
> 1. (Top-side boundary) Neumann boundary<br>
> 2. (Bottom-side boundary) Neumann boundary<br>
> 3. (Left-side boundary) Dirichlet boundary<br>
> 4. (Right-side boundary) Dirichlet boundary<br>
> <br>
> In example ex12.c in SNES, load vector f is defiend as f_0, and neumann boundary condition is defined as f_bd_0. f_1 means the test function gradient term. The example set general type of neumann boundary condition on f_bd_0, but I want to set each boundary
condition individually. Could you give me some advise?<br>
> <br>
> Thansk for reading my email.<br>
> <br>
> Best regard,<br>
> Seungjin Seo<br>
> <br>
> <br>
> <br>
> 보낸 사람: Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>><br>
> 보낸 날짜: 2017년 12월 21일 목요일 오후 3:46:42<br>
> 받는 사람: 서승진 (기계항공및원자력공학부)<br>
> 참조: <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br>
> 제목: Re: [petsc-users] I have a question for PETSc example<br>
> <br>
> <br>
> In deed you can, yes many PETSc examples only feature a single field for simplicity, but multiple fields are fine. You simple provide your FormFunction and FormJacobian to handle both degrees of freedom per cell/node at the same time. You use SNES or TS
depending on if your problem is time dependent and then use -pc_type lu to start. Once you have the physics correct come back and ask us about optimizing the linear solver with preconditioners; but until you have the physics correct it is absurd to waste time
worrying about making the linear solver efficient. What you want to do, many people do and is not a big deal.<br>
> <br>
> Barry<br>
> <br>
> <br>
> > On Dec 21, 2017, at 12:27 AM, 서승진 (기계항공및원자력공학부) <<a href="mailto:mistloin@unist.ac.kr" target="_blank">mistloin@unist.ac.kr</a>> wrote:<br>
> > <br>
> > Dear PETSc-user<br>
> > <br>
> > <br>
> > I am Seungjin Seo, researcher of Korean Advanced Institute of Science and Technology, South Korea.<br>
> > <br>
> > <br>
> > <br>
> > I am trying to solve thermal and fluid equations within a porous structure.<br>
> > The thermal equation includes a non-linear term of pressure and the fluid equation has a boundary condition using temperature gradient.<br>
> > Is there any way I can solve these two equations at the same time, instead of solving temperature first with previous pressure distribution and then pressure and update temperature?<br>
> > Many examples solve only one PDF instead of coupling several physics.<br>
> > Please recommend me which solver is appropriate for my case?<br>
> > <br>
> > Best regards,<br>
> > Seungjin Seo<br>
<br>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></div>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></span></font></div>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></div>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></div>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></div>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></blockquote>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"></font></span></div>
<span class="m_-3039112721082268532x_gmail-HOEnZb"><font color="#888888"><br>
<br clear="all"><span class="HOEnZb"><font color="#888888">
<div><br>
</div>
-- <br>
<div class="m_-3039112721082268532x_gmail-m_-9079360143004781812x_gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener</div>
<div><br>
</div>
<div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~k<wbr>nepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</font></span></font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="m_-3039112721082268532x_gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener</div>
<div><br>
</div>
<div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~<wbr>knepley/</a><br>
</div>
</div>
</div>
</div>
</div>
</font></span></div>
</div>
</div>
</div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>