Page 296 - Process Modelling and Simulation With Finite Element Methods

Coupling Variables Revisited 283

FEMLAB Model

Case 1 of Nicmanis and Hounslow is for aggregation only, characterized by the
assignments P(v,w)=betaO, G=O, and no breakage contributions to b and d
sourcehink terms. The idealized feed is the exponential inlet condition
nin (v ) = exp ( --v ) (7.29)

Launch FEMLAB and in the Model Navigator do the following:

Select 1-D dimension, Geom 1
Select PDE modesdGeneral=+Stationary nonlinear model, weak
form (mode gl). Dependent variable nl, independent variable v>>
Multiphysics Tab. Add Geom 2, Select 2-D dimension
Select PDE modes=+General=Stationary nonlinear form (mode 82)
Dependent variable n2, independent variables v 1, v2>>

Nicmanis and Hounslow [27] suggest that the domain VE [0,2500] has a suitable
ceiling for convergence. Pull down the Options menu and set the grid to
(-50,2550) x (-50,2550) on Geom 2. Pull down the Draw menu.

Draw Mode
~~ l______l____
Geom 1 : Specify geometry; interval [0,2500]
Geom 2: Select Draw Rectangles
Set R1 = [0,2500] ~[0,2500]
Apply
OK

Pull down the options menu and select Add/Edit constants. The Add/Edit
constants dialog box appears.

_s__
Name of constant: tau Expression: 1
Name of constant: nO Expression: 1
Name of constant: beta0 Expression: 1
0
0 Name of constant: GO Expression: 1
Name of constant: v0max Expression: 2500