Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Decay Correction”

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 145 records · Page 8

Some effects of swirl on turbulent mixing and combustion

A general formulation of some effects of swirl on turbulent mixing is given. The basis for the analysis is that momentum transport is enhanced by turbulence resulting from rotational instability of the fluid field. An appropriate form for the turbulent eddy viscosity is obtained by mixing length type arguments. The result takes the form of a corrective factor that is a function of the swirl and acts to increase the eddy viscosity. The factor is based upon the initial mixing conditions implying that the rotational turbulence decays in a manner similar to that of free shear turbulence. Existing experimental data for free jet combustion are adequately matched by using the modifying factor to relate the effects of swirl on eddy viscosity. The model is extended and applied to the supersonic combustion of a ring jet of hydrogen injected into a constant area annular air stream. The computations demonstrate that swirling the flow could: (1) reduce the burning length by one half, (2) result in more uniform burning across the annulus width, and (3) open the possibility of optimization of the combustion characteristics by locating the fuel jet between the inner wall and center of the annulus width.

Rubel, A.↗

Relativistic Radiative and Auger Rates for Fe XXIV

As part of a project to compute improved atomic data for the spectral modeling of iron K lines, we report extensive calculations and comparisons of radiative and Auger rates for transitions involving the K-vacancy states in Fe XXIV. By making use of several computational codes, a detailed study is carried out of orbital representation, configuration interaction, relativistic corrections, cancellation effects, and fine tuning. It is shown that a formal treatment of the Breit interaction is essential to render the important magnetic correlations that take part in the decay pathways of this ion. As a result, the accuracy of the present A-values is firmly ranked at better than 10% while that of the Auger rates at only 15%.

Bautista, M. A.↗

Transient hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation. II. Isomer-resolved unimolecular dynamics

Transient carbon-centered hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation are characterized by their time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) and cyclic ether products. Two distinct •QOOH isomers are examined with radical sites at a primary carbon of one of the methyl groups (β-Me) or a secondary carbon (β-Et) of the ethyl group. Energy-dependent unimolecular rates are obtained from the time-dependent appearance of OH products for the two isomers and compared with statistical microcanonical rates computed using RRKM theory, including heavy-atom tunneling, based on high-level theoretical calculations. A benchmark-corrected approach is utilized to compute high-accuracy stationary-point energies, most importantly, transition-state barriers, for the •QOOH_Me and •QOOH_Et isomers in isopentane oxidation, building on higher-level reference calculations for the oxidation of ethane (C 2 H 5 O 2 ) and propane (C 3 H 7 O 2 ), respectively. Here, the measured rates are compared with RRKM calculations incorporating the benchmark-corrected transition-state parameters, a vibrationally adiabatic multidimensional hindered-rotor treatment of key torsions, and quantum tunneling. Agreement between experiment and theory validates the statistical description and shows faster decay for •QOOH_Et due to its lower barrier. Both β-QOOH isomers decay almost exclusively to OH + cyclic ether products under the conditions studied.

Oxidation processes↗

AutoSAC-Automated Sample Activation Calculation on VULCAN

AutoSAC is an automated python program to conduct sample activation calculation immediately after neutron scattering experiment on VULCAN. AutoSAC takes account of sample dimensions, compositions, forms, sample environments, neutron beam configurations, neutron scattering history, and proton charged resolved beam power etc., and uses a simplified method to rapidly estimate the radioactivity at end of beam or after. It generates a draft report per sample and repository of radioactivity per measurement and provides feedback of current radioactivity and decay time back to a neutron proposal Slack channel to inform users immediately on sample handling. It also utilizes either the sample note field in EPICS or sample composition or thickness correction and takes slack prompt or command to conduct post-beam SAC by demand. The program can be shared to be used in other beamlines at spallation neutron source with modified instrument configurations including instrument specific characteristic beam fluxes.

