FEATool Multiphysics

<h2 id="distinguishing-features">Distinguishing features</h2>
<p>FEATool Multiphysics is a fully integrated physics and PDE simulation environment where the modeling process is subdivided into six steps; preprocessing (<a href="/facts/CAD/5gSXrmmN">CAD</a> and <a href="/facts/Constructive_solid_geometry/Dcs2fzyU">geometry modeling</a>), <a href="/facts/Mesh_generation/4YArjFyH">mesh and grid generation</a>, physics and PDE specification, boundary condition specification, <a href="/facts/Solver/S8LurhlM">solution</a>, and postprocessing and visualization.<a class="footnote-ref" id="fnref:8" href="#fn:8"><sup>8</sup></a>
</p>
<table><tbody><tr><th>1. Geometry Mode</th><th>2. Grid Mode</th><th>3. Multiphysics Mode</th></tr><tr><td></td><td></td><td></td></tr><tr><th>4. Boundary Mode</th><th>5. Solve Mode</th><th>6. Post Mode</th></tr><tr><td></td><td></td><td></td></tr></tbody></table>
<h3>OpenFOAM and SU2 CFD & multi-solver interfaces</h3>
<p>FEATool has introduced a multi-simulation/solver feature whereby integrated interfaces (UI) to popular open-source solvers are available. This enables several solvers to be used from a single unified GUI and CLI without requiring detailed knowledge of the syntax or peculiarities of each solver.
</p><p>The CFD solver interfaces allows fluid dynamics problems to be solved with the finite volume CFD solvers <a href="/facts/OpenFOAM/juBSzlRz">OpenFOAM</a><a class="footnote-ref" id="fnref:9" href="#fn:9"><sup>9</sup></a> and <a href="/facts/SU2_code/IGaothCJ">SU2</a>. Using the SU2 and OpenFOAM GUI<a class="footnote-ref" id="fnref:10" href="#fn:10"><sup>10</sup></a> interfaces automatically converts fluid dynamics models to compatible corresponding mesh, boundary, and control dictionary files, runs simulations, and afterwards imports and interpolates the resulting solutions back into the toolbox. In this way more advanced, larger, and parallel CFD models, for example including turbulence, can be simulated without leaving the FEATool interface.
</p>
<h3>FEniCS multiphysics solver interface</h3>
<p>Similar to the OpenFOAM and SU2 solver interfaces, FEATool also features a fully integrated interface to the <a href="/facts/FEniCS_Project/RNHF2DT9">FEniCS</a> general FEM and multiphysics solver.<a class="footnote-ref" id="fnref:11" href="#fn:11"><sup>11</sup></a> Using the FEATool-FEniCS interface, as both codes feature PDE definition languages, multiphysics problems can automatically be translated and converted to FEniCS <a href="/facts/Python_(programming_language)/YbuGqofa">Python</a> definition files, after which the FEniCS solver is called, and the resulting solution re-imported.
</p>
<h3>Fully scriptable CLI interface</h3>
<p>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 <a href="/facts/MATLAB/qPjLISCk">m-script</a> files.<a class="footnote-ref" id="fnref:12" href="#fn:12"><sup>12</sup></a> The short MATLAB script below illustrates how a complete flow around a cylinder <a href="/facts/Computational_fluid_dynamics/zuWdtNJp">computational fluid dynamics (CFD)</a> 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).<a class="footnote-ref" id="fnref:13" href="#fn:13"><sup>13</sup></a><a class="footnote-ref" id="fnref:14" href="#fn:14"><sup>14</sup></a><a class="footnote-ref" id="fnref:15" href="#fn:15"><sup>15</sup></a> Specifically, custom <a href="/facts/Partial_differential_equation/DyPsBlv7"> partial differential equations (PDE)</a> and expressions can simply be entered and evaluated as string expressions as-is, without need for further compilation or writing custom functions.<a class="footnote-ref" id="fnref:16" href="#fn:16"><sup>16</sup></a>
</p>
% 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 );

