freeboundary-weak.edp

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

//  
//   calcul d'une zone saturation en eau (nappe phr閍tique)
//
verbosity=1;
real weak=1;
string com= "  avec du/dn ";
if(weak) com = " avec du/du = residu  ";
real L=10;        // longueur du domaine					   	
real Q=0.02;      // flux entrant
real K=0.5;	      //permeabilit