Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generative Model”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗
DOE OSTI2011

High-resolution structure of a retroviral protease folded as a monomer

The crystal structure of Mason–Pfizer monkey virus protease folded as a monomer has been solved by molecular replacement using a model generated by players of the online game Foldit. The structure shows at high resolution the details of a retroviral protease folded as a monomer which can guide rational design of protease dimerization inhibitors as retroviral drugs. Mason–Pfizer monkey virus (M-PMV), a D-type retrovirus assembling in the cytoplasm, causes simian acquired immunodeficiency syndrome (SAIDS) in rhesus monkeys. Its pepsin-like aspartic protease (retropepsin) is an integral part of the expressed retroviral polyproteins. As in all retroviral life cycles, release and dimerization of the protease (PR) is strictly required for polyprotein processing and virion maturation. Biophysical and NMR studies have indicated that in the absence of substrates or inhibitors M-PMV PR should fold into a stable monomer, but the crystal structure of this protein could not be solved by molecular replacement despite countless attempts. Ultimately, a solution was obtained in mr-rosetta using a model constructed by players of the online protein-folding game Foldit. The structure indeed shows a monomeric protein, with the N- and C-termini completely disordered. On the other hand, the flap loop, which normally gates access to the active site of homodimeric retropepsins, is clearly traceable in the electron density. The flap has an unusual curled shape and a different orientation from both the open and closed states known from dimeric retropepsins. The overall fold of the protein follows the retropepsin canon, but the C{sup α} deviations are large and the active-site ‘DTG’ loop (here NTG) deviates up to 2.7 Å from the standard conformation. This structure of a monomeric retropepsin determined at high resolution (1.6 Å) provides important extra information for the design of dimerization inhibitors that might be developed as drugs for the treatment of retroviral infections, including AIDS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Assessment of Industrial Load for Demand Response across U.S. Regions of the Western Interconnect

Demand response (DR) has the ability to both increase power grid reliability and potentially reduce operating system costs. Understanding the role of demand response in grid modeling has been difficult due to complex nature of the load characteristics compared to the modeled generation and the variation in load types. This is particularly true of industrial loads, where hundreds of different industries exist with varying availability for demand response. We present a framework considering industrial loads for the development of availability profiles that can provide more regional understanding and can be inserted into analysis software for further study. The developed framework utilizes a number of different informational resources, algorithms, and real-world measurements to perform a bottom-up approach in the development of a new database with representation of the potential demand response resource in the industrial sector across the U.S. This tool houses statistical values of energy and demand response (DR) potential by industrial plant and geospatially locates the information for aggregation for different territories without proprietary information. This report will discuss this framework and the analyzed quantities of demand response for Western Interconnect (WI) in support of evaluation of the cost production modeling with power grid modeling efforts of demand response.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2011

Investigations on Repository Near-Field Thermal Modeling - Repository Science/Thermal Load Management & Design Concepts (M41UF033302)

