SIMULATION OF MAGNETOHYDRODYNAMIC SHOCK WAVE GENERATION, PROPAGATION, AND HEATING IN THE PHOTOSPHERE AND CHROMOSPHERE USING A COMPLETE ELECTRICAL CONDUCTIVITY TENSOR

An electrical conductivity tensor is used in a 1.5D magnetohydrodynamic (MHD) simulation to describe how MHD shock waves may form, propagate, and heat the photosphere and chromosphere by compression and resistive dissipation. The spatial resolution is 1 km. A train of six shock waves is generated by a sinusoidal magnetic field driver in the photosphere with a period T = 30 s, mean of 500 G, and variation of 250 G. The duration of the simulation is 200 s. Waves generated in the photosphere evolve into shock waves at a height z similar to 375 km above the photosphere. The transition of the atmosphere from weakly to strongly magnetized with increasing height causes the Pedersen resistivity eta p to increase to similar to 2000 times the Spitzer resistivity. This transition occurs over a height range of a few hundred kilometers near the temperature minimum of the initial state at z similar to 500 km. The initial state is a model atmosphere derived by Fontenla et al., plus a background magnetic field. The increase eta p is associated with an increase in the resistive heating rate Q. Shock layer thicknesses are similar to 10-20 km. They are nonzero due to the presence of resistive dissipation, so magnetization-induced resistivity plays a role in determining shock structure, and hence the compressive heating rate Q(c). At t = 200 s the solution has the following properties. Within shock layers, Q(maximum) similar to 1.4-7 erg cm(-3) s(-1), and Q(c,maximum) similar to 10-10(3) Q(maximum). Between shock waves, and at some points within shock layers, Q(c) < 0, indicating cooling by rarefaction. The integrals of Q and Q(c) over the shock wave train are F similar to 4.6 x 10(6) erg cm(-2) s(-1) and F-c similar to 1.24 x 10(9) erg cm(-2) s(-1). A method based on the thermal, mechanical, and electromagnetic energy conservation equations is presented for checking the accuracy of the numerical solution, and gaining insight into energy flow and transformation. The method can be applied to higher dimensional simulations. It is suggested that observations be performed to map out the transition region across which the transition from weakly ionized, weakly magnetized plasma to weakly ionized, strongly magnetized plasma occurs, and to correlate it with net radiative loss.