|
7 | 7 | "outputs": [],
|
8 | 8 | "source": [
|
9 | 9 | "%load_ext autoreload\n",
|
10 |
| - "%autoreload 2\n", |
11 |
| - "%matplotlib notebook\n" |
| 10 | + "%autoreload 2" |
12 | 11 | ]
|
13 | 12 | },
|
14 | 13 | {
|
|
29 | 28 | "from josie.bc import Dirichlet, Neumann, NeumannDirichlet\n",
|
30 | 29 | "from josie.geom import Line\n",
|
31 | 30 | "from josie.solver.euler.eos import PerfectGas\n",
|
32 |
| - "from josie.solver.euler.problem import Problem\n", |
33 |
| - "from josie.solver.euler.schemes import Rusanov\n", |
34 | 31 | "from josie.solver.euler.solver import EulerSolver\n",
|
35 | 32 | "from josie.solver.euler.state import Q\n",
|
36 | 33 | "\n",
|
|
41 | 38 | "right = Line([0.25, 0], [.25, 1])\n",
|
42 | 39 | "top = Line([0, 1], [0.25, 1])\n",
|
43 | 40 | "\n",
|
44 |
| - "\n", |
| 41 | + "# EOS\n", |
| 42 | + "eos = PerfectGas(gamma=1.4)\n", |
45 | 43 | "\n",
|
46 | 44 | "# Params\n",
|
47 | 45 | "JET_CENTER = 0.5\n",
|
|
105 | 103 | "metadata": {},
|
106 | 104 | "outputs": [],
|
107 | 105 | "source": [
|
108 |
| - "# Solver\n", |
109 |
| - "eos = PerfectGas(gamma=1.4)\n", |
110 |
| - "scheme = Rusanov(eos)\n", |
| 106 | + "# Scheme\n", |
| 107 | + "from josie.solver.euler.schemes import Rusanov\n", |
| 108 | + "from josie.solver.scheme.time import ExplicitEuler\n", |
| 109 | + "\n", |
| 110 | + "class MyScheme(Rusanov, ExplicitEuler): \n", |
| 111 | + " pass\n", |
| 112 | + "\n", |
| 113 | + "scheme = MyScheme(eos)\n", |
| 114 | + "\n", |
| 115 | + "\n", |
| 116 | + "\n", |
111 | 117 | "solver = EulerSolver(mesh, scheme)\n",
|
112 | 118 | "\n",
|
113 | 119 | "\n",
|
|
143 | 149 | " dt = scheme.CFL(\n",
|
144 | 150 | " solver.values, \n",
|
145 | 151 | " solver.mesh.volumes, \n",
|
146 |
| - " solver.mesh.normals,\n", |
147 | 152 | " solver.mesh.surfaces,\n",
|
148 | 153 | " CFL,\n",
|
149 | 154 | " )\n",
|
|
0 commit comments