Dispersie
Moderator: ArcherBarry
Forumregels
(Middelbare) school-achtige vragen naar het forum "Huiswerk en Practica" a.u.b.
Zie eerst de Huiswerkbijsluiter
(Middelbare) school-achtige vragen naar het forum "Huiswerk en Practica" a.u.b.
Zie eerst de Huiswerkbijsluiter
-
- Berichten: 11
Dispersie
Dit is een wiskundig probleem dat opdook bij een chemisch-fysisch experiment...
Voor een project heb ik dispersiemetingen gedaan. Ik heb concentraties enz. bepaald en een geschikte formule gevonden om dit te verwerken tot een dispersiecoëfficiënte:
C = ( M / sqrt(4 * pi * D * t) ) * exp( -(x - V * t)^2 / (4 * D * t)
MET: C = concentratie (azijnzuur) in mol / liter
M = hoeveelheid (azijnzuur) in mol
D = dispersiecoëfficiënt
x = locatie vanaf inspuitpunt in meter
V = stroomsnelheid in meter / seconde
t = tijd in seconden
Alle gegevens zijn bekend, behalve de dispersiecoëfficiënt. Deze wil ik hiermee berekenen. Daartoe heb ik de formule omgeschreven naar D = ... m.b.v. Maple:
D := -1/2 * (-x + V * t)^2 / (LambertW(-2 * C^2 * pi * (-x + V * t)^2 / (M^2)) * t)
Daarbij is LambertW (x) gedefinieerd als:
LambertW(x) * exp(LambertW(x)) = x
In Maple krijg ik echter geen antwoord uit deze vergelijking nadat ik alle waarden heb ingevuld. Maple rekent het LambertW deel niet uit. Hoe kan ik dit probleemstuk oplossen? Hoeft niet per se met Maple natuurlijk.
Bij voorbaat dank, Ron van Dorp
Voor een project heb ik dispersiemetingen gedaan. Ik heb concentraties enz. bepaald en een geschikte formule gevonden om dit te verwerken tot een dispersiecoëfficiënte:
C = ( M / sqrt(4 * pi * D * t) ) * exp( -(x - V * t)^2 / (4 * D * t)
MET: C = concentratie (azijnzuur) in mol / liter
M = hoeveelheid (azijnzuur) in mol
D = dispersiecoëfficiënt
x = locatie vanaf inspuitpunt in meter
V = stroomsnelheid in meter / seconde
t = tijd in seconden
Alle gegevens zijn bekend, behalve de dispersiecoëfficiënt. Deze wil ik hiermee berekenen. Daartoe heb ik de formule omgeschreven naar D = ... m.b.v. Maple:
D := -1/2 * (-x + V * t)^2 / (LambertW(-2 * C^2 * pi * (-x + V * t)^2 / (M^2)) * t)
Daarbij is LambertW (x) gedefinieerd als:
LambertW(x) * exp(LambertW(x)) = x
In Maple krijg ik echter geen antwoord uit deze vergelijking nadat ik alle waarden heb ingevuld. Maple rekent het LambertW deel niet uit. Hoe kan ik dit probleemstuk oplossen? Hoeft niet per se met Maple natuurlijk.
Bij voorbaat dank, Ron van Dorp
-
- Berichten: 328
Re: Dispersie
Waarom doe je het niet gewoon met een rekenmachine?
C = ( M / sqrt(4 * pi * D * t) ) * exp( -(x - V * t)^2 / (4 * D * t)
( M / sqrt(4 * pi * D * t) ) = C / (exp( -(x - V * t)^2 / (4 * D * t))
sqrt(4 * pi * D * t) = M / C / (exp( -(x - V * t)^2 / (4 * D * t))
(4 * pi * D * t) = (M / C / (exp( -(x - V * t)^2 / (4 * D * t)))2
D = ((M / C / (exp( -(x - V * t)^2 / (4 * D * t)))2)/4 * pi * t (denk ik )
Maar met maple moet je het ook iets anders invullen, zonder spaties bijv. En ik snap niet wat je bedoelt met:
[quote]Daarbij is LambertW (x) gedefinieerd als:
LambertW(x) * exp(LambertW(x)) = x[/quote]
C = ( M / sqrt(4 * pi * D * t) ) * exp( -(x - V * t)^2 / (4 * D * t)
( M / sqrt(4 * pi * D * t) ) = C / (exp( -(x - V * t)^2 / (4 * D * t))
sqrt(4 * pi * D * t) = M / C / (exp( -(x - V * t)^2 / (4 * D * t))
(4 * pi * D * t) = (M / C / (exp( -(x - V * t)^2 / (4 * D * t)))2
D = ((M / C / (exp( -(x - V * t)^2 / (4 * D * t)))2)/4 * pi * t (denk ik )
Maar met maple moet je het ook iets anders invullen, zonder spaties bijv. En ik snap niet wat je bedoelt met:
[quote]Daarbij is LambertW (x) gedefinieerd als:
LambertW(x) * exp(LambertW(x)) = x[/quote]
-
- Berichten: 2.399
Re: Dispersie
[/quote]Roytje schreef: D = ((M / C / (exp( -(x - V * t)^2 / (4 * D * t)))2)/4 * pi * t (denk ik )
Maar met maple moet je het ook iets anders invullen, zonder spaties bijv. En ik snap niet wat je bedoelt met:
[quote]Daarbij is LambertW (x) gedefinieerd als:
LambertW(x) * exp(LambertW(x)) = x
Nee want dan heb je D nog niet geisoleerd. Maar ik snap ook niet wat er met lambertW bedoeld wordt. Verder heb je toch niet ergens een 0 ingevuld he?
-
- Berichten: 311
Re: Dispersie
LambertW(x) * eLambertW(x)=x heeft volgens mij alleen een oplossing als LambertW(x) en x nul zijn of LambertW(x) is -∞ en x is nul.
-
- Berichten: 11
Re: Dispersie
Hey, als ik er ook maar iets van begreep, dan had ik het zelf wel opgelost he! Het probleem is helaas niet zo eenvoudig als voorheen werd geschetst, omdat de D niet gemakkelijk te isoleren is.
- Berichten: 3.507
Re: Dispersie
In ieder geval wel fijn dat je snel reageerde, Ron
I'm not suffering from insanity - I'm enjoying every minute of it!!
-
- Berichten: 11
Re: Dispersie
Eeeh... tjaIn ieder geval wel fijn dat je snel reageerde, Ron
Ik ben een tijdje weggeweest, enne... deze formule moest de werkelijk gemeten waardes bij een experiment simuleren. Uiteindelijk heb ik, omdat alle overige variabelen in de formule wel te bepalen waren, trial en error toegepast. Ik heb net zolang de D veranderd totdat de meetwaardegrafiek en die van de theorie zoveel mogelijk overeenkwamen.
-
- Berichten: 2.035
Re: Dispersie
Ik geloof dat Ron gewoon probeert het record "oude topics omhoog schoppen" te verbeteren!In ieder geval wel fijn dat je snel reageerde, Ron
-
- Berichten: 11
Re: Dispersie
"oude topics omhoog schoppen" staat beter bekend als necro-posting
-
- Berichten: 2.035
Re: Dispersie
Nee, ook al verander je de naam van het record dat je probeert te verbeteren, je krijgt maar 1 prijs!"oude topics omhoog schoppen" staat beter bekend als necro-posting
- Berichten: 6.314
Re: Dispersie
... en weer ontopic.
You can't possibly be a scientist if you mind people thinking that you're a fool. (Douglas Adams)
-
- Berichten: 2.035
Re: Dispersie
Het "ontopic" deel van de discussie was op 24 mei 2005 al sort of afgesloten