Fluid dynamics, finite element matrices from Yousef Saad (SPARSKIT2 collection). Names converted to Univ. of Florida naming convention (.mat changed to .rua). c------------------------------------------------------------------------c c FIDAP MATRICES: c c------------------------------------------------------------------------c c These matrices are extracted from the examples provided by Fidap c c The matrices are stored in Harwell-Boeing format c c For the physic of the problems, see detailed description in Fidap c c c c Matrices extracted by: Abdelkader Baggag (baggag at cs.umn.edu) c with the help of Barry Rackner (bnr at msc.edu) c c c For help or more informations, contact: c c Prof. Yousef Saad (saad at cs.umn.edu) c c c------------------------------------------------------------------------c c GENERAL INFORMATION ON THE PROBLEMS c c------------------------------------------------------------------------c c The elements used are Q2 for the velocity and P1 for the pressure c c for 2D problems and the so-called Brick elements (8 nodes) for 3D c c problems (unless specified) c c c c We are solving the fully coupled Navier-Stokes Equations (unless c c specified), including the Temperature Equation for the Heat Transfer c c Problems and the Chemical Convection-Diffusion Equations for different c c species c c The Mesh is regular structured but not uniform c c c c | Kmax c c | c c | c c | c c |__________________ Jmax c c / c c / c c / c c / c c / Imax c c------------------------------------------------------------------------c c This version is dated: June 09, 1994 c c------------------------------------------------------------------------c c ex1.mat : 2D Transient Poiseuille Flow in a Pipe c Reynolds = 1.0 c Imax = 9, Jmax = 5 c c ex2.mat : 3D Steady Couette Flow c Reynolds = 1.0 c Imax = 9, Jmax = 3 c c ex3.mat : 2D Flow Past a Cylinder in Freestream c Reynolds = 40.0 c Imax = 19, Jmax = 9 c c ex4.mat : 2D Hamel Flow c Reynolds = 0.0 c Imax = 9, Jmax = 5 c c ex5.mat : Non-newtonian Flow in a Channel (Power Law Poiseuille Flow) c Reynolds = Power Law c c ex6.mat : 2D Die Swell Problem c Capillary number = 1.0 c Imax = 43, Jmax = 15 c c ex7.mat : 2D natural Convection in a Square Enclosure c Rayleigh = 1000 c Imax = 9, Jmax = 9 c c ex8.mat : Developing Flow, Vertical Channel, Angle = 0.0 c Rayleigh = 1000 c Imax = 5, Jmax = 13 c c ex9.mat : 2D Jet Impingement Cooling c Reynolds = 0.0 ; (Jet) = 300 c Imax = 21, Jmax = 89 c c ex10.mat : 2D Flow over Multisteps in a Channel c Reynolds = 100 ; Prandtl = 0.6896 c Imax = 13, Jmax = 5 c c ex11.mat : 3D Steady Flow; Cylinder and Flat Plate Heat Exchanger c Imax = 21; Jmax = 7; Kmax = 5 c Density = 1.0; Viscosity = 0.025; SpecifiHeat = 1.0; c Conductivity = 3.760e-3 c c ex12.mat : 2D Flow Stokes Flow in Lid-Driven Wedge (28.5 Degree Angle) c Imax = 7; Jmax = 5 c Density = 1.0; Viscosity = 5.0 c c ex13.mat : Axisymmetric Flow through a Poppet Valve c Imax = 12; Jmax = 13 c Steady, Newtonian, Axi-Symmetric, Isothermal, Nonlinear c Density = 3.122e-2; Viscosity = 6.77e-4 c c ex14.mat : 2D Isothermal Seepage Flow c Imax = 23; Jmax = 23 c Density = 1000.0; Permeability = 100; Porosity = 0.5 c c ex15.mat : 2D Spin up of a Liquid in an Annulus c Imax = 7; Jmax = 7 c Transient, Nonlinear, Newtonian, Cylindrical c Density = 1.0; Viscosity = 0.02 c c ex16.mat : 3D Flow Past a Heated Obstacle c Imax = 13; Jmax = 13; Kmax = 9 c Nonlinear c Density = 1.0; Viscosity = 0.1; Conductivity = 0.1 c c ex17.mat : 3D Flow in a Rectangular Duct Around a 120 Degree Bend c Imax = 5; Jmax = 3; Kmax = 3 c Nonlinear, Turbulent c Density = 1.74e-3; Viscosity = 1.84e-4 1.0 (Mixing) c c ex18.mat : 2D Turbulent Flow Over a backward-Facing Step c Imax = 5; Jmax = 9 c Nonlinear, Turbulent (k-epsilon) c Density = 1.0; Viscosity = 2.2222e-5 c c ex19.mat : 2D Developing Pipe Flow, Turbulent k-epsilon c Imax = 5; Jmax = 5 c Nonlinear, Axi-Symmetric, Turbulent c Density = 1.0; Viscosity = 1.0e-5 c c ex20.mat : 2D Attenuation of a Surface Disturbance c Imax = 7; Jmax = 5 c Transient, Nonlinear, free c Density = 1.0; Viscosity = 0.1 c c ex21.mat : 2D Growth of a Drop from a Nozzle c Imax = 5; Jmax = 3 c Transient, Nonlinear, free, axi-symmetric c Density = 1.0; Viscosity = 1.0 c c ex22.mat : 2D Slot Coater c Imax = 6; Jmax = 3 c Steady, Nonlinear, free c Density = 1.0; Viscosity = 0.1; SurfaceTension = 1.0 c c ex23.mat : 2D Fountain Flow in Injection Molding c Imax = 4; Jmax = 3 c Nonlinear, Free c Density = 1.0; Viscosity = 2706.0 c c ex24.mat : 2D Unsymmetric Forward Roll Coating c Imax = 5; Jmax = 3 c Nonlinear, Free c Viscosity = 1.0; Density = 1.0; SurfaceTension = 10 c c ex25.mat : 2D Drop Growth Resulting from a Moving Piston c Imax = 3; Jmax = 3 c Transient, Nonlinear, Free, AxiSymmetric c Density = 1.0; Viscosity = 1.0; SurfaceTension = 10.0 c c ex26.mat : 2D surface Tension Gradient Driven Thermal Convection c Imax = 3; Jmax = 3 c Steady, Nonlinear, Free, Weakly c Density = 1.0; Viscosity = 1.0; Surface Tension = variable; c Conductivity = 0.565 c c ex27.mat : 2D Crystal Growth Simulation c Imax = 3; Jmax = 3 c Axisymmetric, Transient, Nonlinear, Free, Strongly c Density = 1.0; Viscosity = 1.0; Conductivity = 0.051282 c c ex28.mat : 2D Two Merging Liquids with one External and an interior Interface c Imax = 7; Jmax = 6 c Nonlinear, Steady, Free c Viscosity = 0.50; density = c c ex29.mat : 2D Turbulent flow in AxiSymmetric U-Bend: k-epsilon c Imax = 5; Jmax = 3 c Nonlinear, Turbulent, Axisymmetric c Density = 1.0; Viscosity = 1.0e-5 c c ex30.mat : 3D Crystal Growth: Bridgman Technique c Imax = 11; Jmax = 7; Kmax = 3; MaxPoints = 999 c Nonlinear, Strongly c Density = 1.0; Conductivity = 0.6; SpecifiHeat = 1.0 (Solid region) c density = 306; Viscosity = 4; Enthalpy = 4; Conductivity = 4 (Fluid Region) c c ex31.mat : 2D Dilute Species Deposition on a Titled heated Plate c Imax = 7; Jmax = 3 c Nonlinear, Weakly c Density = 1.0; Viscosity = 0.16; Conductivity = variable; Diffusivity = 0.15 c c ex32.mat : 2D radiation Heat Transfer in an Open Channel c Imax = 5; Jmax = 7 c Nonlinear, weakly c Density = 1.3698; Conductivity = variable; Emissivity = variable c c ex33.mat : 2D Radiation Heat Transfer in a Square Cavity c Imax = 7; Jmax = 7 c Ra = 1.6e7; Rad = 5.88; Density = 1.3698; SpecifiHeat = 0.7300 c Emissivity = 0.45--0.88 c Nonlinear, Strongly c c ex34.mat : 2D Conduction Melting c Imax = 3; Jmax = 3 c Transient, Nonlinear, Strongly c Density = 1.0; Conductivity = 0.1; Diffusivity = 1.0 (NonIsotropic) c Viscosity and Enthalpy definec on curves c c ex35.mat : 2D Turbulent Flow in a Heated Channel c Imax = 41; Jmax = 7 c Reynolds = 25960; k-epsilon c Nonlinear, Turbulent c c ex36.mat : 2D Chemical Vapor Deposition in a Horizontal Reactor c Imax = 7; Jmax = 3 c Nonlinear, Strongly c Density = variable; Viscosity = variable defined on a curve c All the other properties are variable and defined on curves c c ex37.mat : 3D Flow of Plastic in a profile extrusion Die c Imax = 7; Jmax = 9; Kmax = 5 c Nonlinear, NonNewtonian, Weakly c Viscosity = variable; Conductivity = 1.0 c c ex38.mat : Multiple Species Reaction A + 2B -----> 2C + D c Imax = 5; Jmax = 5 c axisymmetric, nonlinear, weakly c all the properties are variable c c ex39.mat : 3D Mixing Tank 2 Blades c Imax = 9; Jmax = 5; Kmax =7 c Nonlinear c Density = 1.0; Viscosity = 1.0 c c ex40.mat : 3D Die Swell Problem on a Square Die c Imax = 5; Jmax = 11; Kmax = 5 c Re = 1.0; Ca = inf c Viscosity = 1.0; Density = 0.5 c c----------------------------------------------------------------------- c-----------------------------------------------------------------------