The various layers of material from the waste package (such as components of the engineered barrier system and the host rock surface) to a given distance within the rock wall at a given distance can be described as concentric circles with varying thermal properties (see Figure 5.1-1). The selected model approach examines the contributions of the waste package, axial waste package neighbors and lateral neighboring emplacement drifts (see Section 5.2.1 and Appendix H, Section 2). In clay and deep borehole media, the peak temperature is driven by the central waste package whereas, in granite and salt, the contribution to the temperature rise by adjacent (lateral) waste packages in drift or emplacement borehole lines is dominant at the time of the peak temperature. Mathematical models generated using Mathcad software provide insight into the effects of changing waste package spacing for six waste forms, namely UOX, MOX, co-extraction, new extraction, E-Chem ceramic and E-Chem metal in four different geologic media (granite, clay, salt and deep borehole). Each scenario includes thermal conductivity and diffusivity for each layer between the waste package and the host rock, dimensions of representative repository designs (such as waste package spacing, drift or emplacement borehole spacing, waste package dimensions and layer thickness), and decay heat curves generated from knowledge of the contents of a given waste form after 10, 50, 100 and 200 years of surface storage. Key results generated for each scenario include rock temperature at a given time calculated at a given radius from the central waste package (Section 5.2.1 and Appendix H, Section 3), the corresponding temperature at the interface of the waste package and EBS material, and at each EBS layer in between (Section 5.2.2 and Appendix H, Section 4). This information is vital to understand the implications of repository design (waste package capacity, surface storage time, waste package spacing, and emplacement drift or borehole spacing) by comparing the peak temperature to the thermal limits of the concentric layers surrounding the waste package; specifically 100 C for the bentonite buffer in granite and clay repositories, 100 C for rock wall in a clay repository and 200 C at the rock wall for a salt repository. These thermal limits are both preliminary and approximate, and serve as a means to evaluate design options rather than determining compliance for licensing situations. The thermal behavior of a salt repository is more difficult to model because it is not a concentric geometry and because the crushed salt backfill initially has a much higher thermal resistance than intact salt. Three models were investigated, namely a waste package in complete contact with crushed salt, secondly a waste package in contact with intact salt, and thirdly a waste package in contact with 75% intact and 25% crushed salt. The latter model best depicts emplacement of a waste package in the corner of an intact salt alcove and subsequently covered with crushed salt backfill to the angle of repose. The most conservative model (crushed salt) had temperatures much higher than the other models and although bounding, is too conservative to use. The most realistic model (75/25) had only a small temperature difference from the simplest (non-conservative, intact salt) model, and is the one chosen in this report (see Section 5.2.3). A trade-study investigating three key variables (surface storage time, waste package capacity and waste package spacing) is important to understand and design a repository. Waste package heat can be reduced by storing for longer periods prior to emplacement, or by reducing the number of assemblies or canisters within that waste package. Waste package spacing can be altered to optimize the thermal load without exceeding the thermal limits of the host rock or EBS components. By examining each of these variables, repository footprint (and therefore cost) can be optimized. For this report, the layout was fixed for each geologic medium based on prior published designs in the international community, but it will be varied in future work. Section 5.2.4 summarizes the conclusions based on varying two of the three parameters (storage time and waste package capacity), and the results are shown in Appendix H, Sections 5 and 6.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

ACCELERATED CREEP TEST (ACT) QUALIFICATION OF CREEP RESISTANCE USING THE WCS CONSTITUTIVE MODEL AND STEPPED ISOSTRESS METHOD (SSM)

In this study, a qualification of accelerated creep-resistance of Inconel 718 is assessed using the novel Wilshire-Cano-Stewart (WCS) model and the stepped isostress method (SSM) and predictions are made to conventional creep data. Conventional creep testing (CCT) is a long-term continuous process, in fact, the ASME B&PV III requires that 10,000+ hours of experiments must be conducted to each heat for materials employed in boilers and/or pressure vessel components. This process is costly and not feasible for rapid development of new materials. As an alternative, accelerated creep testing techniques have been developed to reduce the time needed to characterize the creep resistance of materials. Most techniques are based upon the time-temperature-stress superposition principle (TTSSP) that predicts minimum-creep-strain-rate (MCSR) and stress-rupture behaviors but lack the ability to predict creep deformation and consider deformation mechanisms that occur for experiments of longer duration. The stepped isostress method (SSM) has been developed which enables the prediction of creep deformation response as well as reduce the time needed for qualification of materials. The SSM approach has been successful for polymer, polymeric composites, and recently has been introduced for metals. In this study, the WCS constitutive model, calibrated to SSM test data, qualifies the creep resistance of Inconel 718 at 750°C and predictions are compared to CCT data. The WCS model has proven to make long-term predictions for stress-rupture, minimum-creep-strain-rate (MCSR), creep deformation, and damage in metallic materials. The SSM varies stress levels after time interval adding damage to the material, which can be tracked by the WCS model. The SSM data is calibrated into the model and the WCS model generates realistic predictions of stress-rupture, MSCR, damage, and creep deformation. The calibrated material constants are used to generate predictions of stress-rupture and are post-audit validated using the National Institute of Material Science (NIMS) database. Similarly, the MCSR predictions are compared from previous studies. Finally the creep deformation predictions are compared with real data and is determined that the results are well in between the expected boundaries. Material characterization and mechanical properties can be determined at a faster rate and with a more cost-effective method. This is beneficial for multiple applications such as in additive manufacturing, composites, spacecraft, and Industrial Gas Turbines (IGT).

