Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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 955 records · Page 53

Development of Digital Twin Predictive Model for PWR Components: Updates on Multi Times Series Temperature Prediction Using Recurrent Neural Network, DMW Fatigue Tests, System Level Thermal-Mechanical-Stress Analysis

The long-term operation (LTO) of nuclear power plant (NPP) beyond their original design life of 40 years, can lead to more material damage associated with cyclic fatigue under thermal-mechanical loading cycles and associated long-term exposure of reactor material to the deleterious reactor-coolant environments. However, under this LTO condition the reactor components can still safely operate but may require more frequent Nondestructive Evaluation (NDE) of reactor components. Frequent NDE requirement may lead to frequent shutdown of the NPP. This in turn can lead to power outage and additional NDE-inspection-cost related economic loss. The economic loss can be minimized by reducing uncertainty in life estimation of safety-critical pressure boundary components and by implementing more digital approach such as by using upcoming digital-twin (DT) technology for predicting the structural states (e.g., time and location dependent inside/outside thickness temperature, stress, strain, plastic deformation, etc.) and associated fatigue life of a component in real time. Towards this goal Argonne National Laboratory (ANL) with the sponsorship of DOE Light Water Reactor Sustainability (LWRS) program is working on the development of a DT framework that can be used for real time environmental fatigue prediction of reactor components. The DT framework is based on limited experiment-data, Artificial-intelligence (AI) – Machine-Learning (ML) - Deep-Learning (DL) based techniques and Multiphysics-computational-mechanics such as finite element (FE) based modeling tools. Towards this overall goal, following are some of the major contributions made during the FY21: 1) Multiple 82/182 dissimilar metal weld (DMW) specimens (both solid-weld and joint-weld representing the actual reactor multi-metal nozzles) were fatigue tested. The resulting fatigue lives were compared to the NUREG-6909 based best-fit and design fatigue curves. Additionally, the results of 52/152 DMW fatigue specimens (which were recently tested at Republic of Korea under the sponsorship of International Nuclear Energy Research Initiative - INERI program) were compared to the NUREG-6909 based best-fit and design fatigue curves. From the comparison of 82/182 and 52/152 DMW test data with NUREG-6909 best-fit curve, most of the reported test data fall way away from the NUREG-6909 suggested best-fit or mean curve. The NUREG-6909 suggested best-fit curve is the best-fit curve of austenitic stainless steel and due to lack of enough data on Nickel-based welds, this is currently being used for predicting the life of Nickel-alloy-based welded components. However, the above observation may require higher scaling factor (e.g., ASME suggested factor of 20 on cycles rather than the current NUREG-6909 suggested factor of 12 on cycles) for scaling the austenitic-stainless-steel best-fit-curve for estimating the design or safe-life of a welded component. Accordingly, for example, if a DMW component experience a strain amplitude of 0.6% the PWR-water life of the component would be 52 cycles instead of 85 cycles. However, more DMW tests are required to further ascertain the above-mentioned observations. 2) A system level CAD and finite element model were developed which consists of reactor pressure vessel (RPV), part of steam generator (SG), part of pressurizer (PRZ), hot leg (HL), and surge line (SL). This is with detailed nozzle geometry and thermal-mechanical material properties of different metals to simulate realistic thermal-mechanical stress under connected system global thermal-mechanical boundary conditions. 3) Different system level heat transfer analyses were performed with estimation of relevant heat transfer coefficients. The resulting data were used in subsequent system level thermal-mechanical stress analysis and for generating spatial-temporal training and validation data for a system level digital-twin based temperature predictor. Transient heat transfer analyses were performed considering thermal boundary condition under design-basis (DB) loading and EDF (Électricité de France) data-based grid-load-following (EDF-GLF) loading cycles. 4) System level thermal-mechanical stress analysis was performed for identifying damage-prone hotspots and for future extension of the model for cyclic state prediction. From the system-level model simulation under DB loading cycle it is found that HL and the SL nozzle that connect to the HL can experience significant stress and strain and could be one of the weakest links in the overall reactor coolant system (RCS). 5) An AI/ML based DT model was developed for multi-time-series temperature prediction at any inside/outside thickness locations of PWR pressure boundary components. This is by using Recurrent-neural-network (RNN) and keras machine learning libraries. The RNN model was validated against two laboratory test-based data sets with one obtained through ANL’s in-air fatigue test system and other through PWR-water test loop. The experimentally validated DT model further validated against FE model results to predict thermal scarification related spatialtemporal temperatures at random locations of a component. The well validated DT model was then used for demonstrating spatial-temporal temperature prediction under 100+ years of reactor operation subjected to combined DB, EDF-GLF and randomized grid-load-following (RANDOMGLF) loading Cycles. The expert-elicitation DT model framework was developed assuming field/input/process measurements can be available from a few existing plant sensors and can readily be used by the NPP operators. The above temperature prediction model will feed to the next-step stress analysis model based on which the life of a component can be predicted in realtime, which is one of our future works.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling of Radiolytic Hydrogen Generation of Irradiated Surrogate Aluminum Plates

