Dimensional consistency in fractional differential equations with non singular kernels

Dit artikel presenteert een methode met een eenvoudige variabeletransformatie om de dimensionale consistentie te waarborgen in fractionele differentiaalvergelijkingen met niet-singuliere kernen, en illustreert dit aan de hand van een voorbeeld.

De kern

De "Maatwerk" voor Wiskundige Formules: Een Simpele Uitleg van Gonzálezes Papier

Stel je voor dat je een recept hebt voor een overheerlijke taart (een klassiek natuurkundig model). Dit recept werkt perfect zolang je de ingrediënten in gram en minuten meet. Maar wat als je ineens besluit om de taart niet meer in minuten, maar in "fractie-minuten" te bakken? Of nog gekker: wat als je de tijd niet meer als een rechte lijn ziet, maar als een sliert die soms vastloopt en soms versnelt?

Dat is precies wat er gebeurt in de wereld van fractale differentiaalvergelijkingen. Wetenschappers gebruiken deze wiskunde om complexe systemen te beschrijven, zoals hoe energie verloren gaat in een batterij of hoe een veer trilt. Maar er zit een groot probleem in: als je de simpele "tijd" in je formule vervangt door deze complexe "fractale tijd", breekt de eenheid. Het is alsof je in je taatrecept ineens "gram" vervangt door "een beetje", en dan vraagt de bakker: "Hoeveel gram is dat nou precies?" De formule klopt niet meer; de eenheden zijn in de war.

Het Probleem: De "Gekke" Kernen
In de wiskunde zijn er verschillende manieren om deze "fractale tijd" te definiëren.

1. De oude manier (Riemann-Liouville/Caputo): Dit werkt met een formule die op het begin (tijd = 0) oneindig hoog wordt. Het is alsof je probeert een foto te maken van een explosie op het exacte moment dat het ontploft; het beeld is wazig en onbruikbaar.
2. De nieuwe manier (Caputo-Fabrizio): Dit is de methode waar dit papier over gaat. Hier gebruiken ze een "gladde" formule (een exponentiële kromme) die niet explodeert aan het begin. Het is alsof je een soepel filmpje maakt in plaats van een wazige foto. Dit is veel makkelijker om te gebruiken voor echte natuurkunde, zoals in elektrische schakelingen.

Het Dilemma: De Maatstaf is Verdwijnen
Het probleem met deze nieuwe, gladde methode is dat de wiskundige operator (de tool om de verandering te meten) geen eenheid heeft. Het is alsof je een liniaal hebt die geen centimeters of inches aangeeft, maar gewoon "lengte" meet. Als je deze liniaal gebruikt in een formule voor een elektrische schakeling (zoals een weerstand en een condensator), krijg je een rommeltje: je kunt geen "volt" vergelijken met een "wiskundige grootheid zonder naam".

De Oplossing: Een Nieuwe "Tijds-Bril"
Gabriel González, de auteur van dit papier, komt met een slimme oplossing. Hij zegt: "We moeten de tijd niet zomaar vervangen, we moeten de tijd herschrijven."

Hij introduceert een nieuwe functie, laten we hem $\phi$ (fi) noemen.

• De Analogie: Stel je voor dat je door een bril kijkt die de tijd vertraagt of versnelt, afhankelijk van hoe "fractaal" het systeem is.
• In de gewone wereld (waar $\alpha = 1$) is deze bril leeg; je ziet de tijd zoals hij is.
• In de fractale wereld (waar $0 < \alpha < 1$) is de bril gekleurd. Hij zorgt ervoor dat de eenheden (zoals seconden) weer kloppen, zelfs als je de complexe fractale formule gebruikt.

González zegt eigenlijk: "Vervang niet zomaar 'dt' (een klein stukje tijd) door 'CFD' (de fractale tool). Vervang 'dt' door 'dt gedeeld door een slimme factor $\phi$'."

Door deze factor $\phi$ toe te voegen, wordt de formule weer "dimensioneel consistent". De eenheden vallen weer netjes op hun plek, alsof je de liniaal weer hebt voorzien van een schaalverdeling.

Het Voorbeeld: De RC-Schakeling
Om dit te bewijzen, pakt González een heel bekend voorbeeld: een RC-schakeling (een weerstand en een condensator die samen een batterij opladen).

