Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Performance Rating Method”

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 OSTI2024

Model Quality and Measurement Density Impact on Volt/Volt Ampere Reactive Optimization Performance

The operation of the utility grid is being reshaped by the continuous addition of distributed energy resources and advanced metering infrastructure, which challenge existing grid control strategies. Some utilities deploy advanced distribution management systems (ADMS) to assist with the consolidation of various applications and to augment situational awareness in response to the new power delivery dynamics. An ADMS is an integrated software platform that provides utilities with a way to enhance their reliability, control, and optimization with advanced applications, such as volt/VAR optimization (VVO). A VVO application could serve as a vehicle to deliver cost savings by providing the utility with a method to reduce rates by controlling the voltage and decreasing the energy usage in their service territory. Some utilities are reluctant to integrate an ADMS, because it is a significant investment that requires approval from the public regulatory commission and/or their customers. This paper evaluates the impact on VVO performance when using a lower-quality network model supplemented with additional measurements, which could provide an implementation for cost savings. The results show that a better model quality would provide the highest energy savings; however, some level of telemetry is necessary in all scenarios to prevent voltage exceedances.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Analysis of Thermal and Chemical Effets on Negative Valve Overlap Period Energy Recovery for Low-Temperature Gasoline Combustion

A central challenge for efficient auto-ignition controlled low-temperature gasoline combustion (LTGC) engines has been achieving the combustion phasing needed to reach stable performance over a wide operating regime. The negative valve overlap (NVO) strategy has been explored as a way to improve combustion stability through a combination of charge heating and altered reactivity via a recompression stroke with a pilot fuel injection. The study objective was to analyze the thermal and chemical effects on NVO-period energy recovery. The analysis leveraged experimental gas sampling results obtained from a single-cylinder LTGC engine along with cylinder pressure measurements and custom data reduction methods used to estimate period thermodynamic properties. The engine was fueled by either iso-octane or ethanol, and operated under sweeps of NVO-period oxygen concentration, injection timing, and fueling rate. Gas sampling at the end of the NVO period was performed via a custom dump-valve apparatus, with detailed sample speciation by in-house gas chromatography. The balance of NVO-period input and output energy flows was calculated in terms of fuel energy, work, heat loss, and change in sensible energy. Experiment results were complemented by detailed chemistry single-zone reactor simulations performed at relevant mixing and thermodynamic conditions, with results used to evaluate ignition behavior and expected energy recovery yields. For the intermediate bulk-gas temperatures present during the NVO period (900-1100 K), weak negative temperature coefficient behavior with iso-octane fueling significantly lengthened ignition delays relative to similar ethanol fueled conditions. Faster ethanol ignition chemistry led to lower recovered fuel intermediate yields relative to similar iso-octane fueled conditions due to more complete fuel oxidation. From the energy analysis it was found that increased NVO-period global equivalence ratio, either from lower NVOperiod oxygen concentrations or higher fueling rates, in general led to a greater fraction of net recovered fuel energy and work as heat losses were minimized. These observations were supported by complementary single-zone reactor model results, which further indicated that kinetic time-scales favor chemical energy-consuming exothermic oxidation over slower endothermic reformation. Nonetheless, fuel energy recovery close to the thermodynamic equilibrium solution was achieved for baseline conditions that featured 4% NVO-period oxygen concentration.

ethanol↗
DOE OSTI2022

Water Imbibition and Oil Recovery in Shale: Dynamics and Mechanisms Using Integrated Centimeter-to-Nanometer-Scale Imaging

Water imbibition, and the associated oil displacement, is an important process in shale oil reservoirs after hydraulic fracturing and in water-based enhanced oil recovery (EOR). Current techniques for water imbibition measurement are mostly “black-box”-type methods. A more explicit understanding of the water imbibition/oil recovery dynamics and geological controls is in demand. In this paper, a multiscale imaging technique that covers centimeter to nanometer scale (i.e., core to pore scale), integrating neutron radiography, microcomputed tomography (micro-CT), and scanning electron microscope (SEM) is applied to investigate the water imbibition depth and rate and the cause of heterogeneity of imbibition in shale samples. The dynamic processes of water imbibition in the 1-in. (25.4-mm) core sample were explicitly demonstrated, and the imbibition along the matrix and imbibition through microfractures are distinguished through neutron radiography image analysis. The causes of observed imbibition heterogeneity were further investigated through micro-CT and SEM image analysis for 1.5-mm diameter miniplug samples from different laminas of the 1-in. core samples. Imbibition depth and rate were calculated on the basis of image analysis as well. Estimation of oil recovery through water imbibition in shale matrix was performed for an example shale field. Here this innovative and integrated multiscale imaging technique provides a “white/gray-box” method to understand water imbibition and water-oil displacement in shale. The wide span of the length scale (from centimeter to nanometer) of this technique enables a more comprehensive, accurate, and specific understanding of both the core-scale dynamics and pore-scale mechanisms of water imbibition, oil recovery, and matrix-fracture interaction.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