The hydrogen generation from oxyhydroxide layers under irradiation is an important effect to understand for the long-term storage of aluminum-clad spent nuclear fuel. The original chemical model shows steady production of hydrogen from aluminum oxyhydroxide layers, with no roll-over effect in the curves of the hydrogen concentration until significant pressure has built up within the canister. This study uses a variety of kinetic fits to adjust the old chemical model to have roll-over effects that are consistent with new high-dose experiments. The data from aluminum 6061 coupon tests at high doses were used to fit several kinetic models. The 6061 aluminum samples are representative of the aluminum-clad fuel from the ATR. The fitted kinetic models were then used in a model of a mini-canister with much higher packing density to determine which kinetic rates would fit both sets of experimental tests. The rate constants for the kinetic mechanisms that were close were then tweaked such that the chemical models fit both the experimental data from the capsule tests and the data from the mini-canister tests. The kinetic models which fit the best are either inhibition by hydrogen, or a back reaction wherein H radical is reabsorbed onto the oxyhydroxide layer. The H back reaction rolls over quicker and may fit future collected data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of ATS photographs using a specially designed electronic console

Measurements of cloud motion, cloud size, and growth rate were made on four days in the vicinity of rawin stations in the BOMEX network. These measurements were compared with the rawins and levels of best fit ranging from 980 to 250 mb were found. No correlation was found between cloud size and level of best fit. The magnitude of the vector difference between the cloud motion and wind at the level of best fit was as high as 15 knots for clouds with levels of best fit below the 500-mb level and 8 knots for clouds with levels of best fit above the 500-mb level. When cloud motion vectors with levels of best fit at or above the 500-mb level were compared with winds at the 300-mb level, the vector difference increased to as much as 16 knots or double the greatest minimum vector difference. A comparison between cloud motions and winds at the level of best fit showed that the clouds tended to move slightly faster than the observed winds and to move slightly to the left of the wind direction. Rate of growth appeared to influence whether clouds moved faster or slower than the wind. No relationship between cloud size and amount of deviation from the wind direction was observed.

Blackmer, R. H., Jr.↗

On the cost of approximating and recognizing a noise perturbed straight line or a quadratic curve segment in the plane

Approximation of noisy data in the plane by straight lines or elliptic or single-branch hyperbolic curve segments arises in pattern recognition, data compaction, and other problems. The efficient search for and approximation of data by such curves were examined. Recursive least-squares linear curve-fitting was used, and ellipses and hyperbolas are parameterized as quadratic functions in x and y. The error minimized by the algorithm is interpreted, and central processing unit (CPU) times for estimating parameters for fitting straight lines and quadratic curves were determined and compared. CPU time for data search was also determined for the case of straight line fitting. Quadratic curve fitting is shown to require about six times as much CPU time as does straight line fitting, and curves relating CPU time and fitting error were determined for straight line fitting. Results are derived on early sequential determination of whether or not the underlying curve is a straight line.

Cooper, D. B.↗

Biofilm monitoring coupon system and method of use

