Abstarct: Current project investigate computational modeling of electric discharge created by collisions of particles (ions, electrons, etc.) induced by air and argon. This self-consistent model has been developed to analyze the physical behaviors of high intensity Free- burning arc at atmospheric pressure. The model can be divided into three regions: (1) the cathode boundary laxyer region (2) the ionization region (3) the anode boundary laxyer region. The cathode boundary laxyer has two subsets, 1) the sheath region and 2) pre-sheath. Thin laxyer, which is, formed due to the separation of electrons from the surface of the cathode and the ions (positive) formed goes to the cathode. The beginning of the cathode spot resulting from this laxyer, as well as the pre-sheath laxyer which is plasma, is in Non-Local Thermodynamic Equilibrium (N-LTE) condition due to the degree of ionization and different temperatures. The boundary laxyer of the anode is subject to Non-Local Thermodynamic Equilibrium (N-LTE) in the form of dissipative or cinematic effects, and the anode spot is located in the center of the boundary laxyer of the anode due to the current density and high temperature. The space between the anode and cathode boundary laxyer is the space with a high degree of ionization. That is the plasma column, which is under the Local Chemical Equilibrium (LCE) conditions of the degree of ionization and is equal to the Saha equation. Because Plasma flame is part recombination. And the plasma column is in a state of thermodynamic equilibrium, the number of collisions of particles with themselves creates different amounts of energy. Due to the fact that plasma is a Newtonian fluid, the collisions of particles have a direct effect on transport coefficients, and the plasma can have different physical properties at different pressures. The effect of pressure can be seen even in molar fractions of electrons that with increasing of different pressures, electrons are separated less from the surface of the laxyers of atoms at high pressures. Electric discharge simulations are performed by finite volume analysis in ANSYS Fluent.Using User Defined Function (UDF) Air and argon gas transport coefficients have been used to calculate plasma in magneto-hydrodynamic equations (MHD).The Lorentz force which is affected by hydromagnetic waves such as the magnetic wave and the alpha wave is also written by UDF and considered in the present work.Furthermore, the radiation equation, joule heating and work function is also considered in this work using UDF in our simulation. In modeling done in this work, the stagnation point and the stagnation laxyer have been seen and we obtained the transport coefficients at different pressures, but for modeling in the present work by Fluent we only defined the transport coefficients for pressure of one atmosphere, up to a temperature of 30,000 Kelvin.Our predictions have been compared with experimental results and other numerical predictions reported by previous authors. The analysis of the results showed that our simulation was able to predict the experimental results vey well.