High Throughput Source-less Plasma Deposition of Structured Silicon Anodes for Lithium-Ion Batteries

Amprius developed a manufacturing solution for silicon nanowire anode that relies on an inexpensive, high throughput, and high gas precursor utilization plasma deposition method that uses the anode foils as electrodes for plasma generation. The capacitively couple plasma (CCP) method is used in semiconductor and photovoltaic industry and Amprius modified existing high throughput equipment to use anode foils and to deposit amorphous silicon. The equipment was installed ahead of the program at Amprius site. The rest of the tasks included foil handling and process development. The equipment passed site acceptance tests (SAT) and the process parameter mapping was completed, indicating that the target process window limits produce output materials at the rate and with yield and specifications that meet the manufacturing target criteria. Amprius has hired supporting personnel to optimize processes and run the equipment. A parallel task verified the baseline performance of the silicon anode material, to be used as reference for the new manufacturing method.

25 ENERGY STORAGE↗
DOE OSTI2022

Quantification of Off-Channel Inundated Habitat for Pacific Chinook Salmon ( Oncorhynchus tshawytscha ) along the Sacramento River, California, Using Remote Sensing Imagery

Off-channel areas are one of the most impacted aquatic habitats by humans globally, as extensive agricultural and urban development has limited them to roughly 10% of historical extent. This is also true for California’s Sacramento River Valley, where historically frequent widespread inundation has been reduced to a few off-channel water bodies along the mid-Sacramento River. This remaining shallow-water habitat provides crucial ecological benefits to multiple avian and fish species, but especially to floodplain-adapted species such as Chinook salmon (Oncorhynchus tshawytscha). Characterizing spatiotemporal off-channel dynamics, including inundation extent and residence time, is fundamental to better understanding the intrinsic value of such habitats and their potential to support recovery actions. Remote sensing techniques have been increasingly used to map surface water at regional and local scales, with improved resolutions. As such, this study maps off-channel inundation areas and describes their temporal dynamics by analyzing pixel-based time- series of multiple water indices, modified Normalized Difference Water Index (mNDWI) and the Automated Water Extraction Index (AWEI), generated from LandSat-8 and Sentinel-2 data between 2013–2021. Quantified off-channel area was similar with each water index and method used, but improved performance was associated with Sentinel-2 products and AWEI index to identify wetted areas under lower mainstem discharges. Results indicate an uneven distribution of off-channel habitat in the study area, with limited inundated areas in upstream reaches (<16% of total off-channel area for greater flows). In addition, much less habitat exists for flows under 400 m 3 /s, an important migration cue for endangered winter-run Chinook salmon, limiting juvenile access to areas with enhanced rearing conditions. Off-channel habitat residence times averaged between 7 and 16 days, primarily defined by the rate of receding flows, with rapid flow recession providing marginal off-channel habitat. This study shows reasonable performance of moderate resolution LandSat-8 and Sentinel-2 remote sensing imagery to characterize shallow-water inundated habitat in higher-order rivers, and as a method to inform restoration and native fish recovery efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

An Easy Tool to Predict Survival in Patients Receiving Radiation Therapy for Painful Bone Metastases

