Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Science”

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 541 records · Page 30

Search for Dark Matter with the LZ Experiment

Overwhelming evidence from astrophysical observations indicates that only 20% of the matter in the universe consists of ordinary matter, while the remaining 80% is made of some form of matter called "dark matter''. Although these indirect astrophysical observations convince us that dark matter exists, dark matter has not yet been directly observed. The standard model of particle physics, which has been very successful in explaining the properties of ordinary matter, can neither explain dark matter's presence nor its properties. The discovery and identification of dark matter would have a profound impact on cosmology, astronomy, and particle physics. A leading dark matter candidate consistent with all astrophysical data is a weakly interacting massive particle (WIMP). WIMPs could be detected through their interactions with atomic nuclei in specially designed detectors. Since the collision rates are expected to be very small (a few collisions per year for a 1-ton detector), large detectors with low backgrounds and excellent detection capability for rare collisions are needed to gain knowledge of the nature of dark matter. With dark matter experiments such as the LZ experiment, we are in a position to detect and study this dark matter through sensitive terrestrial direct detection experiments. If no discovery is made with LZ, it will be very difficult to reconcile the sensitivity limit of this experiment with most conventional WIMP scenarios, and alternate paradigms will have to be strongly considered for dark matter. The main focus of the U-M dark matter group over the past two years was on design validation and fabrication for the LZ project, followed by commissioning and operation of the LZ dark matter experiment, and the analysis of the science data during the following funding cycle. The U-M dark matter group had identified two key areas to contribute to the LZ project: the design and fabrication of a radon removal element for the xenon circulation system, and the execution and coordination of prototyping and design validation performed by the small xenon chamber group.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

ESGF Future Architecture Report (V.1.1)

This document reports the findings of the Earth System Grid Federation Future Architecture Face-to-Face Meeting held in the UK 5-7 November 2019, attended by technical representatives from the ESGF partner organisations. ESGF is an International federated archive that provides access to Earth science data. It is fundamentally a distributed file management system with its main function being data dissemination through the provision of services for efficient replication and search. It has been in operation for over ten years and has undergone continuous improvements and enhancements over this time span. The Future Architecture is a major review of the software system including a fundamental reassessment of its requirements with a view to develop a new overall system architecture. This has been initiated in the light of challenges around the maintenance of the existing software legacy and a changing technology landscape and practices in the communities it serves. The work is being carried out under the oversight of the ESGF Executive Committee1,2 and is coordinated activity between ESGF and the European IS-ENES3 project.

54 ENVIRONMENTAL SCIENCES↗

Leveraging graph clustering techniques for cyber‐physical system analysis to enhance disturbance characterisation

Abstract Cyber‐physical systems have behaviour that crosses domain boundaries during events such as planned operational changes and malicious disturbances. Traditionally, the cyber and physical systems are monitored separately and use very different toolsets and analysis paradigms. The security and privacy of these cyber‐physical systems requires improved understanding of the combined cyber‐physical system behaviour and methods for holistic analysis. Therefore, the authors propose leveraging clustering techniques on cyber‐physical data from smart grid systems to analyse differences and similarities in behaviour during cyber‐, physical‐, and cyber‐physical disturbances. Since clustering methods are commonly used in data science to examine statistical similarities in order to sort large datasets, these algorithms can assist in identifying useful relationships in cyber‐physical systems. Through this analysis, deeper insights can be shared with decision‐makers on what cyber and physical components are strongly or weakly linked, what cyber‐physical pathways are most traversed, and the criticality of certain cyber‐physical nodes or edges. This paper presents several types of clustering methods for cyber‐physical graphs of smart grid systems and their application in assessing different types of disturbances for informing cyber‐physical situational awareness. The collection of these clustering techniques provide a foundational basis for cyber‐physical graph interdependency analysis.

97 MATHEMATICS AND COMPUTING↗

Foundational Science to Accelerate Nuclear Energy Innovation [Brochure]