An apparatus and method is disclosed for biofilm monitoring of a water distribution system which includes the mounting of at least one fitting in a wall port of a manifold in the water distribution system with a passage through the fitting in communication. The insertion of a biofilm sampling member is through the fitting with planar sampling surfaces of different surface treatment provided on linearly arrayed sample coupons of the sampling member disposed in the flow stream in edge-on parallel relation to the direction of the flow stream of the manifold under fluid-tight sealed conditions. The sampling member is adapted to be aseptically removed from or inserted in the fitting and manifold under a positive pressure condition and the fitting passage sealed immediately thereafter by appropriate closure means so as to preclude contamination of the water distribution system through the fitting. The apparatus includes means for clamping the sampling member and for establishing electrical continuity between the sampling surfaces and the system for minimizing electropotential effects. The apparatus may also include a plurality of fittings and sampling members mounted on the manifold to permit extraction of the sampling members in a timed sequence throughout the monitoring period.

Sauer, Richard L.↗

Meteoroid and debris special investigation group; status of 3-D crater analysis from binocular imagery

During the 3 month deintegration of the LDEF, the M&D SIG generated approximately 5000 digital color stereo image pairs of impact related features from all space exposed surfaces. Currently, these images are being processed at JSC to yield more accurate feature information. Work is currently underway to determine the minimum number of data points necessary to parametrically define impact crater morphologies in order to minimize the man-hour intensive task of tie point selection. Initial attempts at deriving accurate crater depth and diameter measurements from binocular imagery were based on the assumption that the crater geometries were best defined by paraboloid. We made no assumptions regarding the crater depth/diameter ratios but instead allowed each crater to define its own coefficients by performing a least-squares fit based on user-selected tiepoints. Initial test cases resulted in larger errors than desired, so it was decided to test our basic assumptions that the crater geometries could be parametrically defined as paraboloids. The method for testing this assumption was to carefully slice test craters (experimentally produced in an appropriate aluminum alloy) vertically through the center resulting in a readily visible cross-section of the crater geometry. Initially, five separate craters were cross-sectioned in this fashion. A digital image of each cross-section was then created, and the 2-D crater geometry was then hand-digitized to create a table of XY position for each crater. A 2nd order polynomial (parabolic) was fitted to the data using a least-squares approach. The differences between the fit equation and the actual data were fairly significant, and easily large enough to account for the errors found in the 3-D fits. The differences between the curve fit and the actual data were consistent between the caters. This consistency suggested that the differences were due to the fact that a parabola did not sufficiently define the generic crater geometry. Fourth and 6th order equations were then fitted to each crater cross-section, and significantly better estimates of the crater geometry were obtained with each fit. Work is presently underway to determine the best way to make use of this new parametric crater definition.

Sapp, Clyde A.↗

Modeling ice streams: Derived quantities

The model addressed is a finite-element, map-plane, time-dependent, column-averaged continuity equation solver. The key to the fitting process involves the balance between ice motion dominated by flow in internal layers, and ice motion dominated by sliding at the bed. The fitting process involves an iterative process carried out in the time domain. Beginning with the portion of the ice sheet being modeled identical to the present ice sheet with uniform flow, sliding, and fraction specified at nominal values, the model monitors each nodal point surface elevation. As the calculated surface elevation deviates from the present surface, a correction proportional to the difference is applied to selected parameter sets. This correction is in a sense that would tend to improve the fit at the particular nodal point. A calculated surface elevation that was higher than the present surface would result in an increased fraction, which would tend to lower the calculated surface (if the flow or sliding constant were being used as the fitting parameter, they would be lowered to improve the fit). This process is allowed to proceed as long as is necessary for the situation to stabilize. Typically, this takes tens of thousands of model years, but the rate is dependent on other external forcings such as the accumulation rate. The primary result is that while a typical sample of ice streams from around Antarctica can be fitted quite reasonably using only the fraction of the velocity due to sliding, a different mechanism seems to be in play along the Siple Coast, where reduced sliding constants are required to attain a reasonable fit. Flow is more strongly channelized in this region, and velocities are, in general, higher than are observed in other regions. It is unlikely that the mechanism that controls the ice movement along the Siple Coast is exactly similar to the mechanisms in the other ice streams. The concept of deformable sediments and their contribution to the fast flow along the Siple Coast may have limited applicability to other Antarctic ice streams.

