// -*- c++ -*- // This file is part of ff3d - http://www.freefem.org/ff3d // Copyright (C) 2003 Laboratoire J.-L. Lions UPMC Paris // This program is free software; you can redistribute it and/or modify // it under the terms of the GNU General Public License as published by // the Free Software Foundation; either version 2, or (at your option) // any later version. // This program is distributed in the hope that it will be useful, // but WITHOUT ANY WARRANTY; without even the implied warranty of // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // GNU General Public License for more details. // You should have received a copy of the GNU General Public License // along with this program; if not, write to the Free Software Foundation, // Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. // $Id: hydrostatic.ff,v 1.1 2003/09/26 07:55:59 ArceneL Exp $ // A real-world test case inspired from Stephane Del Pino's PhD Thesis, // paragraph 4.1.2. We solve for the hydrostatic pressure inside a // cavity. vector n = (20,20,20); vector a = (-3,-2,-2); vector b = (3,2,2); scene S = pov("test.pov"); mesh M = structured(n,a,b); domain O = domain(S,outside(<1,0,0>)); solve(H) in O by M{ pde(H) -div(grad(H)) = 0; dnu(H) = 0 on <1,0,0>; H = x*y+20 on M; }; save(opendx,"hydrostatic.dx",H,M); save(medit,"hydrostatic",M); save(medit,"hydrostatic",H,M);