1.3.2 Parabolic PDE’s (Diffusion)

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

#

PDE

description

Mathematica
Maple

hand solved?

Animated?





result

time

result

time

150

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

General initial conditions

62.537

2.083

Yes

151

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Specific initial condition

26.091

1.237

Yes

Yes

152

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

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

27.383

0.931

Yes

Yes

153

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (a)

0.42

3.671

Yes

154

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

IC hat function

35.813

2.606

Yes

Yes

155

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (b)

1.236

4.116

Yes

156

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (c)


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

65.187


handled \(n=3\) case correctly.

15.515

Yes

157

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.3 (d)

545.655

40.644

Yes

158

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

domain from -1 to +1

18.891

1.61

Yes

159

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

domain from -1 to +1

62.92

1.042

Yes

Yes

160

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

with source that depends on space only (general case)

42.268

2.872

Yes

161

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

with source that depends on space only (special case)

32.815

2.741

Yes

Yes

162

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

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

71.123

4.79

Yes

163

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

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

35.558

14.278

Yes

Yes

164

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Math 4567 Exam

18.115

3.305

Yes

165

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

With source

67.959

8.844

166

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

special initial condition

3.189

(Timed out)

600.

Yes

167

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction and Euler-Cauchy Sturm-Liouville

1.816

7.57

168

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction. Using growth form reaction term

93.031

1.13

Yes

169

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction, using logistic form for reaction term

4.709

6.938

170

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction, using Aleee form for reaction term

7.192

11.236

171

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction. Haberman 2.3.8

34.92

3.22

Yes

172

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction

34.502

1.977

Yes

173

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction

0.774

2.005

174

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction

43.242

1.68

175

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion convection (general case)

3.081

2.905

Yes

176

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion convection (special case)

63.055

2.459

Yes

Yes

177

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.2 (general case)

30.562

4.834

Yes

178

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Left end zero, right end insulated, no source

29.39

2.623

Yes

179

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

One end insulated

32.45

2.755

180

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.3.7 (general case)

32.108

2.419

Yes

181

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

specific case

12.122

1.513

Yes

Yes

182

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (a)

(Timed out)

600.

39.665

183

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Both ends insulated, no source

33.1

3.24

Yes

184

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (b) (special case)

1.606

5.717

Yes

185

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (c) (special case)

32.364

19.146

Yes

186

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 2.4.1 (d)

1.387

5.056

187

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

both ends insulated

91.769

12.132

188

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

convection heat loss

32.721

1.529

189

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Pinchover and Rubinstein 6.25

42.65

13.906

190

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

external source

18.726

38.169

191

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction (general case)

91.872

1.319

Yes

Yes

192

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Pinchover and Rubinstein 6.23

64.245

24.371

193

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Pinchover and Rubinstein 6.21

63.132

5.933

194

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Diffusion Reaction. Euler-Cauchy Sturm-Liouville

2.684

38.371

Yes

195

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

convection heat loss

32.892

2.715

Yes

196

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Mixed BC

33.75

5.871

197

Diffusion in 1D

Finite domain (bar), Both ends homogeneous BC

pict

Haberman 8.2.1 (f)

63.62

7.951

198

Diffusion in 1D

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

pict

left end insulated (general case)

33.393

3.734

Yes

199

Diffusion in 1D

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

pict

left end insulated (special case)

7.485

2.226

Yes

Yes

200

Diffusion in 1D

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

pict

right end nonhomogeneous BC (general case)

66.082

3.428

Yes

201

Diffusion in 1D

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

pict

right end nonhomogeneous BC (special case)

0.01

3.907

Yes

202

Diffusion in 1D

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

pict

right end nonhomogeneous BC, special case

10.202

3.405

Yes

Yes

203

Diffusion in 1D

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

pict

convection heat loss

39.799

(Timed out)

600.

Yes

204

Diffusion in 1D

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

pict

nonhomogeneous BC

2.176

10.863

205

Diffusion in 1D

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

pict

nonhomogeneous BC

92.738

2.184

206

Diffusion in 1D

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

pict

nonhomogeneous BC

34.188

20.303

Yes

207

Diffusion in 1D

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

pict

left end BC depends on time (general case)

66.569

15.038

Yes

208

Diffusion in 1D

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

pict

left end BC depends on time (special case)

88.627

15.167

Yes

Yes

209

Diffusion in 1D

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

pict

left end BC depends on time (special case)

86.742

9.942

Yes

Yes

210

Diffusion in 1D

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

pict

Haberman 8.3.6 (special case)

95.024

25.764

Yes

Yes

211

Diffusion in 1D

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

pict

BC depends on time (special case)

12.162

5.486

212

Diffusion in 1D

Finite domain (bar), Periodic BC

pict

Periodic BC (general case). IC given

15.368

9.892

213

Diffusion in 1D

Finite domain (bar), Periodic BC

pict

Periodic BC (general case). No IC given

0.014

1.958

Yes

214

Diffusion in 1D

Finite domain (bar), Periodic BC

pict

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

0.071

2.001

Yes

215

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

both ends non-homogeneous BC (general case)

63.746

6.998

Yes

216

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

