Search NASA⌕ Search

Engineering topics

Lawrence, David

Publications and source records attributed to Lawrence, David.

53 records · Page 3

A.I. for nuclear physics

Nuclear science is concerned with the understanding of the nature of matter, its basic constituents and their interaction to form the elements and the properties we observe. This includes the forms of matter we see around us and also exotic forms such as those that existed in the first moments after the Big Bang and that exist today inside neutron stars. The techniques, tools, and expertise needed for nuclear physics (NP) research are therefore diverse in nature. State-of-the art accelerators are being developed to illuminate the dynamical basis of the core of the atom in terms of the fundamental constituents called quarks and gluons and to increase the number of isotopes with known properties. This scientific infrastructure is reaching scales and complexities that require computational methods for tasks such as anomaly detection in operational data. New methodologies are needed to detect anomalies and to optimize operating parameters, predict failures as well as to discover new optimization algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

JANA2 LDRD (Final Project Report)

The project sought to develop a second generation of the JANA reconstruction framework that leverages modern programming features (C++11 and C++14) while maintaining its basic user centric design. The goals included making the project easy to deploy to large HPC facilities such as NERSC and to include detailed documentation to allow scientists outside of Jefferson Lab to use it at such a facility. This would position JLab to play a key role in the larger Laboratory based computing ecosystem under the ASCR initiative. The experience gained will help inform the design of Exascale class applications for experimental nuclear physics. Developing a next generation JANA provides an advantage over older serial frameworks or ones retro-fitted to provide some multi-threading support. That technique requires significantly more memory than multi-threading which leads to compute jobs being limited by available RAM, and not fully utilizing the available CPU. The smaller memory footprint of multi-threaded programs also leads to increased efficiency due to fewer cache misses. It is for these and other reasons that multi-threading is the more accepted industry standard in the modern era.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Full Implementation of Matrix Approach to Biogeochemistry Module of CLM5

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

58 GEOSCIENCES↗

JANA2 Framework for Event Based and Triggerless Data Processing

Development of the second generation JANA2 multi-threaded event processing framework is ongoing through an LDRD initiative grant at Jefferson Lab. The framework is designed to take full advantage of all cores on modern many-core compute nodes. JANA2 efficiently handles both traditional hardware triggered event data and streaming data in online triggerless environments. Development is being done in conjunction with the Electron Ion Collider development. Anticipated to be the next large scale Nuclear Physics facility constructed. The core framework is written in modern C++ but includes an integrated Python interface. The status of development and summary of the more interesting features are presented.

Lawrence, David↗

ML Track Fitting in Nuclear Physics

Charged particle tracking represents the largest consumer of CPU resources in high data volume Nuclear Physics (NP) experiments. An effort is underway to develop machine learning (ML) networks that will reduce the resources required for charged particle tracking. Tracking in NP experiments represent some unique challenges compared to high energy physics (HEP). In particular, track finding typically represents only a small fraction of the overall tracking problem in NP. This presentation will outline the differences and similarities between NP and HEP charged particle tracking and areas where ML learning may provide a benefit. The status of the specific effort taking place at Jefferson Lab will also be shown.

Britton, Thomas↗

Offsite Data Processing for the GlueX Experiment

The Thomas Jefferson National Accelerator Facility (JLab) 12GeV accelerator upgrade completed in 2015 is now producing data at volumes unprecedented for the lab. The resources required to process this data now exceed the capacity of the onsite farm necessitating the use of offsite computing resources for the first time in the history of JLab. GlueX is now utilizing NERSC and PSC for raw data production. Details of the workflow are presented.

Lawrence, David↗

Cyber-physical Resilient Interoperable Microgrid Networks

The paper presents the results of an industry-academic collaboration between ABB, Duke Energy and the University of Illinois at Urbana-Champaign (UIUC) that has developed and demonstrated a cyberphysical resilient control and communication architecture for deployment of multiple microgrid networks, as well as retrofitting the existing distribution grids by increasing DER penetration levels.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multi-threaded Event Reconstruction with JANA

JANA2 is a multi-threaded event reconstruction framework being developed for Experimental Nuclear Physics. It is an LDRD funded project that will be the successor of the original JANA framework. JANA2 is a near complete rewrite emphasizing C++ language features that have only become available since the C++11 standard. features such as shared pointers, language native threading, and atomics. This paper outlines the status of the project.

Lawrence, David↗

Measuring the Charged Pion Polarizability in the $\gamma\gamma\rightarrow\pi^{+{}}\pi^{-{}}$ Reaction