The foundational science gaps inhibiting the advancement of nuclear energy technologies are identified and tackled in five priority research opportunities. These opportunities pave the way to accelerate the development and ultimately the adoption of new nuclear energy systems. They include the fundamental aspects of ion-electron interactions; novel properties of next-generation coolants and solvents; interfacial dynamics, not only in solids, but in other aspects of nuclear reactors; novel operando and in-situ monitoring and sensing; and artificial intelligence to accelerate condensed phases discovery. Building on the foundation established by previous BES workshops, these opportunities encompass recent advances in fundamental knowledge and focus on the experimental and computational methods needed to resolve major technical challenges for nuclear energy technologies. Through developing fundamental scientific insight as well as pushing the frontiers of modeling complex systems and probing the operation of materials and chemical systems in extreme environments, research motivated by the priorities identified here will further develop the promise, potential, and utilization of nuclear energy for a clean energy future. The PROs are as follows: (1) Master complex electronic structures to tailor thermochemical reactivity, transport, and microstructural evolution; (2) Interrogate and direct the physics and chemistry underpinning next-generation coolants and solvents; (3) Elucidate and control the underlying physics and chemistry of interfaces in complex nuclear environments; (4) Bridge multi-fidelity multi-resolution experiments, computational modeling, and data science to control dynamic behavior; and (5) Harness artificial intelligence to design inherently resilient condensed phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Downscaling atmospheric chemistry simulations with physically consistent deep learning

Abstract. Recent advances in deep convolutional neural network (CNN)-based super resolution can be used to downscale atmospheric chemistry simulations with substantially higher accuracy than conventional downscaling methods. This work both demonstrates the downscaling capabilities of modern CNN-based single image super resolution and video super-resolution schemes and develops modifications to these schemes to ensure they are appropriate for use with physical science data. The CNN-based video super-resolution schemes in particular incur only 39 % to 54 % of the grid-cell-level error of interpolation schemes and generate outputs with extremely realistic small-scale variability based on multiple perceptual quality metrics while performing a large (8×10) increase in resolution in the spatial dimensions. Methods are introduced to strictly enforce physical conservation laws within CNNs, perform large and asymmetric resolution changes between common model grid resolutions, account for non-uniform grid-cell areas, super-resolve lognormally distributed datasets, and leverage additional inputs such as high-resolution climatologies and model state variables. High-resolution chemistry simulations are critical for modeling regional air quality and for understanding future climate, and CNN-based downscaling has the potential to generate these high-resolution simulations and ensembles at a fraction of the computational cost.

54 ENVIRONMENTAL SCIENCES↗

Task I: Dark Matter Search Experiments with Cryogenic Detectors: CDMS-I and CDMS-II Task II: Experimental Study of Neutrino Properties: EXO and KamLAND