An, Ke [Oak Ridge National Laboratory (ORNL), Oak ↗

Theoretical evaluation of rigid baffles in the suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs both two and three dimensional combustor models characterized by concentrated combustion sources at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field for both models. Boundary layer corrections to this unsteady flow are in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral instability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss. The response of the dissipation model to varying baffle blade length, mean flow Mach number, oscillation amplitude, baffle configuration, and oscillation mode is examined.

Combustion instability↗

Observations of stratospheric aerosols associated with the El Chichon eruption

Lidar observations of aerosols were carried out at Aberystwyth between Nov. 1982 and Dec. 1985 using a frequency doubled and frequency tripled Nd/Yag laser and a receiver incorporating a 1 m diameter in a Newtonian telescope configuration. In analyses of the experimental data attention is paid to the magnitude of the coefficient relating extinction and backscatter, the choice being related to the possible presence of aerosols in the upper troposphere and the atmospheric densities employed in the normalisation procedure. The aerosol loading showed marked day to day changes in early months and an overall decay was apparent only after April 1983, this decay being consistent with an e sup -1 time of about 7 months. The general decay was accompanied by a lowering of the layer but layers of aerosols were shown intermittently at heights above the main layer in winter months. The height variations of photon counts corrected for range, or of aerosol backscatter ratio, showed clear signatures of the tropopause. A strong correlation was found between the heights of the tropopause identified from the lidar measurements and from radiosonde-borne temperature measurements. A notable feature of the observations is the appearance of very sharp height gradients of backscatter ratio which seem to be produced by differential advection.

Thomas, L.↗

Heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(s)

The heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) (eq 1) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(S) (eq 2) were investigated over the temperature range 223-296 K in a flow-tube reactor coupled to a quadrupole mass spectrometer. Either a chemical ionization mass spectrometer (CIMS) or an electron-impact ionization mass spectrometer (EIMS) was used to provide suitable detection sensitivity and selectivity. In order to mimic atmospheric conditions, partial pressures of HNO3 and N2O5 in the range 6 x 10(exp -8) - 2 x 10(exp -6) Torr were used. Granule sizes and surface roughness of the solid NaCl substrates were determined by using a scanning electron microscope. For dry NaCl substrates, decay rates of HNO3 were used to obtain gamma(1) = 0.013 +/- 0.004 (1sigma) at 296 K and > 0.008 at 223 K, respectively. The error quoted is the statistical error. After all corrections were made, the overall error, including systematic error, was estimated to be about a factor of 2. HCl was found to be the sole gas-phase product of reaction 1. The mechanism changed from heterogeneous reaction to predominantly physical adsorption when the reactor was cooled from 296 to 223 K. For reaction 2 using dry salts, gamma(2) was found to be less than 1.0 x 10(exp -4) at both 223 and 296 K. The gas-phase reaction product was identified as ClNO2 in previous studies using an infrared spectrometer. An enhancement in reaction probability was observed if water was not completely removed from salt surfaces, probably due to the reaction of N2O5(g) + H2O(s) yields 2HNO3(g). Our results are compared with previous literature values obtained using different experimental techniques and conditions. The implications of the present results for the enhancement of the hydrogen chloride column density in the lower stratosphere after the El Chichon volcanic eruption and for the chemistry of HCl and HNO3 in the marine troposphere are discussed.

Leu, Ming-Taun↗

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗

A theoretical evaluation of rigid baffles in suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs both two and three dimensional combustor models characterized by concentrated combustion sources at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field for both models. Boundary layer corrections to this unsteady flow are used in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral stability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss.

Combustion instability↗

A theoretical evaluation of rigid baffles in suppression of combustion instability

An analytical technique for the prediction of the effects of rigid baffles on the stability of liquid propellant combustors is presented. This analysis employs a three dimensional combustor model characterized by a concentrated combustion source at the chamber injector and a constant Mach number nozzle. An eigenfunction-matching method is used to solve the linearized partial differential equations describing the unsteady flow field. Boundary layer corrections to this unsteady flow are used in a mechanical energy dissipation model to evaluate viscous and turbulence effects within the flow. An integral stability relationship is then employed to predict the decay rate of the oscillations. Results of this analysis agree qualitatively with experimental observations and show that sufficient dissipation exists to indicate that the proper mechanism of baffle damping is a fluid dynamic loss. The response of the dissipation model to varying baffle blade length, mean flow Mach number and oscillation amplitude is examined.

