Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Zeolite membrane cascade for tritium extraction and recovery systems

Membrane separation by zeolite membranes has been proposed as a pre-concentration stage for the tritium extraction from the purge helium of the breeding blanket combined with a final recovery by the catalytic membrane reactor PERMCAT. This fully continuous operation improves the tritium management in fusion machines, minimizing the tritium inventory. For the first time, the permeation measurements for H{sub 2} - He mixtures through a MFI-alumina hollow fibre membrane has been measured for different compositions (0.1 - 20% H{sub 2}) and temperatures. Such a highly permeable membrane, although it shows a limited selectivity, appears attractive for tritium recovery in the blanket. This will imply its operation in a membrane cascade, for which simulation work is ongoing. Mathematically the process is modeled using mass balance equations that can be transformed into the matrix form and solved iteratively assuming a permeate concentration on the first step of iteration, until the separation requirements are fulfilled.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Dynamic tomography with a priori information

We present ''dynamic tomography'' algorithms that allow for the high-resolution, time-resolved imaging of dynamic (i.e., continuously time evolving) complex systems at existing x-ray micro-CT facilities. The behavior of complex systems is constrained by the underlying physics. By exploiting a priori knowledge of the geometry of the physical process being studied to allow the use of sophisticated iterative reconstruction techniques that incorporate constraints, we improve on current frame rates by at least an order of magnitude. This allows time-resolved imaging of previously intractable processes, such as two-phase fluid flow. We present reconstructions from experimental data collected at the Australian National University x-ray micro-CT facility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Information content of and the ability to reconstruct dichroic X-ray tomography and laminography

Dichroic tomography is a 3D imaging technique in which the polarization of the incident beam is used to induce contrast due to the magnetization or orientation of a sample. The aim is to reconstruct not only the optical density but the dichroism of the sample. The theory of dichroic tomographic and laminographic imaging in the parallel-beam case is discussed as well as the problem of reconstruction of the sample’s optical properties. The set of projections resulting from a single tomographic/laminographic measurement is not sufficient to reconstruct the magnetic moment for magnetic circular dichroism unless additional constraints are applied or data are taken at two or more tilt angles. For linear dichroism, three polarizations at a common tilt angle are insufficient for unconstrained reconstruction. However, if one of the measurements is done at a different tilt angle than the other, or the measurements are done at a common polarization but at three distinct tilt angles, then there is enough information to reconstruct without constraints. Possible means of applying constraints are discussed. Furthermore, it is shown that for linear dichroism, the basic assumption that the absorption through a ray path is the integral of the absorption coefficient, defined on the volume of the sample, along the ray path, is not correct when dichroism or birefringence is strong. This assumption is fundamental to tomographic methods. An iterative algorithm for reconstruction of linear dichroism is demonstrated on simulated data.

Marcus, Matthew A. (ORCID:0000000325277586)↗
DOE OSTI2015

Argonne Bubble Experiment Thermal Model Development

This report will describe the Computational Fluid Dynamics (CFD) model that was developed to calculate the temperatures and gas volume fractions in the solution vessel during the irradiation. It is based on the model used to calculate temperatures and volume fractions in an annular vessel containing an aqueous solution of uranium . The experiment was repeated at several electron beam power levels, but the CFD analysis was performed only for the 12 kW irradiation, because this experiment came the closest to reaching a steady-state condition. The aim of the study is to compare results of the calculation with experimental measurements to determine the validity of the CFD model.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2015

Verification and Validation of Digitally Upgraded Control Rooms