Dark Matter Search - During the period of performance, our group continued the search for dark matter in the form of weakly interacting massive particles or WIMPs. As a key member of the CDMS (Cryogenic Dark Matter Search) collaboration, we completed the CDMS II experiment which led the field in sensitivity for more than five years. We fabricated all detectors, and participated in detector testing and verification. In addition, we participated in the construction and operation of the facility at the Soudan Underground Laboratory and played key roles in the data acquisition and analysis. Towards the end of the performance period, we began operating the SuperCDMS Soudan experiment, which consists of 15 advanced Ge (9 kg) detectors. The advanced detector design called iZIP grew out of our earlier DOE Particle Detector R&D program which demonstrated the rejection of surface electrons to levels where they are no longer the dominant source of background. Our group invented this advanced design and these larger detectors were fabricated on the Stanford campus in collaboration with the SLAC CDMS group and the Santa Clara University group. The sensitivity reach is expected to be up to 5 times better than CDMS II after two years of operation. We will check the new limits on WIMPs set by XENON100, and we expect improved sensitivity for light mass WIMPs beyond that of any other existing experiment. Our group includes the Spokesperson for SuperCDMS and continues to make important contributions to improvements in the detector technology which are enabling the very low trigger thresholds used to explore the low mass WIMP region. We are making detailed measurements of the charge transport and trapping within Ge crystals, measuring the diffusive trapping distance of the quasiparticle excitations within the Al phonon collector fins on the detector surface, and we are contributing to the development of much improved detector Monte Carlos which are essential to guide the detector design and optimize the analysis. Neutrino Physics – In the period of performance the neutrino group successfully completed the construction of EXO-200 and commissioned the detector. Science data taking started on Jun 1, 2011. With the discovery of the 2-neutrino double-beta decay in 136-Xe and the first measurement of the 0-neutrino mode resulting in the most stringent limit of Majorana masses, our group continues to be a leading innovator in the field of neutrino physics which is central to DOE-HEP Intensity Frontier program. The phenomenon of neutrino oscillations, in part elucidated by our earlier efforts with the Palo Verde and KamLAND experiments, provides the crucial information that neutrino masses are non-zero and, yet, it contains no information on the value of the neutrino mass scale. In recent times our group has therefore shifted its focus to a high sensitivity 0-neutrino double beta decay program, EXO. The 0-neutrino double beta decay provides the best chance of extending the sensitivity to the neutrino mass scale below 10 meV but, maybe more importantly, it tests the nature of the neutrino wave function, providing the most sensitive probe for Majorana particles and lepton number violation. The EXO program, formulated by our group several years ago, plans to use up to tonnes of the isotope 136-Xe to study the 0-neutrino double beta decay mode. The EXO-200 detector is the first step in this program and it represents the only large US-led and based experiment taking data. The EXO-200 isotope enrichment program broke new grounds for the enterprise of double beta decay. The detector design and material selection program paid off, resulting in a background that is among the very best in the field. The “first light" of EXO-200 was very exciting with the discovery -in the first month of data- of the rarest 2-neutrino double beta decay mode ever observed. The lower limit on the 0-neutrino double beta decay half-life, published in Phys. Rev. Lett. and based on the first 120 days of data is the second best but, when translated into a Majorana mass scale, it is one of the most stringent constraint we have on neutrino masses. Indeed, such a limit was the first result to contradict a claim of discovery in 76-Ge for most nuclear matrix elements calculations. As we continue data taking and plan some modest upgrades to EXO-200 our group is also very active in the design of nEXO, a 5 tonne detector based on the technology demonstrated by EXO-200. Over the years we have made it a tradition to explore the frontier and not to be shy about looking in new directions and re-inventing ourselves to best take advantage of the precious few opportunities provided by Nature. We have also cultivated a number of young people at all levels and, by now, many of the undergraduates, graduate students and postdocs educated by this group have leading positions in academia and industry around the world.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rethinking Visual Analytics for Streaming Data Applications

In the age of data science, the use of interactive information visualization techniques has become increasingly ubiquitous. From online scientific journals to the New York Times graphics desk, the utility of interactive visualization for both storytelling and analysis has become ever more apparent. As these techniques have become more readily accessible, the appeal of combining interactive visualization with computational analysis continues to grow. Arising out of a need for scalable, human-driven analysis, primary objective of visual analytics systems is to capitalize on the complementary strengths of human and machine analysis, using interactive visualization as a medium for communication between the two. These systems leverage developments from the fields of information visualization, computer graphics, machine learning, and human-computer interaction to support insight generation in areas where purely computational analyses fall short. Over the past decade, visual analytics systems have generated remarkable advances in many historically challenging analytical contexts. These include areas such as modeling political systems [Crouser et al. 2012], detecting financial fraud [Chang et al. 2008], and cybersecurity [Harrison et al. 2012]. In each of these contexts, domain expertise and human intuition is a necessary component of the analysis. This intuition is essential to building trust in the analytical products, as well as supporting the translation of evidence into actionable insight. In addition, each of these examples also highlights the need for scalable analysis. In each case, it is infeasible for a human analyst to manually assess the raw information unaided, and the communication overhead to divide the task between a large number of analysts makes simple parallelism intractable. Regardless of the domain, visual analytics tools strive to optimize the allocation of human analytical resources, and to streamline the sensemaking process on data that is massive, complex, incomplete, and uncertain in scenarios requiring human judgment.