% Postprocessing and visualization.
postplot( fea, 'surfexpr', 'sqrt(u^2+v^2)', ...
               'arrowexpr', {'u' 'v'} )

p_cyl_front = evalexpr( 'p', [0.15; 0.2], fea );
p_cyl_back  = evalexpr( 'p', [0.25; 0.2], fea );
delta_p_computed  = p_cyl_front - p_cyl_back
delta_p_reference = 0.117520

<h3>External mesh generator interfaces</h3>
<p>Similar to the external solver interfaces, FEATool features built-in support for the <a href="/facts/Gmsh/sYvz0Xsd">Gmsh</a><a class="footnote-ref" id="fnref:17" href="#fn:17"><sup>17</sup></a> and Triangle<a class="footnote-ref" id="fnref:18" href="#fn:18"><sup>18</sup></a> mesh generators. If requested instead of the built-in mesh generation algorithm,<a class="footnote-ref" id="fnref:19" href="#fn:19"><sup>19</sup></a> FEATool will convert and export appropriate Gridgen2D, Gmsh, or Triangle input data files, call the mesh generators through external system calls, and re-import the resulting grids into FEATool.
</p>
<h3>Other distinguishing features</h3>
<ul><li>Stand-alone operation (without MATLAB) or can be used as a MATLAB toolbox.</li>
<li>Fully cross platform MATLAB interoperability including other <a href="/facts/Toolbox/nVqYVta7">toolboxes</a>.</li>
<li>Extensive FEM basis function library (linear and high order conforming P1-P5, non-conforming, bubble, and vector FEM discretizations).</li>
<li>Support for structured and un-structured line interval, triangles, quadrilaterals, tetrahedral, and hexahedral mesh elements.</li>
<li>28 pre-defined equations and multiphysics modes in 1D, 2D Cartesian and cylindrical coordinates, as well as full 3D.</li>
<li>Support for custom user defined <a href="/facts/PDE5_inhibitor/Ce5r1Qvl">PDE</a> equations.</li>
<li>Mesh and geometry import, export, and conversion between <a href="/facts/OpenFOAM/juBSzlRz">OpenFOAM</a>, <a href="/facts/SU2_code/IGaothCJ">SU2</a>, Dolfin/<a href="/facts/FEniCS_Project/RNHF2DT9">FEniCS</a> XML, GiD,<a class="footnote-ref" id="fnref:20" href="#fn:20"><sup>20</sup></a> <a href="/facts/Gmsh/sYvz0Xsd">Gmsh</a>, GMV,<a class="footnote-ref" id="fnref:21" href="#fn:21"><sup>21</sup></a> Triangle (PSLG), and plain <a href="/facts/ASCII/vGUI33Qu">ASCII</a> grid formats.<a class="footnote-ref" id="fnref:22" href="#fn:22"><sup>22</sup></a></li>
<li>Online postprocssing and image export with <a href="/facts/ParaView/YCvNbd0l">ParaView</a> Glance, <a href="/facts/Plotly/JBJcgC5o">Plotly</a>, and social sharing of results.</li></ul>
<h2 id="see-also">See also</h2>
<ul><li><a href="/facts/Multiphysics/kxUaIqtL">Multiphysics</a></li>
<li><a href="/facts/Computer-aided_engineering/dAsPLr5Y"> Computer-aided engineering (CAE)</a></li>
<li><a href="/facts/Continuum_mechanics/vI7LpFDk">Continuum mechanics</a></li>
<li><a href="/facts/Finite_element_method/eUcfP1vn"> Finite element method (FEM)</a></li></ul>

