Exemplo n.º 1
0
        /// <summary>
        /// Constructor for the Solver
        /// </summary>
        /// <param name="sol_Tol">Convergence Residual</param>
        /// <param name="myMesh">Mesh Object</param>
        /// <param name="myBoundaryConditions">Boundary Condition object</param>
        public Solver(double sol_Tol, Mesh myMesh, int current_time_step, int total_time_steps)
        {
            this.Nodes = myMesh.Node_Array;
            this.sol_Tol = sol_Tol;
            current_dt = current_time_step;
            total_dt = total_time_steps;
            // Assigns passed in mesh object to local mesh object
            Solver_Mesh_Object = myMesh;

            // Calls the Solve() function to use TDMA
            Solve();
        }
Exemplo n.º 2
0
        static void Main(string[] args)
        {
            Console.WriteLine(" ************************************************************************ ");
            Console.WriteLine("                     Welcome to the TEModel Program                       ");
            Console.WriteLine(" ************************************************************************ ");

            double max_Node_Position_Y = 0;

            CSVHandler Get_Current_Profile = new CSVHandler();

            List<double> Current_Profile = new List<double>();

            Current_Profile = Get_Current_Profile.Read_Current_Profile();

             // Can be used to iterate over a current of 1-5 A to the TEM module
             // to construct steady-state temperature plots
            for (int i = 0; i <= Current_Profile.Count - 1; i++)
            {

                // Time in seconds
                double time = 55;

                // Timestep in seconds, set to 9999999 for Steady State Solution
                double timestep = 0.05;
                //double timestep = 99999999;

                // Required iterations, set to 1 for Steady State Solution
                int req_iter = (int)(time / timestep);
                //int req_iter = 1;

                // Current Applied
                //double current = 4;
                //double current = (double)Current_Profile[i];

                // Number of divisions specified in both the X and Y directions
                // in order to form a more dense mesh (30-x, 30-y is mesh independent)
                int divisions_X = 20;
                int divisions_Y = 20;

                // Creates object to write .CSV files for data output
                CSVHandler csv = new CSVHandler();

                // Creates the geometry object and calls the objects
                // constructor which generates a list of Layer objects
                // and a series of X and Y coordinates on which to draw
                // the control volume lines.
                TEMGeometry myGeom_Element = new TEMGeometry();

                // Creating two lists (arrays) of doubles to hold the initial x and y positions
                // of interest from the geometry file
                List<double> InitialX = new List<double>();
                List<double> InitialY = new List<double>();

                // Sets InitialX and InitialY to the array of x and y positions which are
                // output from the geometry class for the TEM (or any other imported class)
                InitialX = myGeom_Element.x_array.ToList();
                InitialY = myGeom_Element.y_array.ToList();

                // Initializes all the material names and thermophysical properties
                Material_Initializer myMaterials = new Material_Initializer();

                // Creates mesh based on the number of divisions and initial coordinates
                // in both the x and y direction.  Also requires a list of layers, materials
                // and an input current for initialization of the source terms
               // Mesh mesh = new Mesh(InitialX, InitialY, divisions_X, divisions_Y, myGeom_Element.Layer_List, myMaterials.Material_List, Current_Profile[0], timestep);

                // Writes the mesh information (to be plotted in MatLAB) to a .CSV file
              //  csv.WriteMesh(mesh.Node_Array);

                Mesh mesh = new Mesh(InitialX, InitialY, divisions_X, divisions_Y, myGeom_Element.Layer_List, myMaterials.Material_List, Current_Profile[0], timestep);

            // Finds highest Y value and assigns it to Max_Node_Position_Y
                foreach (Node node in mesh.Node_Array)
                {
                    if (node.y_Position >= max_Node_Position_Y)
                    {
                        max_Node_Position_Y = node.y_Position;
                    }
                }

                Console.WriteLine("Max Node Position (Y) : " + max_Node_Position_Y.ToString());

                List<double> Temperaturevst = new List<double>();

                foreach (Node node in mesh.Node_Array)
                {
                    node.T = 255.0;
                    node.T_Past = 255.0;
                }

            // t < req_iter
                for (int t = 0; t < Current_Profile.Count; t++)
                {
                    Console.WriteLine("Applied Current for time step " + t.ToString() + " is " + Current_Profile[t].ToString());

                    double current = ((Current_Profile[t]));

                    mesh.Set_Source_Terms(Current_Profile[t]);
                    mesh.Initialize_Influence_Coefficients(timestep);

                    // Creates object which applies the boundary conditions to the node array,
                    // boundary conditions are specified WITHIN the BoundaryCondition.cs file
                    BoundaryCondition myBoundaries = new BoundaryCondition(mesh.Node_Array, timestep);

                    // Solves the temperature field for the given time-step and current
                    Solver mySolver = new Solver(0.0001f, mesh, t, Current_Profile.Count);

                    // Writes the temperature field, and midline temperature through a BiTE element.
                    //csv.Write_Temperature_Field(mesh.Node_Array, ((t * timestep)).ToString("N2") + "_" + current);
                    //csv.Write_Mid_Field(mesh.Node_Array, ((t * timestep)).ToString("N2") + "_" + current);

                    double average_T = 0;
                    int Ceramic_Nodes = 0;

                    foreach (Node node in mesh.Node_Array)
                    {
                        node.T_Past = node.T;
                        node.res = 200;

                        if (node.y_Position > 0.0025 && node.y_Position <= 0.003 && node.Material == "Ceramic")
                        {
                            average_T = average_T + node.T;
                            Ceramic_Nodes++;
                        }
                    }

                    average_T = (average_T / (double)Ceramic_Nodes) - 273.15;

                    Console.WriteLine("Average T: " + average_T + " from an averaged value for " + Ceramic_Nodes + " nodes.");

                    Temperaturevst.Add(average_T);

                }

                csv.Write_Temperature_Versus_Time(Temperaturevst);

            }
            Console.WriteLine("Finished...");

            Console.ReadLine();
        }