97 MATHEMATICS AND COMPUTING↗

Optimal Pulse Processing, Pile-Up Decomposition, and Applications of Silicon Drift Detectors at LCLS

Silicon drift detectors (SDDs) revolutionized spectroscopy in fields as diverse as geology and dentistry. For a subset of experiments at ultrafast, X-ray free-electron lasers (FELs), SDDs can make substantial contributions. Often the unknown spectrum is interesting, carrying science data, or the background measurement is useful to identify unexpected signals. Many measurements involve only several discrete photon energies known a priori, allowing single-event decomposition of pile-up and spectroscopic photon counting. We designed a pulse function and demonstrated that the signal amplitude (i.e., proportional to the detected energy and obtained from fitting with the pulse function), rise time, and pulse height are interrelated, and at short peaking times, the pulse height and pulse area are not optimal estimators for detected energy; instead, the signal amplitude and rise time are obtained for each pulse by fitting, thus removing the need for pulse shaping. By avoiding pulse shaping, rise times of tens of nanoseconds resulted in reduced pulse pile-up and allowed decomposition of remaining pulse pile-up at photon separation times down to hundreds of nanoseconds while yielding time-of-arrival information with the precision of 10 ns. Waveform fitting yields simultaneously high energy resolution and high counting rates (two orders of magnitude higher than current digital pulse processors). At pulsed sources or high photon rates, photon pile-up still occurs. We showed that pile-up spectrum fitting is relatively simple and preferable to pile-up spectrum deconvolution. We then developed a photon pile-up statistical model for constant intensity sources, extended it to variable intensity sources (typical for FELs), and used it to fit a complex pileup spectrum. We subsequently developed a Bayesian pile-up decomposition method that allows decomposing pile-up of single events with up to six photons from six monochromatic lines with 99% accuracy. The usefulness of SDDs will continue into the X-ray FEL era of science. Their successors, the ePixS hybrid pixel detectors, already offer hundreds of pixels, each with a similar performance to an SDD, in a compact, robust and affordable package.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

GLM Observations of the Brightest Lightning in the Americas

Two years of Geostationary Lightning Mapper (GLM) science data are used to document the brightest lighting flashes observed on the Americas continent. The most radiant optical lightning emissions—termed “superbolts”—were first identified by the Vela satellite constellation in the 1970s (Turman, 1977) and are defined in terms of peak optical power. GLM is an integrating sensor that, instead, measures the total optical energy from a lightning pulse. Furthermore, while GLM might not correctly classify short-duration superbolts, its top lightning cases certainly fall in the superbolt category, and the wealth of GLM measurements over its stationary hemispheric field of view provide an unmatched sample of extraordinarily bright lightning. While radiant bolts in excess of 100× the optical energy of typical lightning are ubiquitous across the Americas and result from many types of lightning processes, we find the most radiant cases (>1,000×) are concentrated in the central United States and in the La Plata basin in South America. Coincident Earth Networks Global Lightning Network (ENGLN) observations reveal that these extremely bright emissions usually result from + CG strokes with high peak currents in long horizontal flashes outside of the convective core. Single cases of these megaflashes might produce multiple superbolts over their durations.

54 ENVIRONMENTAL SCIENCES↗

CORPSE model with litter decomposition parameters derived from the LIDET dataset