Combustion instability↗

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗

The Fermi function and the neutron's lifetime

The traditional Fermi function ansatz for nuclear beta decay describes enhanced perturbative effects in the limit of large nuclear charge Z and/or small electron velocity β. We define and compute the quantum field theory object that replaces this ansatz for neutron beta decay, where neither of these limits hold. We present a new factorization formula that applies in the limit of small electron mass, analyze the components of this formula through two loop order, and resum perturbative corrections that are enhanced by large logarithms. We apply our results to the neutron lifetime, supplying the first two-loop input to the long-distance corrections. Our result can be summarized as τ n x |V ud | 2 [1 + 3λ 2 ] [1 + Δ R ] = $\frac{5263.284(17) s}{1 + 27.04(7) x 10^{-3}}$ with |V ud | the up-down quark mixing parameter, τ n the neutron's lifetime, λ the ratio of axial to vector charge, and Δ R the short-distance matching correction. We find a shift in the long-distance radiative corrections compared to previous work, and discuss implications for extractions of |V ud | and tests of the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Calibration of a soft secondary vertex tagger using proton-proton collisions at s = 13 TeV with the ATLAS detector

Several processes studied by the ATLAS experiment at the Large Hadron Collider produce low-momentum b-flavored hadrons in the final state. This paper describes the calibration of a dedicated tagging algorithm that identifies b-flavored hadrons outside of hadronic jets by reconstructing the soft secondary vertices originating from their decays. The calibration is based on a proton-proton collision dataset at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 140 fb -1 . Scale factors used to correct the algorithm’s performance in simulated events are extracted for the b-tagging efficiency and the mistag rate of the algorithm using a data sample enriched in $t\overline{t}$ events. Several orthogonal measurement regions are defined, binned as a function of the multiplicities of soft secondary vertices and jets containing a b-flavored hadron in the event. The mistag rate scale factors are estimated separately for events with low and high average numbers of interactions per bunch crossing. The results, which are derived from events with low missing transverse momentum, are successfully validated in a phase space characterized by high missing transverse momentum and therefore are applicable to new physics searches carried out in either phase space regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Atomic Decay Data for Modeling K Lines of Iron Peak and Light Odd-Z Elements*

Complete data sets of level energies, transition wavelengths, A-values, radiative and Auger widths and fluorescence yields for K-vacancy levels of the F, Na, P, Cl, K, Sc, Ti, V, Cr, Mn, Co, Cu and Zn isonuclear sequences have been computed by a Hartree-Fock method that includes relativistic corrections as implemented in Cowan's atomic structure computer suite. The atomic parameters for more than 3 million fine-structure K lines have been determined. Ions with electron number N greater than 9 are treated for the first time, and detailed comparisons with available measurements and theoretical data for ions with N less than or equal to 9 are carried out in order to estimate reliable accuracy ratings.

modeling K lines↗

Accuracy enhancements for overset grids using a defect correction approach

A defect-correction approach is investigated as a means of enhancing the accuracy of flow computations on overset grids. Typically, overset-grid techniques process and pass information only at grid boundaries. In the current approach, error corrections at all overlapped interior points are injected between grids by using a defect-correction scheme. In some cases this is found to enhance the overall accuracy of the overset-grid method. Locally refined overset grids can be used to provide an efficient solution-adaptation method. The defect correction can also be ultilized as an error-correction technique for a coarse grid by evaluating the residual using a fine base grid, but solving the implicit equations only on the coarse grid. Numerical examples include an accuracy and dissipation study of an unsteady decaying vortex flow, the flow over a NACA 0012 airfoil, and the flow over a mulit-element high-lift airfoil.

Rogers, Stuart E.↗

Long Term Missions at the Sun-Earth Libration Point L1: ACE, SOHO, and WIND

