FLOW5 - A Finite Element Code for Fluid Flow

FLOW5 uses the finite element method to solve for the steady state velocity and pressure in a fluid flow region. The fluid is presumed to be incompressible.

Files you may copy include:

• flow5.f90, the source code;
• cavity.out, the output file;
• elements.txt, the element file, suitable for input to the DISPLAY graphics program;
• uvp.txt, the UVP file, suitable for input to the DISPLAY graphics program;
• flow5.tec, a file suitable for input to the TECPLOT graphics program;

The list of routines includes:

• FLOW5 solves a steady state incompressible flow problem using finite elements.
• BSP evaluates the linear basis functions for pressure.
• DFACT factors a double precision band matrix by elimination.
• DSOLVE solves the linear system A*X=B.
• FP computes the jacobian matrix A of the Navier Stokes residual.
• FX computes the residual RES of the Navier Stokes equations.
• INTPRS interpolates the pressure at the midside nodes.
• NEWTON applies Newton iteration, seeking a solution of FX(G) = 0.
• PLOT writes solution information to two files.
• QBF evaluates a nonisoparametric quadratic basis function.
• REFBSP evaluates a reference element linear basis function.
• REFQBF evaluates a reference element quadratic basis function.
• SETBAN computes the lower band width of the Jacobian matrix.
• SETBAS evaluates the basis functions at each quadrature point.
• SETQUAD sets the abscissas and weights for quadrature in a triangle.
• SETNOD assigns numbers to the nodes and elements.
• SETXY sets the X and Y coordinates of the nodes.
• TRANS calculates the biquadratic reference element transformation.
• UVPNORM returns the "norm" of the computed finite element solution.
• UVALQ evaluates velocities and pressure at a quadrature point in an element.
• VECMAX finds the entry of largest absolute value in a vector.
• WRTEC writes out solution information for use with TECPLOT.
• XOFXSI computes X and Y given XSI and ETA coordinates.

Back to the FORTRAN software page.