As nuclear power plants undertake main control room modernization, a challenge is the lack of a clearly defined human factors process to follow. Verification and validation (V&V) as applied in the nuclear power community has tended to involve efforts such as integrated system validation, which comes at the tail end of the design stage. To fill in guidance gaps and create a step-by-step process for control room modernization, we have developed the Guideline for Operational Nuclear Usability and Knowledge Elicitation (GONUKE). This approach builds on best practices in the software industry, which prescribe an iterative user-centered approach featuring multiple cycles of design and evaluation. Nuclear regulatory guidance for control room design emphasizes summative evaluation—which occurs after the design is complete. In the GONUKE approach, evaluation is also performed at the formative stage of design—early in the design cycle using mockups and prototypes for evaluation. The evaluation may involve expert review (e.g., software heuristic evaluation at the formative stage and design verification against human factors standards like NUREG-0700 at the summative stage). The evaluation may also involve user testing (e.g., usability testing at the formative stage and integrated system validation at the summative stage). An additional, often overlooked component of evaluation is knowledge elicitation, which captures operator insights into the system. In this report we outline these evaluation types across design phases that support the overall modernization process. The objective is to provide industry-suitable guidance for steps to be taken in support of the design and evaluation of a new human-machine interface (HMI) in the control room. We suggest the value of early-stage V&V and highlight how this early-stage V&V can help improve the design process for control room modernization. We argue that there is a need to overcome two shortcomings of V&V in current practice—the propensity for late-stage V&V and the use of increasingly complex psychological assessment measures for V&V.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

MOC Efficiency Improvements Using a Jacobi Inscatter Approximation

In recent weeks, attention has been given to resolving the convergence issues encountered with TCP 0 by trying a Jacobi (J) inscatter approach when group sweeping, where the inscatter source is constructed using the previous iteration flux. This is in contrast to a Gauss-Seidel (GS) approach, which has been the default to-date, where the scattering source uses the most up-to-date flux values. The former is consistent with CASMO, which has no issues with TCP 0 convergence. Testing this out on a variety of problems has demonstrated that the Jacobi approach does indeed provide substantially more stability, though can take more outer iterations to converge. While this is not surprising, there are improvements that can be made to the MOC sweeper to capitalize on the Jacobi approximation and provide substantial speedup. For example, the loop over groups, which has traditionally been the outermost loop in MPACT, can be moved to the interior, avoiding duplicate modular ray trace and coarse ray trace setup (mapping coarse mesh surface indexes), which needs to be performed repeatedly when group is outermost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Progress Update on Iterative Reconstruction of Neutron Tomographic Images

This report satisfies the fiscal year 2016 technical deliverable to report on progress in development of fast iterative reconstruction algorithms for project OR16-3DTomography-PD2Jb, "3D Tomography and Image Processing Using Fast Neutrons." This project has two overall goals. The first of these goals is to extend associated-particle fast neutron transmission and, particularly, induced-reaction tomographic imaging algorithms to three dimensions. The second of these goals is to automatically segment the resultant tomographic images into constituent parts, and then extract information about the parts, such as the class of shape and potentially shape parameters. This report addresses of the component of the project concerned with three-dimensional (3D) image reconstruction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

High-Fidelity Calibration and Characterization of Hyperspectral Computed Tomography System

This document archives the results developed by the Lab Directed Research and Development (LDRD) project sponsored by Sandia National Laboratories (SNL). In this work, a numerical study was performed to show the feasibility of approximating the non-linear operator of SNL's unique high-energy hyperspectral computed tomography system as a sequence of linear operators. The four main results gained from this work include the development of a simulation test-bed using a particle-transport Monte Carlo approach; the demonstration to assemble a linear operator of almost-arbitrary resolution for a given narrow energy window, developing a compression approach to dramatically reduce the size of the linear operator via a spline approach, and the demonstration of using the linear operator to perform processing of x-ray data; in this case, the development of an iterative reconstruction method. This numerical study has indicated that if these results can be replicated on the SNL system, the improved performance could be revolutionary as the method to approximate the nonlinear operator for a hyperspectral CT system is not feasible to perform on a traditional CT system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Modeling and Simulation of Xe-100-type Pebble Bed Gas-Cooled Reactor with SCALE

