Engineering Papers⌕ Search

Engineering topics

Chen, Y

Publications and source records attributed to Chen, Y.

At least 19 records

First-order preemptive Ising-nematic transition in K5Fe4Ag6Te10

Employing inelastic x-ray scattering and neutron scattering techniques, we observed nematic and magnetic phase transitions with distinct characters in K5Fe4Ag6Te10. Upon cooling, the nematic order undergoes a strongly first-order phase transition followed by a magnetic transition at TN≈34.7K that is consistent with a second-order transition. The temperature difference between these two phase transitions is ∼1.0(8) K. The observed phenomenon can be attributed to a distinctive first-order preemptive Ising-nematic transition, a characteristic unique to a quasi-two-dimensional scenario marked by strong out-of-plane spatial anisotropy due to weak coupling. Our studies establish K5Fe4Ag6Te10 as the first material in the family of iron pnictides and chalcogenides that possesses a nematic tricritical point preceding the magnetic one upon decreasing nematic coupling.

Giles-Donovan, N↗

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Center for Frontiers of Subsurface Energy Security (Final Report)

For the past two decades a slowly growing consensus has been that rising temperature is being caused by increases in carbon dioxide (CO 2 ) in the Earth’s atmosphere. What is less consensual is what to do about it. Solutions range from reducing emissions to various ways to capture CO 2 (from effluents of power plants, as plants, etc.) and store it away from the atmosphere. The problem is of such magnitude that it is likely that we will need many methods. Also some methods would cause draconian harm to global economies. This work is about the storage part of the reducing of CO 2 in the atmosphere. There have been many methods proposed for storing CO 2 : in vegetation, at the bottom of oceans, in land farms, and in subsurface formations. The latter is particularly prominent given the large estimated volume in porous formation in the Earth, and the 50-year old experience in injecting CO 2 for enhanced oil recovery. This report is about the science behind subsurface storage, called here geological carbon storage or GCS. Effective GCS requires that billions of metric tons of CO 2 be injected into storage reservoirs annually. GCS feasibility requires balancing the opposing impacts of large emplacement (injection) rates without compromising caprock seal integrity during injection. Once injected, the challenge is to assure the injected CO 2 remains permanently stored. Under most conditions, CO 2 is in a supercritical state, or scCO 2 . scCO 2 is immiscible with the subsurface brine and has a factor of ten smaller viscosity, leading to the development of two-phase fluid flow and scCO 2 bypassing. After the injection of large volumes of scCO 2 , the pore pressure, state of stress, chemical, thermal, and biological steady-state or equilibrium conditions in the subsurface are disrupted, which may lead to unpredictable behavior because the responses to these changes are often non-linear. Current challenges for efficient and reliable scCO 2 storage are: (1) sustaining large injection rates; (2) using available pore space efficiently, and (3) controlling undesired or unexpected behavior, such as scCO 2 leakage. Various aspects of these challenges are addressed by collaborative projects within Center for Frontiers in Subsurface Energy Security (CFSES). Challenge 1 – Sustaining large storage rates. To offset the annual CO 2 emissions of in the United States, ~7 billion metric tons of CO 2 must be injected into the storage reservoirs annually. Injected CO 2 forms a plume that starts spreading, displacing originally resident brine, though inefficiently so. Because of its low viscosity, and buoyant forces, a CO 2 plume rises to the reservoir rock-caprock interface. Injecting large volumes of scCO 2 at sufficiently high rates into subsurface storage formations leads to increases in pore pressures, potential expansion of the reservoir rock, and may result in fracturing of the reservoir rock and/or caprock. CFSES research uses existing and new experimental and modeling approaches to identify chemical-physical controls on the permeability and pore pressure dynamics, anticipating geomechanical (fracturing) events, predicting multi-phase flow patterns and trapping of CO 2 . Challenge 2 – Using pore space with unprecedented efficiency. Since large volumes of scCO 2 are to be stored, the volume available in storage reservoirs must be used efficiently. Storage efficiency is defined as a fraction of pore space occupied by carbon dioxide. Current estimates are that less than 5 percent of available pore volume is available for storage because scCO 2 , having a smaller density and viscosity than reservoir brine, tends to form viscous fingers, and bypass available pore space. This flow regime results in a CO 2 plume spreading over large areas in subsurface. The research goal in Challenge 2 is to advance the fundamental understanding of two-phase flow through porous media, including molecular-to meter- scale experimental and modeling approaches. This work leads to novel strategies to ensure that injected CO 2 occupies more than 5% of the pore space, ideally up to 50%. Challenge 3 – Controlling undesired or unexpected behavior. Following the injection of scCO 2 , the subsurface reservoir re-equilibrates to different physical and chemical states. This re-equilibration in space and time could lead to unexpected and/or undesired behavior: fracturing of reservoir rock and/or caprock, development of preferential flow paths, CO 2 leakage (to the surface or overlaying fresh water reservoirs), induced seismicity, and activation of faults. One of the main research goals of CFSES was to advance the fundamental understanding of how individual chemical and physical processes are coupled. We particularly focused on the chemical-mechanical coupling, and flow-mechanical coupling. We collected experimental and modeling data to determine the temporal and spatial scales of the coupled processes, which may lead to pore collapse caused by creep, activation of fractures and faults, chemically-induced fracture propagation, and self-focusing of scCO 2 flow.

