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.

elcplproj<br />

All map types available for elcplextr can be used also in projection coupling<br />

variables. The common combination of a unit source map and a linear<br />

destination map is not very useful, though. Instead, there is a map type projection<br />

which specifies srcdim+1 vertices in the source geometry and srcdim vertices in the<br />

destination. The basis defined by vectors from the first source vertex to each of the<br />

remaining vertices is mapped onto a right-handed orthogonal system with unit axes.<br />

The basis described by the destination vertices is then mapped onto the first<br />

srcdim-1 dimensions of the same orthogonal basis. This means that the direction<br />

of integration is effectively from the first source vertex to the last.<br />

In addition to the fields present in the elcplextr element, the elcplproj requires<br />

an integration order for the line integrals evaluated for each destination point. The<br />

iorder field specifies the order of polynomials that should be exactly integrated.<br />

Cautionary<br />

Projection coupling is only implemented for simplex meshes. When finding<br />

integration limits, the elcplproj element works directly on the basic polyhedral<br />

mesh. Therefore, results can be inaccurate if the mesh does not properly resolve the<br />

geometry.<br />

Parameter or time dependency in the source transformation is not properly detected<br />

by the solvers, which means that the source transformation will not be updated<br />

between parameter or time steps in that case. Solution dependencies in the<br />

transformation are properly detected, but do not give any Jacobian contributions<br />

from the transformation.<br />

The automatic detection of nonlinear and time or parameter dependent problems<br />

does not work properly in that all problems containing projection coupling variables<br />

are considered to be nonlinear and time dependent.<br />

Examples<br />

Project the diagonal cross section distance on the left and bottom edges of a square.<br />

clear fem<br />

fem.geom = square2;<br />

fem.mesh = meshinit(fem);<br />

fem.elem = {};<br />

el.elem = 'elcplproj';<br />

el.g = {'1'};<br />

el.var = {'d','d'};<br />

clear map1;<br />

map1.type = 'projection';<br />

map1.sg = '1';<br />

map1.sv = {'1','2','3'};<br />

64 | CHAPTER 1: COMMAND REFERENCE

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

Saved successfully!

Ooh no, something went wrong!