The US Department of Energy (DOE) announced the Advanced Reactor Demonstration Program (ARDP) to accelerate the deployment of advanced reactor concepts. Awardees of ARDP funds are expected to demonstrate the operation of an advanced reactor within 7 years of receiving the award. X-Energy’s advanced reactor concept, the Xe-100, was selected as one of two advanced reactor concepts to receive funding to demonstrate the operation of its high-temperature gas-cooled pebble-bed reactor before the end of this decade. As a result of this push to bring advanced reactors to maturation and commercialization, transition and deployment scenario studies are being performed under the Systems Analysis and Integration (SA&I) campaign within the DOE Office of Nuclear Energy (DOE-NE) to evaluate the transition of the current US commercial fleet of light-water reactors (LWRs) to a future fleet of advanced reactors consisting of a mix of ARDP type reactor concepts and advanced LWRs. To accurately evaluate the front- and back-end resource requirements, it is important to perform reactor physics calculations to determine the discharge burnup and isotopic content, fuel residence time, as well as other parameters. For this purpose, a summer project funded by the SA&I campaign allowed for the setup of SCALE models for full-core Xe-100 type high-temperature gas-cooled pebble-bed reactor and a Xe-100 type slice using publicly available information. The core-averaged equilibrium compositions and zone-wise equilibrium compositions for the slice and 3D models, respectively, were obtained following an iterative depletion method developed by Bostelmann et al. using SCALE’s reactor physics sequence TRITON. The slice model was used with TRITON to generate burnup-dependent cross section libraries at different temperatures which can be used with SCALE’s ORIGAMI code to rapidly determine fuel inventory and therefore to perform quick sensitivity studies on parameters such as the pebble location in the core. The SCALE/TRITON transport and depletion calculation for the Xe-100 type slice model indicates that the isotopic concentrations are in good agreement at 1,300 effective full power days (EFPD) for 235 U. An analysis of 236 U results match 239 Pu results would seem to indicate a typographical error in Mulder and Boyes wherein the reported results of 236 U and 239 Pu are reversed. In addition to SCALE/TRITON calculations, a new capability within SCALE/ORIGAMI for the simulation of pebble-bed reactors was used to study the burnup sensitivity with respect to the pebble pathway through the core. The SCALE/ORIGAMI results show that pebbles that travel closer to the reflector for the entire depletion history have a higher burnup than pebbles that travel through the middle of the core because of the higher thermal to fast flux ratio near the reflector. Consequently, a pebble’s burnup is strongly affected by the pebble’s pathway for each pass. Additional phenomena such as temperature distributions in the core and different travel times of the pebbles in the individual radial zones further affect the burnup distribution. The sensitivity of the discharge vector to the pebble pathways taken during each pass can be evaluated in the future using SCALE/ORIGAMI now that the SCALE inputs have been established.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

An Early Investigation of the HHL Quantum Linear Solver for Scientific Applications

In this paper, we explore using the Harrow–Hassidim–Lloyd (HHL) algorithm to address scientific and engineering problems through quantum computing, utilizing the NWQSim simulation package on a high-performance computing platform. Focusing on domains such as power-grid management and climate projection, we demonstrate the correlations of the accuracy of quantum phase estimation, along with various properties of coefficient matrices, on the final solution and quantum resource cost in iterative and non-iterative numerical methods such as the Newton–Raphson method and finite difference method, as well as their impacts on quantum error correction costs using the Microsoft Azure Quantum resource estimator. We summarize the exponential resource cost from quantum phase estimation before and after quantum error correction and illustrate a potential way to reduce the demands on physical qubits. This work lays down a preliminary step for future investigations, urging a closer examination of quantum algorithms’ scalability and efficiency in domain applications.

hybrid software for QC-HPC↗
DOE OSTI2016

CHEMICAL COMPOSITION OF INTERMEDIATE-MASS STAR MEMBERS OF THE M6 (NGC 6405) OPEN CLUSTER