non-homogeneous BC (special case)

16.942

3.002

Yes

217

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Articolo 8.4.1 (special case)

63.188

2.931

Yes

218

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

With source that depends on space only (general case)

52.613

5.634

Yes

219

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

With source that depends on space only (special case)

94.816

3.601

Yes

Yes

220

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

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

69.048

13.045

Yes

221

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time (general case)

84.436

29.074

Yes

222

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Both ends depend on time (special case)

67.824

34.132

Yes

Yes

223

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

both ends nonhomogeneous

63.852

33.68

Yes

224

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Haberman 8.2.1 (a) (general case)

38.949

8.092

Yes

225

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Haberman 8.2.1 (d) (general solution)

54.039

6.713

Yes

226

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)

115.005

(Timed out)

600.

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 (special case)

104.862

37.183

Yes

Yes

228

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)

79.308

18.269

Yes

229

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

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

117.548

74.248

Yes

Yes

230

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Pinchover and Rubinstein 6.17

21.971

12.96

231

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

nonhomogeneous BC

111.737

31.023

232

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Haberman 8.2.2. (a)

122.354

25.014

Yes

233

Diffusion in 1D

Finite domain (bar), Both ends nonhomogeneous BC

pict

Articolo 8.4.3

5.486

18.149

234

Diffusion in 1D

Semi-infinite domain

pict

left end constant (general case)

60.035

2.283

Yes

235

Diffusion in 1D

Semi-infinite domain

pict

left end constant (special case)


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

31.691

3.029

Yes

Yes

236

Diffusion in 1D

Semi-infinite domain

pict

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

0.577

7.363

Yes

237

Diffusion in 1D

Semi-infinite domain

pict

Left end function of time (special case)

60.047

12.418

Yes

Yes

238

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous BC

60.029

2.499

239

Diffusion in 1D

Semi-infinite domain

pict

I.C. not zero


due to IC/BC not zero

13.357

3.183

240

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous BC

60.074

2.483

241

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous BC

60.071

6.834

242

Diffusion in 1D

Semi-infinite domain

pict

nonhomogeneous B.C.

60.033

3.126

243

Diffusion in 1D

Semi-infinite domain

pict

Unit triangle I.C.

61.013

4.208

244

Diffusion in 1D

Semi-infinite domain

pict

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

61.517

6.016

245

Diffusion in 1D

Semi-infinite domain

pict

Diffusion with advection

1.75

5.401

246

Diffusion in 1D

Semi-infinite domain

pict

Practice exam problem

60.108

3.071

247

Diffusion in 1D

Infinite domain

pict

Inverse exponential I.C.

0.869

3.994

Yes

248

Diffusion in 1D

Infinite domain

pict

Advection term

0.067

2.934

249

Diffusion in 1D

Infinite domain

pict

UnitBox I.C.

60.016

6.347

250

Diffusion in 1D

Infinite domain

pict

No source

0.4

6.829

251

Diffusion in 1D

Infinite domain

pict

constant as source

0.425

3.457

252

Diffusion in 1D

Infinite domain

pict

No intial conditions

0.004

0.754

Yes

253

Diffusion in 1D

Infinite domain

pict

piecewise I.C.


due to i.c. not at zero

52.493

5.05

Yes

254

Diffusion in 1D

Infinite domain

pict

Practice exam problem

0.065

1.716

255

Diffusion in 2D

Cartesian coordinates (Rectangle, Square)

pict

No source

4.249

15.159

256

Diffusion in 2D

Cartesian coordinates (Rectangle, Square)

pict

Internal source term

10.246

13.63

257

Diffusion in 2D

Cartesian coordinates (Rectangle, Square)

pict

Articolo 6.6.3

6.641

22.307

258

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

0.97


Cant get series solution

60.266

Yes

259

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

2.743


Do not understand Maple solution

21.679

Yes

Yes

260

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

no \(\theta \) dependency

0.343

6.467

261

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

no \(\theta \) dependency

1.131

1.321

262

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Haberman 8.3.5 (General solution)

19.48

5.786

Yes

263

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Specific example of the above

3.518

4.102

Yes

Yes

264

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Specific example of the above

2.722

6.169

Yes

Yes

265

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Articolo 6.9.1

4.127

3.6

266

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

Articolo 6.9.2

4.582

41.933

267

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

372.358

1.783

Yes

268

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

8.97

5.392

Yes

Yes

269

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

10.372

8.127

Yes

Yes

270

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

4.818

8.399

Yes

271

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

16.146

7.529

Yes

Yes

272

Diffusion in 2D

Polar coordinates (disk, sector, annulus)

pict

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

(Timed out)

600.

7.495

Yes

Yes

273

Diffusion in 3D

Spherical coordinates

pict

No angle dependencies

61.104


Has unresolved Laplace integrals

3.881

274

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.4 (a)

1.883

73.092

Yes

Yes

275

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.4 (b)

6.469

44.451

Yes

Yes

276

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.4 (c)

3.162

43.613

277

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.3 (a)

0.012

4.304

278

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.3 (b)

0.008

4.246

279

Diffusion in 3D

Cylinderical coordinates

pict

Haberman 7.9.3 (c)

0.012

4.19