From: Povolotskyi, M. <mpo...@pu...> - 2019-03-05 19:54:43
|
Thank you, Roy. One more question, please. How can I get dp/dt, if I need? Is there an existing solution for this? Thank you! Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows 10 From: John Peterson<mailto:jwp...@gm...> Sent: Monday, March 4, 2019 12:38 PM To: Povolotskyi, Mykhailo<mailto:mpo...@pu...> Cc: libmesh-users<mailto:lib...@li...> Subject: Re: [Libmesh-users] question about example transient_ex2.C On Mon, Mar 4, 2019 at 11:34 AM Michael Povolotskyi <mpo...@pu...<mailto:mpo...@pu...>> wrote: Thank you. Couple of more questions: 1. You have p and u variables in your equation. Is this a misprint? Yes, sorry there is only one variable, p. 2. It seems there is also a force term. In the example you integrate it over a surface. Can I have a volume force, to solve something like this: (1/c^2) d^2 p / dt^2 - Laplacian(p) + F = 0 ? Sounds reasonable to me. 3. Where cam I read about the Newmark method of solution, that you have implemented? I would assume it something similar to this: https://en.wikipedia.org/wiki/Newmark-beta_method but of course looking at the code to see what it does is your best bet... -- John |