Purpose: Patients with bone metastases have a widely varying survival. A reliable estimation of survival is needed for appropriate treatment strategies. Our goal was to assess the value of simple prognostic factors, namely, patient and tumor characteristics, Karnofsky performance status (KPS), and patient-reported scores of pain and quality of life, to predict survival in patients with painful bone metastases. Methods and Materials: In the Dutch Bone Metastasis Study, 1157 patients were treated with radiation therapy for painful bone metastases. At randomization, physicians determined the KPS; patients rated general health on a visual analogue scale (VAS-gh), valuation of life on a verbal rating scale (VRS-vl) and pain intensity. To assess the predictive value of the variables, we used multivariate Cox proportional hazard analyses and C-statistics for discriminative value. Of the final model, calibration was assessed. External validation was performed on a dataset of 934 patients who were treated with radiation therapy for vertebral metastases. Results: Patients had mainly breast (39%), prostate (23%), or lung cancer (25%). After a maximum of 142 weeks' follow-up, 74% of patients had died. The best predictive model included sex, primary tumor, visceral metastases, KPS, VAS-gh, and VRS-vl (C-statistic = 0.72, 95% CI = 0.70-0.74). A reduced model, with only KPS and primary tumor, showed comparable discriminative capacity (C-statistic = 0.71, 95% CI = 0.69-0.72). External validation showed a C-statistic of 0.72 (95% CI = 0.70-0.73). Calibration of the derivation and the validation dataset showed underestimation of survival. Conclusion: In predicting survival in patients with painful bone metastases, KPS combined with primary tumor was comparable to a more complex model. Considering the amount of variables in complex models and the additional burden on patients, the simple model is preferred for daily use. In addition, a risk table for survival is provided.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Stereotactic Body Radiation Therapy for Locally Advanced Pancreatic Cancer: A Systematic Review and Pooled Analysis of 19 Trials

Purpose: Although surgery is the standard of care for resectable pancreatic cancer (PC), standard-dose chemoradiation therapy and chemotherapy alone are suitable for patients with unresectable disease. Stereotactic body radiation therapy (SBRT) is an alternative, focused local therapy that delivers high radiation doses within a few fractions to the cancer, sparing the surrounding critical tissue. We performed a systematic review and pooled analysis of published trials to evaluate the efficacy and safety of this emerging treatment modality. Methods and Materials: We searched the Cochrane Central Register of Controlled Trials, PubMed, EMBASE, SCOPUS, the Web of Science, and CINAHL for publications regarding SBRT for locally advanced PC. The 1-year overall survival (OS) rate was the primary endpoint, and the median OS, 2-year OS rate, 1-year locoregional control (LRC) rate, and grade 3 to 4 toxicities were the secondary endpoints. A multivariate random-effects meta-analysis was performed to calculate the aggregated OS rates at 1 and 2 years and the 1-year LRC rate. Results: A total of 19 studies, encompassing 1009 patients, were included in the present analysis. The pooled 1-year OS was 51.6% in 13 trials with data available. The median OS ranged from 5.7 to 47 months (median 17). The LRC rate at 1 year was 72.3%. Overall, the occurrence of severe adverse events did not exceed 10%. LRC appeared to correlate with the total SBRT dose and the number of fractions. Conclusions: The advantages of SBRT in terms of treatment time, satisfactory OS, and LRC indicate that it is an effective option for inoperable PC. However, a definitive validation of this treatment modality in large randomized studies is required, owing to the nonrandomized nature of the included studies and the limitations of small single-center series that include mixed populations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Interpreting T-History Data on Supercooling and Phase Change Enthalpy of a Salt-Hydrogel Complex for Building Thermal Management

Phase change materials (PCMs) are promising materials for storing and discharging latent heat in buildings to regulate the thermal environment. Measurement of salt hydrate PCM thermal performance is often challenging because of supercooling and phase segregation. Though differential scanning calorimetry (DSC) has been widely used to measure thermal properties of various phase change materials, DSC uses milligrams of sample and is unrealistic for predicting large-scale materials performance. The T-history method [1], in contrast, handles higher volumes but the analysis typically requires a lumped capacitance assumption. Here, we report and analyze T-History data on a novel salt-hydrogel complex with high viscosity, that renders low Biot number sample preparation extremely challenging. We re-examine the standard T-history method and compared against DSC data in this complex material. We report details of an experimental setup, built with necessary additions to a controlled environmental chamber. We carefully examine the validity of the lumped capacitance assumption and suggest a modified analysis using a computational approach for reliable data. Experimental parameters such as ramp rate and temperature range of the environmental chamber are considered via computational modeling to understand the effect of ramp rate on supercooling [2]. Finally, we report thermal cycling experiments on the salt-hydrogel complex. Our modified approach to measuring supercooling and enthalpy of fusion at large scales is important in better understanding the performance of phase-change materials at scale for building thermal storage. References: [1] Marin, Jose M, Belen Zalba, Luisa F Cabeza, and Harald Mehling. “Determination of Enthalpy Temperature Curves of Phase Change Materials with the Temperature-History Method: Improvement to Temperature Dependent Properties.” Measurement Science and Technology 14, no. 2 (February 1, 2003): 184–89. https://doi.org/10.1088/0957-0233/14/2/305. [2] Safari, A., R. Saidur, F.A. Sulaiman, Yan Xu, and Joe Dong. “A Review on Supercooling of Phase Change Materials in Thermal Energy Storage Systems.” Renewable and Sustainable Energy Reviews 70 (April 2017): 905–19. https://doi.org/10.1016/j.rser.2016.11.272

