Meeting Update

I am implementing the viscosity and resistivity :
The resistivity solver is similar to the talk I gave last week, which implements a scalar resistivity that depends on the local temperature. We can be more sophisticated by allowing the resistivity to be a tensor: remember, the field diffusivity across the flow velocity and along the flow velocity is different by a factor of 2. But I have not discovered yet any MHD code fully implement this complexity. The viscosity is calculated similarly. Currently, the scalar viscosity is considered. We can implement a viscosity tensor later, but it will also be slightly more complicated.
Both the viscosity and resistivity modules are currently finished, which computes the source term induced by these two physics processes at cell centers and face centers respectively. They also calculate a preferred subcycling time based on the CFL number. Although I have not decided how to implement the subcycling since as we discussed before, the energy flow induced by the resistive diffusion and viscous diffusion are likely ignored if we do subcyclilng only in the operator splitted diffusion equation. If we just compare the hydro time step with these time steps, we can get the wanted temporal resolution for energy evolution, but that will mean a slower code.
The next step is to develop a coupling method to integrate these source terms into our MHD solver. I have discussed a bit of this with Jonathan last week. Hopefully we can get some thing running this week.
Also finished editing the paper. It is ready for submission now.

Comments

No comments.