21.01.2015 Views

COMSOL Multiphysics™

COMSOL Multiphysics™

COMSOL Multiphysics™

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

femtime<br />

initial values are consistent (seldom the case since the initial value of the time<br />

derivative is 0). Otherwise, the solver tries to modify the initial values so that they<br />

become consistent. The value Consistent=on can be used for index-1 DAEs. Then<br />

the solver fixes the values of the differential DOFs, and solve for the initial values of<br />

the algebraic DOFs and the time derivative of the differential DOFs. The value<br />

Consistent=bweuler can be used for both index-1 and index-2 DAEs. Then the<br />

solver perturbs the initial values of all DOFs by taking a backward Euler step.<br />

For a DAE system, if Estrat=1, then the algebraic DOFs are excluded from the<br />

error norm of the time discretization error.<br />

You can suggest a size of the initial time step using the property Initialstep.<br />

By default, the solver determines whether the system is differential-algebraic by<br />

looking after zero rows or columns in the mass matrix. If you have a DAE where the<br />

mass matrix has no zero rows or columns, put Masssingular=yes.<br />

The property Maxorder gives the maximum degree of the interpolating polynomial<br />

in the BDF method.<br />

The property Maxstep puts an upper limit on the time step size.<br />

The property Out determines the output arguments and their order. The solution<br />

object fem.sol contains the output times and the corresponding solutions, see<br />

femsol. By default, the time derivatives are not stored in the solution. To store<br />

them, use the Outcomp property, see femsolver. This will also give a more accurate<br />

value in postprocessing of values interpolated in time. The output u is a matrix<br />

whose columns are the solution vectors for the output times. The output tlist is<br />

a row vector containing the output times. The output variable Stop is 0 if a<br />

complete solution was returned, 1 if a partial solution was returned, and 2 if no<br />

solution was returned. For the other outputs, see femlin.<br />

The property Stop makes it possible to return a partial solution when the time<br />

stepping fails at some point. If a failure occurs, the computed time steps are returned<br />

in sol.<br />

Use the property Stopcond to make sure the solver stops before a certain event. You<br />

provide a scalar expression that is evaluated after each time step. The time stepping<br />

is stopped if the real part of the expression is evaluated to something negative. The<br />

corresponding solution, for which the expression is negative is not returned.<br />

You can use the property Keep to tell femtime that certain quantities are constant<br />

in time, which sometimes can speed up the computation, see “Manual Control of<br />

169

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!