36 MATERIALS SCIENCE↗
DOE OSTI2018

Development of Application Function Blocks for Power-Hardware-in-the-Loop Testing of Grid-Connected Inverters: Preprint

Stability and accuracy of the PHIL platform has been a common problem and challenging issue for a long time, thus continuous efforts and research have been devoted to developing interface algorithms to improve upon this. Some of existing works exhibit promising solutions to compensate the errors in the PHIL platform. However, there is a need to harmonize various solutions into a standard framework that provides a general solution for various PHIL test beds and that can be reused with minimum modifications to avoid extra effort. A good way to do this is by using application function blocks (AFBs), which are model-based and can be developed as a standard library component, thus achieving satisfactory reusability and reconfigurability for various platforms. Since PHIL simulation has the flexibility of the software simulation, proper functions can be easily implemented in the digital real-time simulator (DRTS) to preprocess the reference signal before sending it to the power amplifier. Therefore, this paper aims to develop standard AFBs in DRTS for the PHIL interface to test grid-connected inverters. The PHIL interface algorithm is separated into five AFBs that will be presented in detail in this paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Modelling infrasound signal generation from two underground explosions at the Source Physics Experiment using the Rayleigh integral

For this study, we use the Rayleigh integral (RI) as an approximation to the Helmholtz–Kirchoff integral to model infrasound generation and propagation from underground chemical explosions at distances of 250 m out to 5 km as part of the Source Physics Experiment (SPE). Using a sparse network of surface accelerometers installed above ground zero, we are able to accurately create synthetic acoustic waveforms and compare them to the observed data. Although the underground explosive sources were designed to be symmetric, the resulting seismic wave at the surface shows an asymmetric propagation pattern that is stronger to the northeast of the borehole. This asymmetric bias may be attributed to the subsurface geology and faulting of the area and is observed in the acoustic waveforms. We compare observed and modelled results from two of the underground SPE tests with a sensitivity study to evaluate the asymmetry observed in the data. This work shows that it is possible to model infrasound signals from underground explosive sources using the RI and that asymmetries observed in the data can be modelled with this technique.

58 GEOSCIENCES↗
DOE OSTI2013

Distance-redshift relations in an anisotropic cosmological model

In this paper we study an anisotropic model generated from a particular Bianchi type-III metric, which is a generalization of Gödel's metric and an exact solution of Einstein's field equations. We analyse type Ia supernova data, namely the SDSS sample calibrated with the MLCS2k2 fitter, and we verify in which ranges of distances and redshifts the anisotropy could be observed. We also consider, in a joint analysis, the position of the first peak in the CMB anisotropy spectrum, as well as current observational constraints on the Hubble constant. We conclude that a small anisotropy is permitted by the data, and that more accurate measurements of supernova distances above z = 2 might indicate the existence of such anisotropy in the universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Allometric Scaling of Hyporheic Respiration Across Basins in the Pacific Northwest United States

Abstract Hyporheic zones regulate biogeochemical processes in streams and rivers, but high spatiotemporal heterogeneity makes it difficult to predict how these processes scale from individual reaches to river basins. Recent work applying allometric scaling (i.e., power‐law relationships between size and function) to river networks provides a new paradigm for understanding cumulative hyporheic biogeochemical processes. We used previously published model predictions of reach‐scale hyporheic aerobic respiration to explore patterns in allometric scaling across two climatically divergent basins with differing characteristics in the Pacific Northwest, United States. In the model, hydrologic exchange fluxes (HEFs) regulate hyporheic respiration, so we examined how HEFs might influence allometric scaling of respiration. We found consistent scaling behaviors where HEFs were either very low or very high, but differences between basins when HEFs were moderate. Our findings provide initial model‐generated hypotheses for factors influencing allometric scaling of hyporheic respiration. These hypotheses can be used to optimize new data generation efforts aimed at developing predictive understanding of allometries that can, in turn, be used to scale biogeochemical dynamics across watersheds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Avoidable errors in deposited macromolecular structures: an impediment to efficient data mining

