Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

Flexible carbon-based ohmic contacts for organic transistors

The present invention relates to a system and method of organic thin-film transistors (OTFTs). More specifically, the present invention relates to employing a flexible, conductive particle-polymer composite material for ohmic contacts (i.e. drain and source).

Brandon, Erik↗

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗

Electride Mediated Surface Enhanced Raman Scattering (SERS)

An electride may provide surface enhanced Raman scattering (SERS). The electride, a compound where the electrons serve as anions, may be a ceramic electride, such as a conductive ceramic derived from mayenite, or an organic electride, for example. The textured electride surface or electride particles may strongly enhance the Raman scattering of organic or other Raman active analytes. This may also provide a sensitive method for monitoring the chemistry and electronic environment at the electride surface. The results are evidence of a new class of polariton (i.e., a surface electride-polariton resonance mechanism) that is analogous to the surface plasmon-polariton resonance that mediates conventional SERS.

Anderson, Mark S.↗

Radiation, Smoke and Clouds Observed in the Southeastern Atlantic With the Research Scanning Polarimeter During ORACLES

The ObseRvations of Aerosols above Clouds and their interactions (ORACLES) project is making a series of field deployments to the southeastern Atlantic with NASA ER-2 and P3 aircraft to acquire both detailed remote sensing observations and in situ measurements of the aerosols and clouds in that region. This area is home to one of the largest low-level cloud decks on Earth that is seasonally affected by vast plumes of smoke from biomass burning, which in effect provides a natural experiment testing the radiative and microphysical interactions between the smoke and the clouds. The downward solar radiation at the surface, or cloud top, is always reduced by the presence of smoke. However, whether the amount of sunlight reflected back out to space is increased, or decreased by the presence of smoke is sensitively dependent on the brightness of the clouds and the fraction of light that the smoke absorbs each time light hits a smoke particle. In this study we use data from the Research Scanning Polarimeter, an along track scanning instrument, that provides measurements of the Stokes parameters I, Q and U at 410, 470, 555, 670, 865, 960, 1590, 1880 and 2260 nm at 150 viewing angles over a range of +/- 60° from nadir for each contiguous sub-aircraft pixel (~ 300 m in size). A retrieval algorithm is applied to the data acquired with a table look up technique, similar to that of the operational POLDER algorithm, to provide a first guess of the complex refractive index, optical depth and size distribution of the smoke particles together with cloud droplet size and optical depth. A subsequent iterative fitting procedure, where the fact that the doubling/adding method allows the construction of the Green's function for the radiative transfer equation, is used to obtain an efficient and statistically optimal estimate of the aerosol and cloud retrieval parameters. These retrieval parameters are evaluated against in situ observations, when available, and the optical depth and intensive lidar variables that are measured by the High Spectral Resolution Lidar 2. Finally, the aerosol and cloud retrievals are used to evaluate the variations in top of the atmosphere, surface/cloud top shortwave radiative forcing and atmospheric absorption that are caused by variations in the smoke and clouds.

Cairns, Brian↗

Retrieving Aerosol Characteristics From the PACE Mission, Part 1: Ocean Color Instrument

NASA’s Plankton, Aerosol, Clouds, ocean Ecosystem (PACE) satellite mission is scheduled to launch in 2022, with the Ocean Color Instrument (OCI) on board. For the first time reflected sunlight from the Earth across a broad spectrum from the ultraviolet (UV: 350 nm) to the short wave infrared (SWIR: 2260 nm) will be measured from a single instrument at 1 km spatial resolution. While seven discrete bands will represent the SWIR, the spectrum from 350 to 890 nm will be continuously covered with a spectral resolution of 5 nm. OCI will thus combine in a single instrument (and at an enhanced spatial resolution for the UV) the heritage capabilities of the Moderate resolution Imaging Spectroradiometer (MODIS) and the Ozone Monitoring Instrument (OMI), while covering the oxygen A-band (O2A). Designed for ocean color and ocean biology retrievals, OCI also enables continuation of heritage satellite aerosol products and the development of new aerosol characterization from space. In particular the combination of MODIS and OMI characteristics allows deriving aerosol height, absorption and optical depth along with a measure of particle size distribution. This is achieved by using the traditional MODIS visible-to-SWIR wavelengths to constrain spectral aerosol optical depth and particle size. Extrapolating this information to the UV channels allows retrieval of aerosol absorption and layer height. A more direct method to derive aerosol layer height makes use of O2A absorption methods, despite the relative coarseness of the nominal 5 nm spectral resolution of OCI. Altogether the PACE mission with OCI will be an unprecedented opportunity for aerosol characterization that will continue climate data records from the past decades and propel aerosol science forward toward new opportunities.

PACE↗