<h2 id="external-links">External links</h2>
<ul><li><a href="https://www.featool.com/">FEATool Multiphysics website</a></li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1"><p>"FEM Multiphysics Simulation for MATLAB!? (engineer.com)". <a href="https://www.engineering.com/Library/ArticlesPage/tabid/85/ArticleID/12636/categoryId/23/FEM-Multiphysics-Simulation-for-MATLAB.aspx" target="_blank">https://www.engineering.com/Library/ArticlesPage/tabid/85/ArticleID/12636/categoryId/23/FEM-Multiphysics-Simulation-for-MATLAB.aspx</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li>
<li id="fn:2"><p>"Engineering - FEM Multiphysics Simulation for MATLAB (engineering.com)". <a href="https://www.engineering.com/CAM/ArticleID/12636/FEM-Multiphysics-Simulation-for-MATLAB.aspx" target="_blank">https://www.engineering.com/CAM/ArticleID/12636/FEM-Multiphysics-Simulation-for-MATLAB.aspx</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li>
<li id="fn:3"><p>"Modeling the Effects of Increased Glucose Concentration on Intraocular Pressure CSURE 2014 Summer Program" (PDF). <a href="https://www.jics.utk.edu/files/images/csure-reu/2014/Su-Cope/Report.pdf" target="_blank">https://www.jics.utk.edu/files/images/csure-reu/2014/Su-Cope/Report.pdf</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li>
<li id="fn:4"><p>"Multiple Scale Modeling for Predictive Material Deformation Analysis" (PDF). <a href="https://www.math.umass.edu/sites/www.math.umass.edu/files/newsletters/2018_umass_math_newsletter_210.pdf" target="_blank">https://www.math.umass.edu/sites/www.math.umass.edu/files/newsletters/2018_umass_math_newsletter_210.pdf</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li>
<li id="fn:5"><p>"the Department of Civil Engineering at the University of Memphis CIVIL 7117 course notes". <a href="http://www.ce.memphis.edu/7117/notes/notes.html" target="_blank">http://www.ce.memphis.edu/7117/notes/notes.html</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li>
<li id="fn:6"><p>"Department of Mathematics at Lamar University course notes". <a href="http://www.math.lamar.edu/faculty/maesumi/syllabi.html" target="_blank">http://www.math.lamar.edu/faculty/maesumi/syllabi.html</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></p></li>
<li id="fn:7"><p>"Topology Optimization Modeling with MATLAB and FEATool Multiphysics". <a href="https://www.researchgate.net/project/FEATool-Multiphysics/update/5a8adc7db53d2f0bba52d896" target="_blank">https://www.researchgate.net/project/FEATool-Multiphysics/update/5a8adc7db53d2f0bba52d896</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></p></li>
<li id="fn:8"><p>"FEATool Multiphysics online documentation suite". <a href="https://www.featool.com" target="_blank">https://www.featool.com</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></p></li>
<li id="fn:9"><p>OpenCFD. "OpenFOAM® - Official home of The Open Source Computational Fluid Dynamics (CFD) Toolbox". www.openfoam.com. Archived from the original on 22 September 2016. <a href="https://www.openfoam.com/" target="_blank">https://www.openfoam.com/</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></p></li>
<li id="fn:10"><p>"Easy to Use OpenFOAM GUI and CFD Solver Integration with FEATool Multiphysics". <a href="https://www.featool.com/Easy-to-Use-OpenFOAM-GUI" target="_blank">https://www.featool.com/Easy-to-Use-OpenFOAM-GUI</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></p></li>
<li id="fn:11"><p>"The FEniCS Project page". The FEniCS Project. Archived from the original on 5 March 2023. Retrieved 28 July 2016. <a href="https://web.archive.org/web/20230305162911/https://www.fenicsproject.org/" target="_blank">https://web.archive.org/web/20230305162911/https://www.fenicsproject.org/</a> <a href="#fnref:11" class="footnote-back-ref">↩</a></p></li>
<li id="fn:12"><p>"Digital Engineering Editor's Pick: FEATool Multiphysics 1.4 (digitaleng.news)". Archived from the original on 2018-07-24. Retrieved 2018-07-23. <a href="https://web.archive.org/web/20180724032120/https://www.digitaleng.news/de/editors-pick-featool-multiphysics-1-4/" target="_blank">https://web.archive.org/web/20180724032120/https://www.digitaleng.news/de/editors-pick-featool-multiphysics-1-4/</a> <a href="#fnref:12" class="footnote-back-ref">↩</a></p></li>
<li id="fn:13"><p>De Vahl Davis, G. (1996). "Benchmark computations of laminar flow around a cylinder". Flow Simulation with High-Performance Computers II, Notes on Numerical Fluid Dynamics. 52 (3): 547–566. doi:10.1002/fld.1650030305. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:13" class="footnote-back-ref">↩</a></p></li>
<li id="fn:14"><p>On higher order methods for the stationary incompressible Navier-Stokes equations (PhD). University of Heidelberg. 1998. CiteSeerX 10.1.1.38.533. <a href="/wiki/CiteSeerX_(identifier)" target="_blank">/wiki/CiteSeerX_(identifier)</a> <a href="#fnref:14" class="footnote-back-ref">↩</a></p></li>
<li id="fn:15"><p>John, Volker; Matthies, Gunar (2001). "Higher-order finite element discretizations in a benchmark problem for incompressible flows". International Journal for Numerical Methods in Fluids. 37 (8): 885–903. Bibcode:2001IJNMF..37..885J. CiteSeerX 10.1.1.42.8087. doi:10.1002/fld.195. S2CID 15806686. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:15" class="footnote-back-ref">↩</a></p></li>
<li id="fn:16"><p>"Black-Scholes custom equation and PDE modeling tutorial". <a href="https://www.featool.com/doc/tutorials.html#ex_ce1" target="_blank">https://www.featool.com/doc/tutorials.html#ex_ce1</a> <a href="#fnref:16" class="footnote-back-ref">↩</a></p></li>
<li id="fn:17"><p>"Comparison of open source mesh generators (GiD, Gmsh, and Triangle)". 6 March 2018. <a href="https://medium.com/multiphysics/comparison-of-open-source-mesh-generators-distmesh-gmsh-and-triangle-3c1df9c08b75" target="_blank">https://medium.com/multiphysics/comparison-of-open-source-mesh-generators-distmesh-gmsh-and-triangle-3c1df9c08b75</a> <a href="#fnref:17" class="footnote-back-ref">↩</a></p></li>
<li id="fn:18"><p>Shewchuk, Jonathan Richard (1996). "Triangle: Engineering a 2D quality mesh generator and Delaunay triangulator". Applied Computational Geometry Towards Geometric Engineering. Lecture Notes in Computer Science. Vol. 1148. pp. 203–222. CiteSeerX 10.1.1.62.1901. doi:10.1007/BFb0014497. ISBN 978-3-540-61785-3. <a href="978-3-540-61785-3" target="_blank">978-3-540-61785-3</a> <a href="#fnref:18" class="footnote-back-ref">↩</a></p></li>
<li id="fn:19"><p>Persson, Per-Olof; Strang, Gilbert (2004). "A Simple Mesh Generator in MATLAB". SIAM Review. 46 (2): 329–345. Bibcode:2004SIAMR..46..329P. CiteSeerX 10.1.1.84.7905. doi:10.1137/S0036144503429121. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:19" class="footnote-back-ref">↩</a></p></li>
<li id="fn:20"><p>"GiD - The personal pre and post processor homepage". <a href="https://www.gidhome.com" target="_blank">https://www.gidhome.com</a> <a href="#fnref:20" class="footnote-back-ref">↩</a></p></li>
<li id="fn:21"><p>"GMV - The General Mesh Viewer homepage". Archived from the original on 2013-09-26. Retrieved 2018-07-23. <a href="https://web.archive.org/web/20130926190332/http://www-xdiv.lanl.gov/XCM/gmv/GMVHome.html" target="_blank">https://web.archive.org/web/20130926190332/http://www-xdiv.lanl.gov/XCM/gmv/GMVHome.html</a> <a href="#fnref:21" class="footnote-back-ref">↩</a></p></li>
<li id="fn:22"><p>"FEATool Multiphysics technical specifications". <a href="https://www.featool.com/get-featool" target="_blank">https://www.featool.com/get-featool</a> <a href="#fnref:22" class="footnote-back-ref">↩</a></p></li>
</ol>

FEATool Multiphysics open-in-new

FEATool Multiphysics