• Normaal: Als je de batterij aansluit, stijgt de spanning in de condensator soepel tot hij vol is.
• Met de nieuwe methode: González gebruikt zijn nieuwe "tijd-bril" ($\phi$) om de formule aan te passen. Hij lost de vergelijking op en krijgt een nieuw antwoord.
• Het Resultaat: Als je de "fractale knop" (de parameter $\alpha$) op 1 zet, krijg je exact hetzelfde antwoord als de klassieke natuurkunde. Maar als je $\alpha$ lager zet (bijvoorbeeld 0.8), zie je iets interessants: het systeem dissipeert (verliest) energie op een andere manier. Het gedraagt zich alsof er interne wrijving is die niet lineair is.

Waarom is dit belangrijk?
Vroeger moesten wetenschappers "magische getallen" (zoals $\sigma$) toevoegen aan hun formules om de eenheden te laten kloppen, maar niemand wist echt wat die getallen fysiek betekenden.
González' methode is slimmer: hij gebruikt een functie die afhangt van de tijd zelf. Het is alsof je niet een statische maatstaf gebruikt, maar een dynamische liniaal die meebeweegt met het systeem.

Conclusie in het Kort
Dit papier zegt: "Hé, als je die nieuwe, mooie, gladde wiskunde (Caputo-Fabrizio) gebruikt om echte dingen te modelleren, zorg dan dat je de eenheden niet verliest. Doe dit door een slimme 'tijd-functie' toe te voegen die de eenheden in balans houdt."

Het is als het maken van een vertaling van een boek naar een vreemde taal: je kunt niet zomaar woorden vervangen; je moet de grammatica aanpassen zodat de zin nog steeds logisch en correct klinkt. González heeft de grammaticaregels gevonden om de "fractale taal" van de natuurkunde correct te laten spreken.

Technische samenvatting

Samenvatting: Dimensionele consistentie in fractionele differentiaalvergelijkingen met niet-singuliere kernen

Auteur: Gabriel González (Universidad de Guadalajara, Mexico)
Onderwerp: Fractionele Calculus, Caputo-Fabrizio afgeleide, RC-schakeling, Dimensionele consistentie.

---

1. Het Probleem

Fractionele calculus wordt steeds vaker gebruikt om complexe dissipatieve fenomenen te modelleren die geheugen-effecten vertonen. Traditionele definities, zoals de Riemann-Liouville en Caputo-afgeleiden, gebruiken kernen met een machtsfunctie ($t-s)^{-\alpha}$ die singulier zijn bij de oorsprong. Dit kan leiden tot wiskundige moeilijkheden en fysische interpretatieproblemen bij beginvoorwaarden.

Om dit op te lossen, hebben Caputo en Fabrizio een nieuwe definitie voorgesteld met een niet-singuliere exponentiële kern. Een fundamenteel probleem bij het toepassen van deze (en andere) fractionele afgeleiden in fysische modellen is echter de dimensionele inconsistentie:

• De klassieke tijdsafgeleide $d/dt$ heeft de dimensie van tijd$^{-1}$ (bijv. $s^{-1}$).
• De standaard fractionele afgeleide $D^\alpha$ heeft de dimensie tijd$^{-\alpha}$.
• Als men in een fysische vergelijking simpelweg $d/dt$ vervangt door $D^\alpha$, ontstaat er een dimensionele onbalans, wat de fysische interpretatie van het model onmogelijk maakt.

Bestaande oplossingen, zoals het invoeren van een schaalparameter $\sigma$ (met dimensie tijd) om de eenheden te herstellen, worden vaak als willekeurig beschouwd of leiden tot parameters die moeilijk fysisch te interpreteren zijn (bijv. eenheden als $s^{-\alpha}$).

2. Methodologie

De auteur stelt een systematische methode voor om dimensionele consistentie te behouden bij het gebruik van de Caputo-Fabrizio (CF) fractionele afgeleide, zonder willekeurige hulpparameters in te voeren. De aanpak bestaat uit de volgende stappen:

1. Analyse van de CF-afgeleide: De CF-afgeleide wordt gedefinieerd als een dimensieloze operator. Om deze te koppelen aan een fysische grootheid met tijdseenheden, moet een nieuwe variabele worden geïntroduceerd.
2. Invoering van een tijdsfunctie $\tau(t, \alpha)$: In plaats van een constante parameter, wordt een functie $\tau$ geïntroduceerd die afhangt van de tijd $t$ en de fractionele orde $\alpha$. Deze functie wordt gedefinieerd via een integraal:
   $$\tau(t, \alpha) = \int_0^t \frac{dt}{\phi(t, \alpha)}$$
   waarbij $\phi(t, \alpha)$ een functie is met de dimensie van tijd.