A new measurement of the charged pion polarizability ????? is approved to run at Jefferson Lab using the GlueX detector in experimental Hall D. The charged pion polarizability ranks among the most important tests of low-energy QCD presently unresolved by experiment. Analogous to precision measurements of ?0??? that test the intrinsic odd-parity (anomalous) sector of QCD, the pion polarizability tests the intrinsic even-parity sector of QCD. The ????+?? cross section will be accessed via the Primakoff mechanism on a nuclear target. The linearly polarized photon source in Hall D will allow the separation of the Primakoff cross section from coherent ?0(770) production. An additional detector will also need to be constructed to help distinguish ?+?? and ?+?? events. An update on the status of the experiment will be given.

Lawrence, David↗

Matrix approach to land carbon cycle modeling: A case study with the Community Land Model

The terrestrial carbon (C) cycle has been commonly represented by a series of C balance equations to track C influxes into and effluxes out of individual pools in earth system models (ESMs). This representation matches our understanding of C cycle processes well but makes it difficult to track model behaviors. It is also computationally expensive, limiting the ability to conduct comprehensive parametric sensitivity analyses. To overcome these challenges, we have developed a matrix approach, which reorganizes the C balance equations in the original ESM into one matrix equation without changing any modeled C cycle processes and mechanisms. We applied the matrix approach to the Community Land Model (CLM4.5) with vertically-resolved biogeochemistry. The matrix equation exactly reproduces litter and soil organic carbon (SOC) dynamics of the standard CLM4.5 across different spatial-temporal scales. The matrix approach enables effective diagnosis of system properties such as C residence time and attribution of global change impacts to relevant processes. We illustrated, for example, the impacts of CO 2 fertilization on litter and SOC dynamics can be easily decomposed into the relative contributions from C input, allocation of external C into different C pools, nitrogen regulation, altered soil environmental conditions, and vertical mixing along the soil profile. In addition, the matrix tool can accelerate model spin-up, permit thorough parametric sensitivity tests, enable pool-based data assimilation, and facilitate tracking and benchmarking of model behaviors. Altogether, the matrix approach can make a broad range of future modeling activities more efficient and effective.

59 BASIC BIOLOGICAL SCIENCES↗

Ion Mobility Spectrometry - High Resolution LTQ-Orbitrap Mass Spectrometry for Analysis of Homemade Explosives

The detailed chemical characterization of homemade explosives (HMEs) and other chemicals that can mimic or mask the presence of explosives is important for understanding and improving the performance of commercial instrumentation used for explosive detection. To that end, an atmospheric-pressure drift tube ion mobility spectrometry (IMS) instrument has been successfully coupled to a commercial tandem mass spectrometry (MS) system. The tandem MS system is comprised of a linear ion trap and a high resolution Orbitrap analyzer. This IMS-MS combination allows extensive characterization of threat chemical compounds, including HMEs, and complex real-world background chemicals that can interfere with detection. Here, the composition of ion species originating from a specific HME, erythritol tetranitrate, has been elucidated using accurate mass measurements, isotopic ratios, and tandem MS. Gated IMS-MS and high-resolution MS have been used to identify minor impurities that can be indicative of the HME source and/or synthesis route. Comparison between data obtained on the IMS/MS system and on commercial stand-alone IMS instruments used as explosive trace detectors (ETDs) has also been performed. Such analysis allows better signature assignments of threat compounds, modified detection algorithms, and improved overall ETD performance. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Terrestrial ecosystem model performance in simulating productivity and its vulnerability to climate change in the northern permafrost region: Modeled Productivity in Permafrost Regions

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 ± 6 g C m -2 yr -1 ), most models produced higher NPP (309 ± 12 g C m -2 yr -1 ) over the permafrost region during 2000–2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982–2009, there was a twofold discrepancy among models (380 to 800 g C m -2 yr -1 ), which mainly resulted from differences in simulated maximum monthly GPP (GPP max ). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C (Vc max_25 ), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO 2 concentration. In conclusion, these results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPP max as well as their sensitivity to climate change.

54 ENVIRONMENTAL SCIENCES↗

Interactions between land use change and carbon cycle feedbacks: Land Use and Carbon Cycle Feedbacks

