Hello Mohsen,
a general answer is not easy, since it depends on the problem how one
treats the nonlinearity. However in the tutorials you can find two
nonlinear problems, step-15 and step-33.
Best Regards,
Markus
Am 06.07.11 21:12, schrieb Mohsen Ayoobi:
Dear all,
As I see in most of the tutorials when there's a Laplacian for a
variable like U, what is done in assembling the system_matrix is like:
cell_matrix
<http://www.dealii.org/developer/doxygen/deal.II/namespaceLocalIntegrators_1_1Divergence.html#a1a6d2c938fd71efc441cef321ac3baca>(i,j)
+= (fe_values.shape_grad
<http://www.dealii.org/developer/doxygen/deal.II/classFEValuesBase.html#a0a172cb1d5a0e51364b0284288f727a0>
(i, q_point) *
fe_values.shape_grad
<http://www.dealii.org/developer/doxygen/deal.II/classFEValuesBase.html#a0a172cb1d5a0e51364b0284288f727a0>
(j, q_point) *
fe_values.JxW
<http://www.dealii.org/developer/doxygen/deal.II/classFEValuesBase.html#ad097580a2f71878695096cc73b271b9d>
(q_point));
Now I'm wondering how it defers if we have the Laplacian term
multiplied by the same variable we have, ex. U*Laplacian(U).
Thank you in advance.
Mohsen
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
