Electron emission plays a vital role in many modern technologies, from plasma medicine to heavy ion beams for fusion. An accurate theoretical model based upon the physics involved is critical to efficient operation of devices pushing the boundaries of complexity. The interactions between different electron emission mechanisms can severely alter device performance, especially when operating in extreme conditions. This dissertation studies electron emission from the perspectives of increasing geometric and physical mechanism complexities One half of this dissertation derives new relations for space-charge limited emission (SCLE) in non-planar geometries. SCLE is the maximum stable current that may be produced by electron emission before the electric field of the electrons themselves self-limits further emission. In planar devices, this is modeled by the well-established Child-Langmuir (CL) equation. The Langmuir-Blodgett (LB) equations remain the most commonly accepted theory for SCLE for cylindrical and spherical geometries after nearly a century; however, they suffer from being approximations based on a polynomial series expansion fit to a nonlinear differential equation. I derive exact, fully analytic equations for these geometries by using variational calculus to transform the differential equation into a new form that is fully and exactly solvable. This variational approach may be extended to any geometry and offers a full description of the electric field, velocity, and charge density profiles in the diode. SCLE is also an important mechanism for characterizing the operation of devices with an external magnetic field orthogonal to the electric field. This “crossed-field” problem decreases the limiting current as electrons travel longer, curved paths, effectively storing some charge in the gap (moving parallel to the emitter). At a critical magnetic field called the Hull cutoff, electron paths become so tightly curved that the circuit can no longer be completed, a condition called magnetic insulation. Crossed-field SCLE has been accurately modeled in planar devices by Lau and Christenson. Using the variational approach, I replicate their planar results and extend the calculation to cylindrical geometry, a common choice for magnetron devices. Further, I derive additional equations with simplified assumptions that, for the first time, provide an analytic description of experimental results below the Hull cutoff field. Following this I incorporate a series resistor: device resistance (or impedance) changes non-linearly with current and voltage, so I couple Ohm’s Law (OL) to all the models of crossed-field devices. For devices just below the Hull cutoff, I predict analytically and show in simulation novel bi-modal behavior, oscillating between magnetically insulated and non-insulated modes. With crossed-field device assessment, the variational calculus approach to space-charge may be used for numerous applications, including high power microwave sources, relativistic klystron devices, heavy ion beams, Hall thrusters, and plasma processing. The other half of this dissertation derives analytic theories to solve for emission current with three or more electron emission mechanisms simultaneously. In addition to the CL law, SCLE may also occur in neutral, non-vacuum diodes, modeled by the Mott-Gurney (MG) equation. These are the two limiting mechanisms I study; the other major modality of electron emission is direct electron production, the source of current in the device. Electrons are ejected when impelled by high temperature or electric field at the emission surface. These mechanisms are thermionic (or thermal) emission, modeled by the Richardson-Laue-Dushman (RLD) equation, and field emission, modeled by the Fowler-Nordheim (FN) equation, respectively. Additionally, just as I calculated the impedance of devices operating in a crossed-field configuration, all these models can be similarly coupled to OL. I derive models unifying FN, MG, and CL (with an extension linking OL, mentoring an undergraduate) and RLD, FN, and CL. These models are relevant for modern device design, especially as micro- and nano-scale devices seek to eliminate vacuum requirements and as space and military applications require higher temperature tolerances. While multi-physics models, like the ones described above, are important, the single-physics models (FN, RLD, MG, CL, OL) are still valid (and much easier to use) in their respective asymptotic limits. For example, a circuit behaves purely according to OL for very high resistances, according to MG for very high pressures, and so forth. Importantly, when devices operate in transition regions between these asymptotic limits, none of the asymptotic equations match the predictions of multi-physics models. Yet, intersections between the asymptotic equations are easily found, say for a certain set of voltage, gap distance, and pressure, CL=MG. Since both asymptotic equations give the same prediction, we may conclude that both must be inaccurate for those physical parameters! This gives rise to what I term “nexus theory:” solving two or more asymptotic equations simultaneously to rapidly and accurately predict sets of physical parameters at which multi-physics models (specifically, the physics leading to the “nexus point” parameters, points or curves at which nexus conditions are satisfied) are required for accurate device predictions. In fact, I show that multi-physics models are necessary within roughly one to two orders of magnitude from a nexus. In effect, nexus theory provides a simple, powerful tool to determine how complex a model is necessary for a particular device. Both nexus theory and multi-physics results in this dissertation have been successfully used to design devices to operate in specific transition regimes and identify the resulting device behavior.