3. Vervanging van de afgeleide: De gewone tijdsafgeleide wordt vervangen door:
   $$\frac{d}{dt} \rightarrow \frac{1}{\phi(t(\tau), \alpha)} \, {}^{CF}D^\alpha_\tau$$
4. Randvoorwaarde: Er wordt een voorwaarde opgelegd zodat de methode terugkeert naar de klassieke differentiatie wanneer $\alpha \to 1$:
   $$\frac{1}{\phi(t, 1)} \frac{d}{d\tau} = \frac{d}{dt}$$
5. Toepassing op lineaire systemen: De methode wordt toegepast op een eerste-orde differentiaalvergelijking met constante coëfficiënten, die wordt omgezet in een fractionele vergelijking in de variabele $\tau$. De oplossing wordt gevonden via een integrerende factor.

3. Belangrijkste Bijdragen

• Nieuwe aanpak voor dimensionele consistentie: De paper introduceert een methode die geen willekeurige schaalparameters ($\sigma$) vereist, maar in plaats daarvan een functie $\phi(t, \alpha)$ gebruikt die de dimensie van tijd behoudt en afhankelijk is van de fractionele orde.
• Specifiek voor niet-singuliere kernen: De methode is specifiek ontworpen voor de Caputo-Fabrizio operator, waarbij de dimensieloze aard van de operator wordt gecompenseerd door de tijdsfunctie.
• Analytische oplossing: De auteur levert een gesloten vorm oplossing voor de resulterende fractionele differentiaalvergelijking, die geldig is voor elke $\alpha \in (0, 1]$.
• Verificatie: De methode wordt strikt getoetst door te laten zien dat de oplossing convergeert naar de klassieke oplossing wanneer $\alpha \to 1$.

4. Resultaten: Analyse van een RC-schakeling

Als case study wordt een RC-schakeling (weerstand $R$ en condensator $C$) onderzocht die is aangesloten op een DC-batterij.

• Vergelijking: De klassieke vergelijking $dq/dt + \Gamma q = \Gamma q_0$ (waarbij $\Gamma = 1/RC$) wordt omgezet naar een fractionele vorm met de CF-afgeleide.
• Kies van $\phi$: Voor dit specifieke probleem kiest de auteur $\phi(t, \alpha) = e^{-(1-\alpha)\Gamma t} / \Gamma$.
• Oplossing: De lading $q(t, \alpha)$ wordt afgeleid als:
  $$q(t, \alpha) = q_0 \left( 1 - e^{(1-\alpha)\Gamma t} \left( \frac{2-\alpha}{(1-\alpha) + e^{(1-\alpha)\Gamma t}} \right)^{1/(1-\alpha)} \right)$$
• Gedrag:
  • Wanneer $\alpha \to 1$, reduceert de oplossing exact tot de klassieke exponentiële laadkromme: $q(t) = q_0(1 - e^{-\Gamma t})$.
  • Voor $0 < \alpha < 1$ beschrijft de oplossing een dissipatief systeem waarbij de fractionele orde de niet-lokale energie-dissipatie (interne wrijving) modelleert.
  • Grafieken tonen aan dat de spanning over de condensator afhangt van $\alpha$, waarbij lagere waarden van $\alpha$ een langzamere of andersoortige dissipatie vertonen vergeleken met het klassieke geval.

5. Significance en Conclusie

De paper biedt een robuust wiskundig kader voor het modelleren van fysische systemen met fractionele afgeleiden die niet-singuliere kernen gebruiken. De belangrijkste betekenis ligt in het oplossen van het dimensionele probleem zonder de introductie van onmeetbare of willekeurige parameters.

• Fysische interpretatie: De methode zorgt ervoor dat alle termen in de differentiaalvergelijking consistente fysieke eenheden hebben, wat essentieel is voor de interpretatie van experimentele data.
• Generaliseerbaarheid: De auteur stelt dat deze aanpak kan worden uitgebreid naar hogere-orde systemen en andere niet-singuliere fractionele afgeleiden.
• Toepassing: Het biedt een betere basis voor het modelleren van materialen met heterogeniteit en multi-schaal fluctuaties, waar de Caputo-Fabrizio operator vaak wordt toegepast.

Kortom, González demonstreert dat door een slimme vervanging van de tijdsvariabele via een dimensiebehoudende functie, fractionele modellen zowel wiskundig consistent als fysisch interpreteerbaar kunnen worden gemaakt.