Whereas the vast majority of the more than 85 000 crystal structures of macromolecules currently deposited in the Protein Data Bank are of high quality, some suffer from a variety of imperfections. Although this fact has been pointed out in the past, it is still worth periodic updates so that the metadata obtained by global analysis of the available crystal structures, as well as the utilization of the individual structures for tasks such as drug design, should be based on only the most reliable data. Here, selected abnormal deposited structures have been analysed based on the Bayesian reasoning that the correctness of a model must be judged against both the primary evidence as well as prior knowledge. These structures, as well as information gained from the corresponding publications (if available), have emphasized some of the most prevalent types of common problems. The errors are often perfect illustrations of the nature of human cognition, which is frequently influenced by preconceptions that may lead to fanciful results in the absence of proper validation. Common errors can be traced to negligence and a lack of rigorous verification of the models against electron density, creation of non-parsimonious models, generation of improbable numbers, application of incorrect symmetry, illogical presentation of the results, or violation of the rules of chemistry and physics. Paying more attention to such problems, not only in the final validation stages but during the structure-determination process as well, is necessary not only in order to maintain the highest possible quality of the structural repositories and databases but most of all to provide a solid basis for subsequent studies, including large-scale data-mining projects. For many scientists PDB deposition is a rather infrequent event, so the need for proper training and supervision is emphasized, as well as the need for constant alertness of reason and critical judgment as absolutely necessary safeguarding measures against such problems. Ways of identifying more problematic structures are suggested so that their users may be properly alerted to their possible shortcomings.

model validation↗
DOE OSTI2014

Modeling the transmission of community-associated methicillin-resistant Staphylococcus aureus: a dynamic agent-based simulation

Background: Methicillin-resistant Staphylococcus aureus (MRSA) has been a deadly pathogen in healthcare settings since the 1960s, but MRSA epidemiology changed since 1990 with new genetically distinct strain types circulating among previously healthy people outside healthcare settings. Community-associated (CA) MRSA strains primarily cause skin and soft tissue infections, but may also cause life-threatening invasive infections. First seen in Australia and the U.S., it is a growing problem around the world. The U.S. has had the most widespread CA-MRSA epidemic, with strain type USA300 causing the great majority of infections. Individuals with either asymptomatic colonization or infection may transmit CA-MRSA to others, largely by skin-to-skin contact. Control measures have focused on hospital transmission. Limited public health education has focused on care for skin infections. Methods: We developed a fine-grained agent-based model for Chicago to identify where to target interventions to reduce CA-MRSA transmission. An agent-based model allows us to represent heterogeneity in population behavior, locations and contact patterns that are highly relevant for CA-MRSA transmission and control. Drawing on nationally representative survey data, the model represents variation in sociodemographics, locations, behaviors, and physical contact patterns. Transmission probabilities are based on a comprehensive literature review. Results: Over multiple 10-year runs with one-hour ticks, our model generates temporal and geographic trends in CA-MRSA incidence similar to Chicago from 2001 to 2010. On average, a majority of transmission events occurred in households, and colonized rather than infected agents were the source of the great majority (over 95%) of transmission events. The key findings are that infected people are not the primary source of spread. Rather, the far greater number of colonized individuals must be targeted to reduce transmission. Conclusions: Our findings suggest that current paradigms in MRSA control in the United States cannot be very effective in reducing the incidence of CA-MRSA infections. Furthermore, the control measures that have focused on hospitals are unlikely to have much population-wide impact on CA-MRSA rates. New strategies need to be developed, as the incidence of CA-MRSA is likely to continue to grow around the world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Dynamic Simulation Technoeconomic Model for Power Generation