thermodynamics↗
DOE OSTI2016

Analysis of Monolith Cores from an Engineering Scale Demonstration of a Prospective Cast Stone Process

The primary disposition path of Low Activity Waste (LAW) at the DOE Hanford Site is vitrification. A cementitious waste form is one of the alternatives being considered for the supplemental immobilization of the LAW that will not be treated by the primary vitrification facility. Washington River Protection Solutions (WRPS) has been directed to generate and collect data on cementitious or pozzolanic waste forms such as Cast Stone. This report documents the coring and leach testing of monolithic samples cored from an engineering-scale demonstration (ES Demo) with non-radioactive simulants. The ES Demo was performed at SRNL in October of 2013 using the Scaled Continuous Processing Facility (SCPF) to fill an 8.5 ft. diameter x 3.25 ft. high container with simulated Cast Stone grout. The Cast Stone formulation was chosen from the previous screening tests. Legacy salt solution from previous Hanford salt waste testing was adjusted to correspond to the average LAW composition generated from the Hanford Tank Waste Operation Simulator (HTWOS). The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS), were obtained from Lafarge North America in Pasco, WA. In 2014 core samples originally obtained approximately six months after filling the ES Demo were tested along with bench scale molded samples that were collected during the original pour. A latter set of core samples were obtained in late March of 2015, eighteen months after completion of the original ES Demo. Core samples were obtained using a 2” diameter x 11” long coring bit. The ES Demo was sampled in three different regions consisting of an outer ring, a middle ring and an inner core zone. Cores from these three lateral zones were further segregated into upper, middle and lower vertical segments. Monolithic core samples were tested using the Environmental Protection Agency (EPA) Method 1315, which is designed to provide mass transfer rates (release rates) of inorganic analytes contained in monolithic material under diffusion controlled release conditions as a function of leaching time. Compressive strength measurements and drying tests were also performed on the 2015 samples. Leachability indices reported are based on analyte concentrations determined from dissolution of the dried samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Isotopic exchange on solid-phase micro extraction fiber in sediment under stagnant conditions: Implications for field application of performance reference compound calibration

