FEATool Multiphysics ("Finite Element Analysis Toolbox for Multiphysics") is a physics, finite element analysis (FEA), and partial differential equation (PDE) simulation toolbox. FEATool Multiphysics features the ability to model fully coupled heat transfer, fluid dynamics, chemical engineering, structural mechanics, fluid-structure interaction (FSI), electromagnetics, as well as user-defined and custom PDE problems in 1D, 2D (axisymmetry), or 3D, all within a graphical user interface (GUI) or optionally as script files. FEATool has been employed and used in academic research, teaching, and industrial engineering simulation contexts.
Distinguishing features
FEATool Multiphysics is a fully integrated physics and PDE simulation environment where the modeling process is subdivided into six steps; preprocessing (CAD and geometry modeling), mesh and grid generation, physics and PDE specification, boundary condition specification, solution, and postprocessing and visualization.
Easy to use GUI
The FEATool graphical user interface (GUI) shown above was designed for ease of use and minimal prerequisite knowledge. Specifically, because of the full integration of CAD and simulation tools users can switch between pre-processing, analysis, and post-processing modes to change simulation parameters, modify geometry and meshes, to solution evaluation and visualization. In this way FEATool can make up for any performance overhead by saving users time with problem set up and analysis.
OpenFOAM and SU2 CFD solver interfaces
FEATool has introduced a multi-simulation feature whereby interfaces to popular academic and open-source solvers are developed. This feature enables these solvers to be used from the FEATool GUI and CLI without detailed knowledge of the syntax or peculiarities of each solver.
The CFD solver interfaces allows fluid dynamics problems to be solved with the finite volume CFD solvers OpenFOAM and SU2. Using the interfaces automatically converts incompressible Navier-Stokes FEATool models to compatible OpenFOAM/SU2 mesh, boundary, and control dictionary files, runs simulations, and afterwards imports and interpolates the resulting solutions back into FEATool. In this way more advanced, larger, and parallel CFD models, for example including turbulence, can be simulated without leaving the FEATool interface.
FEniCS multiphysics solver interface
Similar to the OpenFOAM and SU2 solver interfaces, FEATool also features a fully integrated interface to the FEniCS general FEM and multiphysics solver. Using the FEATool-FEniCS interface, as both codes feature PDE definition languages, multiphysics problems can automatically be translated and converted to FEniCS Python definition files, after which system calls to the FEniCS solver is made, and the resulting solution re-imported.
Fully scriptable CLI interface
GUI operation is recorded as equivalent function calls, and therefore in addition to binary formats, FEATool simulation models can also be saved and exported as fully scriptable and editable MATLAB compatible m-script files. The short MATLAB script below illustrates how a complete flow around a cylinder computational fluid dynamics (CFD) benchmark problem can be defined and solved with the FEATool m-script functions (including geometry, grid generation, problem definition, solving, and postprocessing all in a few lines of code). Specifically, custom partial differential equations (PDE) and expressions can simply be entered and evaluated as string expressions as-is, without need for further compilation or writing custom functions.
% Geometry and mesh generation. fea.sdim = { 'x' 'y' }; fea.geom.objects = { gobj_rectangle( 0, 2.2, 0, 0.41, 'R1' ), ... gobj_circle( [0.2 0.2], 0.05, 'C1' ) }; fea = geom_apply_formula( fea, 'R1-C1' ); fea.grid = gridgen( fea, 'hmax', 0.02 );
% Problem definition (incompressible Navier-Stokes equations multiphysics mode). fea = addphys( fea, @navierstokes );
% Prescribe fluid viscosity (density is default 1). fea.phys.ns.eqn.coef{ 2,end} = { 0.001 };
% Boundary conditions (Non-specified boundaries are
% per default prescribed no-slip zero velocity walls).
% Inflow (bc type 2) at boundary 4. fea.phys.ns.bdr.sel(4) = 2;
% Outflow (bc type 3, zero pressure) at boundary 2. fea.phys.ns.bdr.sel(2) = 3;
% Parabolic inflow profile x-velocity expression. fea.phys.ns.bdr.coef{ 2,end}{ 1,4} = '4*0.3*y*(0.41-y)/0.41^2';
% Check, parse, and solve problem. fea = parsephys( fea ); fea = parseprob( fea ); fea.sol.u = solvestat( fea );
% Alternatively solve with OpenFOAM or SU2
% fea.sol.u = openfoam( fea ); % fea.sol.u = su2( fea );