Norsk English

Ikke-lineære løsningsmetoder og tidsintegrasjon

Robuste løsningsmetoder er avgjørende når komplekse modeller skal gi pålitelige resultater innen rimelig tid. Vi utvikler algoritmer for ikke-lineære ligningssystemer og tidsavhengige problemer, med vekt på konvergens, nøyaktighet og beregningseffektivitet. Vi hjelper med å forstå hvorfor simuleringer blir krevende, forbedre løsningsstrategier og tilpasse metodene til modellens struktur.

Kontaktpersoner

Hva kan vi hjelpe deg med?

Simuleringer kan stoppe opp eller bruke uforholdsmessig lang tid når modellene er sterkt ikke-lineære, fysiske prosesser virker på ulike tidsskalaer eller koblingene mellom delproblemer blir vanskelige å håndtere. Valg av løsningsmetode og tidssteg er da avgjørende for både robusthet og effektivitet.

Vi har særlig ekspertise i ikke-lineære løsere for stive og sterkt koblede problemer. Arbeidet omfatter Newton-baserte metoder og trust-region-strategier, samt metoder som lokaliserer beregningsarbeidet til de delene av modellen der det trengs mest. Ikke-lineær områdedekomponering og reordning av ukjente gjør det mulig å utnytte modellens struktur og redusere behovet for kostbare globale iterasjoner.

Vi kan hjelpe med å:

  • Undersøke konvergensproblemer og forbedre robustheten i eksisterende simulatorer.
  • Utvikle og tilpasse ikke-lineære løsere som utnytter modellens struktur.
  • Tilpasse tidssteg og løsningsstrategier til utviklingen i simuleringen.
  • Koble fysiske delmodeller gjennom fullt implisitte metoder, sekvensielle formuleringer og operatorsplitting.
  • Vurdere nøyaktighet og beregningskostnad og implementere metodene i effektiv programvare.

Ta kontakt med vår forskningsgruppe for Anvendt beregningsvitenskap dersom du har en modell som er vanskelig å løse, ønsker raskere simuleringer eller trenger nye metoder i eksisterende programvare.

Robust ikke-lineær løsning

Newton-metoder kan gi rask konvergens når iterasjonene er nær løsningen, men kan også stagnere, oscillere eller gi oppdateringer som fører beregningen inn i fysisk eller numerisk problematiske tilstander. Vi arbeider med hvordan modellformulering, initialverdier, skalering og kontroll av oppdateringene kan forbedre løsningsprosessen.

Dette omfatter adaptive relaksasjonsstrategier og trust-region-metoder, der oppdateringene begrenses ut fra hvor den lokale tilnærmingen til problemet er egnet. Vi har blant annet utviklet metoder som utnytter kunnskap om fluksfunksjoner og lokale koblinger for å styre Newton-iterasjonene i flerfaseproblemer.

Vi utvikler også konvergensovervåking som oppdager ugunstig utvikling tidlig. Informasjon fra iterasjonene kan brukes til å justere løsningsstrategien eller avbryte et lite lovende forsøk og prøve på nytt med et kortere tidssteg.

Reordning og lokale løsningsmetoder

Mange diskretiserte modeller har en avhengighetsstruktur som kan utnyttes til å løse ligningene i en gunstig rekkefølge. For transportproblemer kan denne strukturen ofte bestemmes fra strømningsretningen mellom beregningscellene.

Vi har utviklet reordnede ikke-lineære Gauss–Seidel-metoder som følger disse avhengighetene. Når systemet kan ordnes trekantet, kan det løses ved en sekvens av lokale ikke-lineære løsninger, uten gjentatte globale iterasjoner. Gjensidige avhengigheter håndteres ved å samle ukjente i koblede blokker.

Reordningen gjør det mulig å konsentrere det ikke-lineære løsningsarbeidet der det faktisk trengs. Vi har kombinert slike metoder med høyere ordens diskretisering, lokal tidssteging og adaptiv grovinndeling for effektiv simulering av transport i porøse medier.

Ikke-lineær områdedekomponering og prekondisjonering

Konvergensproblemer kan være konsentrert i små deler av en stor modell. Likevel kan de tvinge en global løser til å gjøre mange kostbare iterasjoner over hele beregningsområdet.

Ved ikke-lineær områdedekomponering deles modellen i delområder som behandles med lokale ikke-lineære løsere. Disse løsningene kan brukes som prekondisjonering av det globale problemet, slik at vanskelige lokale forhold håndteres før den samlede modellen korrigeres.

Vi utvikler og undersøker blant annet Schwarz-baserte metoder, adaptive kombinasjoner av Newton-iterasjoner og ikke-lineær prekondisjonering, og strategier for kobling mellom delområder. Arbeidet omfatter både robusthet og effektiv gjennomføring i parallelle simuleringer.

Tidsintegrasjon og adaptiv tidssteging

Et tidssteg må være langt nok til at simuleringen gjennomføres effektivt, men kort nok til å beskrive viktige endringer og gi et håndterbart løsningsproblem. Numerisk stabilitet, tidsdiskretiseringsfeil og ikke-lineær konvergens stiller forskjellige krav som må vurderes i sammenheng.

Vi arbeider med implisitt tidsintegrasjon, adaptiv tidssteging, lokal tidssteging og feilkontroll. Målet er å tilpasse beregningsinnsatsen til hvor og når modellen endrer seg, og å unngå unødvendig nøyaktig løsning av delproblemer når andre feilkilder dominerer.

Vår faglige bakgrunn omfatter også forskning på eksponensielle integratorer og Lie-gruppeintegratorer. Slike metoder utnytter henholdsvis eksponentialfunksjoner av operatorer og geometrisk struktur i differensialligningene. Denne erfaringen inngår i et bredt grunnlag for å vurdere og utvikle tidsintegrasjonsmetoder for ulike typer problemer.

Løsningsstrategier for koblede modeller

Når flere fysiske prosesser virker sammen, kan alle ukjente løses samtidig eller problemet deles i mindre deler. Fullt implisitte metoder behandler koblingene samlet, mens sekvensielle metoder og operatorsplitting kan gi enklere delproblemer og større fleksibilitet.

Vi undersøker hvordan styrken i koblingene, valg av tidssteg og egenskapene til de lokale løserne påvirker nøyaktighet og konvergens. Erfaringen omfatter blant annet koblede strømnings-, transport-, varme- og geomekaniske modeller.

Gjennom samarbeid med industri og programvareleverandører har vi erfaring med å utvikle og teste metodene på krevende simuleringsmodeller. Vi vurderer hele løsningsforløpet: beregningstid, antall iterasjoner, forkastede tidssteg og kvaliteten på resultatene.

Utvalgte publikasjoner

Flere av arbeidene nedenfor har reservoarsimulering som anvendelse. De angriper imidlertid utfordringer som også oppstår i andre simuleringsområder: sterke ikke-lineariteter, lokale konvergensproblemer, koblede prosesser og ulike tidsskalaer. Arbeidene om eksponensielle integratorer viser en bredere bakgrunn innen numerisk tidsintegrasjon.