Sandia National Laboratories has built and successfully tested a dynamic simulation technoeconomic model of the Palo Verde Generating Station that is now being updated to help other US power plants improve operations. Palo Verde, located west of Phoenix, Arizona, is the largest electricity generator in the US at 4 GW. Palo Verde uses — 60 million gallons per day of treated wastewater from Phoenix to cool reactors, and disposes of blowdown in evaporation ponds. The model built for Palo Verde numerically evaluates the economic impact of changing, for example, alternative cooling technologies, water usage and treatment, and influent water chemistry, and is based on detailed accounting of mass, energy, and cash flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Sources and composition of submicron organic mass in marine aerosol particles

Recent studies have proposed a variety of interpretations of the sources and composition of atmospheric marine aerosol particles (aMA) based on a range of physical and chemical measurements collected during open-ocean research cruises. To investigate the processes that affect marine organic particles, this study uses the characteristic functional group composition (from Fourier transform infrared (FTIR) spectroscopy) of aMAP from five ocean regions to show that: (i) The organic functional group composition of aMAP that can be identified as atmospheric primary marine (ocean-derived) aerosol (aPMA) is 65±12% hydroxyl, 21±9% alkane, 6±6% amine, and 7±8% carboxylic acid functional groups. Contributions from photochemical reactions add carboxylic acid groups (15%-25%), shipping effluent in seawater and ship emissions add additional alkane groups (up to 70%), and coastal emissions mix in alkane and carboxylic acid groups from coastal pollution sources. (ii) The organic composition of aPMA is nearly identical to model generated primary marine aerosol particles (gPMA) from bubbled seawater (55% hydroxyl, 32% alkane, and 13% amine functional groups), indicating that its overall functional group composition is the direct consequence of the organic constituents of the seawater source. (iii) While the seawater organic functional group composition was nearly invariant across all three ocean regions studied, the gPMA alkane group fraction increased with chlorophyll-a concentrations (r = 0.79). gPMA from productive seawater had a larger fraction of alkane functional groups (35%) compared to gPMA from non-productive seawater (16%), likely due to the presence of surfactants in productive seawater that stabilize the bubble film and lead to preferential drainage of the more soluble (lower alkane group fraction) organic components. gPMA has a hydroxyl group absorption peak location characteristic of monosaccharides and disaccharides, where the seawater OM hydroxyl group peak location is closer to that of polysaccharides. This may result from the larger saccharides preferentially remaining in the seawater during gPMA and aPMA production

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Analysis of albedo versus cloud fraction relationships in liquid water clouds using heuristic models and large eddy simulation: ALBEDO VERSUS CLOUD FRACTION RELATIONSHIPS

The relationship between the albedo of a cloudy scene $A$ and cloud fraction f c is studied with the aid of heuristic models of stratocumulus and cumulus clouds. Existing work has shown that scene albedo increases monotonically with increasing cloud fraction but that the relationship varies from linear to superlinear. The reasons for these differences in functional dependence are traced to the relationship between cloud deepening and cloud widening. When clouds deepen with no significant increase in f c (e.g., in solid stratocumulus), the relationship between$A$ and f c is linear.When clouds widen as they deepen, as in cumulus cloud fields, the relationship is superlinear. A simple heuristic model of a cumulus cloud field with a power law size distribution shows that the superlinear $A$-f c behavior is traced out either through random variation in cloud size distribution parameters or as the cloud field oscillates between a relative abundance of small clouds (steep slopes on a log-log plot) and a relative abundance of large clouds (flat slopes). Oscillations of this kind manifest in large eddy simulation of trade wind cumulus where the slope and intercept of the power law fit to the cloud size distribution are highly correlated. Further analysis of the large eddy model-generated cloud fields suggests that cumulus clouds grow larger and deeper as their underlying plumes aggregate; this is followed by breakup of large plumes and a tendency to smaller clouds. The cloud and thermal size distributions oscillate back and forth approximately in unison.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Phenome‐to‐genome insights for evaluating root system architecture in field studies of maize

