## Abstract

We develop time-domain approach for simulation of microstrip-connected extremely-high frequency (EHF) solid-state oscillators for close-range radars, including ultrashort-pulse, ultrawide-band (UWB), and noise radars. The circuits utilize high-speed GaN-based active devices such as Gunn diodes (GD) and resonant-tunneling diodes (RTD) capable of operating with enhanced power output. Microstrip interconnects produce time-delay coupling in the system that can create a complicated nonlinear dynamics of oscillations. The circuits can generate self-emerging trains of ultra-short EHF pulses emitted into an open microstrip section for further radiation. The arrays of active devices connected in either parallel (star-case) or series (ladder-case) type of circuits were simulated. Options for generation of chaotic signals in this kind of systems have been considered. An infrared-microwave (IR-EHF) oscillator linked to the resonant antenna was simulated. The oscillator consists of an RTD-driven laser diode (LD) joint to the EHF resonant antenna with a short piece of microstrip section. The oscillator can generate both the EHF pulse radiation and the EHF modulated IR pulses. Both kinds of radiation can be emitted in the free space as the trains of correlated IR-EHF radar pulses. Arrays of oscillators can be used for enhancing the power output of the system.

### Keywords

- time-domain simulations
- solid-state oscillator
- THz
- millimeter wave
- time-delay
- chaos
- distributed systems
- active devices
- Gunn diode

## 1. Introduction

Emerging demands for the EHF oscillators capable of generation of ultra-short pulses and complicated waveforms including chaotic oscillations lead to the development of new approaches to the design and analysis of oscillator systems. The EHF oscillators are of interest for numerous applications. Significant part in these applications belongs to radar systems including, particularly, close-range and noise radars, which require ultrashort-pulse, ultrawide-band, and noise oscillation sources [1, 2]. There are different kinds of the EHF oscillators ranging from microwave power tubes (klystrons, gyrotrons, backward-wave tubes, etc.) to solid-state devices (transistor-based circuits, Gunn diodes, etc.) [3, 4], of which only the latter are discussed in this chapter.

Design of oscillators and circuits is conventionally made in frequency domain. A significant contribution to the design was made by Kurokawa [5, 6] through advancing the negative resistance oscillator concepts and developing stability analysis methods. He developed the impedance approach to the analysis of oscillator systems that makes it possible to design, in particular, multi-device circuits with spatial power combining [6].

Numerous advances to the design and analysis of oscillators have been made in the following years [7, 8, 9]. Significant developments are the extension of the frequency-domain analysis for the account of nonlinear characteristics of active devices [7], the analysis of different impedance and admittance formulations in the oscillator design [8], the application of hybrid harmonic-balance approach [9], etc. A vast literature exists on the design of oscillators with frequency-domain methods.

In this line, time-domain oscillator analysis is not a common practice. In order to deliver essential information about the oscillators and their dynamics, time-domain analysis requires huge amount of numerical simulations of complicated oscillatory systems, which have to be made in a broad range of oscillator parameters.

Despite this difficulty, there are circumstances when such an analysis is a necessity, since no alternative approach can provide adequate information on the oscillator dynamics in the relevant cases. These are the cases when ultrashort-pulse, ultrawide-band, and noise oscillation signals have to be generated [1, 2]. The problem exacerbates when signals should have extremely broad frequency spectrum extended in the EHF and THz bands.

Design of generation and transmission systems for this kind of signals inherently requires the time-domain approach [10]. For passive components like antennas, valuable contributions to mathematics of time-domain modeling that concerns ultrashort-pulse and ultrawide-band signals have been made [11, 12]. The oscillators are, however, much too complicated unstable and nonlinear systems for the efficient simulations. Nonetheless, time-domain modeling is, in fact, the most meaningful approach to the design of oscillators generating this kind of signals [13, 14], though the frequency-domain methods can also be helpful [15].

A practical way of making progress in the analysis of these oscillators is to consider simplified models, which, despite their simplicity, represent essential features of real systems. An important feature of oscillators in the EHF and THz bands is their distributed character. Even though active devices and other discrete elements may be small, their assembly into an operating circuit with extended interconnects, resonators, and antenna components makes the entire system to be comparable to the radiation wavelength.

Thus, time delay arises, essentially, due to the delayed coupling between the components that makes the circuit to operate as a distributed system. Time delay leads to complicated dynamics and, often, to the dynamical chaos in nonlinear systems [16] that makes time-delay oscillators to be attractive devices for numerous applications.

