Search

The thermal, hydrogeological, and chemical processes affecting Kamchatka geothermal reservoirs were studied by using isotope and geochemistry data: (1) The Geysers Valley hydrothermal reservoirs; (2) The Paratunsky low temperature reservoirs; (3) The North-Koryaksky hydrothermal system; (4) The Mutnovsky high temperature geothermal reservoir; (5) The Pauzhetsky geothermal reservoir. In most cases water isotope in combination with Cl- transient data are found to be useful tool to estimate reservoirs natural and disturbed by exploitation recharge conditions, isotopes of carbon-13 (in CO2) data are pointed either active magmatic recharge took place, while SiO2 and Na-K geothermometers shows opposite time transient trends (Paratunsky, Geysers Valley) suggest that it is necessary to use more complicated geochemical systems of water/mineral equilibria.

The presented approach to discretization of functionally defined heterogeneous objects is oriented towards applications associated with numerical simulation procedures, for example, finite element analysis (FEA). Such applications impose specific constraints upon the resulting surface and volume meshes in terms of their topology and metric characteristics, exactness of the geometry approximation, and conformity with initial attributes. The function representation of the initial object is converted into the resulting cellular representation described by a simplicial complex. We consider in detail all phases of the discretization algorithm from initial surface polygonization to final tetrahedral mesh generation and its adaptation to special FEA needs. The initial object attributes are used at all steps both for controlling geometry and topology of the resulting object and for calculating new attributes for the resulting cellular representation.

International audience; in english. The paper presents the results of the theoretical and numerical performance analysis of the experimental scheme for amplification of magnetic flux intensity via its compression by plasma liner. 0D estimations and 2D computations results are compared. The simulations were carried out with the use of RMHD code MARPLE (IMM RAS). The scheme performance affected by the Rayleigh-Taylor instability, developed in the case of initially disturbed plasma shell density, is studied. The possible penetration of the compressor shell plasma from the discharge chamber into the load area results in the nonuniformity of magnetic pressure in it. The simulation proves the possibility of elimination of this unwanted effect by proper selection of the experiment’s parameters. The correlation of the numerical results for this kind of problems using a simplified 0D model and 2D RMHD simulation with the MARPLE code are demonstrated. The prospects of the plasma magnetic flux compression scheme are discussed. Original Russian Text © V.A. Gasilov, S.V. D’yachenko, A.S. Chuvatin, O.G. Ol’khovskaya, A.S. Boldarev, E.L. Kartasheva, G.A. Bagdasarov, 2009, published in Matematicheskoe Modelirovanie, 2009, Vol. 21, No. 11, pp. 57–73.

The work is devoted to some problems in numerical modelling plasma generated in the course of magnetic implosion of multiwire arrays (liners) in Z-pinch discharge mode, including numerical and theoretical analysis of quasispherical imploding conical wire arrays. Experiments were performed at ANGARA-5-1 facility (TRINITI) with the discharge current up to 3 MA. The numerical simulations were carried out on the base of 2D RMHD code MARPLE (IMM RAS). An appropriate description of the early phase of plasma generation during the electrical explosion of wires is a real challenge in MHD modeling of these experiments, in our computations we used the plasma ablation model, which allows to obtain more realistic results as compared with the models where the entire wires are considered to be initially already in plasma state. According to the experimental data we assume that the wire skeleton is stable during the overall time of the plasma generation as its density is significantly more than that of the surrounding plasma. Besides that the skeleton carries essentially less electric current than the plasma. The plasma ablation process lasts until the total wire skeleton evaporation. Experiments show that in the case of current pulse matched load the plasma generation is finished soon after the current maximum. Taking into account small diameter of the wires as compared with interwire gaps and the array dimensions plasma production is simulated by a steady source surface placed inside the computational domain. The mass ablation rate at the unit of the source surface is calculated with respect to the local distribution of the magnetic field induction. The numerical experiments showed that the plasma ablation model has a significant effect on both the entire scheme of plasma dynamics and such values as voltage drop at the load and soft X-ray yield power. The improvement of the ablation model parameters based on experimental and theoretical estimations is an issue of the day in 2D and 3D Z-pinch simulations.

