We propose the architecture of efficiently and flexibly extensible solver system for electromagnetic wave simulations, that can load multi kinds of schemes such as Finite-Difference Time-Domain (FDTD) scheme, Finite Element Method (FEM), and a circuit simulator, with various boundary conditions in the system. Object-oriented approach is a promising method for efficient development of the flexible simulator. The primary object in the architecture is found through our object-oriented analysis as decomposed "region" from whole the simulation space. The decomposed region is considered to be the stage on which the electromagnetic fields play under the local rules. Developers who will extend the functionality of the system can add new classes inherited from the abstract classes in our design depending on the grid structure, the scheme, or the boundary processing method.

We studied a three-wave coupling process occurring in microwave power transmission (MPT) experiment in the ionospheric plasma by performing computer experiments with one-dimensional electromagnetic PIC (Particle-In-Cell) model. In order to examine the spatial variation of the coupling process, we continuously emitted intense electromagnetic wave from an antenna located at a simulation boundary. In the three-wave coupling, a low-frequency electrostatic wave is excited as the consequence of a nonlinear interaction between the forward propagating pump wave and backscattered one. In the computer experiments, low-frequency electrostatic bursts are discontinuously observed in space. The discontinuity of the electrostatic bursts is accounted for by the local electron heating due to the bursts and associated modification of the wave dispersion relation. In a case where the pump wave propagates along the geomagnetic field Bext, several bursts of Langmuir waves are observed. Since the first burst consumes a part of the pump wave energy, the pump wave is weakened and cannot trigger the three-wave coupling beyond the region where the burst occurs. Since the dispersion relation of the Langmuir wave is variable due to the local electron heating by the burst, the coupling condition eventually becomes unsatisfied and the first interaction becomes weak. Another burst of Langmuir waves is observed at a different region beyond the location of the first burst. In the case of perpendicular propagation, the upper hybrid wave, one of the mode branches of the electron cyclotron harmonic waves, is excited. Since the dispersion relation of the upper hybrid wave is less sensitive to the electron temperature, the coupling condition is not easily violated by the temperature increase. As a result, the three-wave coupling periodically takes place in time and eventually the transmission ratio of the microwaves becomes approximately 20% while almost no attenuation of the pump waves is observed after the first electrostatic burst in the parallel case.

For plasma simulations, we developed a one-dimensional (1d) Object-Oriented Particle-in-Cell code for X11-based Unix workstations (XOOPIC) by modifying the current two-dimensional version which was originally developed by PTSG (Plasma theory and simulation group) in the University of California at Berkeley. We implemented a simplified field solve and current deposition in the code. We retained three components of particle velocity, although the spatial variation for particle position and field components is limited to one dimension. To verify the function of the 1d code, we perform simulations with typical models such as the Child-Langmuir current model and electromagnetic wave propagation in plasma. In both cases, the simulation results quantitatively agree with the theory.