We present here the first abundance analysis of 44 late B-, A-, and F-type members of the young open cluster M6 (NGC 6405, age about 75 Myr). Low- and medium-resolution spectra, covering the 4500–5840 Å wavelength range, were obtained using the FLAMES/GIRAFFE spectrograph attached to the ESO Very Large Telescopes. We determined the atmospheric parameters using calibrations of the Geneva photometry and by adjusting the H{sub β} profiles to synthetic ones. The abundances of up to 20 chemical elements, from helium to mercury, were derived for 19 late B, 16 A, and 9 F stars by iteratively adjusting synthetic spectra to the observations. We also derived a mean cluster metallicity of [Fe/H] = 0.07 ± 0.03 dex from the iron abundances of the F-type stars. We find that for most chemical elements, the normal late B- and A-type stars exhibit larger star-to-star abundance variations than the F-type stars probably because of the faster rotation of the B and A stars. The abundances of C, O, Mg, Si, and Sc appear to be anticorrelated with that of Fe, while the opposite holds for the abundances of Ca, Ti, Cr, Mn, Ni, Y, and Ba as expected if radiative diffusion is efficient in the envelopes of these stars. In the course of this analysis, we discovered five new peculiar stars: one mild Am, one Am, and one Fm star (HD 318091, CD-32 13109, GSC 07380-01211, CP1), one HgMn star (HD 318126, CP3), and one He-weak P-rich (HD 318101, CP4) star. We also discovered a new spectroscopic binary, most likely a SB2. We performed a detailed modeling of HD 318101, the new He-weak P-rich CP star, using the Montréal stellar evolution code XEVOL which self-consistently treats all particle transport processes. Although the overall abundance pattern of this star is properly reproduced, we find that detailed abundances (in particular the high P excess) resisted modeling attempts even when a range of turbulence profiles and mass-loss rates were considered. Solutions are proposed which are still under investigation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Spectroscopic measurements of the far-ultraviolet dust attenuation curve at z ∼ 3

We present the first spectroscopic measurements of the shape of the far-ultraviolet (far-UV; λ=950--1500 Å) dust attenuation curve at high redshift (z∼3). Our analysis employs rest-frame UV spectra of 933 galaxies at z∼3, 121 of which have very deep spectroscopic observations (≳7 hr) at λ=850--1300 A-ring , with the Low Resolution Imaging Spectrograph on the Keck Telescope. By using an iterative approach in which we calculate the ratios of composite spectra in different bins of continuum color excess, E(B−V), we derive a dust curve that implies a lower attenuation in the far-UV for a given E(B−V) than those obtained with standard attenuation curves. We demonstrate that the UV composite spectra of z∼3 galaxies can be modeled well by assuming our new attenuation curve, a high covering fraction of H i, and absorption from the Lyman–Werner bands of H{sub 2} with a small (≲20%) covering fraction. The low covering fraction of H{sub 2} relative to that of the H I and dust suggests that most of the dust in the ISM of typical galaxies at z∼3 is unrelated to the catalysis of H{sub 2}, and is associated with other phases of the ISM (i.e., the ionized and neutral gas). The far-UV dust curve implies a factor of ≈2 lower dust attenuation of Lyman continuum (ionizing) photons relative to those inferred from the most commonly assumed attenuation curves for L* galaxies at z∼3. Our results may be utilized to assess the degree to which ionizing photons are attenuated in H ii regions or, more generally, in the ionized or low column density (N(H I)≲10{sup 17.2} cm{sup −2}) neutral ISM of high-redshift galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE KEPLER CATALOG OF STELLAR FLARES

A homogeneous search for stellar flares has been performed using every available Kepler light curve. An iterative light curve de-trending approach was used to filter out both astrophysical and systematic variability to detect flares. The flare recovery completeness has also been computed throughout each light curve using artificial flare injection tests, and the tools for this work have been made publicly available. The final sample contains 851,168 candidate flare events recovered above the 68% completeness threshold, which were detected from 4041 stars, or 1.9% of the stars in the Kepler database. The average flare energy detected is ∼10{sup 35} erg. The net fraction of flare stars increases with g − i color, or decreasing stellar mass. For stars in this sample with previously measured rotation periods, the total relative flare luminosity is compared to the Rossby number. A tentative detection of flare activity saturation for low-mass stars with rapid rotation below a Rossby number of ∼0.03 is found. A power-law decay in flare activity with Rossby number is found with a slope of −1, shallower than typical measurements for X-ray activity decay with Rossby number.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Characterizing and Improving the Data Reduction Pipeline for the Keck OSIRIS Integral Field Spectrograph