This paper presents soft X-rays sources of X-pinch type based on the supersonic gas jets. Different variants of the nozzle design are considered. The high-order TVD finite-difference method on unstructured triangular meshes is used in order to solve gas dynamics equations and find the density distribution in the discharge gap. The results of mathematical modelling of the gas jet formation and of the dynamics of plasma originated from the gas jet for the chosen nozzle design and the results of the experimental study for the pinch current 250-300 kA are presented. The numerical simulations of the electrical explosions of the gas density configurations is performed with the help of radiative magnetohydrodynamics codes RAZRYAD (simplified geometry) and MARPLE (more realistic geometry).

MHD numerical simulations are applied for multiparameter studies of quasi-spherical magnetic compression of plasma liners created by conical multiwire arrays electrical explosion at ANGARA-5-1 facility (TRINITI) with the discharge current 2 to 3 MA and the pulse rise time about 100 ns. The description of the plasma dynamics at different stages of implosion is reproduced as a result of simulation. Numerical and experimental time profiles of voltage drop at the load and soft X-ray yield power are compared. The effect of the geometry changes upon the implosion process is studied. Numerical simulation is based on 2D RMHD code MARPLE (IMM RAS) using unstructured triangular grids. The code implements one-fluid two-temperature MHD model, grid-characteristic method for radiative energy transfer and the model of prolonged plasma ablation to simulate plasma source. The governing MHD system of is completed by electrical equation for the full circuit including the generator itself, leading-in systems and the discharge chamber with the plasma in it. Equations of state, transport and kinetic coefficients, opacity and emissivity coefficients are taken from the tables. The behavior of the discharge is satisfactorily described in general by the above RMHD model. The plasma ablation model appeared to have a significant effect on both the entire scheme of plasma dynamics and such values as voltage drop at the load and soft X-ray yield power. The improvement of this model based on experimental and theoretical estimations is an issue of the day in 2D and 3D Z-pinch simulations. The MARPLE code calibrated against the conical liners simulations proved to be a useful tool for computations aimed to optimization of the experimental setup for 3D implosion of plasma.

Work on Z-pinches, that is under way at the SRC RF TRINITI, is aimed at the resolution of two basic problems: to investigate the problem of physics of X-ray generation in a Z-pinch scheme with the goal of its application for ICF at the "Angara-5-1" facility to develop the engineering aspects of multimegampere facilities intended for using them in a Z-pinch scheme for ICF and making the "Baikal" project. At the "Angara-5-1 facility" the measurements of plasma production rates for the wire and foam liners were taken. It has been shown that for Mo, W, Cu, Fe, and Al a current sheath with a radial thickness much larger than the skin- layer forms. The value of the plasma production rate for these metals was obtained. In the case of the liner from agar-agar foam no formation of such a layer was observed. The spectral measurements taken in hybrid Al-W liners make it possible to estimate the electron temperature and concentration from the spectrum of Al ions. It is also shown that an addition of tungsten into the array compound at a constant mass leads to a decrease of not only the electron temperature but also the final radius of stagnation. The measurements of the hard X-radiation appearing at the moment close to a peak of the X-ray pulse are likely to indicate the formation of an several kiloamperes electron beam. The experiments with conical arrays consisting of tungsten wires demonstrate the effects which arise when the wire array is imploded both in radial and axial directions. For large cone angles there appear two separate radiating regions. The corresponding X-ray pulses have a double-peak structure. Calculations of such liner compression are used to test the programs describing both radial and axial plasma compression. The results obtained at the IMM RAS fairly well provide a quantitative description of the X-ray pulse structure. Such conical liners prove to be a promising tool to test complicated RMHD codes. A pulse generation scheme relevant to the "Baikal" project is under development at the MOL facility. A 30-fold reliable current multiplication in the inductive storage has been obtained. The storage energy was effectively transferred to the subsequent stage of formation, i.e. a magnetic amplifier. The current in the latter was multiplied by a factor of 1.55. At the "Angara-5-1" facility the MEVL code reporting the process of magnetic insulation arising in the in-parallel- connected magnetically insulated vacuum lines has been tested.