A particularly useful simplification arises when making clear distinction between discrete and distributed components and defining the model where discrete units (circuit elements or blocks of elements) are joined by transmission lines (waveguides, microstrips, etc.) in a way that qualitatively represents the actual connectivity of components in the entire system. Then, discrete blocks can be simulated by local equations in time domain and the effects of transmission lines can be accounted by readily available analytic solutions of simplified wave equations. The entire system is then suitable for reasonably efficient time-domain simulations.

We applied the approach to the time-domain analysis of a range of different distributed circuits with active devices where we assumed that the transmission lines are, typically, the microstrip sections and the active devices are either the Gunn diodes or the resonant tunneling diodes in that or another circuit [17, 18, 19, 20, 21, 22, 23, 24, 25]. Transmission lines introduce time delays in the coupling between discrete units that makes the entire circuits to operate as the time-delay oscillators.

Our microstrip-based oscillator models are qualitatively different from other time-delay oscillators usually considered [4]. The difference is that, instead of using pre-defined phase delays in the feedback circuits, we consider time-delays that emerge self-consistently as a result of backward and forward EM wave propagation along the transmission lines with account of their scattering and interference with other process, making the effects particularly complicated. The systems take into account the fact that, at the frequencies of the EHF and THz range, i.e., for the millimeter and sub-millimeter waves, time-delays become unavoidable due to the extended structure of oscillator circuits.

The following Sections present overview of basic results obtained in our time-domain simulations of extended transmission line time-delay oscillators.

## 2. Kinds of circuits and forms of oscillations

We consider solid-state oscillators that can be presented as a combination of both the lumped units (lumped circuit blocks) and distributed microstrip sections (pieces of transmission lines) of different configurations. Microstrip sections provide interconnects between the lumped units and produce time-delay in the coupling between different circuit components.

The lumped units are built up of discrete active and passive devices whose interconnects within each block are of infinitesimal length as compared to the typical wavelength

The distributed sections are the pieces of microwave transmission lines (e.g., microstrip lines as representative elements or any other waveguide structures). Transmission lines (TL) provide time-delay coupling between the lumped units. Time delay appears in the coupling between the lumped units connected by any TL section because of some time needed for the EM wave propagation along the section between the units. The time delay has to be accounted in the analysis when the length

Schematics of a few circuits being considered are presented in Sections 5, 7, and 8 below. Depending on the kind of circuits, different types of nonlinear oscillation can be excited in time-delay systems. The key elements in these distributed systems are the active blocks that contain one or another kind of solid-state active devices.

We consider the EHF solid-state devices such as Gunn diodes or the resonant tunneling diodes (RTD) that can operate in a broad range of frequencies varying from, essentially, 10GHz to about 1 THz and more. The operation of these devices and oscillators is best understood in terms of the negative resistance oscillator concepts [7, 8]. Typically, we consider the Gunn diode circuits in our models, though one example of RTD system is discussed in the Section 8. In practice, the most common are the GaAs Gunn diodes but GaN devices are now of greater interest due to their potential for high-power and high-frequency operation.

The Gunn diodes are simulated using the approximation of limited space-charge accumulation (LSA) mode. In this mode, the strong-field domain in the Gunn diode is bounded to the surface electric contact and can only oscillate near the contact rather than travel through the entire structure. Then, the oscillation frequency of the Gunn diode can vary in a broad range and achieve rather high values.

Using the LSA approximation, the device operation can be described in terms of the given current-voltage characteristics with negative differential resistance (NDR) region. In this model, the current-voltage characteristics of typical Gunn diodes, e.g., GaAs diodes, can be presented in the following form [17, 18, 19, 20, 21, 22, 23, 24].

where

The Gunn diode self-excitation begins when the voltage

The formulation in terms of the current-voltage characteristics, which is typical for the LSA approximation, means that the diode is capable of instant response to any external signal. The operation of such a diode is, formally, not limited from above by any high frequency value. In reality, though, the high-frequency operation is limited by the diode intrinsic capacitance

In our models, we consider both kinds of approximations when the Gunn diode is either not limited in the oscillation frequency or, on the contrary, is characterized by intrinsic capacitance and inductance, which impose the limit on the diode operation frequency. In the latter case, the intrinsic capacitance and inductance are defined as the effective components directly connected to the Gunn diode within the lumped active unit.