Fastook, James↗

Responses to LBNP in men with varying profiles of strength and aerobic capacity: Implications for flight crews

Hemodynamic and hormonal responses to lower-body negative pressure (LBNP) were examined in 24 healthy men to test the hypothesis that responsiveness of reflex control of blood pressure during orthostatic stress is associated with strength and/or aerobic capacity. Subjects underwent treadmill tests to determine peak oxygen uptake (peak VO2) and isokinetic dynamo meter tests to determine leg strength. Based on predetermined criteria, the subjects were classified into one of four fitness profiles of six subjects each matched for age, height, and weight: (1) low strength/low aerobic fitness; (2) low strength/high aerobic fitness; (3) high strength/low aerobic fitness; and (4) high strength/high aerobic fitness. Following 90 min of 6 degree head-down tilt (HDT), each subject underwent graded LBNP through -50 mmHg or presyncope, with maximal duration 15 min. All groups exhibited typical hemodynamic, hormonal, and fluid shift responses during LBNP, with no intergroup differences except for catecholamines. Seven subjects, distributed among the four fitness profiles, became presyncopal. Subjects who showed greatest reduction in mean arterial pressure (MAP) during LBNP had greater elevations in vasopressin and lesser increases in heart rate and peripheral resistance. Peak VO2 nor leg strength were correlated with fall in MAP or with syncopal episodes. We conclude that neither aerobic nor strength fitness characteristics are good predictors of responses to LBNP stress.

Convertino, Victor A.↗

Hemodynamic and hormonal responses to lower body negative pressure in men with varying profiles of strength and aerobic power

Hemodynamic, cardiac, and hormonal responses to lower-body negative pressure (LBNP) were examined in 24 healthy men to test the hypothesis that responsiveness of reflex control of blood pressure during orthostatic challenge is associated with interactions between strength and aerobic power. Subjects underwent treadmill tests to determine peak oxygen uptake (VO2max) and isokinetic dynamometer tests to determine knee extensor strength. Based on predetermined criteria, subjects were classified into one of four fitness profiles of six subjects each, matched for age, height, and body mass: (a) low strength/average aerobic fitness, (b) low strength/high aerobic fitness, (c) high strength/average aerobic fitness, and (d) high strength/high aerobic fitness. Following 90 min of 0.11 rad (6 degrees) head-down tilt (HDT), each subject underwent graded LBNP to -6.7 kPa or presyncope, with maximal duration 15 min, while hemodynamic, cardiac, and hormonal responses were measured. All groups exhibited typical hemodynamic, hormonal, and fluid shift responses during LBNP, with no intergroup differences between high and low strength characteristics. Subjects with high aerobic power exhibited greater (P < 0.05) stroke volume and lower (P < 0.05) heart rate, vascular peripheral resistance, and mean arterial pressure during rest, HDT, and LBNP. Seven subjects, distributed among the four fitness profiles, became presyncopal. These subjects showed greatest reduction in mean arterial pressure during LBNP, had greater elevations in vasopressin, and lesser increases in heart rate and peripheral resistance. Neither VO2max nor leg strength were associated with fall in arterial pressure or with syncopal episodes. We conclude that interactions between aerobic and strength fitness characteristics do not influence responses to LBNP challenge.

Norepinephrine/blood↗

Multilayer Composite Pressure Vessels

