Add Pfem Melting application

applications/PfemMeltingApplication/custom_elements/lagrangian_vms.h

@@ -1006,20 +1006,25 @@ void CalculateLocalVelocityContribution(MatrixType& rDampingMatrix,
 
         double Temperature;
         this->EvaluateInPoint(Temperature, TEMPERATURE, rShapeFunc);
-
+        //Temperature=900.0;
         double Activation_energy;
         this->EvaluateInPoint(Activation_energy, ACTIVATION_ENERGY, rShapeFunc);
 
         double Arrhenius_coefficient;
         this->EvaluateInPoint(Arrhenius_coefficient, ARRHENIUS_COEFFICIENT, rShapeFunc);
-
+	//KRATOS_WATCH(Arrhenius_coefficient)
+	//KRATOS_WATCH(Activation_energy)
+	//KRATOS_WATCH(Temperature)
 	double R=8.31; //universal gas constant
 	//constexpr double E_over_R_polymer = 21877.25;
 	//constexpr double C_polymer = 1.93e12;
-
-	double aux_var_polymer= Arrhenius_coefficient * exp(-Activation_energy/(R*Temperature));
-
-
+	double E_over_R_polymer = Activation_energy / R;
+	//double FFFF = E_over_R_polymer/Temperature;
+	//KRATOS_WATCH(FFFF)
+	//KRATOS_WATCH(E_over_R_polymer)
+	double aux_var_polymer= Arrhenius_coefficient * std::exp(- E_over_R_polymer/Temperature) ;
+        //double aux_var_polymer= 1.0 * std::exp(- 1.0) ; 
+	//KRATOS_WATCH(aux_var_polymer)
         for (unsigned int i = 0; i < TNumNodes; ++i) // iterate over rows
         {
             for (unsigned int j = 0; j < TNumNodes; ++j) // iterate over columns

applications/PfemMeltingApplication/custom_processes/apply_laser_process.hpp

@@ -79,6 +79,10 @@ class ApplyLaserProcess : public Process
 
         mRadius = rParameters["laser_profile"]["radius"].GetDouble();
         mPower = rParameters["laser_profile"]["power"].GetDouble();
+        mDirectionx = rParameters["direction"][0].GetDouble();
+        mDirectiony = rParameters["direction"][1].GetDouble();
+        mDirectionz = rParameters["direction"][2].GetDouble();
+
 
         // Now validate agains defaults -- this also ensures no type mismatch
         //rParameters.ValidateAndAssignDefaults(default_parameters);
@@ -135,28 +139,44 @@ class ApplyLaserProcess : public Process
         Node < 3 > work_point(0, 0.0, 0.0, 0.0);
 
         array_1d<double, 3> direction;
-        direction[0] = 0.0;
-        direction[1] = 0.0;
-        direction[2] = -1.0;
+        /*direction[0] = 0.0;
+        direction[1] = 1.0;
+        direction[2] = 0.0;*/
+
+        direction[0] = mDirectionx;
+        direction[1] = mDirectiony;
+        direction[2] = mDirectionz;
+
+	//KRATOS_WATCH(direction)
         array_1d<double, 3> unitary_dir = direction * (1.0 / MathUtils<double>::Norm3(direction));
 
+        //KRATOS_WATCH(unitary_dir)
         for (size_t k = 0; k < list_of_nodes.size(); k++) {
             const array_1d<double, 3>& coords = list_of_nodes[k]->Coordinates();
             array_1d<double, 3> distance_vector;
             distance_vector[0] = coords[0]-x;
             distance_vector[1] = coords[1]-y;
             distance_vector[2] = coords[2]-z;
+            //KRATOS_WATCH(distance_vector)
+
             const double distance = MathUtils<double>::Norm3(distance_vector); // TODO: squared norm is better here
-            const double distance_projected_to_dir = distance_vector[0]*unitary_dir[0] + distance_vector[1]*unitary_dir[1] + distance_vector[2]*unitary_dir[2];
 
+            //KRATOS_WATCH(distance)  
+            const double distance_projected_to_dir = distance_vector[0]*unitary_dir[0] + distance_vector[1]*unitary_dir[1] + distance_vector[2]*unitary_dir[2];
+            //KRATOS_WATCH(distance_projected_to_dir)   
             const double distance_to_laser_axis = std::sqrt(distance*distance - distance_projected_to_dir*distance_projected_to_dir);
 
+            //KRATOS_WATCH(distance_to_laser_axis) 
+            //KRATOS_THROW_ERROR(std::logic_error, "method not implemented", ""); 
+            //KRATOS_WATCH(mRadius)
             if(distance_to_laser_axis < mRadius) {
                 double& aux= list_of_nodes[k]->FastGetSolutionStepValue(FACE_HEAT_FLUX);
                 aux =  mPower / (Globals::Pi * mRadius * mRadius);
+
+		//KRATOS_WATCH(aux)
             }
         }
-
+        //KRATOS_THROW_ERROR(std::logic_error, "method not implemented", ""); 
         KRATOS_CATCH("");
 	}
 
@@ -210,6 +230,10 @@ class ApplyLaserProcess : public Process
     array_1d<int, 3> mPositionLaserTableId;
     double mPower;
     double mRadius;
+    double mDirectionx;
+    double mDirectiony;
+    double mDirectionz;
+
     /// Assignment operator.
     ApplyLaserProcess& operator=(ApplyLaserProcess const& rOther);
 

applications/PfemMeltingApplication/custom_utilities/heat_source.h

@@ -85,7 +85,7 @@ namespace Kratos
 
             //node_it->FastGetSolutionStepValue(ARRHENIUS_VALUE)=aux_var_polymer;
 
-            node_it->FastGetSolutionStepValue(HEAT_FLUX) = (-1.0) * density * heat_of_vaporization * aux_var_polymer;
+            node_it->FastGetSolutionStepValue(HEAT_FLUX) = 0.0 * (-1.0) * density * heat_of_vaporization * aux_var_polymer;
 
             //KRATOS_WATCH(node_it->FastGetSolutionStepValue(HEAT_FLUX))