Nano-Engineered Catalysts for Direct Methanol Fuel Cells

Nano-engineered catalysts, and a method of fabricating them, have been developed in a continuing effort to improve the performances of direct methanol fuel cells as candidate power sources to supplant primary and secondary batteries in a variety of portable electronic products. In order to realize the potential for high energy densities (as much as 1.5 W h/g) of direct methanol fuel cells, it will be necessary to optimize the chemical compositions and geometric configurations of catalyst layers and electrode structures. High performance can be achieved when catalyst particles and electrode structures have the necessary small feature sizes (typically of the order of nanometers), large surface areas, optimal metal compositions, high porosity, and hydrophobicity. The present method involves electrodeposition of one or more catalytic metal(s) or a catalytic-metal/polytetrafluoroethylene nanocomposite on an alumina nanotemplate. The alumina nanotemplate is then dissolved, leaving the desired metal or metal/polytetrafluoroethylene-composite catalyst layer. Unlike some prior methods of making fine metal catalysts, this method does not involve processing at elevated temperature; all processing can be done at room temperature. In addition, this method involves fewer steps and is more amenable to scaling up for mass production. Alumina nanotemplates are porous alumina membranes that have been fabricated, variously, by anodizing either pure aluminum or aluminum that has been deposited on silicon by electronbeam evaporation. The diameters of the pores (7 to 300 nm), areal densities of pores (as much as 7 x 10(exp 10)sq cm), and lengths of pores (up to about 100 nm) can be tailored by selection of fabrication conditions. In a given case, the catalytic metal, catalytic metal alloy, or catalytic metal/ polytetrafluoroethylene composite is electrodeposited in the pores of the alumina nanotemplate. The dimensions of the pores, together with the electrodeposition conditions, determine the sizes and surface areas of the catalytic particles. Hence, the small features and large surface areas of the porosity translate to the desired small particle size and large surface area of the catalyst (see figure). When polytetrafluoroethylene is included, it is for the purpose of imparting hydrophobicity in order to prevent water from impeding the desired diffusion of gases through the catalyst layer. To incorporate polytetrafluoroethylene into a catalytic-metal/polytetrafluoroethylene nanocomposite, one suspends polytetrafluoroethylene nanoparticles in the electrodeposition solution. The polytetrafluoroethylene content can be varied to obtain the desired degree of hydrophobicity and permeability by gas.

Myung, Nosang↗

Forces Associated with Nonlinear Nonholonomic Constraint Equations

A concise method has been formulated for identifying a set of forces needed to constrain the behavior of a mechanical system, modeled as a set of particles and rigid bodies, when it is subject to motion constraints described by nonholonomic equations that are inherently nonlinear in velocity. An expression in vector form is obtained for each force; a direction is determined, together with the point of application. This result is a consequence of expressing constraint equations in terms of dot products of vectors rather than in the usual way, which is entirely in terms of scalars and matrices. The constraint forces in vector form are used together with two new analytical approaches for deriving equations governing motion of a system subject to such constraints. If constraint forces are of interest they can be brought into evidence in explicit dynamical equations by employing the well-known nonholonomic partial velocities associated with Kane's method; if they are not of interest, equations can be formed instead with the aid of vectors introduced here as nonholonomic partial accelerations. When the analyst requires only the latter, smaller set of equations, they can be formed directly; it is not necessary to expend the labor to form the former, larger set first and subsequently perform matrix multiplications.

Roithmayr, Carlos M.↗

Microgravity and Hypogravity Compatible Methods for the Destruction of Solid Wastes by Magnetically Assisted Gasification

This report summarizes a three-year collaborative effort between researchers at UMPQUA Research Company (URC) and the Chemical Engineering Department at Oregon State University (OSU). The Magnetically Assisted Gasification (MAG) concept was originally conceived as a microgravity and hypogravity compatible means for the decomposition of solid waste materials generated aboard spacecraft, lunar and planetary habitations, and for the recovery of potentially valuable resources. While a number of methods such as supercritical water oxidation (SCW0), fluidized bed incineration, pyrolysis , composting and related biological processes have been demonstrated for the decomposition of solid wastes, none of these methods are particularly well- suited for employment under microgravity or hypogravity conditions. For example, fluidized bed incineration relies upon a balance between drag forces which the flowing gas stream exerts upon the fluidization particles and the opposing force of gravity. In the absence of gravity, conventional fluidization cannot take place. Hypogravity operation can also be problematic for conventional fluidized bed reactors, because the various factors which govern fluidization phenomena do not all scale linearly with gravity. For this reason it may be difficult to design and test fluidized bed reactors in lg, which are intended to operate under different gravitational conditions. However, fluidization can be achieved in microgravity (and hypogravity) if a suitable replacement force to counteract the forces between fluid and particles can be found. Possible alternatives include: centripetal force, electric fields, or magnetic fields. Of these, magnetic forces created by the action of magnetic fields and magnetic field gradients upon ferromagnetic media offer the most practical approach. The goal of this URC-OSU collaborative effort was to develop magnetic hardware and methods to control the degree of fluidization (or conversely consolidation) of granular ferromagnetic media and to employ these innovations in sequential filtration and fluidized bed processes for the segregation and decomposition of solid waste materials, and for the concentration and collection of inorganic residue (ash). This required the development of numerous enabling technologies and tools.

