Suggested Modeling Method

Using the Graphical User Interface

1-31

1

1-32

2 Save the geometry to a model file. The model file is an M-file, so if you want to continue working using the same geometry at your next PDE Toolbox session, simply type the name of the model file at the MATLAB prompt. The pdetool GUI then starts with the model file’s solid geometry loaded. If you save the PDE problem at a later stage of the solution process, the model file also contains commands to recreate the boundary conditions, the PDE coefficients, and the mesh.

3 Move to the next step in the PDE solving process by pressing the button.

The outer boundaries of the decomposed geometry are displayed with the default boundary condition indicated. If the outer boundaries do not match the geometry of your problem, re-enter the Draw mode. You can then correct your CSG model by adding, removing or altering any of the solid objects, or change the set formula used to evaluate the CSG model.

NOTE: The set formula can only be edited while you are in the Draw mode.

If the drawing process resulted in any unwanted subdomain borders, remove them by using the Remove Subdomain Border or Remove All Subdomain Borders option from the Boundary menu.

You can now define your problem’s boundary conditions by selecting the boundary to change and open a dialog box by double-clicking on the boundary or by using the Specify Boundary Conditions. . . option from the Boundary menu.

4 Initialize the triangular mesh. Press the button or use the corresponding Mesh menu option Initialize Mesh. Normally, the mesh algorithm’s default parameters generate a good mesh. If necessary, they can be accessed using the Parameters. . . menu item.

5 If you need a finer mesh, the mesh can be refined by pressing the Refine button. Pressing the button several times causes a successive refinement of the mesh. The cost of a very fine mesh is a significant increase in the number of points where the PDE is solved and, consequently, a significant increase in the time required to compute the solution. Don’t refine unless it is required to achieve the desired accuracy. For each refinement, the number of triangles increases by a factor of four. A better way to increase the

Ω

∂

Using the Graphical User Interface

1-33 accuracy of the solution to elliptic PDE problems is to use the adaptive

solver, which refines the mesh in the areas where the estimated error of the solution is largest. See the entry on adaptmesh in Chapter 5, "Reference" for an example of how the adaptive solver can solve a Laplace equation with an accuracy that requires more than 10 times as many triangles when regular refinement is used.

6 Specify the PDE from the PDE Specification dialog box. You can access that dialog box using the PDE button or the PDE Specification. . . menu item from the PDE menu.

NOTE: This step can be performed at any time prior to solving the PDE since it is independent of the CSG model and the boundaries. If the PDE coefficients are material dependent, they are entered in the PDE mode by double-clicking on the different subdomains.

7 Solve the PDE by pressing the = button or by selecting Solve PDE from the Solve menu. If you don’t want an automatic plot of the solution, or if you want to change the way the solution is presented, you can do that from the Plot Selection dialog box prior to solving the PDE. You open the Plot Selection dialog box by pressing the button with the 3-D solution plot icon or by selecting the Parameters. . . menu item from the Plot menu.

8 Now, from here you can choose one of several alternatives:

• Export the solution and/or the mesh to MATLAB’s main workspace for further analysis.

• Visualize other properties of the solution.

• Change the PDE and recompute the solution.

• Change the mesh and recompute the solution. If you select Initialize Mesh, the mesh is initialized; if you select Refine Mesh, the current mesh is refined. From the Mesh menu, you can also jiggle the mesh and undo previous mesh changes.

• Change the boundary conditions. To return to the mode where you can se-lect boundaries, use the button or the Boundary Mode option from the Boundary menu.

Ω

∂

1

1-34

• Change the CSG model. You can re-enter the draw mode by selecting Draw Mode from the Draw menu or by clicking on one of the Draw mode icons to add another solid object. Back in the Draw mode, you are able to add, change, or delete solid objects and also to alter the set formula.

In addition to the recommended path of actions, there are a number of shortcuts, which allow you to skip over one or more steps. In general, the pdetool GUI adds the necessary steps automatically.

• If you haven’t yet defined a CSG model, and leave the Draw mode with an empty model, pdetool creates an L-shaped geometry with the default boundary condition and then proceeds to the action called for, performing all the steps necessary.

• If you are in the Draw mode and press the ∆ button to initialize the mesh, pdetool first decomposes the geometry using the current set formula and assigns the default boundary condition to the outer boundaries. After that, an initial mesh is created.

• If you press the refine button to refine the mesh before the mesh has been initialized, pdetool first initializes the mesh (and decomposes the geometry, if you were still in the Draw mode).

• If you press the = button to solve the PDE and you have not yet created a mesh, pdetool initializes a mesh before solving the PDE.

• If you select a plot type and choose to plot the solution, pdetool checks to see if there is a solution to the current PDE available. If not, pdetool first solves the current PDE. The solution is then displayed using the selected plot options.

• If you haven’t defined your PDE, pdetool solves the default PDE, which is Poisson’s equation:

-∆ u= 10

(This corresponds to the generic elliptic PDE with c = 1, a = 0, and f = 10.) For the different application modes, different default PDE settings apply.

Using the Graphical User Interface

1-35