Early models, for simplicity, did not account for the diode capacitance and inductance, thus, ignoring the diode frequency limit. The approximation allowed us to significantly simplify the original problems and reduce them to the forms which are more accessible for numerical simulations. In this way, we could consider self-excitation and nonlinear dynamics of the EM field oscillations in a closed two-dimensional (2D) rectangular cavity with an active wall [17, 18] and self-emergence of trains of pulses emitted from one-dimensional (1D) cavity with such a wall on the one side and a dielectric plate as a semi-transparent mirror on the other side that makes the cavity an open resonator [19]. These simulations are discussed in Section 4.

Later models took into account the intrinsic frequencies of active devices. When applied to 1D microstrip circuits and the networks of circuits with either parallel or series connections of branches, we could observe and analyze a series of new effects in self-oscillations of these circuits [20, 21, 22, 23, 24, 25]. The circuits and relevant effects are considered in Sections 5–8.

The dominant effect in these circuits, apart from conventional continuous wave (CW) generation, is the excitation of trains of short radio-frequency (RF) pulses. When using high-speed devices such as GaAs and GaN Gunn diodes in circuits with proper other parameters such as a high resonator frequency and a short length of microstrip sections, one can achieve self-excitation of short trains of ultra-short EHF pulses, which are emitted into an infinite section of another transmission line. Self-developing transitions between either the CW or pulse modes of EHF oscillations are possible, which depend on both the operation conditions and prehistory of oscillations, thus, revealing the hysteresis and bistability effects in the time-delay EHF oscillators being considered.

## 3. Mathematical models and simulation techniques

Time-domain analysis of nonlinear oscillator systems is based on computer simulations for the numerical solutions of oscillatory equations. Typically, ordinary differential equations (ODE) for lumped systems or partial differential equations (PDE) for distributed structures are in use. For more complicated cases, integral-differential, difference-differential, difference-integral-differential, and other kinds of equations may be needed for the adequate modeling of real oscillator systems. The term difference-differential means the equations that account for the finite delay in their arguments that appear in some terms of the equations. The equations of this kind are also called the equations with deviating arguments.

In our models, the delay in the electromagnetic coupling between spatially separated blocks creates multiple time delays in the evolutionary equations that can be derived for the entire oscillator system. The equations are formulated for the electromagnetic field, current, and voltage quantities as functions of time and spatial coordinates in the given circuits. Using various transformations, the equations are reduced to some other forms mentioned above to make their numerical solution more accessible.

Specific examples of equations obtained in this way for some of the models are presented in the following Sections. Typically, the equations are reduced to one or another version of a set of ODEs with time-delay arguments in their terms that describe the EM wave coupling in the distributed microstrip circuits.

For the accurate solution of these equations, we applied a highly efficient and reliable Dormand-Prince method of the (5,3) order of accuracy [26]. A publicly available software code of the method was amended with our extension that provided the dynamical storage of dense output of a solution in a long period of time in the past so as all the time-delayed values were available.

Since the past values are requested by solver at some unknown time nodes in the storage domain that may exceed the formal time-delay value, we created sufficiently large storage domain, applied polynomial interpolation of a high and controllable accuracy that uses a big domain and many nodes around the requested node, and provided a special control that no limitations are broken.

Thus, we obtained a unique, highly accurate and reliable software tool for time domain simulations of complicated nonlinear dynamics of time-delay oscillators that may exhibit dynamical chaos and other unstable transient, oscillatory, and evolutionary developments. The tool appeared to be highly efficient for solving the time-domain problems arising for our models of time-delay oscillatory systems.

More generic software tools could also be used for time-domain simulations of oscillator circuits. At present, SmartSpice simulation software is available for the engineering applications of time-domain circuit analysis [27] though applicability of generic techniques for time-delay problems and complicated circuits is limited.

The other circuit models that assumed an instant response of active devices (excluding the case of 2D cavity) lead to the difference-delay equations rather than time-delay ODEs. These equations are solved by direct iteration process.

In all the problems considered, the initial conditions were spatially uniform steady-state solutions that exist at the given initial parameters including the bias voltage

## 4. Early models, dynamical chaos, and pulses

The early models assumed an instant response of active devices. In this way, they ignored any possible frequency limitations imposed by the limited operation speed of active devices characterized by the relevant intrinsic frequency