A method has been devised to enable the fabrication of lightweight pressure vessels from multilayer composite materials. This method is related to, but not the same as, the method described in gMaking a Metal- Lined Composite-Overwrapped Pressure Vessel h (MFS-31814), NASA Tech Briefs, Vol. 29, No. 3 (March 2005), page 59. The method is flexible in that it poses no major impediment to changes in tank design and is applicable to a wide range of tank sizes. The figure depicts a finished tank fabricated by this method, showing layers added at various stages of the fabrication process. In the first step of the process, a mandrel that defines the size and shape of the interior of the tank is machined from a polyurethane foam or other suitable lightweight tooling material. The mandrel is outfitted with metallic end fittings on a shaft. Each end fitting includes an outer flange that has a small step to accommodate a thin layer of graphite/epoxy or other suitable composite material. The outer surface of the mandrel (but not the fittings) is covered with a suitable release material. The composite material is filament- wound so as to cover the entire surface of the mandrel from the step on one end fitting to the step on the other end fitting. The composite material is then cured in place. The entire workpiece is cut in half in a plane perpendicular to the axis of symmetry at its mid-length point, yielding two composite-material half shells, each containing half of the foam mandrel. The halves of the mandrel are removed from within the composite shells, then the shells are reassembled and bonded together with a belly band of cured composite material. The resulting composite shell becomes a mandrel for the subsequent steps of the fabrication process and remains inside the final tank. The outer surface of the composite shell is covered with a layer of material designed to be impermeable by the pressurized fluid to be contained in the tank. A second step on the outer flange of each end fitting accommodates this layer. Depending on the application, this layer could be, for example, a layer of rubber, a polymer film, or an electrodeposited layer of metal. If the fluid to be contained in the tank is a gas, then the best permeation barrier is electrodeposited metal (typically copper or nickel), which can be effective at a thickness of as little as 0.005 in (.0.13 mm). The electrodeposited metal becomes molecularly bonded to the second step on each metallic end fitting. The permeation-barrier layer is covered with many layers of filament-wound composite material, which could be the same as, or different from, the composite material of the inner shell. Finally, the filament-wound composite material is cured in an ov

DeLay, Tom↗

Modeling Phase-Aligned Gamma-Ray and Radio Millisecond Pulsar Light Curves

Since the discovery of the first eight gamma-ray millisecond pulsars (MSPs) by the Fermi Large Area Telescope, this population has been steadily expanding. Four of the more recent detections, PSR J00340534, PSR J1939+2134 (B1937+21; the first MSP ever discovered), PSR J1959+2048 (B1957+20; the first discovery of a black widow system), and PSR J2214+3000, exhibit a phenomenon not present in the original discoveries: nearly phase-aligned radio and gamma-ray light curves (LCs). To account for the phase alignment, we explore models where both the radio and gamma-ray emission originate either in the outer magnetosphere near the light cylinder or near the polar caps. Using a Markov Chain Monte Carlo technique to search for best-fit model parameters, we obtain reasonable LC fits for the first three of these MSPs in the context of altitude-limited outer gap (alOG) and two-pole caustic (alTPC) geometries (for both gamma-ray and radio emission). These models differ from the standard outer gap (OG)/two-pole caustic (TPC) models in two respects: the radio emission originates in caustics at relatively high altitudes compared to the usual conal radio beams, and we allow both the minimum and maximum altitudes of the gamma-ray and radio emission regions to vary within a limited range (excluding the minimum gamma-ray altitude of the alTPC model, which is kept constant at the stellar radius, and that of the alOG model, which is set to the position-dependent null charge surface altitude). Alternatively, phase-aligned solutions also exist for emission originating near the stellar surface in a slot gap scenario (low-altitude slot gap (laSG) models). We find that the alTPC models provide slightly better LC fits than the alOG models, and both of these give better fits than the laSG models (for the limited range of parameters considered in the case of the laSG models). Thus, our fits imply that the phase-aligned LCs are likely of caustic origin, produced in the outer magnetosphere, and that the radio emission for these pulsars may come from close to the light cylinder. In addition, we were able to constrain the minimum and maximum emission altitudes with typical uncertainties of 30% of the light cylinder radius. Our results therefore describe a third gamma-ray MSP subclass, in addition to the two previously found by Venter et al.: those with LCs fit by standard OG/TPC models and those with LCs fit by pair-starved polar cap models.

Venter, C.↗

Foreground Bias from Parametric Models of Far-IR Dust Emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB plus foreground emission to precision 0.1 percent or better. The Rayleigh-Jeans approximation to the dust spectrum biases the fitted dust spectral index by (Delta)(Beta)(sub d) = 0.2 and the inflationary B-mode amplitude by (Delta)(r) = 0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by (Delta)(r) greater than 0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level (Delta)(r) = 0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

