1.3.3 Parabolic PDE’s (Diffusion)

Table 1.8: Parabolic PDE’s (Diffusion) breakdown of results. Time in seconds

#

PDE

description

Mathematica
Maple

hand solved?

Animated?





result

time

result

time

151

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

General initial conditions

87.065

1.31

Yes

152

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Specific initial condition

29.392

0.619

Yes

Yes

153

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

IC \(u=x(1-x)\)

32.123

0.611

Yes

Yes

154

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (a)

0.402

3.055

Yes

155

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

IC hat function

32.292

1.385

Yes

Yes

156

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (b)

0.963

3.697

Yes

157

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (c)


but \(n=3\) should be special case

94.078


handled \(n=3\) case correctly.

5.432

Yes

158

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (d)

20.406

26.013

Yes

159

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

domain from -1 to +1

24.972

1.119

Yes

160

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

domain from -1 to +1

92.53

0.753

Yes

Yes

161

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

with source that depends on space only (general case)

71.052

1.606

Yes

162

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

with source that depends on space only (special case)

59.534

2.285

Yes

Yes

163

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

with source that depends on space and time (general case)

69.6

3.32

Yes

164

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

with source that depends on space and time (special case)

65.981

10.046

Yes

Yes

165

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Math 4567 Exam

22.702

1.938

Yes

166

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

With source

67.131

6.367

167

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

special initial condition

2.675

(Timed out)

600.

Yes

168

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction and Euler-Cauchy Sturm-Liouville

1.479

8.987

169

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction. Using growth form reaction term

152.58

0.83

Yes

170

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction, using logistic form for reaction term

4.483

27.471

171

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction, using Aleee form for reaction term

7.125

35.232

172

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction. Haberman 2.3.8

94.774

6.826

Yes

173

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction

93.921

1.482

Yes

174

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction

0.622

1.497

175

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction

42.986

1.147

176

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion convection (general case)

2.377

2.102

Yes

177

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion convection (special case)

62.434

1.978

Yes

Yes

178

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.2 (general case)

27.698

2.112

Yes

179

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Left end zero, right end insulated, no source

29.368

2.041

Yes

180

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

One end insulated

32.017

2.24

181

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.7 (general case)

31.697

1.804

Yes

182

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

specific case

10.587

3.334

Yes

Yes

183

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (a)

(Timed out)

600.

26.393

184

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Both ends insulated, no source

32.949

2.065

Yes

185

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (b) (special case)

1.496

4.518

Yes

186

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (c) (special case)

32.276

6.848

Yes

187

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (d)

1.346

4.021

188

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

both ends insulated

91.714

4.442

189

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

convection heat loss

32.722

3.185

190

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Pinchover and Rubinstein 6.25

44.137

8.53

191

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

external source

19.439

19.326

192

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction (general case)

91.768

0.938

Yes

Yes

193

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Pinchover and Rubinstein 6.23

64.376

31.669

194

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Pinchover and Rubinstein 6.21

63.149

3.556

195

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction. Euler-Cauchy Sturm-Liouville

2.654

26.115

Yes

196

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

convection heat loss

32.75

1.936

Yes

197

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Mixed BC

33.604

4.725

198

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 8.2.1 (f)

63.638

4.279

199

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

left end insulated (general case)

33.381

4.627

Yes

200

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

left end insulated (special case)

8.047

1.614

Yes

Yes

201

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

right end nonhomogeneous BC (general case)

96.438

2.418

Yes

202

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

right end nonhomogeneous BC (special case)

0.01

3.152

Yes

203

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

right end nonhomogeneous BC, special case

12.063

2.215

Yes

Yes

204

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

convection heat loss

40.163

(Timed out)

600.

Yes

205

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

nonhomogeneous BC

2.233

9.692

206

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

nonhomogeneous BC

92.576

1.904

207

Diffusion in 1D

Finite domain (bar), left end homogeneous, right end not

pict

nonhomogeneous BC

44.428

15.72

Yes

208

Diffusion in 1D

Finite domain (bar), right end homogeneous, left end not

pict

left end BC depends on time (general case)

66.568

11.76

Yes

209

Diffusion in 1D

Finite domain (bar), right end homogeneous, left end not

pict

left end BC depends on time (special case)

89.752

15.321

Yes

Yes

210

Diffusion in 1D

Finite domain (bar), right end homogeneous, left end not

pict

left end BC depends on time (special case)

88.034

13.846

Yes

Yes

211

Diffusion in 1D

Finite domain (bar), right end homogeneous, left end not

pict

Haberman 8.3.6 (special case)

125.302

18.548

Yes

Yes

212

Diffusion in 1D

Finite domain (bar), right end homogeneous, left end not

pict

BC depends on time (special case)

14.123

5.8

213

Diffusion in 1D

Finite domain (bar), Periodic BC

pict

Periodic BC (general case). IC given

15.983

10.554

214

Diffusion in 1D

Finite domain (bar), Periodic BC

pict

Periodic BC (general case). No IC given

0.015

1.573

Yes

215

Diffusion in 1D

Finite domain (bar), Periodic BC

pict

Periodic BC (general case). Damped heat PDE. No IC given

0.029

1.459

Yes

216

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

both ends non-homogeneous BC (general case)

94.077

4.237

Yes

217

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

non-homogeneous BC (special case)

18.468

2.648

Yes

218

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Articolo 8.4.1 (special case)

93.164

2.548

Yes

219

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

