The featured projects presented below reflect our experience in custom software development and numerical simulations for industry and academia. As software creators we are happy when our programs are used by many people and help to understand complex phenomena in various fields of research.

**Client:**Software development company, Switzerland**Technology:**C++, Matlab**Man hours**: 400**Scope:**The purpose of the project was to develop a C++ code for simulations of a dense gas dispersion over complex terrains. The code is based on shallow layer approach developed in series of papers by Hankin et al [J. Hazard. Mater, A66, p. 211–261, 1999]. The model takes into account the effects of a gas cloud convection and buoyancy, the effect of ground slopes, surface and turbulent shear stresses, top entrainment and the leading edge terms that account for interaction among dense and ambient ﬂuid. The shallow layer equations describing a gas cloud parameters (height, density, velocities) are solved by the flux-corrected transport algorithm of Zalesak, using corner transport scheme as a low order scheme and Zalesak fourth order scheme as a high order scheme. The code is implemented as a C++ library to be incorporated into the existing risk assessment software.**Results:**A typical simulation of a dense gas cloud propagation in a mountain region is presented below.

**Client:**Software development company, Switzerland**Technology:**C++, Matlab**Man hours**: 120**Scope:**The purpose of the project was to develop a C++ code for reconstruction of a 3D wind field from meteorological data in a given region with complex topography (mountains). We applied a modified version of a diagnostic wind model suggested in the paper by Ishikawa [J. of Applied Meteorology, V. 33, p. 733–743, 1994]. The approach ensures a divergence free wind field that has zero normal component at the ground surface and best fits the measured or the given wind parameters. Since a finite-difference discretization of equations leads to rather cumbersome matrix elements we have used Mathematica® to generate some parts of C++ sources (mathematical expressions) automatically to ensure that no errors are present in the code. Huge sparse linear system (matrix size ~ 200000 x 200000) is solved either by open-source (SuperLU) or commercial (Pardiso) library depending on compilation settings. The code is implemented as a C++ library to be incorporated into the existing risk assessment software.**Results:**Streamlines of typical reconstructed and adjusted wind fields are presented below.

**Client:**Graz University of Technology, Austria**Technology:**C/C++, Fortran, Matlab, Mathematica**Man hours**: 2400**Scope:**The purpose of the project was to develop a kinetic linear model of the interaction of helical magnetic perturbations with cylindrical plasmas. The model is then implemented in the numerical code named KiLCA (Kinetic Linear Cylindrical Approximation). For a given equilibrium plasma configuration, KiLCA computes nonlocal plasma conductivity operator and solves Maxwell equations to calculate magnetic perturbation amplitudes inside the plasma. Derived quantities like the dissipated power density and the torques acting on the plasma are evaluated in postprocessing routines within KiLCA. A number of new and important results have been obtained by application of KiLCA to TEXTOR, JET and ITER tokamaks. Those results are well known by the fusion plasma community and are published in top level scientific journals.**Results:**Find attached the recently published papers (+ see references therein).

**Client:**Graz University of Technology, Austria**Technology:**C/C++, Fortran, Matlab, Mathematica**Man hours**: 320**Scope:**The purpose of the project was to upgrade the existing KiLCA code taking into consideration (i) magnetic drifts of particles and (ii) energy conservation in particle collisions. The account of the magnetic drifts in nonlocal plasma conductivity operator is important near the so called "double null" point where a null of background electric field gets close to a resonant magnetic surface.**Results:**The improved plasma response model is used in self-consistent balance simulations to investigate the impact of external helical magnetic perturbations on evolution of plasma parameters in tokamaks.

**Client:**Private person, Italy**Technology:**Fortran, Matlab**Man hours**: 160**Scope:**Laser ablation is used for a growing number of applications, such as pulsed laser deposition, nanoparticle manufacturing, micro-machining, surgery, chemical analysis and many others. Despite of active research for the last decade, the exact mechanisms are not yet fully understood due to many complex phenomena that occur and interact during the laser ablation process. Following a model proposed by A. Bogaerts et al. [Spectrochimica Acta Part B, 58, 1867–1893, 2003] we have developed a code that models the formation of a plasma near the target surface by numerical solving of the coupled Euler and Saha equations. The plasma code calculates the laser beam absorption and is coupled with an ANSYS solver that simulates the heat transfer in the target.**Results:**For high laser irradiances the plasma absorbs considerable amount of the beam energy and its formation has drastic impact on the laser ablation process.

**Client:**Private person, Canada**Technology:**Fortran, Matlab**Man hours**: 100**Scope:**Abyssal ocean currents, being denser than their surrounding fluid, are largely driven by gravity and steered by their local bathymetry. Over large scales they are also strongly influenced by the Coriolis force. The dynamics of these currents becomes complicated close to the equator, where geostrophic balance breaks down as the locally vertical component of the Earth’s rotation vector changes sign. Using a conventional approach based on the shallow water equations with the topography and the Coriolis force terms we have developed a code that models the cross-equatorial abyssal currents.**Results:**Here is a typical simulation of a cross-equatorial abyssal current. Near the equator y = 0 the Coriolis force cannot balance the force caused by the bottom topography slope in a parabolic channel (height(x,y) ~ x^2) and the current crosses the channel axis x = 0.

**Client:**Forschungszentrum Graz, Austria**Technology:**C++, Matlab**Man hours**: 100**Scope:**The purpose of the project was to develop a code for detailed numerical investigation of the so called double-gradient magnetic instability, which is believed to be responsible for the magnetotail flapping oscillations - the fast vertical oscillations of the Earth’s magnetotail plasma sheet with a quasyperiod about 100−200 s. The set of linearized compressible ideal MHD equations in conservative form is solved by means of Lax-Friedrichs finite difference method. The results of numerical simulations confirm predictions of a simplified analytical theory. In particular, the instability growth rate is found to be close to a peak value provided by analytical estimates.**Results:**Find attached one of the recently published papers (+ see references therein).

Adv. in Space Research 48, 1531 (2011)

**Client:**Forschungszentrum Graz, Austria**Technology:**C++, Matlab**Man hours**: 220**Scope:**Magnetic reconnection is a physical process in highly conducting plasmas in which the magnetic topology is rearranged and magnetic energy is converted to kinetic energy, thermal energy, and particle acceleration. We have developed a numerical code for a remote-sensing method of reconstructing the reconnection rate and the location of X-line from single-spacecraft observations. The method is based on the two-dimensional analytical model of time-dependent Petschek-type magnetic reconnection in compressible plasma with asymmetric magnetic field configuration. The reconstruction technique has been successfully applied to series of flux transfer events recorded by spacecrafts in the Earth's magnetosphere.**Results:**Find attached one of the published papers (+ see references therein).

J. Geophys. Res., 112, A10226 (2007)

**Client:**Private person, USA**Technology:**Matlab**Man hours**: 24**Scope:**The purpose of this very small but also very interesting project was to predict the depth of an object towed through water at the end of a cable. The shape of the cable depends on boat velocity, properties of the object and the cable, temperature and salinity of water and is defined by the balance of tension, gravity, buoyancy, drag and lift forces. The resulting ODE is solved numerically and allows to determine all parameters that are important for the successful trolling.**Results:**Here is the trolling calculator GUI.