54 ENVIRONMENTAL SCIENCES↗

The BayREN Integrated Commercial Retrofits (BRICR) Project: An Introduction and Preliminary Results

BayREN Integrated Commercial Retrofits (BRICR) is a DOE-funded project which aims to enhance the capacity of energy efficiency programs to recruit participants, develop retrofits, and measure outcomes in small and medium-sized commercial buildings, a sector notoriously hard to reach and expensive to serve, that accounts for 2/3 of US commercial floor space. To address these barriers, BRICR leverages existing incentives, financing, data, and open source software to facilitate two paths for comprehensive improvements: a deep energy retrofit, or serial upgrades integrated into capital improvement and maintenance cycles. BRICR is developing an integrated workflow for iterative energy modeling of commercial buildings for city energy program managers and auditors - starting with mass building-scale simulation based on public records and proceeding through audit, retrofit, and measurement and verification stages. BRICR builds on existing tools including LBNL’s CBES, NREL’s OpenStudio, PNNL’s Audit Template, and DOE’s BuildingSync and SEED PlatformTM. At each stage, BRICR uses available information to inform simulations (starting with public records but iteratively augmented with observations from energy program staff) to improve the quality of the models that inform decision making. This paper presents initial results from energy models of 1699 office and retail buildings in San Francisco. Building stock data from public records were translated to create a BuildingSync file for each building which was stored in SEED. Each BuildingSync file was then translated to multiple OpenStudio Workflow files for EnergyPlus simulation to estimate energy savings of energy conservation measures (ECM). Energy savings predictions for each ECM were written back to an updated BuildingSync file for each building and re-uploaded to SEED. The distribution of baseline energy simulations was calibrated against the publicly disclosed distribution of energy benchmarking data to increase confidence in results.

Piette, M↗

Study of the wave packet treatment of neutrino oscillation at Daya Bay

The disappearance of reactor $\bar{\nu}_e$ observed by the Daya Bay experiment is examined in the framework of a model in which the neutrino is described by a wave packet with a relative intrinsic momentum dispersion $\sigma_\text{rel}$. Three pairs of nuclear reactors and eight antineutrino detectors, each with good energy resolution, distributed among three experimental halls, supply a high-statistics sample of $\bar{\nu}_e$ acquired at nine different baselines. This provides a unique platform to test the effects which arise from the wave packet treatment of neutrino oscillation. The modified survival probability formula was used to fit Daya Bay data, providing the first experimental limits: $2.38 \cdot 10^{-17} < \sigma_{\rm rel} < 0.23$. Treating the dimensions of the reactor cores and detectors as constraints, the limits are improved: $10^{-14} \lesssim \sigma_{\rm rel} < 0.23$, and an upper limit of $\sigma_{\rm rel} <0.20$ is obtained. All limits correspond to a 95\% C.L. Furthermore, the effect due to the wave packet nature of neutrino oscillation is found to be insignificant for reactor antineutrinos detected by the Daya Bay experiment thus ensuring an unbiased measurement of the oscillation parameters $\sin^22\theta_{13}$ and $\Delta m^2_{32}$ within the plane wave model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-F-P-54: Guidelines to Check Image Registration QA of a Clinical Deformation Registration Software: A Single Institution Preliminary Study