With source that depends on space only (general case)

83.766

6.341

Yes

220

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

With source that depends on space only (special case)

124.882

3.432

Yes

Yes

221

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

With source that depends on space and time only (general case)

69.358

9.333

Yes

222

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time (general case)

85.769

17.345

Yes

223

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time (special case)

98.168

28.194

Yes

Yes

224

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

both ends nonhomogeneous

63.93

14.49

Yes

225

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Haberman 8.2.1 (a) (general case)

39.141

6.474

Yes

226

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Haberman 8.2.1 (d) (general solution)

82.22

5.482

Yes

227

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time with source that depends on space only (general solution)

117.169

(Timed out)

600.

Yes

228

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time with source that depends on space only (special case)

135.484

29.86

Yes

Yes

229

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time with source that depends on time and space (general solution)

80.197

16.293

Yes

230

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time with source present (special case)

148.347

63.454

Yes

Yes

231

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Pinchover and Rubinstein 6.17

26.441

11.401

232

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

nonhomogeneous BC

142.061

16.879

233

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Haberman 8.2.2. (a)

125.581

21.911

Yes

234

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Articolo 8.4.3

5.558

16.523

235

Diffusion in 1D

Semi-infinite domain

pict

left end constant (general case)

42.075

2.036

Yes

236

Diffusion in 1D

Semi-infinite domain

pict

left end constant (special case)


It fail if assumption \(x>0\) is given. A bug

34.038

2.772

Yes

Yes

237

Diffusion in 1D

Semi-infinite domain

pict

Logan p. 76. Left end general function of time (general case)

0.578

6.391

Yes

238

Diffusion in 1D

Semi-infinite domain

pict

Left end function of time (special case)

60.047

10.552

Yes

Yes

239

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous BC

41.502

2.204

240

Diffusion in 1D

Semi-infinite domain

pict

I.C. not zero


due to IC/BC not zero

14.608

2.941

241

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous BC

47.487

2.326

242

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous BC

60.072

6.131

243

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous B.C.

45.447

2.684

244

Diffusion in 1D

Semi-infinite domain

pict

Unit triangle I.C.

35.215

2.839

245

Diffusion in 1D

Semi-infinite domain

pict

I.C. not at \(t=0\)

6.613

4.622

246

Diffusion in 1D

Semi-infinite domain

pict

Diffusion with advection

1.78

4.813

247

Diffusion in 1D

Semi-infinite domain

pict

Practice exam problem

60.111

2.73

248

Diffusion in 1D

Infinite domain

pict

Inverse exponential I.C.

18.765

3.444

Yes

249

Diffusion in 1D

Infinite domain

pict

Advection term

0.072

2.516

250

Diffusion in 1D

Infinite domain

pict

UnitBox I.C.

0.08

4.32

251

Diffusion in 1D

Infinite domain

pict

No source

1.586

5.85

252

Diffusion in 1D

Infinite domain

pict

constant as source

30.487

2.873

253

Diffusion in 1D

Infinite domain

pict

No intial conditions

0.004

0.676

Yes

254

Diffusion in 1D

Infinite domain

pict

piecewise I.C.


due to i.c. not at zero

80.483

3.928

Yes

255

Diffusion in 1D

Infinite domain

pict

Practice exam problem

20.46

1.503

256

Diffusion in 2D

Cartesian coordinates (Rectangle, Square)

pict

No source

4.275

14.339

257

Diffusion in 2D

Cartesian coordinates (Rectangle, Square)

pict

Internal source term

9.849

13.191

258

Diffusion in 2D

Cartesian coordinates (Rectangle, Square)

pict

Articolo 6.6.3

6.824

20.238

259

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

no \(\theta \) dependency, insulated (General solution)

1.016


Cant get series solution

51.553

Yes

260

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

no \(\theta \) dependency, insulated (Specific solution)

2.649


Do not understand Maple solution

17.502

Yes

Yes

261

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

no \(\theta \) dependency

0.359

1.819

262

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

no \(\theta \) dependency

1.181

1.511

263

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Haberman 8.3.5 (General solution)

20.342

5.353

Yes

264

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Specific example of the above

2.901

1.148

Yes

Yes

265

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Specific example of the above

2.482

1.184

Yes

Yes

266

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Inside ring


\(n=1\) causs division by zero

3.69

(Timed out)

600.

267

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Articolo 6.9.1

3.909

0.908

268

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Articolo 6.9.2

4.548

38.487

269

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Haberman 8.2.5 with \(\theta \) dependency (General case)

380.997

0.441

Yes

270

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

With \(\theta \) dependency (Specific example)

8.155

1.062

Yes

Yes

271

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

With \(\theta \) dependency (specific example)

9.736

28.405

Yes

Yes

272

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Insulated with \(\theta \) dependency (General solution)

5.089

5.436

Yes

273

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Insulated with \(\theta \) dependency (Specific example)

15.47

9.737

Yes

Yes

274

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Insulated with \(\theta \) dependency (Specific example)

(Timed out)

600.

9.36

Yes

Yes

275

Diffusion in 3D

Spherical coordinates

pict

No angle dependencies

61.114


Has unresolved Laplace integrals

3.518

276

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.4 (a)

2.161

81.613

Yes

Yes

277

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.4 (b)

7.282

30.911

Yes

Yes

278

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.4 (c)

3.412

27.302

279

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.3 (a)

0.012

4.145

280

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.3 (b)

0.009

4.083

281

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.3 (c)

0.012

4.082