This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of files:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N files have the same columns (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436

Juice, Stephanie↗

Cosmogenic production of 37 Ar in the context of the LUX-ZEPLIN experiment

We estimate the amount of 37 Ar produced in natural xenon via cosmic-ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth’s surface. We then calculate the resulting 37 Ar concentration in a 10-tonne payload (similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage, and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea-level production rate of 37 Ar in natural xenon is estimated to be 0.024 atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1 tonne/month, the average 37 Ar activity after 10 tons are purified and transported underground is 0.058 − 0.090 𝜇⁢Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic 37 Ar will appear as a noticeable background in the early science data, while decaying with a 35-day half-life. This newly noticed production mechanism of 37 Ar should be considered when planning for future liquid-xenon-based experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Molecular Simulations of the Chain Length Dependent Adsorption of C7‐C14 n ‐Alkanes in ZIF‐8

Recent experiments show that the diffusivities of C8 and C10 n-alkanes in ZIF-8 are higher than those of C7 and C9, respectively. We investigated this unusual ‘odd-even' effect by simulating the adsorption of C7-C14 n-alkanes in ZIF-8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest-host energies, isosteric heats, and chain length distributions are analyzed for trends among the n-alkanes. ZIF-8 cages filled with n-alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. In conclusion, results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Methodology for Good Machine Learning with Multi‐Omics Data

In 2020, Novartis Pharmaceuticals Corporation and the U.S. Food and Drug Administration (FDA) started a 4‐year scientific collaboration to approach complex new data modalities and advanced analytics. The scientific question was to find novel radio‐genomics‐based prognostic and predictive factors for HR+/HER− metastatic breast cancer under a Research Collaboration Agreement. This collaboration has been providing valuable insights to help successfully implement future scientific projects, particularly using artificial intelligence and machine learning. This tutorial aims to provide tangible guidelines for a multi‐omics project that includes multidisciplinary expert teams, spanning across different institutions. We cover key ideas, such as “maintaining effective communication” and “following good data science practices,” followed by the four steps of exploratory projects, namely (1) plan, (2) design, (3) develop, and (4) disseminate. We break each step into smaller concepts with strategies for implementation and provide illustrations from our collaboration to further give the readers actionable guidance.

Pharmacology & Pharmacy↗

Uncertainty in inventories for life cycle assessment: State‐of‐the‐art, challenges, and new technologies

Uncertainty is a critical factor that can hinder the quality and potential applications of life cycle assessment (LCA) results. A prominent source of uncertainty stems from the life cycle inventory (LCI) data. Various methodologies exist to estimate the uncertainty associated with LCI data, primarily based on the widely used structured pedigree matrix approach or the computationally intensive Monte Carlo simulation. This perspective review explores how new technologies (e.g., computational algorithms and data collection methods) from data science and related fields can contribute to identifying, quantifying, and reducing uncertainty in LCI modeling. A brief overview of the sources of uncertainty in LCI modeling and how they are addressed in current LCA practice is provided. Additionally, several new technologies are identified, and the potential benefits of their implementation in reducing uncertainties in LCI modeling are discussed. This perspective review concludes by identifying potential areas that require further development for these technologies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

What is (quantitative) system dynamics modeling? Defining characteristics and the opportunities they create

A clear definition of system dynamics modeling can provide shared understanding and clarify the impact of the field. We introduce a set of characteristics that define quantitative system dynamics, selected to capture core philosophy, describe theoretical and practical principles, and apply to historical work but be flexible enough to remain relevant as the field progresses. The defining characteristics are: (1) models are based on causal feedback structure, (2) accumulations and delays are foundational, (3) models are equation-based, (4) concept of time is continuous, and (5) analysis focuses on feedback dynamics. We discuss the implications of these principles and use them to identify research opportunities in which the system dynamics field can advance. These research opportunities include causality, disaggregation, data science and AI, and contributing to scientific advancement. Progress in these areas has the potential to improve both the science and practice of system dynamics.

97 MATHEMATICS AND COMPUTING↗

Efficient proximal subproblem solvers for a nonsmooth trust-region method

In [R. J. Baraldi and D. P. Kouri, Mathematical Programming, (2022), pp. 1-40], we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex and nonsmooth convex function. The principle expense of this method is in computing a trial iterate that satisfies the so-called fraction of Cauchy decrease condition—a bound that ensures the trial iterate produces sufficient decrease of the subproblem model. In this paper, we expound on various proximal trust-region subproblem solvers that generalize traditional trust-region methods for smooth unconstrained and convex-constrained problems. We introduce a simplified spectral proximal gradient solver, a truncated nonlinear conjugate gradient solver, and a dogleg method. Finally, we compare algorithm performance on examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