Purpose: The objective of this study is to verify and analyze the accuracy of a clinical deformable image registration (DIR) software. Methods: To test clinical DIR software qualitatively and quantitatively, we focused on lung radiotherapy and analyzed a single (Lung) patient CT scan. Artificial anatomical changes were applied to account for daily variations during the course of treatment including the planning target volume (PTV) and organs at risk (OAR). The primary CT (pCT) and the structure set (pST) was deformed with commercial tool (ImSimQA-Oncology Systems Limited) and after artificial deformation (dCT and dST) sent to another commercial tool (VelocityAI-Varian Medical Systems). In Velocity, the deformed CT and structures (dCT and dST) were inversely deformed back to original primary CT (dbpCT and dbpST). We compared the dbpST and pST structure sets using similarity metrics. Furthermore, a binary deformation field vector (BDF) was created and sent to ImSimQA software for comparison with known “ground truth” deformation vector fields (DVF). Results: An image similarity comparison was made by using “ground truth” DVF and “deformed output” BDF with an output of normalized “cross correlation (CC)” and “mutual information (MI)” in ImSimQA software. Results for the lung case were MI=0.66 and CC=0.99. The artificial structure deformation in both pST and dbpST was analyzed using DICE coefficient, mean distance to conformity (MDC) and deformation field error volume histogram (DFEVH) by comparing them before and after inverse deformation. We have noticed inadequate structure match for CTV, ITV and PTV due to close proximity of heart and overall affected by lung expansion. Conclusion: We have seen similarity between pCT and dbpCT but not so well between pST and dbpST, because of inadequate structure deformation in clinical DIR system. This system based quality assurance test will prepare us for adopting the guidelines of upcoming AAPM task group 132 protocol.

60 APPLIED LIFE SCIENCES↗

SU-F-I-14: 3D Breast Digital Phantom for XACT Imaging

Purpose: The X-ray induced acoustic computed tomography (XACT) is a new imaging modality which combines X-ray contrast and high ultrasonic resolution in a single modality. Using XACT in breast imaging, a 3D breast volume can be imaged by only one pulsed X-ray radiation, which could dramatically reduce the imaging dose for patients undergoing breast cancer screening and diagnosis. A 3D digital phantom that contains both X-ray properties and acoustic properties of different tissue types is indeed needed for developing and optimizing the XACT system. The purpose of this study is to offer a realistic breast digital phantom as a valuable tool for improving breast XACT imaging techniques and potentially leading to better diagnostic outcomes. Methods: A series of breast CT images along the coronal plane from a patient who has breast calcifications are used as the source images. A HU value based segmentation algorithm is employed to identify breast tissues in five categories, namely the skin tissue, fat tissue, glandular tissue, chest bone and calcifications. For each pixel, the dose related parameters, such as material components and density, and acoustic related parameters, such as frequency-dependent acoustic attenuation coefficient and bandwidth, are assigned based on tissue types. Meanwhile, other parameters which are used in sound propagation, including the sound speed, thermal expansion coefficient, and heat capacity are also assigned to each tissue. Results: A series of 2D tissue type image is acquired first and the 3D digital breast phantom is obtained by using commercial 3D reconstruction software. When giving specific settings including dose depositions and ultrasound center frequency, the X-ray induced initial pressure rise can be calculated accordingly. Conclusion: The proposed 3D breast digital phantom represents a realistic breast anatomic structure and provides a valuable tool for developing and evaluating the system performance for XACT.