The automation of the function-based (FRep) volumetric modeling task is tackled by introducing template parameterized models and a procedure for recovery of constructive models from segmented point-sets. In order to reuse existing models, we propose to parameterize them and to fit the parameters to different point-sets for optimizing and adapting the shape to different objects of the same class of shapes. The automation of the creation of a constructive FRep model is also considered by creating a recovery procedure for a given segmented pointset and a list of corresponding primitives. A genetic algorithm is used to find the best constructive expression for the object with the given set of primitives in the point cloud segmentation and the set of available operations. The proposed approach is illustrated by fitting of different models to point clouds and by the automatic generation of constructive trees from segmented point-sets for real mechanical parts.

In this paper we further develop a novel approach for modelling heterogeneous objects containing entities of various dimensions and representations within a cellular-functional framework based on the implicit complex notion. We provide a brief description for implicit complexes and describe their structure including both the geometry and topology of cells of different types. Then the paper focuses on the development of algorithms for set-theoretic operations on heterogeneous objects represented by implicit complexes. We also describe a step-by-step procedure for the construction of a hybrid model using these operations. Finally, we present a case-study showing how to construct a hybrid model integrating both boundary and function representations. Our examples also illustrate modelling with attributes and dynamic modelling.

This paper describes a technology for modelling and rendering heterogeneous objects containing entities of various dimensionalities within a cellular-functional framework based on the implicit complex notion. Implicit complexes make it possible to combine a cellular representation and a constructive function representation. We describe a formal framework for such a hybrid representation and propose a general structure for implicit complexes. Then, we consider how an implicit complex can be described geometrically and topologically along with its associated attributes. Rendering algorithms for implicit complexes using ray-tracing are also discussed. Finally, we present a case study illustrating the proposed methods and algorithms.

The presented approach to discretization of functionally defined heterogeneous objects is oriented towards applications associated with numerical simulation procedures, for example, finite element analysis (FEA). Such applications impose specific constraints upon the resulting surface and volume meshes in terms of their topology and metric characteristics, exactness of the geometry approximation, and conformity with initial attributes. The function representation of the initial object is converted into the resulting cellular representation described by a simplicial complex. We consider in detail all phases of the discretization algorithm from initial surface polygonization to final tetrahedral mesh generation and its adaptation to special FEA needs. The initial object attributes are used at all steps both for controlling geometry and topology of the resulting object and for calculating new attributes for the resulting cellular representation.

The paper presents an approach to modeling heterogeneous objects as multidimensional point sets with multiple attributes (hypervolumes). A theoretical framework is based on a hybrid model of hypervolumes combining a cellular representation and a constructive representation using real-valued functions. This model allows for independent but unifying representation of geometry and attributes, and makes it possible to represent dimensionally non-homogeneous entities and their cellular decompositions. Hypervolume model components such as objects, operations and relations are introduced and outlined. The framework's inherent multidimensionality allowing, in particular, to deal naturally with time dependence promises to model complex dynamic objects composed of different materials with constructive building of their geometry and attributes. Attributes given at each point can represent properties of arbitrary nature (material, photometric, physical, statistical, etc.). To demonstrate a particular application of the proposed framework, we present an example of multimaterial modeling - a multilayer geological structure with cavities and wells. Another example illustrating the treatment of attributes other than material distributions is concerned with time-dependent adaptive mesh generation where function representation is used to describe object geometry and density of elements in the cellular model of the mesh. The examples have been implemented by using a specialized modeling language and software tools being developed by the authors.