Ik zou de wet van Darcy numeriek moeten oplossen met Matlab.
Nu weet ik dat er zoiets bestaat als pdetools waarbij ik de PDV oplos a.d.h.v. een eindige elementen methode.
Ik gebruik de volgende vorm van de wet van Darcy:
\(\mathbf{u}=-\kappa\cdot\nabla P\)
In een 2D carthesiaans assenstelsel en vroeg me af welke PDV-vorm deze vergelijking is? (je moet keuze maken in pdetools)Ik dacht eerst dat dit een elliptische PDV was maar als je de vorm bekijkt:
\(-\nabla (c\cdot\nabla u) + a\cdot u = f\)
Maar daar zit je dan met 2 nabla's dat uiteindelijk een laplaciaan geeft.Kan er iemand mij verder helpen hoe ik dit (numeriek) wiskundig probleem oplos?
Bedankt !
(PS: is eerste keer dat ik PDV in Matlab oplos)
Eventjes het kader schetsen: Een beginsnelheid is gekend aan een zijde van een ruimte, ik ken kappa ook en ik zou graag de druk 2D schetsen. Mijn grootste probleem is dat ik bij de overstap naar numerieke wiskunde een beetje verward ben :s