An overlooked issue for field application of in situ performance reference compound (PRC) calibration methods is the validity of the assumption that both the sorption of a target compound and desorption of its corresponding PRC follow the first-order kinetics with the same rate constants under stagnant conditions. In the present study, disposable polydimethylsiloxane fibers of 2 sizes (7 and 35 µm) impregnated with 8 13 C-labeled or deuterated PRCs were statically deployed into different marine sediments, from which the kinetics for sorption of the target compounds and desorption of the PRCs were characterized. Nonsymmetrical profiles were observed for exchange of the target analytes and their corresponding PRCs in sediment under stagnant conditions. Further, the hysteretic desorption of PRCs in the kinetic regime may be ascribed to the low chemical potential between the fiber and sediment porewater, which reflects the inability of water molecules to rapidly diffuse through sediment to solvate the PRCs in the aqueous layer around the fiber surface. A moderate correlation (r = 0.77 and r = 0.57, p < 0.05 for both regressions) between the PRC-calibrated equilibrium concentrations of 1,1-dichloro-2,2-bis-(chlorophenyl) ethylene (p,p'-DDE) and polychlorinated biphenyl (PCB)-153 and the lipid normalized levels in worms (Neanthes arenaceodentata) was obtained in co-exposure tests under simulating field conditions, probably resulting from slightly overestimated bioavailability because of the hysteretic desorption of PRCs and toxic effects. Environ Toxicol Chem 2016;35:1978–1985. © 2015 SETAC

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Ionic Liquid–Directed Nanoporous TiNb 2 O 7 Anodes with Superior Performance for Fast–Rechargeable Lithium–Ion Batteries

Nanoporous TiNb 2 O 7 (NPTNO) material is synthesized by a sol–gel method with an ionic liquid (IL) as the nanoporous structure directing template. NPTNO exhibits a high reversible capacity of 210 mAh g –1 even at the charging rate of 50 C and an excellent cyclability of half-cell capacity retention of 74% for 1000 cycles at 5 C and LiNi 0.5 Mn 1.5 O 4 -coupled full-cell capacity retentions of 81% and 87% for 1000 cycles at 1 C and 2 C, respectively. The studies of the 1000 cycled NPTNO electrode illustrate that the IL-directed mesoporous structure can enhance the cyclability of NPTNO cells due to the alleviation of repetitive mechanical stress and volume fluctuation induced by the repetitive Li + insertion-extraction processes. The measured Li + diffusion coefficients from the galvanostatic intermittent titration technique suggest that the IL-templating strategy indeed ensures the fast rechargeability of NPTNO cells based on the fast Li + diffusion kinetics. Benefitting from the nanoporous structure, NPTNO with unhindered Li + diffusion pathways achieves a superior rate capability in the titanium-based oxide materials and the best full-cell cyclability in the TNO materials. Furthermore, the templating potential of IL is demonstrated, and the superb electrochemical performance establishes the IL-directed NPTNO as a promising anode candidate for fast-rechargeable LIBs.

25 ENERGY STORAGE↗
DOE OSTI2018

Uncovering the Cu-driven electrochemical mechanism of transition metal chalcogenides based electrodes

Transition-metal chalcogenides (TMCs) have emerged as attractive anode materials for rechargeable batteries due to their excellent performance and abundant resources. In this work, for the first time, we disclose a unique copper (Cu)-driven conversion process in TMC-based battery systems that involves classic Cu current collector and is considered to be an “activation process”. According to state-of-the-art characterization techniques, Cu was evidenced to gradually replace the transition-metal elements in TMCs to be the active material during cycling. Based on this unique Cu-driven conversion mechanism, we used a facile method to design a new type of sulfur-based battery that presents excellent performance: a reversible capacity of 1.045 mAh cm -2 after 700 cycles at 2Ag -1 , and a good rate capability up to a capacity of 0.33 mAh cm -2 at 20 A g - 1. Lastly, with respect to the large family of TMC compounds, this study introduces a new direction for the design of high-performance energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2013

High-Dose-Rate Prostate Brachytherapy Consistently Results in High Quality Dosimetry

Purpose: We performed a dosimetry analysis to determine how well the goals for clinical target volume coverage, dose homogeneity, and normal tissue dose constraints were achieved with high-dose-rate (HDR) prostate brachytherapy. Methods and Materials: Cumulative dose-volume histograms for 208 consecutively treated HDR prostate brachytherapy implants were analyzed. Planning was based on ultrasound-guided catheter insertion and postoperative CT imaging; the contoured clinical target volume (CTV) was the prostate, a small margin, and the proximal seminal vesicles. Dosimetric parameters analyzed for the CTV were D90, V90, V100, V150, and V200. Dose to the urethra, bladder, bladder balloon, and rectum were evaluated by the dose to 0.1 cm{sup 3}, 1 cm{sup 3}, and 2 cm{sup 3} of each organ, expressed as a percentage of the prescribed dose. Analysis was stratified according to prostate size. Results: The mean prostate ultrasound volume was 38.7 {+-} 13.4 cm{sup 3} (range: 11.7-108.6 cm{sup 3}). The mean CTV was 75.1 {+-} 20.6 cm{sup 3} (range: 33.4-156.5 cm{sup 3}). The mean D90 was 109.2% {+-} 2.6% (range: 102.3%-118.4%). Ninety-three percent of observed D90 values were between 105 and 115%. The mean V90, V100, V150, and V200 were 99.9% {+-} 0.05%, 99.5% {+-} 0.8%, 25.4% {+-} 4.2%, and 7.8% {+-} 1.4%. The mean dose to 0.1 cm{sup 3}, 1 cm{sup 3}, and 2 cm{sup 3} for organs at risk were: Urethra: 107.3% {+-} 3.0%, 101.1% {+-} 14.6%, and 47.9% {+-} 34.8%; bladder wall: 79.5% {+-} 5.1%, 69.8% {+-} 4.9%, and 64.3% {+-} 5.0%; bladder balloon: 70.3% {+-} 6.8%, 59.1% {+-} 6.6%, and 52.3% {+-} 6.2%; rectum: 76.3% {+-} 2.5%, 70.2% {+-} 3.3%, and 66.3% {+-} 3.8%. There was no significant difference between D90 and V100 when stratified by prostate size. Conclusions: HDR brachytherapy allows the physician to consistently achieve complete prostate target coverage and maintain normal tissue dose constraints for organs at risk over a wide range of target volumes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Thermodynamic and Kinetic Mechanisms Governing the Synthesis of Nickel-Poor Cathodes

A deeper understanding of the thermodynamics and kinetics governing the lithiation and layering mechanisms of NMC cathode materials (LiNixMnyCozO 2 , where x + y + z = 1) offers valuable insights for enhancing synthesis methods and improving cathode performance. By employing atomistic and mesoscale approaches informed by in situ powder X-ray diffraction (PXRD) experiments, critical parameters for comprehending lithiation and layering processes and reaction rates were identified. The mesoscale approach captured the evolution of the phases and crystallite size observed in the in situ PXRD, revealing the differences in reaction rates with the use of different lithium salts and starting precursors. Ab initio molecular dynamics (AIMD) underscored the importance of vacancies and structural defects in promoting ion mobility and facilitating the nucleation of a layered domain. This nucleation disrupts the symmetry of disordered phases, ultimately creating a strained phase that serves as a buffer between layered and disordered regions. The lithiation and layering processes reflect a dynamic balance between the thermodynamic drive for a low-energy layered structure and the kinetic of diffusion, which is influenced by temperature and lithium vacancy concentration. Overall, reaction mechanisms are driven by the inherent defects of the intermediate phase that differ for NMC cathode materials. The lithium salts impact the rates of lithiation and layering, with a much slower process for Li 2 CO 3 .

25 ENERGY STORAGE↗
DOE OSTI2012

Direct measurement of the {sup 11}C({alpha},p){sup 14}N reaction at CRIB: A path from pp-chain to CNO

We determined the total reaction rate of the {sup 11}C({alpha},p){sup 14}N reaction relevant to the nucleosynthesis in explosive hydrogen-burning stars. The measurement was performed by means of the thick target method in inverse kinematics with {sup 11}C RI beams. We performed the identification of the ground-state transition and excited-state transitions using time-of-flight information for the first time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

A framework for solving atomistic phonon-structure scattering problems in the frequency domain using perfectly matched layer boundaries

We present a numerical approach to the solution of elastic phonon-interface and phonon-nanostructure scattering problems based on a frequency-domain decomposition of the atomistic equations of motion and the use of perfectly matched layer (PML) boundaries. Unlike molecular dynamic wavepacket analysis, the current approach provides the ability to simulate scattering from individual phonon modes, including wavevectors in highly dispersive regimes. Like the atomistic Green's function method, the technique reduces scattering problems to a system of linear algebraic equations via a sparse, tightly banded matrix regardless of dimensionality. However, the use of PML boundaries enables rapid absorption of scattered wave energies at the boundaries and provides a simple and inexpensive interpretation of the scattered phonon energy flux calculated from the energy dissipation rate in the PML. The accuracy of the method is demonstrated on connected monoatomic chains, for which an analytic solution is known. The parameters defining the PML are found to affect the performance and guidelines for selecting optimal parameters are given. The method is used to study the energy transmission coefficient for connected diatomic chains over all available wavevectors for both optical and longitudinal phonons; it is found that when there is discontinuity between sublattices, even connected chains of equivalent acoustic impedance have near-zero transmission coefficient for short wavelengths. The phonon scattering cross section of an embedded nanocylinder is calculated in 2D for a wide range of frequencies to demonstrate the extension of the method to high dimensions. The calculations match continuum theory for long-wavelength phonons and large cylinder radii, but otherwise show complex physics associated with discreteness of the lattice. Examples include Mie oscillations which terminate when incident phonon frequencies exceed the maximum available frequency in the embedded nanocylinder, and scattering efficiencies larger than two near the Brillouin zone edge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Speeding up N -body simulations of modified gravity: chameleon screening models

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied f ( R ) gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of f ( R ) simulations. For example, a test simulation with 512{sup 3} particles in a box of size 512 Mpc/ h is now 5 times faster than before, while a Millennium-resolution simulation for f ( R ) gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

79 ASTRONOMY AND ASTROPHYSICS↗