OSIRIS is a near-infrared (1.0–2.4 μm) integral field spectrograph operating behind the adaptive optics system at Keck Observatory and one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP, including characterization of the flux assignment algorithm, exploration of spatial rippling in the reduced data cubes, and improvements to several calibration files, including the rectification matrix, bad-pixel mask, and wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

γ DORADUS PULSATIONS IN THE ECLIPSING BINARY STAR KIC 6048106

We present the Kepler photometry of KIC 6048106, which is exhibiting the O’Connell effect and multiperiodic pulsations. Including a starspot on either of the components, light-curve synthesis indicates that this system is a semi-detached Algol with a mass ratio of 0.211, an orbital inclination of 73.°9, and a large temperature difference of 2534 K. To examine in detail both the spot variations and pulsations, we separately analyzed the Kepler time-series data at the interval of an orbital period in an iterative way. The results reveal that the variable asymmetries of the light maxima can be interpreted as the changes with time of a magnetic cool spot on the secondary component. Multiple frequency analyses were performed in the outside-eclipse light residuals after removal of the binarity effects from the observed Kepler data. We detected 30 frequencies with signal to noise amplitude ratios larger than 4.0, of which six ( f {sub 2}– f {sub 6} and f {sub 10}) can be identified as high-order (17 ≤ n ≤ 25) low-degree ( ℓ = 2) gravity-mode pulsations that were stable during the observing run of 200 days. In contrast, the other frequencies may be harmonic and combination terms. For the six frequencies, the pulsation periods and pulsation constants are in the ranges of 0.352–0.506 days and 0.232–0.333 days, respectively. These values and the position on the Hertzsprung–Russell diagram demonstrate that the primary star is a γ Dor variable. The evolutionary status and the pulsation nature of KIC 6048106 are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Cloud Structure of Three Galactic Infrared Dark Star-forming Regions from Combining Ground- and Space-based Bolometric Observations

We have modified the iterative procedure introduced by Lin et al., to systematically combine the submillimeter images taken from ground-based (e.g., CSO, JCMT, APEX) and space (e.g., Herschel, Planck) telescopes. We applied the updated procedure to observations of three well-studied Infrared Dark Clouds (IRDCs): G11.11−0.12, G14.225−0.506, and G28.34+0.06, and then performed single-component, modified blackbody fits to each pixel to derive ∼10″ resolution dust temperature and column density maps. The derived column density maps show that these three IRDCs exhibit complex filamentary structures embedded with rich clumps/cores. We compared the column density probability distribution functions (N-PDFs) and two-point correlation (2PT) functions of the column density field between these IRDCs with several OB-cluster-forming regions. Based on the observed correlation between the luminosity-to-mass ratio and the power-law index of the N-PDF, and complementary hydrodynamical simulations for a 10{sup 4} M{sub ⊙} molecular cloud, we hypothesize that cloud evolution can be better characterized by the evolution of the (column) density distribution function and the relative power of dense structures as a function of spatial scales, rather than merely based on the presence of star-forming activity. An important component of our approach is to provide a model-independent quantification of cloud evolution. Based on the small analyzed sample, we propose four evolutionary stages, namely, cloud integration, stellar assembly, cloud pre-dispersal, and dispersed cloud. The initial cloud integration stage and the final dispersed cloud stage may be distinguished from the two intermediate stages by a steeper than −4 power-law index of the N-PDF. The cloud integration stage and the subsequent stellar assembly stage are further distinguished from each other by the larger luminosity-to-mass ratio (>40 L{sub ⊙}/M{sub ⊙}) of the latter. A future large survey of molecular clouds with high angular resolution may establish more precise evolutionary tracks in the parameter space of N-PDF, 2PT function, and luminosity-to-mass ratio.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Solution of Heliospheric Propagation: Unveiling the Local Interstellar Spectra of Cosmic-ray Species

Local interstellar spectra (LIS) for protons, helium, and antiprotons are built using the most recent experimental results combined with state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic-ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod, providing modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameter optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997 to 2015. The proposed LIS accommodate both the low-energy interstellar CR spectra measured by Voyager 1 and the high-energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

79 ASTRONOMY AND ASTROPHYSICS↗