Three heliophysics missions -- the Advanced Composition Explorer (ACE), Solar Heliospheric Observatory (SOHO), and the Global Geoscience WIND -- have been orbiting the Sun-Earth interior libration point L1 continuously since 1997, 1996, and 2004, respectively. ACE and WIND (both NASA missions) and SOHO (an ESA-NASA joint mission) are all operated from the NASA Goddard Space Flight Center (GSFC). While ACE and SOHO have been dedicated libration point orbiters since their launches, WIND has had also a remarkable 10-year career flying a deep-space, multiple lunar-flyby trajectory prior to 2004. That era featured 36 targeted lunar flybys with excursions to both L1 and L2 before its final insertion in L1 orbit. A figure depicts the orbits of the three spacecraft, showing projections of the orbits onto the orthographic planes of a solar rotating ecliptic frame of reference. The SOHO orbit is a quasi-periodic halo orbit, where the frequencies of the in-plane and out-of-plane motions are practically equal. Such an orbit is seen to repeat itself with a period of approximately 178 days. For ACE and WIND, the frequencies of the in-plane and out-of-plane motions are unequal, giving rise to the characteristic Lissajous motion. ACE's orbit is of moderately small amplitude, whereas WIND's orbit is a large-amplitude Lissajous of dimensions close to those of the SOHO halo orbit. As motion about the collinear points is inherently unstable, stationkeeping maneuvers are necessary to prevent orbital decay and eventual escape from the L1 region. Though the three spacecraft are dissimilar (SOHO is a 3-axis stabilized Sun pointer, WIND is a spin-stabilized ecliptic pole pointer, and ACE is also spin-stabilized with its spin axis maintained between 4 and 20 degrees of the Sun), the stationkeeping technique for the three is fundamentally the same. The technique consists of correcting the energy of the orbit via a delta-V directed parallel or anti-parallel to the Spacecraft-to-Sun line. SOHO achieves this using thrusters oriented in line with the solar direction. WIND achieves the delta-V via pulsing radial thrusters when aligned with the Sun. ACE uses axial thrusters to apply delta-V with a component that is 94% or more aligned with the ACE-Sun line. Sunward thrust adds energy to the orbit preventing decay back toward Earth. Thrust directed anti-Sunward takes energy out of the L1 orbit, thereby preventing escape from the Earth-Moon system into independent heliocentric orbit. Libration point orbit stationkeeping delta-V costs grow exponentially with time elapsed from the last maneuver performed. The doubling time constant is approximately 16 days. For the sake of fuel conservation, and for limiting the absolute magnitude of propulsion performance errors, stationkeeping maneuvers should be performed before the delta-V grows too large; for our purposes 'too large' is considered to be greater than 0.5 m/sec. In practice, the typical interval between burns for this trio is about three months, and the typical delta-V is much smaller than 0.5 m/sec. Typical annual stationkeeping costs have been around 1.0 m/sec for ACE and WIND, and much less than that for SOHO. All three spacecraft have ample fuel remaining; barring contingencies all three could, in principle, be maintained at L1 for decades to come. This paper will review the L1 orbits and the mission history of ACE, WIND, and SOHO, and describe the stationkeeping techniques and orbit maneuver experience. The Lissajous phase control that was practiced for ACE during the period from 1999 to 2001 will also be briefly discussed. The final section will consider the future of these ongoing missions.

Roberts, Craig E.↗

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗

Timelike meson form factors beyond the elastic region from lattice QCD

We present a calculation of the vector-isovector timelike form factors of the pion and the kaon using lattice quantum chromodynamics. We calculate two-point correlation functions with m π ∼ 280 MeV , extracting both the finite-volume spectrum and matrix elements for these states created from the vacuum by a vector current. After determining the coupled-channel π π , K K ¯ scattering amplitudes, we perform the necessary correction for the significant finite-volume effects present in the current matrix elements, leading to the timelike form factors. We find these to be dominated by the presence of the ρ resonance, and we extract its decay constant by an analytic continuation of the amplitudes to the resonance pole. In addition, the spacelike pion form factor is determined on the same lattice configurations, and a dispersive parametrization is used to simultaneously describe the spacelike and elastic timelike regions. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