Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalization method”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dedicated beam position monitor pair for model-independent lattice characterization at NSLS-II

This paper reports recent lattice characterization results obtained at the National Synchrotron Light Source II (NSLS-II) storage ring, conducted without reliance on a lattice model. A pair of beam position monitors (BPMs) with bunch-by-bunch (B×B) resolution, were recently installed in a section of the storage ring free of magnetic fields. The new BPM pair measured the beam, or bunch’s transverse Poincar´e map precisely after the beam was excited. Linear one-turn-matrices (OTM) were then derived, and from these, the 4-dimensional coupled Twiss parameters were extracted at the locations of the BPM pair. By normalizing beam oscillation amplitudes with the Twiss parameters, the global action-variables were obtained. These action-variables facilitated the measurement of the local Twiss parameters observed by other BPMs independent on lattice model. This method is general, and particularly useful in certain scenarios such as a round beam mode in a diffraction-limited light source ring. We applied it to assess both weakly and strongly coupled lattices at the NSLS-II ring. Through analysis of the strongly coupled lattice, the quadrupole tilt errors were estimated to be less than 400 µrad. Utilizing the BPMs’ B×B resolution, for the first time we observed the variations of the linear lattice along a long bunch-train.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

NSRDB (National Solar Radiation Database Software) [SWR-23-77]

This repository contains all of the methods for the NSRDB data processing pipeline. The NSRDB is a serially complete collection of hourly and half-hourly values of meteorological data and the three most common measurements of solar radiation: global horizontal, direct normal and diffuse horizontal irradiance.

Benton, Brandon↗
DOE OSTI2020

A method for characterization of multiple dynamic constitutive parameters of FRCs

We propose a method to measure multiple dynamic material constitutive parameters of unidirectional fiber reinforced composites (FRCs) in a single experiment. Dynamic short-beam shear (DSBS) experiments were performed on a modified Kolsky compression bar, with integration of high-speed imaging and digital image correlation (DIC). The unidirectional FRCs investigated were S-2 glass fiber reinforced matrix of TGDDM-Jeffamine® D230 with monoamine functionalized partially reacted substructures (mPRS) and commercially available SC-15. Analytical solutions of normal and shear strains of a composite beam were derived based on Timoshenko beam theory, assuming material to be transversely isotropic and have different moduli in tension and compression in each principle material orientation. Tensile and compressive moduli were inversely computed through monitoring normal strain slope when specimen was constantly loaded at a speed of ~7.3 m/s within a small deflection. Non-linear shear stress-strain behavior of the composite was described via Ramberg–Osgood equation. Finite element (FE) analysis was conducted in ABAQUS, simultaneously defining via user subroutine UMAT the transverse isotropy of material, bi-modulus constitutive model, and non-linear shear stress-strain relation. The method proposed in this work was validated by comparing strain distributions computed by FE model and DIC measurements. Comparing with traditional dynamic tensile, compressive, and shear experiments on FRCs, this method significantly simplifies the specimen preparation and design of complicated gripping fixtures for multiple experiments. Here, systematic errors resulting from variations of specimen geometry and dimension, loading direction, and instrumentation are reduced, thereby providing compatible data for numerical studies on impact behavior of composites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The pp → W(→ lν) + γ process at next-to-next-to-leading order

We present details of the calculation of the pp → W(→ lν)γ process at next-to-next-to-leading order in QCD, calculated using the jettiness slicing method. The calculation is based entirely on analytic amplitudes. Because of the radiation zero, the NLO QCD contribution from the gq channel is as important as the contribution from the Born $ q\overline{q} $ process, disrupting the normal counting of leading and sub-leading contributions. We also assess the importance of electroweak (EW) corrections, including the EW corrections to both the six-parton channel 0 →$ \overline{u} d\nu {e}^{+}\gamma g $ and the five-parton channel 0 →$ \overline{u} d\nu {e}^{+}\gamma $. Previous experimental results have been shown to agree with theoretical predictions, taking into account the large experimental errors. With the advent of run II data from the LHC, the statistical errors on the data will decrease, and will be competitive with the error on theoretical predictions for the first time. We present numerical results for $ \sqrt{s} $ = 7 and 13 TeV. Analytic results for the one-loop six-parton QCD amplitude and the tree-level seven-parton QCD amplitude are presented in appendices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Dedicated beam position monitor pair for model-independent lattice characterization at NSLS-II