The first model [17, 18] is presented by the 2D rectangular cavity (

The evolution of self-excitation may occur in a different manner depending on the system parameters. An essential parameter is the coupling coefficient

The main result obtained in this model is the emergence of the dynamical chaos of the EM field in the cavity when the coupling coefficient

Most of simulations were made with a square cavity (

Simulations in this 2D model were made using the series expansion of the cavity field in spatial modes with time-dependent expansion coefficients in a way, which is conceptually similar to formulations developed in [11]. In this way, a set of nonlinear ODEs was obtained for the expansion coefficients as unknown functions. The active layer

The second model [19] is formulated as a 1D open-cavity problem where the 1D cavity (

In this case, the problem formulated for the electric field

The main result obtained for these structures is that, typically, a train of the EHF pulses is self-excited in the 1D oscillator. The basic frequency of the EHF oscillations in these pulses is defined by the thickness

Characteristic frequencies of the EHF pulses emerging in these structures would be, e.g., _{2} wafer of thickness

## 5. Active circuit with a remote resonator antenna

The model that removes an essential limitation of instant response of active devices was investigated in detail in [20]. The model is rigorously formulated in terms of the 1D microstrip distributed circuit with spatially separated active and passive components, Figures 1 and 2.

The circuit consists of four basic parts which are (i) an active block

Figures 1 and 2 show a more generic version of the system considered in [20]. The system of a generic kind includes the stub of the length

The governing differential equations and boundary conditions are obtained by applying the wave equations to the transmission line sections and the Kirchhoff circuit equations to the diode and resonator blocks. The wave equations describe the voltage waves

The Kirchhoff equations define boundary conditions imposed at the contact points of both the Gunn diode and resonator circuits. They are formulated in terms of the voltage and current values at the contact points

The initial condition is imposed as the state of no oscillations when the diode voltage is set outside the NDR region by the source voltage

The time and space variables are used in relative units

The equations for this model are reduced to a set of ODEs with time delays

where

## 6. Bistability, hysteresis, and trains of the EHF pulses

Simulations of oscillators with active circuits and remote resonator antennas revealed the existence of two oscillation modes, of which one mode is the CW oscillations and the other is the trains of the EHF pulses. The emergence of one or another oscillation mode depends on the bias voltage of the Gunn diode circuit. The EHF pulse mode arises when the Gunn diode operating voltage

An important feature of the effect is the co-existence of both the pulse and the CW modes of oscillations in some range of operating voltages. The oscillator can generate either the EHF pulses or CW oscillations at the same operating voltage

The effect means bistability of the oscillation modes in the PC-domain. The kind of the mode being excited depends on the history of bias variations, i.e., on the way the oscillator is driven to the operating voltage. In this process, the oscillator reveals an hysteresis in switching between different modes. When the oscillator is driven into the PC-domain through the P-domain starting from small values, the EHF pulses are self-excited. On the contrary, when the oscillator is driven through the C-domain starting from large values, the CW oscillations are generated.

The PC-domain borders also depend on the direction of driving the operating voltage into this domain. So, there are, in fact, two kinds of P, C, and PC domains, which could be labeled, e.g., as P-up and P-down, C-up and C-down, and PC-up and PC-down domains. The effect is illustrated in Figure 3 that shows the I-V curve of the Gunn diode in relative units

Depending on the speed of driving the diode into the relevant stable domain of either the P or C kind, one can excite initially the oscillation mode, which is not intrinsic for that domain, e.g., CW oscillations in the P-domain. This mode is, however, unstable and after a certain period of time it gradually turns into the mode, which is intrinsic for the given domain, e.g., into the pulse mode in P-domain. Similarly, pulses, which could initially be excited in C-domain, gradually transform into the CW oscillations, which remain stable in this domain. The latter process is illustrated in Figure 4.

Let us now consider the properties and the mechanism of emergence of the EHF pulses. The time length of each pulse

The pulse carrier frequency

In relative units, our simulations were made, typically, at the parameter values

The emergence of one or another oscillation mode depends on the length

The formation of trains of the EHF pulses can be explained as follows. If the circuit is designed so that excitation is possible with no resonator at the site

Then, if the conditions are so that oscillations cannot exist with both the active block and the remote resonator engaged, the oscillations cease for the period of time

If the system is so that oscillations can exist in the presence of resonator, the oscillations, once appeared, would continue as a steady-state process. In this case, CW oscillations are excited whose frequency is defined by the Gunn diode and the remote resonator circuits.

The effects described above are tightly connected to a general problem of oscillation quenching and collective behavior of oscillators [29]. They are also related to bifurcations observed in square-wave switching in delay-coupled semiconductor lasers [30]. It is clear that, despite essentially different governing equations, the common feature of time-delay coupling leads to similar consequences in terms of nonlinear dynamics of microwave oscillations and optical polarization in these cases.

The generation of trains of pulses discussed above is obtained in the model that accounts for non-instant response of devices caused by the capacitance and inductance of active units [20]. These simulations confirm the conclusions of a simpler model [19] that the excitation of trains of the EHF pulses is a generic property of those oscillators, which are specified by remote location of their resonator structures.

## 7. Parallel and series connections of active circuits

One can imagine many kinds of active circuits and distributed networks of microstrip connected discrete circuits. There are two basic types of connection of active circuits, which are the parallel and series connections, Figures 6 and 7, respectively.

The parallel connection [21] is conceptually similar to the case of microwave cavity coupling of devices with spatial power combining considered by Kurokawa [6]. Spatial power combining is an important issue in this topic [31]. Our time-domain simulations of microstrip circuits connected in parallel confirm a possibility of increasing the total power output of the system proportional to the number of circuits

The series connection of active circuits was also considered and the effect of nonlinear power combining was demonstrated [22]. The series connection of circuits appears to be less promising then the parallel connection since, due to the self-consistent evolution of the entire system, the basic oscillation frequency, typically, decreases with increasing the number of devices. The active blocks in [22] were different from those in [21] that, partially, could explain the effect. Nonetheless, the increase of the total length of the system in series connection is supposed to be the main reason of reducing the basic oscillation frequency.

Turning back to the parallel connection of active circuits and keeping in mind the explanation of the effect of pulsing presented in Section 6, we can consider a network of

As a test of this possibility, we considered a system of two parallel branches of identical active circuits presented in Figure 2, though of different and, generally, non-commensurable length of time-delay microstrip sections [23]. In this case, despite an apparent simplicity of the system, we observed complicated and, in the lower frequency bands, virtually quasi-chaotic nonlinear oscillations, Figure 8. The effect is similar to the excitation of chaos in the 2D cavity [17, 18] or in the network of dispersive transmission lines where different frequency components take different times for the return of the feedback signal.

Quasi-chaotic character of the radiated wave in the low-frequency band was observed in both the Poincare sections and the auto-correlation functions of the emerging self-oscillations. When comparing auto-correlation functions of trains of the EHF pulses and quasi-chaotic signals arising under the relevant conditions, one can see a revival of correlations over the period of pulse repetition in the train of pulses and, on the contrary, a significant loss of correlation in the quasi-chaotic signal at all the times exceeding the basic period of oscillations, Figure 9.

Poincare sections plotted for the variables

The frequency spectrum of the quasi-chaotic signal shows the presence of chaotic components around the basic oscillation frequency and in the low-frequency band, Figure 11.

When the trains of the EHF pulses are excited, one can obtain rather short length of pulses. Since the process of turning on and off the EHF oscillation pulses is very short, one can obtained pulses that consist of just a few oscillations within each pulse when the microstrip sections between active circuits and resonant structures are sufficiently small, still being of nonzero length.

As an example, in a single-branch circuit of the kind shown in Figures 1 and 2, in case of

One of the oscillators of such a kind was subject to the frequency-domain analysis for the comparison with time-domain simulations [24]. The oscillator had a single active block, though of slightly different design from those above, and a stub of length

The frequency-domain analysis in the form adjusted for the radiating circuits was found to be capable of predicting small-signal oscillation spectra in either single- or multi-frequency cases [24]. At the same time, the approach, naturally, failed when strong nonlinear oscillations of more complicated character have been developed.

## 8. Resonant tunneling diode and laser diode circuit emitting trains of correlated EHF and optical pulses

Finally, we consider a time-delay version of an interesting oscillator that uses the laser diode driven by the resonant tunneling diode. The original form of the oscillator was proposed in [32, 33]. The RTD-LD circuit could generate optical (infrared) LD signals (pulses, oscillations, or chaos) when the RTD was excited by the external radio frequency (RF) bias. The authors made up the oscillator operating at the frequencies up to 2 GHz [33] and analyzed it with numerical simulations as a lumped unit.

The LD used in [33] was an optical communication laser operating at around 1550 nm IR radiation wavelength with an average output power of 5 mW. The RTD-LD hybrid circuit was produced with a minimal length of bonding wires

We considered a distributed version of the circuit that transforms the latter into an open system of the kind as shown in Figures 1 and 2 of Sections 5 and 6, where the Gunn diode is replaced by a monolithic RTD-LD unit and the remote resonator block

Under the oscillation conditions, when using appropriate values of circuit parameters and, particularly, choosing the RTD operating at sufficiently high frequencies, e.g., up to 1 THz [34, 35], one can make the circuit to generate a train of short EHF pulses

The trains of both the EHF and optical pulses can be rapidly turned on and off by the external bias digital signal applied to the RTD so as to produce relatively short pieces of trains of EHF pulses. This property, along with a possibility of synchronous generation of both the EHF and optical pulses, can make the dual EHF-optical oscillators of this kind attractive for potential applications in various forms of close range radar systems. The optical power output of these systems can be increased when making reasonably large chip-on-board LD arrays of, e.g., 20 elements or so in a way similar to conventional high-power LED arrays.

There is an optimal range of the length of microstrip section

In addition to the basic form of dual EHF-optical pulse oscillations presented above, there is another, a more interesting mode of operation. When having a sufficiently large excess in the operation speed of the RTD-LD unit as compared to the duration of separate pulses, each optical pulse can also be modulated in power by the RTD EHF current. This makes both the EHF and optical pulses carrying the same EHF signal, Figure 12.

Figure 12 shows an example of the EHF pulse radiated by the RTD-LD circuit and the relevant EHF-modulated optical pulse. Pulse oscillations have been computed when both the RTD-LD and the resonator circuit intrinsic frequencies are

Detection of the EHF modulation of optical (IR) pulses is yet another issue. Direct detection of the EHF modulation of optical pulses is expected to become possible with an RF optical heterodyne photo-detector proposed in [37].

## 9. Conclusions

We presented time-domain simulations of distributed time-delay oscillators, which show complicated nonlinear dynamics of electromagnetic oscillations generated by these systems.

We considered the models that consist of discrete lumped circuits of active and passive devices and distributed sections of microstrip lines connecting the lumped circuits. The active devices were presented by the Gunn diodes and the resonant tunneling diodes operating in the EHF frequency band. Parallel and series connections of microstrip sections with active devices have been simulated. Other structures like a 2D microwave cavity with a wall covered with active devices and a 1D open resonator made by a similar wall and a thin dielectric mirror have also been investigated.

Simulation models were developed, which rely on the method that reduces the problems for the wave equations in structures with active devices to the problems with time-delay equations of difference-differential kind. A Dormand-Prince method of the (5,3) order of accuracy for ordinary differential equations was applied and extended for solving time-delay equations arising in simulation.

Time-domain simulations revealed a diversity of dynamical effects in time-delay oscillators being considered. A possibility of chaotic or quasi-chaotic oscillations was observed in a 2D cavity with active devices and in some cases of parallel connection of microstrip sections with Gunn diode circuits.

Self-excitation of trains of the EHF pulses in microstrip structures with Gunn diodes and remote resonators have been discovered. A similar kind of trains of either the baseband or the EHF pulses were also observed in a more simplified model of a 1D cavity with an active layer and a dielectric mirror. Bistability in the generation of either the continuous waves or the trains of the EHF pulses in the Gunn diode systems with remote resonator was discovered and hysteresis in switching between these generation modes was observed.

A dual kind of the EHF-pulse and the EHF-modulated optical pulse generator using an RTD-LD time-delay oscillator has been proposed and investigated.

The approach based on the splitting of oscillator systems on discrete parts of active circuits and distributed parts of propagation sections makes it possible time-domain simulations of complex oscillatory structures. The approach has a major advantage over exact modeling with common engineering software. It allows one to address the problems of time-domain simulation of distributed systems with limited computational resources. This makes the analysis of this kind of systems feasible that, otherwise, would be out of reach with conventional simulation tools.

## Acknowledgments

This work was funded in part by NATO research projects NUKR.SFPP 985465 (G5465) and NUKR.SFPP 985395 (G5395) within 37 frames of the Science for Peace and Security (SPS) program.