On Thu, 21 Feb 2013, Manav Bhatia wrote:
Upon further study, it seems like the elem_fixed_solution stores the
element solution in situations where the elem_solution may contain
the velocity or delta solution. Is this correct?
Currently, yes.
---
Roy---
Thanks a lot for the response, Roy.
I did study the code further, and was able to make some progress. I am
currently debugging some AMR issues in my code.
I think some clarity in the API would be helpful. I learnt the hard way
since I assumed that my flux boundary conditions were supposed to go in
On Thu, 28 Feb 2013, Manav Bhatia wrote:
I think some clarity in the API would be helpful. I learnt the hard way since I
assumed that my
flux boundary conditions were supposed to go into side_constraint, only to
realize that it was
not going to get multiplied by dt unless it was in side_resid
Sorry about the delays, but I'm finally finding some time to work
through old emails now.
On Thu, 21 Feb 2013, Manav Bhatia wrote:
I am now sure if I understand the concept of "elem_fixed_solution".
That's my fault. The FEMSystem time integration APIs were
"prematurely optimized", readabili
On Thu, 21 Feb 2013, Manav Bhatia wrote:
What is throwing me off is that the class documentation of Euler2_Solver says:
* Euler solves u' = f(theta*u_new + (1-theta)*u_old),
* Euler2 solves u' = theta*f(u_new) + (1-theta)*f(u_old)
which seems to imply that no mass matrix is accounted for.
On Feb 21, 2013, at 3:25 PM, Roy Stogner wrote:
>>
>> Additionally, within the existing framework, how does one account for
>> *nonlinear* systems that have a mass matrix
>>
>> [M] du/dt = f(u)
>
> Right now, we allow you to request an "elem_fixed_solution" that can
> be used for ei
On Feb 21, 2013, at 3:25 PM, Roy Stogner wrote:
>
>> Additionally, within the existing framework, how does one account for
>> *nonlinear* systems that have a mass matrix
>>
>> [M] du/dt = f(u)
>
> Right now, we allow you to request an "elem_fixed_solution" that can
> be used for either low
Hi Roy,
In the FEM code that I have written, I have implemented Newmark Method
for first and second order nonlinear systems, with and without mass and
damping matrices. While converting second-order to first-order
systems I create a new sparsity patter data structure where the
sub-matrices inhe
On Thu, 21 Feb 2013, Manav Bhatia wrote:
> I would appreciate some help in understanding if the existing System
> classes can handle PDEs (linear/nonlinear) with second order time
> derivative.
IIRC We have one existing system (NewmarkSystem?) designed for second
order time derivatives. I hav
Hi,
I would appreciate some help in understanding if the existing System
classes can handle PDEs (linear/nonlinear) with second order time
derivative. I am going through the FEMSystem documentation and it seems
like these are suitable for a system of the form
du/dt = f(u).
Ofcourse, I
10 matches
Mail list logo