cosmic microwave background↗

Exploring Connectivity in Sequence Space of Functional RNA

Emergence of replicable genetic molecules was one of the marking points in the origin of life, evolution of which can be conceptualized as a walk through the space of all possible sequences. A theoretical concept of fitness landscape helps to understand evolutionary processes through assigning a value of fitness to each genotype. Then, evolution of a phenotype is viewed as a series of consecutive, single-point mutations. Natural selection biases evolution toward peaks of high fitness and away from valleys of low fitness. whereas neutral drift occurs in the sequence space without direction as mutations are introduced at random. Large networks of neutral or near-neutral mutations on a fitness landscape, especially for sufficiently long genomes, are possible or even inevitable. Their detection in experiments, however, has been elusive. Although a few near-neutral evolutionary pathways have been found, recent experimental evidence indicates landscapes consist of largely isolated islands. The generality of these results, however, is not clear, as the genome length or the fraction of functional molecules in the genotypic space might have been insufficient for the emergence of large, neutral networks. Thorough investigation on the structure of the fitness landscape is essential to understand the mechanisms of evolution of early genomes. RNA molecules are commonly assumed to play the pivotal role in the origin of genetic systems. They are widely believed to be early, if not the earliest, genetic and catalytic molecules, with abundant biochemical activities as aptamers and ribozymes, i.e. RNA molecules capable, respectively, to bind small molecules or catalyze chemical reactions. Here, we present results of our recent studies on the structure of the sequence space of RNA ligase ribozymes selected through in vitro evolution. Several hundred thousands of sequences active to a different degree were obtained by way of deep sequencing. Analysis of these sequences revealed several large clusters defined such that every sequence in a cluster can be reached from any other sequence in the same cluster through a series of single point mutations. Sequences in a single cluster appear to adopt more than one secondary structure. The mechanism of refolding within a single cluster was examined. To shed light on possible evolutionary paths in the space of ribozymes, the connectivity between clusters was investigated. The effect of length of RNA molecules on the structure of the fitness landscape and possible evolutionary paths was examined by way of comparing functional sequences of 20 and 80 nucleobases in length. It was found that sequences of different lengths shared secondary structure motifs that were presumed responsible for catalytic activity, with increasing complexity and global structural rearrangements emerging in longer molecules.

Wei, Chenyu↗

Measurement of neutrino oscillations using neutrino and antineutrino beams in the NOvA experiment