This paper reports recent lattice characterization results obtained at the National Synchrotron Light Source II (NSLS-II) storage ring, conducted without reliance on a lattice model. A pair of beam position monitors (BPMs) with bunch-by-bunch (B$\times$B) resolution, were recently installed in a section of the storage ring free of magnetic fields. The new BPM pair measured the beam, or bunch’s transverse Poincaré map precisely after the beam was excited. Linear one-turn-matrices (OTM) were then derived, and from these, the 4-dimensional coupled Twiss parameters were extracted at the locations of the BPM pair. By normalizing beam oscillation amplitudes with the Twiss parameters, the global action-variables were obtained. Additionally, these action-variables facilitated the measurement of the local Twiss parameters observed by other BPMs independent on lattice model. This method is general, and particularly useful in certain scenarios such as a round beam mode in a diffraction-limited light source ring. We applied it to assess both weakly and strongly coupled lattices at the NSLS-II ring. Through analysis of the strongly coupled lattice, the quadrupole tilt errors were estimated to be less than 400 μrad. Utilizing the BPMs’ B$\times$B resolution, for the first time we observed the variations of the linear lattice along a long bunch-train.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vibrational hierarchy leads to dual-phonon transport in low thermal conductivity crystals

Many low-thermal-conductivity ( κ L ) crystals show intriguing temperature ( T ) dependence of κ L : κ L ∝ T -1 (crystal-like) at intermediate temperatures whereas weak T -dependence (glass-like) at high temperatures. It has been in debate whether thermal transport can still be described by phonons at the Ioffe-Regel limit. In this work, we propose that most phonons are still well defined for thermal transport, whereas they carry heat via dual channels: normal phonons described by the Boltzmann transport equation theory, and diffuson-like phonons described by the diffusion theory. Three physics-based criteria are incorporated into first-principles calculations to judge mode-by-mode between the two phonon channels. Case studies on La 2 Zr 2 O 7 and Tl 3 VSe 4 show that normal phonons dominate low temperatures while diffuson-like phonons dominate high temperatures. Our present dual-phonon theory enlightens the physics of hierarchical phonon transport as approaching the Ioffe-Regel limit and provides a numerical method that should be practically applicable to many materials with vibrational hierarchy.

42 ENGINEERING↗
DOE OSTI2021

Comparison of heavy-ion transport simulations: Mean-field dynamics in a box

Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared among each other and to exact calculations. For BUU-like codes, employing the test particle method, the results depend on the combination of the number of test particles and the spread of the profile functions that weight integration over space. These parameters can be properly adapted to give a good reproduction of the analytical zero-sound features. QMD-like codes, using molecular dynamics methods, are characterized by large damping effects, attributable to the fluctuations inherent in their phase-space representation. Moreover, for a given nuclear effective interaction, they generally lead to slower density oscillations, as compared to BUU-like codes. The latter problem is mitigated in the more recent lattice formulation of some of the QMD codes. Furthermore, the significance of these results for the description of real heavy-ion collisions is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantifying dispersity in size and shape of nanoparticles from small-angle scattering data using machine learning based CREASE

Here, we use machine learning (ML) enhanced computational reverse engineering analysis of scattering experiments (CREASE) to interpret small-angle X-ray scattering (SAXS) data obtained from a system of nanoparticles without a priori knowledge of their exact shapes (e.g. spheres or ellipsoids), sizes (0.5–50 nm) and distributions. The SAXS measurements yielded three categories of scattering profiles exhibiting 'strong', 'weak' and 'no' features. Diminishing features (e.g. broadening or disappearing peaks) in scattering profiles have always been attributed to the presence of significant dispersity in the system. Such featureless SAXS data are not suitable for traditional analysis using analytical models. If one were to fit a relevant analytical model (e.g. the lmfit analytical model for polydisperse spheres) to these 'weak' and 'no' SAXS profiles from our nanoparticle systems, one would obtain non-unique interpretations of the data. Relying on electron microscopy to identify the distributions of nanoparticle shapes and sizes is also unfeasible, especially in high-throughput synthesis and characterization loops. In such situations, to identify the distributions of particle sizes and shapes that could be present in the sample, one must rely on methods like ML-CREASE to interpret the data quickly and output all relevant interpretations about the structure present in the system. The ML-CREASE optimization loop takes the experimental scattering profile as input and outputs multiple candidate solutions whose computed scattering profiles match the SAXS profile input. The ML-CREASE method outputs distributions of relevant structural features, such as the volume fraction of the nanoparticles in the system and the mean and standard deviation of the particle size and aspect ratio, assuming a type of distribution (e.g. normal, log-normal) for size and aspect ratio. We find that, for the SAXS profiles analyzed here, accounting for the shape dispersity along with size dispersity of the nanoparticles using ML-CREASE improved the match between the computed scattering profiles and input experimental profiles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Compensatory reserve detects subclinical shock with more expeditious prediction for need of life‐saving interventions compared to systolic blood pressure and blood lactate