Atwater, James E.↗

Overview of Nondestructive Evaluation Capabilities at NASA MSFC

Nondestructive evaluation (NDE) methods examine material integrity without impairing its usefulness. The NDE team at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama is responsible for applying existing methods to new hardware designs and material systems, probability of detection studies to quantify detection capabilities, overseeing NDE requirements for flight hardware, inspecting development articles and flight hardware, engineering support for failure investigations, and investigating emerging NDE methods. NASA MSFC has been on the cutting edge of developing new NDE techniques, first for the space shuttle, then for composite structures, and most recently for additive manufacturing (AM). Standard in-house inspection capabilities include eddy current, magnetic particle, liquid penetrant, ultrasonic, and radiographic testing (including computed tomography). Other advanced inspection capabilities include infrared flash thermography, shearography, acoustic emission, and microwave/millimeter wave testing. Some recent research studies include AM probability of detection testing, investigating in-situ process monitoring methods for AM, and correlating the NDE response of true fatigue cracks with artificial notches.

nondestructive evaluation↗

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗

Plasma Diffusion in Self-Consistent Fluctuations

The problem of particle diffusion in position space, as a consequence ofeleclromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resuiting from an agyrotropic in itiai setting)is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (result ing from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probabi lity distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

Smets, R.↗

Bench-scale synthesis of nanoscale materials

A novel flow-through hydrothermal method used to synthesize nanoscale powders is introduced by Pacific Northwest Laboratory. The process, Rapid Thermal Decomposition of precursors in Solution (RTDS), uniquely combines high-pressure and high-temperature conditions to rapidly form nanoscale particles. The RTDS process was initially demonstrated on a laboratory scale and was subsequently scaled up to accommodate production rates attractive to industry. The process is able to produce a wide variety of metal oxides and oxyhydroxides. The powders are characterized by scanning and transmission electron microscopic methods, surface-area measurements, and x-ray diffraction. Typical crystallite sizes are less than 20 nanometers, with BET surface areas ranging from 100 to 400 sq m/g. A description of the RTDS process is presented along with powder characterization results. In addition, data on the sintering of nanoscale ZrO2 produced by RTDS are included.

Buehler, M. F.↗

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗

Electromagnetic Dissociation Cross Sections using Weisskopf-Ewing Theory

It is important that accurate estimates of crew exposure to radiation are obtained for future long-term space missions. Presently, several space radiation transport codes exist to predict the radiation environment, all of which take as input particle interaction cross sections that describe the nuclear interactions between the particles and the shielding material. The space radiation transport code HZETRN uses the nuclear fragmentation model NUCFRG2 to calculate Electromagnetic Dissociation (EMD) cross sections. Currently, NUCFRG2 employs energy independent branching ratios to calculate these cross sections. Using Weisskopf-Ewing (WE) theory to calculate branching ratios, however, is more advantageous than the method currently employed in NUCFRG2. The WE theory can calculate not only neutron and proton emission, as in the energy independent branching ratio formalism used in NUCFRG2, but also deuteron, triton, helion, and alpha particle emission. These particles can contribute significantly to total exposure estimates. In this work, photonuclear cross sections are calculated using WE theory and the energy independent branching ratios used in NUCFRG2 and then compared to experimental data. It is found that the WE theory gives comparable, but mainly better agreement with data than the energy independent branching ratio. Furthermore, EMD cross sections for single neutron, proton, and alpha particle removal are calculated using WE theory and an energy independent branching ratio used in NUCFRG2 and compared to experimental data.

Adamczyk, Anne M.↗

Predicting the Mineral Composition of Dust Aerosols: Representing Key Processes - Part 1

Soil dust aerosols created by wind erosion are typically assigned globally uniform physical and chemical properties within Earth system models, despite known regional variations in the mineral content of the parent soil. Mineral composition of the aerosol particles is important to their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Here, aerosol mineral composition is derived by extending a method that provides the composition of a wet-sieved soil. The extension accounts for measurements showing significant differences between the mineral fractions of the wetsieved soil and the emitted aerosol concentration. For example, some phyllosilicate aerosols are more prevalent at silt sizes, even though they are nearly absent at these diameters in a soil whose aggregates are dispersed by wet sieving. We calculate the emitted mass of each mineral with respect to size by accounting for the disintegration of soil aggregates during wet sieving. These aggregates are emitted during mobilization and fragmentation of the original undispersed soil that is subject to wind erosion. The emitted aggregates are carried far downwind from their parent soil. The soil mineral fractions used to calculate the aggregates also include larger particles that are suspended only in the vicinity of the source. We calculate the emitted size distribution of these particles using a normalized distribution derived from aerosol measurements. In addition, a method is proposed for mixing minerals with small impurities composed of iron oxides. These mixtures are important for transporting iron far from the dust source, because pure iron oxides are more dense and vulnerable to gravitational removal than most minerals comprising dust aerosols. A limited comparison to measurements from North Africa shows that the model extensions result in better agreement, consistent with a more extensive comparison to global observations as well as measurements of elemental composition downwind of the Sahara, as described in companion articles.

aerosols↗

Examination of Multiple Lithologies Within the Primitive Ordinary Chondrite NWA 5717

Northwest Africa 5717 is a primitive (subtype 3.05) ungrouped ordinary chondrite which contains two apparently distinct lithologies. In large cut meteorite slabs, the darker of these, lithology A, looks to host the second, much lighter in color, lithology B (upper left, Fig. 1). The nature of the boundary between the two is uncertain, ranging from abrupt to gradational and not always following particle boundaries. The distinction between the lithologies, beyond the obvious color differences, has been supported by a discrepancy in oxygen isotopes and an incongruity in the magnesium contents of chondrule olivine. Here, quantitative textural analysis and mineralogical methods have been used to investigate the two apparent lithologies within NWA 5717. Olivine grains contained in a thin section from NWA 7402, thought to be paired to 5717, were also measured to re-examine the distinct compositional range among the light and dark areas. Procedure: Particles from a high-resolution mosaic image of a roughly 13x15cm slice of NWA 5717 were traced in Adobe Photoshop. Due to the large size of the sample, visually representative regions of each lithology were chosen to be analyzed. The resulting layers of digitized particles were imported into ImageJ, which was used to measure their area, along with the axes, the angle from horizontal, and the centroid coordinates of ellipses fitted to each particle following the approach. Resulting 2D pixel areas were converted to spherical diameters employing the unfolding algorithm, which outputs a 3D particle size distribution based on digitized 2D size frequency data. Spatstat was used to create kernel density plots of the centroid coordinates for each region. X-ray compositional maps, microprobe analyses, and Mossbauer spectroscopy was conducted on a thin section of NWA 7402, tentatively paired to NWA 5717.

Cato, M. J.↗

Isentropic phase changes in dissociating gases and the method of sound dispersion for the investigation of homogeneous gas reactions with very high speed : Conclusion

The analytical results of Part I are also applied to sound dispersion by friction and heat conduction, An irreversible change of momentum, energy, and type of particle corresponding to friction, heat conduction, and diffusion effects can appear both in the direction of the sound field and traverse to it. Longitudinal damping, the coupling of longitudinal damping and that due to chemical and physical changes, and coupling of diffusion and compositional changes are treated for a plane sound wave of infinite extent. The same principles are also applied to sound effects in cylindrical tubes. The limitations of the method are discussed in some detail.

SOUND VELOCITY - GASES↗

Charge Resolution of the Silicon Matrix of the ATIC Experiment

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CRT particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. The silicon matrix of the ATIC spectrometer is designed to resolve individual elements from proton to iron. To provide this resolution careful calibration of each pixel of the silicon matrix is required. Firstly, for each electronic channel of the matrix the pedestal value was subtracted taking into account its drift during the flight. The muon calibration made before the flight was used then to convert electric signals (in ADC channel number) to energy deposits in each pixel. However, the preflight muon calibration was not accurate enough for the purpose, because of lack of statistics in each pixel. To improve charge resolution the correction was done for the position of Helium peak in each pixel during the flight . The other way to set electric signals in electronics channels of the Si-matrix to one scale was correction for electric channel gains accurately measured in laboratory. In these measurements it was found that small different nonlinearities for different channels are present in the region of charge Z > 20. The correction for these non-linearities was not done yet. In linear approximation the method provides practically the same resolution as muon calibration plus He-peak correction. For searching a pixel with the signal of primary particle an indication from the cascade in the calorimeter was used. For this purpose a trajectory was reconstructed using weight centers of energy deposits in BGO layers. The point of intersection of this trajectory with Si-matrix and its RMS was determined. The pixel with maximal signal in 3sigma region was taken as sought. The signal in this pixel was corrected by trajectory zenith angle. The preliminary results on charge resolution of the Si-matrix in the range from protons to iron are presented.

Zatsepin, V. I.↗