<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Paul,<div class=""><br class=""></div><div class="">Thank you for this script!! It works very nicely to morph the mesh for the cylinder example.</div><div class="">I modified it to deform the mesh in both x- and y-direction, and made the cylinder into an ellipse (attached picture).</div><div class="">(If someone is interested in the code, I am happy to share it after tidy-up.)</div><div class=""><br class=""></div><div class="">Now I am experimenting with other meshes and blending functions.</div><div class=""><br class=""></div><div class="">Best regards,</div><div class="">Outi</div><div class=""><br class=""></div><div class=""><img height="466" width="423" apple-width="yes" apple-height="yes" apple-inline="yes" id="81C0A7B2-41CC-4237-AC49-5C16ACE182E6" src="cid:6F9E18F7-B6BF-45F9-BD1A-168D9FD2A0E4" class=""></div><div class=""><br class=""></div><div class=""><div><blockquote type="cite" class=""><div class="">On 21 Feb 2015, at 00:34, <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> wrote:</div><br class="Apple-interchange-newline"><div class=""><br class="">Hi Outi,<br class=""><br class="">It's ok... it just means that the solver did not converge all the way... but the<br class="">resulting field is nonetheless sufficiently smooth.<br class=""><br class="">Paul<br class=""><br class="">________________________________________<br class="">From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" class="">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" class="">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>]<br class="">Sent: Friday, February 20, 2015 5:11 PM<br class="">To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="">Subject: Re: [Nek5000-users] Solving mesh Laplace equation (Re: Mesh morphing)<br class=""><br class="">Dear Paul,<br class=""><br class="">Thank you very much, the idea seems exactly like what I would like to<br class="">do. However, there is an error appearing in the Helmholtz solver. From<br class="">logfile:<br class=""><br class=""> call userchk<br class=""> 1.000000000000000E-006 p22 0 1<br class=""> 0 200 **ERROR**: Failed in HMHOLTZ: mshv 1.3206E-02<br class="">2.7147E+03 1.0000E-06<br class=""> done :: userchk<br class=""><br class="">Best regards,<br class="">Outi<br class=""><br class="">Quoting <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>:<br class=""><br class=""><blockquote type="cite" class="">Dear Outi,<br class=""><br class="">I have just added "ocyl2" to the ocyl example directory.<br class=""><br class="">I think it addresses the question you raise. Please advise and let<br class="">us know if you have any difficulties.<br class=""><br class="">Paul<br class=""><br class="">________________________________________<br class="">From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" class="">nek5000-users-bounces@lists.mcs.anl.gov</a><br class="">[<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" class="">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of<br class=""><a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a>]<br class="">Sent: Friday, February 20, 2015 11:39 AM<br class="">To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">nek5000-users@lists.mcs.anl.gov</a><br class="">Subject: [Nek5000-users] Solving mesh Laplace equation (Re: Mesh morphing)<br class=""><br class="">Dear Neks,<br class=""><br class="">What I actually would like to accomplish is to solve the Laplace<br class="">equation for the mesh (with the boundary deformation as a BC). I saw<br class="">that this has been done already at least by Paul (a thread 2009) and<br class="">Matt (for smoothing of a wing mesh, with zero deformation).<br class="">Could any of you please give me a hint on how to solve the Laplace<br class="">equation for the mesh like this? Matt? Paul?<br class="">I would be very grateful.<br class=""><br class="">Best regards,<br class="">Outi<br class=""><br class="">Quoting Outi Tammisola <<a href="mailto:outi@mech.kth.se" class="">outi@mech.kth.se</a>>:<br class=""><br class=""><blockquote type="cite" class="">Dear Neks,<br class=""><br class="">I thought this question might be of interest for someone else as<br class="">well, so posting it here:<br class=""><br class="">I have been looking at the oscillating cylinder example. I would<br class="">like to do something similar, but instead of a prescribed boundary<br class="">velocity, would like to prescribe a given displacement of one<br class="">boundary, and morph my mesh smoothly everywhere to match this<br class="">displacement. For example, if one were to describe a deformation of<br class="">the cylinder rather than a velocity (but without altering the other<br class="">boundaries and without scaling the mesh). This just needs to be done<br class="">once, not at every time step.<br class=""><br class="">Is there an easy way to do this? Can I arrive there somehow by minor<br class="">modifications of the elasticity solver? What equation exactly is the<br class="">elasticity solver solving at every time step?<br class="">(When I outputted the files after every time step for the<br class="">oscillating cylinder, the mesh deformation seemed to match the<br class="">prescribed one only after 2 time steps. Is there a reason for this,<br class="">or did I do something wrong?)<br class=""><br class="">Best regards,<br class="">Outi<br class=""></blockquote><br class=""><br class="">_______________________________________________<br class="">Nek5000-users mailing list<br class=""><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class="">_______________________________________________<br class="">Nek5000-users mailing list<br class="">Nek5000-users@lists.mcs.anl.gov<br class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class=""></blockquote><br class=""><br class="">_______________________________________________<br class="">Nek5000-users mailing list<br class=""><a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class="">_______________________________________________<br class="">Nek5000-users mailing list<br class="">Nek5000-users@lists.mcs.anl.gov<br class="">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class=""></div></blockquote></div><br class=""></div></body></html>