Question on the proper manner to define ...

Simple question.

How should I represent a passive (i.e. gravity drainage) drain in FEFLOW?

Anyone?

You could either use a seepage face condition (h=z, max flux = 0) on all the drainage nodes, or you could use a one-dimensional discrete element with a seepage face condition on one end.

PS: An answer to your first question will follow ASAP