Abstract Understanding the genetic basis of root system architecture (RSA) in crops requires innovative approaches that enable both high‐throughput and precise phenotyping in field conditions. In this study, we evaluated multiple phenotyping and analytical frameworks for quantifying RSA in mature, field‐grown maize in three field experiments. We used forward and reverse genetic approaches to evaluate >1700 maize root crowns, including a diversity panel, a biparental mapping population, and maize mutant and wild‐type alleles at two known RSA genes,DEEPER ROOTING 1(DRO1) andRootless1(Rt1). We show the utility of increasing the dimensionality of traditional two‐dimensional (2D) techniques, referred to as the “2D multi‐view” method, to improve the capture of whole root system information for mapping genetic variation influencing RSA. Comparison of univariate and multivariate genome‐wide association study (GWAS) approaches revealed that multivariate traits were effective at dissecting complex RSA phenotypes and identifying pleiotropic quantitative trait loci (QTLs). Overall, three‐dimensional (3D) root models generated from X‐ray computed tomography and digital phenotyping captured a larger proportion of RSA trait variations compared to other methods of root phenotyping, as evidenced by both genome‐wide and single‐gene analyses. Among the individual root traits, root pulling force emerged as a highly heritable estimate of RSA that identified the largest number of shared QTLs with 3D phenotypes. Our study shows that integrating complementary phenotyping technologies helps to provide a more comprehensive understanding of the genetic architecture of RSA in field‐grown maize.

Genetics & Heredity↗
DOE OSTI2017

Measurement of charged-particle distributions sensitive to the underlying event in $ \sqrt{s}=13 $ TeV proton-proton collisions with the ATLAS detector at the LHC

We present charged-particle distributions sensitive to the underlying event, measured by the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV, in low-luminosity Large Hadron Collider fills corresponding to an integrated luminosity of 1.6 nb –1 . The distributions were constructed using charged particles with absolute pseudorapidity less than 2.5 and with transverse momentum greater than 500 MeV, in events with at least one such charged particle with transverse momentum above 1 GeV. These distributions characterise the angular distribution of energy and particle flows with respect to the charged particle with highest transverse momentum, as a function of both that momentum and of charged-particle multiplicity. The results have been corrected for detector effects and are compared to the predictions of various Monte Carlo event generators, experimentally establishing the level of underlying-event activity at LHC Run 2 energies and providing inputs for the development of event generator modelling. The current models in use for UE modelling typically describe this data to 5% accuracy, compared with data uncertainties of less than 1%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Deformable Registration for Longitudinal Breast MRI Screening

MRI screening of high-risk patients for breast cancer provides very high sensitivity, but with a high recall rate and negative biopsies. Comparing the current exam to prior exams reduces the number of follow-up procedures requested by radiologists. Such comparison, however, can be challenging due to the highly deformable nature of breast tissues. Automated co-registration of multiple scans has the potential to aid diagnosis by providing 3D images for side-by-side comparison and also for use in CAD systems. Although many deformable registration techniques exist, they generally have a large number of parameters that need to be optimized and validated for each new application. Here, we propose a framework for such optimization and also identify the optimal input parameter set for registration of 3D T{sub 1}-weighted MRI of breast using Elastix, a widely used and freely available registration tool. A numerical simulation study was first conducted to model the breast tissue and its deformation through finite element (FE) modeling. This model generated the ground truth for evaluating the registration accuracy by providing the deformation of each voxel in the breast volume. An exhaustive search was performed over various values of 7 registration parameters (4050 different combinations of parameters were assessed) and the optimum parameter set was determined. This study showed that there was a large variation in the registration accuracy of different parameter sets ranging from 0.29 mm to 2.50 mm in median registration error and 3.71 mm to 8.90 mm in 95 percentile of the registration error. Mean registration errors of 0.32 mm, 0.29 mm, and 0.30 mm and 95 percentile errors of 3.71 mm, 5.02 mm, and 4.70 mm were obtained by the three best parameter sets. The optimal parameter set was applied to consecutive breast MRI scans of 13 patients. A radiologist identified 113 landmark pairs (~ 11 per patient) which were used to assess registration accuracy. The results demonstrated that using the optimal registration parameter set, a registration accuracy (in mm) of 3.4 [1.8 6.8] was achieved.

62 RADIOLOGY AND NUCLEAR MEDICINE↗