We explore the role of human land use and land cover change (LULCC) in modifying the terrestrial carbon budget in simulations forced by Representative Concentration Pathway 8.5, extended to year 2300 by using the Community Earth System Model, . Overall, conversion of land (e.g., from forest to croplands via deforestation) results in a model-estimated, cumulative carbon loss of 490 Pg C between 1850 and 2300, larger than the 230 Pg C loss of carbon caused by climate change over this same interval. The LULCC carbon loss is a combination of a direct loss at the time of conversion and an indirect loss from the reduction of potential terrestrial carbon sinks. Approximately 40% of the carbon loss associated with LULCC in the simulations arises from direct human modification of the land surface; the remaining 60% is an indirect consequence of the loss of potential natural carbon sinks. Because of the multicentury carbon cycle legacy of current land use decisions, a globally averaged amplification factor of 2.6 must be applied to 2015 land use carbon losses to adjust for indirect effects. This estimate is 30% higher when considering the carbon cycle evolution after 2100. Most of the terrestrial uptake of anthropogenic carbon in the model occurs from the influence of rising atmospheric CO 2 on photosynthesis in trees, and thus, model-projected carbon feedbacks are especially sensitive to deforestation.

58 GEOSCIENCES↗

Cyber-intrusion Auto-response and Policy Management System (CAPMS)

The Cyber-intrusion Auto-response and Policy Management System (CAPMS) project was funded by a grant from the US Department of Energy (DOE) Cybersecurity for Energy Delivery Systems (CEDS) program with contributions from two partner electric utilities: Southern California Edison (SCE) and Duke Energy. The goal of the project was to demonstrate protecting smart grid assets from a cyber attack in a way that “does not impede critical energy delivery functions.” This report summarizes project goals and activities for the CAPMS project and explores what did and did not work as expected. It concludes with an assessment of possible benefits and value of the system for the future.

97 MATHEMATICS AND COMPUTING↗

A new synoptic scale resolving global climate simulation using the Community Earth System Model

Abstract High‐resolution global climate modeling holds the promise of capturing planetary‐scale climate modes and small‐scale (regional and sometimes extreme) features simultaneously, including their mutual interaction. This paper discusses a new state‐of‐the‐art high‐resolution Community Earth System Model (CESM) simulation that was performed with these goals in mind. The atmospheric component was at 0.25° grid spacing, and ocean component at 0.1°. One hundred years of “present‐day” simulation were completed. Major results were that annual mean sea surface temperature (SST) in the equatorial Pacific and El‐Niño Southern Oscillation variability were well simulated compared to standard resolution models. Tropical and southern Atlantic SST also had much reduced bias compared to previous versions of the model. In addition, the high resolution of the model enabled small‐scale features of the climate system to be represented, such as air‐sea interaction over ocean frontal zones, mesoscale systems generated by the Rockies, and Tropical Cyclones. Associated single component runs and standard resolution coupled runs are used to help attribute the strengths and weaknesses of the fully coupled run. The high‐resolution run employed 23,404 cores, costing 250 thousand processor‐hours per simulated year and made about two simulated years per day on the NCAR‐Wyoming supercomputer “Yellowstone.”

Small, R. Justin↗

GLACE: The Global Land-Atmosphere Coupling Experiment Part 2: Analysis

The twelve weather and climate models participating in the Global Land-Atmosphere Coupling Experiment (GLACE) show both a wide variation in the strength of land-atmosphere coupling and some intriguing commonalities. In this paper, we address the causes of variations in coupling strength - both the geographic variations within a given model and the model-to-model differences. The ability of soil moisture to affect precipitation is examined in two stages, namely, the ability of the soil moisture to affect evaporation, and the ability of evaporation to affect precipitation. Most of the differences between the models and within a given model are found to be associated with the first stage - an evaporation rate that varies strongly and consistently with soil moisture tends to lead to a higher coupling strength. The first stage differences reflect identifiable differences in model parameterization and model climate. Intermodel differences in the evaporation-precipitation connection, however, also play a key role.

Guo, Zhichang↗

"Hot Spots" of Land Atmosphere Coupling

Previous estimates of land-atmosphere interaction (the impact of soil moisture on precipitation) have been limited by a severe paucity of relevant observational data and by the model-dependence of the various computational estimates. To counter this limitation, a dozen climate modeling groups have recently performed the same highly-controlled numerical experiment as part of a coordinated intercomparison project. This allows, for the first time ever, a superior multi-model approach to the estimation of the regions on the globe where precipitation is affected by soil moisture anomalies during Northern Hemisphere summer. Such estimation has many potential benefits; it can contribute, for example, to seasonal rainfall prediction efforts.

Koster, Randal D.↗