Abstract Introduction We conducted a prospective observational study on 205 trauma patients at a level I trauma facility to test the hypothesis that a compensatory reserve measurement (CRM) would identify higher risk for progression to shock and/or need a life‐saving interventions (LSIs) earlier than systolic blood pressure (SBP) and blood lactate (LAC). Methods A composite outcome metric included blood transfusion, procedural LSI, and mortality. Discrete measures assessed as abnormal (ab) were SBP <90 mmHg, CRM <60%, and LAC >2.0. A graded categorization of shock was defined as: no shock (normal [n] SBP [n‐SBP], n‐CRM, n‐LAC); sub‐clinical shock (ab‐CRM, n‐SBP, n‐LAC); occult shock (n‐SBP, ab‐CRM, ab‐LAC); or overt shock (ab‐SBP, ab‐CRM, ab‐LAC). Results Three patients displayed overt shock, 53 displayed sub‐clinical shock, and 149 displayed no shock. After incorporating lactate into the analysis, 86 patients demonstrated no shock, 25 were classified as sub‐clinical shock, 91 were classified as occult shock, and 3 were characterized as overt shock. Each shock subcategory revealed a graded increase requiring LSI and transfusion. Initial CRM was associated with progression to shock (odds ratio = 0.97; p < .001) at an earlier time than SBP or LAC. Conclusions Initial CRM uncovers a clinically relevant subset of patients who are not detected by SBP and LAC. Our results suggest CRM could be used to more expeditiously identify injured patients likely to deteriorate to shock, with requirements for blood transfusion or procedural LSI.

Convertino, Victor A.↗
DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Final Report for Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity

Characterizing in situ stress is essential in Enhanced Geothermal System (EGS) for both risk assessment and operation design. Standard stress measurement techniques such as diagnostic fracture injection tests (DFITs) require creating hydraulic fractures and apply analytical methods such as the G -function analysis to estimate the stress. These classical approaches rely on several assumptions, including but not limited to: (i) the induced fracture is planar and persistent, and (ii) the fracture plane is normal to the minimum horizontal stress. However, these assumptions may not hold in EGS environments like Utah FORGE, where thermo-hydro-mechanical (THM) coupling, wellbore deviation, and rock heterogeneity can alter the local stress condition and consequently lead to fracture complexities in both near-well and far-field regions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Micro-structuration effects on local magneto-transport in [Co/Pd]IrMn thin films

We measured the local magneto-transport (MT) signal with an out-of-plane magnetic field, including magneto-resistance (MR) and Extraordinary Hall effect (EHE), in exchange-biased [Co/Pd]IrMn thin multilayers that are micro-structured with a 100 μm window. We found that when measured locally around the window, the MT signal deviate from the expected behavior. We studied possible causes, including film micro-structuration, electrical contact geometry as well as magnetic field angular tilt. We found that tilting the magnetic field direction with respect to the normal direction does not significantly affect the MT signal, whereas the positioning and geometry of the contacts seem to highly affect the MT signal. For comparison purposes, we carried these MT measurements using the Van-der-Pauw method on a set of four microscopic contacts directly surrounding the window, and on another set of micro-contacts located outside the window, as well as a set of four contacts positioned several millimeters away of each other at the corners of the wafer. If the contacts are sufficiently far apart, the EHE and MR signals have the expected shape and are not significantly affected by the presence of the window. If, on the other hand, the contacts are micro-positioned, the shape of the EHE signal is drastically deformed, and may be modeled as a mix of the standard EHE and MR signals measured on the outer contacts. Furthermore, if the micro-contacts are located directly around the window, the deformation is amplified, and the weight of the MR signal in the mix is further increased by about 40 %. This suggests that the electron path in the Hall geometry is disturbed by both the proximity of the electrodes and by the presence of the window, which both contribute to the deformation for about two-third and one third, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2025

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗
DOE OSTI2024

Radioisotope production using lasers: From basic science to applications

The discovery of chirped pulse amplification has led to great improvements in laser technology, enabling energetic laser beams to be compressed to pulse durations of tens of femtoseconds and focused to a few micrometers. Protons with energies of tens of MeV can be accelerated using, for instance, target normal sheath acceleration and focused on secondary targets. Under such conditions, nuclear reactions can occur, with the production of radioisotopes suitable for medical application. The use of high-repetition lasers to produce such isotopes is competitive with conventional methods mostly based on accelerators. In this paper, we study the production of 67Cu, 63Zn, 18F, and 11C, which are currently used in positron emission tomography and other applications. At the same time, we study the reactions 10B(p,α)7Be and 70Zn(p,4n)67Ga to put further constraints on the proton distributions at different angles, as well as the reaction 11B(p,α)8Be relevant for energy production. The experiment was performed at the 1 PW laser facility at Vega III in Salamanca, Spain. Angular distributions of radioisotopes in the forward (with respect to the laser direction) and backward directions were measured using a high purity germanium detector. Our results are in reasonable agreement with numerical estimates obtained following the approach of Kimura and Bonasera [Nucl. Instrum. Methods Phys. Res., Sect. A 637, 164–170 (2011)].

11C↗