NOvA is a long-baseline accelerator neutrino oscillation experiment using the NuMIneutrino beam from Fermilab. Its main physics goals are to probe the 3-flavour oscillationparameters: neutrino mass hierarchy, CP-violating phase dcp and octant of .23 mixingangle by observing electron neutrino appearance and muon neutrino disappearance. Twofunctionally identical detectors are placed off-axis from the centre of the NuMI beam.The near detector at Fermilab is 100 m underground, and the far detector is locatedon the surface at Ash River, 810 km away from the beam source. The initial neutrinobeam spectra are measured using the near detector data and the oscillation parametersare extracted by fitting the observed data to the predicted neutrino spectrum in the fardetector.This thesis is centered around how to improve the sensitivity of |.m232| and .23 measurementsin the muon neutrino disappearance analysis. NOvA will take data for about12 years. The operation of the NOvA experiment for each year costs tens of millions ofdollars, thus it is valuable to maximise the sensitivity of the analysis. Three samples ofmuon neutrino events are studied in this thesis to improve the analysis sensitivities. First,higher energy muon neutrinos are investigated by extending the energy range in NOvA’scurrent standard analysis. Second, for the sample of events used in NOvA’s existing analysis,a new energy estimator which has been developed to improve the neutrino energy resolution is considered. Furthermore, in addition to binning the events as function ofenergy and hadronic energy fraction, three particle identifiers are introduced to separateneutrino events by signal purity to reduce the effects from backgrounds. Third, an additionallower purity sample of muon neutrino charged current (CC) events that look similarto neutral current events and have not been included in NOvA’s existing analyses havebeen studied.This thesis reanalyses NOvA’s data used in the 2020 analysis, corresponding to anexposure of 13.60×1020 protons on target (POT) in the neutrino beam mode recordedfrom February 6, 2014 to March 20, 2020, and 12.50×1020 protons on target in theantineutrino beam mode recorded between June 29, 2016 to February 26, 2019. Thisthesis has implemented a fit to Asimov fake data, generated where sin2 .23 = 0.59 and.m232 = 2.5 × 10-3 eV2. These sensitivity studies show that the uncertainty range of|.m232| at 1 s in the new analysis is reduced by 5.5% and the significance of maximaldisappearance rejection improves by 7.7%, compared to the standard analysis. This isequivalent to adding 11-16% more data. The best fit values of the oscillation parametersfrom fitting to the far detector (FD) data with the new analysis are found to besin2 .23 = 0.568+0.025-0.043 (sin2 .23 = 0.454+0.046-0.026) and .m232 = 2.399+0.055-0.070 × 10-3 eV2 (.m232= -2.427+0.055-0.067 × 10-3 eV2) for the normal (inverted) hierarchy. The difference in thebest fit for sin2 .23 (.m232) between the new analysis and NOvA’s 2020 analysis is around2% (1.4%). The uncertainty range at 1 s for .m232 decrease by 8% (4%) for the normalhierarchy (inverted hierarchy) compared to the standard analysis. The uncertainty rangefor sin2 .23 is close to the standard analysis. This thesis also implements the fit from combiningelectron neutrino appearance and muon neutrino disappearance. The combinedanalysis shows that the best fit values are very close to the standard analysis. However,the uncertainty range of .m232 at 1 s is reduced by 3.7% using the new analysis. The maximaldisappearance significance is not improved in the new analysis, but the new analysisslightly improves the rejection of the disfavoured octant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimal Stopping Ages for Colorectal Cancer Screening

Importance Prior studies have shown that the benefits, harms, and costs of colorectal cancer (CRC) screening at older ages are associated with a patient’s sex, health, and screening history. However, these studies were hypothetical exercises and not directly informed by data on CRC risk. Objective To identify the optimal stopping ages for CRC screening by sex, comorbidity, and screening history from a cost-effectiveness perspective. Design, Setting, and Participants This economic evaluation first validated the MISCAN-Colon (Microsimulation Screening Analysis–Colon) model against community-based CRC incidence and mortality rates for 2 subcohorts of the PRECISE (Optimizing Colorectal Cancer Screening Precision and Outcomes in Community-Based Populations) cohort. Subsequently, different CRC screening scenarios were simulated in older individuals. Cohorts of US adults aged 76 to 90 years varied by sex and comorbidity status (none, low, moderate, or severe). Statistical and sensitivity analyses were performed from March 2023 to May 2024. Exposures CRC screening histories including fecal immunochemical test (FIT) or colonoscopy, such as a negative colonoscopy result from 10, 15, 20, 25, or 30 years before the index age; 1 to 5 negative FIT results within 5 years of the index age, with different patterns of recency; or a combination of negative colonoscopy and negative FIT results. Main Outcomes and Measures The main outcomes included estimated lifetime clinical outcomes, incremental costs, and quality-adjusted life-years gained (QALYG) associated with 1 additional FIT or colonoscopy. Optimal stopping age for screening, defined as the oldest age for which the incremental cost-effectiveness ratio was still below the willingness-to-pay threshold of $\$$100 000 per QALYG, was evaluated. Results The first of the 2 PRECISE subcohorts used in validating the simulation model included 25 974 adults (15 060 females [58.0%]; 54.7% aged 76 to 80 years) with a negative colonoscopy result 10 years before the index date. The second subcohort consisted of 118 269 adults (67 058 females [56.7%]; 90.5% aged 76 to 80 years) with a negative FIT result 1 year before the index date. Older age, male sex, higher comorbidity levels, and recent CRC screenings were associated with reduced incremental benefit and cost-effectiveness of additional screening. For the reference cohort of 76-year-old females without comorbidities and a negative colonoscopy result 10 years before the index age, 1 additional colonoscopy cost $\$$38 226 per QALYG. For cohorts with otherwise equivalent characteristics, associated costs increased to $\$$1 689 945 per QALYG for females at age 90 years without comorbidities and a negative colonoscopy results 10 years before the index age, $\$$51 604 per QALYG for males at age 76 years without comorbidities and a negative colonoscopy result 10 years before the index age, and $\$$108 480 per QALYG for females at age 76 years with severe comorbidities and a negative colonoscopy result 10 years before the index age and decreased to $\$$16 870 per QALYG for females without comorbidities and a negative colonoscopy result 30 years before the index age. The optimal stopping ages across different cohorts ranged from younger than 76 to 86 years for colonoscopy and younger than 76 to 88 years for FIT. Conclusions and Relevance In this economic evaluation, age, sex, screening history, comorbidity, and future screening modality were associated with the clinical outcomes, cost-effectiveness, and optimal stopping age for CRC screening. These results can inform guideline development and patient-directed informed decision-making.

