A scientific description of the model is available in The Lyon-Fedder-Mobarry (LFM) global MHD magnetospheric simulation code.
We do not recommend running the LFM in a stand-alone configuration. This documentation is still relevant for coupling the LFM with other models (ie. LFM-MIX, CMIT, LFM-RCM, etc).
The Lyon-Fedder-Mobarry (LFM) Code is an integrated simulation model for the global magnetosphere-ionosphere system. The heart of the model is a time-dependent, ideal MHD calculation of the state of the magnetosphere. This magnetospheric model is tightly coupled to a realistic model for the polar ionospheres and is driven by solar wind plasma and magnetic field data upwind of the calculation domain. While lacking important physical processes in both the magnetosphere and ionosphere, the model is the simplest, self-consistent first principles model possible. With current computational capabilities, it is also the only type of model capable of providing reasonable quantitative accuracy in simulating the magnetosphere-ionosphere system.
The architecture of the LFM is based on a number of important design considerations.
The outer boundary of the MHD calculation is handled as either inflow or outflow. The solar wind in the vicinity of the Earth is generally highly supersonic, so inflow is appropriate for those directions where the solar wind passes into the computational domain. The actual inflow conditions are usually described from an appropriately time shifted time series of solar wind data observed by satellite. Other regions of the outer boundary use simple outflow boundaries. The computational domain is large enough that flows have re-accelerated back to supersonic by the time they reach the edge of the grid.
The low altitude boundary of the LFM is handled by incorporating a tightly coupled model for the ionospheric electric field. The MHD calculation is stopped at some altitude between 2 and 3 RE. Field-aligned currents are calculated and mapped along dipole field lines to the ionosphere where they are used as the source term for the height-integrated potential equation. The calculated potential is then mapped back out to the MHD lower boundary where it is used to determine a boundary condition for the velocity and electric field. The algorithms used in the LFM code are described in some detail in Lyon et al. 2004. Details of the calculation of the ionospheric conductances can be found in Fedder et al. 1995.
The LFM grid is a distorted spherical grid with azimuthal symmetry about the polar axis. The polar axis points in the x-direction of the SM coordinate system (roughly toward the Sun). In the (r-?) planes the grid is distorted to place resolution where chosen by the user. The MHD grid usually covers the domain from about 30 RE upwind to 300 RE downwind of the Earth and roughly 100 RE out to the sides. The complementary ionospheric grid is a mapping of the inner surface of the MHD (magnetospheric) domain to the two polar ionospheres. This typically covers the region from the pole to 45-60 degrees latitude.
John Lyon, 6127 Wilder Lab, Dartmouth College, Hanover, NH 03755 ; (603) 646-1242 or Michael Wiltberger, HAO/NCAR, Boulder CO; (303) 497-1532
Email: firstname.lastname@example.org / email@example.com
Lyon, J.G., Numerical methods used in the Lyon-Fedder-Mobarry Global code to model the magnetosphere, Proceedings of ISSS-7, 26-31 March, 2005.
Lyon, J.G., J.A. Fedder and C.M. Mobarry. The Lyon-Fedder-Mobarry (LFM) global MHD magnetospheric simulation code, CISM JASTP Special Issue. PDF here
Fedder, J.A., S.P. Slinker, J.G. Lyon, and R.D. Elphinstone, Global numerical simulation of the growth phase and the expansion onset for substorm observed by Viking, J. Geophys. Res.,100 , 19,083, 1995.