![Page 1: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/1.jpg)
Princeton University
AIChE Centennial – CAST – Philadelphia, 2008
I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people-
Department of Chemical Engineering, PACM & MathematicsPrinceton University, Princeton, NJ 08544
NO EQUATIONS, NO VARIABLES:
An engineering approach to complex systems modeling
Or
Why I fell in love with Newton-Raphson
![Page 2: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/2.jpg)
Princeton University
![Page 3: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/3.jpg)
Princeton University
![Page 4: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/4.jpg)
Princeton University
George Stephanopoulos
Kyriakos Zygourakis
Christos Takoudis
Vassilis Matzouranis
Christos Georgakis
![Page 5: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/5.jpg)
Princeton University
Rutherford Aris Lanny Schmidt
Richard McGehee Don Aronson
Daniel Joseph L.E. “Skip” Scriven
![Page 6: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/6.jpg)
Princeton University
Two ways of taking a model f(x)=0 to the computer
t
x
xo
to
x=x*
simulation
f(x)
xx(0) x(1) x(2)
f(x*)=0
x(0)→x(1)→x(2)→…→x* same x*
algorithm dynamics: unphysicalresult (x*): meaningful
feedback:design of new trial
Extra Bonus: Stability, bifurcation, etc
simulation of ( )dx
f xdt
computer aided analysis
![Page 7: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/7.jpg)
Princeton University
Basil Nicolaenko (Nichols) Ciprian Foias James M. Hyman
![Page 8: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/8.jpg)
Princeton University
![Page 9: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/9.jpg)
Princeton University
![Page 10: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/10.jpg)
Princeton University
![Page 11: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/11.jpg)
Princeton University
“simple”
“complicated” “complex”
J. Ottino, NWU
![Page 12: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/12.jpg)
Princeton University
-1 0 1
No news – drift to Center
SELL BUYNeutral BullishBearish
Bad News Arrives
jump left --Good News Arrives
jump right +
Bad Newsarrival rate(Poisson)
Good Newsarrival rate(Poisson)
If moves to extreme rightagent buys and becomes
neutral
If moves to extreme leftagent sells and becomes
neutral
Sell rate increasesbad news
arrival rate
Buy rate increasesgood newsarrival rate
Influence of other investors
Buyrate
Sellrate
- or lemmings falling off a cliff
![Page 13: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/13.jpg)
Princeton University
SIMULATION RESULTS50000 agents, ε+=0.075, ε - = - 0.072, v0
+=v0-=20
Open loop response, g = 42
![Page 14: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/14.jpg)
Princeton University
SIMULATION RESULTS50000 agents, ε+=0.075, ε - =-0.072, v0
+ = v0- =20
Open loop response, g = 46.5
![Page 15: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/15.jpg)
Princeton University
Multiscale / Complex System Modeling
“Textbook” engineering modeling:macroscopic behavior through macroscopic models(e.g. conservation equations augmented by closures)
Alternative (and increasingly frequent) modeling situation:• Models
– at a FINE / ATOMISTIC / STOCHASTIC level– MD, KMC, BD, LB (also CPMD…)
• Desired Behavior– At a COARSER, Macroscopic Level– E.g. Conservation equations, flow, reaction-diffusion, elasticity
• Seek a bridge– Between Microscopic/Stochastic Simulation– And “Traditional, Continuum” Numerical Analysis– When closed macroscopic equations are not available in closed form
![Page 16: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/16.jpg)
Princeton University
What I told you:
Solve the equations WITHOUT writing them down.
Write “software wrappers” around “fine level” microscopic codes
Top level: all algorithms we know and love Bottom level: MD, kMC, LB, BD, heterogeneous/ discrete media,
CPMD, hybridINTERFACE:
Trade Function/Jacobian Evaluationfor “on demand” experimentation and estimation
“Equation Free” (motivated by “matrix free iterative linear algebra”)
“Variable Free” (use of diffusion maps to discover good reduction coordinates)
Think of the microscopic simulator AS AN EXPERIMENTThat you can set up and run at will
![Page 17: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/17.jpg)
Princeton University
![Page 18: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/18.jpg)
Princeton University
Descriptions
1. Detailed: vi for each variables
2. Moments
Zeroth moment: density ρ
First moment: momentum ρv
Second moment
Third Momentm3
m2
m1
m0
………….
DISTRIBUTIONS & MOMENTS
v
p(v)
vdvvp )(
dvvp )(
![Page 19: Princeton University AIChE Centennial – CAST – Philadelphia, 2008 I.G. Kevrekidis -C. W. Gear, D. Maroudas, R. Coifman and several other good people- Department](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d0d5503460f949e1c28/html5/thumbnails/19.jpg)
Princeton University