Harlass, Matthias [Erasmus Erasmus University Medi↗

Interactions of πK, ππK and KKπ systems at maximal isospin from lattice QCD

We study the interactions of systems of two and three nondegenerate mesons composed of pions and kaons at maximal isospin using lattice QCD, specifically $π^+K^+, π^+π^+K^+$ and $K^+K^+π^+$. Utilizing the stochastic LapH method, we determine the spectrum of these systems on two CLS $N_f$ = 2 + 1 ensembles with pion masses of 200 MeV and 340 MeV, and include many levels in different momentum frames. We constrain the K matrices describing two- and three-particle interactions by fitting the spectrum to the results predicted by the finite-volume formalism, including up to p waves. This requires also results for the $π^+π^+$ and $K^+K^+$ spectrum, which have been obtained previously on the same configurations. We explore different fitting strategies, comparing fits to energy shifts with fits to energies boosted to the rest frame, and also comparing simultaneous global fits to all relevant two- and three-particle channels to those where we first fit twoparticle channels and then add in the three-particle information. We provide the first determination of the three-particle K matrix in $π^+π^+K^+$ and $K^+K^+π^+$ systems, finding statistically significant nonzero results in most cases. We include s and $\textit{p}$ waves in the K matrix for $π^+K^+$ scattering, finding evidence for an attractive p-wave scattering length. We compare our results to Chiral Perturbation Theory, including an investigation of the impact of discretization errors, for which we provide the leading order predictions obtained using Wilson Chiral Perturbation Theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗

Interpreting the Operando XANES of Surface-Supported Subnanometer Clusters: When Fluxionality, Oxidation State, and Size Effect Fight.

X-ray absorption near edge structure (XANES) spectroscopy is widely used for operando catalyst characterization. We show that, for highly fluxional supported nanoclusters, the customary extraction of the oxidation state of the metal from the XANES data by fitting to the bulk standards is highly questionable. The XANES signatures as well as the apparent oxidation state for such clusters arise from a complex combination of many factors, and not only from the chemical composition in reaction conditions (e.g., oxygen content in oxidizing atmosphere). The thermally accessible isomerization and population of several structurally distinct cluster forms, cluster-support interaction, and intrinsic size effects all impact the metal oxidation state and XANES signal. We demonstrate this on copper oxide clusters with different compositions, Cu4Ox(x = 2-5) and Cu5Oy (y = 3, 5), deposited on amorphous alumina and ultrananocrystalline diamond, for which we computed the XANES spectra and compare the results to the experiment. We show in addition that fitting the experimental spectrum to calculated spectra of supported clusters can, in contrast, provide good agreement and insight into the spectrum-composition-structure relation. Experimental XANES interpreted using the proposed fitting scheme shows the partial reduction of Cu oxide clusters at rising temperatures, and pinpoints the specific stoichiometries that dominate in the ensemble of cluster states as the temperature changes.

Zandkarimi, Borna↗