60 APPLIED LIFE SCIENCES↗

SU-F-I-51: CT/MR Image Deformation: The Clinical Assessment QA in Target Delineation

Purpose: To study the deformation effects in CT/MR image registration of head and neck (HN) cancers. We present a clinical indication in guiding and simplifying registration procedures of this process while CT images possessed artifacts. Methods: CT/MR image fusion provides better soft tissue contrast in intracranial GTV definition with artifacts. However, whether the fusion process should include the deformation process is questionable and not recommended. We performed CT/MR image registration of a HN patient with tonsil GTV and nodes delineation on Varian Velocity™ system. Both rigid transformation and deformable registration of the same CT/MR imaging data were processed separately. Physician’s selection of target delineation was implemented to identify the variations. Transformation matrix was shown with visual identification, as well as the deformation QA numbers and figures were assessed. Results: The deformable CT/MR images were traced with the calculated matrix, both translation and rotational parameters were summarized. In deformable quality QA, the calculated Jacobian matrix was analyzed, which the min/mean/max of 0.73/0/99/1.37, respectively. Jacobian matrix of right neck node was 0.84/1.13/1.41, which present dis-similarity of the nodal area. If Jacobian = 1, the deformation is at the optimum situation. In this case, the deformation results have shown better target delineation for CT/MR deformation than rigid transformation. Though the root-mean-square vector difference is 1.48 mm, with similar rotational components, the cord and vertebrae position were aligned much better in the deformable MR images than the rigid transformation. Conclusion: CT/MR with/without image deformation presents similar image registration matrix; there were significant differentiate the anatomical structures in the region of interest by deformable process. Though vendor suggested only rigid transformation between CT/MR assuming the geometry remain similar, our findings indicated with patient positional variations, deformation registration is needed to generate proper GTV coverage, which will be irradiated more accurately in the following boost phase.

60 APPLIED LIFE SCIENCES↗

SU-F-I-69: Absorbed Dose Estimation for a Commercially Available MicroCT Scanner with Various Filtration Techniques

Purpose: To quantify the simulated mean absorbed dose per technique (cGy/mAs) from a commercially available microCT scanner using various filtration techniques. Methods: Monte Carlo simulations using the Geant4 toolkit (version 10) was used utilizing the standard electromagnetic physics model. The Quantum FX microCT scanner (PerkinElmer, Waltham, MA) was modeled incorporating measured energy spectra and spatial dimensions of nominal source-to-object (SOD) distances. The energy distribution was measured using a spectrometer (X-123CdTe, Amptek Inc., Bedford, USA) for the 90 kVp X-ray beams with various filters (including no filter, 1 mm, 2 mm, 3 mm, 4 mm Al and 0.2 mm Cu + 2.5 mm Al). The SOD was set to 154 mm, 104 mm, and 52 mm. A total of 10 million incident particles were processed per simulation. Cutout value was set to 0.1 mm for both photon and electron. The mean dose absorbed (cGy/per incident particle) in a PMMA phantom (length of 2 cm and radius of 3 cm) were recorded. Exposure measurements were taken using a Radcal 9095 system with a protocol of 90 kVp, 200 µA, and ∼12 s beam-On time for the various filters. Results: The mean absorbed dose per mAs for various filtrations and different SOD setups indicated that the dose decreased as the SOD increased and as the amount of filtration increased. For a given SOD, the dose was reduced by as much as ∼13.7% by varying the filter (from 0.2 mm Cu + 2.5 mm Al to no filter). The maximum dose was found to be 0.39 cGy/mAs (SOD of 5.196 cm, no filter) while the minimum dose value was 0.077 cGy/mAs (SOD of 15.4 cm, .2mm Cu + 2mm Al filter). Conclusion: This study estimates easily the mean dose for objects scanned with a microCT scanner with different filtration.

60 APPLIED LIFE SCIENCES↗